当前位置:首页 >> 机械/仪表 >>

ANSYS+APDL+命令集合


ANSYS APDL 命令集合(2) 2007-11-23 16:50:04 /post1 中的几个命令: set, lstep, sbstep, fact, king, time, angle, nset 设定从结果文件 读入的数据 lstep :荷载步数 sbstep:子步数,缺省为最后一步 time: 时间点(如果弧长法则不用) nset: data set number

dscale, wn, dmult 显示变形比例 wn: 窗口号(或 all),缺省为 1 dmult, 0 或 auto : 自动将最大变形图画为构件长的 5% pldisp, kund 显示变形的结构 kund: 0 仅显示变形后的结构 1 显示变形前和变形后的结构 2 显示变形结构和未变形结构的边缘

PRETAB, LAB1, LAB2, ??LAB9 沿线单元长度方向绘单元表数据 LABn : 空: 所有 ETABLE 命令指定的列名 列名: 任何 ETABLE 命令指定的列名 PLLS, LABI, LABJ, FACT, KUND 沿线单元长度方向绘单元表数据 LABI:节点 I 的单元表列名 LABJ:节点 J 的单元表列名 FACT: 显示比例,缺省为 1 kund: 0 不显示未变形的结构 1 变形和未变形重叠 2 变形轮廓和未变形边缘

etable, lab,item,comp 将单元的某项结果制作成表格,以供 pretable 命令输出, lab: 字段名称,自己指定 item: 结果的顶目名称,在每个单元的说明中有(在单元说明表中冒号左 边的 comp, 结果项目名称的分量,在单元说明表中冒号右边的 比如将 plane42 单元的 x 应力分量制成表 etable,sx,x,x LACAL,KCN,KCS,XC,YC,ZC,THXY,THYZ,THZX,PAR1,PAR2 定义区域坐标系统,该命令执行后,ANSYS 坐标系统自动更改为新建立的 坐标系统,故可以定义许多区域坐标系统,以辅助有限元模型的建立。

KCN:该区域坐标系统的确定代号,大于 10 的任何一个号码都可以。 KCS:该区域坐标系统的属性。0,1,2 分别代表卡式坐标,圆柱坐标,球面坐 标。 XC,YC,ZC:该区域坐标系统与整体坐标系统原点的关系。 THXY,THYZ,THZX:该区域坐标系统与整体系统 X,Y,Z 轴的关系 claer,nl1,nl2,lmesh 就是将后面的直线网格化之后的节点和元素都删除 但是共享节点依然存在

mshkey,key 声明是使用自由化网格(key=0) 对应网格(key=1) 或者是混合网格(key=2) 后面两种我因为是新手,所以不大会用,一般都用自由网格~~

关于工作平面: KWPAVE, P1, P2, P3, P4, P5, P6, P7, P8, P9 把工作平面的中心移动到以上几点的平均点 最多 9 如果只选一点,那么就是把工作平面的中心移动到此点

WPOFF, XOFF, YOFF, ZOFF 移动工作平面,注意 xoff,yoff,zoff 是相对当前点的移动量 而不是整体坐标

WPROT, THXY, THYZ, THZX 旋转工作平面 和上面的一样, 是相对当前的工作平面选择一个角度,默认设置是角度为 单位 wpstyl 关闭工作平面显示

Nummrg,label,toler, Gtoler,action,switch 合并相同位置的 item label: 要合并的项目 node: 节点, Elem,单元,kp: 关键点(也合并线,面及点) mat: 材料,type: 单元类型,Real: 实常数 cp:耦合项,CE:约束项,CE: 约束方程,All:所有项 toler: 公差

Gtoler:实体公差 Action: sele 仅选择不合并 空 合并 switch: 较低号还是较高号被保留(low, high) 注意:可以先选择一部分项目,再执行合并。如果多次发生合并命令,一 定要先合并节点, 再合并关键点。 合并节点后, 实体荷载不能转化到单元, 此时可合并关键点解决问题。 我也感觉和 Glue 效果一样,但是它有独到的好处的。 numcmp 是压缩编号,对计算没有影响的。 一个条件命令 *if,val1, oper, val2, base: 条件语句 val1, val2: 待比较的值(也可是字符,用引号括起来) oper: 逻辑操作(当实数比较时,误差为 1e-10) eq, ne, lt, gt, le, ge, ablt, abgt base: 当 oper 结果为逻辑真时的行为 lable: 用户定义的行标志 stop: 将跳出 anasys exit: 跳出当前的 do 循环 cycle: 跳至当前 do 循环的末尾 then: 构成 if-then-else 结构

一个循环命令 *do, par, ival, fval, inc 定义一个 do 循环的开始 par: 循环控制变量 ival, fval, inc:分别为起始值,终值,步长(可正可负) ?? *enddo 定义一个 do 循环的结束

一个网格划分命令 用 desize 定义单元尺寸时单元划分应遵守的级别 高:lesize kesize esize desize 用 smartzing 定义单元尺寸时单元划分应遵守的级别 高:lesize kesize smartsize

定义表、数组等真的很好用哦

*dim, par, type, imax, jmax, kmax, var1, vae2, var3 定义数组 par: 数组名 type: array 数组,如同 fortran,下标最小号为 1,可以多达三维 char 字符串组(每个元素最多 8 个字符) table 表 imax,jmax, kmax 各维的最大下标号 var1,var2,var3 各维变量名,缺省为 row,column,plane 在 ANSYS 帮助系统中关于*SET 命令的注释下列出了 ANSYS 中可以使用的 数学函数。 所有这些数学函数均可以在 ANSYS 环境中使用,这些数学函数 包括: ABS(X) 求绝对值 ACOS(X) 反余弦 ASIN(X) 反正弦 ATAN(X) 反正切 ATAN2(X,Y) 反正切, ArcTangent of (Y/X) , 可以考虑变量 X,Y 的符号 COS(X) 求余弦 COSH(X) 双曲余弦 EXP(X) 指数函数 GDIS(X,Y) 求以 X 为均值,Y 为标准差的高斯分布,在使用蒙地卡罗法研究 随机荷载和随机材料参数时,可以用该函数处理计算结果 LOG(X) 自然对数 LOG10(X) 常用对数(以 10 为基) MOD(X,Y) 求 X/Y 的余数. 如果 Y=0, 函数值为 0 NINT(X) 求最近的整数 RAND(X,Y) 取随机数,其中 X 是下限, Y 是上限 SIGN(X,Y) 取 X 的绝对值并赋予 Y 的符号. Y>=0, 函数值为|X|, Y<0, 函 数值为-|X|,. SIN(X) 正弦 SINH(X) 双曲正弦 SQRT(X) 平方根 TAN(X) 正切 TANH(X) 双曲正切 esel,s,mat,,1 选择材料号为 1 的单元 *get,emin,elem,,num,min 获得最小的单元号 *get,emax,elem,,num,max 获得最大的单元号 *DO,I,emin,emax 作循环 *GET,V1,ELEM,I,VOLU 获得单元的体积存到 V1 的变量中 V=V+V1 求和获得材料 1 的总体积 *enddo

把一个矩阵的一列加起来的方法

提取当前选择集中的结点总数存入变量 aaa1; 提取当前选择集中的结点的最小结点号存入变量 aaa2; 定义 aaa1×2 数组 aaa3; 开始循环: aaa3 数组的第一列存储结点号; aaa3 数组的第二列存储 Sx; 下一个结点号存入变量 aaa2; 循环结束。 /post1 *get,aaa1,node,0,count *get,aaa2,node,0,num,min *dim,aaa3,array,aaa1,2 *do,i,1,aaa1 aaa3(i,1)=aaa2 *get,aaa3(i,2),node,aaa2,s,x aaa2=ndnext(aaa2) *enddo L, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2:在两个关键 点之间定义一条线。 功能:在当前激活坐标系统下,在两个指定关键点之间生成直线或曲线。 P1,P2:线的起点和终点。 NDIV:这条线的单元划分数。一般不用,指定单元划分数推荐用 LESIZE。 这里需要说明一下:如果你的模型相对规则,为了得到高质量的网格,不 妨在划线的时候指定单元划分数,这样,既方便又能按照自己的意愿来分 网。 SPACE:间隔比。通常不用,指定间隔比推荐使用命令 LESIZE。 说明: 线的形状由激活坐标系决定,直角坐标系中将产生一条直线,柱 坐标系中,随关键的坐标不同可能产生直线,圆弧线或螺旋线。 KGEN, ITIME, NP1, NP2, NINC, DX, DY, DZ, KINC, NOELEM, IMOVE:通 过一组关键点生成额外的关键点。 ITIME:生成操作总共执行的次数,如果要生成额外的点,该值必须大于 1。 NP1,NP2,NINC: 被生成的那组关键点的编号为 NP1 至 NP2,编号增量为 NINC (缺省为 1)。 DX,DY,DZ: 关键点在激活坐标系下的位置增量(柱坐标系和球面坐标系下 要注意坐标的变换)。 KINC:生成的点集与原始点集之间的增量值。如果是 0,则指定为最低可 用关键点编号。 NOELEM:指定是否单元和节点也随之生成。0,生成;1,不生成。 IMOVE:指定关键点是否被移除或重新定义。0,按照 ITIME 要求生成额外 关键点;1,移除原始关键点到新的位置,保持编号不变 (ITIME,KINC,NOELEM 被忽略)。

ANSYS 的命令流里经常看到一些相对比较固定的代码,这些代码组合在一 起构成 ANSYS 的一个操作。比如,通过旋转命令将面生成体的操作,命 令流如下: TYPE,2 !指定生成体的单元类型 EXTOPT,ESIZE,18,0, !指定单元划分数 EXTOPT,ACLEAR,1 !清除面网格 VROTAT,ALL,,,,,,15,16, !绕关键点 15,16 构成的轴将所划的面网格旋 转 360 这一组代码在 GUI 方式下通常由一个或几个对话框组成, 还有的命令有时 需要同时使用,否则就会出错。例如:想在自己定义的坐标系下编辑有限 元模型,需要定义坐标系并激活,然后将所有节点移到当前坐标系中。命 令流如下: WPRO,,,90 !将工作平面绕 Y 轴旋转 90 度 CSWPLA,11,1,1,1, !在工作平面原点创建柱坐标系,并激活 NROTAT,ALL !将所有节点旋转到激活坐标系 为了能充分发挥命令流的优势,建议各位把自己常用的代码贡献出来,这 样, 不仅熟悉了命令, 更重要的是, 掌握了一种通过命令流来实现的 ANSYS 操作。 ! DP 材料参数 tb,dp,2 tbdata,,50.0e3,16.7,16.7! ! mp,ex,2,40.0e6 !粉喷桩复合地基特性 mp,dens,2,1770 mp,nuxy,2,0.38 ! ! DP 材料参数 !基座粗砂垫层,20cm tb,dp,3 tbdata,,0.0,25.0,25.0 ! mp,ex,3,46e6 mp,dens,3,1950 mp,nuxy,3,0.30

1.MP,Lab,MAT,C0,C1,C2,C3,C4 定义材料的属性(Material Property),材料属性为固定值时,其值为 C0,当随温度变化时,由后四个参数控制。 MAT:对应 ET 所定义的号码(ITYPE),表示该组属性属于 ITYPE。 Lab:材料属性类别, 任何元素具备何种属性在元素属性表中均有说明。例 如杨氏系数(Lab=EX,EY,EZ),密度(Lab=DENS),泊松比 (Lab=NUXY,NUXYZ,NUZX),剪切模数(Lab=GXY,GYZ,GXZ),热膨胀系数 (Lab=ALPX,ALPY,ALPZ)等。 2./ANTYPE,Antype,Status

声明分析类型,即欲进行哪种分析,系统默认为静力学分析。 Antype=STATIC or 0 静态分析(系统默认) BUCKLE or 1 屈曲分析 MODAL or 2 振动模态分析 HARMIC or 3 调和外力动和系统 TRANS or 4 瞬时动力系统分析 3.SFBEAM, ELEM, LKEY, Lab,VALI, VALJ, VAL2I, VAL2J, IOFFST, JOFFST 定义在梁元素上的分布力。 ELEM:元素号码。 LKEY:建立元素后,依节点顺序梁元素有四个面,该参为分力所施加的面 号。 Lab RES(表示分布压力)。 VALI,VALJ:在 I 点及 J 点分布力的值。 4./pnum,label,key !在有限元模块图形中显示号码。Label=欲显示对象的名称,node 节点, elem 元素,kp 点,line 线,area 面积,volu 体积;key=0 为不显示号码 (系统默认),=1 为显示号码。 5.lesize,nl1,size,angsiz,ndiv,space,kforc,layer1,layer2 ! 定义所 选择线段(nl1,nl1=all 为目前所有的线段)进行元素网格化时元素的大 小(size),元素的大小可用线段的长度(size)或该条线段要分割的元素数 目(ndiv)来确定;space 为间距比(最后一段长与最先一段长的比值,正 值代表以线段方向为基准,负值以中央为基准,系统默认等间距)。 6.plnsol,item,comp !图标节点的解答。以连续的轮廓线表示。 Item 为欲查看何种解答。 Item comp S x,y,z,xy,yz,xz 应力 S 1,2,3 主应力 S eqv,int 等效应力 F x,y,z 结构力 M x,y,z 结构力矩 u x,y,z,sum 位移分量及向量位移 rot x,y,z,sum 旋转位移分量及向量旋转位移 temp 温度

1./UNITS,LABEL 声明单位系统,表示分析时所用的单位,LABEL 表示系统单位,如下所示 LABEL=SI (公制,公尺、公斤、秒) LABEL=CSG (公制,公分、公克、秒) LABEL=BFT (英制,长度=ft) LABEL=BIN (英制,长度=in) 2.节点定义 有限元模型的建立是将机械结构转换为多节点和元素相连接, 所以节点即 为机械结构中一个点的坐标,指定一个号码和坐标位置。在 ANSYS 中所建 立的对象(坐标系、节点、点、线、面、体积等)都有编号。 相关命令 N,NODE,X,Y,Z,THXY,THYZ,THZX 定义节点,若在圆柱坐标系统下 x,y,z 对应 r,θ ,z,在球面系统下对应 r,θ ,?。

NODE:欲建立节点的号码 X,Y,Z:节点在目前坐标系统下的坐标位置 Menu Paths:Main Menu>Preprocessor>Create>Node>In Active CS Menu Paths Main Menu>Preprocessor>Create>Node>On Working Plane NDELE,NODE1,NODE2,NINC 删除在序号在 NODE1 号 NODE2 间隔为 NINC 的所有节点,但若节点已连成 元素,要删除节点必先删除元素。例如: NDELE,1,100,1 !删除从 1 到 100 的所有点 NDELE,1,100,99 !删除 1 和 100 两个点 Menu Paths:Main Menu>Preprocessor>Delete>Nodes NPLOT,KNUM 节点显示, 该命令是将现有卡式坐标系统下节点显示在图形窗口中,以供 使用者参考及查看模块的建立。建构模块的显示为软件的重要功能之一, 以检查建立的对象是否正确。有限元型的建立程中,经常会检查 各个对 象的正确性及相关位置,包含对象视角、对象号码等,所以图形显示为有 限元模型建立过程中不可缺少的步骤。KNUM=0 不显示号码,为 1 显示同 时显示节点号 Menu Paths:Utility Menu>plot>nodes Menu Paths:Utility Menu>plot>Numbering?(选中 NODE 选项) NLIST,NODE1,NODE2,NINC,Lcoord,SORT1,SORT2,SORT3 节点列式, 该命令将现有卡式坐标系统下节点的资料列示于窗口中(会打 开一个新的窗口),使用者可检查建立的坐标点是否正确,并可将资料保 存为一个文件。 如欲在其它坐标系统下显示节点资料,可以先行改变显示 系统,例如圆柱坐标系统,执行命令 DSYS,1。 Menu Paths:Utility Menu>List>Nodes FILL,NODE1,NODE2,NFILL,NSTRT,NINC,ITIME,INC,SPACE 节点的填充命令是自动将两节点在现有的坐标系统下填充许多点, 两节点 间填充的节点个数及分布状态视其参数而定,系统的设定为均分填满。 NODE1,NODE2 为欲填充点的起始节点号码及终结节点号码,例如两节点号 码为 1(NODE1)和 5(NODE2),则平均填充三个节点(2,3,4)介于节 点 1 和 5 之间。 Menu Paths:Main Menu>Preprocessor>Create>Node>Fill between Nds NGEN,ITIME,INC,NODE1,NODE2,NINC,DX,DY,DZ,SPACE 节点复制命令是将一组节点在现有坐标系统下复制到其它位置。 ITIME: 复制的次数,包含自己本身。 INC: 每次复制节点时节点号码的增加量。 NODE1,NODE2,NINC: 选取要复制的节点,即要对哪些节点进行复制。 DX,DY,DZ: 每次复制时在现有坐标系统下,几何位置的改变量。 Menu Paths:Main Menu>Preprocessor>(-Modeling-)Copy>(-Nodes-)Copy ET,ITYPE,Ename,KOPT1,KOPT2,KOPT3,KOPT4,KOPT5,KOPT6,INOPR

元素类型(Element Type)为机械结构系统的含的元素类型种类,例如桌 子可由桌面平面单元各桌脚梁单元构成,故有两个元素类型。ET 命令是 由 ANSYS 元素库中选择某个元素并定义该结构分析所使用的元素类型号 码。 ITYPE:元素类型的号码 Ename:ANSYS 元素库的名称,即使用者所选择的元素。 KOPT1~KOPT6:元素特性编码。 Menu Paths:Main Menu>Preprocessor Element Type>Add/Edit/Delete MP,Lab,MAT,C0,C1,C2,C3,C4 定义材料的属性(Material Property),材料属性为固定值时,其值为 C0,当随温度变化时,由后四个参数控制。 MAT:对应 ET 所定义的号码(ITYPE),表示该组属性属于 ITYPE。 Lab:材料属性类别, 任何元素具备何种属性在元素属性表中均有说明。例 如杨氏系数(Lab=EX,EY,EZ),密度(Lab=DENS),泊松比 (Lab=NUXY,NUXYZ,NUZX),剪切模数(Lab=GXY,GYZ,GXZ),热膨胀系数 (Lab=ALPX,ALPY,ALPZ)等。 Menu paths:Main Menu>Preprocessor>Matial Props>Isotropic R,NSET,R1,R2,R3,R4,R5,R6 定义”实常数”,即某一单元的补充几何特征,如梁单元的面积,壳单元 的厚度。所带的的参数必须与元素表的顺序一致。 Menu paths:Main Menu>Preprocessor>Real Constants E,I,J,K,L,M,N,O,P SOLU 进入解题处理器,当有限元模型建立完以后,便可以进入/SOLU 处理器, 声明各种负载。但大部分负载的载声明也可在/PREP7 中完成,建义全部 负载在/SOLU 处理中进行声明。 /ANTYPE,Antype,Status 声明分析类型,即欲进行哪种分析,系统默认为静力学分析。 Antype=STATIC or 0 静态分析(系统默认) BUCKLE or 1 屈曲分析 MODAL or 2 振动模态分析 HARMIC or 3 调和外力动和系统 TRANS or 4 瞬时动力系统分析 SFBEAM, ELEM, LKEY, Lab,VALI, VALJ, VAL2I, VAL2J, IOFFST, JOFFST 定义在梁元素上的分布力。 ELEM:元素号码。 LKEY:建立元素后,依节点顺序梁元素有四个面,该参为分力所施加的面 号。 Lab RES(表示分布压力)。 VALI,VALJ:在 I 点及 J 点分布力的值。

前处理 /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;mat 对应前面定义的元素类型号码 Itype;c0 为材料特性类 别的值。 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(该 区域坐标系统与整体坐标系统原点关系)。 csys,kcn !声明坐标系统,系统默认为卡式坐标(csys,0)。 k,npt,x,y,z !定义点。npt 为点的号码;x,y,z 为节点在目前坐标系统 下的坐标位置。 kfill,np1,np2,nfill,nstrt,ninc,space ! 点填充。 np1 和 np2 两点间, nfill 为填充点的个数; nstrt,ninc,space 为分布状态。 kgen,itime,np1,np2,ninc,dx,dy,dz,kinc,noelem,imove !点复制。itime 包含本身所复制的次数;knic 为每次复制时点号码增加 量;np1,np2,ninc 点复制范围;dx,dy,dz 每次复制在现有坐标下几何位 置的改变量。 ksymm,ncomp,np1,np2,ninc, kinc,noelem,imove !复制一组(np1,np2,ninc)点对称于某轴(ncomp);knic 为每次复制时 点号码增加量。 kl,nl1,ratio,nk1 !在已知线(nl1)上建立一个点(nk1),该点的位 置由占全线段比例(radio)而定, 比例为 p1 至 nk1 长度与 p1 至 p2 的长度。 kmodif,npt,x,y,z !修改现有点(npt)到新坐标(x,y,z)位置。 knode,npt,node !定义点(npt)于已知节点(node)上。

kdele,np1,np2,ninc !将一组点删除。 ksel,type,item,comp,vmin,vmax,vinc,kabs !选择有效点,type 为选 择方式。 Wpoffs,xoff,yoff,zoff !将工作平面中心点移到另外一点。 Wprota,thxy,thyz,thzx !将工作平面顺时针旋转一个角度。 l,p1,p2,ndiv,space,xv1,yv1,zv1,xv2,yv2,zv2 !由两点定义线段,此 线段的形状可为直线(斜率)为 0,或为曲线(以线段两端斜率 xv1,yv1,zv1,xv2,yv2,zv2 而定);ndiv 为线段在进行网格化时欲分的元 素数目。 Lstr,p1,p2 !用两个点来定义一条直线。 Lcomb,nl1,nl2,keep !将两条线合并为一条线,keep=0 时原线段删除, keep=1 时保留。 Ldiv,nl1,ratio,pdiv,ndiv,keep !将线分割为数条线,nl1 为线段的号 码; ndiv 为线段欲分的段数 (系统默认为两段) 大于 2 时为均分; , ratio 为两段的比例(等于 2 时才作用);keep=0 时原线段删除,keep=1 时保 留。 Lgen,itime,nl1,nl2,ninc,dx,dy,dz,kinc,noelem,imove !线段复制命 令。itime 包含本身所复制的次数;nl1,nl2,ninc 为现有的坐标系统下复 制到其他位置(dx,dy,dz);kinc 为每次复制时线段号码的增加量。 Lfillt,nl1,nl2,rad,pcent !在两相交的线段 nl1,nl2 间产生一条半径 等于 rad 的圆角曲线,同时自动产生三个点,其中两个点在 nl1,nl2 上, 第三个点是新曲线的圆心定(若 pcent=0,则不产生该点)。 Larc,p1,p2,pc,rad !定义两点间的圆弧线,其半径为 rad,pc 为圆弧曲 率中心部分的任何一点,不一定是圆心坐标。 Circle,pcent,rad,paxis,pzero,arc,nseg !产生圆弧线。Pcent 为圆弧 中心坐标点的号码;paxis 定义圆心轴正方向上任意点的号码;Pzero 定 义圆弧线起点轴上的任意点的号码,此点不一定在圆上;rad 圆的半径; nseg 为圆弧线欲划分的段数,完整为 4。 Lang,nl1,p3,ang,phit,locat !产生一新的线段,此新的线段与已存在 的线段 nl1 的夹角为 ang,phit 为新产生点的号码。 L2ang,nl1,nl2,angl,ang2,phit1,phit2 ! 产生新线段。 此新线段与已存 在的直线 nl1 夹角为 ang1,与直线 nl2 的夹角为 ang2。Phit1,Phit2 为 新产生两点的号码。 Ltan,nl1,P3,xv3,yv3,zv3 !产生三次曲线,该曲线方向为 P2 至 P3, 与已知曲线相切于 P2。Xv3,y,v3,zv3 为新线段在终点 P3 处的斜率。 L2tan,nl1,nl2 ! 建立新线段与已知两条相切的方式产生。 若以负值输入, 则相反。 Bspline,p1,p2,p3,p4,p5,p6,xv1,yv1,zv1,xv6,yv6,zv6 !通过 6 点曲 线,并定义两端点的斜率。 spline,p1,p2,p3,p4,p5,p6,xv1,yv1,zv1,xv6,yv6,zv6 ! 通过 6 点曲线, 每点之间形成一新线段,并可以定义两端点的斜率。 Ldele,nl1,nl2,ninc,kswp ! kswp=0 时只删除掉线段本身, =1 时低单元 点一并删除。 Lsel, type,item,comp,vmin,vmax,vinc,kabs !选择有效线段,type 为

选择方式。 A,p1,p2,p3,p4,p5,p6,p7,p8,p9 !由已知点定义面积 Al,l1,l2,l3,l4,l5,l6,l7,l8,l9,l10 !由已知线段定义面积 Agen, itime,na1,na2,ninc,dx,dy,dz,kinc,noelem,imove !面积复制命 令。itime 包含本身所复制的次数;na1,na2,ninc 为现有的坐标系统下复 制到其他位置(dx,dy,dz);kinc 为每次复制时面积号码的增加量。 Arsym,ncomp,na1,na2,ninc,kinc,noelem,imove !复制一组面积 na1,na2,ninc 对称于轴 ncomp;kinc 为每次复制时面积号码的增加量。 Adrag, nl1,nl2,nl3,nl4,nl5,nl6, nlp1,nlp2,nlp3,nlp4,nlp5,nlp6 ! 面积的建立,沿某组线段路径,拉伸而成。 Arotat, nl1,nl2,nl3,nl4,nl5,nl6,pax1,pax2,arc,nseg !建立一组圆 柱形面积,产生方式为绕着某轴(pax1,pax2,为轴上任意两点,并定义轴 的方向)。Nseg 为整个旋转角度方向中欲分段数目。 Aoffst,narea,dist,kinc !复制一块面积,产生方式为平移(offset) 一块面积,以平面法线方向,平移距离为 dist,kinc 为面积号码增加量。 Afillt,na1,na2,rad !建立圆角面积,在两相交平面间产生曲面,rad 为半径。 Askin,nl1,nl2,nl3,nl4,nl5,nl6 !沿已知线建立一个平滑薄层曲面。 Adele,na1,na2,ninc,kswp !kswp=0 时只删除掉面积本身,=1 时低单元 点一并删除。 Asel, type,item,comp,vmin,vmax,vinc,kabs !选择有效面积,type 为 选择方式。 V,p1,p2,p3,p4,p5,p6,p7,p8 !由点定义体积。 Va,a1,a2,a3,a4,a5,a6,a7,a8,a9,a10 !由已知面定义体积 Vgen,itime,nv1,nv2,ninc,dx,dy,dz,kinc,noelem,imove !体积复制。 Vsymm,ncomp,nv1,nv2,ninc,kinc,noelem,imove!对称于轴(ncomp)复 制一组体积 Vdrag,na1,na2,na3,na4,na5,na6,nlp1,?nlp6 ! 体积建立时将一组已知 面积沿着某组线段路径,拉伸而成。 Vrotat,na1,na2,na3,na4,na5,na6,pax1,pax2,arc,nseg ! 建立一组圆柱 形体积,产生方式为绕着某轴(pax1,pax2,为轴上任意两点,并定义轴的 方向)。Nseg 为整个旋转角度方向中欲分段数目。 Vdele,nv1,nv2,ninc,kswp !kswp=0 时只删除掉体积本身,=1 时低单元 点一并删除。 Vsel, type,item,comp,vmin,vmax,vinc,kabs !选择有效体积,type 为 选择方式。 Rectng,x1,x2,y1,y2 !建立长方形面积 Pcirc,rad1,rad2,theta1,theta2 !建立平面圆面积,rad1,rad2 为圆面 积的内径及外径,theta1,theta2 为圆面积的角度范围。系统默认值为 0 度到 360 度,每 90 度分段。 Rpoly,nsides,lside,majrad,minrad !建立一个以工作面中心点为基准 的正多边形面积。边数为 nsides,大小可由边长 lside,或外接圆半径 majard,或内切圆 minrad。 Block,x1,x2,y1,y2,z1,z2 !建立一个长方体区块。

Blc4,xcorner,ycorner,width,height,depth !建立一个长方体区块。 Blc5,xcenter,ycenter,width,height,depth !建立一个长方体区块。区 块体积中心点的 x、y 坐标。 Cylind,rad1,rad2,z1,z2,theta1,theta2 !建立一个圆柱体积,圆柱的 方向为 z 方向并由 z1,z2 为 z 方向长度的范围;rad1,rad2 为圆柱的内外 半径;theat1,theta2 为圆柱的起始、终结角度。 Cyl4,xcenter,ycenter, rad1, theta1, rad2,theta2,depth !建立一 个圆柱体积。以圆柱体积中心点的 x、y 坐标为基准;rad1,rad2 为圆柱 的内外半径;theat1,theta2 为圆柱的起始、终结角度。 Cyl5,xedge1,yedge1,xedge2,yedge2,depth !建立一个圆柱体积。 xedge1,yedge1,xedge2,yedge2 为圆柱上面或下面任一直径的 x、y 起点 坐标与终点坐标。 Cone,rtop,rbot,z1,z2,theta1,theta2 !建立一个圆锥体积。Rtop,z1 为圆锥上平面的半径与长度、rbot,z2 为圆锥下平面的半径与长度; theat1,theta2 为圆锥的起始、终结角度。 Rprism,z1,z2,nsides,lside,majrad,minrad !建立一个正多边形体 积,z1,z2 为 z 方向长度的范围,边数为 nsides;边长 lside;或外接圆半 径 majard;或内切圆 minrad。 !声明元素大小、形状和网格种类 lesize,nl1,size,angsiz,ndiv,space,kforc,layer1,layer2 ! 定义所选 择线段(nl1,nl1=all 为目前所有的线段)进行元素网格化时元素的大小 (size),元素的大小可用线段的长度(size)或该条线段要分割的元素数目 (ndiv)来确定;space 为间距比(最后一段长与最先一段长的比值,正值 代表以线段方向为基准,负值以中央为基准,系统默认等间距)。 kesize,npt,size,fact1,fact2 !定义通过点(npt,npt=all 为通过目 前所有点的线段)的所有线段进行元素网格划分时元素的大小(size), 不含 lesize 所定义的线段。元素的大小仅能用元素的长度(size)输入。 该命令必须成对使用,因为线段基本上含两点。 esize,size,ndiv ! 定义元素网格化时元素的大小。 该命令以目前所有对 象为基准(不含 lesize,kesize 所定义的线段)。元素的大小可用元素的 边长(size)或线段要分成元素数目(ndiv)来确定。 desize,minl,minh,mxel !系统默认元素大小(不含 lesize,kesize,esize 所定义)。 smrtsize,sizval,fac,expnd,trans,angl,angh,gratio,smhlc,smanc,mx itr,sprx !自由网格时,网格大小的高级控制(不含 lesize,kesize,esize 所定义)。一般由 desize 控制元素大小,desize 及 smrtsize 是相互独立的命令,仅能存在一个,执行 smrtsize 命令后 desize 自动无效。 mshkey,key !key=0 自由网格(系统默认);key=1 对应网格;key=2 对 应自由混合(仅适合 2-D 实体)。 mshape,key,dimension !声明网格化时元素的形状。2-D 实体模型采用 四边形(key=0)或全部为三角形(key=1)。 !进行网格化 xatt,mat,real,type,esys !type 元素的形式号码,real 元素的几何参

数属性编号,mat 元素的材料特性属性编号。Esys 为建立元素时所在坐标 系统号码。系统默认值为第一组及卡式坐标。 Xmesh !x 对象网格化后,元素属性由 xatt 决定。 注意:X 表示 L,A,V 即线,面和体。

月 13 日 谁与我同行? 在 www.okok.org 上发表了这个命令流,但实在是“曲高和寡”。箇中没有太先 进的分析技术,为了建模的简便没有使用接触元,仅使用了初始地应力,也没有 考虑注浆后对岩土的加固作用, 对于实际工况的模拟也有所欠缺——没有考虑盾 构的 300 吨的重量,但通过命令流能够建立这样的分析,我想还是比较罕有的。 本想加上加固作用、应力释放以及盾构重量后,再结合工程实测数据,将其作为 论文发表出去, 然而实在没有时间整理。单位花了八千元买了同济曙光盾构分析 模块, 苦于上海与广州两地岩层的差异, 这一软件于我实在强差人意。 想发脾气, 却猛的看见昨日身为研究生的自己也在编写这样的程序,含了金子,却毛躁的没 有筛清其中的沙石, 硌的人生疼,终究坐在电脑前免费为同济大学作了一次程序 调试人员——免费的。 !=====================================================初始应力分析↓ /INPUT,PARAMETER,APDL *IF,CTRL_READ_INITIALSTRESS,GT,0,THEN FINISH /CLEAR /FILENAME,ISANALYSIS /TITLE,Tunnel digging initial stress analysis /UNITS,SI /PREP7 /INPUT,PARAMETER,APDL /INPUT,MODEL,APDL /SOLU ANTYPE,STATIC !分析类型 NROPT,FULL,,ON AUTOTS,ON NLGEOM,OFF DELTIM,0.1,0.05,0.2 EQSLV,AMG !求解器 CNVTOL,F,,0.02,2,0.5 !力收敛准则

LNSRCH,ON ISWRITE,1 !重力加速度 ACEL,0,0,GRAVITY !约束 NSEL,S,LOC,Y,0 NSEL,A,LOC,Y,S_D D,ALL,UY,0 *IF,CTRL_HALFLIFE,GT,0,THEN NSEL,S,LOC,X,0 D,ALL,ROTX,0 D,ALL,ROTY,0 D,ALL,ROTZ,0 *ELSE NSEL,S,LOC,X,S_B *ENDIF NSEL,A,LOC,X,-S_B D,ALL,UX,0 NSEL,S,LOC,Z,S7_Z2 D,ALL,UZ,0 !杀死壳单元 CMSEL,S,T1_SHELL_ELEM,ELEM CMSEL,A,T2_SHELL_ELEM,ELEM EKILL,ALL !开始进行分析 *IF,CTRL_RUNOW,GT,0,THEN ALLSEL /SOLU SOLVE *ENDIF *ENDIF !=====================================================初始应力分析↑

!=====================================================施工分析↓ FINISH /CLEAR /FILENAME,TDANALYSIS /TITLE,Tunnel digging analysis /UNITS,SI /INPUT,PARAMETER,APDL /PREP7 /INPUT,MODEL,APDL

/SOLU ANTYPE,STATIC !分析类型 NROPT,FULL,,ON AUTOTS,ON NLGEOM,OFF !DELTIM,0.1,0.05,0.2 AUTOTS,ON *IF,CTRL_READ_INITIALSTRESS,LT,1,THEN EQSLV,AMG !求解器 *ENDIF CNVTOL,U,,0.025,1,0.001 !力收敛准则 LNSRCH,ON PRED,ON NCNV,2 NEQIT,50 !重力加速度 ACEL,0,0,GRAVITY !约束 NSEL,S,LOC,Y,0 NSEL,A,LOC,Y,S_D D,ALL,UY,0 *IF,CTRL_HALFLIFE,GT,0,THEN NSEL,S,LOC,X,0 D,ALL,ROTX,0 D,ALL,ROTY,0 D,ALL,ROTZ,0 *ELSE NSEL,S,LOC,X,S_B *ENDIF NSEL,A,LOC,X,-S_B D,ALL,UX,0 NSEL,S,LOC,Z,S7_Z2 D,ALL,UZ,0 !读取初始应力 *IF,CTRL_READ_INITIALSTRESS,GT,0,THEN ALLSEL ISFILE,READ,ISANALYSIS,IST *ENDIF !杀死壳单元 CMSEL,S,T1_SHELL_ELEM,ELEM CMSEL,A,T2_SHELL_ELEM,ELEM EKILL,ALL !开始进行分析(第一天什么都不挖) *IF,CTRL_RUNOW,GT,0,THEN

ALLSEL SOLVE *ENDIF !开挖分析 /INPUT,DiggingAnalysis,APDL *IF,CTRL_RUNOW,GT,0,THEN /POST1 CMSEL,S,T1_SHELL_ELEM,ELEM CMSEL,A,T2_SHELL_ELEM,ELEM ESEL,R,LIVE NSLE,S SET,LAST ETABLE,N11,SMISC,1 ETABLE,N22,SMISC,2 ETABLE,MX,SMISC,4 ETABLE,MY,SMISC,5 ESEL,S,LIVE PLNSOL,U,Z,1,1.0 *ELSE ALLSEL EPLOT *ENDIF !=====================================================施工分析↑

!=====================================================后处理↓ !/POST26 !ALLSEL !NSOL,2,66,U,Z,UZ_1 !NSOL,3,71,U,Z,UZ_2 !NSOL,4,70,U,Z,UZ_3 !NSOL,5,69,U,Z,UZ_4 !NSOL,6,68,U,Z,UZ_5 !NSOL,7,67,U,Z,UZ_6 !NSOL,8,64,U,Z,UZ_7 !STORE,MERGE !XVAR,1 !PLVAR,2,3, 4,5,6,7,8 !=====================================================后处理↑

!=====================================================施加荷载并分析 ↓ *IF,CTRL_TRACK_ELEMENT,GT,0,THEN /SHOW,JPEG *ENDIF !开始进行开挖 *IF,CTRL_RUNOW,GT,0,THEN *DO,DAY,1,CTRL_ELEMDIV_DEEP,1 *IF,CTRL_DIGBOT_FIRSE,GT,0,THEN CMSEL,S,T2_SOLID_ELEM,ELEM *ELSE CMSEL,S,T1_SOLID_ELEM,ELEM *ENDIF NSLE,S NSEL,S,LOC,Y,(DAY-1)*CTRL_ELEMSPAC_DEEP NSEL,A,LOC,Y,DAY*CTRL_ELEMSPAC_DEEP ESLN,R,1 EKILL,ALL *IF,CTRL_DIGBOT_FIRSE,GT,0,THEN CMSEL,S,T2_SOLID_ELEM,ELEM *ELSE CMSEL,S,T1_SOLID_ELEM,ELEM *ENDIF NSEL,S,LOC,Y,(DAY+1)*CTRL_ELEMSPAC_DEEP NSEL,A,LOC,Y,DAY*CTRL_ELEMSPAC_DEEP ESLN,R,1 NSLE,S NSEL,R,LOC,Y,DAY*CTRL_ELEMSPAC_DEEP *IF,CTRL_DIGBOT_FIRSE,GT,0,THEN SF,ALL,PRES,T2_SHELLMACHINE_PRES *ELSE SF,ALL,PRES,T1_SHELLMACHINE_PRES *ENDIF *IF,CTRL_DIGBOT_FIRSE,GT,0,THEN CMSEL,S,T2_SHELL_ELEM,ELEM *ELSE CMSEL,S,T1_SHELL_ELEM,ELEM *ENDIF NSEL,S,LOC,Y,(DAY-1)*CTRL_ELEMSPAC_DEEP NSEL,A,LOC,Y,DAY*CTRL_ELEMSPAC_DEEP ESLN,R,1 EALIVE,ALL NSLE,S !释放围岩压力

*IF,CTRL_DIGBOT_FIRSE,GT,0,THEN RELEASEPRES=T2_RELEASE_PRES *ELSE RELEASEPRES=T1_RELEASE_PRES *ENDIF !水压力 *IF,CTRL_WATERLOAD_APPLIED,GT,0,THEN ENUM=ELNEXT(0) *DOWHILE,ENUM,GT,0 NNUM1=NELEM(ENUM,1) NNUM2=NELEM(ENUM,2) NNUM3=NELEM(ENUM,3) NNUM4=NELEM(ENUM,4) NPRESS1=(WATER_LOAD_UPPOS-NZ(NNUM1))*(WATER_LOAD_DOWN-WATER_LOAD_UP)/ (WATER_LOAD_UPPOS-WATER_LOAD_DOWNPOS) NPRESS2=(WATER_LOAD_UPPOS-NZ(NNUM2))*(WATER_LOAD_DOWN-WATER_LOAD_UP)/ (WATER_LOAD_UPPOS-WATER_LOAD_DOWNPOS) NPRESS3=(WATER_LOAD_UPPOS-NZ(NNUM3))*(WATER_LOAD_DOWN-WATER_LOAD_UP)/ (WATER_LOAD_UPPOS-WATER_LOAD_DOWNPOS) NPRESS4=(WATER_LOAD_UPPOS-NZ(NNUM4))*(WATER_LOAD_DOWN-WATER_LOAD_UP)/ (WATER_LOAD_UPPOS-WATER_LOAD_DOWNPOS) SFE,ENUM,2,PRES,1,NPRESS1+RELEASEPRES,NPRESS2+RELEASEPRES,NPRESS3+REL EASEPRES,NPRESS4+RELEASEPRES ENUM=ELNEXT(ENUM) *ENDDO *ENDIF ALLSEL SOLVE *IF,CTRL_TRACK_ELEMENT,GT,0,THEN ESEL,S,LIVE EPLOT *ENDIF *ENDDO *IF,CTRL_DIG_ALL,GT,0,THEN *DO,DAY,1,CTRL_ELEMDIV_DEEP,1 *IF,CTRL_DIGBOT_FIRSE,GT,0,THEN CMSEL,S,T1_SOLID_ELEM,ELEM *ELSE CMSEL,S,T2_SOLID_ELEM,ELEM *ENDIF NSLE,S NSEL,S,LOC,Y,(DAY-1)*CTRL_ELEMSPAC_DEEP NSEL,A,LOC,Y,DAY*CTRL_ELEMSPAC_DEEP ESLN,R,1

EKILL,ALL *IF,CTRL_RELEASE_PRES,GT,0,THEN ALLSEL SOLVE *ENDIF *IF,CTRL_DIGBOT_FIRSE,GT,0,THEN CMSEL,S,T1_SOLID_ELEM,ELEM *ELSE CMSEL,S,T2_SOLID_ELEM,ELEM *ENDIF NSEL,S,LOC,Y,(DAY+1)*CTRL_ELEMSPAC_DEEP NSEL,A,LOC,Y,DAY*CTRL_ELEMSPAC_DEEP ESLN,R,1 NSLE,S NSEL,R,LOC,Y,DAY*CTRL_ELEMSPAC_DEEP *IF,CTRL_DIGBOT_FIRSE,GT,0,THEN SF,ALL,PRES,T1_SHELLMACHINE_PRES *ELSE SF,ALL,PRES,T2_SHELLMACHINE_PRES *ENDIF *IF,CTRL_DIGBOT_FIRSE,GT,0,THEN CMSEL,S,T1_SHELL_ELEM,ELEM *ELSE CMSEL,S,T2_SHELL_ELEM,ELEM *ENDIF NSEL,S,LOC,Y,(DAY-1)*CTRL_ELEMSPAC_DEEP NSEL,A,LOC,Y,DAY*CTRL_ELEMSPAC_DEEP ESLN,R,1 EALIVE,ALL NSLE,S !释放围岩压力 *IF,CTRL_DIGBOT_FIRSE,GT,0,THEN RELEASEPRES=T1_RELEASE_PRES *ELSE RELEASEPRES=T2_RELEASE_PRES *ENDIF !水压力 *IF,CTRL_WATERLOAD_APPLIED,GT,0,THEN ENUM=ELNEXT(0) *DOWHILE,ENUM,GT,0 NNUM1=NELEM(ENUM,1) NNUM2=NELEM(ENUM,2) NNUM3=NELEM(ENUM,3) NNUM4=NELEM(ENUM,4)

NPRESS1=(WATER_LOAD_UPPOS-NZ(NNUM1))*(WATER_LOAD_DOWN-WATER_LOAD_UP)/ (WATER_LOAD_UPPOS-WATER_LOAD_DOWNPOS) NPRESS2=(WATER_LOAD_UPPOS-NZ(NNUM2))*(WATER_LOAD_DOWN-WATER_LOAD_UP)/ (WATER_LOAD_UPPOS-WATER_LOAD_DOWNPOS) NPRESS3=(WATER_LOAD_UPPOS-NZ(NNUM3))*(WATER_LOAD_DOWN-WATER_LOAD_UP)/ (WATER_LOAD_UPPOS-WATER_LOAD_DOWNPOS) NPRESS4=(WATER_LOAD_UPPOS-NZ(NNUM4))*(WATER_LOAD_DOWN-WATER_LOAD_UP)/ (WATER_LOAD_UPPOS-WATER_LOAD_DOWNPOS) SFE,ENUM,2,PRES,1,NPRESS1+RELEASEPRES,NPRESS2+RELEASEPRES,NPRESS3+REL EASEPRES,NPRESS4+RELEASEPRES ENUM=ELNEXT(ENUM) *ENDDO *ENDIF ALLSEL SOLVE *IF,CTRL_TRACK_ELEMENT,GT,0,THEN ESEL,S,LIVE EPLOT *ENDIF *ENDDO *ENDIF *ENDIF *IF,CTRL_TRACK_ELEMENT,GT,0,THEN /SHOW,CLOSE *ENDIF !=====================================================施加荷载并分析 ↑

!=====================================================定义单元↓ ET,1,SHELL181 KEYOPT,1,3,2 KEYOPT,1,8,2 ET,2,SOLID45 T1_SHELL_TYPENUM=1 T2_SHELL_TYPENUM=1 S1_SOLID_TYPENUM=2 S2_SOLID_TYPENUM=2 S3_SOLID_TYPENUM=2 S4_SOLID_TYPENUM=2 S5_SOLID_TYPENUM=2 S6_SOLID_TYPENUM=2

S7_SOLID_TYPENUM=2 !=====================================================定义单元↑

!=====================================================定义实常数↓ R,1 R,2 R,5 R,6 R,7 R,8 R,9 R,10 R,11 T1_SHELL_REALNUM=1 T2_SHELL_REALNUM=2 S1_SOLID_REALNUM=5 S2_SOLID_REALNUM=6 S3_SOLID_REALNUM=7 S4_SOLID_REALNUM=8 S5_SOLID_REALNUM=9 S6_SOLID_REALNUM=10 S7_SOLID_REALNUM=11 !=====================================================定义实常数↑

!=====================================================定义材料↓ MP,PRXY,1,T1_CON_PRXY MP,EX,1,T1_CON_EX MP,DENS,1,T1_CON_DENS MP,PRXY,2,T2_CON_PRXY MP,EX,2,T2_CON_EX MP,DENS,2,T2_CON_DENS MP,EX,5,S1_EX MP,PRXY,5,S1_PRXY MP,DENS,5,S1_DENS TB,DP,5 TBDATA,1,S1_C,S1_FI,S1_FIF MP,EX,6,S2_EX MP,PRXY,6,S2_PRXY MP,DENS,6,S1_DENS TB,DP,6 TBDATA,1,S2_C,S2_FI,S2_FIF

MP,EX,7,S3_EX MP,PRXY,7,S3_PRXY MP,DENS,7,S1_DENS TB,DP,7 TBDATA,1,S3_C,S3_FI,S3_FIF MP,EX,8,S4_EX MP,PRXY,8,S4_PRXY MP,DENS,8,S1_DENS TB,DP,8 TBDATA,1,S4_C,S4_FI,S4_FIF MP,EX,9,S5_EX MP,PRXY,9,S5_PRXY MP,DENS,9,S1_DENS TB,DP,9 TBDATA,1,S5_C,S5_FI,S5_FIF MP,EX,10,S6_EX MP,PRXY,10,S6_PRXY MP,DENS,10,S1_DENS TB,DP,10 TBDATA,1,S6_C,S6_FI,S6_FIF MP,EX,11,S7_EX MP,PRXY,11,S7_PRXY MP,DENS,11,S1_DENS TB,DP,11 TBDATA,1,S7_C,S7_FI,S7_FIF T1_SHELL_MATNUM=1 T2_SHELL_MATNUM=2 S1_SOLID_MATNUM=5 S2_SOLID_MATNUM=6 S3_SOLID_MATNUM=7 S4_SOLID_MATNUM=8 S5_SOLID_MATNUM=9 S6_SOLID_MATNUM=10 S7_SOLID_MATNUM=11 !=====================================================定义材料↑

!=====================================================定义截面↓ SECTYPE,1,SHELL,,T1_SECTION,3 SECDATA,T1_THICKNESS/3,1 SECDATA,T1_THICKNESS/3,1 SECDATA,T1_THICKNESS/3,1 SECOFFSET,BOT SECTYPE,2,SHELL,,T2_SECTION,3

SECDATA,T2_THICKNESS/3,1 SECDATA,T2_THICKNESS/3,1 SECDATA,T2_THICKNESS/3,1 SECOFFSET,BOT T1_SHELL_SECNUM=1 T2_SHELL_SECNUM=2 !=====================================================定义截面↑

!=====================================================建立几何模型↓ /PREP7 !建立隧道管线上 KSN=1000 LSN=1000 NUMSTR,KP,KSN NUMSTR,LINE,LSN KSEL,NONE LSEL,NONE K,,T1_X,0,T1_Z K,,T1_X,S_D,T1_Z L,1000,1001 !建立隧道管线下 KSN=1100 LSN=1100 NUMSTR,KP,KSN NUMSTR,LINE,LSN KSEL,NONE LSEL,NONE K,,T2_X,0,T2_Z K,,T2_X,S_D,T2_Z L,1100,1101 !建立土层 1 线面 KSN=2100 LSN=2100 ASN=2100 NUMSTR,KP,KSN NUMSTR,LINE,LSN NUMSTR,AREA,ASN KSEL,ALL LSEL,NONE K,2100,-S_B,0,S1_Z1 K,2101,S_B,0,S1_Z1 K,2110,T1_X-T1_R1-T1_SOL_DIST,0,S1_Z1 K,2111,T1_X+T1_R1+T1_SOL_DIST,0,S1_Z1

K,2120,T1_X,0,S1_Z1 K,2200,-S_B,0,S2_Z1 K,2201,S_B,0,S2_Z1 K,2210,T1_X-T1_R1-T1_SOL_DIST,0,S2_Z1 K,2211,T1_X+T1_R1+T1_SOL_DIST,0,S2_Z1 K,2220,T1_X,0,S2_Z1 L,2101,2201 L,2200,2100 L,2110,2210 L,2111,2211 L,2120,2220 LESIZE,ALL,,,CTRL_ELEMDIV_S1 LSEL,NONE L,2110,2120 L,2120,2111 L,2211,2220 L,2220,2210 LESIZE,ALL,,,CTRL_ELEMDIV_GENERAL LSEL,NONE L,2100,2110 L,2111,2101 L,2201,2211 L,2210,2200 LESIZE,ALL,,,CTRL_ELEMDIV_SIDE LSEL,ALL ASEL,NONE A,2100,2110,2210,2200 A,2110,2120,2220,2210 A,2120,2111,2211,2220 A,2111,2101,2201,2211 CM,S1_AREA,AREA !建立土层 2 线面 KSN=2200 LSN=2200 ASN=2200 NUMSTR,KP,KSN NUMSTR,LINE,LSN NUMSTR,AREA,ASN KSEL,ALL LSEL,NONE K,2300,-S_B,0,S3_Z1 K,2301,S_B,0,S3_Z1 K,2310,T1_X-T1_R1-T1_SOL_DIST,0,S3_Z1 K,2311,T1_X+T1_R1+T1_SOL_DIST,0,S3_Z1

CIRCLE,1000,T1_R1,1001 LDIV,ALL,0.5 CM,T1_SOLID_LINE,LINE L,2300,2200 L,2201,2301 L,2311,2211 L,2210,2310 L,2204,1000 L,2208,1000 L,2205,1000 L,2209,1000 L,2202,1000 LESIZE,ALL,,,CTRL_ELEMDIV_GENERAL LSEL,NONE L,2310,2300 L,2301,2311 LESIZE,ALL,,,CTRL_ELEMDIV_SIDE L,2310,2204 L,2210,2208 L,2220,2205 L,2211,2209 L,2311,2202 LESIZE,ALL,,,CTRL_ELEMDIV_45 LSEL,ALL ASEL,NONE A,2200,2210,2310,2300 A,2210,2208,2204,2310 A,2210,2220,2205,2208 A,2220,2211,2209,2205 A,2211,2311,2202,2209 A,2211,2201,2301,2311 CM,S2_AREA,AREA ASEL,NONE A,2204,2208,1000 A,2208,2205,1000 A,2205,2209,1000 A,2209,2202,1000 CM,S2_TUNNEL_AREA,AREA !建立土层 3 线面 KSN=2300 LSN=2300 ASN=2300 NUMSTR,KP,KSN NUMSTR,LINE,LSN

NUMSTR,AREA,ASN KSEL,ALL LSEL,NONE K,2400,-S_B,0,S4_Z1 K,2401,S_B,0,S4_Z1 K,2410,T1_X-T1_R1-T1_SOL_DIST,0,S4_Z1 K,2411,T1_X+T1_R1+T1_SOL_DIST,0,S4_Z1 K,2420,T1_X,0,S4_Z1 L,2300,2400 L,2310,2410 L,2311,2411 L,2301,2401 L,1000,2207 L,1000,2203 L,1000,2206 LESIZE,ALL,,,CTRL_ELEMDIV_GENERAL LSEL,NONE L,2400,2410 L,2411,2401 LESIZE,ALL,,,CTRL_ELEMDIV_SIDE L,2410,2207 L,2420,2203 L,2411,2206 LESIZE,ALL,,,CTRL_ELEMDIV_45 LSEL,ALL ASEL,NONE A,2300,2310,2410,2400 A,2310,2204,2207,2410 A,2207,2203,2420,2410 A,2203,2206,2411,2420 A,2206,2202,2311,2411 A,2202,2311,2411,2206 A,2311,2301,2401,2411 CM,S3_AREA,AREA ASEL,NONE A,2204,1000,2207 A,2207,1000,2203 A,2203,1000,2206 A,2206,1000,2202 CM,S3_TUNNEL_AREA,AREA !建立土层 4 线面 KSN=2400 LSN=2400 ASN=2400

NUMSTR,KP,KSN NUMSTR,LINE,LSN NUMSTR,AREA,ASN KSEL,ALL LSEL,NONE K,2500,-S_B,0,S5_Z1 K,2501,S_B,0,S5_Z1 K,2510,T2_X-T2_R1-T2_SOL_DIST,0,S5_Z1 K,2511,T2_X+T2_R1+T2_SOL_DIST,0,S5_Z1 K,2520,T2_X,0,S5_Z1 L,2510,2520 L,2520,2511 LESIZE,ALL,,,CTRL_ELEMDIV_GENERAL LSEL,NONE L,2500,2510 L,2511,2501 LESIZE,ALL,,,CTRL_ELEMDIV_SIDE L,2400,2500 L,2401,2501 L,2410,2510 L,2420,2520 L,2411,2511 LESIZE,ALL,,,CTRL_ELEMDIV_S4 LSEL,ALL ASEL,NONE A,2400,2410,2510,2500 A,2410,2420,2520,2510 A,2420,2411,2511,2520 A,2411,2401,2501,2511 CM,S4_AREA,AREA !建立土层 5 线面 KSN=2500 LSN=2500 ASN=2500 NUMSTR,KP,KSN NUMSTR,LINE,LSN NUMSTR,AREA,ASN KSEL,ALL LSEL,NONE K,2600,-S_B,0,S6_Z1 K,2601,S_B,0,S6_Z1 K,2610,T2_X-T2_R1-T2_SOL_DIST,0,S6_Z1 K,2611,T2_X+T2_R1+T2_SOL_DIST,0,S6_Z1 CIRCLE,1100,T2_R1,1101

LDIV,ALL,0.5 CM,T2_SOLID_LINE,LINE L,2600,2500 L,2501,2601 L,2611,2511 L,2510,2610 L,2504,1100 L,2508,1100 L,2505,1100 L,2509,1100 L,2502,1100 LESIZE,ALL,,,CTRL_ELEMDIV_GENERAL LSEL,NONE L,2610,2600 L,2601,2611 LESIZE,ALL,,,CTRL_ELEMDIV_SIDE L,2610,2504 L,2510,2508 L,2520,2505 L,2511,2509 L,2611,2502 LESIZE,ALL,,,CTRL_ELEMDIV_45 LSEL,ALL ASEL,NONE A,2500,2510,2610,2600 A,2510,2508,2504,2610 A,2510,2520,2505,2508 A,2520,2511,2509,2505 A,2511,2611,2502,2509 A,2511,2501,2601,2611 CM,S5_AREA,AREA ASEL,NONE A,2504,2508,1100 A,2508,2505,1100 A,2505,2509,1100 A,2509,2502,1100 CM,S5_TUNNEL_AREA,AREA !建立土层 S6 线面 KSN=2600 LSN=2600 ASN=2600 NUMSTR,KP,KSN NUMSTR,LINE,LSN NUMSTR,AREA,ASN

KSEL,ALL LSEL,NONE K,2700,-S_B,0,S7_Z1 K,2701,S_B,0,S7_Z1 K,2710,T2_X-T2_R1-T2_SOL_DIST,0,S7_Z1 K,2711,T2_X+T2_R1+T2_SOL_DIST,0,S7_Z1 K,2720,T2_X,0,S7_Z1 L,2600,2700 L,2610,2710 L,2611,2711 L,2601,2701 L,1100,2507 L,1100,2503 L,1100,2506 LESIZE,ALL,,,CTRL_ELEMDIV_GENERAL LSEL,NONE L,2700,2710 L,2711,2701 LESIZE,ALL,,,CTRL_ELEMDIV_SIDE L,2710,2507 L,2720,2503 L,2711,2506 LESIZE,ALL,,,CTRL_ELEMDIV_45 LSEL,ALL ASEL,NONE A,2600,2610,2710,2700 A,2610,2504,2507,2710 A,2507,2503,2720,2710 A,2503,2506,2711,2720 A,2506,2502,2611,2711 A,2502,2611,2711,2506 A,2611,2601,2701,2711 CM,S6_AREA,AREA ASEL,NONE A,2504,1100,2507 A,2507,1100,2503 A,2503,1100,2506 A,2506,1100,2502 CM,S6_TUNNEL_AREA,AREA !建立土层 7 线面 KSN=2700 LSN=2700 ASN=2700 NUMSTR,KP,KSN

NUMSTR,LINE,LSN NUMSTR,AREA,ASN KSEL,ALL LSEL,NONE K,2800,-S_B,0,S7_Z2 K,2801,S_B,0,S7_Z2 K,2810,T2_X-T2_R1-T2_SOL_DIST,0,S7_Z2 K,2811,T2_X+T2_R1+T2_SOL_DIST,0,S7_Z2 K,2820,T2_X,0,S7_Z2 L,2800,2801 LESIZE,ALL,,,CTRL_ELEMDIV_TNB LSEL,NONE L,2700,2800 L,2701,2801 L,2710,2810 L,2720,2820 L,2711,2811 LESIZE,ALL,,,CTRL_ELEMDIV_S7 LSEL,ALL ASEL,NONE A,2700,2710,2810,2800 A,2710,2720,2820,2810 A,2720,2711,2811,2820 A,2711,2701,2801,2811 CM,S7_AREA,AREA !拉伸围岩 ALLSEL KSN=3000 LSN=3000 ASN=3000 VSN=3000 NUMSTR,KP,KSN NUMSTR,LINE,LSN NUMSTR,AREA,ASN NUMSTR,VOLU,VSN VDRAG,ALL,,,,,,1000 CMSEL,S,T1_SOLID_LINE,LINE ASLL,S LSLA,S LSEL,R,LOC,Y,S_D/2 KSLL,S LSLK,S,1 ASLL,S,1 CM,T1_AREA,AREA

CMSEL,S,T2_SOLID_LINE,LINE ASLL,S LSLA,S LSEL,R,LOC,Y,S_D/2 KSLL,S LSLK,S,1 ASLL,S,1 CM,T2_AREA,AREA !设置纵向直线划分细度 LSEL,S,LOC,Y,S_D/2 LESIZE,ALL,,,CTRL_ELEMDIV_DEEP !=====================================================建立几何模型↑

!=====================================================划分单元↓ ESN=1000 NUMSTR,ELEM,ESN *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,S1_AREA,AREA VSLA,S ESEL,NONE VATT,S1_SOLID_MATNUM,S1_SOLID_REALNUM,S1_SOLID_TYPENUM VMESH,ALL CM,S1_ELEM,ELEM ESN=2000 NUMSTR,ELEM,ESN *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,S2_AREA,AREA VSLA,S ESEL,NONE VATT,S2_SOLID_MATNUM,S2_SOLID_REALNUM,S2_SOLID_TYPENUM VMESH,ALL CM,S2_ELEM,ELEM ESN=3000 NUMSTR,ELEM,ESN *IF,CTRL_HALFLIFE,GT,0,THEN

KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,S3_AREA,AREA VSLA,S ESEL,NONE VATT,S3_SOLID_MATNUM,S3_SOLID_REALNUM,S3_SOLID_TYPENUM VMESH,ALL CM,S3_ELEM,ELEM ESN=4000 NUMSTR,ELEM,ESN *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,S4_AREA,AREA VSLA,S ESEL,NONE VATT,S4_SOLID_MATNUM,S4_SOLID_REALNUM,S4_SOLID_TYPENUM VMESH,ALL CM,S4_ELEM,ELEM ESN=5000 NUMSTR,ELEM,ESN *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,S5_AREA,AREA VSLA,S ESEL,NONE VATT,S5_SOLID_MATNUM,S5_SOLID_REALNUM,S5_SOLID_TYPENUM VMESH,ALL CM,S5_ELEM,ELEM ESN=6000 NUMSTR,ELEM,ESN *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,S6_AREA,AREA

VSLA,S ESEL,NONE VATT,S6_SOLID_MATNUM,S6_SOLID_REALNUM,S6_SOLID_TYPENUM VMESH,ALL CM,S6_ELEM,ELEM ESN=7000 NUMSTR,ELEM,ESN *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,S7_AREA,AREA VSLA,S ESEL,NONE VATT,S7_SOLID_MATNUM,S7_SOLID_REALNUM,S7_SOLID_TYPENUM VMESH,ALL CM,S7_ELEM,ELEM ESN=8000 NUMSTR,ELEM,ESN *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,S2_TUNNEL_AREA,AREA VSLA,S ESEL,NONE VATT,S2_SOLID_MATNUM,S2_SOLID_REALNUM,S2_SOLID_TYPENUM VMESH,ALL CM,S2_TUNNEL_ELEM,ELEM *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,S3_TUNNEL_AREA,AREA VSLA,S ESEL,NONE VATT,S3_SOLID_MATNUM,S3_SOLID_REALNUM,S3_SOLID_TYPENUM VMESH,ALL CM,S3_TUNNEL_ELEM,ELEM CMSEL,A,S2_TUNNEL_ELEM,ELEM CM,T1_SOLID_ELEM,ELEM

ESN=9000 NUMSTR,ELEM,ESN *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,S5_TUNNEL_AREA,AREA VSLA,S ESEL,NONE VATT,S5_SOLID_MATNUM,S5_SOLID_REALNUM,S5_SOLID_TYPENUM VMESH,ALL CM,S5_TUNNEL_ELEM,ELEM *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,S6_TUNNEL_AREA,AREA VSLA,S ESEL,NONE VATT,S6_SOLID_MATNUM,S6_SOLID_REALNUM,S6_SOLID_TYPENUM VMESH,ALL CM,S6_TUNNEL_ELEM,ELEM CMSEL,A,S5_TUNNEL_ELEM,ELEM CM,T2_SOLID_ELEM,ELEM !建立壳单元 ESN=10000 NUMSTR,ELEM,ESN *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,T1_AREA,AREA ESEL,NONE AATT,T1_SHELL_MATNUM,T1_SHELL_REALNUM,T1_SHELL_TYPENUM,,T1_SHELL_SECN UM AMESH,ALL CM,T1_SHELL_ELEM,ELEM ESN=11000 NUMSTR,ELEM,ESN *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B

LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,T2_AREA,AREA ESEL,NONE AATT,T2_SHELL_MATNUM,T2_SHELL_REALNUM,T2_SHELL_TYPENUM,,T2_SHELL_SECN UM AMESH,ALL CM,T2_SHELL_ELEM,ELEM !=====================================================划分单元↑

!=====================================================施加荷载并分析 ↓ !初始应力场分析 /SOLU ANTYPE,STATIC !分析类型 NROPT,FULL,,ON AUTOTS,ON NLGEOM,OFF EQSLV,AMG !求解器 NSUB,2 !子步数 NEQIT,50 !每个子步的平衡迭代次数 CNVTOL,F,5000,0.02,2 !力收敛准则 LNSRCH,ON OUTRES,ALL,ALL PSTRES,ON !布置结构约束以及进行初应力的分析 !重力加速度 ACEL,0,0,GRAVITY NSEL,S,LOC,Y,0 NSEL,A,LOC,Y,S_D D,ALL,UY,0 *IF,CTRL_HALFLIFE,GT,0,THEN NSEL,S,LOC,X,0 D,ALL,ROTX,0 D,ALL,ROTY,0 D,ALL,ROTZ,0 *ELSE NSEL,S,LOC,X,S_B *ENDIF NSEL,A,LOC,X,-S_B D,ALL,UX,0 NSEL,S,LOC,Z,S7_Z2

D,ALL,UZ,0 !杀死壳单元 CMSEL,S,T1_SHELL_ELEM,ELEM CMSEL,A,T2_SHELL_ELEM,ELEM EKILL,ALL !开始进行分析 *IF,CTRL_RUNOW,GT,0,THEN ALLSEL /SOLU SOLVE /PREP7 UPGEOM,-1,,,MERIN,RST /SOLU SOLVE *ENDIF !=====================================================施加荷载并分析 ↑ *IF,CTRL_RUNOW,GT,0,THEN /POST1 CMSEL,S,T1_SHELL_ELEM,ELEM CMSEL,A,T2_SHELL_ELEM,ELEM ESEL,R,LIVE NSLE,S SET,LAST SET,LAST ETABLE,N11,SMISC,1 ETABLE,N22,SMISC,2 ETABLE,MX,SMISC,4 ETABLE,MY,SMISC,5 ESEL,S,LIVE PLNSOL,U,Z,1,1.0 *ELSE ALLSEL EPLOT *ENDIF

!=====================================================改变视角↓ /VIEW,1,0.393639727536,-0.893266354468,0.217078292054 /ANG,1,-63.6619259293 /AUTO,1

/ESHAPE,1 !=====================================================改变视角↑

!=====================================================定义单元↓ ET,1,SHELL181 KEYOPT,1,3,2 KEYOPT,1,8,2 ET,2,SOLID45 T1_SHELL_TYPENUM=1 T2_SHELL_TYPENUM=1 S1_SOLID_TYPENUM=2 S2_SOLID_TYPENUM=2 S3_SOLID_TYPENUM=2 S4_SOLID_TYPENUM=2 S5_SOLID_TYPENUM=2 S6_SOLID_TYPENUM=2 S7_SOLID_TYPENUM=2 !=====================================================定义单元↑

!=====================================================定义实常数↓ R,1 R,2 R,5 R,6 R,7 R,8 R,9 R,10 R,11 T1_SHELL_REALNUM=1 T2_SHELL_REALNUM=2 S1_SOLID_REALNUM=5 S2_SOLID_REALNUM=6 S3_SOLID_REALNUM=7 S4_SOLID_REALNUM=8 S5_SOLID_REALNUM=9 S6_SOLID_REALNUM=10 S7_SOLID_REALNUM=11 !=====================================================定义实常数↑

!=====================================================定义材料↓ MP,PRXY,1,T1_CON_PRXY MP,EX,1,T1_CON_EX MP,DENS,1,T1_CON_DENS MP,PRXY,2,T2_CON_PRXY MP,EX,2,T2_CON_EX MP,DENS,2,T2_CON_DENS MP,EX,5,S1_EX MP,PRXY,5,S1_PRXY MP,DENS,5,S1_DENS TB,DP,5 TBDATA,1,S1_C,S1_FI,S1_FIF MP,EX,6,S2_EX MP,PRXY,6,S2_PRXY MP,DENS,6,S1_DENS TB,DP,6 TBDATA,1,S2_C,S2_FI,S2_FIF MP,EX,7,S3_EX MP,PRXY,7,S3_PRXY MP,DENS,7,S1_DENS TB,DP,7 TBDATA,1,S3_C,S3_FI,S3_FIF MP,EX,8,S4_EX MP,PRXY,8,S4_PRXY MP,DENS,8,S1_DENS TB,DP,8 TBDATA,1,S4_C,S4_FI,S4_FIF MP,EX,9,S5_EX MP,PRXY,9,S5_PRXY MP,DENS,9,S1_DENS TB,DP,9 TBDATA,1,S5_C,S5_FI,S5_FIF MP,EX,10,S6_EX MP,PRXY,10,S6_PRXY MP,DENS,10,S1_DENS TB,DP,10 TBDATA,1,S6_C,S6_FI,S6_FIF MP,EX,11,S7_EX MP,PRXY,11,S7_PRXY MP,DENS,11,S1_DENS TB,DP,11 TBDATA,1,S7_C,S7_FI,S7_FIF T1_SHELL_MATNUM=1 T2_SHELL_MATNUM=2

S1_SOLID_MATNUM=5 S2_SOLID_MATNUM=6 S3_SOLID_MATNUM=7 S4_SOLID_MATNUM=8 S5_SOLID_MATNUM=9 S6_SOLID_MATNUM=10 S7_SOLID_MATNUM=11 !=====================================================定义材料↑

!=====================================================定义截面↓ SECTYPE,1,SHELL,,T1_SECTION,3 SECDATA,T1_THICKNESS/3,1 SECDATA,T1_THICKNESS/3,1 SECDATA,T1_THICKNESS/3,1 SECOFFSET,BOT SECTYPE,2,SHELL,,T2_SECTION,3 SECDATA,T2_THICKNESS/3,1 SECDATA,T2_THICKNESS/3,1 SECDATA,T2_THICKNESS/3,1 SECOFFSET,BOT T1_SHELL_SECNUM=1 T2_SHELL_SECNUM=2 !=====================================================定义截面↑

!=====================================================建立几何模型↓ /PREP7 !建立隧道管线上 KSN=1000 LSN=1000 NUMSTR,KP,KSN NUMSTR,LINE,LSN KSEL,NONE LSEL,NONE K,,T1_X,0,T1_Z K,,T1_X,S_D,T1_Z L,1000,1001 !建立隧道管线下 KSN=1100 LSN=1100 NUMSTR,KP,KSN NUMSTR,LINE,LSN KSEL,NONE

LSEL,NONE K,,T2_X,0,T2_Z K,,T2_X,S_D,T2_Z L,1100,1101 !建立土层 1 线面 KSN=2100 LSN=2100 ASN=2100 NUMSTR,KP,KSN NUMSTR,LINE,LSN NUMSTR,AREA,ASN KSEL,ALL LSEL,NONE K,2100,-S_B,0,S1_Z1 K,2101,S_B,0,S1_Z1 K,2110,T1_X-T1_R1-T1_SOL_DIST,0,S1_Z1 K,2111,T1_X+T1_R1+T1_SOL_DIST,0,S1_Z1 K,2120,T1_X,0,S1_Z1 K,2200,-S_B,0,S2_Z1 K,2201,S_B,0,S2_Z1 K,2210,T1_X-T1_R1-T1_SOL_DIST,0,S2_Z1 K,2211,T1_X+T1_R1+T1_SOL_DIST,0,S2_Z1 K,2220,T1_X,0,S2_Z1 L,2101,2201 L,2200,2100 L,2110,2210 L,2111,2211 L,2120,2220 LESIZE,ALL,,,CTRL_ELEMDIV_S1 LSEL,NONE L,2110,2120 L,2120,2111 L,2211,2220 L,2220,2210 LESIZE,ALL,,,CTRL_ELEMDIV_GENERAL LSEL,NONE L,2100,2110 L,2111,2101 L,2201,2211 L,2210,2200 LESIZE,ALL,,,CTRL_ELEMDIV_SIDE LSEL,ALL ASEL,NONE A,2100,2110,2210,2200

A,2110,2120,2220,2210 A,2120,2111,2211,2220 A,2111,2101,2201,2211 CM,S1_AREA,AREA !建立土层 2 线面 KSN=2200 LSN=2200 ASN=2200 NUMSTR,KP,KSN NUMSTR,LINE,LSN NUMSTR,AREA,ASN KSEL,ALL LSEL,NONE K,2300,-S_B,0,S3_Z1 K,2301,S_B,0,S3_Z1 K,2310,T1_X-T1_R1-T1_SOL_DIST,0,S3_Z1 K,2311,T1_X+T1_R1+T1_SOL_DIST,0,S3_Z1 CIRCLE,1000,T1_R1,1001 LDIV,ALL,0.5 CM,T1_SOLID_LINE,LINE L,2300,2200 L,2201,2301 L,2311,2211 L,2210,2310 L,2204,1000 L,2208,1000 L,2205,1000 L,2209,1000 L,2202,1000 LESIZE,ALL,,,CTRL_ELEMDIV_GENERAL LSEL,NONE L,2310,2300 L,2301,2311 LESIZE,ALL,,,CTRL_ELEMDIV_SIDE L,2310,2204 L,2210,2208 L,2220,2205 L,2211,2209 L,2311,2202 LESIZE,ALL,,,CTRL_ELEMDIV_45 LSEL,ALL ASEL,NONE A,2200,2210,2310,2300 A,2210,2208,2204,2310

A,2210,2220,2205,2208 A,2220,2211,2209,2205 A,2211,2311,2202,2209 A,2211,2201,2301,2311 CM,S2_AREA,AREA ASEL,NONE A,2204,2208,1000 A,2208,2205,1000 A,2205,2209,1000 A,2209,2202,1000 CM,S2_TUNNEL_AREA,AREA !建立土层 3 线面 KSN=2300 LSN=2300 ASN=2300 NUMSTR,KP,KSN NUMSTR,LINE,LSN NUMSTR,AREA,ASN KSEL,ALL LSEL,NONE K,2400,-S_B,0,S4_Z1 K,2401,S_B,0,S4_Z1 K,2410,T1_X-T1_R1-T1_SOL_DIST,0,S4_Z1 K,2411,T1_X+T1_R1+T1_SOL_DIST,0,S4_Z1 K,2420,T1_X,0,S4_Z1 L,2300,2400 L,2310,2410 L,2311,2411 L,2301,2401 L,1000,2207 L,1000,2203 L,1000,2206 LESIZE,ALL,,,CTRL_ELEMDIV_GENERAL LSEL,NONE L,2400,2410 L,2411,2401 LESIZE,ALL,,,CTRL_ELEMDIV_SIDE L,2410,2207 L,2420,2203 L,2411,2206 LESIZE,ALL,,,CTRL_ELEMDIV_45 LSEL,ALL ASEL,NONE A,2300,2310,2410,2400

A,2310,2204,2207,2410 A,2207,2203,2420,2410 A,2203,2206,2411,2420 A,2206,2202,2311,2411 A,2202,2311,2411,2206 A,2311,2301,2401,2411 CM,S3_AREA,AREA ASEL,NONE A,2204,1000,2207 A,2207,1000,2203 A,2203,1000,2206 A,2206,1000,2202 CM,S3_TUNNEL_AREA,AREA !建立土层 4 线面 KSN=2400 LSN=2400 ASN=2400 NUMSTR,KP,KSN NUMSTR,LINE,LSN NUMSTR,AREA,ASN KSEL,ALL LSEL,NONE K,2500,-S_B,0,S5_Z1 K,2501,S_B,0,S5_Z1 K,2510,T2_X-T2_R1-T2_SOL_DIST,0,S5_Z1 K,2511,T2_X+T2_R1+T2_SOL_DIST,0,S5_Z1 K,2520,T2_X,0,S5_Z1 L,2510,2520 L,2520,2511 LESIZE,ALL,,,CTRL_ELEMDIV_GENERAL LSEL,NONE L,2500,2510 L,2511,2501 LESIZE,ALL,,,CTRL_ELEMDIV_SIDE L,2400,2500 L,2401,2501 L,2410,2510 L,2420,2520 L,2411,2511 LESIZE,ALL,,,CTRL_ELEMDIV_S4 LSEL,ALL ASEL,NONE A,2400,2410,2510,2500 A,2410,2420,2520,2510

A,2420,2411,2511,2520 A,2411,2401,2501,2511 CM,S4_AREA,AREA !建立土层 5 线面 KSN=2500 LSN=2500 ASN=2500 NUMSTR,KP,KSN NUMSTR,LINE,LSN NUMSTR,AREA,ASN KSEL,ALL LSEL,NONE K,2600,-S_B,0,S6_Z1 K,2601,S_B,0,S6_Z1 K,2610,T2_X-T2_R1-T2_SOL_DIST,0,S6_Z1 K,2611,T2_X+T2_R1+T2_SOL_DIST,0,S6_Z1 CIRCLE,1100,T2_R1,1101 LDIV,ALL,0.5 CM,T2_SOLID_LINE,LINE L,2600,2500 L,2501,2601 L,2611,2511 L,2510,2610 L,2504,1100 L,2508,1100 L,2505,1100 L,2509,1100 L,2502,1100 LESIZE,ALL,,,CTRL_ELEMDIV_GENERAL LSEL,NONE L,2610,2600 L,2601,2611 LESIZE,ALL,,,CTRL_ELEMDIV_SIDE L,2610,2504 L,2510,2508 L,2520,2505 L,2511,2509 L,2611,2502 LESIZE,ALL,,,CTRL_ELEMDIV_45 LSEL,ALL ASEL,NONE A,2500,2510,2610,2600 A,2510,2508,2504,2610 A,2510,2520,2505,2508

A,2520,2511,2509,2505 A,2511,2611,2502,2509 A,2511,2501,2601,2611 CM,S5_AREA,AREA ASEL,NONE A,2504,2508,1100 A,2508,2505,1100 A,2505,2509,1100 A,2509,2502,1100 CM,S5_TUNNEL_AREA,AREA !建立土层 S6 线面 KSN=2600 LSN=2600 ASN=2600 NUMSTR,KP,KSN NUMSTR,LINE,LSN NUMSTR,AREA,ASN KSEL,ALL LSEL,NONE K,2700,-S_B,0,S7_Z1 K,2701,S_B,0,S7_Z1 K,2710,T2_X-T2_R1-T2_SOL_DIST,0,S7_Z1 K,2711,T2_X+T2_R1+T2_SOL_DIST,0,S7_Z1 K,2720,T2_X,0,S7_Z1 L,2600,2700 L,2610,2710 L,2611,2711 L,2601,2701 L,1100,2507 L,1100,2503 L,1100,2506 LESIZE,ALL,,,CTRL_ELEMDIV_GENERAL LSEL,NONE L,2700,2710 L,2711,2701 LESIZE,ALL,,,CTRL_ELEMDIV_SIDE L,2710,2507 L,2720,2503 L,2711,2506 LESIZE,ALL,,,CTRL_ELEMDIV_45 LSEL,ALL ASEL,NONE A,2600,2610,2710,2700 A,2610,2504,2507,2710

A,2507,2503,2720,2710 A,2503,2506,2711,2720 A,2506,2502,2611,2711 A,2502,2611,2711,2506 A,2611,2601,2701,2711 CM,S6_AREA,AREA ASEL,NONE A,2504,1100,2507 A,2507,1100,2503 A,2503,1100,2506 A,2506,1100,2502 CM,S6_TUNNEL_AREA,AREA !建立土层 7 线面 KSN=2700 LSN=2700 ASN=2700 NUMSTR,KP,KSN NUMSTR,LINE,LSN NUMSTR,AREA,ASN KSEL,ALL LSEL,NONE K,2800,-S_B,0,S7_Z2 K,2801,S_B,0,S7_Z2 K,2810,T2_X-T2_R1-T2_SOL_DIST,0,S7_Z2 K,2811,T2_X+T2_R1+T2_SOL_DIST,0,S7_Z2 K,2820,T2_X,0,S7_Z2 L,2800,2801 LESIZE,ALL,,,CTRL_ELEMDIV_TNB LSEL,NONE L,2700,2800 L,2701,2801 L,2710,2810 L,2720,2820 L,2711,2811 LESIZE,ALL,,,CTRL_ELEMDIV_S7 LSEL,ALL ASEL,NONE A,2700,2710,2810,2800 A,2710,2720,2820,2810 A,2720,2711,2811,2820 A,2711,2701,2801,2811 CM,S7_AREA,AREA !拉伸围岩 ALLSEL

KSN=3000 LSN=3000 ASN=3000 VSN=3000 NUMSTR,KP,KSN NUMSTR,LINE,LSN NUMSTR,AREA,ASN NUMSTR,VOLU,VSN VDRAG,ALL,,,,,,1000 CMSEL,S,T1_SOLID_LINE,LINE ASLL,S LSLA,S LSEL,R,LOC,Y,S_D/2 KSLL,S LSLK,S,1 ASLL,S,1 CM,T1_AREA,AREA CMSEL,S,T2_SOLID_LINE,LINE ASLL,S LSLA,S LSEL,R,LOC,Y,S_D/2 KSLL,S LSLK,S,1 ASLL,S,1 CM,T2_AREA,AREA !设置纵向直线划分细度 LSEL,S,LOC,Y,S_D/2 LESIZE,ALL,,,CTRL_ELEMDIV_DEEP !=====================================================建立几何模型↑

!=====================================================划分单元↓ ESN=1000 NUMSTR,ELEM,ESN ALLSEL *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,S1_AREA,AREA VSLA,S ESEL,NONE VATT,S1_SOLID_MATNUM,S1_SOLID_REALNUM,S1_SOLID_TYPENUM

VMESH,ALL CM,S1_ELEM,ELEM ESN=2000 NUMSTR,ELEM,ESN ALLSEL *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,S2_AREA,AREA VSLA,S ESEL,NONE VATT,S2_SOLID_MATNUM,S2_SOLID_REALNUM,S2_SOLID_TYPENUM VMESH,ALL CM,S2_ELEM,ELEM ESN=3000 NUMSTR,ELEM,ESN ALLSEL *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,S3_AREA,AREA VSLA,S ESEL,NONE VATT,S3_SOLID_MATNUM,S3_SOLID_REALNUM,S3_SOLID_TYPENUM VMESH,ALL CM,S3_ELEM,ELEM ESN=4000 NUMSTR,ELEM,ESN ALLSEL *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,S4_AREA,AREA VSLA,S ESEL,NONE VATT,S4_SOLID_MATNUM,S4_SOLID_REALNUM,S4_SOLID_TYPENUM VMESH,ALL CM,S4_ELEM,ELEM

ESN=5000 NUMSTR,ELEM,ESN ALLSEL *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,S5_AREA,AREA VSLA,S ESEL,NONE VATT,S5_SOLID_MATNUM,S5_SOLID_REALNUM,S5_SOLID_TYPENUM VMESH,ALL CM,S5_ELEM,ELEM ESN=6000 NUMSTR,ELEM,ESN ALLSEL *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,S6_AREA,AREA VSLA,S ESEL,NONE VATT,S6_SOLID_MATNUM,S6_SOLID_REALNUM,S6_SOLID_TYPENUM VMESH,ALL CM,S6_ELEM,ELEM ESN=7000 NUMSTR,ELEM,ESN ALLSEL *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,S7_AREA,AREA VSLA,S ESEL,NONE VATT,S7_SOLID_MATNUM,S7_SOLID_REALNUM,S7_SOLID_TYPENUM VMESH,ALL CM,S7_ELEM,ELEM ESN=8000 NUMSTR,ELEM,ESN

ALLSEL *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,S2_TUNNEL_AREA,AREA VSLA,S ESEL,NONE VATT,S2_SOLID_MATNUM,S2_SOLID_REALNUM,S2_SOLID_TYPENUM VMESH,ALL CM,S2_TUNNEL_ELEM,ELEM ALLSEL *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,S3_TUNNEL_AREA,AREA VSLA,S ESEL,NONE VATT,S3_SOLID_MATNUM,S3_SOLID_REALNUM,S3_SOLID_TYPENUM VMESH,ALL CM,S3_TUNNEL_ELEM,ELEM CMSEL,A,S2_TUNNEL_ELEM,ELEM CM,T1_SOLID_ELEM,ELEM ESN=9000 NUMSTR,ELEM,ESN ALLSEL *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,S5_TUNNEL_AREA,AREA VSLA,S ESEL,NONE VATT,S5_SOLID_MATNUM,S5_SOLID_REALNUM,S5_SOLID_TYPENUM VMESH,ALL CM,S5_TUNNEL_ELEM,ELEM ALLSEL *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1

ASLL,S,1 *ENDIF CMSEL,R,S6_TUNNEL_AREA,AREA VSLA,S ESEL,NONE VATT,S6_SOLID_MATNUM,S6_SOLID_REALNUM,S6_SOLID_TYPENUM VMESH,ALL CM,S6_TUNNEL_ELEM,ELEM CMSEL,A,S5_TUNNEL_ELEM,ELEM CM,T2_SOLID_ELEM,ELEM !建立壳单元 ESN=10000 NUMSTR,ELEM,ESN ALLSEL *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,T1_AREA,AREA ESEL,NONE AATT,T1_SHELL_MATNUM,T1_SHELL_REALNUM,T1_SHELL_TYPENUM,,T1_SHELL_SECN UM AMESH,ALL CM,T1_SHELL_ELEM,ELEM ESN=11000 NUMSTR,ELEM,ESN ALLSEL *IF,CTRL_HALFLIFE,GT,0,THEN KSEL,S,LOC,X,0,-S_B LSLK,S,1 ASLL,S,1 *ENDIF CMSEL,R,T2_AREA,AREA ESEL,NONE AATT,T2_SHELL_MATNUM,T2_SHELL_REALNUM,T2_SHELL_TYPENUM,,T2_SHELL_SECN UM AMESH,ALL CM,T2_SHELL_ELEM,ELEM !=====================================================划分单元↑

!=====================================================主要参数↓ CTRL_RUNOW=1 CTRL_HALFLIFE=1 CTRL_ELEMDIV_GENERAL=2 CTRL_ELEMDIV_SIDE=2 CTRL_ELEMDIV_DEEP=12 CTRL_ELEMDIV_45=2 CTRL_ELEMDIV_S1=2 CTRL_ELEMDIV_S4=1 CTRL_ELEMDIV_S7=2 CTRL_ELEMDIV_TNB=CTRL_ELEMDIV_GENERAL*2+CTRL_ELEMDIV_SIDE*2 CTRL_DIGBOT_FIRSE=1 CTRL_WATERLOAD_APPLIED=1 CTRL_TRACK_ELEMENT=0 CTRL_DIG_ALL=1 CTRL_READ_INITIALSTRESS=1 T1_X=0 T1_Z=-16.301 T2_X=0 T2_Z=-24.701 T1_R1=3 T1_R2=2.7 T2_R1=3 T2_R2=2.7 T1_R=(T1_R1+T1_R2)/2 T2_R=(T2_R1+T2_R2)/2 T1_SOL_DIST=0.6 T2_SOL_DIST=0.6 T1_THICKNESS=T1_R1-T1_R2 T2_THICKNESS=T2_R1-T2_R2 SHELL_DEEP=1.2 S_B=(T1_R1+T2_R)*2 S_D=SHELL_DEEP*CTRL_ELEMDIV_DEEP CTRL_ELEMSPAC_DEEP=S_D/CTRL_ELEMDIV_DEEP S1_DELTA_H=0.00 S2_DELTA_H=0.00 S3_DELTA_H=0.00 S4_DELTA_H=0.00 S5_DELTA_H=0.00 S6_DELTA_H=0.00 S7_DELTA_H=0.00 S1_H=S1_DELTA_H+11 S2_H=S2_DELTA_H+T1_R+T1_SOL_DIST S3_H=S3_DELTA_H+T1_R+T1_SOL_DIST

S5_H=S5_DELTA_H+T2_R+T2_SOL_DIST S6_H=S6_DELTA_H+T2_R+T2_SOL_DIST S7_H=S7_DELTA_H+12.6 S4_H=T1_Z-T2_Z-S3_H-S5_H S1_Z1=T1_Z+S2_H+S1_H S1_Z2=T1_Z+S2_H S2_Z1=T1_Z+S2_H S2_Z2=T1_Z S3_Z1=T1_Z S3_Z2=T1_Z-S3_H S4_Z1=T1_Z-S3_H S4_Z2=T2_Z+S5_H S5_Z1=T2_Z+S5_H S5_Z2=T2_Z S6_Z1=T2_Z S6_Z2=T2_Z-S5_H S7_Z1=T2_Z-S6_H S7_Z2=T2_Z-S6_H-S7_H T1_CON_EX=3.01E10 T2_CON_EX=3.01E10 T1_CON_DENS=2500 T2_CON_DENS=2500 T1_CON_PRXY=0.35 T2_CON_PRXY=0.35 LC1_RATIO=1.35 LC2_RATIO=1.35 LC3_RATIO=1.35 S1_C=8E3 S2_C=25E3 S3_C=25E3 S4_C=29e3 S5_C=50E3 S6_C=50E3 S7_C=300E3 S1_FI=25 S2_FI=28 S3_FI=28 S4_FI=28 S5_FI=28 S6_FI=28 S7_FI=28 S1_FIF=S1_FI S2_FIF=S2_FI S3_FIF=S3_FI

S4_FIF=S4_FI S5_FIF=S5_FI S6_FIF=S6_FI S7_FIF=S7_FI S1_PRXY=0.25 S2_PRXY=0.2 S3_PRXY=0.2 S4_PRXY=0.2 S5_PRXY=0.2 S6_PRXY=0.2 S7_PRXY=0.2 S1_EX=30E6 S2_EX=60E6 S3_EX=60E6 S4_EX=60E6 S5_EX=60E6 S6_EX=60E6 S7_EX=130E6 S1_DENS=1990 S2_DENS=2002 S3_DENS=2002 S4_DENS=2002 S5_DENS=2002 S6_DENS=2002 S7_DENS=2200 GRAVITY=9.818 WATER_DEEPTH=1 WATER_DENS_N=9.800 WATER_LOAD_UPPOS=S1_Z1-WATER_DEEPTH WATER_LOAD_DOWNPOS=S7_Z2 WATER_LOAD_UP=0 WATER_LOAD_DOWN=-(WATER_LOAD_UPPOS-WATER_LOAD_DOWNPOS)*WATER_DENS_N T1_SHELLMACHINE_PRES=18000.0 T2_SHELLMACHINE_PRES=42000.0 T1_RELEASE_PRES=6000.00 T2_RELEASE_PRES=7000.00 SIN45=0.70710678118654752440084436210485 COS45=0.70710678118654752440084436210485 !=====================================================主要参数↑ 12:01 | 添加评论 | 固定链接 | 写入日志 6 月 23 日 我做的地下工程分析

!******************************************************************** ** ! 车陂南盾构井 !******************************************************************** ** FINISH /CLEAR !******************************************************************** ** /PREP7 !******************************************************************** ** !----------------------------------------------------------------------! 全局分析定义 !----------------------------------------------------------------------! 分析类型 ANTYPE,STATIC ! 开关变量 MBAS=1 ! 中板分析开关, 只要大于 0 都会进行分析 HWAS=1 ! 半水分析开关, 只要大于 0 都会进行分析, 否则按照全水进行分析 QZHF=5 ! 强制划分细度 ANOW=1 ! 现在马上分析 ! 荷载定义 WATER_LOAD_UP=21E3 ! 上部水荷载 WATER_LOAD_DOWN=247.21E3 ! 下部水荷载 HALFWATER_LOAD_DOWN=21E3 ! 半水下部荷载 SOIL_LOAD_UP=8.4E3 ! 上部土荷载 SOIL_LOAD_UP2=23.101E3 ! 上部土荷载 2 SOIL_LOAD_DOWN=134.71E3 ! 下部土荷载 SOIL_K_RATIO=5E7 ! 土壤机床系数 SOIL_LOAD_COVER=21E3 ! 覆土荷载 LIVE_LOAD_BIG=33.2E3 ! 活荷载大 LIVE_LOAD_SMA=4E3 ! 活荷载 OVER_LOAD=20E3 ! 超载 SPRING_LINK_AREA=1 ! 用作弹簧的杆件截面大小 SPRING_LINK_LENGTH=1 ! 用作弹簧的杆件长度 SPRING_LINK_AS=9E7 ! 用作弹簧的刚度几近无穷的杆件

! 荷载系数 CASERATIO1=1.35 CASERATIO2=1.35 CASERATIO3=1.35 CASERATIO4=1.4 CASERATIO5=1.4 CASEON1=1 CASEON2=1 CASEON3=1 CASEON4=1 CASEON5=0 ! 荷载与其系数相乘 SOIL_LOAD_COVER=SOIL_LOAD_COVER*CASERATIO1 SOIL_LOAD_UP=SOIL_LOAD_UP*CASERATIO2 SOIL_LOAD_UP2=SOIL_LOAD_UP2*CASERATIO2 SOIL_LOAD_DOWN=SOIL_LOAD_DOWN*CASERATIO2 WATER_LOAD_UP=WATER_LOAD_UP*CASERATIO3 WATER_LOAD_DOWN=WATER_LOAD_DOWN*CASERATIO3 HALFWATER_LOAD_DOWN=HALFWATER_LOAD_DOWN*CASERATIO3 OVER_LOAD=OVER_LOAD*CASERATIO4 LIVE_LOAD_BIG=LIVE_LOAD_BIG*CASERATIO5 LIVE_LOAD_SMA=LIVE_LOAD_SMA*CASERATIO5 !----------------------------------------------------------------------! 单元定义 !----------------------------------------------------------------------ET,1,BEAM189 KEYOPT,1,1,1 ET,2,SHELL181 KEYOPT,2,3,2 KEYOPT,2,8,2 ET,3,COMBIN14 KEYOPT,3,2,3

ET,4,COMBIN39 KEYOPT,4,1,0 KEYOPT,4,2,0 KEYOPT,4,3,3 ET,5,LINK10 KEYOPT,5,3,1 ET,6,LINK10 KEYOPT,6,3,1 ET,7,MPC184 KEYOPT,7,1,1 !----------------------------------------------------------------------! 实常数定义 !----------------------------------------------------------------------R,1, R,2, R,3,3.14E20 R,4,-1,9E9,1,1E-20 R,5,SPRING_LINK_AREA R,6,20 R,7 !----------------------------------------------------------------------! 材料定义 !----------------------------------------------------------------------! C30 混凝土 MP,PRXY,1,0.2 MP,EX,1,3E10 MP,DENS,1,2.7E3 ! C40 混凝土 MP,PRXY,2,0.2 MP,EX,2,3.25E10 MP,DENS,2,2.7E3 ! 用作弹簧的杆件刚度 MP,PRXY,5,0.01 MP,EX,5,3.25E10

! 用作弹簧的杆件刚度 MP,PRXY,6,0.01 MP,EX,6,SPRING_LINK_AS ! 超级无敌刚性杆 MP,PRXY,7,0.01 MP,EX,7,SPRING_LINK_AS !----------------------------------------------------------------------! 截面定义 !----------------------------------------------------------------------! 0.8*0.8 立柱或梁 SECTYPE,1,BEAM,RECT,S1,5 SECDATA,0.8,0.8,4,4 SECOFFSET,CENT ! 0.8*1 中间梁 SECTYPE,2,BEAM,RECT,S2,5 SECDATA,0.8,1,4,4 SECOFFSET,CENT ! 0.8*1.7 顶梁 SECTYPE,3,BEAM,RECT,S3,5 SECDATA,0.8,1.7,4,4 SECOFFSET,USER,0,0.45 ! 1*1.8 底梁 SECTYPE,4,BEAM,RECT,S4,5 SECDATA,1,1.8,4,4 SECOFFSET,USER,0,0.4 ! 0.4 板 SECTYPE,5,SHELL,,B1,5 SECDATA,0.4,1 SECOFFSET,MID ! 0.8 板 SECTYPE,6,SHELL,,B2,5 SECDATA,0.8,1 SECOFFSET,MID ! 1板 SECTYPE,7,SHELL,,B3,5

SECDATA,1,1 SECOFFSET,MID ! 1.6*0.8 底梁 SECTYPE,8,BEAM,RECT,S5,5 SECDATA,1.6,0.8,4,4 SECOFFSET,CENT !----------------------------------------------------------------------! 几何定义 !----------------------------------------------------------------------! 特征点定义 X0=-12.3 X1=-11.5 X2=-4.15 X3=4.15 X4=11.5 X5=12.3 Y0=-7.4 Y1=-6.6 Y2=0 Y3=6.6 Y4=7.4 Z0=-2.5 Z1=0 Z2=8.1 Z3=14.72 Z4=21.37 ! 底层 K,101,X4,Y1,Z1 K,102,X4,Y2,Z1 K,103,X4,Y3,Z1 K,104,X3,Y3,Z1 K,105,X2,Y3,Z1 K,106,X1,Y3,Z1 K,107,X1,Y2,Z1 K,108,X1,Y1,Z1 K,109,X2,Y1,Z1 K,110,X3,Y1,Z1 K,111,X3,Y2,Z1 K,112,X2,Y2,Z1 A,101,102,111,110

A,102,103,104,111 A,111,104,105,112 A,110,111,112,109 A,112,105,106,107 A,109,112,107,108 CM,L1_KP,KP CM,L1_AREA,AREA CM,L1_LINE,LINE KSEL,S,LOC,X,X2 KSEL,A,LOC,X,X3 LSLK,S,1 ASLL,S,1 CM,L1_MIDBOARD_AREA,AREA ASEL,INVE CM,L1_SIDEBOARD_AREA,AREA ! 中下层 KSEL,NONE LSEL,NONE ASEL,NONE CMSEL,S,L1_SIDEBOARD_AREA,AREA AGEN,2,ALL,,,0,0,Z2-Z1,100 CMSEL,U,L1_KP,KP CMSEL,U,L1_LINE,LINE CMSEL,U,L1_AREA,AREA CM,L2_KP,KP CM,L2_LINE,LINE CM,L2_AREA,AREA ! 中上层 KSEL,NONE LSEL,NONE ASEL,NONE CMSEL,S,L1_SIDEBOARD_AREA,AREA AGEN,2,ALL,,,0,0,Z3-Z1,200 CMSEL,U,L1_KP,KP CMSEL,U,L1_LINE,LINE CMSEL,U,L1_AREA,AREA CM,L3_KP,KP CM,L3_LINE,LINE CM,L3_AREA,AREA ! 顶层 KSEL,NONE LSEL,NONE

ASEL,NONE CMSEL,S,L1_AREA,AREA AGEN,2,ALL,,,0,0,Z4-Z1,300 CMSEL,U,L1_KP,KP CMSEL,U,L1_LINE,LINE CMSEL,U,L1_AREA,AREA CM,L4_KP,KP CM,L4_LINE,LINE CM,L4_AREA,AREA ! 侧墙底层 KSEL,ALL LSEL,ALL ASEL,NONE A,101,102,202,201 A,102,103,203,202 A,103,104,204,203 A,104,105,205,204 A,105,106,206,205 A,106,107,207,206 A,107,108,208,207 A,108,109,209,208 A,109,110,210,209 A,110,101,201,210 CM,L1_Q_AREA,AREA ! 侧墙中下层 KSEL,ALL LSEL,ALL ASEL,NONE A,201,202,302,301 A,202,203,303,302 A,203,204,304,303 A,204,205,305,304 A,205,206,306,305 A,206,207,307,306 A,207,208,308,307 A,208,209,309,308 A,209,210,310,309 A,210,201,301,310 CM,L2_Q_AREA,AREA ! 侧墙中上层 KSEL,ALL LSEL,ALL

ASEL,NONE A,301,302,402,401 A,302,303,403,402 A,303,304,404,403 A,304,305,405,404 A,305,306,406,405 A,306,307,407,406 A,307,308,408,407 A,308,309,409,408 A,309,310,410,409 A,310,301,401,410 CM,L3_Q_AREA,AREA ! 中柱 KSEL,ALL LSEL,NONE *DO,I,100,300,100 L,I+11,I+111 L,I+12,I+112 *ENDDO CM,CENT_COL_LINE,LINE ! 底梁 KSEL,S,LOC,X,X2 KSEL,R,LOC,Z,Z1 LSLK,S,1 CM,L1_BEAM1_LINE,LINE KSEL,S,LOC,X,X3 KSEL,R,LOC,Z,Z1 LSLK,S,1 CM,L1_BEAM2_LINE,LINE CMSEL,A,L1_BEAM1_LINE,LINE CM,L1_BEAM_LINE,LINE ! 中下梁 KSEL,S,LOC,X,X2 KSEL,R,LOC,Z,Z2 LSLK,S,1 CM,L2_BEAM1_LINE,LINE KSEL,S,LOC,X,X3 KSEL,R,LOC,Z,Z2 LSLK,S,1 CM,L2_BEAM2_LINE,LINE CMSEL,A,L2_BEAM1_LINE,LINE CM,L2_BEAM_LINE,LINE

! 中上梁 KSEL,S,LOC,X,X2 KSEL,R,LOC,Z,Z3 LSLK,S,1 CM,L3_BEAM1_LINE,LINE KSEL,S,LOC,X,X3 KSEL,R,LOC,Z,Z3 LSLK,S,1 CM,L3_BEAM2_LINE,LINE CMSEL,A,L3_BEAM1_LINE,LINE CM,L3_BEAM_LINE,LINE ! 顶梁 KSEL,S,LOC,X,X2 KSEL,R,LOC,Z,Z4 LSLK,S,1 CM,L4_BEAM1_LINE,LINE KSEL,S,LOC,X,X3 KSEL,R,LOC,Z,Z4 LSLK,S,1 CM,L4_BEAM2_LINE,LINE CMSEL,A,L4_BEAM1_LINE,LINE CM,L4_BEAM_LINE,LINE ! 组东墙 KSEL,S,LOC,X,X4 LSLK,S,1 ASLL,S,1 CM,WALL_EAST_AREA,AREA ! 组西墙 KSEL,S,LOC,X,X1 LSLK,S,1 ASLL,S,1 CM,WALL_WEST_AREA,AREA ! 组南墙 KSEL,S,LOC,Y,Y1 LSLK,S,1 ASLL,S,1 CM,WALL_SOUTH_AREA,AREA ! 组北墙 KSEL,S,LOC,Y,Y3 LSLK,S,1

ASLL,S,1 CM,WALL_NORTH_AREA,AREA ! 选定盾构预留楼板 KSEL,S,LOC,X,X2 KSEL,A,LOC,X,X3 KSEL,U,LOC,Z,Z1 LSLK,S,1 ASLL,S,1 CM,DUNGOU_BROAD_AREA,AREA ! 硬性规定所有直线的划分细度 LSEL,ALL LESIZE,ALL,,,QZHF ! 给定柱属性并划分单元 CMSEL,S,CENT_COL_LINE,LINE LATT,2,1,1,,,,1 LMESH,ALL ! 给定顶梁属性并划分单元 CMSEL,S,L4_BEAM_LINE,LINE LATT,1,1,1,,,,3 LMESH,ALL ! 给定中间梁属性并划分单元 CMSEL,S,L2_BEAM_LINE,LINE CMSEL,A,L3_BEAM_LINE,LINE LATT,1,1,1,,,,2 LMESH,ALL ! 给定底梁属性并划分单元 CMSEL,S,L1_BEAM_LINE,LINE LATT,1,1,1,,,,4 LMESH,ALL ! 给定底板属性并划分单元 CMSEL,S,L1_AREA,AREA AATT,1,2,2,,7 AMESH,ALL ! 给定中间板属性并划分单元 CMSEL,S,L2_AREA,AREA CMSEL,A,L3_AREA,AREA *IF,MBAS,LE,0,THEN

CMSEL,U,DUNGOU_BROAD_AREA,AREA *ENDIF AATT,1,2,2,,5 AMESH,ALL ! 给定顶板属性并划分单元 CMSEL,S,L4_AREA,AREA *IF,MBAS,LE,0,THEN CMSEL,U,DUNGOU_BROAD_AREA,AREA *ENDIF AATT,1,2,2,,6 AMESH,ALL ! 给定属性并划分墙 CMSEL,S,L1_Q_AREA,AREA CMSEL,A,L2_Q_AREA,AREA CMSEL,A,L3_Q_AREA,AREA AATT,1,2,2,,6 AMESH,ALL ! 复制东墙为挡土墙 ALLSEL,ALL NMAX_EAST=0 *GET,NMAX_EAST,NODE,0,NUM,MAXD CMSEL,S,WALL_EAST_AREA,AREA AGEN,2,ALL,,,X5-X4,0,0,1000 CMSEL,U,WALL_EAST_AREA,AREA CM,SRWALL_EAST_AREA,AREA ! 复制西墙为挡土墙 ALLSEL,ALL NMAX_WEST=0 *GET,NMAX_WEST,NODE,0,NUM,MAXD CMSEL,S,WALL_WEST_AREA,AREA AGEN,2,ALL,,,X0-X1,0,0,1000 CMSEL,U,WALL_WEST_AREA,AREA CM,SRWALL_WEST_AREA,AREA ! 复制南墙为挡土墙 ALLSEL,ALL NMAX_SOUTH=0 *GET,NMAX_SOUTH,NODE,0,NUM,MAXD CMSEL,S,WALL_SOUTH_AREA,AREA AGEN,2,ALL,,,0,Y0-Y1,0,2000

CMSEL,U,WALL_SOUTH_AREA,AREA CM,SRWALL_SOUTH_AREA,AREA ! 复制北墙为挡土墙 ALLSEL,ALL NMAX_NORTH=0 *GET,NMAX_NORTH,NODE,0,NUM,MAXD CMSEL,S,WALL_NORTH_AREA,AREA AGEN,2,ALL,,,0,Y4-Y3,0,2000 CMSEL,U,WALL_NORTH_AREA,AREA CM,SRWALL_NORTH_AREA,AREA ! 抗拔桩 NSEL,S,LOC,X,X2 NSEL,A,LOC,X,X3 NSEL,R,LOC,Y,Y2 NSEL,R,LOC,Z,Z1 D,ALL,UZ,0 ! 挡土墙 CMSEL,S,SRWALL_EAST_AREA,AREA CMSEL,A,SRWALL_WEST_AREA,AREA CMSEL,A,SRWALL_SOUTH_AREA,AREA CMSEL,A,SRWALL_NORTH_AREA,AREA NSLA,S,1 NSEL,R,LOC,Z,Z1 D,ALL,ALL,0 ! 制造基础节点 NSEL,ALL NNUMMAX=5000 *GET,NNUMMAX,NODE,0,NUM,MAX CMSEL,S,L1_AREA,AREA NSLA,S,1 INTERDIST=(Y2-Y1)/(QZHF*2) REACT=INTERDIST/100 *DO,YPOS,Y1+INTERDIST,Y2-INTERDIST,INTERDIST*2 NSEL,U,LOC,Y,YPOS-REACT,YPOS+REACT *ENDDO INTERDIST=(Y3-Y2)/(QZHF*2) *DO,YPOS,Y2+INTERDIST,Y3-INTERDIST,INTERDIST*2 NSEL,U,LOC,Y,YPOS-REACT,YPOS+REACT *ENDDO CM,L1_ROOT_NODE,NODE NGEN,2,NNUMMAX,ALL,,,0,0,-SPRING_LINK_LENGTH

CMSEL,U,L1_ROOT_NODE,NODE CM,L1_BASE_NODE,NODE ! 制造底座弹簧 CMSEL,S,L1_ROOT_NODE,NODE ESEL,NONE NNUM=NDNEXT(0) NCOUNT=10 *GET,NCOUNT,NODE,0,COUNT SINGLESOIL_K_RATIO=((X5-X0)*(Y4-Y0)*SOIL_K_RATIO)/NCOUNT R,4,0,0,1,SINGLESOIL_K_RATIO MP,EX,5,SINGLESOIL_K_RATIO TYPE,5 REAL,5 MAT,5 *DOWHILE,NNUM,GT,0 E,NNUM,NNUM+NNUMMAX NNUM=NDNEXT(NNUM) *ENDDO CM,L1_SPRING_ELEM,ELEM ! 组织结构主体单元 CMSEL,S,SRWALL_EAST_AREA CMSEL,A,SRWALL_WEST_AREA CMSEL,A,SRWALL_SOUTH_AREA CMSEL,A,SRWALL_NORTH_AREA ESLA,S ESEL,INVE CM,MAINSTRUCTURE_ELEM,ELEM ! 建立挡土墙其余部分 K,3101,X5,Y0,Z1 K,3103,X5,Y4,Z1 K,3106,X0,Y4,Z1 K,3108,X0,Y0,Z1 K,3201,X5,Y0,Z2 K,3203,X5,Y4,Z2 K,3206,X0,Y4,Z2 K,3208,X0,Y0,Z2 K,3301,X5,Y0,Z3 K,3303,X5,Y4,Z3 K,3306,X0,Y4,Z3 K,3308,X0,Y0,Z3 K,3401,X5,Y0,Z4 K,3403,X5,Y4,Z4

K,3406,X0,Y4,Z4 K,3408,X0,Y0,Z4 A,3101,1101,1201,3201 A,3201,1201,1301,3301 A,3301,1301,1401,3401 A,1103,3103,3203,1203 A,1203,3203,3303,1303 A,1303,3303,3403,1403 A,3103,2103,2203,3203 A,3203,2203,2303,3303 A,3303,2303,2403,3403 A,2106,3106,3206,2206 A,2206,3206,3306,2306 A,2306,3306,3406,2406 A,3106,1106,1206,3206 A,3206,1206,1306,3306 A,3306,1306,1406,3406 A,1108,3108,3208,1208 A,1208,3208,3308,1308 A,1308,3308,3408,1408 A,3108,2108,2208,3208 A,3208,2208,2308,3308 A,3308,2308,2408,3408 A,2101,3101,3201,2201 A,2201,3201,3301,2301 A,2301,3301,3401,2401 LESIZE,ALL,,,QZHF AATT,1,2,2,,6 AMESH,ALL ! 连接挡土墙与侧墙 HNNUM=0 CMSEL,S,WALL_EAST_AREA,AREA NSLA,S,1 NNUM=NDNEXT(0) *DOWHILE,NNUM,GT,0 NSEL,ALL TNX=NX(NNUM) TNY=NY(NNUM) TNZ=NZ(NNUM) HNNUM=NODE(X5,TNY,TNZ) *IF,TNZ,EQ,Z4,THEN TYPE,7 REAL,7

MAT,7 *ELSE TYPE,6 REAL,6 MAT,6 *ENDIF *IF,HNNUM,GT,0,THEN E,NNUM,HNNUM *ENDIF NSLA,S,1 NNUM=NDNEXT(NNUM) *ENDDO CMSEL,S,WALL_WEST_AREA,AREA NSLA,S,1 NNUM=NDNEXT(0) *DOWHILE,NNUM,GT,0 NSEL,ALL TNX=NX(NNUM) TNY=NY(NNUM) TNZ=NZ(NNUM) HNNUM=NODE(X0,TNY,TNZ) *IF,TNZ,EQ,Z4,THEN TYPE,7 REAL,7 MAT,7 *ELSE TYPE,6 REAL,6 MAT,6 *ENDIF *IF,HNNUM,GT,0,THEN E,NNUM,HNNUM *ENDIF NSLA,S,1 NNUM=NDNEXT(NNUM) *ENDDO CMSEL,S,WALL_SOUTH_AREA,AREA NSLA,S,1 NNUM=NDNEXT(0) *DOWHILE,NNUM,GT,0 NSEL,ALL TNX=NX(NNUM)

TNY=NY(NNUM) TNZ=NZ(NNUM) HNNUM=NODE(TNX,Y0,TNZ) *IF,TNZ,EQ,Z4,THEN TYPE,7 REAL,7 MAT,7 *ELSE TYPE,6 REAL,6 MAT,6 *ENDIF *IF,HNNUM,GT,0,THEN E,NNUM,HNNUM *ENDIF NSLA,S,1 NNUM=NDNEXT(NNUM) *ENDDO CMSEL,S,WALL_NORTH_AREA,AREA NSLA,S,1 NNUM=NDNEXT(0) *DOWHILE,NNUM,GT,0 NSEL,ALL TNX=NX(NNUM) TNY=NY(NNUM) TNZ=NZ(NNUM) HNNUM=NODE(TNX,Y4,TNZ) *IF,TNZ,EQ,Z4,THEN TYPE,7 REAL,7 MAT,7 *ELSE TYPE,6 REAL,6 MAT,6 *ENDIF *IF,HNNUM,GT,0,THEN E,NNUM,HNNUM *ENDIF NSLA,S,1 NNUM=NDNEXT(NNUM) *ENDDO CM,SRWNW_SPRING_ELEM,ELEM

!******************************************************************** ** /SOLU !******************************************************************** ** ! 设置荷载叠加 !SFCUM,PRES,ADD ! 给定重力加速度 ACEL,0,0,9.8 ! 顶板荷载 UPBOARDLOAD=0 *IF,CASEON1,GT,0,THEN UPBOARDLOAD=UPBOARDLOAD+SOIL_LOAD_UP2 *ENDIF *IF,CASEON3,GT,0,THEN UPBOARDLOAD=UPBOARDLOAD+WATER_LOAD_UP *ENDIF *IF,CASEON4,GT,0,THEN UPBOARDLOAD=UPBOARDLOAD+OVER_LOAD *ENDIF *IF,UPBOARDLOAD,NE,0,THEN CMSEL,S,L4_AREA,AREA *IF,MBAS,LE,0,THEN CMSEL,U,DUNGOU_BROAD_AREA,AREA *ENDIF ESLA,S SFE,ALL,2,PRES,1,UPBOARDLOAD,UPBOARDLOAD,UPBOARDLOAD,UPBOARDLOAD *ENDIF ! 侧墙荷载 *IF,CASEON3,GT,0,THEN CMSEL,S,L1_Q_AREA,AREA CMSEL,A,L2_Q_AREA,AREA CMSEL,A,L3_Q_AREA,AREA ESLA,S ENUM=ELNEXT(0) *DOWHILE,ENUM,GT,0 NNUM1=NELEM(ENUM,1) NNUM2=NELEM(ENUM,2) NNUM3=NELEM(ENUM,3) NNUM4=NELEM(ENUM,4) NLOC1=NZ(NNUM1) NLOC2=NZ(NNUM2)

NLOC3=NZ(NNUM3) NLOC4=NZ(NNUM4) NPRES1=(Z4-NLOC1)*(WATER_LOAD_DOWN-WATER_LOAD_UP)/(Z4-Z1)+WATER_LOAD_ UP NPRES2=(Z4-NLOC2)*(WATER_LOAD_DOWN-WATER_LOAD_UP)/(Z4-Z1)+WATER_LOAD_ UP NPRES3=(Z4-NLOC3)*(WATER_LOAD_DOWN-WATER_LOAD_UP)/(Z4-Z1)+WATER_LOAD_ UP NPRES4=(Z4-NLOC4)*(WATER_LOAD_DOWN-WATER_LOAD_UP)/(Z4-Z1)+WATER_LOAD_ UP ESEL,S,,,ENUM SFE,ALL,2,PRES,1,NPRES1,NPRES2,NPRES3,NPRES4 ESLA,S ENUM=ELNEXT(ENUM) *ENDDO *ENDIF ! 挡土墙荷载 *IF,CASEON2,GT,0,THEN CMSEL,S,SRWALL_EAST_AREA,AREA CMSEL,A,SRWALL_WEST_AREA,AREA CMSEL,A,SRWALL_NORTH_AREA,AREA CMSEL,A,SRWALL_SOUTH_AREA,AREA ESLA,S ENUM=ELNEXT(0) *DOWHILE,ENUM,GT,0 NNUM1=NELEM(ENUM,1) NNUM2=NELEM(ENUM,2) NNUM3=NELEM(ENUM,3) NNUM4=NELEM(ENUM,4) NLOC1=NZ(NNUM1) NLOC2=NZ(NNUM2) NLOC3=NZ(NNUM3) NLOC4=NZ(NNUM4) NPRES1=(Z4-NLOC1)*(SOIL_LOAD_DOWN-SOIL_LOAD_UP)/(Z4-Z1)+SOIL_LOAD_UP NPRES2=(Z4-NLOC2)*(SOIL_LOAD_DOWN-SOIL_LOAD_UP)/(Z4-Z1)+SOIL_LOAD_UP NPRES3=(Z4-NLOC3)*(SOIL_LOAD_DOWN-SOIL_LOAD_UP)/(Z4-Z1)+SOIL_LOAD_UP NPRES4=(Z4-NLOC4)*(SOIL_LOAD_DOWN-SOIL_LOAD_UP)/(Z4-Z1)+SOIL_LOAD_UP ESEL,S,,,ENUM SFE,ALL,2,PRES,1,NPRES1,NPRES2,NPRES3,NPRES4 ESLA,S ENUM=ELNEXT(ENUM) *ENDDO *ENDIF

! 底板荷载 LOWBOARDLOAD1=0 LOWBOARDLOAD2=0 *IF,CASEON3,GT,0,THEN LOWBOARDLOAD1=LOWBOARDLOAD1+WATER_LOAD_DOWN LOWBOARDLOAD2=LOWBOARDLOAD2+WATER_LOAD_DOWN *ENDIF *IF,CASEON5,GT,0,THEN LOWBOARDLOAD1=LOWBOARDLOAD1-LIVE_LOAD_BIG LOWBOARDLOAD2=LOWBOARDLOAD2-LIVE_LOAD_SMA *ENDIF *IF,LOWBOARDLOAD1,NE,0,THEN CMSEL,S,L1_SIDEBOARD_AREA,AREA ESLA,S SFE,ALL,1,PRES,1,LOWBOARDLOAD1,LOWBOARDLOAD1,LOWBOARDLOAD1,LOWBOARDLO AD1 *ENDIF *IF,LOWBOARDLOAD2,NE,0,THEN CMSEL,S,L1_MIDBOARD_AREA,AREA ESLA,S SFE,ALL,1,PRES,1,LOWBOARDLOAD2,LOWBOARDLOAD2,LOWBOARDLOAD2,LOWBOARDLO AD2 *ENDIF !----------------------------------------------------------------------! 求解 !----------------------------------------------------------------------*IF,ANOW,GT,0,THEN CMSEL,S,L1_BASE_NODE,NODE D,ALL,ALL,0 ALLSEL,ALL ! 分析选项 NSUBST,10 TIME,1 AUTO,ON NLGEOM,ON LNSRCH,ON SOLVE *ENDIF

!******************************************************************** ** /POST1 !******************************************************************** ** *IF,ANOW,GT,0,THEN SET,FIRST LSEL,ALL ESLL,S ETABLE,A,SMISC,2 ETABLE,B,SMISC,15 PLLS,A,B ETABLE,C,SMISC,3 ETABLE,D,SMISC,16 PLLS,C,D ASEL,ALL ESLA,S ETABLE,MX,SMISC,4 ETABLE,MY,SMISC,5 ETABLE,MZ,SMISC,6 PRETAB,MX,MY,MZ CMSEL,S,MAINSTRUCTURE_ELEM,ELEM !/PSF,PRES,NORM,1,0,ON *ENDIF

ANSYS 移动热源 apdl 程序

2009-8-26 16:39

finish /clear /prep7 /nopr et,1,solid70 mp,dens,1,5250 mp,EX,1,2e11 mp,kxx,1,30 mp,kyy,1,30 mp,kzz,1,15 mp,PRXY,1,0.3!no printf as devo mp,pryz,1,0.3

mp,prxz,1,0.3 mp,c,1,450 mp,prxy,1,0.3 save csys,1 cylind,,50,,2000 vsweep,all /solu !*afun,deg *dim,flux2,table,24,16,16,y,z,time *do,k,1,16,1 *do,i,1,24,1 flux2(i,0,k)=(i-1)*15 *enddo *do,j,1,16,1 flux2(0,j,k)=(j-1)*2000/15 *enddo *do,i,1,24,1 *do,j,1,16,1 dis=sqrt(((j-1)*2000/15-(k-1)*2000/15)**2+(((i-1)*15-45)/180)**2) *if,dis,eq,0,then flux2(i,j,k)=8000000000 *else flux2(i,j,k)=0 *endif *enddo *enddo flux2(0,0,k)=(k-1)*0.1 *enddo antype,trans tunif,25 kbc,1 *set,delt,1 timint,on autots,off tt=1.6

*do,tm,1e-6,tt,0.1 /view,1,1,1,1 /ang,1 /rep,fast /go time,tm

autots,0 deltim,,,,1 kbc,0 tsres,erase nsel,s,loc,x,50 !nsel,r,loc,y,0,245 allsel,all SF,all,HFLUX,%flux2% /status,solu solve outres,all,all /psf,hflux,,2 /replot *enddo outres,all,all finish

/SOLU LSCLEAR,SOLID /PREP7 ETCHG,TTS FINISH /SOLU asel,s,loc,z,0 D,all,uz,0 *DO,TM,1E-6,1.6,0.1 TIME,TM LDREAD,TEMP,,,TIME,0,'file','RTH' DELTIM,,,,1 SOLVE *ENDDO FINISH


相关文章:
ANSYS+APDL+命令集合
ANSYS+APDL+命令集合_工学_高等教育_教育专区。ANSYS+APDL+命令集合ANSYS APDL 命令集合(2) 2007-11-23 16:50:04 /post1 中的几个命令: set, lstep, sb...
ANSYS+APDL+命令集合
ANSYS+APDL+命令集合_机械/仪表_工程科技_专业资料。APDL命令ANSYS APDL 命令集合(2) 2007-11-23 16:50:04 /post1 中的几个命令: set, lstep, sbstep, ...
ANSYSAPDL命令汇总
ANSYSAPDL命令汇总_计算机软件及应用_IT/计算机_专业资料。ANSYS 命令大全 A A,...Lab6 该命令与 D 命令相对应, 定义约束, KPOI 为受限点的号码, VALUE 为受...
ANSYS APDL命令汇总
ANSYS APDL 命令汇总 A A,P1,P2,P3,P4,P5,P6,P7,P8,P9 此命令用已知的...,cname8 将组元分组形成组元集合 aname: 组元集名称 cname1……cname8: 已...
经典APDL-ANSYS
经典APDL-ANSYS_计算机软件及应用_IT/计算机_专业资料。今日推荐 180...八天学会ansys apdl 10页 免费 ANSYS+APDL+命令集合 78页 2下载券 基于ANSYS...
Ansys 常用命令全集
Assembly 组元/ 组元/集合 Plot/Replot 画图/ 画图/重新画图 Pan,Zoom,Rotate...APDL ANSYS 参数化设计 语言 Line Fillet 在两条线的过渡生 成线 Arbitrary ...
ANSYS命令汇总
ansys中典型案例命令流汇... 3页 1下载券 ANSYS APDL命令汇总 44页 1下载券...NL1 也可以是集合名。 PAX1, PAX2 定义旋转轴的关键点。 ARC 弧长(度)。 ...
ANSYS命令大全
ANSYS命令大全_计算机软件及应用_IT/计算机_专业资料 暂无评价|0人阅读|0次下载|举报文档ANSYS命令大全_计算机软件及应用_IT/计算机_专业资料。APDL 换行与续行APDL...
ANSYS中的APDL命令
ANSYS中的APDL命令_工学_高等教育_教育专区。结合自身经验, 命令(二 结合自身...,cname8 将组元分组形成组元集合 aname: 组元集名称 cname1……cname8: 已...
ANSYS中的APDL命令
荷载步文件范围 ANSYS 中的 APDL 命令 二) 命令(二 1.1 /prep7(进入前...,cname8 将组元分组形成组元集合 aname: 组元集名称 cname1……cname8: 已...
更多相关标签:
ansys apdl命令流 | ansys apdl 命令 | ansys apdl教程 | ansys apdl实例 | ansys apdl | ansys apdl打不开 | ansys apdl视频教程 | ansys17.0 apdl打不开 |