隧道强度检验ANSYS命令流

合集下载

有限元分析软件ANSYS命令流中文说明

有限元分析软件ANSYS命令流中文说明有限元分析软件ANSYS命令流中文说明（1）CommandVSBV, NV1, NV2, SEPO, KEEP1, KEEP2 —Subtracts volumes from volumes，用于2个solid相减操作，最终目的是要nv1-nv2=?通过后面的参数设置，可以得到很多种情况：sepo项是2个体的边界情况，当缺省的时候，是表示2个体相减后，其边界是公用的，当为sepo的时候，表示相减后，2个体有各自的独立边界。

keep1与keep2是询问相减后，保留哪个体？当第一个为keep时，保留nv1,都缺省的时候，操作结果最终只有一个体，比如：vsbv,1,2,sepo,,keep,表示执行1-2的操作，结果是保留体2，体1被删除，还有一个1-2的结果体，现在一共是2个体（即1-2与2），且都各自有自己的边界。

如vsbv,1,2,,ke ep,,则为1-2后，剩下体1和体1-2，且2个体在边界处公用。

同理，将v换成a及l是对面和线进行减操作！mp,lab, mat, co, c1,…….c4 定义材料号及特性lab: 待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens）ex: 弹性模量nuxy: 小泊松比alpx: 热膨胀系数reft: 参考温度reft: 参考温度prxy: 主泊松比gxy: 剪切模量mu: 摩擦系数dens: 质量密度mat: 材料编号（缺省为当前材料号）co: 材料特性值，或材料之特性，温度曲线中的常数项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,degVSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWPType，是选择的方式，有选择（s）,补选（a），不选(u)，全选(all)、反选（inv）等,其余方式不常用Item, Comp 是选取的原则以及下面的子项如volu 就是根据实体编号选择，loc 就是根据坐标选取，它的comp就可以是实体的某方向坐标！其余还有材料类型、实常数等MIN, VMAX, VINC，这个就不必说了吧！,例：vsel,s,volu,,14vsel,a,volu,,17,23,2上面的命令选中了实体编号为14，17，19，21，23的五个实体VDELE, NV1, NV2, NINC, KSWP: 删除未分网格的体nv1:初始体号nv2:最终的体号ninc:体号之间的间隔kswp=0:只删除体kswp=1:删除体及组成关键点,线面如果nv1=all,则nv2,ninc不起作用其后面常常跟着一条显示命令VPLO,或aplo,nplo,这个湿没有参数的命令，输入后直接回车，就可以显示刚刚选择了的体、面或节点，很实用的哦！Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点为下一步做准备Type: S: 选择一组新节点（缺省）R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0” 使用正负号“1”仅用绝对值下面是单元生死第一个载荷步中命令输入示例：!第一个载荷步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命令流.doc

八天学会 Ansys 命令流为方便大家的交流和学习, 特推出 " 跟我学命令流 " 课程本课程分为三部分: 前处理 , 加载求解 , 后处理每部分的学习时间:10 天 , 共计 30 天每天学习大约10 个命令希望本课程对大家能有所帮助第一天目标 : 熟悉 ANSYS基本关键字的含义k --> Keypoints 关键点l --> Lines 线a --> Area 面v --> Volumes 体e --> Elements 单元n --> Nodes 节点cm --> component 组元et --> element type 单元类型mp --> material property 材料属性r --> real constant 实常数d --> DOF constraint约束f --> Force Load 集中力sf --> Surface load on nodes 表面载荷bf --> Body Force on Nodes 体载荷ic --> Initial Conditions 初始条件第二天目标 : 了解命令流的整体结构, 掌握每个模块的标识!文件说明段/BATCH/TITILE,test analysis !定义工作标题/FILENAME,test !定义工作文件名/PREP7 ! 进入前处理模块标识!定义单元 , 材料属性 , 实常数段ET,1,SHELL63 ! 指定单元类型ET,2,SOLID45 !指定体单元MP,EX,1,2E8 !指定弹性模量MP,PRXY,1,0.3 !入泊松比MP,DENS,1,7.8E3 ! 入材料密度R,1,0.001 !指定壳元常数- 厚度......!建立模型K,1,0,0,, !定关点K,2,50,0,,K,3,50,10,,K,4,10,10,,K,5,10,50,,K,6,0,50,,A,1,2,3,4,5,6, !由关点生成面......!划分网格ESIZE,1,0,AMESH,1......FINISH ! 前理束/SOLU ! 入求解模!施加束和荷DL,5,,ALLSFL,3,PRES,1000SFL,2,PRES,1000......SOLVE !求解FINISH ! 求解模束/POST1 ! 入通用后理器....../POST26 ! 入程后理器⋯⋯/EXIT,SAVE ! 退出并存以下是日志文件中常出的一些命令的明, 希望能大家在整理LOG文件有所帮助/ANGLE ! 指定绕轴旋转视图/DIST !说明对视图进行缩放/DEVICE ! 设置图例的显示, 如: 风格 , 字体等/REPLOT ! 重新显示当前图例/RESET ! 恢复缺省的图形设置/VIEW ! 设置观察方向/ZOOM !对图形显示窗口的某一区域进行缩放第三天生成关键点和线部分1.生成关键点K, 关键点编号 ,X 坐标 ,Y 坐标 ,Z 坐标例:K,1,0,0,02.在激活坐标系生成直线LSTR,关键点 P1, 关键点 P2例LSTR,1,23.在两个关键点之间连线L, 关键点 P1, 关键点 P2例L,1,2注: 此命令会随当前的激活坐标系不同而生成直线或弧线4.由三个关键点生成弧线LARC,关键点 P1, 关键点 P2, 关键点 PC,半径 RAD例LARC,1,3,2,0.05注: 关键点 PC是用来控制弧线的凹向5.通过圆心半径生成圆弧CIRCLE,关键点圆心 , 半径 RAD,,,, 圆弧段数NSEG例:CIRCLE,1,0.05,,,,46.通过关键点生成样条线BSPLIN,关键点 P1, 关键点 P2,关键点 P3, 关键点 P4, 关键点 P5, 关键点 P6 例:BSPLIN,1,2,3,4,5,67.生成倒角线LFILLT, 线 NL1, 线 NL2,倒角半径RAD例LFILLT,1,2,0.0058.通过关键点生成面A, 关键点 P1, 关键点 P2, 关键点 P3, 关键点 P4, 关键点 P5, 关键点 P6,P7,P8...例:A,1,2,3,49.通过线生成面AL, 线 L1, 线 L2, 线 L3, 线 L4, 线 L5, 线 L6, 线 L7, 线 L8, 线 L9, 线 L10例:AL,5,6,7,810.通过线的滑移生成面ASKIN, 线 NL1,线 NL2,线 NL3,线 NL4, 线 NL5, 线 NL6, 线 NL7,线 NL8, 线 NL9例:ASKIN,1,4,5,6,7,8注:线 1 为滑移的导向线第四天目标 : 掌握常用的实体 - 面的生成生成矩形面1.通过矩形角上定位点生成面BLC4,定位点 X 方向坐标 XCORNER,定位点 Y 方向坐标 YCORNER,矩形宽度 WIDTH,矩形高度HEIGHT,矩形深度 DEPTH例:BLC4,0,0,5,3,02.通过矩形中心定位点生成面BLC5,定位点 X 方向坐标 XCENTER,定位点 Y 方向坐标 YCENTER,矩形宽度 WIDTH,矩形高度HEIGHT,矩形深度 DEPTH注:与上条命令的不同就在于矩形的定位点不一样例:BLC5,2.5,1.5,5,3,03. 通过在工作平面定义矩形X.Y 坐标生成面RECTNG,矩形左边界 X 坐标 X1,矩形右边界X 坐标 X2, 矩形下边界Y坐标 Y1, 矩形上边界Y 坐标Y2例:RECTNG,0,5,0,3生成圆面4.通过中心定位点生成实心圆面CYL4,定位点 X 方向坐标 XCENTER,定位点 Y 方向坐标 YCENTER,圆面的内半径 RAD1,内圆面旋转角度 THETA1,圆面的外半径 RAD2,外圆面旋转角度THETA2,圆面的深度 DEPTH注:如要实心的圆面则不用 RAD2,THETA2,DEPTH例:CYL4,0,0,5,3605.生成扇形圆面命令介绍如上例 1 实心扇形 :CYL4,0,0,5,60例2 扇形圆环 :CYL4,0,0,5,60,10,60例3 整的圆环 :CYL4,0,0,5,360,10,360注: 同时可通过定义圆面的深度以生成柱体6.通过在工作平面定义起始点生成圆面CYL5,开始点 X坐标 XEDGE1,开始点 Y坐标 YEDGE1,结束点 X坐标 XEDGE2,结束点 Y坐标YEDGE2, 圆面深度 DEPTH例:CYL5,0,0,2,2,7.通过在工作平面定义内外半径和起始角度来生成圆面PCIRC,内半径 RAD1,外半径 RAD2,起始角度THETA1,结束角度THETA2例LCIRC,2,5,30,1808.生成面与面的倒角AFILLT, 面 1 的编号 NA1, 面 2 的编号 NA2,倒角半径RAD例:AFILLT,2,5,2下一讲 : 多边形面的生成第五天目标 : 掌握多边形面和体的生成1.生成多边形面命令 :RPR4, 多边形的边数NSIDES,中心定位点X 坐标 XCENTER,中心定位点Y 坐标 YCENTER, 中心定位点距各边顶点的距离RADIUS,多边形旋转角度THETA例:RPR4,4,0,0,0.15,30注: 这条命令可通过定义不同的NSIDES生成三边形 , 四边形 ,...,八边形2.生成多边形体命令 :RPR4, 多边形的边数NSIDES,中心定位点X 坐标 XCENTER,中心定位点Y 坐标 YCENTER, 中心定位点距各边顶点的距离RADIUS,多边形旋转角度THETA,多边形的深度DEPTH例:RPR4,4,0,0,0.15,30,0.1注: 多边形体和面命令唯一的不同就在于深度DEPTH的定义到此 , 关键点 , 线 , 面的生成讲解已结束, 下一讲 : 体的生成第六天目标 : 掌握体的生成命令1.通过关键点生成体命令 :V, 关键点 P1, 关键点 P2, P3, P4, P5, P6, P7, P8例:V,4,5,6,7,15,24,252.通过面生成体命令 :VA, 面 A1, 面 A2, A3, A4, A5, A6, A7, A8, A9, A10例:VA,3,4,5,8,103.通过长方形角上定位点生成体命令 :BLC4该命令前面在讲生成面的时候已作介绍, 唯一的不同在于深度DEPTH的定义 .4.通过长方形中心定位点生成面命令 :BLC55.通过定义长方体起始位置生成体命令 :BLOCK,开始点 X 坐标 X1,结束点 X 坐标 X2, Y1, Y2, Z1, Z2例:BLOCK,2,5,0,2,1,36.生成圆柱体基本命令通生成圆形面, 不同在于DEPTH的定义基本命令 :CYL4基本命令 :CYL5基本命令 :CYLIND7.生成棱柱基本命令通生成多边形, 不同在于DEPTH的定义基本命令 :RPR48.通过球心半径生成球体命令 :SPH4, 球心 X 坐标 XCENTER,球心 Y 坐标 YCENTER,半径 RAD1,半径 RAD2例:SPH4,1,1,2,59.通过直径上起始点坐标生成球体命令 :SPH5, 起点 X 坐标 XEDGE1,起点 Y坐标 YEDGE1,结束点 X 坐标 XEDGE2,结束点 Y 坐标YEDGE2例:SPH5,2,5,7,610.在工作平面起点通过半径和转动角度生成球体命令 :SPHERE,半径 RAD1,半径 RAD2,转动角度THETA1,转动角度THETA2例:SPHERE,2,5,0,6011.生成圆锥体命令 :CONE,底面半径 RBOT,顶面半径 RTOP,底面高 Z1, 顶面高 Z2, 转动角度 THETA1,转动角度THETA2例:CONE,10,20,0,50,0,180下一讲 : 布尔操作第七天目标 : 掌握常用的布尔操作命令1.沿法向延伸面生成体命令 :VOFFST,面的编号NAREA,面拉伸的长度DIST, 关键点增量KINC例:VOFFST,1,2,,2.通过坐标的增量延伸面生成体命令 :VEXT, 面 1 的编号 NA1,面 2 的编号 NA2,增量 NINC,X 方向的增量 DX,Y 方向的增量 DY,Z 方向的增量 DZ, RX, RY, RZ例:VEXT,1,5,1,1,2,2,3. 面绕轴旋转生成体命令 :VROTAT,面 1 的编号PAX1,定位轴关键点 2 编号NA1,面 2 的编号PAX2,旋转角度NA2,NA3, NA4, NA5, NA6,ARC,生成体的段数 NSEG定位轴关键点 1 编号例:VROTAT,1,2,,,,,4,5,360,44. 沿线延伸面生成体命令 :VDRAG,面 1 的编号 NA1,面 2 的编号 NA2, NA3, NA4, NA5, NA6,导引线 2 的编号 NLP2, NLP3, NLP4, NLP5, NLP6例:VDRAG,2,3,,,,,8,导引线 1 的编号NLP1, 5. 线绕轴旋转生成面命令 :AROTAT,线 1 的编号 NL1, NL2, NL3, NL4, NL5, NL6, 位轴关键点 2 的编号 PAX2,旋转角度ARC,生成面的段数例:AROTAT,3,4,,,,,6,8,360,4定位轴关键点NSEG1 的编号PAX1,定6.沿线延伸线生成面命令 :ADRAG,线 1 的编号 NL1,NL2, NL3, NL4, NL5, NL6, 导引线 1 的编号 NLP1, NLP2, NLP3, NLP4, NLP5, NLP6例:ADRAG,3,,,,,,87.同理可以延伸关键点 , 相应的命令如下 :LROTAT, NK1, NK2, NK3, NK4, NK5, NK6, PAX1, PAX2, ARC, NSEGLDRAG, NK1, NK2, NK3, NK4, NK5, NK6, NL1, NL2, NL3, NL4, NL5, NL6各选项的含义雷同于上 .8.延伸一条线命令 LEXTND,线的编号 NL1, 定位关键点编号 NK1,延伸的距离 DIST, 原有线是否保留控制项KEEP例LEXTND,5,2,1.5,09.布尔操作 : 加命令 LCOMB,线编号 NL1,线编号 NL2, 是否修改控制项KEEP例LCOMB,2,5注: 对面和体的相应为:VADD,AADD.选项的含义都类似10.布尔操作 : 粘接和搭接搭接的核心关键字为:OVLAP,随实体的不同略有不同, 如 :对体为 VOVLAP对面为 AOVLAP对线为 LOVLAP粘接的核心关键字为:GLUE, 随实体的不同略有不同, 如 :对体为 VGLUE对面为 AGLUE对线为 LGLUE但其他的选项的含义是类似的, 这里就不再累述.下一讲 : 移动 , 复制 , 映射 , 删除 ...第八天目标 : 掌握体素的移动, 复制 , 删除 , 映射一. 移动关键点命令 :KMODIF,关键点编号NPT,移动后的坐标X, 移动后的坐标Y, 移动后的坐标Z例:KMODIF,5,0,0,2二. 移动复制关键点命令 :KGEN,复制次数选项 ITIME, 起始关键点编号 NP1,结束关键点编号 NP2,增量 NINC,偏移DX,偏移 DY,偏移 DZ,关键点编号增量 KINC, 生成节点单元控制项 NOELEM,原关键点是否被修改选项 IMOVE例:KGEN,2,1,10,1,2,2,2,,,,注:IMOVE 选项说明 , 设置为 0 时, 不修改原关键点 , 即为复制 , 设置为 1 时 , 修改原关键点 , 即为移动 , 从而通过控制 IMOVE选项实现移动或复制 .三. 移动复制线命LGEN,ITIME,NL1,NL2,NINC,DX,DY,DZ,KINC,NOELEM,IMOVE各选项的含义同上四. 移动复制面命:AGEN,ITIME,NA1,NA2,NINC,DX,DY,DZ,KINC,NOELEM,IMOVE各选项的含义同上五. 移动复制体命令 :VGEN,ITIME,NV1,NV2,NINC,DX,DY,DZ,KINC,NOELEM,IMOVE各选项的含义同上六. 修改面的法向方向命令 :ANORM,面的编号 ANUM,单元的法向方向是否修改选项NOEFLIP例:ANORM,2七.体素的删除基本的命令为 :*DELE组合不同的关键字形成不同的命令如:KDELE,LDELE,ADELE,VDELE基本的命令格式为 :*DELE, 起始体素编号 N*1, 结束体素编号 N*2, 增量 NINC,是否删除体素下层的元素选项 KSWP如LDELE,2,5,1,1八.体素的映射基本的命令为 :*SYMM组合不同的关键字形成不同的命令如:KSYMM,LSYMM,ARSYM,VSYMM基本的命令格式为 :*SYMM,映射轴选项 NCOMP,起始体素编号 N*1, 结束体素编号 N*2, 增量NINC,关键点编号增量 KINC,NOELEM, IMOVE如:VSYMM,X,1,10,1,,,,。

ANSYS命令流(入门必备)

ANSYS命令集/EXIT,Slab,Fname,Ext,Dir Slab=ALL 保存所有资料Slab=NOSA VE所有更改资料不保存Slab=MODEL保存实体模型，有限元模型，负载的资料（系统默认）例：/EXIT，ALL-------------------------------------------------------- /FILNAM，Fname Fname=工作文件名称，不要扩展名例：/FILNAM，Sanpangzi--------------------------------------------------------/SA VE,Fname,Ext,Dir 保存目前所有的Datebase资料，即更新Jobname.db--------------------------------------------------------/RESUME,Fname,Ext,Dir,NOPAR 回到最后SA VE时的Datebase状态--------------------------------------------------------/CLEAR 清除所有Datebase资料-------------------------------------------------------- LOCAL，KCN，KCS，XC，YC，ZC，THXY，THYZ，THZX，PAR1，PAR2 定义区域坐标系统KCN 区域坐标系统代号，大于10的任何号码KCS=0，1，20=笛卡儿坐标1=圆柱坐标2=球面坐标XC，YC，ZC 该区域坐标原点与整体坐标原点的关系THXY，THYZ，THZX 该区域坐标与整体坐标XYZ轴的关系例：LOCAL，11，1，1，1，0-------------------------------------------------------- CSYS，0，1，2声明当前坐标系统例：CSYS，0-------------------------------------------------------- /UNITS，LABEL 声明系统分析时所用的单位LABEL=SI (米，千克，秒)LABEL=CGS (厘米，克，秒)LABEL=BFT (英尺)LABEL=BIN (英寸)例：LABEL，SI-------------------------------------------------------- /PREP7进入通用前处理器-------------------------------------------------------- N，NODE，X，Y，Z，THXY，THYZ，THZX 定义节点NODE 节点号码X，Y，Z 节点在当前坐标系中位置例：N，1，2，3，4-------------------------------------------------------- NDELE，NODE1，NODE2，NINC 删除已建立的节点NODE1，NODE2 删除从NODE1到NODE2的节点，如1到100 NINC 间隔号码，1为1到100全删，2为1，3，5 (99)例：NDELE，1，100，2--------------------------------------------------------NPLOT，KNUM 将节点显示在图形窗口中KNUM=0不显示节点号码KNUM=1显示节点号码--------------------------------------------------------NLIST，NODE1，NODE2，NINC将节点资料列在窗口中例：NLIST--------------------------------------------------------NGEN，ITIME，INC，NODE1，NODE2，NINC，DX，DY，DZ，SPACE 复制节点ITIME 复制次数，包括本身INC复制时节点号码增量NODE1，NODE2，NINC 要复制的节点DX，DY，DZ 复制出的节点的位置改变量例：NGEN，4，5，1，5，1，1，2，3 将节点1到5复制4次，每次复制X，Y，Z方向分别移动1，2，3单位长度-------------------------------------------------------- FILL，ITIME，INC，NODE1，NODE2，NINC，DX，DY，DZ，SPACE 填充节点（默认为均分填充）例：FILL，1，100 在节点1到100之间填充2，3 (99)-------------------------------------------------------- ET，ITYPE，Ename，KOPT1……KOPT6，INOPR 定义元素ITYPE 元素类型编号Ename 所使用元素名称KOPT1-KOPT6 元素特性编码例：ET，1，LINK1 第1类元素为LINK1单元-------------------------------------------------------- MP，Lab，MAT，C0，C1，C2，C3，C4 定义材料特性材料特性为固定值，其值为C0材料特性随温度变化，由C1-C4控制Lab 材料特性类别MAT 对应ET所定义的元素类型编号ITYPELab=EX，EY，EZ 杨氏系数Lab=DENS 密度Lab=PRXY，PRYZ，PRZX 泊松比Lab=GXY，GYZ，GZX 剪力模数例：MP，EX，1，207E9 第一类元素的杨氏系数为207E9-------------------------------------------------------- R，NSET，R1……R6 定义元素类型几何特性NSET 属性组别号码（系统默认值1）R1-R6 所定义元素类型几何特性值例：R，1，1E-4，2.09E-10 ，0.005-------------------------------------------------------- E，I，J，K，L，M，N，O，P定义元素连接方式I-P 定义元素节点顺序的号码例：E，1，2，5，7 四节点元素的节点顺序为1，2，5，7-------------------------------------------------------- EPLOT，KNUM 将元素显示在图形窗口中ENUM=0 不显示元素ENUM=1 显示元素-------------------------------------------------------- ELIST 将元素资料列在窗口中-------------------------------------------------------- EDELE，IEL1，IEL2，INC 删除已建立的元素IEL1，IEL2，INC 欲删除元素的范围例：EDELE，1，10，1-------------------------------------------------------- EGEN，ITIME，NINC，IEL1，IEL2，IEINC，MINC，IINC，RINC，CINC，DX，DY，DZ 复制元素ITIME 复制次数，包括本身NINC复制时节点号码增量IEL1，IEL2，IEINC 欲复制的元素范围DX，DY，DZ 复制出的元素的位置改变量例：EGEN，6，12，1，4，1 将元素1到4复制6次-------------------------------------------------------- /PNUM，Label，KEY 在图形中显示号码Label=NODE，ELEM，KP，LINE，AREA，VOLUKEY=0 不显示号码KEY=1 显示号码例：/PNUM，ELEM，1-------------------------------------------------------- /SOLU 进入解题处理器--------------------------------------------------------ANTYPE，Antype，Status 声明分析类型Antype=STATIC or 0 静态分析（系统默认）Antype=BUCKLE or 1 屈曲分析Antype=MODAL or 2 振动模态分析Antype=HARMIC or 3 调和外力动力系统分析Antype=TRANS or 4 瞬时动力系统分析例：ANTYPE，STATIC-------------------------------------------------------- F，NODE，Lab，V ALUE，V ALUE2，NEND，NINC定义节点上的集中力NODE 节点号码Lab 外力形式Lab=FX，FY，FZ，MX，MY，MZ 结构力学Lab=HEAT 热学的热流量Lab=AMP，CHRG 电学的电流，电荷Lab=FLUX 磁学的磁通量V ALUE 外力大小NODE，NEND，NINC 施力节点范围例：F，1，FY，，-200，5，1 =F，ALL，FY，-200 节点1-5 的Y方向定义集中力-200（注意FY，，表明V ALUE2默认）-------------------------------------------------------- FDELE，NODE，Lab，NEND，NINC 删除节点集中力例：FDELE，1，FY，5，1 = FDELE，ALL-------------------------------------------------------- D，NODE，Lab，V ALUE，V ALUE2，NEND，NINC，Lab2，Lab3，Lab4，Lab5，Lab6 定义节点自由度的限制NODE，NEND，NINC 选取自由度约束节点的范围Lab 相对元素每一个节点受自由度约束的形式结构力学Lab=UX，UY，UZ（直线位移）Lab=ROTX，ROTY，ROTZ（旋转位移）例：D，1，UX，，，5，1 节点1到5 X方向约束D，1，UX，，，5，1，UY 节点1到5 X Y方向约束D，1，ALL，，，5，1 节点1到5 全部自由度约束注意：使用命令前要先定义元素--------------------------------------------------------DDELE，NODE，Lab，NEND，NINC 取消节点自由度约束例：DDELE，ALL-------------------------------------------------------- DLIST，NODE1，NODE2，NINC 列出节点自由度约束--------------------------------------------------------DL，LINE，AREA，Lab，V ALUE1，V ALUE2定义线自由度限制LINE 线号AREA 线所属面积例：DL，8，3，ALL 定义面积3上面线8的约束注意：同时有DLLIST，DLDELE命令--------------------------------------------------------SFBEAM，ELEM，LKEY，Lab，V ALI，V ALJ，V AL2I，V AL2J，IOFFST，JOFFST定义分布力作用于梁元素的方式及大小ELEM 分布力所作用的元素编号LKEY 梁元素的4个面中分布力所在面号码Lab=PRES （表示分布压力）V ALI，V ALJ 在I，J点分布力的值例：SFBEAM，1，1，PRES，60，30 元素1上1号面作用分布力SFBEAM，1，1，PRES，-30，60--------------------------------------------------------SFE，ELEM，LKEY，Lab，KV AL，V AL1，V AL2，V AL3，V AL4定义分布力作用于元素上的方式及大小ELEM 分布力所作用的元素编号LKEY 分布力作用边，面的号码Lab=PRES （表示分布压力）V AL1-V AL4 分布力在元素边，面上节点的值例：SFE，4，2，PRES，，20，60元素4的第2边，面作用分布力--------------------------------------------------------SF，NLIST，Lab，V ALUE，V ALUE2 定义节点间分布力NLIST 分布力作用边或面上所有节点Lab=PRES例：SF，ALL，PRES，10注意：SFE适用于非均匀分布力，作用在元素的边或面上SF适用于均匀分布力，作用在节点之间-------------------------------------------------------- SFLIST，NODE，Lab 显示分布力-------------------------------------------------------- NSEL，Type，Item，COMP，VMIN，VMAX，VINC，KABS节点选择命令Type 选择方式Type=S 选择一组节点为ACTIVE点Type=R 在现有ACTIVE点中选出一部分作为ACTIVE点Type=A 在选择一部分节点，加入ACTIVE点中Type=U 在现有ACTIVE点中，排除某些节点Type=ALL 选择全部节点作为ACTIVE点Item =NODE 用节点号码选取┇Item =LOC 用节点坐标选取COMP=无┇COMP=X 以节点X坐标为准VMIN，VMAX，VINC 节点选取范围例：NSEL，S，NODE，，1，13，1选1-13中奇数点为ACTIVE点NSEL，A，NODE，，14，20，1 选14-20加入ACTIVE点--------------------------------------------------------OUTPR，Item，FREQ，Cname 分析结果是否显示于输出窗口Item=ALL 所有结果Item=NSOL 节点自由度结果FREQ 负载的次数FREQ=ALL 最后负载例：OUTPR，ALL，ALL注意：仅用于小程序--------------------------------------------------------SOLVE 开始解题--------------------------------------------------------/POST1 进入后处理器--------------------------------------------------------PLDISP，KUND 显示结构变形结构KUND=0 显示变形后结构形状KUND=1 显示变形前后结构形状KUND=0 显示变形前后结构形状，但仅显示外观--------------------------------------------------------PLESOL，Item，Comp显示元素的解答Item（何种解答）Comp（Item分量）S X，Y，Z，XY，YZ，ZX 应力S 1，2，3 主应力S EQV，INT 等效应力F X，Y，Z 结构力M X，Y，Z 结构力矩例：PLESOL，S，X，Y，Z--------------------------------------------------------PLNSOL，Item，Comp显示节点的解答Item（何种解答）Comp（Item分量）S X，Y，Z，XY，YZ，ZX 应力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 旋转位移例：PLNSOL，S，Y PLNSOL，U，X--------------------------------------------------------PRESOL，Item，Comp 打印元素解答Item（何种解答）Comp（Item分量）S X，Y，Z，XY，YZ，ZX 应力F X，Y，Z 结构力M X，Y，Z 结构力矩例：PRESOL，S，X--------------------------------------------------------PRNSOL，Item，Comp 打印节点解答Item（何种解答）Comp（Item分量）U X，Y，Z 位移U X，Y，Z方向及总向量方向位移S COMP 应力S PRIN 主应立，等效应力例：PRNSOL，U PRNSOL，S注意：查看结果通常使用PLDISP，1 PLNSOL，U，Y……PRNSOL，S其中PLNSOL中Comp不能省略，PRNSOL中可以省略--------------------------------------------------------TYPE，Itype 声明建立元素时，元素形式号码（对应ET的Itype）例：ET，1，LINK1 ET，2，PLANE42TYPE，1建立LINK1元素TYPE，2建立PLANE42元素--------------------------------------------------------REAL，NSET 声明建立元素时，元素几何参数属性编号（对应R，NSET）--------------------------------------------------------MAT，MAT 声明建立元素时，元素材料特性属性编号（对应MP，MAT）--------------------------------------------------------LSWRITE，LSNUM 多重负载资料保存至文件Jobname.S0i 例：LSWRITE 自动编号，不用输入--------------------------------------------------------LSSOLVE，SLMIN，LSMX，LSINC读取所定义的多重负载并解答SLMIN，LSMX，LSINC 读取负载的范围例：LSSOLVE，1，4 获得1-4负载的解答--------------------------------------------------------DDELE，NODE，Lab，NEND，NINC删除定义的约束条件NODE，NEND，NINC 删除约束的节点范围Lab删除约束的节点方向例：DDELE，1，UY，5，2 删除1，3，5节点的Y向约束--------------------------------------------------------FDELE，NODE，Lab，NEND，NINC删除定义的集中力NODE，NEND，NINC 删除集中力的范围Lab删除外力的方向例：FDELE，1，FY，5，2 删除1，3，5节点的Y向集中力--------------------------------------------------------SFDELE，Nlist，Lab 将已定义的面载荷删除Nlist 面负载所含的节点（由NSEL选择，设Nlist=ALL）Lab=PRES （结构力学）--------------------------------------------------------SFEDELE，ELEM，LKEY，Lab 将负载从元素上删除ELEM 元素编号LKEY 负载作用边，面的号码Lab=PRES （结构力学）--------------------------------------------------------SET，Lstep，SBSTEP，FACT，KIMG，TIME，NGLE，NSET检查负载结果Lstep=负载编号例：SET，2 检查第二负载的结果PLDISP，1……注意：此命令要在/POST1中使用-------------------------------------------------------- FILE，Fname，Ext，Dir 读取分析后的结果文件例：FILE，TEST，RST--------------------------------------------------------ANTYPE，Antype，Status 声明分析类型例：ANTYPE=MODAL or 2 模态分析-------------------------------------------------------- MODOPT，Method，NMODE，FREQB，FREQE，PRMODE，NUMKEY 选择模态分析方法Method=REDUC 降阶法Method=SUBSP 次空间法NMODE 欲求模态个数（降阶法小于主自由度一半）FREQB，FREQE 欲探讨振动频率范围（默认全部）PRMODE 分析后，模态结果保存到结果文件的个数例：MODOPT，SUBSP，5-------------------------------------------------------- M，NODE，Lab1，NEND，NINC，Lab2。

ansys命令流使用方法

ansys命令流使用方法
在ANSYS中，命令流是一种用于执行特定操作的自动化工具。

以下是ANSYS命令流使用的一般步骤：
1. 打开ANSYS软件并加载您要使用的工程文件。

2. 在ANSYS Graphical User Interface (GUI) 中，将鼠标指针放
在工具栏上。

在“Run”下拉菜单中选择“Command Line”。

3. 在命令行窗口中，输入和编辑您想执行的命令。

您可以使用ANSYS的命令语言以及相关命令进行模型操作、网格生成、
求解等。

4. 您可以通过多种方式输入命令：直接在命令行中输入、从脚本文件中读取、从ANSYS GUI中的日志文件中复制粘贴等。

5. 您可以使用命令流中的参数和变量来进行自动化操作。

使用“!VARIABLE”语句定义变量，并通过“!VARIABLE = value”语
句赋值。

6. 使用ANSYS的各种功能命令对模型进行操作。

例如，在预
处理阶段，您可以使用命令生成几何体、定义材料属性、设定网格、添加边界条件等。

7. 在求解阶段，使用命令启动求解器，设置求解器选项，运行求解器，并监视求解器的输出。

8. 在结果后处理阶段，使用命令读取并处理结果数据，生成图形、报告等。

9. 执行命令流，您可以一次性执行整个命令流，或者逐个执行命令。

10. 您还可以将命令流保存为脚本文件，以便将来再次使用。

以上是ANSYS命令流的一般用法，具体的命令和语法取决于您的特定需求和ANSYS的版本。

建议您参考ANSYS的官方文档和教程，以获得更详细和准确的使用说明。

ansys隧道内力命令流

wprota,,,-6 !工作平面绕Y轴旋转-6度
asbw,all !用工作平面切割所选择的所有面
*enddo
asel,s,,,14,15,1 !选择面
esize,
amesh,11 !划分11号面
!划分围岩网格
!设置网格份数
lsel,s,,,5,11,1 !选择线
lsel,a,,,13,15,1
adele,1 !删除面
NUMCMP,AREA !压缩面编号
!生成锚杆
asel,s,,,14 !选择面14
l,77,74
l,75,68
l,74,69
l,73,70
l,71,64
l,70,65
l,69,66
l,67,60
l,66,61
l,65,62
a,60,63,79,76 !通过4个关键点生成一个面积区域
asbl,all,all !通过线分割面生成新面
*set,t1,0.25
*set,r_karst,3.6 !溶洞半径
*set,r1,2.5E-2
*afun,deg
!定义单元类型
ET,1,BEAM3 !定义衬砌支护单元
type,3
lmesh,all
allsel
!划分开挖掉土体单元网格
mat,3
type,2
mshkey,0 !设定自由网格划分
mshape,0 !设定四边形网格划分
asel,a,,,19,20,1
asel,a,,,22,25,3
asel,a,,,26,28,2
asel,a,,,31,32,1
asel,a,,,34,50,2

ANSYS实例分析连拱隧道开挖命令流

!III级围岩，二维，全断面开挖/cle/filname,ch-three-2D,1!几何参数!计算外围区域*set,x1,-60*set,x2,60*set,y1,-30*set,y2,50*set,l1,3.5 !锚杆长度*set,r1,0.0125 !锚杆半径*set,th1,0.15 !初次衬砌厚度*set,th2,0.05 !中导洞喷层厚!材料参数!中隔墙*set,ec25,29.5e9*set,uc25,0.2*set,dens25,2500!喷射混凝土*set,ec20,21e9 *set,uc20,0.2*set,dens20,2200!锚杆*set,el,170e9*set,ul,0.3*set,densl,7800!围岩*set,er,6e9*set,ur,0.28*set,densr,2500 *set,cr,0.8e6*set,fir,40/prep7!定义单元类型et,1,plane42keyopt,1,3,2 !平面应变问题et,2,beam3keyopt,2,6,1 !打开输出内力选项et,3,link1!定义实常数r,1,th1,th1/12,th1r,2,th2,th2/12,th2r,3,3.1415926*r1**2!几何建模!建中隔墙k,1,-3,-1.69k,2,-3,-0.89k,3,-3,-0.37k,4,-1.33552,-0.37k,5,-0.95134,1.57593k,6,-1.25121,3.53662k,7,1.25121,3.53662k,8,0.95134,1.57593k,9,1.33552,-0.37k,10,3,-0.37k,11,3,-0.89k,12,3,-1.69l,1,2l,2,3l,3,4larc,4,6,5larc,6,7,5,3larc,7,9,8l,9,10l,10,11l,11,12l,12,1!左侧围岩及锚杆圈k,13,-3,0.81k,14,-2.52524,2.42963 k,15,-1.483023,4.11778 k,16,-1.77650,4.67036 k,17,-2.058364,5.09399k,18,-2.37658,5.49104 k,19,-2.734606,5.864085 k,20,-3.12567,6.20235 k,21,-5.88951,7.72342 k,22,-9,8.25k,23,-12.11049,7.72342 k,24,-14.8743,6.20235 k,25,-16.57763,3.97894 k,26,-17.02812,1.21456 k,27,-16.9,-0.89k,28,-20.4,-0.89k,29,-20.51567,0.91966 k,30,-19.7906,5.36705 k,31,-17.05,8.944k,32,-13.26253,11.02838 k,33,-9,11.75k,34,-4.146154,10.80594 k,35,0,8.11142k,36,1.143324,6.85081!导洞左侧圈l,3,13 larc,13,6,14 !左侧开挖圈larc,6,16,15 larc,16,18,17 larc,18,20,19 larc,20,22,21 larc,22,24,23 larc,24,26,25 l,26,27!锚杆圈l,28,29 larc,29,31,30 larc,31,33,32 larc,33,35,34 l,35,36!内圈线合并lcomb,15,16 lcomb,15,17 lcomb,15,18 lcomb,15,19!外圈线合并lcomb,20,21lcomb,20,22lcomb,20,23numcmp,allldiv,14,,,2ldiv,15,,,36ldiv,16,,,36ldiv,17,,,2kdele,19,22,1kdele,25,27,1numcmp,allkdele,17numcmp,all!锚杆线（40根，含建3维时可用的另外20根锚杆）l,16,21l,22,93l,17,20*do,i,23,57,1*enddol,18,19l,18,2numcmp,all!内层细化区域k,100,-16.9,-1.69 k,101,-20.4,-1.69 k,102,-24.9,-1.69 k,103,-24.9,-0.89 k,104,-24.9,16.25625 k,105,,16.25625k,106,x1,-1.69k,107,x1,16.25625 k,108,-24.9,y2k,109,-24.9,y1l,1,100l,100,101l,101,102l,102,103l,104,105l,19,103l,70,104l,18,100l,19,101l,102,109l,102,106l,104,107l,104,108!左侧线镜像至右并建立各面allslsymm,x,11,143,1rectng,x1,x2,y1,y2asbl,1,allnummrg,all,,,low numcmp,alllcomb,287,288lcomb,289,290lcomb,294,295 lcomb,296,297 lcomb,130,292 numcmp,allalls!定义材料属性!岩体*do,i,1,6,1mp,ex,i,ermp,prxy,i,ur mp,dens,i,densr tb,dp,itbdata,1,cr,fir*enddo!中导洞及隧道初支*do,i,7,10,1mp,ex,i,ec20mp,prxy,i,uc20 mp,dens,i,dens20*enddo!锚杆左右侧*do,i,11,12,1 mp,ex,i,elmp,prxy,i,ul mp,dens,i,densl*enddo!中隔墙mp,ex,13,ec25 mp,prxy,13,uc25 mp,dens,13,dens25!中导洞初支!左侧type,2mat,7real,2lsel,s,line,,1,2,1lesize,1,,,1 lesize,2,,,1 lesize,11,,,2 lesize,12,,,6 lmesh,all!右侧mat,8real,2lsel,s,line,,8,9,1 lsel,a,line,,259,260,1 lsel,a,line,,5 lesize,9,,,1 lesize,8,,,1 lesize,260,,,2 lesize,259,,,6 lesize,5,,,4 lmesh,all!左侧隧道初支mat,9real,1lsel,s,line,,13,15,1lsel,a,line,,288 lesize,13,,,3 lesize,all,,,1 lmesh,all!右侧隧道初支mat,10real,1lsel,s,line,,139,140,1 lsel,a,line,,142,176,1 lsel,a,line,,286 lsel,a,line,,292 lesize,286,,,3 lesize,all,,,1 lmesh,all!左隧道锚杆type,3mat,11real,3lsel,s,line,,89,123,2 lsel,a,line,,263,269,3lsel,a,line,,210 lsel,a,line,,266 lesize,263,,,3 lesize,269,,,3 lesize,271,,,1 lesize,279,,,1 lesize,210,,,1 lesize,266,,,1 lesize,all,,,5,2 lmesh,all!右隧道锚杆type,3mat,12real,3lsel,s,line,,212,246,2 lsel,a,line,,277,280,3 lsel,a,line,,265 lsel,a,line,,270lsel,a,line,,87 lsel,a,line,,17 lesize,277,,,3 lesize,280,,,3 lesize,265,,,1 lesize,270,,,1 lesize,87,,,1 lesize,17,,,1 lesize,all,,,5,2 lmesh,all!中墙type,1mat,1lsel,s,line,,1,10,1 lesize,3,,,2 lesize,4,,,8 lesize,6,,,8 lesize,7,,,2 lesize,10,,,8 amesh,97!左中导洞mat,2 amesh,2!右中导洞mat,3 amesh,3!左隧道内岩体mat,4lsel,s,line,,124 lesize,124,,,20 amesh,101!右隧道内岩体mat,5lsel,s,line,,261 lesize,261,,,20 amesh,102!导洞上方围岩mat,6lsel,s,line,,264,267,3lsel,a,line,,281lsel,a,line,,268,278,10lsel,a,line,,272lesize,264,,,3lesize,278,,,3lesize,all,,,1!左侧锚杆加固圈lsel,s,line,,88,122,2 !还有264,267,281为扩展到3维的锚杆lesize,all,,,5,2lsel,s,line,,53,86,1lsel,a,line,,16lsel,a,line,,287lesize,all,,,1!右侧锚杆加固圈lsel,s,line,,211,245,2 !还有278,87,17为扩展到3维的锚杆lesize,all,,,5,2 lsel,s,line,,177,209,1 lsel,a,line,,291,294,3 lsel,a,line,,141 lesize,all,,,1!小外围上部lsel,s,line,,64,86,1 lsel,a,line,,188,209,1 lsel,a,line,,294 lccat,alllsel,a,line,,130,132,1 lsel,a,line,,252 lesize,130,,,46 lesize,132,,,8,2 lesize,252,,,8,2 amesh,103!小外围左lsel,s,line,,53,63,1 lsel,a,line,,287lccat,alllsel,s,line,,129 lsel,a,line,,131 lesize,129,,,13 lesize,131,,,8,2 amesh,99!小外围右lsel,s,line,,291 lsel,a,line,,177,187,1 lsel,a,line,,141 lccat,alllsel,a,line,,250 lesize,250,,,13 lsel,a,line,,251 lesize,251,,,8,2 amesh,100!小外围下lsel,s,line,,128lsel,s,line,,127 lesize,127,,,8,2lsel,s,line,,133,134,1 lesize,all,,,1lsel,s,line,,126 lesize,126,,,5,2lsel,s,line,,125 lesize,125,,,20lsel,s,line,,262 lesize,262,,,20lsel,s,line,,253,254,1 lesize,all,,,1lsel,s,line,,247 lesize,247,,,5,2lsel,s,line,,249 lesize,249,,,1lsel,s,line,,248 lesize,248,,,8,2 amesh,42,43 amesh,1,4,3 amesh,82,83!大外围网格划分lsel,s,line,,274,276,2 lsel,a,line,,136,137,1 lsel,a,line,,256,257,1 lsel,a,line,,282,285,3 lesize,all,,,18,2.5lsel,s,line,,273,283,10 lsel,a,line,,135,255,120 lesize,all,,,15,2.5lsel,s,line,,275,284,9 lsel,a,line,,138,258,120 lesize,all,,,18,2.5lsel,s,line,,290,295,5lesize,290,,,46 lesize,295,,,74lsel,s,line,,289,293 lesize,all,,,14lsel,s,line,,128,129 lccat,alllsel,s,line,,249,250 lccat,alllsel,s,line,,125,127,1 lsel,a,line,,10lsel,a,line,,262 lsel,a,line,,247,248,1 lccat,allamesh,allallslsel,s,lccatldele,all numcmp,allallselsave,data!施加约束allsnsel,s,loc,x,x1 nsel,a,loc,x,x2 d,all,uxnsel,s,loc,y,y1d,all,uyallsacel,,10fini/soluantype,static pred,onlnsrch,on nlgeom,on nropt,full outress,all,all deltim,0.1,0.05,0.2!荷载步1time,1esel,s,type,,2,3,1ekill,allesel,allesel,s,livensle,snsel,inved,all,all,0nsel,allesel,all! iswrite,1 !生成初应力文件solvesave,s1,db!荷载步2time,2esel,s,mat,,2,3,1ekill,allesel,s,mat,,7,8,1ealive,all nsle,sddele,all,all esel,allesel,s,live nsle,snsel,inved,all,allesel,s,mat,,1 mpchg,13,all allssolvesave,s2,db!荷载步3 time,3esel,s,mat,,4,7,3 ekill,allesel,s,mat,,9,11,2 esel,a,mat,,3 ealive,all nsle,sddele,all,all esel,allesel,s,livensle,snsel,inved,all,allallssolvesave,s3,db!荷载步4time,4esel,s,mat,,5,8,3 esel,a,mat,,3 ekill,allesel,s,mat,,10,12,2 ealive,allnsle,sddele,all,all esel,allesel,s,livensle,snsel,inved,all,allallssolvesave,s4,dbfini/post1!载荷步1结果resume,s1,dbfile,ch-three-2D,rst/pbc,all,,0 !不显示边界约束等条件set,1,lastesel,s,liveplnsol,u,sum,0,1.0plnsol,u,xplnsol,u,yplnsol,s,1plnsol,s,3plnsol,s,eqvplnsol,eppl,eqv,0,1!载荷步2结果resume,s2,db/pbc,all,,0esel,s,liveesel,u,mat,,13plnsol,u,sum,0,1 !未减去自重引起的合位移云图plnsol,u,x !未减去自重引起的位移plnsol,u,y !未减去自重引起的位移plnsol,s,1plnsol,s,3plnsol,s,eqvplnsol,eppl,eqv,0,1esel,s,mat,,13plnsol,u,yplnsol,s,1plnsol,s,3plnsol,s,eqv!中导洞初支内力esel,s,liveesel,r,ename,,3 !选择beam3单元类型etable,fx_i,smisc,1etable,fx_j,smisc,7etable,fy_i,smisc,2etable,fy_j,smisc,8etable,mz_i,smisc,6etable,mz_j,smisc,12plls,mz_i,mz_j,-0.2 !显示弯矩plls,fx_i,fx_j,0.2 !显示轴力plls,fy_i,fy_j,0.2 !显示剪力!查看位移场file,ch-three-2D,rstset,2,lastesel,s,liveesel,u,mat,,13lcdef,1,1 !定义第1个载荷步为工况1lcoper,sub,1 !当前载荷步的结果减去工况1的结果(位移结果有问题）esel,s,liveplnsol,u,sum !显示合位移云图plnsol,u,xplnsol,u,y!载荷步3结果resume,s3,db/pbc,all,,0esel,s,liveesel,u,mat,,13plnsol,s,eqvplnsol,s,1plnsol,s,3plnsol,eppl,eqv,0,1!中墙等效应力场esel,s,mat,,13plnsol,s,eqvplnsol,s,1plnsol,s,3!中导洞初支及左隧道衬砌内力esel,s,liveesel,r,ename,,3 !选择beam3单元类型etable,fx_i,smisc,1etable,fx_j,smisc,7etable,fy_i,smisc,2etable,fy_j,smisc,8etable,mz_i,smisc,6etable,mz_j,smisc,12plls,mz_i,mz_j,-0.2 !显示弯矩（0.2为显示比例) plls,fx_i,fx_j,0.2 !显示轴力plls,fy_i,fy_j,0.2 !显示剪力!查看锚杆轴力esel,s,liveesel,r,ename,,1etable,forx,smisc,1plls,forx,forx,0.05!查看位移场file,ch-three-2D,rstset,3,lastlcdef,1,1 !定义第1个载荷步为工况1 lcoper,sub,1 !当前载荷步的结果减去工况1的结果esel,s,liveesel,u,mat,,3plnsol,u,sum !显示合位移云图plnsol,u,xplnsol,u,y!载荷步4结果resume,s4,db/pbc,all,,0esel,s,liveesel,u,mat,,13plnsol,s,eqvplnsol,s,1plnsol,s,3plnsol,eppl,eqv,0,1!中墙等效应力场esel,s,mat,,13plnsol,s,eqvplnsol,s,1plnsol,s,3plnsol,u,y!左及右隧道衬砌内力esel,s,liveesel,r,ename,,3 !选择beam3单元类型etable,fx_i,smisc,1etable,fx_j,smisc,7etable,fy_i,smisc,2etable,fy_j,smisc,8etable,mz_i,smisc,6etable,mz_j,smisc,12plls,mz_i,mz_j,-0.2 !显示弯矩（0.2为显示比例) plls,fx_i,fx_j,0.2 !显示轴力plls,fy_i,fy_j,0.4 !显示剪力!查看锚杆轴力esel,allesel,r,ename,,1etable,forx,smisc,1plls,forx,forx,0.05!查看位移场file,ch-three-2D,rstset,4,lastlcdef,1,1 !定义第1个载荷步为工况1 lcoper,sub,1 !当前载荷步的结果减去工况1的结果esel,s,liveesel,u,mat,,3plnsol,u,sum !显示合位移云图plnsol,u,xplnsol,u,y。

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为半径。

隧道开挖ansys模拟分析

隧道台阶法开挖的有限元模拟分析1.力学模型的建立岩体的性质是十分复杂的，在地下岩体的力学分析中，要全面考虑岩体的所有性质几乎是不可能的。

建立岩体力学模型，是将一些影响岩石性质的次要因素略去，抓住问题的主要矛盾，即着眼于岩体的最主要的性质。

在模型中，简化的岩体性质有强度、变形、还有岩体的连续性、各项同性及均匀性等。

考虑岩石的性质和变形特性，以及外界因素的影响，采用的模型有弹性、塑性、弹塑性、粘弹性、粘弹塑性等。

根据对隧道的现场调查及试验结果分析，围岩具有明显的弹塑性性质。

因此，根据隧道的实际情况，考虑岩体的弹塑性性质，在符合真实施工工序和支护措施的基础上，在数值模拟过程中将计算模型简化成弹塑性平面应变问题，采用Drucker—Prager屈服准则来模拟围岩的非线性并且不考虑其体积膨胀，混凝土材料为线弹性且不计其非线性变形。

对地下工程开挖进行分析，一般有两种计算模型：（1）“先开洞，后加载”在加入初始地应力场前，首先将开挖掉的单元从整体刚度矩阵中删除，然后对剩余的单元加入初始地应力场进行有限元计算。

（2）“先加载，后开洞”这种方法是首先在整个计算区域内作用地应力场，然后在开挖边界上施加反转力，经过有限元计算得到所需要的应力、位移等物理量。

两种方法对线弹性分析而言，所得到的应力场是相同的，而位移场是不同的，模型（2）（即：“先加载，后开洞”）更接近实际情况。

在实际地下工程开挖中部分岩体已进入塑性状态，必须用弹塑性有限元进行计算分析，而塑性变形与加载的路径有关，所以模拟计算必须按真实的施工过程进行，即在对地下工程开挖进行弹塑性数值模拟过程中，必须遵循“先加载，后开洞”的原则。

在有限元法中，求解非线性问题最常采用的方法是常刚度初应力法。

对于弹塑性问题，由于塑性变形不可恢复，应力和应变不再是一一对应的关系，即应力状态与加载路径有关，因此应该用增量法求解。

弹塑性应力增量与应变增量之间的关系可近似地表示为}{}]{[}]){[]([}{][}{0σεεεσd d D d D D d D d p ep +=-== （1）式中，][D —弹性矩阵，][p D —塑性矩阵。

ansys命令流中文说明

ansys命令流中文说明KB、KE: 待划分线的定向关键点起始、终止号SECNUM: 截面类型号u SECPLOT,SECID,MESHKEY 画梁截面的几何形状及网格划分SECID:由SECTYPE命令分配的截面编号MESHKEY：0：不显示网格划分1：显示网格划分u /ESHAPE, SCALE 按看似固体化分的形式显示线、面单元SCALE: 0:简单显示线、面单元1：使用实常数显示单元形状u esurf, xnode, tlab, shape 在已存在的选中单元的自由表面覆盖产生单元xnode: 仅为产生surf151 或surf152单元时使用tlab: 仅用来生成接触元或目标元top 产生单元且法线方向与所覆盖的单元相同，仅对梁或壳有效，对实体单元无效Bottom产生单元且法线方向与所覆盖的单元相反，仅对梁或壳有效，对实体单元无效Reverse 将已产生单元反向Shape: 空与所覆盖单元形状相同Tri 产生三角形表面的目标元注意：选中的单元是由所选节点决定的，而不是选单元，如同将压力加在节点上而不是单元上u Nummrg,label,toler, Gtoler,action,switch 合并相同位置的itemlabel: 要合并的项目node: 节点， Elem,单元，kp: 关键点（也合并线，面及点）mat: 材料，type: 单元类型，Real: 实常数cp：耦合项，CE：约束项，CE: 约束方程，All：所有项toler: 公差Gtoler：实体公差Action: sele 仅选择不合并空合并switch: 较低号还是较高号被保留（low, high）注意：可以先选择一部分项目，再执行合并。

如果多次发生合并命令，一定要先合并节点，再合并关键点。

合并节点后，实体荷载不能转化到单元，此时可合并关键点解决问题。

u Lsel, type, item, comp, vmin, vmax, vinc, kswp 选择线type: s 从全部线中选一组线r 从当前选中线中选一组线a 再选一部线附加给当前选中组aunoneu(unselect)inve: 反向选择item: line 线号loc 坐标length 线长comp: x,y,zkswp: 0 只选线1 选择线及相关关键点、节点和单元u Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点为下一步做准备Type: S: 选择一组新节点（缺省）R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0” 使用正负号“1”仅用绝对值u NSLL,type, nkey 选择与所选线相联系的节点u nsla, type, nkey: 选择与选中面相关的节点type：s 选一套新节点r 从已选节点中再选a 附加一部分节点到已选节点u 从已选节点中去除一部分nkey: 0 仅选面内的节点1 选所有和面相联系的节点（如面内线，关键点处的节点）u esel, type, item, comp, vmin, vmax, vinc, kabs 选择一组单元Type: S: 选择一组单元（缺省）R: 在当前组中再选一部分作为一组A: 为当前组附加单元U: 在当前组中不选一部分单元All: 选所有单元None: 全不选Inve: 反向选择当前组（？）Stat: 显示当前选择状态Item： Elem: 单元号Type: 单元类型号Mat: 材料号Real: 实常数号Esys: 单元坐标系号u ALLSEL, LABT, ENTITY 选中所有项目LABT: ALL: 选所有项目及其低级项目BELOW: 选指定项目的直接下属及更低级项目ENTITY: ALL: 所有项目（缺省）VOLU:体高级AREA:面LINE :线KP:关键点ELEM:单元NODE:节点低级u Tshap,shape 定义接触目标面为2D、3D的简单图形Shape: line:直线Arc:顺时针弧Tria:3点三角形Quad:4点四边形………….2.6 根据需要耦合某些节点自由度u cp, nset, lab,,node1,node2,……node17nset: 耦合组编号lab: ux,uy,uz,rotx,roty,rotznode1-node17: 待耦合的节点号。

史上最全的ANSYS命令流大全

《史上最全的ANSYS命令流查询与解释》【1】*************************************************************************************对ansys主要命令的解释1，/PREP7 ! 加载前处理模块2，/CLEAR,NOSTART ! 清除已有的数据, 不读入启动文件的设置(不加载初始化文件)初始化文件是用于记录用户和系统选项设置的文本文件/CLEAR, START !清除系统中的所有数据,读入启动文件的设置/FILENAME, 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 ! 申明分析类型是静力分析(STA TIC或者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隧道命令流一本书上的隧道命令流，我试了几次总有问题，发出来请大家帮忙看下哪里错了。

某明挖隧道由围岩、支护、锚杆几个部分组成。

fini/cle/filename,exm4/title,the analysis of the arch tunnel*afun,deg!隧道分析范围*set,x_zs,-50*set,y_zs,30*set,x_yx,50*set,y_yx,-45!加锚区设定*set,num,30*set,th_s,0*set,th_f,180-th_s*set,th,(th_f-th_s)/num!隧道断面*set,yx_1_x,0*set,yx_1_y,0*set,yx_2_x,0*set,yx_2_y,11.37*set,yx_3_x,3.34*set,yx_3_y,0*set,yx_4_x,-3.34*set,yx_4_y,0*set,r_1,5.8*set,r_2,14.4*set,r_3,2.67 *set,r_4,2.67 *set,r_jg_1,8.8 *set,r_jg_2,17.4 *set,r_jg_3,5.67 *set,r_jg_4,5.67 !**********set,t1,0.25*set,r1,25e-3 e_r=1.5e9u_r=0.4dens_r=1826.5 c_r=0.16e6fi_r=28!*******u_c=0.2dens_c=2449 !***********8 e_a=170e9u_a=0.3dens_a=7959 !********/prep7et,1,3 keyopt,1,6,1 et,2,42 keyopt,2,3,2 et,3,1!衬砌mp,ex,1,e_cmp,prxy,1,u_cmp,dens,1,dens_c!围岩mp,ex,2,e_rmp,prxy,2,u_rmp,dens,2,dens_r tb,dp,2tbdata,1,c_r,fi_r!挖去的土体mp,ex,3,e_rmp,prxy,3,u_rmp,dens,3,dens_r tb,dp,3tbdata,1,c_r,fi_r!锚杆mp,ex,4,e_amp,prxy,4,u_amp,dens,4,dens_a!r,1,t1,t1/12,t1r,2,3.14*(r1/2)**2 !********k,1,yx_1_x,yx_1_y,0 k,2,yx_2_x,yx_2_y,0 k,3,yx_3_x,yx_3_y,0 k,4,yx_4_x,yx_4_y,0k,101,yx_1_x,yx_1_y,10 k,102,yx_2_x,yx_2_y,10 k,103,yx_3_x,yx_3_y,10 k,104,yx_4_x,yx_4_y,10 circle,1,r_1,101 circle,2,r_2,102circle,3,r_3,103circle,4,r_4,104lcsl,alllsel,s,,,23,26,1lsel,a,,,29,42,1lsel,a,,,5,6,1lsel,a,,,10,13,3lsel,a,,,44,45,1lsel,invecm,zf,linecmsel,s,zf,linelsel,inveldel,all,,,lallselcircle,1,r_jg_1,101circle,2,r_jg_2,102circle,3,r_jg_3,103circle,4,r_jg_4,104 cmsel,s,zf,line lsel,invelcsl,alllsel,s,,,49,58,1lsel,a,,,31,34,1lsel,a,,,37,42,1lsel,a,,,60,65,5lsel,a,,,61,64,3lsel,a,,,44,45,1lsel,a,,,5,6,1ldel,all,,,1allsel!隧道分析范围rectng,x_zs,x_yx,y_zs,y_yx asbl,all,all !******生成锚杆asel,s,,,4wprota,,-90 wprota,,,-th_s*do,i,1,num,1 wprota,,,-th asbw,allwpcsysasel,s,loc,y,th_s,th_f asel,inveaadd,allallsel!********定义两集合asel,s,,,2lsla,slsel,r,loc,y,th_s,th_f cm,zf_mq,line asel,s,,,2lsla,slsel,r,loc,y,th_f,th_s+360 cm,zf_wmq,line allsel!利用工作平面切割面csyswpcsyswpof,-1.5*r_jg_1，wprota,,,-90 asbw,allwpcsyswpof,1.5*r_jg_1，wprota,,,-90 asbw,allwpcsyswpof,,1.2*r_jg_1，wprota,,90asbw,allwpcsyswpof,,-0.6*r_jg_1，wprota,,90 asbw,allwpcsys!划分支护单元cmsel,s,zf_mq,linecmsel,a,zf_wmq,linelmesh,allallsel!/efacet,1!***锚杆type,3mat,4real,2wpcsyscsys,1lsel,s,loc,y,th_s*do,i,1,num,1 lsel,a,loc,y,th_s+th*i, *enddolsel,u,,,49,59,10 lsel,u,,,28,72,44 lmesh,allallsel!******挖去土体type,2mat,3mshkey,0 mshape,0esize,1amesh,2lsel,s,,,68,70,1 lsel,a,,,72lsel,a,,,25,65,40 lsel,a,,,28,59,31lesize,all,1allselesize,3!划分围岩的有限元网格type,2mat,2asel,s,loc,x,1.5*r_jg_1,100*r_jg_1allselasel,s,loc,x,1.5*r_jg_1,100*r_jg_1 asel,inve asel,u,,,2mshape,0amesh,allallsel!边界条件施加csys,0nsel,s,loc,x,x_zsnsel,a,loc,x,x_yxd,all,uxallselnsel,s,loc,y,y_yxd,all,uy/soluantype,staticpred,onlnsrch,onnlgeom,onnropt,fulloutres,all,allnsub,7,10!step1time,1acel,,,9.8esel,s,type,,1,3,2ekill,allesel,allesel,s,livensle,snsel,inved,all,allallselsolvesave,step1,db!step2time,2esel,s,mat,,3ekill,allesel,s,type,,1,3,2nsle,sddele,all,allesel,allesel,s,livensle,snsel,inved,all,all,0nsel,allesel,allsolvesave,step2,db!*********************** /post1 resume,'step1','db'file,'exm4','rst'set,1,lastesel,s,liveplnsol,u,sum,0,1.0plnsol,u,xplnsol,u,yplnsol,s,eqvresume,'step2','db'file,'exm4','rst'set,2,last!************esel,s,type,,1etable,fx_i,smisc,1etable,fx_j,smisc,7etable,fy_i,smisc,2etable,fx_j,smisc,8etable,mz_i,smisc,6etable,mz_j,smisc,12plls,mz_i,mz_j,-0.5plls,fx_i,fx_j,0.2plls,fy_i,fy_j,0.5esel,s,type,,3etable,forx,smisc,1plls,forx,forx,0.05etable,ybx,lepel,1plls,ybx,ybx,0.05resume,'step2','db' file,'exm4','rst' set,2,lastlcdef,1,1lcoper,sub,1 esel,s,live plnsol,u,sum plnsol,u,x plnsol,u,y。

(完整版)ANSYS命令流总结(全)

ANSYS结构解析单元功能与特征/POST1/可以构成一一些命令，一般是一种整体命令（ session),三十也有特别，比方是办理 ! 是说明说明符号，，与其余软件的说明是相同的， 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, V ALUELab 为参数名称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杆单元单元名称简称节点数节点自由度特征备注LINK12D杆2Ux,Uy EPCSDGB常用杆元LINK83D杆Ux,Uy,Uz EPCSDGBLINK103D仅受拉EDGB模拟缆索的废弛及或仅受压杆缝隙LINK113D线性调理EGB模拟液压缸和大转器动LINK1803D有限应变杆EPCDFGB另可考虑粘弹塑性E- 弹性 (Elasticity)，P-塑性(Plasticity)，C-蠕变(Creep)，S-膨胀(Swelling)，D-大变形或大挠度deflection)， F- 大应变 (Large strain)或有限应变(Finite strain)，B-单元存亡(Birth and dead),G-化 (Stress stiffness)或几何刚度(Geometric stiffening)，A-自适应降落(Adaptive descent)等。

ANSYS工程结构数值分析---命令流实例

!------------------------------------------------------!EX8.26 钢筋混凝土简支梁数值分析!分离式模型,关闭压碎,keyopt(1)=0,keyopt(7)=1 !力加载,位移收敛准则，误差1.5%,1/4模型分析!--------------------------------------------finish/clear/config,nres,2000/prep7!1.定义单元与材料性质--------------------et,1,solid65,,,,,,,1 !K1=0,k7=1et,2,link8mp,ex,1,13585mp,prxy,1,0.2fc=14.3ft=1.43tb,concr,1tbdata,,0.5,0.95,ft,-1tb,miso,1,,11tbpt,,0.0002,fc*0.19tbpt,,0.0004,fc*0.36tbpt,,0.0006,fc*0.51tbpt,,0.0008,fc*0.64tbpt,,0.001,fc*0.75tbpt,,0.0012,fc*0.84tbpt,,0.0014,fc*0.91tbpt,,0.0016,fc*0.96tbpt,,0.0018,fc*0.99tbpt,,0.002,fctbpt,,0.0033,fc*0.85mp,ex,2,2.0e5mp,prxy,2,0.3tb,biso,2tbdata,,300,0pi=acos(-1)r,1,0.25*pi*22*22r,2,0.25*pi*22*22/2r,3,0.25*pi*10*10r,4,0.25*pi*10*10/2!2.创建几何模型blc4,,,150/2,300,2000/2*do,i,1,9wpoff,,,100vsbw,all*enddo wpcsys,-1 wpoff,,,50 vsbw,all wpcsys,-1 wprota,,-90 wpoff,,,30 vsbw,all wpoff,,,240 vsbw,all wpcsys,-1 wpoff,30 wprota,,,90 vsbw,all wpcsys,-1!3.划分钢筋网格elemsiz=50 lsel,s,loc,x,30 lsel,r,loc,y,30 cm,zj,linelatt,2,1,2 lesize,all,elemsiz lsel,s,loc,x,75 lsel,r,loc,y,30 cm,zjb,line latt,2,2,2 lesize,all,elemsiz lsel,s,loc,x,30 lsel,r,loc,y,270 cm,jlj,linelatt,2,3,2 lesize,all,elemsiz lsel,s,tan1,z lsel,r,loc,y,30,270 lsel,r,loc,x,30,70 lsel,u,loc,z,50 cm,gj,linelatt,2,3,2 lesize,all,elemsiz lsel,s,loc,z,0 lsel,r,loc,y,30,270 lsel,r,loc,x,30,70 cm,gjb,line latt,2,4,2lesize,all,elemsizlsel,allcmsel,s,zjcmsel,a,zjbcmsel,a,jljcmsel,a,gjcmsel,a,gjbcm,gj,linelmesh,alllsel,all!4.划分混凝土网格vatt,1,,1mshkey,1esize,elemsizvmesh,allallsel,all!5.施加荷载和约束lsel,s,loc,y,0lsel,r,loc,z,900dl,all,,uyasel,s,loc,z,0da,all,symmasel,s,loc,x,75da,all,symm!-----------------p0=180000q0=p0/150/100asel,s,loc,z,0,50asel,r,loc,y,300sfa,all,1,pres,q0allsel,all!6.求解控制设置/soluantype,0nsubst,80outres,all,allautos,onneqit,50cnvtol,u,,0.015solve!7.进入POST1查看结果/post1set,lastesel,s,type,,2etable,saxl,ls,1plls,saxl,saxlesel,s,type,,1/device,vector,onplcrackeall!8.进入时程后处理/post26nsol,2,205,u,yprod,3,1,,,,,,p0/1000prod,4,2,,,,,,-1xvar,4plvar,3!EX2.6D 60kg/m钢轨实体建模FINISH/CLEAR/PREP7!1.定义参数H1=30.5H2=48.5H3=48.5B1=45.75B2=29.25B3=16.5B4=36.5B5=10XD1=1/3XD2=1/9XD3=1/3XD4=1/20R0=400R1=2R2=4R3=40R4=20R5=25R7=13R8=80R9=300R10=2FAI=43H=H1+2*H2+H3!2.创建轨底主要组成部分的线K,1K,2,B1+B2K,3,B1+B2,H1-B1*XD1-B2*XD2K,4,B1,H1-B1*XD1K,5,0,H1*DO,I,1,4L,I,I+1*ENDDO!3.计算确定轨头下三个圆弧的位置CTA=A TAN(XD3)BTA=(90-FAI)*ACOS(-1)/180-CTAA=1/TAN(CTA)C1=B3/2+R0-R6*SIN(CTA)-(R5-R6)*COS(BTA)C2=R5*SIN(BTA)+R6*COS(CTA)-R6*SIN(BTA)-H2C=(C1-C2/TAN(CTA))/(R0-R5)C3=A*C-SQRT(A*A-C*C+1)C3=C3/(A*A+1)REFA=ASIN(C3)C1=R6*SIN(CTA)+(R5-R6)*COS(BTA)+(R0-R5)*COS(REFA) !4.创建轨腰和轨头结构线ROUX=B3/2+R0-C1ROUY=ROUX*TAN(CTA)ANGI=ASIN((H1+H2)/R0)K,6,B3/2+R0-R0*COS(ANGI)K,7,B3/2+R0-R0*COS(REFA),H1+H2+R0*SIN(REFA) LARC,6,7,2,R0DX1=R6*(COS(BTA)-SIN(CTA))DY1=R6*(COS(CTA)-SIN(BTA))K,8,ROUX-DX1,H1+2*H2+ROUY-DY1K,9,ROUX,H1+2*H2+ROUYLARC,7,8,2,R5LARC,8,9,2,R6K,10,B4,H1+2*H2+B4*XD3K,11,B4-(H3-B4*XD3)*XD4,HL,9,10L,10,11ANGI=ASIN(B5/R9)K,12,0,HK,13,B5,H-(R9-R9*COS(ANGI)) LARC,12,13,1,R9DX1=B4-B5+R8*SIN(ANGI)ANGI1=ASIN(DX1/R8)Y1=R8*COS(ANGI)-R8*COS(ANGI1) K,14,B4,KY(13)-Y1LARC,13,14,1,R8!5.倒角或弧线连接LFILLT,1,2,R1LFILLT,2,3,R2LFILLT,3,4,R3LPTN,4,5LFILLT,17,18,R4LFILLT,8,9,R10LPTN,9,11LFILLT,20,22,R7LDELE,15,16,1,1LDELE,19,21,2,1LSYMM,X,ALLNUMMRG,ALLNUMCMP,ALLAL,ALLASUMVOFFST,1,300!EX2.9B 斜向圆台建模finish/clear/prep7r0=3r1=2h=4cyl4,,,r0wpoff,r1,,hcyl4,,,r1*do,i,1,4askin,i,i+4*enddova,all!斜向延伸创建圆台vgen,1,1,,,,3*r0,,,,1wpcsyscyl4,,,r0vext,7,,,r1,,h,r1/r0,r1/r0!EX2.9B 已知函数方程时的建模!玫瑰线1==================================== !极坐标方程为ρ=asin(4Φ/3)finish/clear/prep7a=20csys,1*do,i,0,360*3fei=i*acos(-1)/180r0=a*sin(fei*4/3)k,,r0,i*enddo*do,i,1,360*3l,i,i+1*enddo!玫瑰线2=================================== !极坐标方程为ρ=asin(4Φ)finish/clear/prep7a=20csys,1*do,i,0,360fei=i*acos(-1)/180r0=a*sin(fei*4)k,,r0,i*enddo*do,i,1,360l,i,i+1*enddo!长辐圆内旋轮线============================ !参数方程x=(a-b)cost+λcos(a-b)/bt!参数方程y=(a-b)sint-λsin(a-b)/btfinish/clear/prep7a=50b=5lmda=15*afun,deg*do,i,0,360x1=(a-b)*cos(i)+lmda*cos((a-b)/b*i)y1=(a-b)*sin(i)-lmda*sin((a-b)/b*i)k,,x1,y1*enddo*do,i,1,360l,i,i+1*enddo!椭圆抛物面================================ !方程为z=x*x/a/a+y*y/b/b!使用蒙皮方法finish/clear/prep7a=20b=10n=20*do,i,1,niz=i/10lsel,nonewpoff,,,izcyl4,,,1adele,alllsscale,all,,,a*sqrt(iz),b*sqrt(iz),,,,1*enddoallsel*do,i,1,4*(n-1),4askin,i,i+4askin,i+1,i+5askin,i+2,i+6askin,i+3,i+7*enddo!EX2.9C 旋转图素建模!类花瓶建模--------------------finish/clear/prep7K,1K,2,100K,3,300,400K,4,150,550k,5,130,880K,6,300,1150k,7,400,1100k,8,500,1150bsplin,allk,1001,,1000arotat,all,,,,,,1,1001,360,5!类轮建模----------------------finish/clear/prep7!先创建类轮结构的一个断面，然后旋转之blc4,,,4,6cyl4,2,6,1.5asba,1,2k,100,-3k,101,-3,10vrotat,all,,,,,,100,101!EX3.1 自定义箱形截面finish/clear/prep7k,1k,2,2k,3,2,2.2k,4,3,2.3k,5,3,2.5k,6,0,2.5a,1,2,3,4,5,6blc4,,0.2,1.7,2asba,1,2wpoff,1.7wprota,,,90asbw,allwpoff,,,0.3asbw,allwpoff,,0.2wprota,,90asbw,allwpoff,,,-2asbw,allarsym,x,allwpcsysaglue,allet,1,plane82lesize,all,,,1amesh,allsecwrite,myboxfinish/clear/prep7et,1,beam189mp,ex,1,3.0e10mp,prxy,1,0.167sectype,1,beam,mesh secread,mybox,,,meshk,1k,2,10k,100,5,5l,1,2lesize,all,,,20latt,1,,1,,100,,1lmesh,all/eshape,1Eplot!EX3.2 自定义多种材料截面finish/clear/prep7Ro=1.5Ri=1.0csys,1cyl4,,,ricyl4,,,roaptn,allwprota,,90asbwa,allwprota,,,90asbw,allwpcsyset,1,plane82mymat1=4mymat2=7mp,ex,mymat1,1.0mp,ex,mymat2,2.0asel,s,loc,x,0,riaatt,mymat1,,1asel,s,loc,x,ri,roaatt,mymat2,,1allselesize,0.25mshape,0,2dmshkey,1amesh,allsecwrite,mycsolid,sectfinish/clear/prep7et,1,beam189mym1=4mym2=7mp,ex,mym1,3.0e10mp,prxy,mym1,0.167 mp,ex,mym2,2.1e11mp,prxy,mym2,0.3 sectype,1,beam,mesh secread,mycsolid,sect,,mesh k,1k,2,,,10l,1,2lesize,all,,,20latt,,,1,,,,1lmesh,all/eshape,1/pnum,mat,1eplot!EX3.3A 单个方位关键点示例FINISH/CLEAR/PREP7ET,1,BEAM189MP,EX,1,2.1E5MP,PRXY,1,0.3SECTYPE,1,BEAM,ISECDATA,100,40,160,10,10,8K,1K,2,,,1000L,1,2K,100,,500,500LATT,1,,1,,100,,1LGEN,4,1,,,500LESIZE,ALL,,,10LMESH,ALL/ESHAPE,1EPLOT!EX3.3B 同时使用KB和KE方位关键点示例FINISH/CLEAR/PREP7ET,1,BEAM189MP,EX,1,2.1E5MP,PRXY,1,0.3SECTYPE,1,BEAM,ISECDATA,100,40,160,10,10,8L0=1000DL=500DXC=400K,1K,2,,,L0L,1,2K,100,,DLK,200,DXC,-DLK,300,2*DXC,DLK,301,2*DXC+DLK,400K,500,8*DXCLGEN,5,1,,,DXCLATT,1,,1,,100,,1LSEL,S,,,2LATT,1,,1,,200,,1LSEL,S,,,3LATT,1,,1,,300,301,1LSEL,S,,,4LATT,1,,1,,400,,1LSEL,S,,,5LATT,1,,1,,500,,1LSEL,ALLLESIZE,ALL,,,50LMESH,ALL/ESHAPE,1EPLOT!EX3.3C 双方位关键点翻转示例FINISH/CLEAR/PREP7ET,1,BEAM189MP,EX,1,2.1E5MP,PRXY,1,0.3SECTYPE,1,BEAM,I SECDATA,100,40,160,10,10,8L0=1000DL=500DXC=400K,1K,2,,,L0K,3,DXCK,4,DXC,,L0/2K,5,DXC,,L0L,1,2L,3,4L,4,5K,100,,DLk,101,,-DLK,200,DXC,DLk,201,2*DXCK,202,DXC,-DLLATT,1,,1,,100,101,1lesize,all,,,100LSEL,S,,,2LATT,1,,1,,200,201,1lesize,all,,,50LSEL,S,,,3LATT,1,,1,,201,202,1lesize,all,,,50LSEL,ALLLMESH,ALL/ESHAPE,1EPLOT!EX3.4A 两种单元形状和两种网格划分比较finish/clear/prep7et,1,plane82k,1k,2,8k,3,7,6k,4,1,6a,1,2,3,4esize,1mshape,0mshkey,0!mshape,0!mshkey,1!mshape,1!mshkey,0!mshape,1!mshkey,1amesh,all!EX3.4B 中间节点位置控制网格划分比较finish/clear/prep7cyl4,,,4,,8,60lesize,all,,,2mshape,0mshkey,1mshmid,0!mshmid,1!mshmid,2amesh,all!EX3.5 线上单元尺寸设置示例!下边密上边稀finish/clear/prep7et,1,plane82blc4,,,10,10lsel,s,tan1,ylesize,all,,,10lsel,s,loc,x,0lesize,all,,,9,1/8lsel,s,loc,x,10lesize,all,,,9,8lsel,allmshape,0mshkey,1amesh,all!中间密外边稀finish/clear/prep7et,1,plane82blc4,,,10,10lsel,s,tan1,ylesize,all,,,10,-1/5lsel,s,tan1,xlesize,all,,,9,-1/8lsel,allmshape,0mshkey,1amesh,all!EX3.6 合并线和连接线以进行映射网格划分finish/clear/prep7et,1,plane82k,1,5k,2,10k,3,11,6k,4,6,15k,5,-1,8k,6,,4l,1,2l,2,3l,3,4larc,4,5,3,10l,5,6l,6,1al,allesize,3mshape,0mshkey,1lccat,1,2lcomb,4,5amesh,all!EX3.7 合并线或连接线的网格划分设置finish/clear/prep7et,1,plane82k,1,5k,2,10k,3,11,6k,4,6,15k,5,-1,8k,6,,4l,1,2l,2,3l,3,4larc,4,5,3,10l,5,6l,6,1al,allesize,,10lesize,6,,,8lesize,1,,,4lesize,2,,,3lesize,4,,,5lesize,5,,,2mshape,0mshkey,1lcomb,1,2!lccat,1,2!如采用lcatt,1,2则采用7个划分数lcomb,4,5amesh,all!EX3.7B 简化映射网格划分finish/clear/prep7et,1,plane82k,1,5k,2,10k,3,11,6k,4,6,15k,5,-1,8k,6,,4l,1,2l,2,3l,3,4larc,4,5,3,10l,5,6l,6,1al,allesize,,6mshape,0amap,1,2,5,3,4!EX3.8A 过渡四边形映射网格finish/clear/prep7k,1k,2,10,-1k,3,8,6k,4,1,3a,1,2,3,4lesize,1,,,8lesize,3,,,3lesize,4,,,7lesize,2,,,2mshape,0,2dmshkey,1amesh,all!EX3.8B 过渡四边形映射网格finish/clear/prep7et,1,plane42k,1k,2,10,-1k,3,8,6k,4,1,3a,1,2,3,4lesize,1,,,11lesize,3,,,3lesize,4,,,2lesize,2,,,2mshape,0,2dmshkey,1amesh,all!EX3.9A 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lesize,all,,,3lesize,5,,,15allselmshape,0,3dvmesh,all!EX3.9B 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lesize,all,,,4lesize,7,,,12mshape,0,3dmshkey,1vmesh,all!EX3.9C 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lsel,s,,,1,8,7lesize,all,,,6lsel,s,,,2,7,5lesize,all,,,3lsel,s,,,4,5lesize,all,,,7lsel,s,,,3,6,3lesize,all,,,2lsel,s,,,9,12,1lesize,all,,,5allselmshape,0,3dmshkey,1vmesh,all!EX3.9D 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lesize,2,,,2lesize,4,,,2lesize,5,,,2lesize,1,,,3lesize,3,,,3lesize,7,,,6lesize,6,,,7lesize,11,,,7lesize,9,,,3lesize,10,,,3lesize,12,,,3allselmshape,0,3dmshkey,1vmesh,all!EX3.10 带两孔的长方体的扫掠网格划分finish/clear/prep7a=10r=2et,1,mesh200,6et,2,solid45blc4,,,2*a,a,acyl4,a/2,a/2,r,,,,awprota,,90cyl4,1.5*a,a/2,r,,,,-avsbv,1,2vsbv,4,3wprota,,,90wpoff,,,a/2vsbw,allwpoff,,,a/2vsbw,allwpoff,,,a/2vsbw,allwpcsyswpoff,,a/2,a/2vsbw,allwprota,,90vsbw,allwpcsysesize,1amap,105,15,16,26,63amap,107,16,13,60,26amap,108,13,14,28,60amap,103,15,14,28,63lesize,94,,,4lesize,79,,,5vsel,s,loc,x,0,avsweep,allasel,s,loc,y,aasel,r,loc,x,a,2*alccat,2,45lccat,57,71lccat,68,78lccat,65,72mshape,0,2dmshkey,1amesh,allvsel,s,loc,x,a,2*avsweep,allallsel/view,1,1,2,3Eplot!EX3.11 拉伸类命令生成体单元网格finish/clear/prep7et,1,82et,2,95blc4,,,4,4blc4,6,,4,4blc4,12,,4,4blc4,18,,4,4esize,1amesh,allesize,,8vrotat,1,,,,,,1,4,90vext,2,,,,,10,0.5,0.5voffst,3,10vdrag,4,,,,,,35!EX3.12 圆的网格划分finish/prep7et,1,plane82r0=10cyl4,,,r0cyl4,3*r0,,,,r0 wprota,,90asbw,allwprota,,,90asbw,allwpoff,,,3*r0asbw,allwpcsys,-1asel,s,loc,x,-r0,r0lsla,slesize,all,,,8mshape,0,2dmshkey,1amesh,allallselasel,s,loc,x,2*r0,4*r0 lesize,all,,,8lsel,r,length,,r0lesize,all,,,8,0.1,1 amesh,allallsel!EX3.13 圆环的网格划分finish/clear/prep7et,1,plane82r0=10cyl4,,,r0/3,,r0,90cyl4,2*r0,,r0/10,,r0,90 asel,s,loc,x,-r0,r0lsla,slesize,all,,,8lsel,r,length,,r0*2/3 lesize,all,,,3,,1 mshape,0,2dmshkey,1amesh,allasel,s,loc,x,2*r0,4*r0lesize,5,,,12lesize,7,,,6lsel,s,length,,r0*9/10lesize,all,,,7amesh,all!EX3.14 圆柱面和圆柱体的网格划分!圆柱面finish/clear/prep7r0=10h0=50et,1,shell63cyl4,,,r0adele,1cm,l1cm,linek,50k,51,,,h0l,50,51adrag,l1cm,,,,,,5lsel,s,loc,z,0lesize,all,,,6lsel,s,length,,h0lesize,all,,,8mshape,0,2dmshkey,1amesh,all!圆柱体finish/clear/prep7r0=10h0=50et,1,solid95cyl4,,,r0,,,,h0wprota,,90vsbw,allwprota,,,90vsbw,allmshape,0,3dmshkey,1lsel,s,loc,z,0lsel,s,length,,h0lesize,all,,,8vmesh,all!EX3.15 圆锥的网格划分!扫掠网格划分finish/clear/prep7cone,10,,,15,,90et,1,200,7et,2,95lsel,s,,,5,6lesize,all,,,12,0.5lsel,alllesize,3,,,6mshape,0,2dmshkey,1amesh,3vsweep,1,3,4vsymm,x,allvsymm,y,allvglue,all!六面体映射网格划分finish/clear/prep7cone,10,,,15,,90et,2,95lesize,all,,,8mshape,0,3dmshkey,1vmesh,allvsymm,x,allvsymm,y,allvglue,all!EX3.16 正多边形面的通用网格划分finish/clear/prep7ns=11rpr4,ns,,,10kp0=100+nsk,kp0*if,mod(ns,2),eq,0,then*do,i,1,ns/2l,kp0,2*i-1*enddo*else*do,i,1,nsl,kp0,i*enddo*endiflsel,s,,,ns+1,2*nscm,l1cm,lineallselasbl,1,l1cm*if,mod(ns,2),eq,0,then cmsel,s,l1cmlesize,all,,,8lsel,invelesize,all,,,4*elselesize,all,,,8*endifallselmshape,0mshkey,1amesh,all!EX3.17 球体及球面网格划分finish/clear/prep7r0=10sphere,,r0,,90vsbw,allvdele,2,,,1esize,2et,1,solid95mshape,0,3dmshkey,1vsymm,x,all vsymm,y,all vsymm,z,allvglue,allfinish/clear/prep7r0=10sphere,,r0,,90 vsbw,allvdele,2,,,1vdele,allasel,s,loc,x,0asel,a,loc,y,0asel,a,loc,z,0adele,all,,,1allselesize,2et,1,shell63 mshape,0,2d mshkey,1amesh,allarsym,x,allarsym,y,allarsym,z,all nummrg,all!EX3.18 孔板网格划分finish/clear/prep7a0=300b0=800r0=15blc4,,,a0,b0cyl4,a0/4,b0/8,r0 agen,2,2,,,a0/2 agen,2,2,3,1,,b0/8 agen,2,2,5,1,,b0*5/8 asel,s,,,2,9,1cm,a2cm,areaallselwprota,,-90*do,i,1,5wpoff,,,b0/16asbw,all*enddowpoff,,,b0*5/16*do,i,1,5wpoff,,,b0/16asbw,all*enddowprota,,,90*do,i,1,3wpoff,,,a0/4asbw,all*enddowpcsys,-1numcmp,alllsel,s,radius,,r0lesize,all,,,8lsel,invelesize,all,,,4lsel,allet,1,82mshape,0,2dmshkey,1asel,u,loc,y,b0/16,b0*5/16 asel,u,loc,y,b0*11/16,b0*15/16 lsla,slsel,r,tan1,xlesize,all,50,,,,1amesh,allallselamap,21,30,31,54,62 amap,22,31,32,54,64 amap,32,29,30,62,76 amap,33,29,32,64,76amap,42,34,35,76,82 amap,43,35,36,76,84 amap,41,33,34,53,82 amap,44,33,36,53,84amap,19,22,23,49,64 amap,20,23,24,49,60 amap,9,21,22,64,74 amap,31,21,24,60,74amap,38,26,27,74,84 amap,39,27,28,74,80 amap,37,25,26,50,84 amap,40,25,28,50,80amap,17,14,15,43,59 amap,18,15,16,43,63 amap,7,13,14,59,73 amap,8,13,16,63,73amap,34,18,19,73,79 amap,30,17,18,79,44 amap,35,19,20,73,83 amap,36,17,20,44,83amap,15,6,7,39,63amap,16,7,8,39,58amap,4,6,5,63,71amap,5,5,8,58,71amap,27,10,11,71,83 amap,28,11,12,71,78 amap,26,9,10,83,40 amap,29,9,12,40,78!EX3.19 角支架的网格划分finish/clear/prep7blc4,,,150,50blc4,100,,50,-50cyl4,,25,25cyl4,125,-50,25aadd,allnumcmp,allcyl4,,25,10cyl4,125,-50,10 asel,s,,,2,3cm,a1cm,area asel,allasba,1,a1cm lcomb,1,6lfillt,1,2,20 asbl,4,6adele,1,,,1lfillt,3,4,20al,18,19,20 aadd,all numcmp,all!---------------------- wprota,,90 wpoff,,,-25 asbw,all wpoff,,,75 asbw,all wprota,,,90 asbw,all wpoff,,,125 asbw,all wpcsys,-1 wpoff,25 wprota,,,90 asbw,all kwpave,18 asbw,all kwpave,3 wprota,,90 asbw,all kwpave,21 wprota,,-45 asbw,8 wpcsys,-1!-------------------- et,1,plane82 mshape,0,2d mshkey,1 esize,6lesize,33,,,6 lesize,37,,,6 lesize,42,,,6amesh,3,5,2amap,6,9,10,4,23amap,7,9,12,1,23amesh,11,13,2amesh,1,9,8amesh,2,4,2amap,16,3,18,26,28amap,12,14,15,5,28amap,15,13,14,6,28lcomb,19,22lcomb,27,46amesh,10,14,4!========================= !EX3.20 大板小孔的网格划分A finish/clear/prep7a0=100blc4,,,a0,a0cyl4,,,a0/100asba,1,2csys,1k,50,a0/20k,51,a0/20,90l,50,51asbl,all,1wprota,,90,45asbw,allet,1,plane82mshape,0,2dmshkey,1lsel,s,loc,x,0,a0/20lesize,all,,,8lsel,alllesize,7,,,12,0.1lesize,8,,,12,0.1lesize,9,,,12,10amesh,allcsys,0arsym,x,allarsym,y,allnummrg,all!============================= !EX3.20 大板小孔的网格划分Bfinish/clear/prep7a0=100blc4,,,a0,a0cyl4,,,a0/100asba,1,2csys,1k,50,a0/20k,51,a0/20,90l,50,51asbl,all,1et,1,plane82mshape,0,2dmshkey,1lesize,5,,,8lesize,1,,,10lesize,4,,,8,5lesize,6,,,8,5amesh,1lesize,7,,,20,0.1lesize,8,,,20,0.1amap,2,50,51,2,4csys,0arsym,x,allarsym,y,allnummrg,all!EX3.21 Z形面网格划分及局部细分finish/clear/prep7blc4,,,15,10blc4,10,6,14,12aadd,allwprota,,-90wpoff,,,6asbw,allwpoff,,,4wprota,,,90wpoff,,,10asbw,allwpoff,,,5asbw,allwpcsys,-1et,1,82esize,2mshape,0,2dmshkey,1amesh,allkrefine,9,10,1,1,,off !不光滑处理（节点位置不变）aclear,allamesh,allkrefine,9,10,1,1,,smooth !光滑处理，节点位置可改变，但既有单元不删除aclear,allamesh,allkrefine,9,10,1,1 !光滑处理，可删除既有单元重细分，节点位置改变!EX3.22 分布小孔环形面的网格划分finish/clear/prep7r1=15r2=25r3=3n=8cyl4,,,r1,,r2,180/ncyl4,0.5*(r1+r2),,r3asba,1,2kl,2larc,6,8,4,0.5*(r1+r2)asbl,all,4lesize,4,,,4lesize,10,,,4lesize,6,,,8lesize,3,,,6lesize,8,,,4lesize,5,,,8lesize,9,,,4lesize,7,,,4et,1,plane82mshape,0,2damap,2,6,7,4,8amap,1,5,6,8,1!lrefine,5,6,1,1,smootharsym,y,allcsys,1agen,n,all,,,,360/nnummrg,all!EX3.23 同心圆轴体的网格划分finish/clear/prep7cylind,50,0,100,150cylind,40,0,200,250cylind,30,0,100,250cylind,20,0,50,250cylind,10,0,0,250vptn,allwprota,,,90vsbw,allwprota,,90vsbw,allet,1,95csys,1lsel,s,length,,50lesize,all,,,5lsel,invelsel,u,radius,,10,50lesize,all,,,4n1=6n2=6n3=8n4=10n5=12lsel,r,loc,x,0,10lesize,all,,,n1,,1lsel,s,radius,,10lesize,all,,,n1lsel,s,radius,,20lsel,s,radius,,30lesize,all,,,n3lsel,s,radius,,40lesize,all,,,n4lsel,s,radius,,50lesize,all,,,n5allselmshape,0,3dmshkey,1vmesh,all!EX3.24A 圆柱与长方体组合模型的网格划分!采用六面体映射网格划分finish/clear/prep7a0=30h1=15h2=25r0=7blc4,,,a0/2,a0/2,h1cyl4,,,r0,,,90,h1+h2vptn,allaccat,4,6et,1,95esize,3mshape,0,3dmshkey,1vmesh,allasel,s,accaadele,alllsel,s,lccaldele,allallselvsymm,x,allvsymm,y,allnummrg,all!EX3.24B 圆柱与长方体组合模型的网格划分!采用过渡六面体映射网格划分finish/cleara0=30h1=15h2=25r0=7blc4,,,a0/2,a0/2,h1cyl4,,,r0,,,90,h1+h2vptn,alllsel,s,radius,,r0lsel,a,length,,r0lesize,all,,,6lsel,s,loc,z,0lsel,a,loc,z,h1lesize,all,,,7asel,s,loc,x,a0/2asel,a,loc,y,a0/2accat,allallselet,1,95esize,3mshape,0,3dmshkey,1vmesh,allasel,s,accaadele,alllsel,s,lccaldele,allallselvsymm,x,allvsymm,y,allnummrg,all!EX3.25-任意圈数圆柱形螺旋弹簧的扫略网格划分FINISH/CLEAR/PREP7!1.定义弹簧参数（同上）D=4C=8DZ=C*DT=DZ/6N=4.7*IF,T,LT,D,THEN*ENDIFTKPD=90!2.创建全部螺旋线CSYS,1TDEG=N*360TDEG1=MOD(TDEG,TKPD) N0=(TDEG-TDEG1)/TKPD+1 *AFUN,DEG*DO,I,1,N0CTA=(I-1)*TKPDZ=T/360*CTAK,I,DZ/2,CTA,Z*ENDDO*IF,TDEG1,LT,1.0E-2,THEN *ELSEN0=N0+1CTA=CTA+TDEG1Z=T/360*CTAK,N0,DZ/2,CTA,Z*ENDIF*DO,I,1,N0-1L,I,I+1*ENDDOCM,L1,LINE!3.在螺旋线端部创建簧丝截面CSYS,0WPOFF,DZ/2 WPROTA,,90CYL4,,,D/2wprota,,90asbw,allwprota,,,90asbw,allcm,a1cm,areaVDRAG,a1cm,,,,,,L1et,1,mesh200,7et,2,solid95cmsel,s,a1cmlsla,slesize,all,,,6amesh,allallselesize,d/2vsweep,all!EX3.26 具坑缺陷圆柱的网格划分finish/clear/prep7l0=50r0=10r1=5r2=1h0=5et,1,solid45cyl4,,,r0,,,,l0wpave,0,r0,l0/2wprota,,90con4,,,r1,r2,h0vsbv,1,2wprota,,90vsbw,allwpoff,,,r1+3vsbw,allwpoff,,,-2*(r1+3)vsbw,allwprota,,,90vsbw,allwpoff,,,r1+1vsbw,allwpoff,,,-2*(r1+1)vsbw,allwpcsys,-1wprota,,90vsbw,allwpoff,,,-r1vsbw,allwpoff,,,2*r1vsbw,allwpcsys,-1numcmp,allksel,s,loc,z,l0/2ksel,r,loc,x,0ksel,u,loc,y,0*get,kp1,kp,0,num,min kp2=kpnext(kp1) ksel,s,loc,z,l0/2+r2*get,kp3,kp,0,num,min allselk,1000,,-r0,l0/2+5l,kp3,1000asel,none*get,l1,line,0,num,max arotat,l1,,,,,,kp1,kp2 cm,a1cm,areaallselvsba,all,a1cm!vsel,s,loc,z,17,33!vsel,r,loc,x,-6,6!aslv,s!aplotaccat,128,61accat,160,138 accat,224,202 accat,183,35accat,31,134accat,151,164 accat,215,231 accat,82,195accat,130,5accat,162,145 accat,229,209 accat,192,89accat,12,123accat,117,155 accat,109,219 accat,14,178mshape,0,3d mshkey,1esize,1vmesh,all!EX3.27长方体开圆柱槽的网格划分finish/clear/prep7blc5,,,20,10,30wpoff,,5,5cyl4,,,8,,,,10wpoff,,,10cyl4,,,4,,,,12vsel,s,,,2,3cm,v1cm,voluallselwpcsys,-1vsbv,1,v1cmvoffst,15,5voffst,16,3voffst,17,15vptn,allwprota,,,90vsbw,allwpcsys,-1wpoff,,,5vsbw,allwpoff,,,10vsbw,allwpoff,,,12vsbw,allet,1,mesh200,7esize,2mshape,0,2dmshkey,1asel,s,loc,z,30lccat,6,74lccat,8,68amesh,allasel,s,loc,z,0lccat,3,73lccat,1,67amesh,allallselet,2,solid95mshape,0,3dmshkey,1vsweep,all!ex4.1累加节点自由度约束值finish/clear/prep7et,1,beam3k,1k,2,10l,1,2esize,1lmesh,alld,1,ux,1e-2,,,,uyd,1,rotzd,2,alldlistnsel,s,d,u,0,0.1dofsel,s,uydscale,2.5dlistdofsel,alldcum,addd,1,ux,2e-2dlistdcum,ignod,2,uy,1.0dlist!ex4.2 对线施加约束并转换finish/clear/prep7et,1,95blc4,,,10,10,10dl,7,,ux,0.1dl,5,,alldl,11,6,symmdl,10,6,asymdl,6,,symmdllistesize,2vmesh,alldtrandlist!ex4.3 施加集中荷载与节点坐标系finish/clear/prep7et,1,beam4k,1k,2,5k,3,10l,1,2l,2,3local,12,0,,,,90nrotat,alldk,1,allfk,2,fy,-1000esize,1lmesh,allnrotat,alllplotfk,3,fy,1000f,6,fx,-1000sbctraneplot!ex4.4a 3d单元sf加载示例finish/clear/prep7et,1,95blc4,,,10,10,20esize,,4vmesh,allasel,s,loc,y,10sf,all,pres,1000asel,s,loc,z,20nsla,s,1sf,all,pres,1000!ex4.4b 2d单元sf加载示例finish/clear/prep7et,1,82blc4,,,100,200blc4,30,60,40,80asba,1,2wprota,,-90wpoff,,,60asbw,allwpoff,,,80asbw,allwprota,,,90wpoff,,,30asbw,allwpoff,,,40asbw,allwpcsys,-1esize,5amesh,all/psf,pres,norm,2sf,all,pres,100sfdele,all,presnsel,s,loc,x,0sf,all,pres,100nsel,s,loc,x,15,20esln,s,1nsel,r,loc,x,15sf,all,pres,110nsel,s,loc,x,40,60nsel,r,loc,y,10,30esln,s,1sf,all,pres,100lsel,s,loc,x,100nsll,s,1esln,snsel,s,loc,x,95sf,all,pres,-100eplot!ex4.5 节点号及其荷载函数finish/clear/prep7et,1,45blc4,,,10,10,20esize,5vmesh,all*dim,mypres,,100*do,i,1,100mypres(i)=i*10.0*enddosffun,pres,mypres(1)nsel,s,loc,y,10sf,all,pres,10sflist*do,i,1,100mypres(i)=i*50.0*enddonsel,s,loc,z,20sf,all,pres,0allselsflist!ex4.6a 2d平面单元plane82 finish/clear/prep7et,1,82blc4,,,10,50esize,2amesh,allsfe,1,4,pres,,100,50sfe,6,4,pres,,200,40nsel,s,loc,x,10esln,ssfe,all,2,pres,,-100/psf,pres,norm,2,0,1eplot!ex4.6b 3d体单元solid95 finish/clear/prep7et,1,95blc4,,,10,10,30esize,5vmesh,all/psf,pres,norm,2,0,1sfe,22,3,pres,,100sfe,22,6,pres,,100,110,120,130 sfelist!ex4.6c 3d壳单元shell63 finish/clear/prep7et,1,63wprota,,90blc4,,,10,10esize,5amesh,all/psf,pres,norm,2,0,1sfe,1,1,pres,,100sfe,1,3,pres,,100sfe,1,6,pres,,100sfe,4,1,pres,,100,110,120,130!ex4.7 在梁单元上施加荷载finish/clear/prep7et,1,beam3k,1k,2,10l,1,2esize,,10lmesh,all/pnum,elem,1sfbeam,3,1,pres,50,100 sfbeam,5,1,pres,100 sfbeam,7,1,pres,50,100,,,0.2,0.1sfbeam,9,1,pres,100,,,,0.4,-1sfbeam,3,2,pres,50,100,,,0.2,0.1!ex4.8a 在线上施加面荷载finish/clear/prep7et,1,82blc4,,,10,30esize,5amesh,all/psf,pres,norm,2sfl,4,pres,10,60sfl,2,pres,60sftraneplot!ex4.8b 在线上施加面荷载finish/clear/prep7et,1,shell63wprota,,90blc4,,,10,30esize,5amesh,allsfl,3,pres,100sftraneplot!ex4.9 利用荷载梯度在直角坐标系下的施加方法finish/clear/prep7et,1,82blc4,,,10,60esize,2amesh,all/psf,pres,norm,2sfgrad,pres,,y,0,-5nsel,s,loc,x,0nsel,r,loc,y,0,40sf,all,pres,600sfgrad,pres,,y,30,-20nsel,s,loc,x,10sf,all,pres,0allseleplot!ex4.10利用荷载梯度在柱坐标系下的施加方法finish/clear/prep7csys,1k,1,10,-90k,2,10,90k,3,10,90,30l,1,2l,2,3adrag,1,,,,,,2ldele,2,,,1numcmp,allet,1,63esize,2amesh,allsfgrad,pres,1,y,-90,1sfa,all,2,pres,400sftran!可考察下述命令及其施加方法!①!sfgrad,pres,1,y,270,1sfa,all,2,pres,400sftran!②!local,12,1cscir,12,1sfgrad,pres,12,y,270,1sfa,all,2,pres,400sftran!ex4.11 ldread命令的示例finish/clear/filname,ldtest1/prep7et,1,plane82mp,ex,1,2.1e5mp,prxy,1,0.3blc4,,,10,40esize,2amesh,alld,all,allsfl,3,pres,100/solusolve/post1prrsolfinish/filename,ldtest2/solulsclear,allnsel,s,loc,y,0d,all,allallselldread,reac,,,,,ldtest1,rst solve/post1plnsol,s,y!ex4.12 初应力荷载finish/clear/filname,colu1/prep7et,1,plane82mp,ex,1,2e5mp,nuxy,1,0.3blc4,,,1,10esize,2amesh,allnsel,s,loc,y,0d,all,uyd,1,uxnsel,s,loc,y,10sf,all,pres,-10。

ansys 命令流

1. A，P1，P2，…，P17，P18（以点定义面）2. AADD，NA1，NA2，…NA8，NA9（面相加）3. AATT，MAT，REAL，TYPE，ESYS，SECN（指定面的单元属性）【注】ESYS为坐标系统号、SECN为截面类型号。

4. *ABBR，Abbr，String（定义一个缩略词）5. ABBRES，Lab，Fname，Ext（从文件中读取缩略词）6. ABBSAVE，Lab，Fname，Ext（将当前定义的缩略词写入文件）7. ABS，IR，IA，--，--，Name，--，--，FACTA（取绝对值）【注】*************8. ACCAT，NA1，NA2（连接面）9. ACEL，ACEX，ACEY，ACEZ（定义结构的线性加速度）10. ACLEAR，NA1，NA2，NINC（清除面单元网格）11. ADAMS，NMODES，KSTRESS，KSHELL【注】*************12. ADAPT, NSOLN, STARGT, TTARGT, FACMN, FACMX, KYKPS, KYMAC【注】*************13. ADD，IR, IA, IB, IC, Name, --,-- , FACTA, FACTB, FACTC（变量加运算）14. ADELE，NA1，NA2，NINC，KSWP（删除面）【注】KSWP =0删除面但保留面上关键点、1删除面及面上关键点。

15. ADRAG，NL1，NL2，…，NL6，NLP1，NLP2，…，NLP6（将既有线沿一定路径拖拉成面）16. AESIZE，ANUM，SIZE（指定面上划分单元大小）17. AFILLT，NA1，NA1，RAD（两面之间生成倒角面）18. AFSURF，SAREA，TLINE（在既有面单元上生成重叠的表面单元）19. *AFUN, Lab（指定参数表达式中角度单位）20. AGEN, ITIME, NA1, NA2, NINC, DX, DY, DZ, KINC, NOELEM, IMOVE（复制面）21. AGLUE，NA1，NA2，…，NA8，NA9（面间相互粘接）22. AINA，NA1，NA2，…，NA8，NA9（被选面的交集）23. AINP，NA1，NA2，…，NA8，NA9（面集两两相交）24. AINV，NA，NV（面体相交）25. AL，L1，L2，…，L9，L10（以线定义面）26. ALIST，NA1，NA2，NINC，Lab（列表显示面的信息）【注】Lab=HPT时，显示面上硬点信息，默认为空。

ansys隧道荷载结构法命令流及设计师对隧道设计的经验总结

/prep7*set,shuxiang,276170 !竖向荷载*set,zuoce,125500 !水平荷载*set,youce,-125500*set,jizhongli,-100000 !地震集中力*set,uxishu,136.966e6*0.5 !水平弹性抗力系数*set,vxishu,87.4348e6*0.5 !竖向弹性抗力系数!定义单元类型、实常数及材料性质et,1,beam3r,1,0.3,0.00225,0.3mp,ex,1,3.1e10 !弹性模量mp,dens,1,2500 !密度mp,prxy,1,0.167 !泊松比!创建几何模型画点k,1,0,2.55k,2,-2.04,1.53k,3,-2.6,-0.15k,4,-2.319721,-2.023237k,5,-1.648852,-2.71152k,6,0,-2.98k,7,1.648852,-2.71152k,8,2.319721,-2.023237k,9,2.6,-0.15k,10,2.04,1.53k,1000,0,0 !中心点!画隧道轮廓线larc,1,2,1000,2.55larc,2,3,1000,2.8larc,3,4,1000,6.4larc,4,5,1000,1.04969larc,5,6,1000,5.2larc,6,7,1000,5.2larc,7,8,1000,1.04969larc,8,9,1000,6.4larc,9,10,1000,2.8larc,10,1,1000,2.55!选择左边的所有线allsel !选择所有的实体lsel,s,loc,x,0,-5 !选择左边的所有线:x从0到－5 lcomb,all !合并所选的线lesize,all,,,30 !把上面合并的线等分30段!选择右边的所有线allsellsel,s,loc,x,0,5lcomb,alllesize,all,,,30allsel !全选!划分单元lmesh,allnplot !显示节点!加弹簧!y方向*do,i,23,40PSPRNG,i,TRAN,vxishu,,-0.3,, , !tran-直线的*enddoPSPRNG,2,TRAN,vxishu,,-0.3,, , !固定结构用的!x方向*do,i,15,26PSPRNG,i,TRAN,uxishu,-0.3,,, ,*enddo*do,i,37,48PSPRNG,i,TRAN,uxishu,0.3,,, ,*enddoallselfinish!进入求解层，施加荷载，定义荷载步等/solu!施加约束nsel,s,,,2 !选择约束的节点nsel,s,loc,x,0dnsel,r,loc,y,-2.98d,all,ux !施加水平方向的约束allselfcum,add,, !一定要，使荷载能叠加!竖向荷载*do,i,3,15a=-shuxiang*0.5*abs(nx(i)-nx(i+1)) f,i+1,fy,af,i,fy,a*enddo*do,i,47,59a=-shuxiang*0.5*abs(nx(i)-nx(i+1)) f,i+1,fy,af,i,fy,a*enddo!对未循环的节点施加竖向荷载a=-shuxiang*0.5*abs(nx(1)-nx(3)) f,3,fy,af,1,fy,2*af,60,fy,a!水平荷载*do,i,3,30a=zuoce*0.5*abs(ny(i)-ny(i+1))f,i+1,fx,af,i,fx,a*enddo*do,i,32,59a=youce*0.5*abs(ny(i)-ny(i+1))f,i+1,fx,af,i,fx,a*enddof,1,fx,jizhongli !加集中力!对未循环的节点施加水平荷载a=youce*0.5*abs(ny(31)-ny(2))f,31,fx,-af,32,fx,a!施加重力acel,,9.8solve !求解finish!进入后处理，定义荷载工况并组合，输出图片和文本文件/post1!显示弯矩图etable,mi,smisc,6etable,mj,smisc,12plls,mi,mj,-1!显示轴力图etable,fi,smisc,1etable,fj,smisc,7plls,fi,fj,1!显示变形图pldisp,1综合该讨论的各家建议：隧道设计流程各家观点：1、我们做设计的时候也是这样做的,先用同济曙光采用地层结构法模拟开挖过程,对喷锚支护及开挖方法进行分析,选择一个合适的开挖顺序,然后用荷载结构法计算初衬的变形, 和二衬的承载能力和裂缝宽度。

ansys命令流操作大全.

ansys——ANSYS命令流（Ⅰ）1. A，P1，P2，…，P17，P18（以点定义面）2. AADD，NA1，NA2，…NA8，NA9（面相加）3. AATT，MAT，REAL，TYPE，ESYS，SECN（指定面的单元属性）【注】ESYS为坐标系统号、SECN为截面类型号。

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 EPCSDGBLINK10 3D仅受拉或仅受压杆EDGB 模拟缆索的松弛及间隙LINK11 3D线性调节器EGB 模拟液压缸和大转动LINK180 3D有限应变杆EPCDFGB 另可考虑粘弹塑性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实例教学

山岭隧道受力ANSYS有限元分析实例教学目录一、问题重述 (1)二、模型的建立 (3)2.1模型绘制 (3)2.2模型参数选取 (3)2.3模型网格划分 (3)2.4计算外荷载（计算DK5+632断面） (4)2.5施加荷载与约束 (7)三、求解模型与受拉地基弹簧的修正 (8)四、求解结果 (10)4.1弯矩、轴力应力云图 (10)4.2关键节点内力 (10)五、附录 (11)附录1 全部节点等效荷载表 (11)附录2 全部节点内力表 (13)附录3 剪力图 (16)一、问题重述隧道起讫里程为DK4+843.5~DK6+430，全长1586.5m ，DK5+632处采用暗挖法施工，该断面的地层及结构等信息见下图。

(a) 纵断面图（单位：m ）(b) 横断面图（单位：cm ）里程 D K 5+632300°∠65°根据地质资料得：围岩级别为Ⅳ级，隧道上方土体重度依次从上往下取γ1 =18 kN/m3，γ2=23kN/m3。

请采用荷载-结构模式对该断面衬砌结构（仅二次衬砌）进行受力分析：（1）试求隧道围岩压力和有限元模型的等效节点力（不考虑重力）？（要求：单元长度取0.3m，画出单元和节点图，编制表格列出各节点的等效节点力）。

（2）采用有限元软件计算结构内力，绘制弯矩图和轴力图，列出特征部位的内力二、模型的建立2.1模型绘制在ANSYS建模，以二次衬砌中轴线为轮廓，隧道断面模型如下图：图2.1 隧道断面尺寸示意图（cm）二次衬砌采用Beam188梁单元模拟，地基弹簧采用Combin14弹簧单元模拟。

隧道纵向计算长度取为1m，二次衬砌参数选取如下表：2.3模型网格划分单元长度取0.3m，网格划分后，单元图、节点图分别如下：图2.2 离散化-节点图图2.3离散化-单元图2.4外荷载的计算（计算DK5+632断面） 2.4.1 验算坑道高度与跨度之比1232644.3032120.902 1.71232723.553212H B ++++==++++＜式中，H 表示坑道高度，B 表示坑道跨度根据我国《铁路隧道设计规范》，可以采用统计法计算。