最新opensees总结

1、精品文档1、定义梁柱单元局部坐标轴的命令流为：geomTransf Linear $transfTag $vecxzX $vecxzY $vecxzZ其中，$transfTag代表局部坐标轴矢量的编号，$vecxzX $vecxzY $vecxzZ表示局部坐标轴的方向矢量值。2、OPENSEE的刚性隔板假定命令流格式为：rigidDiaphragm $perpDirn $masterNodeTag $slaveNodeTag1 $slaveNodeTag2 .其中，$perpDirn表示刚性隔板的方法，如实例中楼板的刚性隔板的平移方向为U1（X方向）与U2（丫方向），即1-2平面，该值应为 3

2、。$masterNodeTag为主结点，$slaveNodeTag1为从结点。主结点一般为 刚性隔板刚心。实例中：rigidDiaphragm 3 35 2，表示刚性隔板平动方向为1-2平面，刚心主节点为 35点，2号结点为从结点。3、弹性梁柱单元的命令流：element elasticBeamColumn $eleTag $iNode $jNode $A $E $G $J $Iy $Iz $transfTag需要提供截面的截面积 A、截面Y轴惯性矩Iy，截面Z轴惯性矩Iz，截面扭转矩，截面材料的弹性模量 E及 剪切模量G。其中：$transfTag与$eleTag是一致的，表示一个单元有自已

3、特定的坐标轴向量，为了编程 的方便。陈：例题三4、非线性材料模型的定义（1）uniaxialMaterial Steel01 1 335 200000 0.00001表示，钢筋的屈服强度为 335MPa，弹性模量为200000MPa，硬化系数为0.00001，即屈服平台基本上为水 平段。将混凝土材料本构 C40改为非线性混凝土本构【Concrete01】，命令流如下：（2）uniaxialMaterial Concrete01 2 -26.8 -0.002 -10 -0.0033材料参数意见参考图所示。注意：混凝土本构 Concrete01是最简单的混凝土本构，注意数值是负数， 即表示受压段。

4、该本构没有受拉段，即受拉强度为0,表示结构一分析即进入弹塑性。5、 采用纤维单元，需要定义纤维截面，纤维截面的定义如下面代码所示section Fiber 1 fiber -1.125E+002 -2.700E+002 4.500E+003 2 fiber 1.150E+002 -2.650E+002 4.900E+002 1以上命令流表示，纤维截面编号为1，内部为子命令流，表示每一个纤维的信息，每一个纤维的定义格式如下：fiber $Y $Z $Area $Mat命令中，$Y表示每个纤维的截面 Y坐标（截面中心为原点0）; $Z表示每个纤维的截面Z坐标（截面中 心为原点0）； $ Area表

5、示每个纤维的贡献面积；$ Mat表示每个纤维使用的非线性材料本构的编号。注意：纤维的坐标与材料切线模量可以组装成截面的刚度，而纤维的坐标与材料的应力可以组装成截面的 内力（抗力），那么每个纤维的应变可以通过截面的变形与坐标求出。采用纤维截面的单元，即为基于平 截面假定。截面变形求解应变是基于平截面假定的。6 采用的单元为非线性梁柱单元，即基于柔度法的纤维单元（Nonlinear BeamColumn Element or ForceBeam Column Element），需要输入命令流如下：element nonlinearBeamColumn $eleTag $iNode $jNode $

6、numlntgrPts $secTag $transfTag 其中，$eleTag为单元编号；$iNode为开始结点；$jNode为结束结点；$numIntgrPts为积分点数量； $secTag为纤维截面编号，$transfTag为局部坐标轴编号。积分点数量，也就纤维单元的计算截面数量，纤维单元的刚度与抗力是由截面刚度与抗力沿杆件长度积分所得，显然，不能将全部截面积分，只能采用精品文档离散积分，OPENSEES默认的积分方法是高斯-洛贝塔积分(Guass-Lobotto)，各阶积分点分布及权函数 如下图所示。(OPENSEES后期版本可提供其它积分方法：Legendre, Radau,New

7、tonCotes, Trapezoidal )。梁柱单元的局部坐标的计算与编号在上章节已描述。9=0.0(0i= 1 / 3&=0.5<0：=4 3=1.0to-=l ;，3*51=0.0&=0.2763<4=0.7237£5=1.001=1/6o；=5/6c>s=l 6*"1VVTi=0.0&=0.1726&=0.5=0.8273 i=rOC0j=0J o>=O.5444 ©3=0.7111 ©4=0.5444 ®5=0J7、截面变形记录：本实例采用纤维单元，即可输出每段积分点处截面的变

8、形，记录截面变形的命令如下：recorder Element -file $File -time -eleRange 1 6 section $SecPosdeformation其中，$File为记录的文本文件名； -eleRange 1 6,表示记录单元为16号；$SecPos为输出积分点号，如 下图所示，开始结点为1号，如果积分点数为N，结束结点的截面编号为 N，deformation表示记录内容为 截面变形。截面变形输出结果中，输出三个数值分别是轴向应变，绕Z轴弯曲曲率，绕Y轴弯曲曲率。精品文档8、均布荷载的定义。 本实例采用均布线荷载，即施加于单元的均布荷载，OPENSEES的线性荷载

9、有两种，一种是均布荷载，另一种是线上点荷载。以均布荷载为例，命令流格式如下：eleLoad -ele $eleTag1 v$eleTag2 .-type -beamUniform $Wy $Wz v$Wx>其中，$eleTag1为施加荷载的单元编号，-type -beamUniform 表示均布线荷载，$Wy $Wz v$Wx>为三个轴方向的均布荷载值。$Wy表示截面局部坐标轴垂直方向的均布荷载，$Wz表示截面局部坐标轴方向的均布荷载，$Wx表示沿单无长度方向的均布荷载，如下图所示。梁荷载向量示意图陈：例题四9、截面组装。纤维截面是依据纤维截面(Section Fiber)命令流形

10、成，命令流如下所示:#NC400X400 section Fiber 1 fiber -1.667E+002 -1.667E+002 4.444E+003 2 fiber 1.650E+002 0.000E+000 4.900E+002 1可见，纤维截面的编号为1 (纤维截面算1个截面本构)，对应该截面的抗剪与抗扭本构编号分别为201, 精品文档301 及 401，通过截面组装的方式将这四个截面本构组装在一起，命令流如下： section Aggregator 1001 201 Vy 301 Vz 401 T -section 1其中，1001为组装后的截面编号，201 Vy代表Vy （Y方向

11、抗剪）采用单轴本构 201号；301 Vz代表Vz （Z方向抗剪）采用单轴本构 301号；401 T代表T （抗扭）采用单轴本构 401号；-section 1表示参与组 装的纤维截面为 1 号截面；就这样将三个材料本构与一个纤维截面组装在一组，那么组装后的纤维截面编 号即改为 1001 号。那么，在纤维单元定义时，单元采用的纤维截面编号即不为1 号，为 1001 号，如下命令流所示。element nonlinearBeamColumn 1 1 2 4 1001 1上述 6 ：element nonlinearBeamColumn $eleTag $iNode $jNode $numIntg

12、rPts $secTag $transfTag其中， $eleTag 为单元编号； $iNode 为开始结点； $jNode 为结束结点； $numIntgrPts 为积分点数量 ； $secTag 为纤维截面编号， $transfTag 为局部坐标轴编号。10、Push-over分析设置。Push-over是高层建筑结构弹塑性分析方法的一种，其分析过程分两个部分，一个就是静力弹塑性分析部分，另一部分为 Push-Over 曲线转化为能力需求曲线进行抗震评估。本实例只 完成第一部分。一般情况下， Push-over 过程分重力荷载加载过程（力控制），水平位移加载过程（位移 控制）。命令流如下：

13、重力加载的命令流如下：puts "gravity"# Load Case = DEADpattern Plain 1 Linear eleLoad -ele 5 -type -beamUniform 0 -3.797E+000 0 eleLoad -ele 46 -type -beamUniform 0 -3.797E+000 0puts "analysis"constraints Plainnumberer Plainsystem BandGeneraltest EnergyIncr 1.0e-6 200algorithm Newtonintegrat

14、or LoadControl 1.000E-001analysis Staticanalyze 10与上述章节介绍的是一样的，是一个简单的力控制加载，假定总施加的重力荷载为G,总荷载步为10步,每步施加 0.1G。注意：施加的荷载工况( pattern Plain )与分析设置，一般连在一起，表明后面的分析工况采用的荷载就是上面的命令流描述的荷载。所以力控制分析( loadControl)的分析工况所采用的荷载为上面的荷载( pattern Plain 1 )。位移控制加载的命令流如下：loadConst 0.0puts "pushover"# Load Case = PU

15、SHpattern Plain 2 Linear load 6 1.350E+004 0.000E+000 0.000E+000 0.000E+000 0.000E+000 0.000E+000 load 20 4.500E+003 0.000E+000 0.000E+000 0.000E+000 0.000E+000 0.000E+000 精品文档 puts "analysis" constraints Plain numberer Plain system BandGeneral test EnergyIncr 1.0e-6 200 algorithm Newton integrator DisplacementControl 8 1 2.000E-001 analysis Static analyze 100 上述命令流中