opensees命令解释.docxVIP

1、定义梁柱单元局部坐标轴的命令流为：?geomTransf?Linear?$transfTag?$vecxzX?$vecxzY?$vecxzZ?其中，$transfTag?代表局部坐标轴矢量的编号，$vecxzX?$vecxzY?$vecxzZ?表示局部坐标轴的方向矢量值。?2、OPENSEES?的刚性隔板假定命令流格式为：?rigidDiaphragm?$perpDirn?$masterNodeTag?$slaveNodeTag1?$slaveNodeTag2?...其中，$perpDirn?表示刚性隔板的方法，如实例中楼板的刚性隔板的平移方向为U1（X?方向）与U2（Y?方向），即1-2?平面，该值应为3。$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?是一致的，表示一个单元有自已特定的坐标轴向量，为了编程的方便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?是最简单的混凝土本构，注意数值是负数，即表示受压段。该本构没有受拉段，即受拉强度为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?表示每个纤维的贡献面积；$?Mat?表示每个纤维使用的非线性材料本构的编号。?注意：纤维的坐标与材料切线模量可以组装成截面的刚度，而纤维的坐标与材料的应力可以组装成截面的内力（抗力），那么每个纤维的应变可以通过截面的变形与坐标求出。采用纤维截面的单元，即为基于平截面假定。截面变形求解应变是基于平截面假定的。?6、采用的单元为非线性梁柱单元，即基于柔度法的纤维单元（Nonlinear?BeamColumn?Element?or?Force?Beam?Column?Element），需要输入命令流如下：?element?nonlinearBeamColumn?$eleTag?$iNode?$jNode?$numIntgrPts?$secTag?$transfTag?其中，$eleTag?为单元编号；$iNode?为开始结点；$jNode?为结束结点；$numIntgrPts?为积分点数量；$secTag?为纤维截面编号，$transfTag?为局部坐标轴编号。积分点数量，也就纤维单元的计算截面数量，纤维单元的刚度与抗力是由截面刚度与抗力沿杆件长度积分所得，显然，不能将全部截面积分，只能采用离散积分，OPENSEES?默认的积分方法是高斯-洛贝塔积分（Guass-Lobotto），各阶积分点分布及权函数如下图所示。（OPENSEES?后期版本可提供其它积分方法：Legendre,?Radau,NewtonCotes,?Trapezoidal）。梁柱单元的局部坐标的计算与编号在上章节已描述。??7、截面变形记录：本实例采用纤维单元，即可输出每段积分点处截面的变形，记录截?面变形的命令如下：?recorder?Element?-file?$File?-time?-eleRange?1?6?section?$SecPos?def