[01]. [Tool][软件] OSMatTest: OpenSees Material Test [OpenSees单轴材料测试工具]

[02]. OpenSees Steel01 Material Test

[03]. OpenSees Steel02 Material Test

[04]. OpenSees RambergOsgoodSteel Material Test

[05]. OpenSees Concrete01 Material Test

[06]. OpenSees Concrete02 Material Test

[07]. OpenSees Concrete04 Material Test

[08]. OpenSees ConcreteD Material Test

[09]. OpenSees Concrete06 Material Test

[10]. [Tool][OSFSV] OSFSV v2017: OpenSees Fiber Section Viewer [OpenSees纤维截面可视化工具]

[11]. [科研][OpenSees] BWBN Material (Pinching Hysteretic Bouc-Wen Material) Test

[12]. [科研][OpenSees] BoucWen Material Test

[13]. [Tool][工具][科研][更新] OSFSV v2021: OpenSees Fiber Section Viewer [OpenSees纤维截面可视化工具]

[14]. [OSFSV][案例][OpenSees] OSFSV: OSFiberSectionViewer提示“sqrt: DOMAIN error”错误的解决办法!!

[15]. [笔记][非线性][OpenSees] 多个单轴二折线非线性本构并联是什么样子的?

[16]. [笔记] OpenSees中的瑞丽阻尼( Rayleigh Damping in OpenSees)

[17]. [科研][软件][OpenSees] OSBWTest: OpenSees BoucWen Material Test ( OpenSees BoucWen材料测试)

70 thoughts on “OpenSees 软件笔记

  1. yihong Reply

    崔济东博士您好,我用openseespy建立了一个平面框架,边柱与中柱强轴方向不一致,惯性矩不同,施加地震波得到的时程曲线后端是震荡的,请问这可能是什么原因呢

  2. Jinggang Reply

    崔博士,我想咨询下您知不知道opensees输出截面变形结果的各项内容代表的具体内容?因为我采用如下命令:
    recorder Element -file $dataDir/deformationPier1.txt -time -eleRange 157 157 section 1 deformation
    输出157号单元1#积分点的截面变形数据,但是输出结果有七项内容,有些搞不懂每一项具体代表什么内容,所以向您请教下,多谢!

  3. MR。宋 Reply

    想请教一下崔博士,以下的不收敛问题可能是什么原因,我检查了一下模型,可能是约束条件、但并没有发现什么可以修改的地方。(以下用了多个算法尝试收敛,但均无效,模型是梁柱节点,荷载分析是往复加载分析)
    WARNING: CTestEnergyIncr::test() – failed to converge
    after: 6 iterations
    current EnergyIncr: 395.293 (max: 1e-06) Norm deltaX: 0.122196, Norm deltaR: 3.15897e+06
    NewtonRaphson::solveCurrentStep() -the ConvergenceTest object failed in test()
    StaticAnalysis::analyze() – the Algorithm failed at iteration: 0 with domain at load factor -0.326601
    OpenSees > analyze failed, returned: -3 error flag
    Trying Newton with Initial Tangent ..
    WARNING: CTestNormDispIncr::test() – failed to converge
    after: 3000 iterations current Norm: 0.274599 (max: 1e-06, Norm deltaR: 953445)
    NewtonRaphson::solveCurrentStep() -the ConvergenceTest object failed in test()
    StaticAnalysis::analyze() – the Algorithm failed at iteration: 0 with domain at load factor 0.0949372
    OpenSees > analyze failed, returned: -3 error flag

    • CJD Post authorReply

      如果是约速条件,可以试试把本构设置为弹性 试试是否有同样错误。

    • qiao wang Reply

      老师好,我问个问题,如果想在opensees中去了解一下材料的本构,就是在opensees材料的文件中学习吗?比如说:我想了解ViscousDamper这种材料,是在ViscousDamper.cpp文件里查看吗?​
      然后如果是在在材料的基础上进行二次开发,是编写新的.cpp和.h文件,从而得到新的材料吗?”

  4. 四个火 Reply

    你好,我想请问一下,我用opensees建单墩模型,轴压比取的0.1,但是滞回曲线捏拢的很厉害,有什么办法能解决呢?

    • CJD Post authorReply

      用了什么模型,什么加载条件?常规纤维模型好像不容易做出显著的捏拢哦!

      • 四个火 Reply

        一个简单的单墩,在上面施加力,单元用的nonlinearBeamColumn
        这里是加载条件:
        pattern Plain 1 Linear {
        for {set i 1} {$i<=$N_PierNode} {incr i 1} {
        load $PierNode($i) 0.0 0.0 [expr -$PierNodalMass($i)*$g] 0.0 0.0 0.0 ;
        }
        load $DeckNode 0.0 0.0 [expr -$DeckMass*$g*$xs] 0.0 0.0 0.0;
        }
        source gravity.tcl
        puts "Gravity analysis is done"

        #pushover
        puts "set pushover analysis"
        loadConst -time 0.0;
        pattern Plain 2 Linear {
        load $DeckNode 1 0.0 0.0 0.0 0.0 0.0;
        }
        constraints Plain
        numberer Plain
        system BandGeneral
        test EnergyIncr 1.0e-6 400
        algorithm KrylovNewton
        integrator DisplacementControl $DeckNode 1 1.0e-3
        analysis Static

        for {set i 1} {$i<=400} {incr i 10} {
        set ok 0
        set step 1
        while {$ok == 0 & $step <= $i} {
        set ok [analyze 1]
        set step [expr $step+1]
        }
        integrator DisplacementControl $DeckNode 1 -1.0e-3
        set ok 0
        set step 1
        while {$ok == 0 & $step <= [expr 2*$i]} {
        set ok [analyze 1]
        set step [expr $step+1]
        }
        integrator DisplacementControl $DeckNode 1 1.0e-3
        set ok 0
        set step 1
        while {$ok == 0 & $step <= [expr $i]} {
        set ok [analyze 1]
        set step [expr $step+1]
        }
        }

        puts "step $step"
        puts "pushover analysis is done"

    • 王梓杨 Reply

      请问解决了吗?我也遇到这个问题了 方便交流一下吗

  5. Reply

    崔博士,您好!我跑opensees出现了下面的问题:
    Gravity
    goGravity
    WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 10
    WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
    StaticAnalysis::analyze() – the Algorithm failed at iteration: 0 with domain at load factor 0.1
    OpenSees > analyze failed, returned: -3 error flag
    Gravity Done.
    Recorder
    Analysis
    Rayleigh
    solve
    UMD2FA: matrix is singular. Only 972 pivots found.
    WARNING UmfpackGenLinSolver::solve(void)- 4 returned in factorization UMD2FA()
    WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
    DirectIntegrationAnalysis::analyze() – the Algorithm failed at time 0.02
    OpenSees > analyze failed, returned: -3 error flag
    Ground Motion Done. End Time: 0.000000
    OpenSees >
    这是为什么呢请问?

  6. Reply

    Gravity
    goGravity
    WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 10
    WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
    StaticAnalysis::analyze() – the Algorithm failed at iteration: 0 with domain at load factor 0.1
    OpenSees > analyze failed, returned: -3 error flag
    Gravity Done.
    Recorder
    Analysis
    Rayleigh
    solve
    UMD2FA: matrix is singular. Only 972 pivots found.
    WARNING UmfpackGenLinSolver::solve(void)- 4 returned in factorization UMD2FA()
    WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
    DirectIntegrationAnalysis::analyze() – the Algorithm failed at time 0.02
    OpenSees > analyze failed, returned: -3 error flag
    Ground Motion Done. End Time: 0.000000
    OpenSees >

  7. 张静思 Reply

    崔博士您好,有一个opensees的问题想请教您。
    做了一个双柱桥的推覆分析2维模型可以运行,但改成三维模型运行时直接闪退,不知道是什么问题。
    希望能得到您的解答,感谢!!

  8. Mr YUE Reply

    崔博士你好。想请教建模的问题。opensees桩土先建立PySimple1、TzSimple1、QzSimple1材料,然后通过零长单元建模,但是提交后总是到这个地方就中断了,无法进行下去。分析问题应该在桩土采用PySimple1、TzSimple1、QzSimple1材料这个地方,但是都是按照软件给的参数进行设置的,又无法看出问题。请您指教一下

    • CJD Post authorReply

      具体也不太清楚是什么问题。可以尝试把模型简化,看问题是否还存在。然后逐步查。另外有时候如果是一些参数或者算法问题,可以根据提示看看源代码,也许会找到一些线索。

  9. ling Reply

    崔博士好,我有一个问题想问问您,结构的pushover曲线一直往上走不屈服,该怎么调整呢,和选用的钢筋或者混凝土材料的opensees标号有关吗,谢谢

    • CJD Post authorReply

      可能的情况很多,这个得具体检查才知道,先检查是不是材料非线性参数设置有问题。

      • ling Reply

        谢谢老师,我之前说的不太对,就是其实是有屈服点的但是后面那个曲线一直往上跑,我把推的位移增大之后,出来的曲线往上升的更明显了,您能告诉我是怎么回事吗

  10. ZY Reply

    崔博士您好,请问FSAM 材料可用于三维壳单元建模吗?

  11. qi Reply

    进行截面分析 弯矩曲率分析时 输出纤维应力应变有数据 ,但输出的是弯矩
    recorder Element -file Element1.out -time -ele 1 2 Section 1 fiber 5.6 4.8 strainstress

  12. KingsLey Reply

    崔博士,想请教一下,如何用Opensees做一个简化的核电安全壳模型?可以用什么软件把模型建立,然后再导入Opensees进行计算呢?恳请崔博士指导一下!

  13. tf朱 Reply

    崔博士,您好,请问在opensees中可以直接提取动力分析时的截面弯曲-曲率关系吗。

  14. Zhao Tianrun Reply

    崔博士您好!
    冒昧请教您一个问题。我用series命令定义双折线弹性与塑性材料结合的本构:
    uniaxialMaterial Steel01 1 1.95e6 2e11 0.001
    uniaxialMaterial ElasticBilin 2 2500 10500 150
    uniaxialMaterial Series 3 1 2
    在进行pushover分析时可以得到相应的骨架曲线,但是在进行地震时程分析时,结果显示材料超过屈服力与屈服位移时仍处于线性状态,即滞回曲线为一条双折线,无塑性滞回现象。请问为什么时程分析时无法体现塑性特性呢?

      • Zhao Tianrun Reply

        我建立了一个简单的2维模型,采用element zeroLength 并赋予了X向的材料本构,pushover分析得到了我想要的骨架曲线,但是一旦进行时程分析就只有弹性材料发挥作用,底部节点剪力与顶部位移关系的曲线为双折线且没有屈服段;
        另外采用普通的弹塑性材料,如steel01进行时程分析则可以得到正常的滞回曲线。是series组合材料无法有效进行时程分析吗?

        • CJD Post authorReply

          我具体也没用过series,但按道理应该不会才对。不过从理论上来说,串联材料的确比并联材料相对复杂一点。因为从材料本构状态确定层面上来说,是应变、曲率驱动的。并联的话,由于里面各个材料应变是一致的,理论上就和单个材料的状态确定一样,这个相对较好操作,给同样一个应变增量,起码可以确定内部所有材料的力增量,进而获得总体力增量,再迭代。但对于串联材料来说,应该处理上比较麻烦,涉及到应变如何分配,尤其是如果串联很多个材料的话,理论上会很麻烦。你可以试试这样,分两步走(1)不要测试框架了,直接做一根truss构件,材料series,单自由度往复加载及时程分析,看看是不是没办法滞回。这个是直接测试是不是因为太复杂,series就不存在滞回。(2)第二个,最简单的方法,就是看series类的源代码,看看里面状态确定时候是如何处理,立刻就知道了。如果你会看源代码,那第一步可以省了。

          • Zhao Tianrun

            好的。感谢崔博士解答! 😀

  15. RW Reply

    崔博士,您好!请问OpenSees模型如果出现拟静力收敛而动力分析(PGA很小也不收敛,只能计算1%左右)不收敛的情况应该检查哪一部分呢?

    • CJD Post authorReply

      看提示是什么,看是哪部分不收敛,最好能通过提示,翻翻源代码找原因。

  16. Reply

    请问一下输出纤维截面中的单根纤维的应力应变没数据会是有啥么问题吗

  17. 郑跃 Reply

    崔博士您好,请问用OpenSees进行截面分析时,用哪个命令记录截面中单根纤维的应力、应变?

    • CJD Post authorReply

      看看 Element Recorder命令,可以记录指定局部坐标的纤维响应

    • CJD Post authorReply

      比如:recorder Element -file Element1.out -time -ele 1 3 section 1 fiber 0.10 0.10 stressStrain

  18. 豆沙包 Reply

    想问博士一个问题。我想提取出预应力钢筋中的预应力,应该用recorder中的哪一个呀?axialfroce么?

  19. YP Deng Reply

    崔博士您好,冒昧请教您一个问题,Opensees里如何批量输出弧形纤维的应力呢?

  20. hhhhhapo Reply

    想问老师一个问题。
    我在模拟一个节段梁,用什么material和element模拟节段中间的开合,不知道用什么。

    • CJD Post authorReply

      不知道你是不是要 gap 和 hook,看看这两类单元或材料的特性。

  21. 进击的人类 Reply

    崔博士您好,我有一个问题想请教您,pinching4材料的rDisp,rForce,uForce,K1,K2,K3,K4,Klim是怎么确定的呢,我查了很多资料都是说根据实验数据确定,但都没有说用什么方法。

    • CJD Post authorReply

      一般是根据实验拟合了,具体可以参考网站的 BMMF 二折线拟合程序:http://www.jdcui.com/?p=12043 及 Opensees Mat Test 程序:http://www.jdcui.com/?p=2571

  22. RW Reply

    崔博士您好,想请教您一下,如果想在opensees中定义一种只能承受拉力的材料应该如何考虑呢?

  23. fengmengyao Reply

    崔博士您好,我想请教您一个问题,关于ReinforcingSteel材料本构的Cf、Cd、α等参数的取值与设置该怎么取呢?

  24. YangXF Reply

    崔博士,你好。利用opensees做pushdown分析,梁柱采用纤维模型,模拟空间框架节点,类似与Joint2D单元,请问这个节点该怎么做?

  25. 萨卡斯基 Reply

    崔老师,您好。
    我想请教一下,在做桥墩的延性分析时,需要检测受压区核心混凝土的应变,这个应变检测位置具体应该检测那里,或是有没有具体受压区核心混凝土的应变的检测规定,请您给我解惑。
    非常感谢。

    • CJD Post authorReply

      直接测量核心区的混凝土应变可能不行了,我的感觉是,能否测量边缘钢筋的应变,间接评估混凝土的应变。

    • CJD Post authorReply

      我刚才回复的是构件试验。你是做试验吧。还是说数值模拟?

      • 萨卡斯基 Reply

        就是用OpenSees做数值模拟的时候,屈服应变一般是用受拉区外侧钢筋初次达到屈服时的应变,极限应变就是受压区核心混凝土达到极限压应变的时候,但是受压核心区我监测不同的点,得到的结果也不一样,所以很困惑。

      • 萨卡斯基 Reply

        我使用OpenSees做数值模拟,构件的屈服应变我取的是受拉区最外侧钢筋初次屈服时的应变,极限应变我取的是受压区核心混凝土达到极限压应变的值,但是受压区核心混凝土的监测点不同,算出的值也不同,所以很困惑,到底要取哪个区域的值算合适。

        • CJD Post authorReply

          那取在不同点,应变直肯定不同了。而且就算同一个点,单元划分不同等,计算的结果有差异。这些问题你具体分析吧。

          • 萨卡斯基

            非常感谢您的答复。

  26. jianzhang Reply

    麻烦请教下,对于复合材料的shell结构,opensees有没有较好的定义材料属性和单元的方法,本来想采用梁单元进行简单模拟,但是结构截面两个方向的材料属性(关于弹性对称轴)并不一致,如下图所示结果

  27. ChengQian Reply

    崔博士,您好!用OpenSees对一个超高层结构进行弹塑性时程分析,怎样得到各构件的耗能情况?求指导,谢谢

  28. WHY Reply

    崔博士 想问一下想给selfcentering定义初始预紧力,使用InitStressMaterial为何计算无变化。新手使用,还希望得到大神讲解

    • Li Yan Reply

      崔博士您好 请问OpenSees可以输出结构的耗能情况吗?

  29. Jiang linwei Reply

    崔博士,请问一下在OPENSEES如果想定义C50混凝土的数据应该咋办?

  30. NJUST李 Reply

    好网站,挺有帮助的。应该早一点知道的。很多免费软件,非常方便,用这里的软件算等效屈服弯矩和延性一下子就算出来了,比导入cad求面域方便的多,适合做学位论文时大规模处理数据使用。

    • CJD Post authorReply

      我也是之前用到就做出来,给大家,猜想应该很多学生做论文会用到。谢谢支持。

Leave a Reply

Your email address will not be published.

This site uses Akismet to reduce spam. Learn how your comment data is processed.