ANSYS单元特性之命令流算例
ANSYS命令流教程
引用ANSYS 入门教程(26) - 网格划分高级技术(b)四、单元有效性检查不良的单元形状会导致不准确的结果，然而并没有判别单元形状好坏的通用标准，也就是说一种单元形状对一个分析可能导致不准确的结果，但可能对另一种分析的结果又是可接受的。
在计算过程中，ANSYS可能不出现单元形状警告信息，也可能会出现很多个单元形状警告信息，这都不能说明单元形状就一定会导致准确或不准确的结果，因此单元形状的好坏和结果的准确性完全依赖用户的判断和分析。
1. 单元形状参数限值设置命令：SHPP, Lab, VALUE1, VALUE2ANSYS 单元形状检查是缺省的，但控制单元形状检查的参数可以修改。
⑴Lab = ON：激活单元形状检查。
VALUE1 可取:ANGD：SHELL28 单元角度检查。
ASPECT：单元纵横比检查。
如四边形单元警告限值为20，错误限值为1E6；PARAL：对边平行度检查。
如无中间节点的四边形该项的警告限值为70°，如超过150°则给出错误信息。
MAXANG：最大角度检查。
无中间节点的四边形单元该项警告限值为155°，而其错误限值为179.9°；JACRAT：雅可比率检查。
简单地说，雅可比率表达了“单元”模拟“实际”的计算可靠性，比率越高越不可靠。
如h 单元的警告限值为30，超过30 单元形状就很不理想（与母单元形状相差甚远）。
WARP：歪曲率检查。
对于四边形面单元、壳单元或体单元的面等，当其严重歪曲时造成不好的单元形状，此值越高表示单元歪曲越严重。
也可用ALL 关闭或激活所有选项。
⑵Lab = WARN：仅激活警告模式，对超过错误限制的单元只给出警告信息而不致网格划分失败。
而Lab=ON 则一旦超过错误限制时将导致网格划分失败。
⑶Lab = OFF：完全关闭单元形状检查，可通过设置VALUE1 的值而关闭个别形状检查。
如VALUE1 可取ANGD、ASPECT、PARAL、MAXANG、JACRAT、WARP 及ALL 等。
ANSYS命令流
VDRAG,NA1,NA2,NA3,NA4,NA5,NA6,NLP1,NLP2,NLP3,NLP4,NLP5,NLP6体积（V olume）的建立是由一组面积（NA1~NA6），延某组线段（NL1~NL6）为路径，拉伸而成。
Menu Paths:Main Menu>Operate>Extrude/Sweep>Along LinesVROTAT,NA1,NA2,NA3,NA4,NA5,NA6,P AX1,P AX2,ARC,NSEG建立柱形体积，即将一组面（NA1~NA6）绕轴PAX1,PAX2旋转而成，以已知面为起点，ARC为旋转的角度，NSEG为整个旋转角度中欲分的数目。
Menu Paths:Main Menu>Operate>Extrude/Sweep>About Axis【例4-4】综合点、线、面、体积练习/PREP7/PNUM,KP,1 $/PUM,LINE,1/PNUM,AREA,1 $/PNUM,VOLU,1K,1,2 ！建立点1，坐标（2，0）K,2,3,4 ！建立点2，体坛（3，4）K,3,-0.5,3 ！建立点3，体坛（-0.5，3）K,4,-2,0.5 ！建立点4，体坛（-2，0.5）A,1,2,3,4 ！由点1，2，3，4生成面CSYS,1 ! 转为圆柱坐标系统K,5,4 ！建立点5，半径=4，角度=0K,6,4,45 ！建立点6，半径=4，角度=45K,7,4,45,5 ！建立点7，半径=4，角度=45，高度=5K,8,4,,5 ！建立点7，半径=4，角度=0，高度=5A,5,6,7,8 ！生成一个中心角为45度的圆柱面ADEL,1 ！删除面1LDEL,3 ！删除线1A,1,2,3,4 ！生成面ADEL,2A,5,6,8A,8,5,6CSYS ！转到卡式坐标K,,7,1,3/PNUM,KP,1V,8,5,6,94.8 用体素创建ANSYS对象这里先引入体素（Primitive）的概念，ANSYS中，体素指预先定义好的具有共同形状的面或体。
(完整版)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命令流使用方法
在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命令流学习笔记5-workbench中命令流的一些应用
ANSYS命令流学习笔记5-workbench中命令流的⼀些应⽤!ANSYS命令流学习笔记5workbench中命令流的⼀些应⽤学习重点：1. 定义单元类型2. 使⽤各向异性材料时，定义其单元为圆柱坐标系3. 有⾓度吊装时，定义吊装约束4. workbench、APDL的联合仿真案例如下：如下图模型，四个顶点通过杆件连接，进⾏吊装时的有限元分析。
1. 建⽴模型，设定必要的坐标系。
分成两个solid，内部的圆柱，剩余的矩形部分。
建⽴两个坐标系，分别⽤于指定各向异性材料的属性、吊装的固定点。
下图，建⽴圆柱坐标系，编号100，⽤于指定各向异性材料。
下图，建⽴直⾓坐标系，编号12，⽤于指定吊装固定点。
2. 建⽴named selection，⽅便在命令流中选择必要的元素。
下图，将四个吊装点，中间的圆柱，分别定义为任何名称，必须是英⽂才能⽤于APDL 命令中。
3. 定义边界条件，施加重⼒加速度，在static structural 下插⼊command(APDL) ，内容如下/prep7allscmsel,s,c1,elem !选择c1单元所有节点,既圆柱体的所有单元emodif,all,esys,100 !其坐标系转换为100坐标系，因为缠绕的各向异性材料必须在圆柱坐标系下定义单元的坐标。
!完成对各向异性材料的坐标系设定。
et,10,10 !定义编号为10的，link10单元r,10,0.01 !定义编号为10的实常数0.01，⽤于定义link10单元的截⾯积0.01mm^2*get,nmax,node,,num,max !获取node的最⼤数值，储存在nmax名称的变量⾥csys,12 !调⽤csys12坐标系n,nmax+1, !csys坐标原点建⽴node，后⾯会将其固定，既吊装的固定点mat,1type,10 !选取编号10的单元类型real,10 !选取编号10的实常数cmsel,s,kk1,node !选择kk1点，kk1已经定义为named selection*get,k1,node,,num,max !获取已选节点中的节点数最⼤值，既kk1的节点编号，取值为k1cmsel,s,kk2,node ! k2,k3,k4⽅法类型k1*get,k2,node,,num,maxcmsel,s,kk3,node*get,k3,node,,num,maxcmsel,s,kk4,node*get,k4,node,,num,maxalls !全选e,nmax+1,k1e,nmax+1,k2e,nmax+1,k3e,nmax+1,k4 !建⽴四个link单元d,nmax+1,allalls/sol此外，下图中单位必须保持⼀致，不然计算很容易出问题。
ansys命令流说明
24. AINV，NA，NV（面体相交）
25. AL，L1，L2，…，L9，L10（以线定义面）
26. ALIST，NA1，NA2，NINC，Lab（列表显示面的信息）
【注】Lab=HPT时，显示面上硬点信息，默认为空。
【注】*************
13. ADD，IR, IA, IB, IC, Name, --,-- , FACTA, FACTB, FACTC（变量加运算）
14. ADELE，NA1，NA2，NINC，KSWP（删除面）
【注】KSWP =0删除面但保留面上关键点、1删除面及面上关键点。
RETAIN=ON（对于全是四边形的网格，细化不会将三角形引入）、OFF（可能将三角形引入）
51. AREVERSE，ANUM，NOEFLIP（将面的法线方向反向）
52. AROTAT，NL1，NL2，NL3，NL4，NL5，NL6，PAX1，PAX2，ARC，NSEG（绕轴旋转生成面）
78. BFK，KPOI，Lab，VAL1，VAL2，VAL3，PHASE（在关键上施加体载荷）
79. BFKDELE，KPOI，Lab（删除关键点上的体载荷）
80. BFKLIST，KPOI，Lab（列表显示线关键点上的体载荷）
81. BFL，LINE，Lab，VAL1，VAL2，VAL3，PHASE（在线上施加体载荷）
15. ADRAG，NL1，NL2，…，NL6，NLP1，NLP2，…，NLP6（将既有线沿一定路径拖拉成面）
16. AESIZE，ANUM，SIZE（指定面上划分单元大小）
17. AFILLT，NA1，NA1，RAD（两面之间生成倒角面）
ANSYS命令流
5.3 元素形状定义元素形状在2-D结构中可分为四边形和三角形，在3-D结构中可分为六面体和角锥体。
当实体模型进行对应网格划分时，2-D及3-D结构所产生的元素必为四边形及六面体，当无法进行对应网格化时，程序会自动用自由网格化，所以2-D结构将自行以四边形和三角形的混合方式进行，3-D结构以角锥体方式进行。
网格化，有默认尺寸大小，也就是说不给定线段和网格数目，仍然可以进行网格划分，但不一定能满足设计者的要求。
元素大小基本上在线段上定义，可用线段数目和线段长度来划分，通常以线段数目分割比较方便。
分割时可采用均分或不均分，不均分以线段方向或中间为准，根据数定义可得到渐增或渐减的效果。
除此之外，也可以以整体对象为基准，确定网格的大小。
此外在自由网格化一般不需要定义线段的数目及大小，程序将提供智能化控制（SMRTSIZE）；而指定线段进行元素数目及大小的声明，大多用于对应网格化。
5.4网格划分工具网格划分工具是网格控制的一种快捷方式，它能方便的实现单元属性控制、智能网格划分控制、尺寸控制、自由网格划分和对应网格划分、执行网格划分、清除网格划分以及局部细分。
参见图5-2。
程序默认为自由网格划分，元素形状以四边形、六面体为准优先，三角形、角锥次之，网格化时，如果实体模型能够对应网格化，而且相对应边长度不是差的很多，则必以对应网格化优先考虑进行。
网格划分工具中，我们一般只用它的一两组功能，即可达到要求。
这里有必要知道尺寸控制的优先级。
缺省单元尺寸控制：●对线划分的指定被最先考虑●关键点附近的单元尺寸作为第二级考虑对象●总体单元尺寸作为第三级考虑对象●缺省尺寸最后考虑智能单元尺寸的优先顺序●对线划分的指定被最先考虑●关键点附近的单元尺寸作为第二级考虑对象，当考虑到小的几何特征和曲率时，可以忽略它●总体单元尺寸作为第三级考虑对象，当考虑到小的几何特征和曲率时，可以忽略它●智能单元尺寸设置最后考虑【例5-1】综合练习，如图5-3所示/PREP7ET,1,PLANE42K,1,,-2.5 $K,2,6,-2.5 $K,3,,2.5 $K,4,6,2.5CSYS,1K,5,10,-30 $K,6,10,30CSYS,0SA VEL,1,2 $L,4,3CSYS,1,L,2,4 $L,5,6CSYS,0A,1,2,4,3A,2,5,6,4SA VE ! 可以有RESUME命令回复到当前点AMESH,ALL ! 此时网格以系统默认的尺寸进行自由划分，但由于面积符合对应网格化的要求，所以会进行对应网划分。
ANSYS命令流
ANSYS命令流一、建模/PREP7 进入前处理器。
（执行命令前要先输这个，否则会悲剧的提示错误。
）ET,1,SOLID45 定义单元类型1为SOLID45.MP,EX,1,3.45E8 定义材料属性1的弹性模量为3.45e8MP,PRXY,1,0.2 定义材料属性1的泊松比为0.2MP,DENS,1,2.653061E-3 定义材料属性1的密度为2.653061E-3R,2,5.4438,0.006868 定义实常数2，面积为5.4438，初应变为0.006868 （这个实常数定义对应于LINK8单元，其他的单元要查相应的格式。
）K,601,30,-94.4835,-1298.9 建立关键点601，坐标为后面的三个值。
KGEN,2,601,610,1,40,0,0,20 复制节点601至610一次，每个节点x坐标增量为40，y、z增量为0，节点编号增量为20。
L,902,903 连接节点901、902成线。
LGEN,13,127,,,50,,,3 复制线的命令流。
将线127复制12次，x 向位置增量为50，y、z向为0，线编号增量为3。
LSEL,S,LINE,,127,152 从所有图素中选择编号125至152间所有的线。
LSEL,S,LINE,,125 从所有图素中选择编号为125的线。
LSEL,A,LINE,,126 在上面已经选择线125的基础上，添加选择编号为126的线。
A,101,102,103,104,105,106,107 将上述7个关键点（最多10个）练成面。
VOFFST,113,175,400 将面113沿着垂直于面的方向偏移175，节点编号增量为400。
V,301,3301,306,307, 401, 3401, 406, 407 将上述八个关键点连结成体。
（注意前四个点和后四个点的选择顺序需一致，即同为顺时针或逆时针）-注意：关键点个数最多8个。
V,304,3304,305, 404,3404,405 将上述六个关键点连结成体。
ANSYS单元特性之命令流算例
EX1.1 （LINK1）（1）进入后处理模块,显示节点位移和杆件内力MID_NO DE = NODE (A/2,-B,0 )! 寻找距离位置（A/2,-B,0）最近的点，存入MID_NODE*GET,DISP,NODE,MID_NO DE,U,Y！提取节点MI D_NOD E上的位移UY，若果已知要求的节点，直接提取即可。
LEFT_E L = ENEARN (MID_NO DE)! 需找距离节点MID_N ODE最近的单元，存入LEFT_ELETABLE,STRS,L S,1! 用轴向应力S AXL的编号“LS,1”定义单元表S TRS*GET,STRSS,ELEM,LEFT_E L,ETAB,STRS! 从单元表STRS中提取LEFT_EL单元的应力结果，存入变量STRSS。
注意：提取的轴向应力结果具体到指定的单元。
（2）申明数组,提取计算结果,并比较计算误差*DIM,LABEL,CHAR,2！定义2个元素的字符型数组LABE L*DIM,VALUE,,2,3！定义2*3的数值型数组V ALU ELABEL(1) = 'STRS_M Pa','DEF_mm' ! 给字符型数组的第1个元素赋值*VFILL,VALUE(1,1),DATA,1,-0.05498! 给其他数值型数组中的元素赋值*VFILL,VALUE(1,2),DATA,STRSS,DISP*VFILL,VALUE(1,3),DATA,ABS(STRSS/1 ) ,ABS( DISP /0.05498)/OUT,EX1_1,out ！将输出内容重定向到文件EX1_1.out/COM ! 以注释形式输出内容/COM,------------------- EX1.1 RESULT S COMPAR ISON---------------------/COM,/COM, | TARGET | ANSYS | RATIO/COM,*VWRITE,LABEL(1),VALUE(1,1),VALUE(1,2),VALUE(1,3)(1X,A8,' ',F10.3,' ',F10.3,' ',1F5.3)/COM,----------------------------------------------------------------/OUT ! 结束数据重定向，关闭输出文件FINISH*LIST,EX1_1,out ! 列表显示文件EX1_1.out的内容EX1.2 （LINK1）/PNUM, NODE,1！打开节点编号显示/NUMBER, 2！只显示编号，不使用色彩列表显示节点位移和单元的计算结果PRDISP! 列表显示节点位移值计算结果ETABLE, MFORX,SMISC,1！以杆单元的轴力为内容，建立单元表M F ORXETABLE, SAXL, LS, 1 ！以杆的轴向应力为内容，建立单元表S A XLETABLE, EPELAX L, LEPEL, 1! 以杆单元的轴向应变为内容，建立单元表E PELAX LPRETAB! 显示单元表中的计算结果/NUMBER, 0！显示编号，并使用彩色PLETAB, MFORX！用色度图显示单元表MF ORX中杆件轴力图EX1.3 （LINK1）NSEL,S,LOC,Y,1.0 ！选择所有位于Y=1.0位置上的节点FSUM！累计叠加选择集中所有节点上的反力*GET,REAC_1,FSUM,,ITEM,FY ！将累加结果中的FY（Y方向的力）保存到变量R E AC_1中EX1.4 （LINK1）R,1,65e-6！定义第1类实常数，杆件截面面积为65mm^2，在转化为国际单位制时操作TREF,70 ! 设定参考温度为70度BFUNIF,TEMP,80 ! 温度从原来的70度均匀上升到80度（TREF+10）EX1.5 （PLANE42 AND CONT AC26）ETABLE,STRSX,S,X！定义X方向的应力为单元表S TRS X*GET,STRSSX,ELEM,3,ETAB,STRSX！从单元表STRS X中提取3号单元的X向应力，存入STRS S X。
(完整版)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, 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杆单元(Large 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命令流实例
ansys命令流实例1.结构静力分析－－梁分析/PREP7 !进入前处理器ET,1,BEAM3 !定义单元类型R,1,0.25,0.0052,0.5 !定义实常数MP,EX,1,210E6MP,PRXY,1,0.3 !定义材料属性N,1,0N,2,1N,3,5N,4,7N,5,7N,6,9N,7,11N,8,11N,9,13N,10,14 !定义节点E,1,2E,2,3E,3,4E,5,6E,6,7E,8,9E,9,10 !生成单元CP,1,UX,4,5CP,2,UY,4,5CP,3,UX,7,8CP,4,UY,7,8 !耦合节点FINISH/SOLU !进入求解器D,2,UXD,2,UYD,3,UYD,6,UYD,9,UY !施加位移约束F,10,FY,-4 !施加集中约束SFBEAM,1,1,PRES,4,4SFBEAM,2,1,PRES,4,4 !施加均布力SOLVE !求解/POST1 !进入后处理器PLDISP !绘制结构变形图PRDISP !列出各节点的位移ETABLE,IMOMENT,SMISC,6ETABLE,JMOMENT,SMISC,12ETABLE,ISHEAR,SMISC,2ETABLE,JSHEAR,SMISC,8 !将节点弯矩、剪力制表PRETAB !列表显示单元的弯矩、剪力/TITLE,SHEAR FORCE DISTRIBUTION !设置剪力分布图的标题PLLS,ISHEAR,JSHEAR !绘制剪力分布图/TITLE,BENDING MOMENT IDSTRIBUTION !设置弯矩分布图的标题PLLS,IMOMENT,JMOMENT !绘制弯矩分布图流固耦合实例RAD=0.8 !底面半径H=1G=9.8OMEGAR=2ROU=1000 !定义参数变量/PREP7 !进入前处理器ET,1,FLUID79 !选择单元类型KEYOPT,1,3,1 !设置单元关键字MP,EX,1,2E9 !设置杨氏模量MP,DENS,1,ROU !设置材料密度K,1K,2,RADK,3,RAD,HK,4,,H !生成关键点A,1,2,3,4 !连接关键点生成面积LESIZE,ALL,,,10 !设置网格划分精度AMESH,ALL !将面积划分网络/SOLU !进入求解器DL,2,,UXDL,1,,UYNSEL,S,LOC,XDSYM,SYMM,XD,ALL,UXNSEL,ALL !施加位移约束ACEL,,GOMEGA,,OMEGAR !施加惯性力SOLVE !求解/POST1SET,LAST !进入通用后处理器PLNSOL,U,X,0,1 !绘制应力云图UCENT=UY(22)UEDGE=UY(12)UELEV=UEDGE-UCENT !提取节点位移结构静力分析－－壳结构内力分析LENGTH=100YOUNG=200000THICKNESS=2FORCE=1000DENSITY=9E-6 !将材质、载荷、板的几何尺寸等参数化/PREP7 !进入前处理器MP,EX,1,YOUNGMP,NUXY,1,0.3MP,DENS,1,DENSITY !定义材质ET,1,SHELL63 !定义单元类型R,1,THICKNESS,THICKNESS,THICKNESS,THICKNESS !定义实常数!构建结构的几何模型K,1,0,0K,2,LENGTH,0K,3,LENGTH,LENGTHK,4,0,LENGTH !定义关键点A,1,2,3,4!通过关键点生成面LSEL,ALLLESIZE,ALL,,,16AMESH,ALL !设定网格划分参数，划分网格FINISH/SOLU !进入求解器NSEL,S,LOC,X,0,0D,ALL,ALL,0 !选择X＝0的节点将其固定NSEL,S,LOC,X,LENGTH,LENGTHD,ALL,ALL,0 !选择X＝LENGTH的节点将其固定NSEL,S,LOC,X,0.5*LENGTH,0.5*LENGTHNSEL,R,LOC,Y,0.5*LENGTH,0.5*LENGTHF,ALL,FZ,FORCEALLSEL !捕捉板的中心点并在中心点处施加集中力荷载SOLVE !求解FINISH/POST1 !进入后处理器/DSC,,10PLNSOL,U,Z,0,1 !绘图显示板的竖向变形NSEL,ALL !提取板的最大竖向变形NSORT,U,Z,1,1 !将节点的位移绝对值以升序排序*GET,MAXDEFLECTION,SORT,0,MAX !提取位移最大值并赋给变量NSEL,S,LOC,X,0NSEL,A,LOC,X,LENGTH,LENGTH !选择固定边节点NSORT,S,EQV,1,1 !将节点等效应力的绝对值以升序排序*GET,MAXSTRESS,SORT,0,MAX !提取等效应力最大值并赋给变量结构静力分析－－平面桁架分析/PREP7 !进入前处理器ET,1,LINK1 !选择单元R,1,0.1 !定义实常数MP,EX,1,30E6MP,PRXY,1,0.3 !定义材料属性N,1,0N,2,4N,3,8N,4,12N,5,0,3N,6,4,3N,7,8,3 !生成节点E,1,2E,2,3E,3,4E,4,7E,3,7E,2,7E,2,6E,2,5E,1,5E,5,6E,6,7 !生成单元FINISH/SOLU !进入求解器D,1,UXD,1,UYD,5,UX !实加位移约束F,2,FY,-15F,3,FY,-15F,4,FY,-15 !施加集中力SOLVE !求解/POST1 !进入能用后处理器PRESOL,FORC !列表显示反力。
ANSYS命令流详解（超全）
ANSYS命令流详解（超全）一、定义材料号及特性mp,lab, mat, co, c1,…….c4lab: 待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens）ex: 弹性模量nuxy: 小泊松比alpx: 热膨胀系数reft: 参考温度reft: 参考温度prxy: 主泊松比gxy: 剪切模量mu: 摩擦系数dens: 质量密度mat: 材料编号（缺省为当前材料号）c 材料特性值，或材料之特性，温度曲线中的常数项c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数二、定义DP材料：首先要定义EX和泊松比：MP，EX，MAT，……MP，NUXY，MAT，……定义DP材料单元表（这里不考虑温度）：TB，DP，MAT进入单元表并编辑添加单元表：TBDATA，1，CTBDATA，2，ψTBDATA，3，……如定义：EX=1E8，NUXY=0.3，C=27，ψ=45的命令如下：MP，EX，1，1E8MP，NUXY，1，0.3TB，DP，1TBDATA，1，27TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,deg三、单元生死载荷步!第一个载荷步TIME,... !设定时间值（静力分析选项）NLGEOM,ON !打开大位移效果NROPT,FULL !设定牛顿－拉夫森选项ESTIF,... !设定非缺省缩减因子（可选）ESEL,... !选择在本载荷步中将不激活的单元EKILL,... !不激活选择的单元ESEL,S,LIVE !选择所有活动单元NSLE,S !选择所有活动结点NSEL,INVE !选择所有非活动结点（不与活动单元相连的结点）D,ALL,ALL,0 !约束所有不活动的结点自由度（可选）NSEL,ALL !选择所有结点ESEL,ALL !选择所有单元D,... !施加合适的约束F,... !施加合适的活动结点自由度载荷SF,... !施加合适的单元载荷BF,... !施加合适的体载荷SAVESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D,F,SF和BF命令得到更详细的解释。
ANSYS命令流
令查看。
建立元素前必须先行定义使用者欲选择的元素型号、元素材料特性、元素几何特性等，为了程序的协调性一般在/PREP7后，就定义元素型号及相关资料，只要在建立元素前说明使用哪种元素即可。
相关命令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/DeleteMP,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>IsotropicR,NSET,R1,R2,R3,R4,R5,R6定义”实常数”，即某一单元的补充几何特征，如梁单元的面积，壳单元的厚度。
所带的的参数必须与元素表的顺序一致。
Menu paths:Main Menu>Preprocessor>Real ConstantsE,I,J,K,L,M,N,O,P定义元素的连接方式，元素表已对该元素连接顺序作出了说明，通常2-D平面元素节点顺序采用顺时针逆时针均可以，但结构中的所有元素并不一定全采用顺时针或逆时针顺序。
ansys 温度荷载实例分析命令流
!XXXXXXXXXXXXXXXXXX!定义单元类型!XXXXXXXXXXXXXXXXXXFINI/CLE/TITLE,QSHang Pylon Test/UNITS,SICSYS,0!*afun,deg/PREP7/triad,on/view,1,1,1,1 !定义ISO查看/NERR,1000000et,1,plane55 !定义二维热单元et,2,surf153mp,kxx,1,3.3mp,kyy,1,3.3 !设定砼的导热系数mp,kzz,1,3.3!mp,dens,1,2400 !定义砼的密度mp,c,1,925 !定义砼的比热mp,prxy,1,0.2 !泊松比C50mp,ex,1,3.45e10 !混凝土的弹性模量C50 mp,alpx,1,1e-5!mptemp,1,0,2643,2750,2875 !定义铸钢的热性能!mpdata,kxx,2,1.44,1.54,1.22,1.22!mpdata,enth,2,0,128.1,163.8,174.2!mpplot,kxx,2!mpplot,enth,2!XXXXXXXXXXXXXXXXXXX!建立混凝土几何模型!XXXXXXXXXXXXXXXXXXX!************plane1K,1,0,0,0k,2,0,3.5,0k,3,6.885,3.5k,4,6.885,4k,5,21.2,4k,6,23.7,0l,1,2l,2,3l,3,4l,4,5local,11,1,21.2,0,0,,,,4/2.5,1 l,5,6CSYS,0k,7,18.9,0k,8,18.9,2k,9,19.4,2.5k,10,21,2.5k,11,22.2,0l,1,7l,7,8l,8,9l,9,10local,12,1,21,0,0,,,,2.5/1.2,1 l,10,11CSYS,0l,11,6al,all!************plane2 asel,nonelsel,nonek,12,0.8,0.25k,13,0.3,0.75k,14,0.3,2k,15,0.8,2.5k,16,5.2,2.5k,17,5.7,2k,18,5.7,0.75k,19,5.2,0.25a,12,13,14,15,16,17,18,19k,20,6.8,0.25k,21,6.3,0.75k,22,6.3,2k,23,6.8,2.5k,24,11.2,2.5k,25,11.7,2k,26,11.7,0.75k,27,11.2,0.25a,20,21,22,23,24,25,26,27k,28,12.8,0.25k,29,12.3,0.75k,30,12.3,2k,31,12.8,2.5k,32,15.3,2.5k,33,15.8,2k,34,15.8,0.75k,35,15.3,0.25a,28,29,30,31,32,33,34,35asel,scm,plane,areaallsasba,1,plane!************画网格aatt,1,1,AESIZE,all,0.02 !将所有面网格划分时的单元尺寸设置为0.02m MSHAPE,0,2D !利用四边形单元进行网格划分MSHKEY,2 !采用映射网格划分单元AMESH,ALL !对所有网格进行划分单元allscm,area1,areaarsym,x,area1,,,,0,0allsarsym,y,all,,,,0,0asel,allaglue,all!Nummrg,elem!Nummrg,node!Nummrg,kpnummrg,allnumcmp,all!**************定义组件allslsel,s,length,,6.885lsel,a,length,,0.5lsel,a,length,,14.32lsel,r,loc,y,0,100cm,Eligne,linensll,scm,Enode,nodeTYPE,2ESURFallslsel,s,length,,6.885lsel,a,length,,0.5lsel,a,length,,14.32lsel,r,loc,y,-100,0cm,Wligne,linensll,scm,Wnode,nodeTYPE,2ESURFallslsel,s,length,,5.173lsel,r,loc,x,0,100cm,Sligne,linensll,scm,Snode,nodeTYPE,2ESURFallslsel,s,length,,5.173lsel,r,loc,x,-100,0cm,Nligne,linensll,scm,Nnode,nodeALLSTYPE,2ESURF!**************计算日照时间内，即14个小时内的空气的温度值*DIM,t_outdoor,,13 !定义室外空气温度数组为t_outdoor*DO,t,6,18 !给室外空气温度数组t_outdoor赋值*SET,t_outdoor(t-5),31.5+5.5*sin(3.1415926535898*(t-9)/12)*ENDDO!**************日照强度*DIM,E_insolation,,13 !东边日照*DIM,S_insolation,,13 !东边日照*DIM,W_insolation,,13 !西边日照*DIM,N_insolation,,13 !北边日照E_insolation(1) =240 $S_insolation(1) =33 $ W_insolation(1) = 28 $ N_insolation(1) =89E_insolation(2) =521 $S_insolation(2) =80 $W_insolation(2) =69 $ N_insolation(2) = 160E_insolation(3) =629 $S_insolation(3) = 115 $W_insolation(3) =99 $ N_insolation(3) = 166E_insolation(4) =640 $S_insolation(4) =148 $W_insolation(4) = 127 $ N_insolation(4) = 152E_insolation(5) =536 $S_insolation(5) =184 $W_insolation(5) = 142 $ N_insolation(5) = 142E_insolation(6) =374 $S_insolation(6) =229 $W_insolation(6) = 153 $ N_insolation(6) = 153E_insolation(7) =178 $S_insolation(7) =240 $W_insolation(7) = 153 $ N_insolation(7) = 153E_insolation(8) =178 $S_insolation(8) =229 $W_insolation(8) = 322 $ N_insolation(8) = 153E_insolation(9) =165 $S_insolation(9) =184 $W_insolation(9) = 461 $ N_insolation(9) = 142E_insolation(10) =148 $S_insolation(10) =148 $W_insolation(10) =550 $ N_insolation(10) = 152E_insolation(11)=115 $S_insolation(11)=115 $W_insolation(11) =541 $ N_insolation(11) =166E_insolation(12)=80 $S_insolation(12)=80 $W_insolation(12)= 448 $ N_insolation(12) = 160E_insolation(13)=33 $S_insolation(13)=33 $W_insolation(13)= 206 $ N_insolation(13) = 89!**************综合换热系数h=12.47+3.33*1.5!**************综合大气温度*DIM,Et_synthetiser,,13*DIM,St_synthetiser,,13*DIM,Wt_synthetiser,,13*DIM,Nt_synthetiser,,13*DO,x,1,13,1*SET,Et_synthetiser(x),t_outdoor(x)+(E_insolation(x)/17.465)*0.65*SET,St_synthetiser(x),t_outdoor(x)+(S_insolation(x)/17.465)*0.65*SET,Wt_synthetiser(x),t_outdoor(x)+(W_insolation(x)/17.465)*0.65*SET,Nt_synthetiser(x),t_outdoor(x)+(N_insolation(x)/17.465)*0.65*ENDDO!**************求解/soluantype,trans !设置求解类型瞬态热分析trnopt,full ! 指定瞬态分析的求解方法为完全法timint,on ! 打开时间积分效应tunif,27.6 !设置构件的初始温度为20摄氏度!tref,20outres,all ! 将除SV AR和LOCI以外的所有计算数据写入数据库和文件中!antype,4!autots,on!trnopt,full!lnsrch,on!outres,basic,last!outpr,basic,lastl=0*do,m,1,13,1*do,r,3600,3600,3600time,(m-1)*3600+rl=l+1nsubst,1,100,1 !设置每个荷载的子步数数为1，最大值为100，最小值为1 autots,on ! 打开自动时间步长跟踪eqslv,JCG !指定方程求解器为JCGkbc,0 !使用递增方式加载线性荷载!施加对流荷载!alls!nsel,s,,,Enode!d,all,temp,Et_synthetiser(m)!alls!nsel,s,,,Wnode!d,ALL,temp,Wt_synthetiser(m)!alls!nsel,s,,,Snode!d,ALL,temp,St_synthetiser(m)!alls!nsel,s,,,Nnode!d,ALL,temp,Nt_synthetiser(m)!bf,Enode,hgen,E_insolation(m)!bf,Wnode,hgen,W_insolation(m)!bf,Snode,hgen,S_insolation(m)!bf,Nnode,hgen,N_insolation(m)allsSFl,Eligne,conv,17.456,,Et_synthetiser(m) SFL,Wligne,conv,17.456,,Wt_synthetiser(m) SFL,Sligne,conv,17.456,,St_synthetiser(m) SFL,Nligne,conv,17.456,,Nt_synthetiser(m)allslsel,s,,,Eligneesll,sesel,r,type,,1sfe,all,,HFLUX,,0.65*E_insolation(m)allslsel,s,,,Wligneesll,sesel,r,type,,1sfe,all,,HFLUX,,0.65*W_insolation(m)allslsel,s,,,Sligneesll,sesel,r,type,,1sfe,all,,HFLUX,,0.65*S_insolation(m)allslsel,s,,,Nligneesll,sesel,r,type,,1sfe,all,,HFLUX,,0.65*N_insolation(m)!alls!nsel,s,,,Enode!esln,s!esel,r,type,,1!sfe,all,,conv,,0.65*E_insolation(m)!alls!nsel,s,,,Wnode!esln,s!esel,r,type,,2!sfe,all,,conv,,0.65*W_insolation(m)!alls!nsel,s,,,Snode!esln,s!esel,r,type,,2!sfe,all,,conv,,0.65*S_insolation(m)!alls!nsel,s,,,Nnode!esln,s!esel,r,type,,2!sfe,all,,conv,,0.65*N_insolation(m)!sfa,5,,conv,10,t_outdoor(m) !给面5定义当前荷载步的对流换热系数及周围环境温度!sfa,21,,conv,10,t_outdoor(m) !给面21定义当前荷载步的对流换热系数及周围环境温度!sfa,18,,conv,10,t_outdoor(m) !给面18定义当前荷载步的对流换热系数及周围环境温度!sfa,10,,conv,10,t_outdoor(m) !给面10定义当前荷载步的对流换热系数及周围环境温度!施加太阳辐射荷载!bfv,1,HGEN,heat_eq(1,m)+long_wave(1) !给体1施加当前荷载步的等效生热率!bfv,4,HGEN,heat_eq(2,m)+long_wave(2) !给体4施加当前荷载步的等效生热率!bfv,3,HGEN,heat_eq(3,m)+long_wave(3) !给体3施加当前荷载步的等效生热率!bfv,2,HGEN,heat_eq(4,m)+long_wave(4) !给体2施加当前荷载步的等效生热率!bfe,all,hgen,1,1e10allssolve*enddo*enddo!*************************进入热应力求解/PREP7ETCHG,TTS !热到结构分析转换/SOLantype,4 !瞬态timint,1,struct !Turns on transient effectstimint,0,thermtimint,0,magtimint,0,elecautots,on !打开自动时间步长，有利于非线性收敛tref,20trnopt,full !完全瞬态nropt,full !设置牛顿-拉普森选项kbc,0 !渐变荷载nlgeom,on !打开大变形效应LUMPM,0PSTRES,ONNROPT,FULL, ,on!定义热应力计算参考温度TREF,20,!*************************边界条件DSYM,SYMM,Z,0csys,0LSEL,ALLDL,ALL,,UX,0DL,ALL,,UY,0DL,ALL,,UZ,0*DO,I,1,1355LDREAD,TEMP,,,1*i, ,'model_BD1','rth',' ' !读入热分析的计算结果OUTRES,ALL,ALL,TIME,1*iDELTIM,1,1,1SOLVE*ENDDO*DO,I,1,500LDREAD,TEMP,,,1355+10*i, ,'model_BD1','rth',' ' !读入热分析的计算结果OUTRES,ALL,ALL,TIME,1355+10*iDELTIM,10,10,10SOLVE*ENDDO*DO,I,1,300LDREAD,TEMP,,,6355+50*i, ,'model_BD1','rth',' ' !读入热分析的计算结果OUTRES,ALL,ALL,TIME,6355+50*iDELTIM,50,50,50SOLVE*ENDDO*DO,I,1,300LDREAD,TEMP,,,21355+100*i, ,'model_BD1','rth',' ' !读入热分析的计算结果OUTRES,ALL,ALL,TIME,21355+100*iDELTIM,100,100,100 SOLVE*ENDDOSET,LASTCSYS,0PATH,zdirec,2,50PPA TH,1,,0,0,0.17PPA TH,2,,0,0.7,0.17 PDEF,,S,xPLPATH,Sx/IMAGE,SA VE,zpath1,BMPSET,LASTCSYS,0PATH,zdirec,2,50PPA TH,1,,-1.72,0.4,0.17 PPA TH,2,,1,0.4,0.17 PDEF,,S,xPLPATH,Sx/IMAGE,SA VE,zpath2,BMPSET,LASTCSYS,0PATH,zdirec,2,50PPA TH,1,,0,0.7,0PPA TH,2,,0,0.7,0.35 PDEF,,S,zPLPATH,Sz/IMAGE,SA VE,zpath3,BMPSET,LASTCSYS,0PATH,zdirec,2,50PPA TH,1,,-1.72,0.7,0.17 PPA TH,2,,1,0.7,0.17 PDEF,,S,zPLPATH,Sz/IMAGE,SA VE,zpath4,BMPvsel,s,,,1 !选择编号为1的体元素eslv,r !选择当前所选体元素上的所有单元nsle,r !选择当前所选单元上的所有节点*get,node_total,node,0,count ! 提取当前所选节点的个数*get,node_num1(1),node,0,num,min ! 提取当前所选节点中节点的最小编号，并赋值给变量node_num1(1)*do,n,2,node_total,1*get,node_num1(n),node,node_num1(n-1),nxth ! 将当前所选择的节点的编号赋值给数组node_num1*enddo*do,f,1,node_total,1*get,node_temp1(f),node,node_num1(f),temp ! 提取当前所选节点的节点温度值，并赋值给数组node_temp1*enddototal=0*do,f,1,node_total,1total=total+node_temp1(f) ! 计算当前所选节点的温度值之和，并赋值给变量total*enddotarea=total/node_total+273 ! 计算当前所选节点温度值的平均值，并赋值给tareatsky=t_outdoor(m)-6+273 !计算当前时刻的天空温度值，并赋值给tsky*if,fuhao(1),eq,100,thenvalue_4=value_3*((tsky**4-tarea**4)*(1+cos_ph(1))/2) !计算与天空之间值长波辐射净强度*elsevalue_4=value_3*((t_earth(m)**4-tarea**4)*(1-cos_ph(1))/2) !计算与地面之间值长波辐射净强度*endif*set,long_wave(1,l),value_4/0.03 !计算长波辐射净强度的等效生热率allselvsel,s,,,4 !选择编号为4的体元素eslv,r !选择当前所选体元素上的所有单元nsle,r !选择当前所选单元上的所有节点*get,node_total,node,0,count ! 提取当前所选节点的个数*get,node_num1(1),node,0,num,min ! 提取当前所选节点中节点的最小编号，并赋值给变量node_num1(1)*do,n,2,node_total,1*get,node_num1(n),node,node_num1(n-1),nxth ! 将当前所选择的节点的编号赋值给数组node_num1*enddo*do,f,1,node_total,1*get,node_temp1(f),node,node_num1(f),temp ! 提取当前所选节点的节点温度值，并赋值给数组node_temp1*enddototal=0*do,f,1,node_total,1total=total+node_temp1(f) ! 计算当前所选节点的温度值之和，并赋值给变量total*enddotarea=total/node_total+273 ! 计算当前所选节点温度值的平均值，并赋值给tareatsky=t_outdoor(m)-6+273 !计算当前时刻的天空温度值，并赋值给tsky*if,fuhao(2),eq,100,thenvalue_4=value_3*((tsky**4-tarea**4)*(1+cos_ph(2))/2) !计算与天空之间值长波辐射净强度*elsevalue_4=value_3*((t_earth(m)**4-tarea**4)*(1-cos_ph(2))/2) !计算与地面之间值长波辐射净强度*endif*set,long_wave(2,l),value_4/0.03 !计算长波辐射净强度的等效生热率allselvsel,s,,,3 !选择编号为3的体元素eslv,r !选择当前所选体元素上的所有单元nsle,r !选择当前所选单元上的所有节点*get,node_total,node,0,count ! 提取当前所选节点的个数*get,node_num1(1),node,0,num,min ! 提取当前所选节点中节点的最小编号，并赋值给变量node_num1(1)*do,n,2,node_total,1*get,node_num1(n),node,node_num1(n-1),nxth ! 将当前所选择的节点的编号赋值给数组node_num1*enddo*do,f,1,node_total,1*get,node_temp1(f),node,node_num1(f),temp ! 提取当前所选节点的节点温度值，并赋值给数组node_temp1*enddototal=0*do,f,1,node_total,1total=total+node_temp1(f) ! 计算当前所选节点的温度值之和，并赋值给变量total*enddotarea=total/node_total+273 ! 计算当前所选节点温度值的平均值，并赋值给tareatsky=t_outdoor(m)-6+273 !计算当前时刻的天空温度值，并赋值给tsky*if,fuhao(3),eq,100,thenvalue_4=value_3*((tsky**4-tarea**4)*(1+cos_ph(3))/2) !计算与天空之间值长波辐射净强度*elsevalue_4=value_3*((t_earth(m)**4-tarea**4)*(1-cos_ph(3))/2) !计算与地面之间值长波辐射净强度*endif*set,long_wave(3,l),value_4/0.03 !计算长波辐射净强度的等效生热率allselvsel,s,,,2 !选择编号为2的体元素eslv,r !选择当前所选体元素上的所有单元nsle,r !选择当前所选单元上的所有节点*get,node_total,node,0,count ! 提取当前所选节点的个数*get,node_num1(1),node,0,num,min ! 提取当前所选节点中节点的最小编号，并赋值给变量node_num1(1)*do,n,2,node_total,1*get,node_num1(n),node,node_num1(n-1),nxth ! 将当前所选择的节点的编号赋值给数组node_num1*enddo*do,f,1,node_total,1*get,node_temp1(f),node,node_num1(f),temp ! 提取当前所选节点的节点温度值，并赋值给数组node_temp1*enddototal=0*do,f,1,node_total,1total=total+node_temp1(f) ! 计算当前所选节点的温度值之和，并赋值给变量total*enddotarea=total/node_total+273 ! 计算当前所选节点温度值的平均值，并赋值给tareatsky=t_outdoor(m)-6+273 !计算当前时刻的天空温度值，并赋值给tsky*if,fuhao(4),eq,100,thenvalue_4=value_3*((tsky**4-tarea**4)*(1+cos_ph(4))/2) !计算与天空之间值长波辐射净强度*elsevalue_4=value_3*((t_earth(m)**4-tarea**4)*(1-cos_ph(4))/2) !计算与地面之间值长波辐射净强度*endif*set,long_wave(4,l),value_4/0.03 !计算长波辐射净强度的等效生热率allsel*enddo*enddo。
ansys命令流语法
ansys命令流语法ANSYS命令流语法是使用ANSYS软件进行仿真分析的关键部分。
它是一种将命令以特定顺序组合在一起的方式，以实现特定的分析目标。
本文将介绍ANSYS命令流语法的基本语法规则和常用命令，以及如何使用它们进行仿真分析。
一、ANSYS命令流语法的基本语法规则1. 命令的基本格式：命令[选项] [参数1, 参数2, ...]2. 命令的执行顺序：ANSYS命令流是按照命令的顺序逐条执行的。
如果需要改变执行顺序，可以使用条件语句、循环语句等控制结构。
3. 注释：可以在命令流中添加注释，以"!"开头。
注释部分不会被执行，可以用于解释命令的用途或添加说明。
4. 变量和参数：可以使用变量和参数来存储和传递数据。
变量以"$"开头，参数以"%"开头。
二、常用命令1. Preprocessor命令：用于定义和准备分析模型的预处理操作。
- /PREP7：进入预处理器界面。
- ET，MP，REAL等：定义单元类型、材料属性、实数等。
- K，L，A等：创建节点、单元、区域等。
2. Solution命令：用于设置和运行分析求解器。
- /SOLU：进入求解器界面。
- SOLVE，ANTYPE等：设置分析类型、求解选项等。
- D，S等：定义边界条件、加载条件等。
3. Postprocessor命令：用于后处理和分析结果的可视化。
- /POST1：进入后处理器界面。
- PLOT，PDEF等：绘制图形、定义图形属性等。
- PRINT，*VWRITE等：输出结果数据。
三、使用ANSYS命令流语法进行仿真分析使用ANSYS命令流语法进行仿真分析的一般步骤如下：1. 导入几何模型：使用CAD软件创建几何模型，并将其导入ANSYS 中。
2. 定义材料属性：根据实际材料的物理特性，使用MP命令定义材料属性。
3. 网格划分：使用网格划分命令划分几何模型，生成有限元网格。
ansys命令流 esel的用法
ansys命令流 esel的用法ESEL是ANSYS中的一个命令流（Command Flow），用于选择特定实体。
它的语法如下：ESEL,TYPE,type_name,OPER,operation_name其中：- TYPE表示实体的类型，可以是NODE（节点）、ELEM （单元）、KP（关键点）等。
- type_name是要选择的实体的名称或编号。
名称可以是节点或单元的标签，编号是其在实体列表中的位置。
- OPER表示选择操作的类型，可以是LOC（通过LOCATION 选择）、NUM（通过编号选择）等。
- operation_name是选择操作的具体名称或参数。
以下是几个常用的ESEL命令流示例：1. 选择单元：ESEL,TYPE,ELEM,REAL,属性名,操作符,阈值例如：ESEL,TYPE,ELEM,REAL,VONM,GT,50表示选择属性“VONM”值大于50的单元。
2. 选择节点：ESEL,TYPE,NODE,LOC,X,操作符,值例如：ESEL,TYPE,NODE,LOC,X,GT,0表示选择X坐标大于0的节点。
3. 选择关键点：ESEL,TYPE,KP,NUM,ID例如：ESEL,TYPE,KP,NUM,1表示选择关键点列表中的第一个关键点。
4. 选择实体列表中的一部分：ESEL,TYPE,ENTITY_LIST,OPR,operation_name,start_num,end_ num,inc_num例如：ESEL,TYPE,ENTITY_LIST,OPR,RANGE,5,10,1表示选择实体列表中第5到第10个实体，包括每个实体。
以上只是ESEL命令流的一些常见用法，实际使用时可以根据具体的需求和ANSYS版本来选择合适的操作。
对于更详细的命令流用法可以参考ANSYS的官方文档或者执行“help esel”命令获取帮助信息。
ANSYS热分析教程及命令流算例
＝＝＝＝＝【热力耦合分析单元简介】＝＝＝＝＝＝SOLID5－三维耦合场实体具有三维磁场、温度场、电场、压电场和结构场之间有限耦合的功能。
本单元由8个节点定义，每个节点有6个自由度。
在静态磁场分析中，可以使用标量势公式（对于简化的RSP，微分的DSP，通用的GSP）。
在结构和压电分析中，具有大变形的应力钢化功能。
与其相似的耦合场单元有PLANE13、SOLID62和SOLID98。
INFIN9－二维无限边界用于模拟一个二维无界问题的开放边界。
具有两个节点，每个节点上带有磁向量势或温度自由度。
所依附的单元类型可以为PLANE13和PLANE53磁单元，或PLANE55和PLANE77和PLANE35热单元。
使用磁自由度（ＡＺ）时，分析可以是线性的也可以是非线性的，静态的或动态的。
使用热自由度时，只能进行线性稳态分析。
PLANE13－二维耦合场实体具有二维磁场、温度场、电场和结构场之间有限耦合的功能。
由4个节点定义，每个节点可达到4个自由度。
具有非线性磁场功能，可用于模拟B－H曲线和永久磁铁去磁曲线。
具有大变形和应力钢化功能。
当用于纯结构分析时，具有大变形功能，相似的耦合场单元有SOLID5、SOLID98和SOLID62。
LINK31－辐射线单元用于模拟空间两点间辐射热流率的单轴单元。
每个节点有一个自由度。
可用于二维（平面或轴对称）或三维的、稳态的或瞬态的热分析问题。
允许形状因子和面积分别乘以温度的经验公式是有效的。
发射率可与温度相关。
如果包含热辐射单元的模型还需要进行结构分析，辐射单元应当被一个等效的或（空）结构单元所代替。
LINK32－二维传导杆用于两节点间热传导的单轴单元。
该单元每个节点只有一个温度自由度。
可用于二维（平面或轴对称）稳态或瞬态的热分析问题。
如果包含热传导杆单元的模型还需进行结构分析，该单元可被一个等效的结构单元所代替。
LINK33－三维传导杆用于节点间热传导的单轴单元。
该单元每个节点只有一个温度自由度。
ANSYS算例4722基于4节点四边形单元的矩形薄板分析GUI及命令流
【ANSYS 算例】4.7.2(2) 基于4节点四边形单元的矩形薄板分析针对【MATLAB 算例】4.7.2(1)的模型，即如图4-21所示的一个薄平板，在右端部受集中力F 作用，其中的参数为：75110Pa,=1/3,=0.1m,=110N E t F μ=⨯⨯。
在ANSYS 平台上，完成相应的力学分析。
(a) 问题描述 (b) 有限元分析模型图4-21 右端部受集中力作用的薄平板解答 下面对该平面结构进行整体建模和分析，按照图4-21(b)的模型，该结构被划分为左右两个4节点单元。
1 基于图形界面(GUI)的交互式操作(step by step)(1) 进入ANSYS(设定工作目录和工作文件)程序 →ANSYS →ANSYS product launcher →Filemanagement →Working directory (设置工作目录) →Jobname (设置工作文件名): rectangular →Run → OK (2) 设置计算类型ANSYS Main Menu : Preferences … → Structural → OK(3) 选择单元类型ANSYS Main Menu : Preprocessor →Element Type →Add/Edit/Delete … →Add … →Solid ：Quad 4node 42 →OK (返回到Element Types 窗口) → Options … →K3: Plane Strs w/thk(带厚度的平面应力问题) →OK →Close(4) 定义材料参数ANSYS Main Menu : Preprocessor →Material Props →Material Models →Structural →Linear →Elastic → Isotropic: EX:1e7 (弹性模量)，PRXY: 0.3 (泊松比) → OK → 鼠标点击该窗口右上角的“ ”来关闭该窗口(5) 定义实常数以确定平面问题的厚度ANSYS Main Menu: Preprocessor →Real Constants … →Add/Edit/Delete →Add →Type 1→ OK →Real Constant Set No: 1 (第1号实常数), THK: 0.1 (平面问题的厚度) →OK →Close(6) 建立节点ANSYS Main Menu: Preprocessor →Modeling → Create → Nodes → In Active CS → NODE:(节点编号) → X,Y, Z:(节点坐标) → THXY ，THYZ ，THZX ：（旋转角度，不填，默认为零）→ Apply → 最后一个节点坐标输入→ OK(7) 建立单元ANSYS Main Menu: Preprocessor →Modeling → Create → Elements → User Numbered → Thru Nodes → Number to assign to element: 1 → Apply → 用鼠标依次选取点3，5，6，4 → Apply → Number to assign to element: 2 → Apply → 用鼠标依次选取点1，3，4，2 → OK (8) 模型施加约束和外载左边两节点加X和Y两方向的约束ANSYS Main Menu: Solution→ Define Loads → Apply →Structural → Displacement→ On Nodes→用鼠标选择左侧边上两个节点(节点5，6) → OK→ Lab2 DOFs to be constrained : All DOF → VALUE → 0(默认值为零)→ OK右边加Y方向的外载ANSYS Main Menu: Solution→ Define Loads → Apply →Structural → Force/Moment→ On Nodes→用鼠标选择右侧边上的两个节点(节点1，2)→ OK → Lab Direction of force/mom: FY ，VALUE：-5.0e4→ OK(9) 分析计算ANSYS Main Menu: Solution → Solve → Current LS→OK→Close (Solution is done! ) →关闭文字窗口(10) 结果显示ANSYS Main Menu: General Postproc → List Results→ Nodel Solution→ DOF Solution →Displacement vector sum→ OK(可得到节点位移结果如下)ANSYS Main Menu: General Postproc → List Results→ Reaction Solu → All items → OK , (可得到5，6节点支反力)ANSYS Main Menu: General Postproc →List Results→Element Solution→Stress →X-component of stress→ OK(得到单元X方向应力结果如下，依照相同方法可以得到其它方向结果)(11) 退出系统ANSYS Utility Menu: File→ Exit …→ Save Everything→OK2 完整的命令流!%%%%%%%%% [ANSYS算例] 4.7.2(2) %%%%% begin %%%%/prep7 !进入前处理et,1,plane42 !定义单元类型(no.1 plane42 )KEYOPT,1,3,3 !设置带厚度的平面应力mp,ex,1,1e7 !设置弹性模量mp,prxy,1,0.33333 !设置泊松比R,1,0.1, ! 定义厚度(0.1)!------定义6个节点N,1,2,1,,,,, ! 节点1,坐标(2,1,0)N,2,2,,,,,, ! 节点2,坐标(2,0,0)N,3,1,1,,,,, ! 节点3,坐标(1,1,0)N,4,1,0,,,,, ! 节点4,坐标(1,0,0)N,5,0,1,,,,, ! 节点5,坐标(0,1,0)N,6,0,0,,,,, ! 节点6,坐标(0,0,0)E,5,6,4,3 ! 由4个节点生成一个单元E,3,4,2,1 ! 由4个节点生成另一个单元d,5,all ! 对5号节点, 完全位移约束d,6,all ! 对6号节点, 完全位移约束F,1,FY,-50000 ! 对1号节点, 施加FY=-50000F,2,FY,-50000 ! 对2号节点, 施加FY=-50000!=====进入求解模块/solu !求解模块solve !求解finish !退出所在模块!=====进入一般的后处理模块/POST1 !进入后处理PLDISP,1 !计算的变形位移显示(变形前与后的对照)finish !退出!%%%%%%%%% [ANSYS算例] 4.7.2(2) %%%%% end %%%%。
ANSYS单元特性之命令流算例
EX1.1 （LINK1）（1）进入后处理模块,显示节点位移和杆件内力MID_NODE = NODE (A/2,-B,0 )! 寻找距离位置（A/2,-B,0）最近的点，存入MID_NODE*GET,DISP,NODE,MID_NODE,U,Y！提取节点MID_NODE上的位移UY，若果已知要求的节点，直接提取即可。
LEFT_EL = ENEARN (MID_NODE)! 需找距离节点MID_NODE最近的单元，存入LEFT_EL ETABLE,STRS,LS,1! 用轴向应力SAXL的编号“LS,1”定义单元表STRS*GET,STRSS,ELEM,LEFT_EL,ETAB,STRS! 从单元表STRS中提取LEFT_EL单元的应力结果，存入变量STRSS。
注意：提取的轴向应力结果具体到指定的单元。
（2）申明数组,提取计算结果,并比较计算误差*DIM,LABEL,CHAR,2！定义2个元素的字符型数组LABEL*DIM,V ALUE,,2,3！定义2*3的数值型数组V ALUELABEL(1) = 'STRS_MPa','DEF_mm' ! 给字符型数组的第1个元素赋值*VFILL,V ALUE(1,1),DATA,1,-0.05498 ! 给其他数值型数组中的元素赋值*VFILL,V ALUE(1,2),DATA,STRSS,DISP*VFILL,V ALUE(1,3),DATA,ABS(STRSS /1 ) ,ABS( DISP /0.05498 )/OUT,EX1_1,out ！将输出内容重定向到文件EX1_1.out/COM ! 以注释形式输出内容/COM,------------------- EX1.1 RESULTS COMPARISON ---------------------/COM,/COM, | TARGET | ANSYS | RATIO/COM,*VWRITE,LABEL(1),V ALUE(1,1),V ALUE(1,2),V ALUE(1,3)(1X,A8,' ',F10.3,' ',F10.3,' ',1F5.3)/COM,----------------------------------------------------------------/OUT ! 结束数据重定向，关闭输出文件FINISH*LIST,EX1_1,out ! 列表显示文件EX1_1.out的内容EX1.2 （LINK1）/PNUM, NODE,1！打开节点编号显示/NUMBER, 2！只显示编号，不使用色彩列表显示节点位移和单元的计算结果PRDISP! 列表显示节点位移值计算结果ETABLE, MFORX,SMISC,1！以杆单元的轴力为内容，建立单元表MFORXETABLE, SAXL, LS, 1 ！以杆的轴向应力为内容，建立单元表SAXLETABLE, EPELAXL, LEPEL, 1! 以杆单元的轴向应变为内容，建立单元表EPELAXLPRETAB! 显示单元表中的计算结果/NUMBER, 0！显示编号，并使用彩色PLETAB, MFORX ！用色度图显示单元表MFORX中杆件轴力图EX1.3 （LINK1）NSEL,S,LOC,Y,1.0 ！选择所有位于Y=1.0位置上的节点FSUM！累计叠加选择集中所有节点上的反力*GET,REAC_1,FSUM,,ITEM,FY ！将累加结果中的FY（Y方向的力）保存到变量REAC_1中EX1.4 （LINK1）R,1,65e-6！定义第1类实常数，杆件截面面积为65mm^2，在转化为国际单位制时操作TREF,70 ! 设定参考温度为70度BFUNIF,TEMP,80 ! 温度从原来的70度均匀上升到80度（TREF+10）EX1.5 （PLANE42 AND CONTAC26）ETABLE,STRSX,S,X！定义X方向的应力为单元表STRSX*GET,STRSSX,ELEM,3,ETAB,STRSX！从单元表STRSX中提取3号单元的X向应力，存入STRSSX。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
EX1.1 （LINK1）（1）进入后处理模块,显示节点位移和杆件内力MID_NODE = NODE (A/2,-B,0 )! 寻找距离位置（A/2,-B,0）最近的点，存入MID_NODE*GET,DISP,NODE,MID_NODE,U,Y！提取节点MID_NODE上的位移UY，若果已知要求的节点，直接提取即可。
LEFT_EL = ENEARN (MID_NODE)! 需找距离节点MID_NODE最近的单元，存入LEFT_EL ETABLE,STRS,LS,1! 用轴向应力SAXL的编号“LS,1”定义单元表STRS*GET,STRSS,ELEM,LEFT_EL,ETAB,STRS! 从单元表STRS中提取LEFT_EL单元的应力结果，存入变量STRSS。
注意：提取的轴向应力结果具体到指定的单元。
（2）申明数组,提取计算结果,并比较计算误差*DIM,LABEL,CHAR,2！定义2个元素的字符型数组LABEL*DIM,V ALUE,,2,3！定义2*3的数值型数组V ALUELABEL(1) = 'STRS_MPa','DEF_mm' ! 给字符型数组的第1个元素赋值*VFILL,V ALUE(1,1),DATA,1,-0.05498 ! 给其他数值型数组中的元素赋值*VFILL,V ALUE(1,2),DATA,STRSS,DISP*VFILL,V ALUE(1,3),DATA,ABS(STRSS /1 ) ,ABS( DISP /0.05498 )/OUT,EX1_1,out ！将输出内容重定向到文件EX1_1.out/COM ! 以注释形式输出内容/COM,------------------- EX1.1 RESULTS COMPARISON ---------------------/COM,/COM, | TARGET | ANSYS | RATIO/COM,*VWRITE,LABEL(1),V ALUE(1,1),V ALUE(1,2),V ALUE(1,3)(1X,A8,' ',F10.3,' ',F10.3,' ',1F5.3)/COM,----------------------------------------------------------------/OUT ! 结束数据重定向，关闭输出文件FINISH*LIST,EX1_1,out ! 列表显示文件EX1_1.out的内容EX1.2 （LINK1）/PNUM, NODE,1！打开节点编号显示/NUMBER, 2！只显示编号，不使用色彩列表显示节点位移和单元的计算结果PRDISP! 列表显示节点位移值计算结果ETABLE, MFORX,SMISC,1！以杆单元的轴力为内容，建立单元表MFORXETABLE, SAXL, LS, 1 ！以杆的轴向应力为内容，建立单元表SAXLETABLE, EPELAXL, LEPEL, 1! 以杆单元的轴向应变为内容，建立单元表EPELAXLPRETAB! 显示单元表中的计算结果/NUMBER, 0！显示编号，并使用彩色PLETAB, MFORX ！用色度图显示单元表MFORX中杆件轴力图EX1.3 （LINK1）NSEL,S,LOC,Y,1.0 ！选择所有位于Y=1.0位置上的节点FSUM！累计叠加选择集中所有节点上的反力*GET,REAC_1,FSUM,,ITEM,FY ！将累加结果中的FY（Y方向的力）保存到变量REAC_1中EX1.4 （LINK1）R,1,65e-6！定义第1类实常数，杆件截面面积为65mm^2，在转化为国际单位制时操作TREF,70 ! 设定参考温度为70度BFUNIF,TEMP,80 ! 温度从原来的70度均匀上升到80度（TREF+10）EX1.5 （PLANE42 AND CONTAC26）ETABLE,STRSX,S,X！定义X方向的应力为单元表STRSX*GET,STRSSX,ELEM,3,ETAB,STRSX！从单元表STRSX中提取3号单元的X向应力，存入STRSSX。
注意：已知某一具体单元，在定义完单元表后，就可以直接用*GET语句提取这一具体单元的表中定义的项目。
EX2.1 （BEAM3）NPLOT! 显示节点位置图形，但是不显示节点号码NPLOT, 1！显示节点位置图形，同时显示节点号码NLIST！列表显示节点在直角坐标系下的坐标值DSYS, 1 ！改变显示坐标系，列出资料时，转换到圆柱坐标系NLIST！列表显示节点在圆柱坐标系下的坐标值定义以两端弯矩和剪力为内容的单元表，并列出单元的单元表数据ETABLE,IMOMENT, SMISC, 6！以单元I节点弯矩为内容，定义单元表ETABLE,JMOMENT, SMISC, 12！以单元J节点弯矩为内容，定义单元表ETABLE, ISHEAR, SMISC, 2！以单元I节点剪力为内容，定义单元表ETABLE, JSHEAR, SMISC, 8 ！以单元J节点剪力为内容，定义单元表PRETAB！列表显示单元表中单元的计算结果用单元表数据绘制剪力图/TITLE, SHEAR FORCE DIAGRAM! 定义剪力图窗口标题PLLS, ISHEAR, JSHEAR！结构剪力分布图用单元表数据绘制弯矩图/TITLE, BENDING MOMENT DIAGRAM！定义弯矩图窗口标题PLLS, IMOMENT, JMOMENT！结构弯矩分布图EX2.2 （BEAM3）E,1,2！过1,2节点定义第1个单元EGEN,4,1,1 ！按照前面的单元模式，生成4次，每次两端节点编号增加1，生成其余3个单元。
SET,1,1 ！读入第1荷载步第1子步的计算结果PRNSOL,U,COMP！列表显示节点线位移和阶巅峰线位移的矢量和PRNSOL,ROT,COMP！列表显示节点转角和节点角位移的矢量和EX2.3 （BEAM4）ET,1,BEAM4,,,,,,1 ！定义单元类型为BEAM4，并控制在单元坐标系中输出弯矩NSEL,S,,,5！选择编号不小于5的所有节点DSYM,SYMM,X！定义刚选定的所有节点关于X轴的对称位移边界条件SFBEAM,ALL,1,PRES,314！在所有梁单元表面施加均布荷载进入后处理模块，显示位移计算结果，提取左边支座节点的转角和跨中挠度/POST1NSEL,S,,,1,5,4！选中1号和5号节点PRNSOL,U,Z！列表显示选择集中所有节点的Z向位移PRNSOL,ROT,Y ！列表显示选择集中的所有节点的绕Y轴的转角ROTYNSEL,ALL！重新选中所有节点PRRSOL！列表显示被约束的节点的支反力（Print Reaction Solution）RGHT_END = NODE (2.54,0,0)！选择距离(2.54,0,0)最近的节点，并存入RGHT_ENDLFT_END = NODE (0,0,0)！选择距离(0,0,0)最近的节点，并存入LFT_END*GET,UZ_MX_C2,NODE,RGHT_END,U,Z！获取节点RGHT_END的线位移UZ*GET,SLOPE_C2,NODE,LFT_END,ROT,Y ！获取节点LFT_END的转角ROTYFINISH在时间历程后处理器中，叠加处理5号节点的弯矩/POST26RFORCE,2,RGHT_END,M,Y ！叠加RGHT_END节点的弯矩MY，存入第2个POST26变量STORE！保存变量*GET,M_MX_C2,V ARI,2,EXTREM,VMAX ！获取第2个POST26变量的最大值FINISHEX2.4 （SOLID5，SOLID95，SOLID98，SOLID147）smrt,off ! 关闭智能化功能LESIZE,ALL,,,10! 定义所有线的划分段数为10段ESIZE,25.4! 定义默认的单元尺寸为25.4V,1,2,4,3,5,6,8,7! 用关键点连接成体VMESH,1 ! 执行体单元划分操作FK,2,FX,4.4483/4! 在2号关键点上施加沿着X方向的，大小为4.4483/4的轴向力*REPEAT,4,2 ! 将前面的命令执行4次，每次操作对象的编号增加2，即完成对2、4、6关键点的荷载定义FKDELE,ALL,FY ! 删除所有关键点上的轴向力定义宏程序，完成自由端位移计算结果的处理*CREATE,MAC! 定义宏程序计算这些结果SET,ARG1,1 ! 定义参数ARG1=1LCDEF,ARG1,ARG1! 基于计算结果定义荷载工况NSEL,S,LOC,X,152.4! 选择节点PRNSOL,U,COMP! 打印位移结果及其节点位移矢量和LCFACT,ARG1,ARG2! 将缩放因子“AGR2”应用给荷载工况1LCASE,ARG1! 将工况1读入到数据库PRNSOL,U,COMP! 打印位移结果及其节点位移矢量和*END! 结束宏程序的定义使用上面定义的宏程序处理所有的3个工况/COM*** !*USE,MAC,1,(1/7.6E-4) ! 使用宏程序MAC，使用参数1，(1/7.6E-4)SET,1,1 ! 读入第1荷载步的计算结果LCSEL,S,0,0 ! 选择荷载步的子集N1=NODE(152.4 , 0, 0) ! 选择距离(152.4 , 0, 0)位置最近的节点，命名为N1*GET,UX1,NODE,N1,U,X ! 获取节点N1的位移UX处理第2个荷载步的计算结果*USE,MAC,2,(1/2.4712) ! 使用宏MAC，使用参数2，(1/2.4712)SET,2,1 !LCSEL,S,0,0 !*GET,UY1,NODE,N1,U,Y !处理第3个荷载步的计算结果*USE,MAC,3,(1/10.9646) !SET,3,1 !LCSEL,S,0,0 !*GET,UZ1,NODE,N1,U,Z !定义字符串变量和相应的理论计算结果*DIM,LABEL,CHAR,3,2 !*DIM,V ALUE,,3,3 !LABEL(1,1) = 'DEFL X ','DEFL Y ','DEFL Z '!LABEL(1,2) = 'mm','mm','mm' !*VFILL,V ALUE(1,1),DATA,7.6E-4,2.471,10.965 !*VFILL,V ALUE(1,2),DATA,UX1,UY1,UZ1!*VFILL,V ALUE(1,3),DATA,ABS(UX1/(7.6E-4)),ABS(UY1/2.471),ABS(UZ1/10.965) SA VE, TABLE_1 ! 保存数据到数据库TABLE_1FINISHNSEL,S,LOC,X,152.4NSEL,R,LOC,Y,2.54NSEL,R,LOC,Z,1.27*GET,MIDD,NDMX ! 获取当前最大节点号，并存入MIDD从第1个表中恢复数据，绘制单元应力云图RESUME,TABLE_1 ! 从数据库中恢复数据/COM *** CLIPPED AND CAPPED DISPLAY OF STRESS CONTOURS ***NSEL,ALL ! 选择所有节点/VIEW,1,2,1,1 ! 设置观察方向为从（2,1,1）指向原点EPLOT ! 绘制单元图/TYPE,1,CAP ! 显示梁的剖面图/DIST,1,5.08 ! 设置观察距离为5.08mm/FOCUS,1,.3,.15,.09 ! 设置截开截面的位置/Auto ! 以最佳比例显示PLNSOL,S,X ! 显示X方向的应力删除宏文件，临时表文件TABLE_*/DELETE,MAC ! 删除宏文件/DELETE,TABLE_1 ! 删除临时性数据库EX2.5 （PIPE16）F,2,MZ,1112*914! 定义作用在节点2上的绕Z轴大小为M=FD=1112*914N.MM的力偶F,2,FX,-1112 ! 定义作用在2号节点的沿X轴大小为1112N的力EX2.6 （SHELL63，BEAM4）定义耦合自由度，以保证变形对称CP,1,UZ,2,12 ! 定义2号节点和12号节点的UZ自由度为耦合自由度，定义为集合1CP,2,ROTY,2,12! 定义2号节点和12号节点的ROTY自由度为耦合自由度，定义为集合1 CPSGEN,6,1,1,2! 按照前面定义的耦合自由度集合1到2，循环6次生成其余的耦合自由度ETABLE,STRS,S,1 ！用壳单元SHELL63的最大拉应力S1为内容定义单元表STRS ESORT,STRS！按最大拉应力对单元表STRS的排序*GET,SMAX1,SORT,,MAX！获取排序后的最大拉应力值，存入变量SMAX1PRNSOL,DOF ！列表显示节点位移计算结果*DO,I,1,10 ！循环10次定义不同位置的单元不同的截面几何参数R,I,1,1,(I-1)*76.2*12.7*12.7*12.7/120,6.35,1 ！REAL,1,AREA1,IZ1，IY1，TKZB1，TKYB1,IX1 RMORE,1,1,I*76.2*12.7*12.7*12.7/120,6.35,1！AREA2,IZ2，IY2，TKZB2，TKYB2,IX2*ENDDORMODIF,1,3,76.2*12.7*12.7*12.7/120/2 ! 修改实常数RLIST！列出实常数/NOPR ！关闭输入显示RESUME,TABLE_2 从数据库中恢复数据/GOPR！打开输入显示EX2.7 （BEAM54）*GET,UY,NODE,1,U,Y ! 提取1号节点的Y方向的线位移UY，并保存到变量UY*status,parm！显示内存变量的内容EX2.8 （BEAM54）CNVTOL,F,,1E-4！设置荷载力的收敛误差为1E-4NCNVTOL,M,,1E-1! 设置集中力偶的收敛误差为0.1N.MEND_NODE = NODE (0,3048/2,0)! 选择距离(0,3048/2,0)最近的节点，命名为END_NODE*GET,DEF,NODE,END_NODE,U,X！提取节点END_NODE的位移UX，存入变量DEFETABLE,TENS,NMISC,1 ！以BEAM54的最大拉应力（NMISC,1）为内容定义单元表TENS ETABLE,COMP,NMISC,2 ！以BEAM54的最大拉应力（NMISC,2）为内容定义单元表COMP*GET,STS_TENS,ELEM,1,ETAB,TENS ！提取1号单元的最大拉应力，存入变量STS_TENS*GET,STS_COMP,ELEM,1,ETAB,COMP ！提取1号单元的最大压应力，存入变量STS_COMPEX3.1 （BEAM3）*GET,FCR,MODE,1,FREQ！提取第1阶频率计算结果，并保存到变量FCR.。