移动体热源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命令流大全（整理）
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.建模尺寸为5mm×5mm×0.5mm的SiCpAl材料（此处特别注意坐标系的位置，建议用DM建模）2.网格的划分实现移动热源路径网格划分密集，其他部分网格划分尺寸正常，方便计算，提高效率。
下图为网格划分示意图及网格划分设置。
3.分析计算（1）入射激光功率密度分布为高斯分布：I(r)=(1−R)2Pπω2exp(−2r2ω2)，（1）式中R为反射率，P为激光输出功率，ω为入射激光光斑半径，r为横向柱坐标，r=0处为光斑中心。
对模型施加负载条件时，将高斯热流密度用ANSYS APDL编写命令流加载到SiCp/Al上表面。
已知条件：功率P，行走速度v=5mm/s，光斑半径ω=0.05mm，时间t=1s。
则有，高斯移动热源：Q=2Pπ0.000052⋅exp(−2×(x−0.0025)2+(y−0.005×t)20.000052)×90%步骤：①首先在Mechanical APDL经典界面施加创建高斯热源函数的命令流；将上述公式带入，如下图所示。
②完成好函数输入之后，保存函数；然后读入刚刚保存的函数，命名为HFLUX，如下图所示：③提取命令流，命令流在文末（2）选择要加载的面，重命名为A1（3）载荷步设置在Step Controls（步长控制）中，Step End Time（结束时间）设置为1s，打开Auto Time Stepping(自动划分时间)，在Define By中设置为Substeps（子步），设置为100子步。
如下图所示。
（4）边界设置取模型下表面为绝缘边界，其他表面为对流换热面，对流换热系数取10W/(m3·K)，设置实验环境温度为22℃。
4.温度场仿真结果及分析命令流*DEL,_FNCNAME*DEL,_FNCMTID*DEL,_FNCCSYS*SET,_FNCNAME,'HFLUX'*SET,_FNCCSYS,0! /INPUT,1.func,,,1*DIM,%_FNCNAME%,TABLE,6,24,1,,,,%_FNCCSYS%!! Begin of equation: 4e8*exp(-2*(({X}-0.0025)^2+({Y}-0.005* ! {TIME})^2)/0.00005^2)*SET,%_FNCNAME%(0,0,1), 0.0, -999*SET,%_FNCNAME%(2,0,1), 0.0*SET,%_FNCNAME%(3,0,1), 0.0*SET,%_FNCNAME%(4,0,1), 0.0*SET,%_FNCNAME%(5,0,1), 0.0*SET,%_FNCNAME%(6,0,1), 0.0*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0, 0, 0, 0*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 0, 0, -1*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 2, -2*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 2, 0, 0, -3*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 3, -1*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0.0025, 0, 0, 2*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 2, 2, -1*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 2, 0, 0, -3*SET,%_FNCNAME%(0,9,1), 0.0, -4, 0, 1, -3, 17, -1*SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 0.005, 0, 0, 1*SET,%_FNCNAME%(0,11,1), 0.0, -3, 0, 1, -1, 3, 1*SET,%_FNCNAME%(0,12,1), 0.0, -1, 0, 1, 3, 2, -3*SET,%_FNCNAME%(0,13,1), 0.0, -3, 0, 2, 0, 0, -1*SET,%_FNCNAME%(0,14,1), 0.0, -5, 0, 1, -1, 17, -3*SET,%_FNCNAME%(0,15,1), 0.0, -1, 0, 1, -4, 1, -5*SET,%_FNCNAME%(0,16,1), 0.0, -3, 0, 1, -2, 3, -1*SET,%_FNCNAME%(0,17,1), 0.0, -1, 0, 0.00005, 0, 0, 0*SET,%_FNCNAME%(0,18,1), 0.0, -2, 0, 2, 0, 0, -1*SET,%_FNCNAME%(0,19,1), 0.0, -4, 0, 1, -1, 17, -2*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 1, -3, 4, -4*SET,%_FNCNAME%(0,21,1), 0.0, -1, 7, 1, -1, 0, 0*SET,%_FNCNAME%(0,22,1), 0.0, -2, 0, 4e8, 0, 0, -1*SET,%_FNCNAME%(0,23,1), 0.0, -3, 0, 1, -2, 3, -1*SET,%_FNCNAME%(0,24,1), 0.0, 99, 0, 1, -3, 0, 0! End of equation: 4e8*exp(-2*(({X}-0.0025)^2+({Y}-0.005*{TIME})^2)/0.00005^2) !-->SF,A1,HFLUX,%HFLUX%。
(完整版)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命令流大全
《史上最全的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定”义单元表STRS ETABLE, MFORX,SMISC,1 ! 以杆单元的轴力为内容, 建立单元表MFORX ETABLE, SAXL, LS, 1 ! 以杆单元的轴向应力为内容, 建立单元表SAXLETABLE, EPELAXL, LEPEL, 1 ! 以杆单元的轴向应变为内容, 建立单元表EPELAXLETABLE,STRS_ST,LS,1 ! 以杆件的轴向应力“LS,1为”内容定义单元表STRS_STETABLE, STRS_CO, LS,1 ! 以杆件的轴向应力“LS,1定”义单元表STRS_CO ETABLE,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,STR”S_CO 从单元表STRS_CO 中提取COPPER_E 单元的应力结果，存入变量STRSS_CO10 FINISH ! 退出以前的模块11, /CLEAR, START ! 清除系统中的所有数据,读入启动文件的设置12 /UNITS, SI ! 申明采用国际单位制14 /NUMBER, 2 ! 只显示编号, 不使用彩色/NUMBER, 0 ! 显示编号, 并使用彩色15 /SOLU ! 进入求解模块:定义力和位移边界条件,并求解ANTYPE, STATIC ! 申明分析类型是静力分析（STATIC 或者0） OUTPR, BASIC, ALL ! 在输出结果中, 列出所有荷载步的基本计算结果OUTPR,BASIC,ALL ! 指定输出所有节点的基本数据OUTPR,BASIC,LAST ! 选择基本输出选项,直到最后一个荷载步OUTPR,,1 ! 输出第 1 个荷载步的基本计算结果OUTPR,BASIC,1 ! 选择第 1 荷载步的基本输出项目OUTPR,NLOAD,1 ! 指定输出第 1 荷载步的内容OUTRES,ALL,0 ! 设置将所有数据不记录到数据库。
史上最全的ANSYS命令流大全
《史上最全的ANSYS命令流查询与解释》[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 !设置将所有数据不记录到数据库。
move source
ANSYS Workbench移动热源施加本篇博文主要介绍如何在ANSYS WORKBENCH里面如何施加移动热源，本人也是通过借鉴网上资料、论坛和请教交流，做出的一个移动热源初级实例。
1.问题描述如下图所示，尺寸为0.1x0.1x0.005m长方体，在长方体中间沿着Y方向施加一个移动热源，热源的速度为0.1m/s，热源为热流密度，值为时间位移函数，如下图所示。
2.分析思路（1）首先在APDL经典界面施加创建高斯热源函数的命令流；（2）在WB中创建瞬态分析模块，创建有限元模型；（3）将APDL命令流插入到WB中；（4）计算求解查看后处理。
3.步骤（1）创建高斯热源函数命令流打开ANSYS经典界面，在函数编辑器下创建如下函数：4e7*exp(-3*(({X}-0.05)^2+({Y}-0.01*{TIME})^2)/0.005^2) 如下图所示：完成好函数输入之后，保存函数；然后读入刚刚保存的函数，命名为HFLUX，如下图所示：到此，高斯热源函数即完成创建，只需要将以上操作的命令流提取出来即可，命令流件文章末尾。
（2）在WB中创建瞬态热分析模块，创建几何模型、材料属性和划分网格，注意中间的网格要细化，如下图所示：在几何体上表面创建一个Named Selection，命名为A1，如下图所示：求解设置，设置仿真时间为10s，子步为50，如下图所示；创建对流换热，选择除上表面之外的其余5个面。
（3）插入命令流。
将（1）中的命令流粘贴到WB中，最后再加上如下语句：SF,A1,HFLUX, %HFLUX%（4）进行求解，求解完成之后，查看温度变化，如下图所示：4.总结对于WB热源的分析，本文只为提供一种基于WB平台下移动热源施加的方法，本人也是刚学习，整个过程肯定存在不足，请大家积极留言，给出好的建议和意见，一起学习。
这在有限元分析中是一个非常有意义的课题，希望大家给出宝贵意见，一起对该问题进行深一步的研究，谢谢！附件：命令流*DEL,_FNCNAME*DEL,_FNCMTID*DEL,_FNCCSYS*SET,_FNCNAME,'HFLUX'*SET,_FNCCSYS,0! /INPUT,HFLUX.func,,,1*DIM,%_FNCNAME%,TABLE,6,24,1,,,,%_FNCCSYS%! Begin of equation:4e7*exp(-3*(({X}-0.05)^2+({Y}-0.01*{TIME})^2)/0.005^2) *SET,%_FNCNAME%(0,0,1), 0.0, -999*SET,%_FNCNAME%(2,0,1), 0.0*SET,%_FNCNAME%(3,0,1), 0.0*SET,%_FNCNAME%(4,0,1), 0.0*SET,%_FNCNAME%(5,0,1), 0.0*SET,%_FNCNAME%(6,0,1), 0.0*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0, 0, 0, 0*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 0, 0, -1*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 2, -2 *SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 3, 0, 0, -3*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 3, -1*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0.05, 0, 0, 2 *SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 2, 2, -1*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 2, 0, 0, -3*SET,%_FNCNAME%(0,9,1), 0.0, -4, 0, 1, -3, 17, -1 *SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 0.01, 0, 0, 1 *SET,%_FNCNAME%(0,11,1), 0.0, -3, 0, 1, -1, 3, 1*SET,%_FNCNAME%(0,12,1), 0.0, -1, 0, 1, 3, 2, -3*SET,%_FNCNAME%(0,13,1), 0.0, -3, 0, 2, 0, 0, -1*SET,%_FNCNAME%(0,14,1), 0.0, -5, 0, 1, -1, 17, -3 *SET,%_FNCNAME%(0,15,1), 0.0, -1, 0, 1, -4, 1, -5 *SET,%_FNCNAME%(0,16,1), 0.0, -3, 0, 1, -2, 3, -1 *SET,%_FNCNAME%(0,17,1), 0.0, -1, 0, 0.005, 0, 0, 0 *SET,%_FNCNAME%(0,18,1), 0.0, -2, 0, 2, 0, 0, -1*SET,%_FNCNAME%(0,19,1), 0.0, -4, 0, 1, -1, 17, -2 *SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 1, -3, 4, -4 *SET,%_FNCNAME%(0,21,1), 0.0, -1, 7, 1, -1, 0, 0*SET,%_FNCNAME%(0,22,1), 0.0, -2, 0, 4e7, 0, 0, -1 *SET,%_FNCNAME%(0,23,1), 0.0, -3, 0, 1, -2, 3, -1 *SET,%_FNCNAME%(0,24,1), 0.0, 99, 0, 1, -3, 0, 0! End of equation:4e7*exp(-3*(({X}-0.05)^2+({Y}-0.01*{TIME})^2)/0.005^2) SF,a1,HFLUX, %HFLUX%。
(完整版)史上最全的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结构解析单元功能与特征/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命令流解释大全
一、定义材料号及特性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=，C=27，ψ=45的命令如下：MP，EX，1，1E8MP，NUXY，1，TB，DP，1TBDATA，1，27TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,deg三、单元生死载荷步!第一个载荷步TIME,... !设定时间值（静力分析选项）NLGEOM,ON !打开大位移效果NROPT,FULL !设定牛顿－拉夫森选项ESTIF,... !设定非缺省缩减因子（可选）ESEL,... !选择在本载荷步中将不激活的单元EKILL,... !不激活选择的单元ESEL,S,LIVE !选择所有活动单元NSLE,S !选择所有活动结点NSEL,INVE !选择所有非活动结点（不与活动单元相连的结点）D,ALL,ALL,0 !约束所有不活动的结点自由度（可选）NSEL,ALL !选择所有结点ESEL,ALL !选择所有单元D,... !施加合适的约束F,... !施加合适的活动结点自由度载荷SF,... !施加合适的单元载荷BF,... !施加合适的体载荷SAVESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D,F,SF和BF命令得到更详细的解释。
(完整版)ANSYS最常用命令流+中文注释(超级大全)
ANSYS最常用命令流+中文注释VSBV, 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,,keep,,则为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，MA T，……MP，NUXY，MAT，……定义DP材料单元表（这里不考虑温度）：TB，DP，MA T进入单元表并编辑添加单元表：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, KSWP Type，是选择的方式，有选择（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,... !施加合适的体载荷SA VESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D, F,SF和BF命令得到更详细的解释。
ANSYS命令流使用方法(中文)修改
Finish(退出四大模块，回到BEGIN层)/clear (清空内存，开始新的计算)1．定义参数、数组，并赋值.2．/prep7(进入前处理)定义几何图形：关键点、线、面、体定义几个所关心的节点，以备后处理时调用节点号。
设材料线弹性、非线性特性设置单元类型及相应KEYOPT设置实常数设置网格划分，划分网格根据需要耦合某些节点自由度定义单元表3．/solu加边界条件设置求解选项定义载荷步求解载荷步4./post1（通用后处理）5./post26 （时间历程后处理）6.PLOTCONTROL菜单命令7.参数化设计语言8.理论手册Finish(退出四大模块，回到BEGIN层)/clear (清空内存，开始新的计算)1.定义参数、数组，并赋值.dim, par, type, imax, jmax, kmax, var1, vae2, var3 定义数组par: 数组名type：array 数组，如同fortran,下标最小号为1，可以多达三维（缺省）char 字符串组（每个元素最多8个字符）tableimax,jmax, kmax 各维的最大下标号var1,var2,var3 各维变量名，缺省为row,column,plane(当type为table时) 2./prep7(进入前处理)2.1 设置单元类型及相应KEYOPTET, itype, ename, kop1……kop6, inopr 设定当前单元类型Itype：单元号Ename：单元名设置实常数Keyopt, itype, knum, valueitype: 已定义的单元类型号knum: 单元的关键字号value: 数值注意：如果,则必须使用keyopt命令，否则也可在ET命令中输入2.2 定义几个所关心的节点，以备后处理时调用节点号。
n,node,x,y,z,thxy, thyz, thzx 根据坐标定义节点号如果已有此节点，则原节点被重新定义，一般为最大节点号。
ansys命令流详解超全)
2． /prep7(进入前处理)
定义几何图形：关键点、线、面、体
定义几个所关心的节点，以备后处理时调用节点号。
设材料线弹性、非线性特性
设置单元类型及相应KEYOPT
设置实常数
设置网格划分，划分网格
根据需要耦合某些节点自由度
定义单元表
存盘
3．/solu
加边界条件
设置求解选项
u vdele, nv1, nv2, ninc, kswp 删除体
kswp: 0 只删除体
1 删除体及面、关键点（非公用）
u vgen, itime, nv1, nv2, ninc, dx, dy, dz, kinc, noelem, imove 移动或拷贝体
itime: 份数
nv1, nv2, ninc：拷贝对象编号
ex: 弹性模量
nuxy: 小泊松比
alpx: 热膨胀系数
reft: 参考温度
reft: 参考温度
prxy: 主泊松比
gxy: 剪切模量
mu: 摩擦系数
dens: 质量密度
mat: 材料编号（缺省为当前材料号）
c 材料特性值，或材料之特性，温度曲线中的常数项
c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数
,例：vsel,s,volu,,14
vsel,a,volu,,17,23,2
上面的命令选中了实体编号为 14，17，19，21，23的五个实体
u rforce, nvar, node, item, comp, name 指定待存储的节点力数据
nvar: 变量号
node: 节点号
item comp
F x, y.z
移动体热源ANSYS命令流
!非线性材料求解成功GAUSS分布Z向扫描!!考虑了和基板对流情况下加载热流密度的热分析rb=0.15mm!粉床7.8X3.3mm，烧结区域6X1.5mm!考虑粉床表面的辐射、对流!考虑粉末状、熔化状（液态）、凝固态的边界条件变化!考虑金属蒸发会带走一部分能量，故效率减半eff＝0.7!r=0.2mmFINISH/CLEAR,START/BATCH !设置程序模式为“批处理”/CWD,'d:\Ansyswork\EBM FEM Simulation for Metals\3D Temperture distribution' !更改当前工作目录/FILNAME,EBM_db,0/TITLE,3D T emperture distrbution!选择分析问题的类型：热/NOPR/PMETH,OFF,0KEYW,PR_SET,1KEYW,PR_STRUC,0KEYW,PR_THERM,1KEYW,PR_FLUID,0KEYW,PR_ELMAG,0KEYW,MAGNOD,0KEYW,MAGEDG,0KEYW,MAGHFE,0KEYW,MAGELC,0KEYW,PR_MULTI,0KEYW,PR_CFD,0/GO/UNITS,SI !Standard Units!************************************************************************************** ***********!!开始，进入前处理器，定义下列属性－!0参数定义!1定义单元类型!2实常数!3材料性能!4几何模型!5网格划分控制!6其他!************************************************************************************** ***********!/PREP7!0参数定义!模型参数!几何模型尺寸me-3=mm*SET,Block_x,7.8e-3*SET,Block_x1,0.9e-3*SET,Block_x2,6.9e-3*SET,Block_y,3.3e-3*SET,Block_y1,0.9e-3*SET,Block_y2,2.4e-3*SET,Block_z,0.3e-3*SET,Block_z1,0.9e-3!划分网格的大小*SET,LMeshx,3e-4 !line1*SET,LMeshy,3e-4 !line4*SET,LMeshz,3e-4 !line9!边界条件参数*SET,Tinit,30 !初始温度!1定义单元类型ET,1,SOLID70ET,2,surf152 !表面效应单元r,2KEYOPT,2,4,1 !表面效应单元设置KEYOPT,2,8,4ET,3,SOLID87!2实常数!Nothing!3材料性能MPTEMPMPTEMP,1,25,400,1200,1600,2000,2500MPDATA,DENS,1,1,3983,4688.4,5406.8,5948.8,6908,6331 !密度MPDATA,KXX,1,1,1.6278,2.3812,3.6396,29.8,30.6,32 !热导率MPDATA,c,1,1,472,535,679,859,945,1034 !比热MPTEMP,,,,,,,,MPTEMP,1,0,200,400,800,1200MPTEMP,6,1400,1454,1800,2000MPDATA,ENTH,1,1,0,7.9729e+8,1.6400e+9,3.4621e+9,5.3897e+9 MPDATA,ENTH,1,6,6.3415e+9,8.2595e+9,9.8798e+9,1.0816e+10 !热焓MPTEMP,,,,,,,,MPTEMP,1,20,300,600,900,1200,1500,1800,2000,3000,4000 MPDATA,HF,1,1,6,50,120,200,250,378,700,850,1746,3799 !对流系数MPTEMP,,,,,,,,MPTEMP,1,20,200,400,600,800,900,1000,1100,1200,1420,1460MPDATA,DENS,2,1,7966,7893,7814,7724,7630,7583,7535,7486,7436,7320,7320MPDATA,KXX,2,1,13.31,16.33,19.47,22.38,25.07,26.33,27.53,28.67,29.76,31.95,32 !热导率mpdata,c,2,1,470,508,550,592,634,655,676,698,719,765,765MPTEMP,,,,,,,,MPTEMP,1,0,200,400,800,1200MPTEMP,6,1400,1454,1800,2000MPDATA,ENTH,2,1,0,7.9729e+8,1.6400e+9,3.4621e+9,5.3897e+9MPDATA,ENTH,2,6,6.3415e+9,8.2595e+9,9.8798e+9,1.0816e+10!4几何模型BLOCK,0,Block_x,0,Block_y,0,Block_z,BLOCK,Block_x1,Block_x2,Block_y1,Block_y2,0,Block_z,vsel,allVOVLAP,allBLOCK,0,Block_x,0,Block_y,0,-Block_z1,vglue,allnumcmp,all!划分网格lesize,8,LMeshx lesize,5,LMeshy lesize,13,LMeshz vsel,s,,,1 type,1mat,1 mshkey,1 vmesh,allallselesize,0.3e-3 mshkey,0 mshape,1,3d vsel,s,,,2 type,3mat,2 vmesh,allallselesize,0.3e-3 mshkey,0 mshape,1,3dvsel,s,,,3type,3mat,1vmesh,allallsel!成形件轨迹单元质心排序width=Block_y2-Block_y1!wpoffs,Block_x1,Block_y1,0!local,11,1,Block_x1,Block_y1,0 ALLSEL,ALL*GET,EMAX,ELEM,,COUNT,MAXVSEL,S, , , 1ESLV,S*GET,AEMAX,ELEM,,COUNT1,MAX *DIM,ANE,ARRAY,AEMAX,*DIM,ANE1,ARRAY,AEMAX,*DIM,ANE2,ARRAY,AEMAX,*DIM,ANE3,ARRAY,AEMAX,*DIM,ANE4,ARRAY,AEMAX,*DIM,ANE5,ARRAY,AEMAX,*DIM,ANEX,ARRAY,AEMAX*DIM,ANEY,ARRAY,AEMAX*DIM,ANEZ,ARRAY,AEMAX*DIM,ANEORDER,ARRAY,AEMAX,,1!形心排序VSEL,S, , , 1ESLV,S*GET,ANSEL,ELEM,,COUNT1II=0*DO,I,1,AEMAX*IF,ESEL(I),EQ,1,THENII=II+1ANE(II)=I*ENDIF*ENDDO*DO,I,1,ANSEL*GET,ANEZ(I),ELEM,ANE(I),CENT,Z*GET,ANEY(I),ELEM,ANE(I),CENT,Y*GET,ANEX(I),ELEM,ANE(I),CENT,X *ENDDOe=0*do,i,1,AEMAX*if,ANEY(I),ge,4*width/5+Block_y1,then e=e+1ANE5(e)=ANE(i)*endif*enddod=0*do,i,1,AEMAX*if,ANEY(I),ge,3*width/5+Block_y1,then *if,ANEY(I),lt,4*width/5+Block_y1,then d=d+1ANE4(d)=ANE(i)*endif*endif*enddoc=0*do,i,1,AEMAX*if,ANEY(I),ge,2*width/5+Block_y1,then *if,ANEY(I),lt,3*width/5+Block_y1,then c=c+1ANE3(c)=ANE(i)*endif*endif*enddob=0*do,i,1,AEMAX*if,ANEY(I),ge,width/5+Block_y1,then *if,ANEY(I),lt,2*width/5+Block_y1,then b=b+1ANE2(b)=ANE(i)*endif*endif*enddoa=0*do,i,1,AEMAX*if,ANEY(I),lt,width/5+Block_y1,thena=a+1ANE1(a)=ANE(i)*endif*enddo!*************************对1道焊缝区的单元按质心X坐标进行从低到高排序****************************esel,none*do,i,1,a*if,ANE1(i),NE,0,thenesel,a,elem,,ANE1(i)*endif*enddoMINE=0*DIM,ANEX1,ARRAY,a*DIM,ANEY1,ARRAY,a*DIM,ANEZ1,ARRAY,a*DO,I1,1, aesel,u,elem,,mine*GET,ANSEL1,ELEM,,COUNT1II=0*DO,I,1,AEMAX*IF,ESEL(I),EQ,1,THENII=II+1ANE1(II)=I*ENDIF*ENDDO*DO,I,1,ANSEL1*GET,ANEZ1(I),ELEM,ANE1(I),CENT,Z*GET,ANEY1(I),ELEM,ANE1(I),CENT,Y*GET,ANEX1(I),ELEM,ANE1(I),CENT,X *ENDDOMINZ=1E20MINY=1E20MINX=1E20*DO,I,1,ANSEL1*IF,ANEZ1(I),LT,MINZ,THEN MINZ=ANEZ1(I)MINY=ANEY1(I)MINX=ANEX1(I)MINE=ANE1(I)*ELSE*IF,ANEZ1(I),EQ,MINZ,THEN *IF,ANEY1(I),LT,MINY,THEN MINZ=ANEZ1(I)MINY=ANEY1(I)MINX=ANEX1(I)MINE=ANE1(I)*ELSE*IF,ANEY1(I),EQ,MINY,THEN *IF,ANEX1(I),LT,MINX,THEN MINZ=ANEZ1(I)MINY=ANEY1(I)MINX=ANEX1(I)MINE=ANE1(I)*ENDIF*ENDIF*ENDIF*ENDIF*ENDIF*ENDDOANEORDER(I1)=MINE*ENDDO!*************************对2道焊缝区的单元按质心X坐标进行从高到低排序****************************MAXE=1esel,none*do,j,1,b*if,ANE2(j),NE,0,thenesel,a,elem,,ANE2(j)*endif*enddo*DIM,ANEX2,ARRAY,b*DIM,ANEY2,ARRAY,b*DIM,ANEZ2,ARRAY,b*DO,J1,1, besel,u,elem,,MAXE*GET,ANSEL2,ELEM,,COUNT1JJ=0*DO,J,1,AEMAX*IF,ESEL(I),EQ,1,THENJJ=JJ+1ANE2(JJ)=J*ENDIF*ENDDO*DO,J,1,ANSEL2*GET,ANEZ2(J),ELEM,ANE2(J),CENT,Z*GET,ANEY2(J),ELEM,ANE2(J),CENT,Y*GET,ANEX2(J),ELEM,ANE2(J),CENT,X*ENDDOMAXZ=1E20MAXY=1E20MAXX=-1E20*DO,J,1,ANSEL2*IF,ANEZ2(J),LT,MAXZ,THEN MAXZ=ANEZ2(J)MAXY=ANEY2(J)MAXX=ANEX2(J)MAXE=ANE2(J)*ELSE*IF,ANEZ2(J),EQ,MAXZ,THEN *IF,ANEY2(J),LT,MAXY,THEN MAXZ=ANEZ2(J)MAXY=ANEY2(J)MAXX=ANEX2(J)MAXE=ANE2(J)*ELSE*IF,ANEY2(J),EQ,MAXY,THEN *IF,ANEX2(J),GT,MAXX,THENMAXZ=ANEZ2(J)MAXY=ANEY2(J)MAXX=ANEX2(J)MAXE=ANE2(J)*ENDIF*ENDIF*ENDIF*ENDIF*ENDIF*ENDDOANEORDER(J1+a)=MAXE*ENDDO!*************************对3道焊缝区的单元按质心X坐标进行从低到高排序****************************MINE=1esel,none*do,k,1,c*if,ANE3(k),NE,0,thenesel,a,elem,,ANE3(k)*endif*enddo*DIM,ANEX3,ARRAY,c*DIM,ANEY3,ARRAY,c*DIM,ANEZ3,ARRAY,c*DO,K1,1, cesel,u,elem,,mine*GET,ANSEL3,ELEM,,COUNT1KK=0*DO,K,1,AEMAX*IF,ESEL(K),EQ,1,THENKK=KK+1ANE3(KK)=K*ENDIF*ENDDO*DO,K,1,ANSEL3*GET,ANEZ3(K),ELEM,ANE3(K),CENT,Z*GET,ANEY3(K),ELEM,ANE3(K),CENT,Y*GET,ANEX3(K),ELEM,ANE3(K),CENT,X *ENDDOMINZ=1E20MINY=1E20MINX=1E20*DO,K,1,ANSEL3*IF,ANEZ3(K),LT,MINZ,THENMINZ=ANEZ3(K)MINY=ANEY3(K)MINX=ANEX3(K)MINE=ANE3(K)*ELSE*IF,ANEZ3(K),EQ,MINZ,THEN*IF,ANEY3(K),LT,MINY,THENMINZ=ANEZ3(K)MINY=ANEY3(K)MINX=ANEX3(K)MINE=ANE3(K)*ELSE*IF,ANEY3(K),EQ,MINY,THEN*IF,ANEX3(K),LT,MINX,THENMINZ=ANEZ3(K)MINY=ANEY3(K)MINX=ANEX3(K)MINE=ANE3(K)*ENDIF*ENDIF*ENDIF*ENDIF*ENDIF*ENDDOANEORDER(K1+a+b)=MINE*ENDDO!*************************对4道焊缝区的单元按质心X坐标进行从高到低排序****************************MAXE=1esel,none*do,j,1,d*if,ANE4(j),NE,0,thenesel,a,elem,,ANE4(j)*endif*enddo*DIM,ANEX4,ARRAY,d*DIM,ANEY4,ARRAY,d*DIM,ANEZ4,ARRAY,d*DO,J1,1, desel,u,elem,,MAXE*GET,ANSEL4,ELEM,,COUNT1JJ=0*DO,J,1,AEMAX*IF,ESEL(I),EQ,1,THENJJ=JJ+1ANE2(JJ)=J*ENDIF*ENDDO*DO,J,1,ANSEL4*GET,ANEZ4(J),ELEM,ANE4(J),CENT,Z*GET,ANEY4(J),ELEM,ANE4(J),CENT,Y *GET,ANEX4(J),ELEM,ANE4(J),CENT,X *ENDDOMAXZ=1E20MAXY=1E20MAXX=-1E20*DO,J,1,ANSEL4*IF,ANEZ4(J),LT,MAXZ,THENMAXZ=ANEZ4(J)MAXY=ANEY4(J)MAXX=ANEX4(J)MAXE=ANE4(J)*ELSE*IF,ANEZ4(J),EQ,MAXZ,THEN*IF,ANEY4(J),LT,MAXY,THENMAXZ=ANEZ4(J)MAXY=ANEY4(J)MAXX=ANEX4(J)MAXE=ANE4(J)*ELSE*IF,ANEY4(J),EQ,MAXY,THEN*IF,ANEX4(J),GT,MAXX,THENMAXZ=ANEZ4(J)MAXY=ANEY4(J)MAXX=ANEX4(J)MAXE=ANE4(J)*ENDIF*ENDIF*ENDIF*ENDIF*ENDIF*ENDDOANEORDER(J1+a+b+c)=MAXE*ENDDO!*************************对5道焊缝区的单元按质心X坐标进行从低到高排序****************************esel,none*do,i,1,e*if,ANE5(i),NE,0,thenesel,a,elem,,ANE5(i)*endif*enddoMINE=1*DIM,ANEX5,ARRAY,e*DIM,ANEY5,ARRAY,e*DIM,ANEZ5,ARRAY,e*DO,I1,1, eesel,u,elem,,mine*GET,ANSEL5,ELEM,,COUNT1II=0*DO,I,1,AEMAX*IF,ESEL(I),EQ,1,THENII=II+1ANE5(II)=I*ENDIF*ENDDO*DO,I,1,ANSEL5*GET,ANEZ5(I),ELEM,ANE5(I),CENT,Z*GET,ANEY5(I),ELEM,ANE5(I),CENT,Y *GET,ANEX5(I),ELEM,ANE5(I),CENT,X *ENDDOMINZ=1E20MINY=1E20MINX=1E20*DO,I,1,ANSEL1*IF,ANEZ5(I),LT,MINZ,THENMINZ=ANEZ5(I)MINY=ANEY5(I)MINX=ANEX5(I)MINE=ANE5(I)*ELSE*IF,ANEZ5(I),EQ,MINZ,THEN*IF,ANEY5(I),LT,MINY,THENMINZ=ANEZ5(I)MINY=ANEY5(I)MINX=ANEX5(I)MINE=ANE5(I)*ELSE*IF,ANEY5(I),EQ,MINY,THEN*IF,ANEX5(I),LT,MINX,THENMINZ=ANEZ5(I)MINY=ANEY5(I)MINX=ANEX5(I)MINE=ANE5(I)*ENDIF*ENDIF*ENDIF*ENDIF*ENDIF*ENDDOANEORDER(I1+a+b+c+d)=MINE*ENDDO!成形件轨迹单元质心排序FINISH!*********************体热源加载求解**************************************/SOLANTYPE,4 !瞬态分析TRNOPT,FULL !指定瞬态分析选项eqslv !指定求解器solcontrol,on !求解控制(用优化求解器)ESEL,ALLEPLOT/AUTO,1/REPLOTallsel,alloutres,all,allTOFFST,273tunif,24kbc,1timint,onnsub1=2nsub2=20asel,s,loc,z,Block_z !加对流条件asel,a,loc,x,0asel,a,loc,x,Block_x,. asel,a,loc,y,0asel,a,loc,y,Block_ynslasf,all,conv,-1,24allseldt=3e-2 ! 定义每次成形所需要的时间tim=0!生死单元控制!*DO,I1,1,AEMAX! ESEL,S,,,ANEORDER(I1)! EKILL,ALL! ESEL,S,LIVE! EPLOT!*enddoj=aemax !j小于AEmax*DO,I1,1,jnsubst,,,,1esel,s,,,ANEORDER(I1)! EALIVE,ANEORDER(I1-1)! ESEL,S,LIVEEPLOT! *ENDIFesel,allEPLOT!***************体载荷加载*********tim=tim+dttime,tim!bfe,ANEORDER(I1),hgen,1,3.91e12bfe,ANEORDER(I1),hgen,1,4.63e12*0.7!0.6 T_end=300 !*************************************SOLVEbfedele,ANEORDER(I1),hgen*ENDDO/soluANTYPE,,restbfedele,all,hgentime,4 !step101开始完全冷却NSUBST,15,30,5solvesavetime,5 !step102开始完全冷却NSUBST,15,30,5solvesavetime,6 !step103开始完全冷却NSUBST,15,30,5solvesavetime,7 !step104开始完全冷却NSUBST,15,30,5solvesavetime,8 !step105开始完全冷却NSUBST,15,30,5solvesavetime,9 !step106开始完全冷却NSUBST,15,30,5solvesavetime,10 !step107开始完全冷却NSUBST,15,30,5solvesavetime,30 !step108开始完全冷却NSUBST,15,30,5solvesavetime,300 !step109开始完全冷却NSUBST,15,30,5solvesave!******************************生成节点温度分布云图动画文件的程序段********************************/VIEW,1,1,2,3/ANG,1/REP,FAST/post1/seg,dele!/cont,1,15,0,3500/16,3500/dscale,1,1.0avprin,0,0avres,1/seg,multi,temp1,0.1esel,all*do,i,1,37set,(i-1)*3+1plnsol,temp*enddo/seg,off,temp37,1/anfile,save,temp37,avi/image,save,temp37,jpeganim,1,1,1finish!***********************************进入后处理*****************************************************/post26 ! 进入时间历程后处理器nsol,2,138,temp,,temp138 ! 取出(3.9，0.9,0.3)的温度值并赋给变量2号nsol,3,163,temp,,temp163 ! 取出(3.9，2.4,0.3)的温度值并赋给变量3号nsol,4,249,temp,,temp249 ! 取出(6.6，1.2,0.3)的温度值并赋给变量3号nsol,5,154,temp,,temp154 ! 取出(6.6，2.4,0.3)的温度值并赋给变量3号!nsol,4,149,temp,,temp149 ! 取出(6.9，1.2,0.3)的温度值并赋给变量3号!nsol,5,148,temp,,temp148 ! 取出(6.9，2.4,0.3)的温度值并赋给变量3号/gropt,axnsc,1.5, ! 缩放图标字体的大小plvar,2,3 ! 图形显示变量2,3号!prvar,2,3/image,save,temp23,jpeg ! 抓图plvar,4,5 ! 图形显示变量2,3号!prvar,4,5/image,save,temp45,jpeg ! 抓图plvar,2,4 ! 图形显示变量2,3号!prvar,2,4/image,save,temp24,jpeg ! 抓图finish。
(完整版)ANSYS命令流使用方法(中文)修改
Finish(退出四大模块，回到BEGIN层)/clear (清空内存，开始新的计算)1．定义参数、数组，并赋值.2．/prep7(进入前处理)定义几何图形：关键点、线、面、体定义几个所关心的节点，以备后处理时调用节点号。
设材料线弹性、非线性特性设置单元类型及相应KEYOPT设置实常数设置网格划分，划分网格根据需要耦合某些节点自由度定义单元表3．/solu加边界条件设置求解选项定义载荷步求解载荷步4./post1（通用后处理）5./post26 （时间历程后处理）6.PLOTCONTROL菜单命令7.参数化设计语言8.理论手册Finish(退出四大模块，回到BEGIN层)/clear (清空内存，开始新的计算)1.定义参数、数组，并赋值.dim, par, type, imax, jmax, kmax, var1, vae2, var3 定义数组par: 数组名type：array 数组，如同fortran,下标最小号为1，可以多达三维（缺省）char 字符串组（每个元素最多8个字符）tableimax,jmax, kmax 各维的最大下标号var1,var2,var3 各维变量名，缺省为row,column,plane(当type为table时) 2./prep7(进入前处理)2.1 设置单元类型及相应KEYOPTET, itype, ename, kop1……kop6, inopr 设定当前单元类型Itype：单元号Ename：单元名设置实常数Keyopt, itype, knum, valueitype: 已定义的单元类型号knum: 单元的关键字号value: 数值注意：如果,则必须使用keyopt命令，否则也可在ET命令中输入2.2 定义几个所关心的节点，以备后处理时调用节点号。
n,node,x,y,z,thxy, thyz, thzx 根据坐标定义节点号如果已有此节点，则原节点被重新定义，一般为最大节点号。
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为半径。
9、*AFUN，Lab在参数表达式中，为角度函数指定单位．Lab：指定将要使用的角度单位．有３个选项．RAD：在角度函数的输入与输出中使用弧度单位（默认）DEG：在角度函数的输入与输出中使用度单位．STAT：显示该命令当前的设置（即是度还是弧度）．10、Agen, itime,na1,na2,ninc,dx,dy,dz,kinc,noelem,imove ！面积复制命令。
Ansys命令流大全
Ansys命令流大全1、A,P1,P2,P3,P4,P5,P6,P7,P8,P9此命令用已知的一组关键点点（P1~P9）来定义面（Area），最少使用三个点才能围成面，同时产生转围绕些面的线。
点要依次序输入，输入的顺序会决定面的法线方向。
如果超过四个点，则这些点必须在同一个平面上。
Menu Paths:MainMenu>Preprocessor>Create>Arbitrary>Through KPs2、*ABBR，Abbr,String－－定义一个缩略语．Abbr:用来表示字符串＂String＂的缩略语，长度不超过８个字符．String：将由＂Abbr＂表示的字符串，长度不超过６０个字符．3、ABBRES，Lab，Fname，Ext－从一个编码文件中读出缩略语．Lab：指定读操作的标题，NEW：用这些读出的缩略语重新取代当前的缩略语（默认）CHANGE：将读出的缩略语添加到当前缩略语阵列，并替代现存同名的缩略语．Ext：如果＂Fname＂是空的，则缺省的扩展命是＂ABBR＂．4、ABBSAV，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热分析教程及命令流算例
＝＝＝＝＝【热力耦合分析单元简介】＝＝＝＝＝＝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高斯热源加载步骤及命令流
1.建模尺寸为5mm×5mm×0.5mm的SiCpAl材料（此处特别注意坐标系的位置，建议用DM建模）2.网格的划分实现移动热源路径网格划分密集，其他部分网格划分尺寸正常，方便计算，提高效率。
下图为网格划分示意图及网格划分设置。
3.分析计算（1）入射激光功率密度分布为高斯分布：I(r)=(1−R)2Pπω2exp(−2r2ω2)，（1）式中R为反射率，P为激光输出功率，ω为入射激光光斑半径，r为横向柱坐标，r=0处为光斑中心。
对模型施加负载条件时，将高斯热流密度用ANSYS APDL编写命令流加载到SiCp/Al上表面。
已知条件：功率P，行走速度v=5mm/s，光斑半径ω=0.05mm，时间t=1s。
则有，高斯移动热源：Q=2Pπ0.000052⋅exp(−2×(x−0.0025)2+(y−0.005×t)20.000052)×90%步骤：①首先在Mechanical APDL经典界面施加创建高斯热源函数的命令流；将上述公式带入，如下图所示。
②完成好函数输入之后，保存函数；然后读入刚刚保存的函数，命名为HFLUX，如下图所示：③提取命令流，命令流在文末（2）选择要加载的面，重命名为A1（3）载荷步设置在Step Controls（步长控制）中，Step End Time（结束时间）设置为1s，打开Auto Time Stepping(自动划分时间)，在Define By中设置为Substeps（子步），设置为100子步。
如下图所示。
（4）边界设置取模型下表面为绝缘边界，其他表面为对流换热面，对流换热系数取10W/(m3·K)，设置实验环境温度为22℃。
4.温度场仿真结果及分析命令流*DEL,_FNCNAME*DEL,_FNCMTID*DEL,_FNCCSYS*SET,_FNCNAME,'HFLUX'*SET,_FNCCSYS,0! /INPUT,1.func,,,1*DIM,%_FNCNAME%,TABLE,6,24,1,,,,%_FNCCSYS%!! Begin of equation: 4e8*exp(-2*(({X}-0.0025)^2+({Y}-0.005* ! {TIME})^2)/0.00005^2)*SET,%_FNCNAME%(0,0,1), 0.0, -999*SET,%_FNCNAME%(2,0,1), 0.0*SET,%_FNCNAME%(3,0,1), 0.0*SET,%_FNCNAME%(4,0,1), 0.0*SET,%_FNCNAME%(5,0,1), 0.0*SET,%_FNCNAME%(6,0,1), 0.0*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0, 0, 0, 0*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 0, 0, -1*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 2, -2*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 2, 0, 0, -3*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 3, -1*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0.0025, 0, 0, 2*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 2, 2, -1*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 2, 0, 0, -3*SET,%_FNCNAME%(0,9,1), 0.0, -4, 0, 1, -3, 17, -1*SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 0.005, 0, 0, 1*SET,%_FNCNAME%(0,11,1), 0.0, -3, 0, 1, -1, 3, 1*SET,%_FNCNAME%(0,12,1), 0.0, -1, 0, 1, 3, 2, -3*SET,%_FNCNAME%(0,13,1), 0.0, -3, 0, 2, 0, 0, -1*SET,%_FNCNAME%(0,14,1), 0.0, -5, 0, 1, -1, 17, -3*SET,%_FNCNAME%(0,15,1), 0.0, -1, 0, 1, -4, 1, -5*SET,%_FNCNAME%(0,16,1), 0.0, -3, 0, 1, -2, 3, -1*SET,%_FNCNAME%(0,17,1), 0.0, -1, 0, 0.00005, 0, 0, 0*SET,%_FNCNAME%(0,18,1), 0.0, -2, 0, 2, 0, 0, -1*SET,%_FNCNAME%(0,19,1), 0.0, -4, 0, 1, -1, 17, -2*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 1, -3, 4, -4*SET,%_FNCNAME%(0,21,1), 0.0, -1, 7, 1, -1, 0, 0*SET,%_FNCNAME%(0,22,1), 0.0, -2, 0, 4e8, 0, 0, -1*SET,%_FNCNAME%(0,23,1), 0.0, -3, 0, 1, -2, 3, -1*SET,%_FNCNAME%(0,24,1), 0.0, 99, 0, 1, -3, 0, 0! End of equation: 4e8*exp(-2*(({X}-0.0025)^2+({Y}-0.005*{TIME})^2)/0.00005^2) !-->SF,A1,HFLUX,%HFLUX%。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
!非线性材料求解成功GAUSS分布Z向扫描!!考虑了和基板对流情况下加载热流密度的热分析rb=0.15mm!粉床7.8X3.3mm，烧结区域6X1.5mm!考虑粉床表面的辐射、对流!考虑粉末状、熔化状（液态）、凝固态的边界条件变化!考虑金属蒸发会带走一部分能量，故效率减半eff＝0.7!r=0.2mmFINISH/CLEAR,START/BATCH !设置程序模式为“批处理”/CWD,'d:\Ansyswork\EBM FEM Simulation for Metals\3D Temperture distribution' !更改当前工作目录/FILNAME,EBM_db,0/TITLE,3D Temperture distrbution!选择分析问题的类型：热/NOPR/PMETH,OFF,0KEYW,PR_SET,1KEYW,PR_STRUC,0KEYW,PR_THERM,1KEYW,PR_FLUID,0KEYW,PR_ELMAG,0KEYW,MAGNOD,0KEYW,MAGEDG,0KEYW,MAGHFE,0KEYW,MAGELC,0KEYW,PR_MULTI,0KEYW,PR_CFD,0/GO/UNITS,SI !Standard Units!****************************************************************************** *******************!!开始，进入前处理器，定义下列属性－!0参数定义!1定义单元类型!2实常数!3材料性能!4几何模型!5网格划分控制!6其他!****************************************************************************** *******************!/PREP7!0参数定义!模型参数!几何模型尺寸me-3=mm*SET,Block_x,7.8e-3*SET,Block_x1,0.9e-3*SET,Block_x2,6.9e-3*SET,Block_y,3.3e-3*SET,Block_y1,0.9e-3*SET,Block_y2,2.4e-3*SET,Block_z,0.3e-3*SET,Block_z1,0.9e-3!划分网格的大小*SET,LMeshx,3e-4 !line1*SET,LMeshy,3e-4 !line4*SET,LMeshz,3e-4 !line9!边界条件参数*SET,Tinit,30 !初始温度!1定义单元类型ET,1,SOLID70ET,2,surf152 !表面效应单元r,2KEYOPT,2,4,1 !表面效应单元设置KEYOPT,2,8,4ET,3,SOLID87!2实常数!Nothing!3材料性能MPTEMPMPTEMP,1,25,400,1200,1600,2000,2500MPDATA,DENS,1,1,3983,4688.4,5406.8,5948.8,6908,6331 !密度MPDATA,KXX,1,1,1.6278,2.3812,3.6396,29.8,30.6,32 !热导率MPDATA,c,1,1,472,535,679,859,945,1034 !比热MPTEMP,,,,,,,,MPTEMP,1,0,200,400,800,1200MPTEMP,6,1400,1454,1800,2000MPDATA,ENTH,1,1,0,7.9729e+8,1.6400e+9,3.4621e+9,5.3897e+9 MPDATA,ENTH,1,6,6.3415e+9,8.2595e+9,9.8798e+9,1.0816e+10 !热焓MPTEMP,,,,,,,,MPTEMP,1,20,300,600,900,1200,1500,1800,2000,3000,4000 MPDATA,HF,1,1,6,50,120,200,250,378,700,850,1746,3799 !对流系数MPTEMP,,,,,,,,MPTEMP,1,20,200,400,600,800,900,1000,1100,1200,1420,1460MPDATA,DENS,2,1,7966,7893,7814,7724,7630,7583,7535,7486,7436,7320,7320MPDATA,KXX,2,1,13.31,16.33,19.47,22.38,25.07,26.33,27.53,28.67,29.76,31.95,32 !热导率mpdata,c,2,1,470,508,550,592,634,655,676,698,719,765,765MPTEMP,,,,,,,,MPTEMP,1,0,200,400,800,1200MPTEMP,6,1400,1454,1800,2000MPDATA,ENTH,2,1,0,7.9729e+8,1.6400e+9,3.4621e+9,5.3897e+9MPDATA,ENTH,2,6,6.3415e+9,8.2595e+9,9.8798e+9,1.0816e+10!4几何模型BLOCK,0,Block_x,0,Block_y,0,Block_z,BLOCK,Block_x1,Block_x2,Block_y1,Block_y2,0,Block_z,vsel,allVOVLAP,allBLOCK,0,Block_x,0,Block_y,0,-Block_z1,vglue,allnumcmp,all!划分网格lesize,8,LMeshxlesize,5,LMeshylesize,13,LMeshzvsel,s,,,1type,1mat,1mshkey,1vmesh,allallselesize,0.3e-3mshkey,0mshape,1,3dvsel,s,,,2type,3mat,2vmesh,allallselesize,0.3e-3mshkey,0mshape,1,3dvsel,s,,,3type,3mat,1vmesh,allallsel!成形件轨迹单元质心排序width=Block_y2-Block_y1!wpoffs,Block_x1,Block_y1,0!local,11,1,Block_x1,Block_y1,0 ALLSEL,ALL*GET,EMAX,ELEM,,COUNT,MAXVSEL,S, , , 1ESLV,S*GET,AEMAX,ELEM,,COUNT1,MAX*DIM,ANE,ARRAY,AEMAX,*DIM,ANE1,ARRAY,AEMAX,*DIM,ANE2,ARRAY,AEMAX,*DIM,ANE3,ARRAY,AEMAX,*DIM,ANE4,ARRAY,AEMAX,*DIM,ANE5,ARRAY,AEMAX,*DIM,ANEX,ARRAY,AEMAX*DIM,ANEY,ARRAY,AEMAX*DIM,ANEZ,ARRAY,AEMAX*DIM,ANEORDER,ARRAY,AEMAX,,1!形心排序VSEL,S, , , 1ESLV,S*GET,ANSEL,ELEM,,COUNT1II=0*DO,I,1,AEMAX*IF,ESEL(I),EQ,1,THENII=II+1ANE(II)=I*ENDIF*ENDDO*DO,I,1,ANSEL*GET,ANEZ(I),ELEM,ANE(I),CENT,Z*GET,ANEY(I),ELEM,ANE(I),CENT,Y*GET,ANEX(I),ELEM,ANE(I),CENT,X *ENDDOe=0*do,i,1,AEMAX*if,ANEY(I),ge,4*width/5+Block_y1,then e=e+1ANE5(e)=ANE(i)*endif*enddod=0*do,i,1,AEMAX*if,ANEY(I),ge,3*width/5+Block_y1,then *if,ANEY(I),lt,4*width/5+Block_y1,then d=d+1ANE4(d)=ANE(i)*endif*endif*enddoc=0*do,i,1,AEMAX*if,ANEY(I),ge,2*width/5+Block_y1,then *if,ANEY(I),lt,3*width/5+Block_y1,then c=c+1ANE3(c)=ANE(i)*endif*endif*enddob=0*do,i,1,AEMAX*if,ANEY(I),ge,width/5+Block_y1,then*if,ANEY(I),lt,2*width/5+Block_y1,then b=b+1ANE2(b)=ANE(i)*endif*endif*enddoa=0*do,i,1,AEMAX*if,ANEY(I),lt,width/5+Block_y1,thena=a+1ANE1(a)=ANE(i)*endif*enddo!*************************对1道焊缝区的单元按质心X坐标进行从低到高排序****************************esel,none*do,i,1,a*if,ANE1(i),NE,0,thenesel,a,elem,,ANE1(i)*endif*enddoMINE=0*DIM,ANEX1,ARRAY,a*DIM,ANEY1,ARRAY,a*DIM,ANEZ1,ARRAY,a*DO,I1,1, aesel,u,elem,,mine*GET,ANSEL1,ELEM,,COUNT1II=0*DO,I,1,AEMAX*IF,ESEL(I),EQ,1,THENII=II+1ANE1(II)=I*ENDIF*ENDDO*DO,I,1,ANSEL1*GET,ANEZ1(I),ELEM,ANE1(I),CENT,Z*GET,ANEY1(I),ELEM,ANE1(I),CENT,Y*GET,ANEX1(I),ELEM,ANE1(I),CENT,X*ENDDOMINZ=1E20MINY=1E20MINX=1E20*DO,I,1,ANSEL1*IF,ANEZ1(I),LT,MINZ,THENMINZ=ANEZ1(I)MINY=ANEY1(I)MINX=ANEX1(I)MINE=ANE1(I)*ELSE*IF,ANEZ1(I),EQ,MINZ,THEN*IF,ANEY1(I),LT,MINY,THENMINZ=ANEZ1(I)MINY=ANEY1(I)MINX=ANEX1(I)MINE=ANE1(I)*ELSE*IF,ANEY1(I),EQ,MINY,THEN*IF,ANEX1(I),LT,MINX,THENMINZ=ANEZ1(I)MINY=ANEY1(I)MINX=ANEX1(I)MINE=ANE1(I)*ENDIF*ENDIF*ENDIF*ENDIF*ENDIF*ENDDOANEORDER(I1)=MINE*ENDDO!*************************对2道焊缝区的单元按质心X坐标进行从高到低排序****************************MAXE=1esel,none*do,j,1,b*if,ANE2(j),NE,0,thenesel,a,elem,,ANE2(j)*endif*enddo*DIM,ANEX2,ARRAY,b*DIM,ANEY2,ARRAY,b*DIM,ANEZ2,ARRAY,b*DO,J1,1, besel,u,elem,,MAXE*GET,ANSEL2,ELEM,,COUNT1JJ=0*DO,J,1,AEMAX*IF,ESEL(I),EQ,1,THENJJ=JJ+1ANE2(JJ)=J*ENDIF*ENDDO*DO,J,1,ANSEL2*GET,ANEZ2(J),ELEM,ANE2(J),CENT,Z*GET,ANEY2(J),ELEM,ANE2(J),CENT,Y*GET,ANEX2(J),ELEM,ANE2(J),CENT,X *ENDDOMAXZ=1E20MAXY=1E20MAXX=-1E20*DO,J,1,ANSEL2*IF,ANEZ2(J),LT,MAXZ,THENMAXZ=ANEZ2(J)MAXY=ANEY2(J)MAXX=ANEX2(J)MAXE=ANE2(J)*ELSE*IF,ANEZ2(J),EQ,MAXZ,THEN*IF,ANEY2(J),LT,MAXY,THENMAXZ=ANEZ2(J)MAXY=ANEY2(J)MAXX=ANEX2(J)MAXE=ANE2(J)*ELSE*IF,ANEY2(J),EQ,MAXY,THEN*IF,ANEX2(J),GT,MAXX,THENMAXZ=ANEZ2(J)MAXY=ANEY2(J)MAXX=ANEX2(J)MAXE=ANE2(J)*ENDIF*ENDIF*ENDIF*ENDIF*ENDIF*ENDDOANEORDER(J1+a)=MAXE*ENDDO!*************************对3道焊缝区的单元按质心X坐标进行从低到高排序****************************MINE=1esel,none*do,k,1,c*if,ANE3(k),NE,0,thenesel,a,elem,,ANE3(k)*endif*enddo*DIM,ANEX3,ARRAY,c*DIM,ANEY3,ARRAY,c*DIM,ANEZ3,ARRAY,c*DO,K1,1, cesel,u,elem,,mine*GET,ANSEL3,ELEM,,COUNT1KK=0*DO,K,1,AEMAX*IF,ESEL(K),EQ,1,THENKK=KK+1ANE3(KK)=K*ENDIF*ENDDO*DO,K,1,ANSEL3*GET,ANEZ3(K),ELEM,ANE3(K),CENT,Z*GET,ANEY3(K),ELEM,ANE3(K),CENT,Y*GET,ANEX3(K),ELEM,ANE3(K),CENT,X*ENDDOMINZ=1E20MINY=1E20MINX=1E20*DO,K,1,ANSEL3*IF,ANEZ3(K),LT,MINZ,THENMINZ=ANEZ3(K)MINY=ANEY3(K)MINX=ANEX3(K)MINE=ANE3(K)*ELSE*IF,ANEZ3(K),EQ,MINZ,THEN*IF,ANEY3(K),LT,MINY,THENMINZ=ANEZ3(K)MINY=ANEY3(K)MINX=ANEX3(K)MINE=ANE3(K)*ELSE*IF,ANEY3(K),EQ,MINY,THEN*IF,ANEX3(K),LT,MINX,THENMINZ=ANEZ3(K)MINY=ANEY3(K)MINX=ANEX3(K)MINE=ANE3(K)*ENDIF*ENDIF*ENDIF*ENDIF*ENDIF*ENDDOANEORDER(K1+a+b)=MINE*ENDDO!*************************对4道焊缝区的单元按质心X坐标进行从高到低排序****************************MAXE=1esel,none*do,j,1,d*if,ANE4(j),NE,0,thenesel,a,elem,,ANE4(j)*endif*enddo*DIM,ANEX4,ARRAY,d*DIM,ANEY4,ARRAY,d*DIM,ANEZ4,ARRAY,d*DO,J1,1, desel,u,elem,,MAXE*GET,ANSEL4,ELEM,,COUNT1JJ=0*DO,J,1,AEMAX*IF,ESEL(I),EQ,1,THENJJ=JJ+1ANE2(JJ)=J*ENDIF*ENDDO*DO,J,1,ANSEL4*GET,ANEZ4(J),ELEM,ANE4(J),CENT,Z*GET,ANEY4(J),ELEM,ANE4(J),CENT,Y*GET,ANEX4(J),ELEM,ANE4(J),CENT,X *ENDDOMAXZ=1E20MAXY=1E20MAXX=-1E20*DO,J,1,ANSEL4*IF,ANEZ4(J),LT,MAXZ,THENMAXZ=ANEZ4(J)MAXY=ANEY4(J)MAXX=ANEX4(J)MAXE=ANE4(J)*ELSE*IF,ANEZ4(J),EQ,MAXZ,THEN*IF,ANEY4(J),LT,MAXY,THENMAXZ=ANEZ4(J)MAXY=ANEY4(J)MAXX=ANEX4(J)MAXE=ANE4(J)*ELSE*IF,ANEY4(J),EQ,MAXY,THEN*IF,ANEX4(J),GT,MAXX,THENMAXZ=ANEZ4(J)MAXY=ANEY4(J)MAXX=ANEX4(J)MAXE=ANE4(J)*ENDIF*ENDIF*ENDIF*ENDIF*ENDIF*ENDDOANEORDER(J1+a+b+c)=MAXE*ENDDO!*************************对5道焊缝区的单元按质心X坐标进行从低到高排序****************************esel,none*do,i,1,e*if,ANE5(i),NE,0,thenesel,a,elem,,ANE5(i)*endif*enddoMINE=1*DIM,ANEX5,ARRAY,e*DIM,ANEY5,ARRAY,e*DIM,ANEZ5,ARRAY,e*DO,I1,1, eesel,u,elem,,mine*GET,ANSEL5,ELEM,,COUNT1II=0*DO,I,1,AEMAX*IF,ESEL(I),EQ,1,THENII=II+1ANE5(II)=I*ENDIF*ENDDO*DO,I,1,ANSEL5*GET,ANEZ5(I),ELEM,ANE5(I),CENT,Z*GET,ANEY5(I),ELEM,ANE5(I),CENT,Y*GET,ANEX5(I),ELEM,ANE5(I),CENT,X*ENDDOMINZ=1E20MINY=1E20MINX=1E20*DO,I,1,ANSEL1*IF,ANEZ5(I),LT,MINZ,THENMINZ=ANEZ5(I)MINY=ANEY5(I)MINX=ANEX5(I)MINE=ANE5(I)*ELSE*IF,ANEZ5(I),EQ,MINZ,THEN*IF,ANEY5(I),LT,MINY,THENMINZ=ANEZ5(I)MINY=ANEY5(I)MINX=ANEX5(I)MINE=ANE5(I)*ELSE*IF,ANEY5(I),EQ,MINY,THEN*IF,ANEX5(I),LT,MINX,THENMINZ=ANEZ5(I)MINY=ANEY5(I)MINX=ANEX5(I)MINE=ANE5(I)*ENDIF*ENDIF*ENDIF*ENDIF*ENDIF*ENDDOANEORDER(I1+a+b+c+d)=MINE*ENDDO!成形件轨迹单元质心排序FINISH!*********************体热源加载求解**************************************/SOLANTYPE,4 !瞬态分析TRNOPT,FULL !指定瞬态分析选项eqslv !指定求解器solcontrol,on !求解控制(用优化求解器)ESEL,ALLEPLOT/AUTO,1/REPLOTallsel,alloutres,all,allTOFFST,273tunif,24kbc,1timint,onnsub1=2nsub2=20asel,s,loc,z,Block_z !加对流条件asel,a,loc,x,0asel,a,loc,x,Block_xasel,a,loc,y,0asel,a,loc,y,Block_ynslasf,all,conv,-1,24allseldt=3e-2 ! 定义每次成形所需要的时间tim=0!生死单元控制!*DO,I1,1,AEMAX! ESEL,S,,,ANEORDER(I1)! EKILL,ALL! ESEL,S,LIVE! EPLOT!*enddoj=aemax !j小于AEmax*DO,I1,1,jnsubst,,,,1esel,s,,,ANEORDER(I1)! EALIVE,ANEORDER(I1-1)! ESEL,S,LIVEEPLOT! *ENDIFesel,allEPLOT!***************体载荷加载*********tim=tim+dttime,tim!bfe,ANEORDER(I1),hgen,1,3.91e12bfe,ANEORDER(I1),hgen,1,4.63e12*0.7!0.6 T_end=300 !************************************* SOLVEbfedele,ANEORDER(I1),hgen*ENDDO/soluANTYPE,,restbfedele,all,hgentime,4 !step101开始完全冷却NSUBST,15,30,5solvesavetime,5 !step102开始完全冷却NSUBST,15,30,5solvesavetime,6 !step103开始完全冷却NSUBST,15,30,5solvesavetime,7 !step104开始完全冷却NSUBST,15,30,5solvesavetime,8 !step105开始完全冷却NSUBST,15,30,5solvesavetime,9 !step106开始完全冷却NSUBST,15,30,5solvesavetime,10 !step107开始完全冷却NSUBST,15,30,5solvesavetime,30 !step108开始完全冷却NSUBST,15,30,5solvesavetime,300 !step109开始完全冷却NSUBST,15,30,5solvesave!******************************生成节点温度分布云图动画文件的程序段********************************/VIEW,1,1,2,3/ANG,1/REP,FAST/post1/seg,dele!/cont,1,15,0,3500/16,3500/dscale,1,1.0avprin,0,0avres,1/seg,multi,temp1,0.1esel,all*do,i,1,37set,(i-1)*3+1plnsol,temp*enddo/seg,off,temp37,1/anfile,save,temp37,avi/image,save,temp37,jpeganim,1,1,1finish!***********************************进入后处理*****************************************************/post26 ! 进入时间历程后处理器nsol,2,138,temp,,temp138 ! 取出(3.9，0.9,0.3)的温度值并赋给变量2号nsol,3,163,temp,,temp163 ! 取出(3.9，2.4,0.3)的温度值并赋给变量3号nsol,4,249,temp,,temp249 ! 取出(6.6，1.2,0.3)的温度值并赋给变量3号nsol,5,154,temp,,temp154 ! 取出(6.6，2.4,0.3)的温度值并赋给变量3号!nsol,4,149,temp,,temp149 ! 取出(6.9，1.2,0.3)的温度值并赋给变量3号!nsol,5,148,temp,,temp148 ! 取出(6.9，2.4,0.3)的温度值并赋给变量3号/gropt,axnsc,1.5, ! 缩放图标字体的大小plvar,2,3 ! 图形显示变量2,3号!prvar,2,3/image,save,temp23,jpeg ! 抓图plvar,4,5 ! 图形显示变量2,3号!prvar,4,5/image,save,temp45,jpeg ! 抓图plvar,2,4 ! 图形显示变量2,3号!prvar,2,4/image,save,temp24,jpeg ! 抓图finish。