ANSYS应用——m法计算桩基时模拟桩与土作用的命令流
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法是一种常用于估算桩基承载力的分析方法。
用ANSYS软件分析桩土相互作用
3.2.2 建模并划分网格 本模型简化为三维空间实体结构，桩采用钢管桩实际尺寸，土体取长为20m，宽20m（相当于10倍桩
径，以减小边界效应），深100.58m 的长方体块来模拟。
126
第十三届中国海洋（岸）工程学术讨论会论文集
用 ANSYS 软件分析桩土相互作用
陈建民，肖 花
（中国石油大学石油工程学院，山东 东营 257061）
摘要：用 ANSYS 程序对导管架平台的桩与土之间的相互作用进行了研究。结果表明，在平台载荷作用下，桩上的最大应
力点出现在距海底一定深度的范围内，土中的最大应力也出现在离海底一定距离的土层中，而不是在海底处。桩周土层中应 力分布显著不对称，一侧的土层受力大，另一侧土层的受力小。随着深度的增加，桩身水平方向位移呈衰减波的形式变化。 衰减波的衰减速度很快，波峰的绝对值在第二个峰处就已减小到 0.001m。
159.91
159.91
1.44
13
Very dense sandy SILT
—82.75
9.60
Nq=20
81.39
4.79
14
Very stiff lean CLAY
qmax=4.79 MPa
—88.85
9.80
169.97
169.97
1.53
15
Very stiff lean CLAY
—92.66
刚性桩、弹性桩、弹性长桩的判断标准，与桩打入土中的深度 h 有关。桩打入土中的深度 h 同相对刚
度 T 的比值 Z max 称为相对桩长， Z max 可用于判断一桩为刚性桩、弹性桩或弹性长桩的判据。实验和分析
ansys桩基础算例
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⾥做地震分析时，需要对结构施加地震惯性荷载，地震惯性⼒是通过加速度的⽅式输⼊进结构的，然后与结构的质量⼀起形成动⼒计算时的惯性荷载，下⾯说⼀下在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命令流大全
《史上最全的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命令流总结(全)
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法计算桩基时模拟桩与土作用的命令流
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命令流大全（整理）
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为半径。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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。