ansys抗震计算命令流:ANSYS地震波命令流求助

dim,d,array,1,1000,0这个最好不要用一个字母建立数组名字再说了，也要看你求解的报错是什么才好下定论是那一块出问题

在ansys中输入 help,time 然后就会跳出help来告诉你所用的答案。

单或二层框架弹性分析，需要考虑四种负载恒载，活荷载，风荷载和起重机负载 1，几何（BEAM3和beam54）成立，定义必需的元素表中，包括杆端力和最大应力和最小应力。 然后保存到数据库中。负载的四个标准值（未乘以偏系数），和四个分别保存为施加载荷步文件。 2，采用溶液 - >从LS文件，解决了四种负载 3，荷载组合，订单流程如下 / POST1 LCDEF，1,1 BR> LCDEF，2,2 LCDEF，3,3 LCDEF，4,4！定义四个条件分别计算 lcfact 4负载下，1,1.2 lcfact，2,1.4 lcfact，3,1.19 lcfact，4,1.4！指定每个组合系数 LCASE，1的条件！读条件1，数据库= 1 sumtype，首席！指定对象以及操作 LCOPER，加2！荷载组合，数据库=数据库+ 2 LCOPER，加，4！荷载组合，数据库=数据库+ 4 LCOPER，lprin！计算线性主应力 LCWRITE，11！该数据库并将结果写入11条件的恒载+活载+吊车荷载结果，点击看LCASE，1 lcfact，2,1.19 lcfact，4 1.19！改变组合系数 sumtype，首席 LCOPER，添加2 LCOPER，添加，3 LCOPER，添加4 LCOPER，lprin LCWRITE，12！ 12个结果数据库写入条件是恒载+活载+吊车荷载+风荷载结果，点击看！ ......其他的荷载组合 ！使用LCASE之后，正的条件可以被转移到n，并查看其变形和内力 ！您可以使用下面的命令来获得流动条件11和12，13，99，以较高者为准，然后查看最大应力 LCASE，11 LCASE，分钟，12 LCASE，分钟， 13 LCWRITE，98 LCASE 98 ！查看98条件下的应力分布...... LCASE，11 LCASE，最大，12 LCASE，最大，13 LCWRITE，99 LCASE 99 ！检查应力分布状况99 ...... 使用下面的命令来定义并阅读一些载荷工况 LCDEF_ LCDEF生成负载情况下的结果文件的结果的一列，LCNO，LSTEP，SBSTEP，KIMG LCNO指针的任意数（1-99），被分配到LSTEP，SBSTEP和FILE命令中指定的负载条件。缺省值被加到一个前一个。 LLSTEP被定义作为负载的情况下的负载的步骤号。默认值为1 SBSTEP副荷载步数。默认值是负载阶跃负载阶跃的一个孩子。 KIMG只为1- 0-多个分析用虚部注意的复杂的分析的实数部分通过建立一个链接到文件指针的结果，以产生一个结果的列表载荷工况。该指针（LCNO）可以用在LCASE或LCOPER命令来读取负载情况下的数据到数据库中。 LCDEF，ERASE删除所有载荷工况指针（和所有负载情况下的文件）。随着LCDEF，LCNO，ERASE删除指定的荷载工况指针LCNO（以及相应的文档）。当选项为ERASE，所有指针都将被删除，但仅限于默认的文件扩展名（LCWRITE）将被删除。写LCDEF，STAT看到所有选定的工况（LCSEL）状态，写LCDEF，STAT，ALL看所有负载情况下的状态。 STAT命令可以用来列出所有荷载工况。 LCFILE看看如何建立一个指向荷载工况文件（书面LCWRITE）导致列。从文件单元谐波负载的情况下的结果被读入存储在所述零位置的数据。 加法，ANSYS是相当困难的研究深入，对我来说，是没有办法学一两年，达到熟练程度

我看到的一个例子是进入塑性后（屈服后）用21个（自定）应力-应变点拟合，或者你就直接输入屈服强度～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～MP,EX,2,115e3 MP,PRXY,2,0.3TB, MISO,2,1,21,0TBPT,,0.008260869565,950TBPT,,0.01030599278,966.25TBPT,,0.01148279762,982.5TBPT,,0.01319060099,998.75TBPT,,0.01568580449,1015TBPT,,0.01934123137,1031.25TBPT,,0.02469575862,1047.5TBPT,,0.03252378039,1063.75TBPT,,0.04393191531,1080TBPT,,0.06049295633,1096.25TBPT,,0.08443048782,1112.5TBPT,,0.1188721142,1128.75TBPT,,0.1681951823,1145TBPT,,0.2384966466,1161.25TBPT,,0.3382288502,1177.5TBPT,,0.4790561384,1193.75TBPT,,0.677004215,1210TBPT,,0.9539960686,1226.25TBPT,,1.339896413,1242.5TBPT,,1.875222619,1258.75TBPT,,2.614726016,1275 ～～～～～～～～～～～～～～～～～～～～～～～你又加上温度非线性，屈服后的非线性行为下弹模不是固定的吧，你这样也能算，定量上差一些吧

!载荷步为n，假设等于100，只需要修改参数n就可以跑后面的命令流!每一步的最大位移储存在 UMAX+载荷步 这个参数里面!每一步的最大应力储存在 SMAX+载荷步 这个参数里面n=100set,first!最大总位移，若是某个方向的位移可以修改下面命令中的sum为x,y,z等方向nsort,u,sum,0,0,all GET,NODEMAX,SORT,0,IMAXGET,UMAX1,node,NODEMAX,u,sum!MISES应力，若是某个方向的应力可以修改下面命令中的eqv为x,y,z,xy,yz,zx,int等nsort,s,eqv,0,0,all GET,NODEMAX,SORT,0,IMAXGET,SMAX1,node,NODEMAX,s,eqvdo,ii,2,n,1set,next!最大总位移，若是某个方向的位移可以修改下面命令中的sum为x,y,z等方向nsort,u,sum,0,0,all GET,NODEMAX,SORT,0,IMAXGET,UMAX%ii%,node,NODEMAX,u,sum!MISES应力，若是某个方向的应力可以修改下面命令中的eqv为x,y,z,xy,yz,zx,int等nsort,s,eqv,0,0,all GET,NODEMAX,SORT,0,IMAXGET,SMAX%ii%,node,NODEMAX,s,eqvenddostat,all

bed_l=20 ! 长度mbed_w=6 ! 宽度bed_t1=4e-2 ! 厚度1bed_t2=6e-2bed_t3=18e-2bed_t4=15e-2bed_t5=15e-2bed_t6=20e-2bed_t7=2.5ex_1=1500e6 ! 弹性模量paex_2=1300e6ex_3=1600e6ex_4=800e6ex_5=600e6ex_6=90e6ex_7=90e6PRXY_1=0.35PRXY_2=0.35PRXY_3=0.25PRXY_4=0.35PRXY_5=0.35PRXY_6=0.4PRXY_7=0.4finish/prep7Et,1,SOLID191KEYOPT, 1, 8, 1mp,ex,1,ex_1mp,ex,2,ex_2mp,ex,3,ex_3mp,ex,4,ex_4mp,ex,5,ex_5mp,ex,6,ex_6mp,ex,7,ex_7mp,prxy,1,PRXY_1mp,prxy,2,PRXY_2mp,prxy,3,PRXY_3mp,prxy,4,PRXY_4mp,prxy,5,PRXY_5mp,prxy,6,PRXY_6mp,prxy,7,PRXY_7r,1,7RMORE,RMORE,1,,bed_t1,2,,bed_t2RMORE,3,,bed_t3,4,,bed_t4RMORE,5,,bed_t5,6,,bed_t6RMORE,7,,bed_t7,BLOCK, 0, bed_l, 0, bed_w,0,bed_t1+bed_t2+bed_t3+bed_t4+bed_t5+bed_t6+bed_t7LSEL, s, line, , 9, 12LESIZE,all, , , 1esize,0.7vMESH,1allsel/ESHAPE,1.0 ! 这个命令来看单元的层次，后面改成0就是没有单元层了，与下面的一个命令一起用就可以看到了/REPLOTasel,s,area,,2,6nsla,s,1d,all,all,0allsel/VIEW, 1, -0.202926483459 , 0.833744595865 , -0.513508316560/ANG, 1, -172.396473039 /REPLO PRES_bed=100e3/400e-4 ! 压强，paSFA, 1, 1, pres,PRES_bed ! 出现警告不用管allselfinish/solusolvefinish/post1LAYER,7PRESOL, s, int这个是一种建模求解方式！

好费劲啊！呵呵，你参考下！关于屈曲分析的。 /solu ! 进入求解器antype,0 ! 静力求解eqslv,spar ! 选择spar求解器pstres,on ! 打开预应力选项solve ! 求解finish /solu ! 再一次求解，前一次的主要是考虑预应力效应，这个屈曲分析中比较常见antype,1 ! 屈曲分析bucopt,lanb,1,0,0 ! 指定屈曲分析的分析选项，后面的三个数字是默认值mxpand,1,0,0,1,0.001 ! 模态扩展分析扩展分析，并制定计算单元和反作用力求解信息，后面的0.001是默认值solve finish /post1 ! 进入后处理器set,first ! 选择第一载荷步set,list ! 列出第一载荷步的求解信息结果pldisp,1 ! 显示模态结果，也是位移，包括单元和轮廓线finish /config,nres,100 ! 允许存在结果文件中的最大号码（一般指的是载荷步）/prep7 ! 进入预处理器tb,biso,1,1,2 ! 指定属性数据表，用于同向双线型分析，塑性的分布情况等，指定1好材料属性，再1指定温度号码tbtemp,0 ! 为数据表定义温度tbdata,,2.0e8,0,,,, ! 为刚才定义的属性数据表定义数据，从默认开始，这个要结合所选择单元才能看明白，第一个省略说明是默认选项一次定义6个数据upgeom,0.01,1,1,'qiangdu','rst'," ! 第一个是比例因子，向后面定义模型，小于1，说明缩小，两个1是载荷步，和子步；后面的是新生成的文件的名字和后缀将分析所得的位移加到有限元模型的 节点上并更新有限元模型的几何形状finish /solu antype,0 ! 稳态分析nlgeom,1 ! 指出包含大变形outres,all,all ! 输出所有的结点和单元的求解信息arclen,1,0 ! 激活非线性分析的弧长法，不过最大半径设为0arctrm,u,0.5,78,uy ! 弧长求解终止条件，当一次位移等于或者大于设定值时停止分析，位移设定为0.5，以78结点的位移为标准，Uy方向有效，也就是78结点爱UY方向有0.5位移是就停止分析！nsubst,200,,,1 ! 指定本载荷步为200，并把当前载荷步的时间步作为起始时间步solve finish /post26 ! 进入时间后处理器nsol,2,78,u,y,deflection ! 把78结点的在Y方向的位移结果储存在deflection标识名下，2为代号（特征号，这个求解的）rforce,3,1861,f,x,reactionf ! 把1861结点在X的结构力储存在reactionf标识下，3为代号/axlab,x,deflection ! 表明图形显示中X轴标签的显示内容，下面为Y轴的/axlab,y,reactionf xvar,2 ! 指定X轴的2代号的信息，图形，plvar,3 ! 显示3代号的信息，这里也就是在Y轴的显示图形，最多可以显示10个变量，一次finish ! 推出处理器

dim,d,array,0这个最好不要用一个字母建立数组名字再说了。

time 然后就会跳出help来告诉你所用的答案。

单或二层框架弹性分析，需要考虑四种负载恒载，活荷载，风荷载和起重机负载 1，几何（BEAM3和beam54）成立，定义必需的元素表中，包括杆端力和最大应力和最小应力。然后保存到数据库中。负载的四个标准值（未乘以偏系数），和四个分别保存为施加载荷步文件。采用溶液 - >从LS文件，解决了四种负载 3，订单流程如下/1 BR>定义四个条件分别计算 lcfact 4负载下，指定每个组合系数 LCASE，1的条件！读条件1，首席！指定对象以及操作 LCOPER，加，荷载组合，数据库=数据库+ 4 LCOPER，计算线性主应力 LCWRITE，该数据库并将结果写入11条件的恒载+活载+吊车荷载结果，改变组合系数 sumtype，添加，12个结果数据库写入条件是恒载+活载+吊车荷载+风荷载结果，点击看！......其他的荷载组合！使用LCASE之后，正的条件可以被转移到n，并查看其变形和内力！您可以使用下面的命令来获得流动条件11和12，以较高者为准，然后查看最大应力 LCASE，分钟，查看98条件下的应力分布...... LCASE，最大，检查应力分布状况99 ...... 使用下面的命令来定义并阅读一些载荷工况LCDEF_ LCDEF生成负载情况下的结果文件的结果的一列，指针的任意数（1-99），被分配到LSTEP，SBSTEP和FILE命令中指定的负载条件。缺省值被加到一个前一个。被定义作为负载的情况下的负载的步骤号。默认值为1 SBSTEP副荷载步数。默认值是负载阶跃负载阶跃的一个孩子。只为1- 0-多个分析用虚部注意的复杂的分析的实数部分通过建立一个链接到文件指针的结果，以产生一个结果的列表载荷工况。该指针（LCNO）可以用在LCASE或LCOPER命令来读取负载情况下的数据到数据库中。ERASE删除所有载荷工况指针（和所有负载情况下的文件）。ERASE删除指定的荷载工况指针LCNO（以及相应的文档）。当选项为ERASE，所有指针都将被删除，但仅限于默认的文件扩展名（LCWRITE）将被删除。STAT看到所有选定的工况（LCSEL）状态，写LCDEF，STAT，ALL看所有负载情况下的状态。STAT命令可以用来列出所有荷载工况。

载荷步为n!假设等于100，只需要修改参数n就可以跑后面的命令流，每一步的最大位移储存在 UMAX+载荷步 这个参数里面!每一步的最大应力储存在 SMAX+载荷步 这个参数里面n=100set!最大总位移!若是某个方向的位移可以修改下面命令中的sum为x，yz,ii,n,最大总位移!若是某个方向的位移可以修改下面命令中的sum为x，z等方向nsort,node,u,sum,MISES应力!若是某个方向的应力可以修改下面命令中的eqv为x，y,z,xy,yz,zx,int等nsort,s,node,s,eqvenddostat,all,

!载荷步为n，假设等于100，只需要修改参数n就可以跑后面的命令流!每一步的最大位移储存在 UMAX+载荷步 这个参数里面!每一步的最大应力储存在 SMAX+载荷步 这个参数里面n=100set,first!最大总位移，若是某个方向的位移可以修改下面命令中的sum为x,y,z等方向nsort,u,sum,0,0,all GET,NODEMAX,SORT,0,IMAXGET,UMAX1,node,NODEMAX,u,sum!MISES应力，若是某个方向的应力可以修改下面命令中的eqv为x,y,z,xy,yz,zx,int等nsort,s,eqv,0,0,all GET,NODEMAX,SORT,0,IMAXGET,SMAX1,node,NODEMAX,s,eqvdo,ii,2,n,1set,next!最大总位移，若是某个方向的位移可以修改下面命令中的sum为x,y,z等方向nsort,u,sum,0,0,all GET,NODEMAX,SORT,0,IMAXGET,UMAX%ii%,node,NODEMAX,u,sum!MISES应力，若是某个方向的应力可以修改下面命令中的eqv为x,y,z,xy,yz,zx,int等nsort,s,eqv,0,0,all GET,NODEMAX,SORT,0,IMAXGET,SMAX%ii%,node,NODEMAX,s,eqvenddostat,all

bed_l=20!长度mbed_w=6!宽度bed_t1=4e-2!厚度1bed_t2=6e-2bed_t3=18e-2bed_t4=15e-2bed_t5=15e-2bed_t6=20e-2bed_t7=2.5ex_1=1500e6!弹性模量paex_2=1300e6ex_3=1600e6ex_4=800e6ex_5=600e6ex_6=90e6ex_7=90e6PRXY_1=0.35PRXY_2=0.35PRXY_3=0.25PRXY_4=0.35PRXY_5=0.35PRXY_6=0.4PRXY_7=0.4finish/prep7Et,ex,ex_7mp,PRXY_7r,bed_t1,bed_t7,BLOCK,bed_l,bed_w,bed_t1+bed_t2+bed_t3+bed_t4+bed_t5+bed_t6+bed_t7LSEL,line,这个命令来看单元的层次,后面改成0就是没有单元层了,与下面的一个命令一起用就可以看到了/,area,VIEW,ANG,-172.396473039 /,REPLO PRES_bed=100e3/,压强;pres,出现警告不用管allselfinish/,solusolvefinish/,7PRESOL;s!int这个是一种建模求解方式。

好费劲啊！你参考下！/进入求解器antype,静力求解eqslv,打开预应力选项solve!再一次求解，前一次的主要是考虑预应力效应，这个屈曲分析中比较常见antype,屈曲分析bucopt,指定屈曲分析的分析选项，后面的三个数字是默认值mxpand,模态扩展分析扩展分析，并制定计算单元和反作用力求解信息，后面的0.001是默认值solve finish /进入后处理器set,选择第一载荷步set,list!列出第一载荷步的求解信息结果pldisp,显示模态结果，也是位移，包括单元和轮廓线finish /config,允许存在结果文件中的最大号码（一般指的是载荷步）/进入预处理器tb,指定属性数据表，用于同向双线型分析，塑性的分布情况等，指定1好材料属性，再1指定温度号码tbtemp,为数据表定义温度tbdata,为刚才定义的属性数据表定义数据,从默认开始,这个要结合所选择单元才能看明白,第一个省略说明是默认选项一次定义6个数据upgeom,',",第一个是比例因子;向后面定义模型;小于1,说明缩小;两个1是载荷步;和子步,后面的是新生成的文件的名字和后缀将分析所得的位移加到有限元模型的 节点上并更新有限元模型的几何形状finish /稳态分析nlgeom，指出包含大变形outres，输出所有的结点和单元的求解信息arclen;激活非线性分析的弧长法,不过最大半径设为0arctrm!u,uy,弧长求解终止条件,当一次位移等于或者大于设定值时停止分析!位移设定为0.5，以78结点的位移为标准,Uy方向有效,也就是78结点爱UY方向有0.5位移是就停止分析,指定本载荷步为200，并把当前载荷步的时间步作为起始时间步solve finish /，进入时间后处理器nsol，u,把78结点的在Y方向的位移结果储存在deflection标识名下!2为代号（特征号，这个求解的）rforce;f,把1861结点在X的结构力储存在reactionf标识下,3为代号/!x，表明图形显示中X轴标签的显示内容,下面为Y轴的/,y,指定X轴的2代号的信息;图形,