ANSYS应用——m法计算桩基时模拟桩与土作用的命令流

M-Method Calculation of Pile Foundation ExampleBackground:The M-method is a commonly used analytical approach for estimating the capacity of pile foundations. It considers both the axial and lateral loads acting on the pile, as well as the soil-pile interaction. This method is based on the assumption that the pile-soil interaction can be represented by a spring-dashpot system.Example:Let's consider a pile foundation installed in a soil deposit with known geotechnical properties. The pile has a diameter of 1.5 meters and a length of 20 meters. The axial load acting on the pile is 1000 kN, and the lateral load is 200 kN at a depth of 5 meters.Calculation Steps:Soil Properties:Determine the soil's modulus of elasticity (E) and shear modulus (G).Determine the soil's cohesion (c) and friction angle (φ).Pile Properties:Determine the pile's cross-sectional area (A) and moment of inertia (I).Determine the pile's modulus of elasticity (Ep).Spring-Dashpot Model:Calculate the axial stiffness (Kz) and lateral stiffness (Ky) of the pile-soil system using the M-method formulas.Determine the damping coefficients (Cz and Cy) if required for dynamic analysis.Load-Displacement Analysis:Apply the axial and lateral loads to the pile.Solve the governing equations of motion to determine the displacement and internal forces in the pile.Capacity Estimation:Evaluate the pile's ultimate axial and lateral capacities based on the displacement and internal force results.Compare the estimated capacities with the applied loads to assess the pile's safety.Conclusion:Using the M-method, we can estimate the axial and lateral capacities of pile foundations. This method provides a practical and efficient way to analyze pile-soil interaction and assess the pile's performance under various loading conditions.背景：M法是一种常用于估算桩基承载力的分析方法。

1问题描述桩基础是桥梁工程中广泛应用的重要基础形式之一。

如果场地浅层土的承载力低，无法满足桥梁结构对地基变形和承载力的要求时，需要考虑采用柱基础。

此次课程设计模拟了混泥土桩基（摩擦型）在竖向均布荷载作用下的反应。

具体设计资料如下：1.1柱基础假定场地的软弱土层较厚，桩端达不到坚硬土层或岩层上，桩顶的荷载主要靠桩身与土体之间的摩擦力来支承，桩尖处土层反力很小，可以忽略不计。

桩身采用C20混泥土，混泥E=3.2×1010N/m2，混泥土密度2500 KG/m3，混泥土泊松比0.167。

土抗压弹性模量C1.2土体由于桩基对周围土体的影响随着深度和影响半径的增大而逐渐减小，因此土体按照有限E=2.6×108N/m2，土体密度体积来考虑。

假设桩身周围的土体均质，土体的抗压弹性模量C1900 KG/m3，土体的泊松比0.42，桩基与周围土体的摩擦系数取0.2。

1.3荷载状况桥跨上部结构传递下来的荷载简化成竖向均布荷载，直接作用于桩基础顶部，不考虑水平力和弯矩的影响。

竖向均布荷载设计值为50×104Pa。

2单元的选择2.1桩基础混凝土桩基础，采用SOLID45单元。

SOLID45单元是八节点三维实体单元,每一个节点具有三个自由度。

单元的几何形状、结点位置和单元坐标系如图1所示。

该单元具有塑性、蠕变、膨胀、应力强化以及大变形大应变和模拟各向异性等功能,所以模型中的桩基础混凝土单元采用SOLID45实体单元。

图1 SOLID45单元2.2土体土体单元选择170，土体与桩基的接触单元选择173。

ANSYS中能用于岩土材料的模型只有DP模型。

DP模型是理想弹塑性模型，理想弹塑性即应力达到屈服极限以后，应力不再增大，但是应变会一直增大。

ANSYS中设定DP模型需要输入3个参数，粘聚力，内摩擦角，膨胀角，其中的膨胀角是用来控制体积膨胀的大小的。

在岩土工程中，一般密实的砂土和超强固结土在发生剪切的时候会出现体积膨胀，因为颗粒重新排列了；而一般的砂土或者正常固结的土体，只会发生剪缩。

ANSYS加载命令流⼀，在ANSYS⾥施加地震惯性⼒的⽅法在ANSYS⾥做地震分析时，需要对结构施加地震惯性荷载，地震惯性⼒是通过加速度的⽅式输⼊进结构的，然后与结构的质量⼀起形成动⼒计算时的惯性荷载，下⾯说⼀下在ANSYS ⾥施加地震惯性⼒的⽅法。

⾸先，将三个⽅向的地震加速度放到⼀个⽂本⽂件⾥，如accexyz.txt，在这个数据⽂件⾥共放三列数据，每列为⼀个⽅向的地震加速度值，这⾥仅给出数据⽂件中前⼏⾏的数据：-0.227109E-02 -0.209046E+00 0.467072E+01-0.413893E-02 -0.168195E+00 0.261523E+01-0.574753E-02 -0.157890E+00 0.809014E-01-0.731227E-02 -0.152996E+00 0.119975E+01-0.876865E-02 -0.138102E+00 0.130902E+01-0.101067E-01 -0.131582E+00 0.143611E+00 .......................然后，再建⼀个⽂本⽂件⽤来存放三个⽅向的地震加速度时间点，如time.txt，在这个数据⽂件⾥仅⼀列数据，对应于加速度数据⽂件⾥每⼀⾏的时间点，这⾥给出数据⽂件中前⼏⾏数据：0.100000E-010.200000E-010.300000E-010.400000E-010.500000E-010.600000E-01.......................编写如下的命令流⽂件，并命名为acce.inp*dim,ACCEXYZ,TABLE,2000,3 !01⾏*vread,ACCEXYZ(1,1),accexyz,txt,,JIK,3,2000 !02⾏(3e16.6) !03⾏*vread,ACCEXYZ(1,0),time,txt !04⾏(e16.6) !05⾏ACCEXYZ(0,1)=1 !06⾏ACCEXYZ(0,2)=2 !07⾏，同上ACCEXYZ(0,3)=3 !08⾏，同上finish/SOLUANTYPE,transbtime=0.01 !定义计算起始时间etime=15.00 !定义计算结束时间dtime=0.01 !定义计算时间步长*DO,itime,btime,etime,dtimetime,itimeAUTOTS,0NSUBST,1, , ,1KBC,1acel,ACCEXYZ(itime,1),ACCEXYZ(itime,2),ACCEXYZ(itime,3) !施加三个⽅向的地震加速度SOLVE*ENDDO最后，在命令窗⼝⾥输⼊/input,acce,inp即可对结构进⾏地震动⼒分析。

《史上最全的ANSYS命令流查询与解释》【1】*************************************************************************************对ansys主要命令的解释1， /PREP7 ! 加载前处理模块2， /CLEAR,NOSTART ! 清除已有的数据, 不读入启动文件的设置(不加载初始化文件)初始化文件是用于记录用户和系统选项设置的文本文件/CLEAR, START !清除系统中的所有数据,读入启动文件的设置/, EX10.5 ! 定义工程文件名称/TITLE, EX10.5 SOLID MODEL OF AN AXIAL BEARING ! 指定标题4， F,2,FY,-1000 ! 在2号节点上施加沿着-Y方向大小为1000N的集中力6， FINISH ! 退出模块命令7， /POST1 ! 加载后处理模块8， PLDISP,2 ! 显示结构变形图,参数“2”表示用虚线绘制出原来结构的轮廓9， ETABLE,STRS,LS,1 ! 用轴向应力SAXL的编号”LS,1”定义单元表STRSETABLE, MFORX,SMISC,1 ! 以杆单元的轴力为内容, 建立单元表MFORXETABLE, SAXL, LS, 1 ! 以杆单元的轴向应力为内容, 建立单元表SAXLETABLE, EPELAXL, LEPEL, 1 ! 以杆单元的轴向应变为内容, 建立单元表EPELAXLETABLE,STRS_ST,LS,1 !以杆件的轴向应力“LS,1”为内容定义单元表STRS_STETABLE, STRS_CO, LS,1 !以杆件的轴向应力“LS,1”定义单元表STRS_COETABLE,STRSX,S,X ! 定义X方向的应力为单元表STRSXETABLE,STRSY,S,Y ! 定义Y方向的应力为单元表STRSY*GET,STRSS_ST,ELEM,STEEL_E, ETAB, STRS_ST !从单元表STRS_ST中提取STEEL_E单元的应力结果，存入变量STRSS_ST;*GET, STRSS_CO,ELEM,COPPER_E,ETAB,STRS_CO”从单元表STRS_CO中提取COPPER_E单元的应力结果，存入变量STRSS_CO10 FINISH !退出以前的模块11, /CLEAR, START ! 清除系统中的所有数据,读入启动文件的设置12 /UNITS, SI !申明采用国际单位制14 /NUMBER, 2 !只显示编号, 不使用彩色/NUMBER, 0 ! 显示编号, 并使用彩色15 /SOLU ! 进入求解模块:定义力和位移边界条件,并求解ANTYPE, STATIC ! 申明分析类型是静力分析(STATIC或者0)OUTPR, BASIC, ALL ! 在输出结果中, 列出所有荷载步的基本计算结果OUTPR,BASIC,ALL !指定输出所有节点的基本数据OUTPR,BASIC,LAST ! 选择基本输出选项,直到最后一个荷载步OUTPR,,1 ! 输出第1个荷载步的基本计算结果OUTPR,BASIC,1 ! 选择第1荷载步的基本输出项目OUTPR,NLOAD,1 ! 指定输出第1荷载步的内容OUTRES,ALL,0 !设置将所有数据不记录到数据库。

ANSYS结构分析单元功能与特性/可以组成一一些命令，一般是一种总体命令（session),三十也有特殊，比如是处理/POST1! 是注释说明符号，，与其他软件的说明是一样的，ansys不作为命令读取，* 此符号一般是APDL的标识符，也就是ansys的参数化语言，如*do ,,,*enddo等等NSEL的意思是node select，即选择节点。

s就是select，选择。

DIM 是定义数组的意思。

array 数组。

MP命令用来定义材料参数。

K是建立关键点命令。

K,关键点编号,x坐标,y坐标，z坐标。

K, NPT, X, Y, Z是定义关键点，K是命令，NPT是关键点编号，XYZ是坐标。

NUMMRG, keypoint 用这个命令，要保证关键点的位置完全一样，只是关键点号不一样的才行。

这个命令对于重复的线面都可以用。

这个很简单，压缩关键。

Ngen 复制节点e,节点号码：这个命令式通过节点来形成单元NUMCMP,ALL：压缩所有编号，这样你所有的线都会按次序重新编号~你要是需要固定的线固定的标号NSUBST,100,500,50 ：通过指定子步数来设置载荷步的子步LNSRCH线性搜索是求解非线性代数方程组的一种技巧，此法会在一段区间内，以一定的步长逐步搜索根，相比常用的牛顿迭代法所要耗费的计算量大得多，但它可以避免在一些情况下牛顿迭代法出现的跳跃现象。

LNSRCH 激活线性搜索PRED 激活自由度求解预测NEQIT指定一个荷载步中的最大子步数AUTOTS 自动求解控制打开自动时间步长.KBC -指定阶段状或者用跳板装载里面一个负荷步骤。

SPLINE:P1，P2，P3，P4，P5，P6，XV1，YV1，ZV1，XV6，YV6，ZV6（生成分段样条曲线）*DIM，Par，Type，IMAX，JMAX，KMAX，Var1，Var2，Var3（定义载荷数组的名称）【注】Par: 数组名Type： array 数组，如同fortran,下标最小号为1，可以多达三维（缺省）char 字符串组（每个元素最多8个字符）tableIMAX，JMAX，KMAX 各维的最大下标号Var1，Var2，Var3 各维变量名，缺省为row,column,plane(当type 为table时)/config是设置ansys配置参数的命令格式为/CONFIG, Lab, VALUELab为参数名称 value为参数值例如：/config，MXEL，10000的意思是最大单元数为10000杆单元:LINK1、8、10、11、180梁单元：BEAM3、4、23、24，44，54，188，189管单元:PIPE16，17，18，20，59，602D实体元:PLANE2，25，42，82，83，145，146，182，1833D实体元:SOLID45，46，64，65，72，73，92，95，147，148，185，186，187，191壳单元:SHELL28，41，43，51，61，63，91，93，99，143，150，181，208，209弹簧单元:COMBIN7，14，37，39，40质量单元:MASS21接触单元:CONTAC12，52，TARGE169，170，CONTA171，172，173，174，175，178矩阵单元:MATRIX27，50表面效应元:SURF153，154粘弹实体元:VISCO88，89，106，107，108，超弹实体元:HYPER56，58，74，84，86，158耦合场单元:SOLID5，PLANE13，FLUID29，30，38，SOLID62，FLUID79，FLUID80，81，SOLID98，FLUID129，INFIN110，111，FLUID116，130界面单元:INTER192，193，194，195显式动力分析单元:LINK160，BEAM161，PLANE162，SHELL163，SOLID164，COMBI16杆单元单元名称简称节点数节点自由度特性备注LINK1 2D杆 2 Ux,Uy EPCSDGB常用杆元LINK8 3D杆Ux,Uy,Uz EPCSDGBLINK103D仅受拉或仅受压杆EDGB模拟缆索的松弛及间隙LINK11 3D线性调节器EGB模拟液压缸和大转动LINK183D有限应EPCDFG另可考虑粘0 变杆 B 弹塑性E-弹性(Elasticity)，P-塑性(Plasticity)，C-蠕变(Creep)，S-膨胀(Swelling)，D-大变形或大挠度(Large deflection)，F-大应变(Large strain)或有限应变(Finite strain)，B-单元生死(Birth and dead),G-应力刚化(Stress stiffness)或几何刚度(Geometric stiffening)，A-自适应下降(Adaptive descent)等。

ANSYS应用——m法计算桩基时模拟桩与土作用的命令流!z1spring1.mac/PREP7WPCSYS,-1 !工作平面恢复到默认状态CSYS,0dpile1=1.5!桩直径*****allsel,allcmsel,s,znsll,s,1NSEL,r,LOC,Z,-5-5-5.36,-6-5-5.36!建立泥面以下桩的节点集合CM,SPNODE1,NODEallsel,all,allcmsel,s,znsll,s,1NSEL,r,LOC,Z,-6-5-5.36,-20-5-5.36CM,SPNODE2,NODEallsel,all,allcmsel,s,znsll,s,1NSEL,r,LOC,Z,-20-5-5.36,-25-5-5.36CM,SPNODE3,NODEallsel,all,allcmsel,s,znsll,s,1NSEL,r,LOC,Z,-25-5-5.36,-27-5-5.36CM,SPNODE4,NODEallsel,all,allcmsel,s,znsll,s,1NSEL,r,LOC,Z,-27-5-5.36,-31-5-5.36CM,SPNODE5,NODEallsel,all,allspring1/PREP7CMSEL,S,SPNODE1 !选择泥面以下,桩的节点集合*GET,enum,ELEM,,COUNT, , , ,*GET, NNUM, NODE, 0, COUNT,*GET, nmin, NODE, 0, num, min,*DiM,AA1,ARRAY,NNUM*DiM,AA2,ARRAY,NNUM*DiM,AA3,ARRAY,NNUMCMSEL,S,SPNODE1*do,i,1,nnumAA1(i)=NMINnmin=ndnext(nmin)*enddo !给集合AA1赋值*do,i,1,nnumNSEL,S,,,AA1(i)ESLN,S*GET,emin,ELEM,0,num,min !获取最小的单元编号n1=NELEM(emin,1)n2=NELEM(emin,2)Z12=NZ(N1)-NZ(N2)X12=NX(N1)-NX(N2)*IF,(Z12*X12),GE,0,THENSIG=1*ELSESIG=-1*ENDIFZ12=ABS(Z12)L12=DISTND(N1,N2)CLI=TAN(ACOS(Z12/L12))NGEN,2,100000,AA1(I), , ,-1,,1*CLI*SIGNGEN,2,200000,AA1(I), , ,,1AA2(i)=AA1(i)+100000AA3(i)=AA1(i)+200000TYPE,3 !将单元赋为弹簧单元,为etsp单元类型,见begin!MAT,mpsp !mpsp为材料号,见beginESYS,0R,enum+i,8000*(-10-5.36-NZ(AA1(I)))*(L12)*dpile1, !根据M法给弹簧赋予刚度,其中zlsize为桩单元长!!REAL,enum+ie,AA1(I),AA2(I)e,AA1(I),AA3(I)*enddo/SOLALLSEL,ALL*do,i,1,NNUMD,AA2(i),allD,AA3(i),all*ENDDOspring1spring2/PREP7CMSEL,S,SPNODE2 !选择泥面以下,桩的节点集合*GET,enum,ELEM,,COUNT, , , ,*GET, NNUM, NODE, 0, COUNT,*GET, nmin, NODE, 0, num, min,*DiM,AA1,ARRAY,NNUM*DiM,AA2,ARRAY,NNUM*DiM,AA3,ARRAY,NNUMCMSEL,S,SPNODE2*do,i,1,nnumAA1(i)=NMINnmin=ndnext(nmin)*enddo !给集合AA1赋值*do,i,1,nnumNSEL,S,,,AA1(i)ESLN,S*GET,emin,ELEM,0,num,min !获取最小的单元编号n1=NELEM(emin,1)n2=NELEM(emin,2)Z12=NZ(N1)-NZ(N2)X12=NX(N1)-NX(N2)*IF,(Z12*X12),GE,0,THENSIG=1*ELSESIG=-1*ENDIFZ12=ABS(Z12)L12=DISTND(N1,N2)CLI=TAN(ACOS(Z12/L12))NGEN,2,100000,AA1(I), , ,-1,,1*CLI*SIGNGEN,2,200000,AA1(I), , ,,1AA2(i)=AA1(i)+100000AA3(i)=AA1(i)+200000TYPE,3 !将单元赋为弹簧单元,为etsp单元类型,见begin!MAT,mpsp !mpsp为材料号,见beginESYS,0R,enum+i,4500*(-10-5.36-NZ(AA1(I)))*(L12)*dpile1, !根据M法给弹簧赋予刚度,其中zlsize为桩单元长!!REAL,enum+ie,AA1(I),AA2(I)e,AA1(I),AA3(I)*enddo/SOLALLSEL,ALL*do,i,1,NNUMD,AA2(i),allD,AA3(i),all*ENDDOspring2spring3/PREP7CMSEL,S,SPNODE3 !选择泥面以下,桩的节点集合*GET,enum,ELEM,,COUNT, , , ,*GET, NNUM, NODE, 0, COUNT,*GET, nmin, NODE, 0, num, min,*DiM,AA1,ARRAY,NNUM*DiM,AA2,ARRAY,NNUM*DiM,AA3,ARRAY,NNUMCMSEL,S,SPNODE3*do,i,1,nnumAA1(i)=NMINnmin=ndnext(nmin)*enddo !给集合AA1赋值*do,i,1,nnumNSEL,S,,,AA1(i)ESLN,S*GET,emin,ELEM,0,num,min !获取最小的单元编号n1=NELEM(emin,1)n2=NELEM(emin,2)Z12=NZ(N1)-NZ(N2)X12=NX(N1)-NX(N2)*IF,(Z12*X12),GE,0,THENSIG=1*ELSESIG=-1*ENDIFZ12=ABS(Z12)L12=DISTND(N1,N2)CLI=TAN(ACOS(Z12/L12))NGEN,2,100000,AA1(I), , ,-1,,1*CLI*SIGNGEN,2,200000,AA1(I), , ,,1AA2(i)=AA1(i)+100000AA3(i)=AA1(i)+200000TYPE,3 !将单元赋为弹簧单元,为etsp单元类型,见begin!MAT,mpsp !mpsp为材料号,见beginESYS,0R,enum+i,10000*(-10-5.36-NZ(AA1(I)))*(L12)*dpile1, !根据M法给弹簧赋予刚度,其中zlsize为桩单元长!!REAL,enum+ie,AA1(I),AA2(I)e,AA1(I),AA3(I)*enddo/SOLALLSEL,ALL*do,i,1,NNUMD,AA2(i),allD,AA3(i),all*ENDDOspring3spring4/PREP7CMSEL,S,SPNODE4 !选择泥面以下,桩的节点集合*GET,enum,ELEM,,COUNT, , , ,*GET, NNUM, NODE, 0, COUNT,*GET, nmin, NODE, 0, num, min,*DiM,AA1,ARRAY,NNUM*DiM,AA2,ARRAY,NNUM*DiM,AA3,ARRAY,NNUMCMSEL,S,SPNODE3*do,i,1,nnumAA1(i)=NMINnmin=ndnext(nmin)*enddo !给集合AA1赋值*do,i,1,nnumNSEL,S,,,AA1(i)ESLN,S*GET,emin,ELEM,0,num,min !获取最小的单元编号n1=NELEM(emin,1)n2=NELEM(emin,2)Z12=NZ(N1)-NZ(N2)X12=NX(N1)-NX(N2)*IF,(Z12*X12),GE,0,THENSIG=1*ELSESIG=-1*ENDIFZ12=ABS(Z12)L12=DISTND(N1,N2)CLI=TAN(ACOS(Z12/L12))NGEN,2,100000,AA1(I), , ,-1,,1*CLI*SIGNGEN,2,200000,AA1(I), , ,,1AA2(i)=AA1(i)+100000AA3(i)=AA1(i)+200000TYPE,3 !将单元赋为弹簧单元,为etsp单元类型,见begin!MAT,mpsp !mpsp为材料号,见beginESYS,0R,enum+i,20000*(-10-5.36-NZ(AA1(I)))*(L12)*dpile1, !根据M法给弹簧赋予刚度,其中zlsize为桩单元长!!REAL,enum+ie,AA1(I),AA2(I)e,AA1(I),AA3(I)*enddo/SOLALLSEL,ALL*do,i,1,NNUMD,AA2(i),allD,AA3(i),all*ENDDOspring4spring5/PREP7CMSEL,S,SPNODE5 !选择泥面以下,桩的节点集合*GET,enum,ELEM,,COUNT, , , ,*GET, NNUM, NODE, 0, COUNT,*GET, nmin, NODE, 0, num, min,*DiM,AA1,ARRAY,NNUM*DiM,AA2,ARRAY,NNUM*DiM,AA3,ARRAY,NNUMCMSEL,S,SPNODE3*do,i,1,nnumAA1(i)=NMINnmin=ndnext(nmin)*enddo !给集合AA1赋值*do,i,1,nnumNSEL,S,,,AA1(i)ESLN,S*GET,emin,ELEM,0,num,min !获取最小的单元编号n1=NELEM(emin,1)n2=NELEM(emin,2)Z12=NZ(N1)-NZ(N2)X12=NX(N1)-NX(N2)*IF,(Z12*X12),GE,0,THENSIG=1*ELSESIG=-1*ENDIFZ12=ABS(Z12)L12=DISTND(N1,N2)CLI=TAN(ACOS(Z12/L12))NGEN,2,100000,AA1(I), , ,-1,,1*CLI*SIGNGEN,2,200000,AA1(I), , ,,1AA2(i)=AA1(i)+100000AA3(i)=AA1(i)+200000TYPE,3 !将单元赋为弹簧单元,为etsp单元类型,见begin!MAT,mpsp !mpsp为材料号,见beginESYS,0R,enum+i,30000*(-10-5.36-NZ(AA1(I)))*(L12)*dpile1, !根据M法给弹簧赋予刚度,其中zlsize为桩单元长!!REAL,enum+ie,AA1(I),AA2(I)e,AA1(I),AA3(I)*enddo/SOLALLSEL,ALL*do,i,1,NNUMD,AA2(i),allD,AA3(i),all*ENDDOspring5。

Ansys命令流大全（整理）1、A,P1,P2,P3,P4,P5,P6,P7,P8,P9此命令用已知的一组关键点点（P1~P9）来定义面（Area），最少使用三个点才能围成面，同时产生转围绕些面的线。

点要依次序输入，输入的顺序会决定面的法线方向。

如果超过四个点，则这些点必须在同一个平面上。

Menu Paths:Main Menu>Preprocessor>Create>Arbitrary>Through KPs2、*ABBR，Abbr,String－－定义一个缩略语．Abbr:用来表示字符串＂String＂的缩略语，长度不超过８个字符．String：将由＂Abbr＂表示的字符串，长度不超过６０个字符．3、ABBRES，Lab，Fname，Ext－从一个编码文件中读出缩略语．Lab：指定读操作的标题，NEW：用这些读出的缩略语重新取代当前的缩略语（默认）CHANGE：将读出的缩略语添加到当前缩略语阵列，并替代现存同名的缩略语．Ext：如果＂Fname＂是空的，则缺省的扩展命是＂ABBR＂．4、ABBSA V，Lab，Fname，Ext－将当前的缩略语写入一个文本文件里Lab：指定写操作的标题，若为ALL，表示将所有的缩略语都写入文件（默认）5、add, ir, ia,ib,ic,name,--,--,facta, factb, factc将ia,ib,ic变量相加赋给ir变量ir, ia,ib,ic：变量号name: 变量的名称6、Adele,na1,na2,ninc,kswp ！kswp=0时只删除掉面积本身，＝1时低单元点一并删除。

7、Adrag, nl1,nl2,nl3,nl4,nl5,nl6, nlp1,nlp2,nlp3,nlp4,nlp5,nlp6 ！面积的建立，沿某组线段路径，拉伸而成。

8、Afillt,na1,na2,rad ！建立圆角面积，在两相交平面间产生曲面，rad为半径。

ANSYS中两种模拟桩土共同作用的方法及对比作者：孙浩马孟启来源：《科学与财富》2011年第01期[摘要] ANSYS中模拟桩土共同作用的常用方法有两种，分别是运用接触问题模拟和运用弹簧单元模拟，本文介绍了这两种模拟方法，并通过一个工程实例作了对比。

[关键词] 非线性共同作用接触问题弹簧单元1、桩土共同作用常用的两种模拟方法在桥梁结构的整体计算中，有时候在荷载作用下的变位是必须考虑的。

桥梁的基础以桩基最常见，模拟桩土作用的方法大致有两种：一是将桩、承台和土甚至上部结构整体建模，土看成弹塑性材料，用适用于D-P材料的solid45模拟，桩、承台和土间用接触单元连接；二是将土看成弹簧，用弹簧或者LINK单元模拟。

第二种方法又可以分为两类，一类是算出整个桩基在土面处的水平、竖向和转动刚度，然后直接在土面处施加弹簧单元，这样就没有桩了；二类是将一定厚度的土层的刚度系数计算出来，然后再这层土的中间处施加弹簧，这时可以用梁单元模拟桩，但单元的抗弯刚度取实际值的0.8倍。

2、用D-P材料模拟2.1接触问题的基本常识此方法主要运用到接触问题，接触问题是一种高度非线性行为，它存在两大难点：其一，在求解问题之前，不知道接触区域，表面之间是接触或分开是未知的、突然变化的，这些随载荷、材料边界条件和其他因素而定；其二，大多数的接触问题需要计算摩擦，有几种摩擦的模型可供挑选，它们都是非线性的，摩擦使问题的收敛性变得困难。

接触问题分为两个基本类型：刚体-柔体的接触和柔体-柔体接触。

ANSYS支持三种接触方式：点-点、点-面和面-面。

每种接触方式使用的接触单元适用于某类问题。

2.2运用接触问题模拟桩土作用的方法桩土的共同作用属于面-面接触问题。

在模拟桩土作用这样的3D的接触单元中，ANSYS 支持刚体-柔体的面-面接触单元，刚性面被当作“目标”面，分别用Target169和Target170来模拟2-D和3-D的“目标”面，柔性体的表面被当作“接触”面，用Conta171、Conta172、Conta173、Conta174来模拟。

第27卷第6期 岩 土 力 学 V ol.27 No.6 2006年6月 Rock and Soil Mechanics Jun. 2006收稿日期：2004-09-08作者简介：唐世栋，男，1952年生，博士，副教授，从事岩土工程教学和研究工作。

文章编号：1000－7598－(2006) 06―0973―04基于ANSYS 软件模拟桩的挤入过程唐世栋，李 阳（同济大学 地下建筑与工程系，上海 200092）摘 要：基于ANSYS 软件分析了桩土之间的相互作用，模拟了桩打入时土中的应力、应变情况。

通过结合ANSYS 中的接触分析和生死单元，以DP 材料来模拟土体，采用循环命令的方式来分析桩土接触时复杂的应力状态。

模拟结果得到了圆孔扩张理论和极限平衡法的验证。

关 键 词：ANSYS ；桩；桩土作用；Drucker-Prager 屈服准则 中图分类号：O 245 TU 473.1+2 文献标识码：AAnalysis of a driven pile by ANSYSTANG Shi-dong, LI Yang(Department of Geotechnical Engineering, Tongji University, Shanghai 200092, China)Abstract: The interaction of soil and pile is analyzed by ANSYS; and the state of stress and strain in the soil after the pile was driven in is simulated. By combining the contact analysis with birth - death element of ANSYS, and proposed the soil material comply with the Drucker-Prager criterion, the circulation command is used to analyze complicated soil stress state. These results are verified by limit equilibrium method and theory of expansion of a cavity.Key words: ANSYS; pile; interaction of soil and pile; Drucker-Prager criterion1 概 述挤土桩在打入的过程中，对周围土体的挤密会使地表隆起、桩侧产生位移应力分区、桩尖下形成应力泡等现象[1]。

有限元分析软件ANSYS命令流中文说明4 4有限元分析软件ANSYS命令流中文说明4/42010-05-23 21：151设置分析类型ANTYPE，Antype,status,ldstep,action其中antype表示分析类型STATIC：静态分析MODAL：模态分析TRANS：瞬态分析SPECTR：谱分析2 KBC，KEY制定载荷为阶跃载荷还是递增载荷EKY=0递增方式KEY=1阶跃方式3 SOLVE开始一个求解运算4 LSSOLVE读入并求解多个载荷步5 TIME，time设置求解时间有时在分析中需要进入后处理，然后在保持进入后处理之前的状态的情况下接着算下去，可以使用以下的方法：PARSAV,ALL,PAR,TXT！PARSAV命令是储存ANSYS的参数，ALL代表所有参数，PAR是文件名，TXT是扩展名/SOLU ANTYPE,REST,CruStep-1,,CONTINUE！ANTYPE是定义分析类型的命令，REST代表重启动，CruStep代表本载荷步的编号PARRES,NEW,PAR,TXT！PARRES是恢复参数的命令，NEW表示参数是以刷新状态恢复，PAR和TXT 代表了储存了参数的文件名和扩展名如果有单元生死的问题，可以这样处理：ALLSEL,ALL*GET,E_SUM_MAX,ELEM,NUM,MAX！得到单元的最大编号，即单元的总数ESEL,S,LIVE！选中"生"的单元*GET,E_SUM_AL,ELEM,COUNT*DIM,E_POT_AL,E_SUM_MAX！单元选择的指示*DIM,E_NUM_AL,E_SUM_AL！单元编号的数组J=0！读出所选单元号*DO,I,1,E_SUM_MAX*VGET,E_POT_AL(I),ELEM,I,ESEL！对所有单元做循环，被选中的单元标志为"1"*IF,E_POT_AL(I),EQ,1,THEN J=J+1 E_NUM_AL(J)=I*ENDIF*ENDDO ALLSEL,ALL在重启动之后恢复单元生死状态*if,E_SUM_AL,ne,0,then*do,i,1,Num_Alive esel,a,E_NUM_AL(i)*enddo ealive,all allsel*endif/WINDOW,WN,XMIN,XMAX,YMIN,YMAX,NCOPY注意x的坐标是-1到1.67，y坐标是-1到1 Xmin=off on，FULL,LEFT,RIGH,TOP,BOT,LTOP,LBOT,RTOP,RBOT注意一个问题，除了1号窗口外，其他的不能用鼠标操作，只用先发/view 和/dist，然后用/replot。

ANSYS命令流详解（超全）一、定义材料号及特性mp,lab, mat, co, c1,…….c4lab: 待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens）ex: 弹性模量nuxy: 小泊松比alpx: 热膨胀系数reft: 参考温度reft: 参考温度prxy: 主泊松比gxy: 剪切模量mu: 摩擦系数dens: 质量密度mat: 材料编号（缺省为当前材料号）c 材料特性值，或材料之特性，温度曲线中的常数项c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数二、定义DP材料：首先要定义EX和泊松比：MP，EX，MAT，……MP，NUXY，MAT，……定义DP材料单元表（这里不考虑温度）：TB，DP，MAT进入单元表并编辑添加单元表：TBDATA，1，CTBDATA，2，ψTBDATA，3，……如定义：EX=1E8，NUXY=0.3，C=27，ψ=45的命令如下：MP，EX，1，1E8MP，NUXY，1，0.3TB，DP，1TBDATA，1，27TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,deg三、单元生死载荷步!第一个载荷步TIME,... !设定时间值（静力分析选项）NLGEOM,ON !打开大位移效果NROPT,FULL !设定牛顿－拉夫森选项ESTIF,... !设定非缺省缩减因子（可选）ESEL,... !选择在本载荷步中将不激活的单元EKILL,... !不激活选择的单元ESEL,S,LIVE !选择所有活动单元NSLE,S !选择所有活动结点NSEL,INVE !选择所有非活动结点（不与活动单元相连的结点）D,ALL,ALL,0 !约束所有不活动的结点自由度（可选）NSEL,ALL !选择所有结点ESEL,ALL !选择所有单元D,... !施加合适的约束F,... !施加合适的活动结点自由度载荷SF,... !施加合适的单元载荷BF,... !施加合适的体载荷SAVESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D,F,SF和BF命令得到更详细的解释。

桩 -土荷载传递规律 ANSYS模型分析摘要：本文主要介绍了岩土工程中桩土体系荷载传递的相关理论，重点描述了桩基础在ANSYS模型中位移和应力的变化特点，根据ANSYS模型木桩数据分析的结果，模拟分析桩土体系的荷载传递原理。

关键词：ANSYS建模桩基础荷载传递有限元引言目前国内外很多科学工作者采用了实验模型对桩基进行了研究，并获得大量有益的成果，同时推动了岩土工程理论研究上的进展，但是工程建设中，桩基工程施工时实验模拟和现场测定均会带来更多的人为误差，根据ANSYS模型对桩土作用研究即可减少这种误差，只要正确的设定参数，将会方便快捷的掌握桩基础荷载传递的规律。

1桩-土体系ANSYS模型的建立1.1桩-土体系的基本参数设定本文以木桩的ANSYS10.0模型进行建模分析，通过模拟木桩的实际荷载传递规律来分析桩土中荷载传递机理以及相关特性。

在ANSYS建模过程中，木桩的尺寸采用如下数值：桩径34.12mm，桩高度30cm的圆柱形形状，地基土采用天然黄土，假设木桩埋深220cm。

接触单元的弹性模量采取2E+4 ；泊松比采取0.499，摩擦系数0.2，木桩的弹性模量采取8E+3 ；泊松比采用0.38。

桩在竖向荷载下的工作过程中非常复杂，所以其特性取决于多方面的变通因素，在桩土的作用分析中，一系列的简化是不可避免的。

因此在本文在创建ANSYS模型创建的过程中，必须对其真实模型还需要作进一步的简化和处理，具体如下：(1)考虑到土弹塑性的本构关系，一般土层采用Drucker-Prager(DP)材料模型。

DP材料的材料特性值包括粘聚力C、内摩擦角φ和膨胀角φf 。

膨胀角φf被用来控制体积膨胀的大小，对压实的颗粒状材料，当材料受剪时，颗粒将会膨胀，如果膨胀角φf=0，则不会发生体积膨胀，其塑性行为被假定为理想弹塑性。

另外，此种材料还考虑了由于屈服而引起的体积膨胀，但不考虑温度变化的影响，适用于混凝土、岩石和土壤等颗粒状材料。

用MIDAS模拟桩-土相互作用（“m法”确定土弹簧刚度）北京迈达斯技术有限公司2009年05月1、引言土与结构相互作用的研究已有近60～70年的历史，待别是近30年来，计算机技术的发展为其提供了有力的分析手段。

桩基础是土建工程中广泛采用的基础形式之一，许多建于软土地基上的大型桥梁结构往往都采用桩基础，桩-土动力相互作用又是土-结构相互作用问题中较复杂的课题之一。

至今已有不少关于桩基动力特性的研究报告，国内外研究人员也提出了许多不同的桩-土动力相互作用计算方法。

从研究成果的归类来看，理论上主要有离散理论和连续理论及两者的结合，解决的方法一般有集中质量法、有限元法、边界元法和波动场法。

60～70年代，美国学者J.penzien等在解决泥沼地上大桥动力分析时提出了集中质量法，目前已在国内外得到了广泛的应用。

集中质量法将桥梁上部结构多质点体系和桩一土体系的质量联合作为一个整体，来建立整体耦联的地震振动微分方程组进行求解。

该模型假定桩侧土是Winkler连续介质。

以半空间的Mindlin静力基本解为基础，将桩-土体系的质量按一定的厚度简化并集中为一系列质点，离散成一理想化的参数系统。

并用弹簧和阻尼器模拟土介质的动力性质，形成一个包括地下部分的多质点体系。

土弹簧刚度的确定，除考虑使用较为精确的有限元或边界元方法外，较为简便的方法是采用Penzien模型中提供的土弹簧计算方法或参照现行规范中土弹簧的计算方法。

我国公路桥涵地基与基础设计规范(JTG D63-2007)用的“m法”计算方法和参数选取方面比Penzien 的方法要简单和方便，且为国内广大工程师所熟.“m法”的基本原理是将桩作为弹性地基梁，按Winkler假定（梁身任一点的土抗力和该点的位移成正比）求解。

但是，由于桩-土相互作用的实验数据不足，土的物性取值有时亦缺乏合理性，在确定土弹簧的刚度时，仍有不少问题未能很好解决。

特别是，“m法”中m的取值对弹簧刚度的计算结果影响很大，且不能反映地震波的频率特性和强度带来的影响。

1．ANSYS SOLID65环向布置钢筋的例子 (3)2.混凝土非线性计算实例（1）- MISO单压 (5)3.混凝土非线性计算实例（2）-MISO约束压 (6)4.混凝土非线性计算实例（3）- KINH滞回 (9)5.混凝土非线性计算实例（4）- KINH压-拉裂 (11)6.混凝土非线性计算实例（5） (12)7.混凝土非线性计算实例（6） (14)8.混凝土非线性计算实例（7）-MISO滞回 (16)9.混凝土非线性计算实例（8） (18)10.混凝土非线性计算实例（9）-梁平面应力 (20)11.四层弹簧－质点模型的地震分析 (22)12.悬臂梁地震分析 (48)13.用beam 54单元描述变截面梁的例子 (72)14.变截面梁实例 (73)15.拱桥浇筑过程分析-单元生死应用实例 (74)16.简支梁实体与预应力钢筋分析实例 (75)17. 简单的二维焊接分析－单元生死实例 (77)18.隧道开挖（三维）的命令流 (84)19.岩土接触分析实例 (101)20.钢筋混凝土管的动力响应特性分析实例 (109)21.隧道模拟开挖命令流（入门） (116)22.螺栓连接的模拟实现问题 (119)23.道路的基层、垫层模量与应力之间的关系 (129)23.滞回分析 (151)24.模拟某楼层浇注 (153)25.在面上施加移动的面力 (155)27.在任意面施加任意方向任意变化的压力 (159)28.预紧分析 (160)29.几何非线性+塑性+接触+蠕变 (162)30.埋设在地下的排水管道 (167)32.幕墙企业玻璃简化计算 (172)33.等截面杆单元生死应用实例 (188)34.梁板建模联系 (189)36.简单的例子－如何对结构的振动控制分析 (192)37.模态分析结果的输出实例 (194)38.火车过桥动态加载实例（部分） (196)39.悬索结构的找形和计算的例题 (213)40.陶瓷杆撞击铝板的例子 (218)41.求反作用力的APDL命令法 (221)42.LS-DYNA实例（部分） (222)43.路面分层填筑对路基的影响 (223)44.一个例子（含地震影响，求振兴与频率） (227)45.接触面上的压力总和 (231)46.施加位置函数荷载 (235)47.非线性分析考虑刚度退化 (236)48.一个圆形水池的静力分析 (237)49.ANSYS中混凝土模式预应力模拟的算例 (238)50.悬臂梁受重力作用发生大变形求其固有频率 (240)51.循环对称结构模态分析 (242)52.三角平台受谐波载荷作用的结构响应 (244)53.三角平台受一地震谱激励的应力分布和支反力 (246)54.三角平台受时程载荷作用的应力分布和变形过程 (248)55.经典层合板理论 (250)56．定易圆轨迹的例子 (257)57．模拟门式刚架施工－单元生死 (257)58．钢筋混凝土整体式模型例子 (260)59．在荷载步之间改变材料属性例子 (262)60．含预应力的特征值屈曲计算 (263)61．振型叠加计算及工况组合例子 (265)62．柱子稳定分析算（预应力，特征值屈曲，初始缺陷） (268)63．module MConcrete !混凝土模板 (271)64.混凝土开裂实例 (279)65．螺栓网格划分 (280)66．自由液面的土石坝平面渗流分析 (281)67．导出刚度矩阵 (285)68．某混凝土拱坝工程施工期及运行期温度场仿真分析 (286)69．移动温度荷载计算 (293)70．SHSD用于壳-实体装配实例An (295)71．ansys显示－隐式－回弹分析实例 (299)72．工况组合的经典例子 (314)1．ANSYS SOLID65环向布置钢筋的例子! 一个管道，环向配筋率为1%，纵向配筋率为0.5%，径向配筋率为0.1%!FINISH/CLEAR/PREP7!*! 单元属性ET,1,SOLID65!*KEYOPT,1,1,0KEYOPT,1,5,0KEYOPT,1,6,0KEYOPT,1,7,1!*!实参数1：不同方向配筋R,1,2,.001, , ,2, .01,RMORE, 90, ,2,.005 ,90 ,90 ,!Adds real constants to a set.!材料属性!混凝土基本材料属性MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,30E2MPDATA,PRXY,1,,.2!屈服准则TB,MISO,1,1,5,TBTEMP,0TBPT,,0.0005,15TBPT,,0.001,21TBPT,,0.0015,24TBPT,,0.002,27TBPT,,0.003,24!破坏准则TB,CONC,1,1,9,TBTEMP,0TBDATA,,.5,.9,3,30,,!抗拉和抗压是10 倍的关系TBDATA,,,,1,,,MPTEMP,,,,,,,,MPTEMP,1,0!钢材基本属性MPDATA,EX,2,,200E3MPDATA,PRXY,2,,.27!屈服准则TB,BISO,2,1,2,TBTEMP,0TBDATA,,310,2E3,,,,! 管道内径和外径CYL4,0,0,3000, , , ,10000CYL4,0,0,2000, , , ,10000VSBV, 1, 2! 定义局部柱坐标CSWPLA,11,1,1,1,! Defines a local coordinate system at the origin of the working plane. KWPAVE, 11!move to keypoint11WPRO,,-90.000000,VSBW, 3WPCSYS,-1,0!Defines the working plane location based on a coordinate system. KWPAVE, 1WPRO,,,-90.000000VSBW,ALLESIZE,500,0,! 注意：设定单元局部坐标VATT, 1, 1, 1, 11!*VSWEEP,ALL/DEVICE,VECTOR,1/ESHAPE,1.0!Displays elements with shapes determined from the real constants or section/REPLO! 注意：红色代表最大配筋方向，绿色代表其次，蓝色表示最小配筋方向2.混凝土非线性计算实例（1）- MISO单压!MISO单压FINISH/CLEAR/PREP7ET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,MISO,1,1,15TBPT,, 100E-6, 3.0TBPT,, 300E-6, 8.3TBPT,, 600E-6,14.6TBPT,, 900E-6,19.1TBPT,,1100E-6,21.0TBPT,,1250E-6,22.0TBPT,,1400E-6,22.6TBPT,,1550E-6,22.8TBPT,,1650E-6,22.7TBPT,,1800E-6,22.3TBPT,,2000E-6,21.4TBPT,,2800E-6,16.8TBPT,,3200E-6,14.7TBPT,,3800E-6,12.3TBPT,,4600E-6,9.9!TB,MKIN,1!TBTEMP,,STRAIN!TBDATA,,600E-6,1100E-6,1600E-6,3000E-6,4500E-6 !TBTEMP,0!TBDATA,,15,21.5,23,16,9BLOCK,0,50,0,50,0,50MSHAPE,0,3DMSHKEY,1ESIZE,10VMESH,ALLNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50CP,1,UY,ALLNSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUANTYPE,STATICAUTOTS,OFFOUTRES,ALL,1TIME,50NSUBST,50D,1,UY,-.25NSEL,ALLSOLVEFINISH/POST26NSOL,2,1,U,Y,UYRFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50ADD,5,3,,,STRESS,,,-1/2500/AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH3.混凝土非线性计算实例（2）-MISO约束压!MISO约束压FINISH/CLEAR/PREP7ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,MISO,1,1,15TBPT,, 100E-6, 3.0TBPT,, 300E-6, 8.3TBPT,, 600E-6,14.6TBPT,, 900E-6,19.1TBPT,,1100E-6,21.0TBPT,,1250E-6,22.0TBPT,,1400E-6,22.6TBPT,,1550E-6,22.8TBPT,,1650E-6,22.7TBPT,,1800E-6,22.3TBPT,,2000E-6,21.4TBPT,,2800E-6,16.8TBPT,,3200E-6,14.7TBPT,,3800E-6,12.3TBPT,,4600E-6,9.9!TB,MKIN,1!TBTEMP,,STRAIN!TBDATA,,600E-6,1100E-6,1600E-6,3000E-6,4500E-6 !TBTEMP,0!TBDATA,,15,21.5,23,16,9BLOCK,0,50,0,50,0,50ESIZE,10VMESH,ALLMSHAPE,0,3DNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50CP,1,UY,ALLD,ALL,UX,0D,ALL,UZ,0NSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUOUTRES,ALL,1TIME,20NSUBST,20,0,20D,1,UY,-.04LSWRITE,1TIME,30NSUBST,20,0,20D,1,UY,-.06LSWRITE,2TIME,40NSUBST,100D,1,UY,-.1LSWRITE,3LSSOLVE,1,3FINISH/POST26NSOL,2,1,U,Y,UY RFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50 ADD,5,3,,,STRESS,,,-1/2500 /AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH4.混凝土非线性计算实例（3）- KINH滞回!KINH滞回FINISH/CLEAR/PREP7ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,KINH,1,1,10TBPT,, 150E-6, 4.5TBPT,, 600E-6,14.8TBPT,,1000E-6,20.25TBPT,,1300E-6,22.3TBPT,,1480E-6,22.8TBPT,,1620E-6,22.8TBPT,,1800E-6,22.3TBPT,,2000E-6,21.4TBPT,,3500E-6,12.8TBPT,,5000E-6, 9.0BLOCK,0,50,0,50,0,50ESIZE,10VMESH,ALLMSHAPE,0,3DNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50CP,1,UY,ALLNSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUOUTRES,ALL,ALLTIME,20NSUBST,20,0,20D,1,UY,-.04LSWRITE,1TIME,35NSUBST,15,0,15D,1,UY,-.01LSWRITE,2TIME,90NSUBST,55,0,55D,1,UY,-.12LSWRITE,3TIME,125NSUBST,35,0,35D,1,UY,-.05LSWRITE,4TIME,175NSUBST,50,0,50D,1,UY,-.15LSWRITE,5LSSOLVE,1,5FINISH/POST26NSOL,2,1,U,Y,UY RFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50 ADD,5,3,,,STRESS,,,-1/2500 /AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH5.混凝土非线性计算实例（4）- KINH压-拉裂!KINH压——拉裂压到峰值，泄载、反向加载到拉裂FINISH/CLEAR/PREP7ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,KINH,1,1,10TBPT,, 150E-6, 4.5TBPT,, 600E-6,14.8TBPT,,1000E-6,20.25TBPT,,1300E-6,22.3TBPT,,1480E-6,22.8TBPT,,1620E-6,22.8TBPT,,1800E-6,22.3TBPT,,2000E-6,21.4TBPT,,3500E-6,12.8TBPT,,5000E-6, 9.0BLOCK,0,50,0,50,0,50ESIZE,10VMESH,ALLMSHAPE,0,3DNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50CP,1,UY,ALLNSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUOUTRES,ALL,ALLTIME,20NSUBST,20,0,20D,1,UY,-.04LSWRITE,1TIME,45NSUBST,25,0,25D,1,UY,.01LSWRITE,2LSSOLVE,1,2FINISH/POST26NSOL,2,1,U,Y,UYRFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50ADD,5,3,,,STRESS,,,-1/2500/AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH6.混凝土非线性计算实例（5）FINISH/CLEAR/PREP7ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1 TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1 TB,KINH,1,1,18 TBPT,, 100E-6, 3.0 TBPT,, 300E-6, 8.3 TBPT,, 600E-6,14.6 TBPT,, 900E-6,19.1 TBPT,,1100E-6,21.0 TBPT,,1250E-6,22.0 TBPT,,1400E-6,22.6 TBPT,,1550E-6,22.8 TBPT,,1650E-6,22.7 TBPT,,1800E-6,22.3 TBPT,,2000E-6,21.4 TBPT,,2400E-6,19.1 TBPT,,2800E-6,16.8 TBPT,,3200E-6,14.7 TBPT,,3600E-6,13.0 TBPT,,4100E-6,11.3 TBPT,,4600E-6,9.9BLOCK,0,50,0,50,0,50 ESIZE,10VMESH,ALLMSHAPE,0,3DNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50CP,1,UY,ALLNSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUOUTRES,ALL,ALLTIME,50NSUBST,50,0,50D,1,UY,-.25SOLVEFINISH/POST26NSOL,2,1,U,Y,UYRFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50ADD,5,3,,,STRESS,,,-1/2500/AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH7.混凝土非线性计算实例（6）!MISO约束压FINISH/CLEAR/PREP7ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,KINH,1,1,15TBPT,, 100E-6, 3.0TBPT,, 300E-6, 8.3TBPT,, 600E-6,14.6TBPT,, 900E-6,19.1TBPT,,1100E-6,21.0TBPT,,1250E-6,22.0TBPT,,1400E-6,22.6TBPT,,1550E-6,22.8TBPT,,1650E-6,22.7TBPT,,1800E-6,22.3TBPT,,2000E-6,21.4TBPT,,2800E-6,16.8TBPT,,3200E-6,14.7TBPT,,3800E-6,12.3TBPT,,4600E-6,9.9!TB,MKIN,1!TBTEMP,,STRAIN!TBDATA,,600E-6,1100E-6,1600E-6,3000E-6,4500E-6 !TBTEMP,0!TBDATA,,15,21.5,23,16,9BLOCK,0,50,0,50,0,50ESIZE,10VMESH,ALLMSHAPE,0,3DNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50CP,1,UY,ALLD,ALL,UX,0D,ALL,UZ,0NSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUOUTRES,ALL,1TIME,20NSUBST,20,0,20D,1,UY,-.04LSWRITE,1TIME,30NSUBST,20,0,20D,1,UY,-.06LSWRITE,2TIME,40NSUBST,100D,1,UY,-.1LSWRITE,3LSSOLVE,1,3FINISH/POST26NSOL,2,1,U,Y,UYRFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50ADD,5,3,,,STRESS,,,-1/2500/AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH8.混凝土非线性计算实例（7）-MISO滞回!MISO滞回FINISH/CLEAR/PREP7ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,MISO,1,1,15 TBPT,, 100E-6, 3.0 TBPT,, 300E-6, 8.3 TBPT,, 600E-6,14.6 TBPT,, 900E-6,19.1 TBPT,,1100E-6,21.0 TBPT,,1250E-6,22.0 TBPT,,1400E-6,22.6 TBPT,,1550E-6,22.8 TBPT,,1650E-6,22.7 TBPT,,1800E-6,22.3 TBPT,,2000E-6,21.4 TBPT,,2800E-6,16.8 TBPT,,3200E-6,14.7 TBPT,,3800E-6,12.3 TBPT,,4600E-6,9.9BLOCK,0,50,0,50,0,50 ESIZE,10VMESH,ALL MSHAPE,0,3DNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50 CP,1,UY,ALLNSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUOUTRES,ALL,ALL TIME,20NSUBST,20,0,20D,1,UY,-.04 LSWRITE,1TIME,35NSUBST,15,0,15D,1,UY,-.01LSWRITE,2TIME,60NSUBST,25,0,25D,1,UY,-.06LSWRITE,3TIME,80NSUBST,20,0,20D,1,UY,-.02LSWRITE,4TIME,110NSUBST,30,0,30D,1,UY,-.08LSWRITE,5LSSOLVE,1,5FINISH/POST26NSOL,2,1,U,Y,UYRFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50ADD,5,3,,,STRESS,,,-1/2500/AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH9.混凝土非线性计算实例（8）!MISO单压FINISH/CLEAR/PREP7ANTYPE,STATICET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,MISO,1,1,10TBPT,, 150E-6, 4.5TBPT,, 600E-6,14.8TBPT,,1000E-6,20.25TBPT,,1300E-6,22.3TBPT,,1480E-6,22.8TBPT,,1620E-6,22.8TBPT,,1800E-6,22.3TBPT,,2000E-6,21.4TBPT,,3500E-6,12.8TBPT,,5000E-6, 9.0!TB,MKIN,1!TBTEMP,,STRAIN!TBDATA,,600E-6,1100E-6,1600E-6,3000E-6,4500E-6 !TBTEMP,0!TBDATA,,15,21.5,23,16,9BLOCK,0,50,0,50,0,50ESIZE,10VMESH,ALLMSHAPE,0,3DNSEL,S,LOC,X,0D,ALL,UX,0NSEL,S,LOC,Y,0D,ALL,UY,0NSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,50CP,1,UY,ALLNSEL,ALLFINISH/VIEW,1,1,1,1/REPLOT/SOLUOUTRES,ALL,1TIME,50NSUBST,100,0,100D,1,UY,-.25SOLVEFINISH/POST26NSOL,2,1,U,Y,UYRFORCE,3,1,F,Y,FYADD,4,2,,,STRAIN,,,-1/50ADD,5,3,,,STRESS,,,-1/2500/AXLAB,X,STRAIN/AXLAB,Y,STRESSXVAR,4PLVAR,5FINISH10.混凝土非线性计算实例（9）-梁平面应力!梁平面应力FINISH/CLEAR/PREP7ET,1,SOLID65R,1MP,EX,1,30E3MP,NUXY,1,0.2TB,CONCR,1,1TBTEMP,0TBDATA,1,0.3,0.5,2.5,-1TB,MISO,1,1,15TBPT,, 100E-6, 3.0TBPT,, 300E-6, 8.3TBPT,, 600E-6,14.6TBPT,, 900E-6,19.1TBPT,,1100E-6,21.0TBPT,,1250E-6,22.0TBPT,,1400E-6,22.6TBPT,,1550E-6,22.8TBPT,,1650E-6,22.7TBPT,,1800E-6,22.3TBPT,,2000E-6,21.4TBPT,,2800E-6,16.8 TBPT,,3200E-6,14.7 TBPT,,3800E-6,12.3 TBPT,,4600E-6, 9.9ET,2,LINK8R,2,30R,3,6MP,EX,2,200E3MP,NUXY,2,0.3TB,BKIN,2,1TBTEMP,0TBDATA,,300,0/ESHAPE,1/TRIAD,LBOT/VIEW,1,1,1,1N,1,0,0,0N,2,20,0,0NGEN,11,2,1,2,1,0,20 NGEN,43,22,1,22,1,0,0,20 E,1,2,4,3,23,24,26,25 EGEN,10,2,1TYPE,2REAL,2MAT,2E,3,25E,4,26REAL,3E,19,41E,20,42EGEN,42,22,1,14E,47,63E,48,64EGEN,10,88,589,590/NUMBER,1/PNUM,TYPE,1/REPLOTNSEL,S,LOC,Z,0D,ALL,UZ,0NSEL,S,LOC,Y,0NSEL,R,LOC,Z,740D,ALL,UY,0NSEL,ALLFINISH/SOLUANTYPE,STATICNSEL,S,LOC,Y,200NSEL,R,LOC,Z,200D,ALL,UY,-2TIME,200NSUBST,10LSWRITE,1D,ALL,UY,-0.21TIME,21NSUBST,100LSWRITE,2D,ALL,UY,-1TIME,100NSUBST,100LSWRITE,3NSEL,ALLAUTOTS,1OUTRES,ALL,1PRED,ONLSSOLVE,1,3FINISH11.四层弹簧－质点模型的地震分析包括频率分析，谱分析和时程分析作者：陆新征，清华大学土木工程系FINI/CLEAR/UNITS,SI/PREP7DAMPRATIO=0.02 !振型阻尼系数都为0.02ET,1,BEAM188KEYOPT,1,7,1KEYOPT,1,8,1ET,2,MASS21,,,4MP,EX,1,210E9MP,NUXY,1,0.3MP,DENS,1,7800MPDATA,DAMP,1,,DAMPRATIO SECTYPE,3,BEAM,ISECDATA,0.2,0.2,0.5,0.01,0.01,0.008 R,1,1.6E2,1.6E2R,2,1.2E2,1.2E2K,1,0,0,0K,2,0,3,0K,3,0,6,0K,4,0,9,0K,5,0,12,0K,6,0,0,10000L,1,2L,2,3L,3,4L,4,5LSEL,,,,1,4LATT,1,0,1,,6,,3ALLSLESIZE,ALL, , ,3 , , , , ,1 LMESH,ALLNUMMRG,ALLNUMCMP,ALLTYPE,2REAL,1E,2TYPE,2REAL,2E,8E,14E,20*SET,NT,1001*SET,DT,0.02*DIM,AC,,NT*VREAD,AC(1),RECORD,TXT(F8.3)/SOLUD,1,ALL,ANTYPE,0SOLVEFINISH/SOLU !模态分析ANTYPE,2MODOPT,SUBSP,8MXPAND,8, , ,1SOLVEFINI!得到自振频率1*GET,FREQ1,MODE,1,FREQ/ESHAPE,1.0/SOLUANTYPE,SPECTR !谱分析SPOPT,SPRSMDAMP,1,DAMPRATIO, , , , , ,SED,1SVTYPE,2,1FREQ,0.167,0.25,0.333,0.44444,0.5,0.667,1,1.25,1.667SV,,0.154350625,0.191590625,0.210210625,0.224175625,0.250716714,0.329 514922,0.484352764,0.598723486,0.786897371FREQ,2.222,10,100000SV,,1.034212766,1.034212766,0.3528!按7度多遇地震，地震影响系数0.08 第一组III场地TG＝0.45SRSS,0.01,DISPSOLVEFINISH/POST1/INP,,MCOMNSORT,U,SUM*GET,DMAX_SPE,SORT, ,MAX !反应谱最大变形!开始地震计算FINISH/CONFIG,NRES,20000/SOLUANTYPE,TRANSTRNOPT,FULLALPHAD,2*DAMPRATIO*FREQ1*2*3.1415926BETAD,2*DAMPRATIO/(FREQ1*2*3.1415926)*DO,I,1,1001ACEL,AC(I),0,0TIME,I*0.02OUTRES,ALL,ALLSOLVE*ENDDOFINISH/POST26FILENSOL,2,20,U,X,PLVAR,2, , , , , , , , , ,STORE,MERGEXVAR,1PLVAR,2,*GET,DMAX_TRANS,VARI,2,EXTREM,VMAX, , ! 时程分析最大位移数据文件record.txt-0.0014-0.0108-0.0101-0.0088-0.0095-0.012-0.0142-0.0128-0.011-0.0085-0.0085-0.0131-0.0176-0.0194-0.0162-0.0144-0.0108-0.0082-0.0042-0.0066-0.0131-0.019-0.0196-0.00660.0030.0141-0.0049-0.0144 -0.0203 -0.026 -0.0325 -0.0306 -0.0172 -0.0197 -0.0163 -0.0164 -0.0067 0.0025 0.015 0.0236 0.0252 0.0336 0.0463 0.0492 0.0419 0.0359 0.0271 0.0235 0.0339 0.0412 0.053 0.0639 0.0732 0.0652 0.0599 0.04 0.04 0.0063 -0.0515 -0.0787 -0.0603 -0.0484 -0.025 -0.0059 0.0134 0.0308 0.0499 0.071 0.0995 0.12190.1449 0.1155 0.0935 0.0892 0.0926 0.0839 0.0901 0.0993 0.1209 0.0328 -0.1475 -0.2066 -0.1989 -0.2034 -0.1816 -0.1725 -0.1752 -0.1753 -0.1805 -0.163 -0.1347 -0.1087 -0.0782 -0.0429 -0.0017 0.036 0.0785 0.1164 0.1598 0.196 0.2412 0.2729 0.3036 0.32 0.3417 0.2821 0.2324 -0.1198 -0.2373 -0.164 -0.1865 -0.1095 -0.07530.0113 0.0533 0.0895 0.1186 0.1757 0.0576 -0.2631 -0.1547 -0.1729 -0.1012 -0.0579 0.0237 -0.067 -0.198 -0.1641 -0.1685 -0.1481 -0.1231 -0.1001 -0.0751 -0.0523 -0.0271 -0.0044 0.0188 -0.0095 -0.0433 -0.0838 -0.0951 -0.0716 -0.0599 -0.0334 -0.0108 0.0185 0.042 0.0673 -0.0097 -0.0372 -0.004 0.0011 0.0344 0.0565 0.0883 0.1130.0219 0.0241 0.0683 0.0689 0.1318 0.1353 0.204 -0.0931 -0.1308 -0.0692 -0.0546 0.0072 0.0675 -0.1067 -0.1488 -0.1071 -0.1162 -0.0762 -0.0559 -0.0215 -0.0126 -0.0674 -0.0324 -0.0337 -0.0109 0.0017 0.0299 0.0488 0.0608 0.0222 -0.0032 -0.0245 0.0077 0.0211 0.0568 0.0826 0.1206 0.1478 0.1737 0.0421 0.0029 0.0259 0.0293-0.0147 0.0143 0.0206 0.0499 0.0645 0.0957 0.1128 0.1447 0.1629 0.1945 0.1856 0.1984 0.1769 0.125 -0.1207 -0.0542 -0.0384 -0.0311 -0.1118 -0.1661 -0.2464 -0.2025 -0.1835 -0.1317 -0.096 -0.0325 0.0154 0.0816 0.1319 0.1818 -0.0058 -0.0169 0.0285 0.0447 0.0983 0.1424 0.1853 0.2456 0.1685 -0.138 -0.0999 -0.1089 -0.0907-0.125 -0.2111 -0.1617 -0.1692 -0.1306 -0.1111 -0.0773 -0.051 -0.0544 -0.12 -0.1209 -0.1158 -0.1145 -0.0717 -0.0546 0.0064 -0.0804 -0.1634 -0.0859 -0.0961 -0.0396 -0.0147 0.0319 0.0648 0.0876 0.0472 0.0198 -0.0027 0.0292 0.0445 0.0785 0.1033 0.1352 0.1606 0.1861 0.1281 0.064 0.0204 0.0314 0.0373 0.0496 0.0235 -0.0084-0.0113 -0.0229 -0.0248 -0.0157 -0.0069 0.0147 0.0379 0.0579 0.0255 -0.0041 -0.0428 -0.0133 0.0095 0.023 -0.0129 -0.005 0.008 0.021 0.038 0.051 0.0157 -0.0032 -0.0111 0.0005 0.0076 0.0035 -0.0095 -0.0036 -0.0016 0.0038 0.0085 -0.0056 -0.0304 -0.0421 -0.0244 -0.0236 -0.0177 -0.0129 -0.0018 0.0203 -0.0108 -0.0091 -0.0034-0.0111 -0.0099 -0.0002 0.0073 0.0235 0.0355 0.0705 0.0779 0.0184 -0.0263 -0.0124 -0.0042 0.0159 0.0048 -0.0219 -0.0467 -0.0428 -0.0216 -0.0043 0.0159 0.032 0.0419 0.0123 -0.016 -0.0204 -0.0082 -0.0206 -0.0137 -0.0055 0.0053 0.0134 0.0266 0.0232 0.0079 -0.0008 0.02 0.0435 0.0492 0.0191 0.0092 -0.0022 -0.0021 0.00520.0255 0.0368 0.0525 0.0541 0.0425 0.0398 0.0559 0.0756 0.0365 0.0411 0.0098 -0.0204 -0.0249 -0.0405 -0.0413 -0.0471 -0.0433 -0.0458 -0.0057 0.0178 -0.0208 -0.0492 -0.053 -0.0362 -0.0405 -0.0308 -0.0316 -0.0265 -0.0265 -0.0269 -0.0345 -0.0309 -0.0217 -0.0078 0.0087 0.0281 0.031 0.0358 0.0341 0.0358 0.0287 0.0305 0.01120.0136 0.0384 -0.0861 -0.1349 -0.1342 -0.1354 -0.1193 -0.1042 -0.0829 -0.0651 -0.0444 -0.0258 -0.006 -0.0091 -0.0182 -0.0147 0.0085 0.0163 0.005 0.0264 0.0582 0.0867 0.12 0.1695 0.1111 -0.11 -0.0366 -0.0445 -0.0236 -0.096 -0.0656 -0.0597 -0.067 -0.0552 -0.0027 0.0378 0.1072 0.1669 0.0947 0.0408 0.0667 0.0132 -0.0095-0.0827 -0.1152 -0.115 -0.0803 -0.0369 0.0029 0.0545 0.1178 0.161 -0.027 0.0034 -0.0056 0.002 0.0146 0.0537 0.0798 -0.0205 -0.059 -0.0169 -0.0175 -0.0028 0.0074 0.0382 0.0567 0.0753 0.0801 0.0592 0.0304 0.0023 0.0064 -0.0406 -0.0451 -0.0079 0.0168 0.0567 0.0093 -0.0055 0.0044 -0.0123 -0.0282 -0.0437 -0.0352 -0.02550.0205 0.0519 0.0854 0.1144 0.0733 0.0237 -0.0368 -0.0271 -0.0217 -0.0873 -0.0973 -0.0589 -0.0336 0.0077 0.0259 0.0508 0.0361 0.0081 -0.0056 -0.0209 -0.0317 -0.0238 -0.0376 -0.055 -0.0722 -0.0803 -0.0523 -0.034 -0.0011 0.0065 -0.0037 -0.0005 -0.0168 -0.041 -0.008 0.0079 0.0374 0.0615 0.0665 0.0254 -0.0057 -0.0474 -0.0356-0.0048 0.0126 0.0379 0.0241 -0.0227 -0.0428 -0.0679 -0.0661 -0.059 -0.0513 -0.0408 -0.0309 -0.0266 -0.0541 -0.0628 -0.0908 -0.1107 -0.0881 -0.077 -0.0582 -0.0473 -0.0333 -0.0199 0.002 0.0211 0.0432 0.0613 0.0767 0.0933 0.1066 0.113 0.1187 0.1247 0.1334 0.1594 0.1797 0.2037 0.1236 0.0442 -0.014 -0.0666 -0.0555 -0.0693-0.1246 -0.1179 -0.105 -0.092 -0.0743 -0.0809 -0.085 -0.086 -0.0863 -0.0873 -0.0868 -0.0885 -0.0537 0.0052 0.0215 0.0245 0.058 0.0314 0.0236 0.0485 0.0589 0.0525 0.0355 0.0197 0.0199 0.0492 0.0343 0.0288 0.0432 0.0239 0.0088 0.0077 -0.0148 -0.0077 -0.0019 0.0075 0.0044 -0.0145 -0.0316 -0.0241 -0.0028 0.0182 0.04260.0512 0.0466 0.0479 0.0193 0.0222 0.0274 0.0393 0.0504 0.0577 0.0588 0.0822 0.0797 0.0949 0.0345 0.0045 -0.0123 -0.0347 -0.0426 -0.0416 -0.0275 -0.027 0.0074 0.0428 -0.0231 -0.0387 -0.0083 0.0139 0.0445 0.0027 -0.0697 -0.0796 -0.0251 -0.0135 0.0079 -0.0115 -0.0251 -0.0333 -0.0269 -0.0301 -0.02 -0.0067 -0.0038 0.01050.0344 0.0957 0.0898 0.0179 -0.0362 -0.0994 -0.0807 -0.0744 -0.0539 -0.033 -0.0128 0.0031 0.0148 0.0508 -0.0022 -0.0489 -0.0358 -0.0691 -0.0516 -0.0371 0.0088 0.0632 0.0841 0.1276 0.1388 0.1193 0.0751 0.0225 -0.0088 -0.0227 0.0074 0.0181 0.0544 0.0399 0.0045 -0.0082 -0.0185 -0.002 0.0006 -0.0117 -0.021 -0.0303 -0.0512-0.0579 -0.0266 -0.0178 0.004 0.0098 0.0137 0.0221 0.0437 0.0091 -0.0548 -0.0555 -0.0243 -0.0081 0.025 0.041 0.0182 -0.0027 -0.0243 -0.0015 0.0247 0.0482 0.0783 0.0622 0.0331 -0.0014 -0.0195 -0.0247 -0.0212 -0.011 0.005 0.0241 -0.0034 -0.0216 -0.0471 -0.0363 -0.0195 -0.0018 0.017 -0.008 0.0005 0.023 0.0374 0.06010.0432 0.0344 0.0505 0.0653 0.0683 0.0172 -0.017 -0.0527 -0.0664 -0.0387 -0.0222 -0.0033 0.0119 -0.0128 -0.0351 -0.0514 -0.0335 -0.0218 -0.0012 0.0142 0.007 -0.0063 -0.012 -0.0322 -0.0346 -0.0091 0.0073 0.0309 0.0472 0.0603 0.0576 0.033 -0.0073 -0.0777 -0.0608 -0.0438 -0.0209 0.0031 0.035 0.0293 0.0121 0.0338 0.03170.0206 0.0198 0.0174 0.0021 -0.0144 -0.0343 -0.0339 -0.0145 -0.0028 0.017 -0.0096 -0.0255 -0.0279 -0.0388 -0.0242 -0.0215 -0.0182 -0.0174 -0.0038 -0.0027 -0.0185 -0.0123 0.0087 0.0343 0.0695 0.091 0.0853 0.076 0.0513 0.0186 0.0015 -0.019 -0.0151 -0.0073 0.0021 0.0129 0.0215 0.0024 -0.0124 -0.0329 -0.0519 -0.0708 -0.0579-0.0307 -0.0145 -0.0009 -0.018 -0.0318 -0.0465 -0.0391 -0.0345 -0.0316 -0.0435 -0.0491 -0.0475 -0.042 -0.0361 -0.0277 -0.0258 -0.0139 -0.0068 0.0507 0.0722 0.0878 0.0782 0.0765 0.0439 0.008 0.0013 -0.0126 -0.0015 0.003 0.0104 0.0104 0.0193 0.0205 0.0074 -0.0056 -0.0072 0.007 0.0106 0.0147 -0.0009 -0.0159 -0.0187 -0.00070.0105 -0.0115 -0.0302 -0.0309 -0.0095 -0.0058 0.0004 0.002 0.005 0.0057 0.0097 0.0134 0.0177 0.0218 0.0261 0.0302 0.0346 0.0386 0.0474 0.0393 0.0238 0.0115 -0.0079 -0.0124 0.0054 0.0027 -0.025 -0.0566 -0.063 -0.0591 -0.0413 -0.0068 0.0272 0.0277 -0.0021 -0.006 -0.011 -0.0221 -0.0416 -0.0519 -0.0222 0.003 0.00790.0171 0.0253 0.0323 0.0391 0.0164 -0.0136 -0.0323 -0.0291 -0.0287 -0.0304 -0.0339 -0.0245 -0.0076 0.0125 0.0376 0.0402 0.0245 0.0156 -0.004 -0.0153 -0.0289 -0.0316 -0.0111 0.0094 0.0335 0.0576 0.0424 0.0143 -0.0007 -0.0135 -0.027 -0.0341 -0.0357 -0.0396 -0.0402 -0.0488 -0.048 -0.0406 -0.0407 -0.0351 -0.0187 -0.0057 0.0044-0.0072-0.0169-0.01150.0126finish/clear12.悬臂梁地震分析!定义单元类型/prep7ET,1,BEAM3R,1,0.09,0.000675,0.3, !0.3*0.3m截面MP,EX,1,2.06e11MP,NUXY,1,0.3MP,DAMP,1,0.02 !阻尼比0.02MP,dens,1,7.85e3!--------------------------------------------------!定义节点单元和位移约束*do,i,1,11n,i,(i-1)*0.5,0,0 !5m长,分10段*enddo*do,i,1,10e,i,i+1*enddod,1,ALL,0!--------------------------------------------------!定义和读入时程曲线NT=1000 !时程曲线有NT个点DT=0.01 !时间间隔*dim,ac,,NT !读入数据,这个数据文件可以用excel等软件来写/input,tianjin,txt !天津波东西向共10秒!--------------------------------------------------/SOLUNSUBST,1, , ,1 !1个子步OUTRES,ALL,1 !输出每个子步的结果ANTYPE,TRANS !时程分析*do,i,1,NTACEL,0,ac,0solve*enddo!--------------------------------------------------/POST26NSOL,2,11,U,Y,PLVAR,2,3, , ,下面是文件tianjin.txt，这个文件可以用excel等软件制作或编写很短的程序就能写出来先贴前500个点ac( 1)= -0.06334598ac( 2)= -0.04417088ac( 3)= -0.02188456ac( 4)= -0.00622243ac( 5)= 0.01599961ac( 6)= 0.03431334ac( 7)= 0.05332774ac( 8)= 0.07603510ac( 9)= 0.09322070ac( 10)= 0.11568701ac( 11)= 0.13706695ac( 12)= 0.15446468ac( 13)= 0.17966937ac( 14)= 0.16881546ac( 15)= 0.10252541ac( 16)= 0.03081323ac( 17)= -0.02835768ac( 18)= -0.04577148ac( 19)= 0.00371867ac( 20)= 0.05160179ac( 21)= 0.08542017ac( 22)= 0.14456214ac( 23)= 0.15453540ac( 24)= 0.07974091ac( 25)= -0.00708701ac( 26)= -0.08872422ac( 27)= -0.15415934ac( 28)= -0.16546641ac( 29)= -0.15415934ac( 30)= -0.15538390ac( 31)= -0.15533569ac( 32)= -0.15030567ac( 33)= -0.14931896ac( 34)= -0.14984284ac( 35)= -0.14375542 ac( 36)= -0.14322831 ac( 37)= -0.14130631 ac( 38)= -0.07721145 ac( 39)= 0.02666387 ac( 40)= 0.06092257 ac( 41)= 0.03247169 ac( 42)= 0.01391368 ac( 43)= -0.00049497 ac( 44)= -0.02224453 ac( 45)= -0.04070290 ac( 46)= -0.05431446 ac( 47)= -0.05096541 ac( 48)= -0.03585931 ac( 49)= -0.02533004 ac( 50)= -0.01427366 ac( 51)= -0.00180630 ac( 52)= 0.00995074 ac( 53)= 0.02199384 ac( 54)= 0.03315949 ac( 55)= 0.04807596 ac( 56)= 0.05553902 ac( 57)= 0.03931121 ac( 58)= -0.00056568 ac( 59)= -0.05451052 ac( 60)= -0.07338993 ac( 61)= -0.02214490 ac( 62)= 0.05202605 ac( 63)= 0.11738081 ac( 64)= 0.17217417 ac( 65)= 0.18113498 ac( 66)= 0.15603314 ac( 67)= 0.14236373 ac( 68)= 0.12819293 ac( 69)= 0.10743007 ac( 70)= 0.09416241 ac( 71)= 0.07652363 ac( 72)= 0.05719426 ac( 73)= 0.04319058 ac( 74)= 0.02332124 ac( 75)= 0.00492073 ac( 76)= -0.00940756 ac( 77)= -0.04927159 ac( 78)= -0.13574918。

/config,nres,2000000/units,si !国际单位制fk1=1!法向接触刚度ftoln=0.1!初始渗透因子/prep7et,1,solid45A=3.14*0.3*0.3cylind,0,0.192,0,11,0,90 !旋喷桩上部内侧v1cylind,0.192,0. 2,0,11,0,90!钢管v2cylind,0.2,0.3,0,11,0,90 !旋喷桩上部外侧v3cylind,0,0.3,11,22,0,90 !旋喷桩下部v4cylind,0.3,1.5,0,2.6,0,90 !2.6m淤泥质粉质粘土v5cylind,0.3,1.5,2.6,7.92,0,90 !5.32m粉土v6cylind,0.3,1.5,7.92,16.3,0,90 !8.38m淤泥v7cylind,0.3,1.5,16.3,18.67,0,90!2.37m粉质粘土v8cylind,0.3,1.5, 18.67,20.47 ,0,90!1.8m粉质粘土v9cylind,0.3,1.5, 20.47,25.3 ,0,90!4.83m粉土v10cylind,0,0.3,22,42,0,90 !桩底内侧20m粉质粘土v11cylind,0.3,1.5,25.3,42,0,90 !桩底外侧20m粉质粘土v12allsel,allvsel,s,,,2mp,ex,2,2e11 !钢管mp,nuxy,2,0.3mp,dens,2,7850allsel,allvsel,s,,,1,4vsel,u,,,2mp,ex,3,1.5e8mp,nuxy,3,0.2mp,dens,3,2131!tb,concr,3 !水泥土破坏模型，材料属性!tbdata,,0.7,0.9,0.3e6,2e6!张开裂缝的剪力传递系数,闭合裂缝的剪力传递系数,单轴抗拉强度, 单轴抗压强度allsel,allvsel,s,,,5vplotmp,ex,5,3e6 !淤泥质粉质黏土mp,nuxy,5,0.35mp,dens,5,1950!mp,mu,5,0.633 !设置摩擦系数为0.633tb,dp,5 !土层5的D-P参数tbdata,1,13.2e3,10.1,0 !土的粘聚力，摩擦角，膨胀角allsel,allvsel,s,,,6vplotmp,ex,6,9e6 !粉土mp,nuxy,6,0.3mp,dens,6,1970!mp,mu,6,0.633 !设置摩擦系数为0.633tb,dp,6 !土层6的D-P参数tbdata,1,15.2e3,21.2,30 !土的粘聚力，摩擦角，膨胀角allsel,allvsel,s,,,7vplotmp,ex,7,1.5e6 !淤泥mp,nuxy,7,0.4mp,dens,7,1670!mp,mu,6,0.633 !设置摩擦系数为0.633tb,dp,7 !土层7的D-P参数tbdata,1,12.9e3,9,0 !土的粘聚力，摩擦角，膨胀角allsel,allvsel,s,,,8vplotmp,ex,8,5e6 !粉质粘土mp,nuxy,8,0.35mp,dens,8,1880!mp,mu,8,0.633 !设置摩擦系数为0.633tb,dp,8 !土层8的D-P参数tbdata,1,15.8e3,15.2,0 !土的粘聚力，摩擦角，膨胀角allsel,allvsel,s,,,9vplotmp,ex,9,8e6 !粉质粘土mp,nuxy,9,0.35mp,dens,9,2000!mp,mu,9,0.633 !设置摩擦系数为0.633tb,dp,9 !土层9的D-P参数tbdata,1,37.1e3,19.8,0 !土的粘聚力，摩擦角，膨胀角allsel,allvsel,s,,,10mp,ex,10,10e6 !粉土mp,nuxy,10,0.3mp,dens,10,1980!mp,mu,10,0.633 !设置摩擦系数为0.633tb,dp,10 !土层10的D-P参数tbdata,1,14.3e3,20.1,10 !土的粘聚力，摩擦角，膨胀角allsel,allvsel,s,,,11,12vplotmp,ex,11,6e6 !粉质粘土mp,nuxy,11,0.35mp,dens,11,1910!mp,mu,10,0.633 !设置摩擦系数为0.633tb,dp,11 !土层11的D-P参数tbdata,1,23.4e3,15.9,0 !土的粘聚力，摩擦角，膨胀角!布尔运算粘结体所有的除桩之外的部分ALLSEL,ALLvsel,s,,,5,12,1vplotvglue,allvsel,s,,,1,4vplotvglue,allallsel,all/pnum,line,1saveallsel,allallsel,alllsel,s,,,3,4lsel,a,,,13,14lesize,all,,,6lsel,s,,,7,8lsel,a,,,18,19lesize,all,0.1lsel,s,,,55,58lesize,all,,,1type,1mat,2 !钢管vmesh,6vsel,s,,,1,7,6vsel,a,,,8vplottype,1mat,3 !水泥土lsel,s,,,1,2lesize,all,,,6lsel,s,,,55,56 lesize,all,,,1lsel,s,,,22,24,2 lesize,all,,,10,0.8,1 lsel,s,,,3lesize,all,,,6lsel,s,,,13lesize,13,,,6lsel,s,,,25lesize,25,,,6lsel,s,,,7,9lsel,a,,,18,19lsel,a,,,30,31lsel,a,,,67,69 lesize,all,0.1lsel,s,,,4,6 !体8 lsel,a,,,14,26,12 lesize,all,,,6lsel,s,,,57,58 lesize,all,,,1lsel,s,,,65,66 lesize,all,,,10,0.8,1 aadd,35,31,2 aslv,slsla,sFLST,2,3,4,ORDE,3 FITEM,2,6 FITEM,2,57 FITEM,2,65 LCCAT,P51X FLST,2,3,4,ORDE,3 FITEM,2,5 FITEM,2,58 FITEM,2,66 LCCAT,P51X vsweep,1vsweep,8 allsel,allvsel,u,,,1,6,5 vsel,u,,,7,8 vsel,u,,,18vplotlsel,s,,,44,46,2 lesize,all,,,6 lsel,s,,,43,45,2 lesize,all,,,5,0.8,1 lsel,s,,,51,54, lesize,all,0.1 lsel,s,,,47,49,2 lesize,all,,,6 lsel,s,,,48,50,2 lesize,all,,,5,0.8,1 type,1mat,5 vsweep,5 vsel,s,,,13 lsel,s,,,136,139 lesize,all,0.1 lsel,s,,,59,61,2 lesize,all,,,6 lsel,s,,,60,62,2 lesize,all,,,5,0.8,1 type,1mat,6 vsweep,13 vsel,s,,,14 lsel,s,,,140,143 lesize,all,0.1 lsel,s,,,71,73,2 lesize,all,,,6 lsel,s,,,72,74,2 lesize,all,,,5,0.8,1 type,1mat,7 vsweep,14lsel,s,,,144,147 lesize,all,0.1 lsel,s,,,83,85,2 lesize,all,,,6 lsel,s,,,84,86,2 lesize,all,,,5,0.8,1 type,1mat,8 vsweep,15 vsel,s,,,16lsel,s,,,148,151 lesize,all,0.1 lsel,s,,,95,97,2 lesize,all,,,6 lsel,s,,,96,98,2 lesize,all,,,5,0.8,1 type,1mat,9 vsweep,16 vsel,s,,,19lsel,s,,,158,163 lesize,all,0.1 lsel,s,,,107,109,2 lesize,all,,,6 lsel,s,,,108,110,2 lesize,all,,,5,0.8,1 type,1mat,10 vsweep,19 vsel,s,,,17lsel,s,,,152,156 lesize,all,0.1 lsel,s,,,118,128,10 lesize,all,,,6 lsel,s,,,155,157,2 lesize,all,,,5,0.8,1 type,1mat,11 vsweep,17lsel,s,,,154,158,2lsel,a,,,159lsel,a,,,123lesize,all,0.1lsel,s,,,118,120,lesize,all,,,6lsel,s,,,115,117lesize,all,,,6type,1mat,11vsweep,18!用接触向导建立接触对et,2,170et,3,173keyopt,3,9,0! 包括由几何形状和接触表面偏移所引起的穿透keyopt,3,12,2 !不分开接触，允许滑动，法向不分开!桩建立目标面r,20 ,,,fk1,ftoln, , ,mp,mu,100,0.85allsel,allvsel,s,,,1,6,5vsel,a,,,7,8aslv,saplotasel,s,,,14asel,a,,,19,39,20aplotcm,target,areatype,2real,20mat,100nsla,s,1esurf,allesel,s,type,,2eplotr,21,,,fk1,ftoln, , ,mp,mu,101,0.3allsel,allvsel,s,,,1,6,5aslv,saplotasel,s,,,14 asel,a,,,19,39,20 aplotcm,target,area type,2real,21mat,101nsla,s,1 esurf,allesel,s,type,,2 eplotr,22,,,fk1,ftoln, , , mp,mu,102,0.15 allsel,allvsel,s,,,1,6,5 vsel,a,,,7,8 aslv,saplotasel,s,,,14 asel,a,,,19,39,20 aplotcm,target,area type,2real,22mat,102nsla,s,1 esurf,allesel,s,type,,2 eplotr,23,,,fk1,ftoln, , , mp,mu,103,0.3 allsel,allvsel,s,,,1,6,5 vsel,a,,,7,8 aslv,saplotasel,s,,,14 asel,a,,,19,39,20 aplotcm,target,area type,2real,23mat,103nsla,s,1 esurf,allesel,s,type,,2 eplotr,24,,,fk1,ftoln, , , mp,mu,104,0.35 allsel,allvsel,s,,,1,6,5 vsel,a,,,7,8 aslv,saplotasel,s,,,14 asel,a,,,19,39,20 aplotcm,target,area type,2real,24mat,104nsla,s,1 esurf,allesel,s,type,,2 eplotr,25,,,fk1,ftoln, , , mp,mu,105,0.3 allsel,allvsel,s,,,1,6,5 vsel,a,,,7,8 aslv,saplotasel,s,,,14 asel,a,,,19,39,20 aplotcm,target,area type,2real,25mat,105nsla,s,1 esurf,allesel,s,type,,2 eplot!桩建立接触面allsel,allvsel,u,,,1,6,5 vsel,u,,,7,8 vplotasel,s,,,26 asel,a,,,72 asel,a,,,76 asel,a,,,80 asel,a,,,84 asel,a,,,95 asel,a,,,59 aplottype,3real,20mat,100 nsla,s,1 esurf,all allsel,allvsel,u,,,1,6,5 vsel,u,,,7,8 vplotasel,s,,,72 aplottype,3real,21mp,mu,101,0.3 mat,101 nsla,s,1 esurf,all allsel,allvsel,u,,,1,6,5 vsel,u,,,7,8 vplotasel,s,,,76 aplottype,3real,22mp,mu,102,0.15mat,102 nsla,s,1 esurf,all allsel,allvsel,u,,,1,6,5 vsel,u,,,7,8 vplotasel,s,,,80 aplottype,3real,23mp,mu,103,0.3 mat,103 nsla,s,1 esurf,all allsel,allvsel,u,,,1,6,5 vsel,u,,,7,8 vplotasel,s,,,84 aplottype,3real,24mp,mu,104,0.35 mat,104 nsla,s,1 esurf,all allsel,allvsel,u,,,1,6,5 vsel,u,,,7,8 vplotasel,s,,,59,95,36 aplottype,3real,25mp,mu,105,0.3 mat,105 nsla,s,1 esurf,all/psymb,esys,1allselgplotsavefinish!求解过程/solu!定义面约束条件asel,s,loc,x,0 !侧面施加对称约束da,all,symmasel,s,loc,y,0da,all,symmasel,s,loc,z,42 !底面试加约束da,all,allallsel,all !土外层施加约束asel,s,,,86asel,a,,,94asel,a,,,82asel,a,,,78asel,a,,,74asel,a,,,70asel,a,,,25aplotda,all,allALLSEL,ALLantype,staticnlgeom,onsolcontrol,onautot,onpred,onneqit,50!每个子步允许的最大迭代次数outpr,all,allOUTRES,all,allFirstF=840e3/ALastF=840e3/A!考虑注浆会增加30％的部分Increment=60e3/A*do,Force,FirstF,LastF,Incrementallsel,allvsel,s,,,1vsel,a,,,6vsel,a,,,7aslv,s,asel,s,,,1,asel,a,,,29asel,a,,,34aplotsfa,all,,pres,Forceallsel,allnsubst,10,50,5solve*enddofinish【下载本文档，可以自由复制内容或自由编辑修改内容，更多精彩文章，期待你的好评和关注，我将一如既往为您服务】。

第三章回弹/filname,huitan/units,si/prep7et,1,plane42keyopt,1,3,2mp,ex,1,2.160687961e11 mp,nuxy,1,0.3TB,MISO,1,1,25, TBTEMP,0TBPT,,0.000814,175.88e6, TBPT,,0.001531,324.98e6, TBPT,,0.002819,591.26e6, TBPT,,0.003846,801.457e6, TBPT,,0.004239,873.22e6, TBPT,,0.004594,924.399e6, TBPT,,0.004955,956.88e6, TBPT,,0.00531,976.99e6, TBPT,,0.005699,990.2138e6, TBPT,,0.006261,1000e6, TBPT,,0.006957,1008.34e6, TBPT,,0.007804,1015.5e6, TBPT,,0.008779,1021.185e6, TBPT,,0.009827,1026.47e6, TBPT,,0.010929,1031.289e6, TBPT,,0.012066,1035.444e6, TBPT,,0.0132,1039.693e6, TBPT,,0.014354,1044.886e6, TBPT,,0.015494,1049.796e6, TBPT,,0.016631,1053.951e6, TBPT,,0.017785,1058.106e6, TBPT,,0.019039,1062.827e6, TBPT,,0.020276,1067.643e6, csys,1k,1,0.00512k,2,0.0064kgen,36,1,2,1,,10*do,j,1,69,2a,j,j+1,j+3,j+2*enddoa,71,72,2,1savelesize,all,,,4amesh,allnummrg,allnumcmp,allsavecp,1,uy,166,170,171,172,173finish/solunsel,s,loc,y,270d,all,allnplotallselnplotnsel,s,loc,y,90d,all,uxnplotallseloutres,all,alltime,1NSUBST,100nlgeom,ond,166,uy,-0.00128solve*get,def1,node,166,U,Ynsel,s,loc,y,90ddele,all,uyallself,166,fy,0solve*get,def2,node,166,U,Yfinish/post1 !进入后处理器set,list !显示所有计算的结果set,1,last,1 !读取第一个加载步最后子步的结果pldisp,1 !同时显示变形前与变形后结构的形状plnsol,s,eqv,1,1 !显示节点的等效应力set,2,last,1 !读取第二个加载步最后子步的结果pldisp,1 !同时显示变形前与变形后结构的形状。

ANSYS命令流使用方法（中文）修改Finish(退出四大模块，回到BEGIN层)/clear (清空内存，开始新的计算)1．定义参数、数组，并赋值.2．/prep7(进入前处理)定义几何图形：关键点、线、面、体定义几个所关心的节点，以备后处理时调用节点号。

设材料线弹性、非线性特性设置单元类型及相应KEYOPT设置实常数设置网格划分，划分网格根据需要耦合某些节点自由度定义单元表3．/solu加边界条件设置求解选项定义载荷步求解载荷步4./post1（通用后处理）5./post26 （时间历程后处理）6.PLOTCONTROL菜单命令7.参数化设计语言8.理论手册Finish(退出四大模块，回到BEGIN层)/clear (清空内存，开始新的计算)1.定义参数、数组，并赋值.dim, par, type, imax, jmax, kmax, var1, vae2, var3 定义数组par: 数组名type：array 数组，如同fortran,下标最小号为1，可以多达三维（缺省）char 字符串组（每个元素最多8个字符）tableimax,jmax, kmax 各维的最大下标号var1,var2,var3 各维变量名，缺省为row,column,plane(当type 为table时) 2./prep7(进入前处理)2.1 设置单元类型及相应KEYOPTET, itype, ename, kop1……kop6, inopr 设定当前单元类型Itype：单元号Ename：单元名设置实常数Keyopt, itype, knum, valueitype: 已定义的单元类型号knum: 单元的关键字号value: 数值注意：如果,则必须使用keyopt命令，否则也可在ET命令中输入2.2 定义几个所关心的节点，以备后处理时调用节点号。

n,node,x,y,z,thxy, thyz, thzx 根据坐标定义节点号如果已有此节点，则原节点被重新定义，一般为最大节点号。

prep7/pnum,label,key!在有限元模块图形中显示号码。

Label＝欲显示对象的名称，node节点，elem元素，kp点，line线，area面积，volu体积；key=0为不显示号码（系统默认），＝1为显示号码。

et,itype,ename,kopt1, kopt2, kopt3, kopt4, kopt5, kopt6,inopr!元素类型定义。

Itype为元素类型号码，通常由1开始；ename为ANSYS元素库的名称，如beam3,plane42,solid45等；kopt1～kopt6为元素特性编码，如beam3的kopt6＝1时，表示分析后的结果可输出节点的力及力矩，link1无需任何元素特性编码。

mp,lab,mat,c0,c1,c2,c3,c4！定义材料特性。

Lab为材料特性类别，如杨氏系数lab=ex、ey、ez，密度lab=dens，泊松比lab=nuxy、nuyz、nuzx，剪力模数lab=gxy、gyz、gxz，热膨胀系数lab=alpx、alpy、alpz，热传导系数lab=kxx、kyy、kzz，比热lab=c,摩擦系数lab=mu, 参考温度reft；mat对应前面定义的元素类型号码Itype；c0为材料特性类别的值; c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项, 4次项的系数。

r,nset,r1,r2,r3,r4,r5,r6！元素几何特性。

nset通常由1开始；r1～r6几何特性的值。

注：solid45元素不需要此命令，beam3单元有area截面积，惯性矩izz，高度height等。

例如：r,1,3e-4（截面积）,2.5e-9（惯性矩）,0.01（高度）local,kcn,kcs,xc,yc,zc,thxy,thyz,thzx,par1,par2！定义区域坐标系统。

kcn区域坐标系统代号（大于10）；kcs区域坐标系统属性（0为卡式坐标，1为圆柱坐标，2为球面坐标）；xc,yc,zc（该区域坐标系统与整体坐标系统原点关系）。