ANSYS热应力分析命令流
(最新整理)ANSYS热分析详解
(完整)ANSYS热分析详解编辑整理：尊敬的读者朋友们：这里是精品文档编辑中心，本文档内容是由我和我的同事精心编辑整理后发布的，发布之前我们对文中内容进行仔细校对，但是难免会有疏漏的地方，但是任然希望（(完整)ANSYS热分析详解）的内容能够给您的工作和学习带来便利。
同时也真诚的希望收到您的建议和反馈，这将是我们进步的源泉，前进的动力。
本文可编辑可修改，如果觉得对您有帮助请收藏以便随时查阅，最后祝您生活愉快业绩进步，以下为(完整)ANSYS热分析详解的全部内容。
第一章简介一、热分析的目的热分析用于计算一个系统或部件的温度分布及其它热物理参数，如热量的获取或损失、热梯度、热流密度（热通量〕等。
热分析在许多工程应用中扮演重要角色,如内燃机、涡轮机、换热器、管路系统、电子元件等。
二、ANSYS的热分析•在ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Thermal、ANSYS/FLOTRAN、ANSYS/ED五种产品中包含热分析功能,其中ANSYS/FLOTRAN不含相变热分析。
•ANSYS热分析基于能量守恒原理的热平衡方程，用有限元法计算各节点的温度，并导出其它热物理参数。
•ANSYS热分析包括热传导、热对流及热辐射三种热传递方式.此外，还可以分析相变、有内热源、接触热阻等问题。
三、ANSYS 热分析分类•稳态传热:系统的温度场不随时间变化•瞬态传热：系统的温度场随时间明显变化四、耦合分析•热－结构耦合•热－流体耦合•热－电耦合•热－磁耦合•热－电－磁－结构耦合等第二章 基础知识一、符号与单位 W/m 2—℃ 二、传热学经典理论回顾热分析遵循热力学第一定律,即能量守恒定律:● 对于一个封闭的系统（没有质量的流入或流出〕PE KE U W Q ∆+∆+∆=-式中： Q —— 热量；W -- 作功；∆U ——系统内能；∆KE ——系统动能;∆PE —-系统势能；●对于大多数工程传热问题：0＝＝PE KE ∆∆； ●通常考虑没有做功:0=W ， 则：U Q ∆=； ● 对于稳态热分析：0=∆=U Q ，即流入系统的热量等于流出的热量;●对于瞬态热分析：dt dU q =,即流入或流出的热传递速率q 等于系统内能的变化。
ANSYS耦合场分析_热应力
例如: 如果结构网格包括在热模型中不存在的圆角 时，许多节点将落在热模型的外面。如果圆角足 够大而且热模型足够细致，圆角区域的载荷将不 能写出。
10-15
基本过程
在热-应力分析中，由温度求解得到的节点温度 将在结构分析中用作体载荷。 当在顺序求解使用手工方法时将热节点温度施加到结构单元上有两种选项。选择 的原则在于结构模型和热模型是否有相似的网格划分:
如果热和结构的单元有相同 的节点号码. . .
1
• 热模型自动转换为结构模型，使用 ETCHG 命令(见相应单元表格)。 • 温度可以直接从热分析结果文件读出 并使用LDREAD 命令施加到结构模型 上。
10-3
直接方法 - 例题
在第七章对流部分中，介绍了FLUID66和FLUID116热—流单元。该单元具有 热和压力自由度，因此是直接耦合场单元。
ANSYS有一些其他的耦合单元，具有结构，热，电，磁等自由度。绝大多数 的实际问题只涉及到少数几个物理场的耦合。这里提供了几个涉及到热现象 的直接耦合场分析。
• 不同场之间使用统一的单位制。例如，在热-电分析中，如果电瓦单位使 用瓦(焦耳/秒)，热单位就不能使用Btu/s。
• 由于需要迭代计算，热耦合场单元不能使用子结构。
10-6
直接方法 - 加载, 求解, 后处理
在直接方法的加载，求解，后处理中注意以下方面： • 如果对带有温度自由度的耦合场单元选择 瞬态 分析类型的话： – 瞬态温度效果可以在所有耦合场单元中使用。 – 瞬态电效果(电容，电感)不能包括在热-电分析中(除非只是TEMP和VOLT自由度 被 激活)。 – 带有磁向量势自由度的耦合场单元可以用来对瞬态磁场问题建模(如,SOLID62). 带 有标量势自由度的单元只能模拟静态现象(SOLID5)。 • 学习每种单元的自由度和允许的载荷。耦合场单元允许的相同位置(节点,单元面等)施加多 种类型的载荷 (D, F, SF, BF) 。 • 耦合场分析可以使高度非线性的。考虑使用Predictor 和 Line Search 功能改善收敛性。 • 考虑使用Multi-Plots功能将不同场的结果同时输出到多个窗口中。
ANSYS工程应用教程_热与电磁学篇
ANSYS 工程应用教程_热与电磁学篇随着ANSYS 版本的不断更新，ANSYS 的应用领域也日益广泛。
作为融结构、热、流体、电磁、声学为一体的大型通用有限元分析软件，可广泛应用于核工业、铁道、石油化工、航空航天、机械制造、能源、汽车交通、国防军工、电子、土木工程、造船、生物医学、轻工、地矿、水利、日用家电、等一般工业及科学研究领域。
热分析包括稳态热分析、瞬态热分析、热辐射、相变、热应力等，电磁场分析包括二维静态、谐性、瞬态磁场分析，三维静态、谐性、瞬态磁场分析，高频电磁场分析和电场分析等。
ANSYS 热分析简介：图形用户界面方式（GUI ）或命令流方式进行计算。
ANSYS 如何进行热分析：实际上，其基本原理是先将所处理的对象划分成有限个单元（包含若干节点），然后根据能量守恒原理求解一定边界条件和初始条件下每一节点处的热平衡方程，由此计算出各节点温度，继而进一步求解出其他相关量。
耦合场分析：这类涉及两个和多个物理场相互作用的问题为耦合场分析。
主要方法有直接耦合和间接耦合。
直接耦合解法的耦合单元包含所有的自由度，仅仅通过一次求解就能得出耦合场分析结果。
这种方法实际上是通过计算包含所有必须项的单元矩阵或单元载荷向量来实现的。
间接耦合法又称为序贯耦合法，通过把第一磁场分析的结果作为第二次场分析的载荷来实现良种场的耦合。
三种基本传热方式：传导：当物理内部存在温度差时，热量将从高温部分传递到低温部分；而且不同温度的物体相互接触时热量会从高温物体传递到低温物体。
傅立叶定律，又称导热基本定律hot cold A(T T )t dQ κ-=，Q 为时间t 内的传热量，κ为热传导率，T 为温度，A 为面积，d 为两平面之间的距离。
对流：温度不同的各部分流体之间发生相对运动所引起的热量传递方式。
流体被加热时：w f q h(t t )=-流体被冷却时：f w q h(t t )=-，w t 和f t 分别为壁面温度和流体温度，h 为对流热系数。
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杆单元单元名称简称节点数节点自由度特性备注LINK1 2D杆 2 Ux,Uy EPCSDGB常用杆元LINK8 3D杆Ux,Uy,Uz EPCSDGBLINK103D仅受拉或仅受压杆EDGB模拟缆索的松弛及间隙LINK11 3D线性调节器EGB模拟液压缸和大转动LINK180 3D有限应变杆EPCDFGB 另可考虑粘弹塑性E-弹性(Elasticity)，P-塑性(Plasticity)，C-蠕变(Creep)，S-膨胀(Swelling)，D-大变形或大挠度(Large deflection)，F-大应变(Large strain)或有限应变(Finite strain)，B-单元生死(Birth and dead),G-应力刚化(Stress stiffness)或几何刚度(Geometric stiffening)，A-自适应下降(Adaptive descent)等。
换热器热应力分析ansys命令流
! ***************环境设置***************FINISH/CLEAR/UNITS,SI/FILN,E102/TITLE, MECHANICAL AND THERMAL STRESS ANALYSIS OF FIXED TUBE-SHEET EXCHANGER! ********* 参数设定**********SET,TSHELL,17 !壳程筒体壁厚*SET,HSHELL,2500 !壳程筒体分析长度*SET,DSHELL,806 !壳程筒体内径*SET,HPSH,20 !壳程侧管板锻造凸肩高度*SET,RSHELL,15 !壳程侧管板锻造圆角半径*SET,SHIGH,100 !管板厚度*SET,TCHANNEL,18.5 !管程筒体壁厚*SET,HCHANNEL,350 !管程筒体分析长度*SET,HPCH,20 !管程侧管板锻造凸肩高度*SET,RCHANNEL,15 !管程侧管板锻造圆角半径*SET,DP,25 !换热管外直径*SET,SP,38 !换热管管间距*SET,TPIPE,2.5 !换热管壁厚*SET,TINSULAT,50 !保温层厚度*SET,NY,20 !换热管行数*DIM,NXZ,,NY !以行数定义一维数组个数，每个数组中存放该行换热管根数*SET,NXZ(1),6,6,6,6,5,6,5,6,5,5*SET,NXZ(11),5,5,4,5,4,4,3,3,2,1! ***************剖分参数设定****************SET,NGH,8 !换热管周向剖分数*SET,NGMESH,40 !换热管轴向剖分数*SET,DENSITYMESHG,2 !换热管轴向剖分加密倍数*SET,NR0,60 !管板环向剖分数*SET,NHOU,2 !壳程和管程筒体壁厚剖分数*SET,NINSULAT,3 !保温层壁厚剖分数*SET,NHSHELL,40 !壳程筒体分析长度剖分数*SET,NCHANNEL,8 !管程筒体分析长度剖分数*SET,NR,4 !管板锻造圆角周向剖分数*SET,NSHIGH,4 !管板厚度剖分数*SET,NHPSH,2 !壳程侧管板锻造凸肩高度剖分数*SET,NHPCH,2 !管程侧管板锻造凸肩高度剖分数! ***************载荷条件****************SET,PS,0.58 !壳程设计压力*SET,PT,2 !管程设计压力! ***************物性参数****************SET,ASHELL,0.001764 !壳程筒体内介质传热系数*SET,ACHANNEL,0.000191 !管程筒体内介质传热系数*SET,ATUBE,0.000553 !换热管内介质传热系数*SET,AENVIR,0.000006 !环境空气传热系数*SET,TEPTUBE,250 !管程介质温度*SET,TEPSHEL,140.5 !壳程介质温度*SET,TEPENVIR,20 !环境温度*SET,KSHELL,0.0371 !壳程筒体金属热导率*SET,KCHANNEL,0.0188 !管程筒体金属热导率*SET,KSHIGH,0.01618 !管板金属热导率*SET,KTUBE,0.0189 !换热管金属热导率*SET,KINSULAT,0.00025 !保温层岩棉热导率! ***************前处理***************/PREP7ET,1,SHELL57 !定义辅助建模的壳单元类型ET,2,SOLID70 !进行热分析的体单元类型MP,KXX,1,KSHELL !定义材料导热系数，1号材料对应筒体钢板MP,KXX,2,KCHANNEL !定义材料导热系数，2号材料对应管箱钢板MP,KXX,3,KSHIGH !定义材料导热系数，3号材料对应管板锻件MP,KXX,4,KTUBE !定义材料导热系数，4号材料对应换热管钢管MP,KXX,5,KINSULAT !定义材料导热系数，5号材料对应保温层岩棉! ***************建立模型***************管板下表面,四分之一*IF,DSHELL/2.0D0-RSHELL,LE,DSHELL/2.0D0+TSHELL-TCHANNEL-RCHANNEL,THEN !判断并确定一管板圆平面半径，其距两侧圆角边缘最小为5MMROL=DSHELL/2.0D0-RSHELL-5*ELSEROL=DSHELL/2.0D0+TSHELL-TCHANNEL-RCHANNEL-5*ENDIFPCIRC,,ROL,,90 !生成管板圆平面，其面号为1*AFUN,DEG !转换成角度制NTOL=0 !换热管总数计数器清零*DO,J,1,NY !利用循环产生以换热管外直径为直径的圆平面NTEMP=J*DO,I,1,NXZ(NTEMP)*IF,NTEMP,NE,NINT(NTEMP/2)*2,THEN !根据当前换热管行数的奇偶性确定工作面偏移的圆心坐标WPOFFS,SP*SIN(60)+2*SP*SIN(60)*(I-1),SP*(J-1)*COS(60) !行数J为奇数时该行第I根换热管圆心坐标*ELSEWPOFFS,2*SP*SIN(60)*(I-1),SP*(J-1)*COS(60) !行数J为偶数时该行第I根换热管圆心坐标*ENDIFPCIRC,,DP/2.0D0 !产生以换热管外直径为直径的圆平面NTOL=NTOL+1 !计数器加1WPCSYS,1,0 !工作面圆心坐标回零*ENDDO*ENDDOASEL,U,AREA,,1 !选中除管板圆平面的所有面CM,BASE1,AREA !将上一行选中的所有以换热管外直径为直径的圆平面命名为BASE1ALLSEL !全选ASBA,1,BASE1 !布尔减运算，产生布管区管板的下表面换热管下表面,四分之一*DO,J,1,NY !利用循环产生换热管下表面的圆环面，方法与上一个循环相同NTEMPJ=J*DO,I,1,NXZ(NTEMPJ)NTEMPI=I*IF,NTEMPJ,NE,NINT(NTEMPJ/2)*2,THENWPOFFS,SP*SIN(60)+2*SP*SIN(60)*(I-1),SP*(J-1)*COS(60)*ELSEWPOFFS,2*SP*SIN(60)*(I-1),SP*(J-1)*COS(60)*ENDIF*IF,NTEMPJ,EQ,1,THENPCIRC,DP/2.0D0-TPIPE,DP/2.0D0,,90PCIRC,DP/2.0D0-TPIPE,DP/2.0D0,90,180*ELSE*IF,NTEMPJ,EQ,NINT(NTEMPJ/2)*2,THEN*IF,NTEMPI,EQ,1,THENPCIRC,DP/2.0D0-TPIPE,DP/2.0D0,270,360PCIRC,DP/2.0D0-TPIPE,DP/2.0D0,,90*ELSEPCIRC,DP/2.0D0-TPIPE,DP/2.0D0*ENDIF*ELSEPCIRC,DP/2.0D0-TPIPE,DP/2.0D0*ENDIF*ENDIFWPCSYS,1,0*ENDDO*ENDDOAGLUE,ALL !将新产生的换热管圆环面与管板平面粘和为一体NUMCMP,AREA !使所有面的面号从小到大重新排序NUMCMP,LINE !使所有线的线号从小到大重新排序LSEL,S,RADIUS,,ROL !选中管板圆平面的外边缘线ASLL,S,0 !选中与此线相关的面，即布管区管板的下表面CM,BASPLT,AREA !给布管区管板的下表面命名为BASPLTALLSEL !全选ASEL,U,AREA,,BASPLT !选中所有换热管下表面的圆环面CM,BASTUB,AREA !给换热管下表面的圆环面命名为BASTUBALLSEL !全选网格剖分参数LSEL,S,RADIUS,,ROL !选中管板圆平面的外边缘线LESIZE,ALL,,,NR0 !指定该线的剖分数ALLSEL !全选LSEL,S,RADIUS,,DP/2.0D0-TPIPE,DP/2.0D0 !选中构成换热管下表面的圆环面的圆环线LESIZE,ALL,,,NGH/4 !指定这些线的剖分数ALLSEL !全选LSEL,S,LENGTH,,TPIPE !选中长度为换热管壁厚的线LESIZE,ALL,,,1 !换热管壁厚方向剖分1份ALLSEL !全选布管区管板体MSHAPE,0,2D !指定面单元形状为四边形TYPE,1 !指定面单元类型ESIZE,SP-DP !定义单元尺寸MAT,3 !设定材料属性AMESH,BASPLT !对布管区管板的下表面进行网格划分K,10001,0,0,0 !定义拉伸布管区管板的下表面所需的关键点K,10002,0,0,TCHANNEL !定义拉伸布管区管板的下表面所需的关键点K,10003,0,0,SHIGH-TSHELL!定义拉伸布管区管板的下表面所需的关键点K,10004,0,0,SHIGH !定义拉伸布管区管板的下表面所需的关键点NUMCMP,LINE !使所有线的线号从小到大重新排序L,10001,10002,NHOU !生成布管区管板的下表面的拉伸辅助线，并指定剖分数L,10002,10003,NSHIGH !生成布管区管板的下表面的拉伸辅助线，并指定剖分数L,10003,10004,NHOU !生成布管区管板的下表面的拉伸辅助线，并指定剖分数*GET,LINE0,LINE,,NUM,MAX!得到辅助线的线号TYPE,2 !指定体单元类型VDRAG,BASPLT,,,,,,LINE0-2,LINE0-1,LINE0 !产生布管区管板的体单元ACLEAR,BASPLT !清除布管区管板的下表面的面单元换热管体MSHAPE,0,2D !指定面单元形状为四边形TYPE,1 !指定面单元类型MAT,4 !设定材料属性AMESH,BASTUB !对换热管下表面的圆环面进行网格划分K,10005,0,0,SHIGH+RSHELL+HPSH+HSHELL !定义拉伸换热管下表面的圆环面所需的关键点L,10004,10005 !生成换热管下表面的圆环面的拉伸辅助线*GET,LINE1,LINE,,NUM,MAX!得到辅助线的线号LESIZE,LINE1,,,NGMESH,DENSITYMESHG !指定辅助线剖分份数和加密倍数TYPE,2 !指定体单元类型VDRAG,BASTUB,,,,,,LINE0-2,LINE0-1,LINE0,LINE1 !产生换热管的体单元ACLEAR,BASTUB !清除换热管下表面的圆环面的面单元管程筒体零度剖面WPROTA,,90 !旋转工作面NUMCMP,AREA !使所有面的面号从小到大重新排序*GET,NAREA1,AREA,,NUM,MAX !得到当前最大面号RECTANG,DSHELL/2.0D0+TSHELL-TCHANNEL,DSHELL/2.0D0+TSHELL,-(RCHANNEL+HPCH+HCHAN NEL),-(RCHANNEL+HPCH) !生成管程筒体零度剖面管板与壳程筒体和管程筒体连接区的零度剖面RECTANG,DSHELL/2.0D0+TSHELL-TCHANNEL,DSHELL/2.0D0+TSHELL,-(RCHANNEL+HPCH),-RCHAN NEL !生成管板管程侧凸肩零度剖面NUMCMP,LINE !使所有线的线号从小到大重新排序*GET,LINE2,LINE,,NUM,MAX !得到当前最大线号K,200001,DSHELL/2.0D0+TSHELL-TCHANNEL-RCHANNEL,0,-RCHANNEL !管板管程侧锻造圆角圆心坐标K,400001,DSHELL/2.0D0-RSHELL,,SHIGH+RSHELLK,300001,DSHELL/2.0D0+TSHELL-TCHANNEL,0,-RCHANNEL !管板管程侧锻造圆角弧线起点K,300002,DSHELL/2.0D0+TSHELL-TCHANNEL-RCHANNEL+RCHANNEL*COS(45),,-RCHANNEL+RCHA NNEL*COS(45)!管板管程侧锻造圆角弧线中点K,300003,DSHELL/2.0D0+TSHELL-TCHANNEL-RCHANNEL,,, !管板管程侧锻造圆角弧线终点K,300004,ROL,,, !形成此区零度剖面的关键点K,300005,ROL,,TCHANNEL !形成此区零度剖面的关键点K,300006,ROL,,SHIGH-TSHELL !形成此区零度剖面的关键点K,300007,ROL,,SHIGH !形成此区零度剖面的关键点K,300008,DSHELL/2.0D0-RSHELL,,SHIGH !形成此区零度剖面的关键点K,300009,DSHELL/2.0D0-RSHELL+RSHELL*COS(45),,SHIGH+(RSHELL-RSHELL*COS(45)) !形成此区零度剖面的关键点K,300010,DSHELL/2.0D0,,SHIGH+RSHELL !形成此区零度剖面的关键点K,300011,DSHELL/2.0D0+TSHELL,,SHIGH+RSHELL !形成此区零度剖面的关键点K,300012,DSHELL/2.0D0+TSHELL,,SHIGH-TSHELL !形成此区零度剖面的关键点K,300013,DSHELL/2.0D0+TSHELL,,TCHANNEL !形成此区零度剖面的关键点K,300014,DSHELL/2.0D0+TSHELL,,-RCHANNEL !形成此区零度剖面的关键点K,300015,DSHELL/2.0D0+TSHELL-TCHANNEL-RCHANNEL,,TCHANNEL !形成此区零度剖面的关键点K,300016,DSHELL/2.0D0-TSHELL,,SHIGH-TSHELL !形成此区零度剖面的关键点LARC,300001,300002,200001,RCHANNEL !生成此区零度剖面的外边缘线LINE2+1 LESIZE,LINE2+1,,,NR/2 !指定弧线剖分数L,300001,300014,NHOU !生成网格划分辅助线LINE2+2L,300002,300013,NHOU !生成网格划分辅助线LINE2+3L,300013,300014,NR/2 !生成此区零度剖面的外边缘线LINE2+4 AL,LINE2+1,LINE2+2,LINE2+3,LINE2+4 !生成此区零度剖面LARC,300002,300003,200001,RCHANNEL !生成此区零度剖面的外边缘线LINE2+5 LESIZE,LINE2+5,,,NR/2 !指定弧线剖分数L,300003,300015,NHOU !生成网格划分辅助线LINE2+6L,300013,300015,NR/2 !生成网格划分辅助线LINE2+7AL,LINE2+3,LINE2+5,LINE2+6,LINE2+7 !生成此区零度剖面L,300003,300004,1 !生成此区零度剖面的外边缘线LINE2+8 L,300004,300005,NHOU !生成网格划分辅助线LINE2+9L,300005,300006,NSHIGH !生成网格划分辅助线LINE2+10L,300006,300007,NHOU !生成网格划分辅助线LINE2+11L,300007,300008,1 !生成此区零度剖面的外边缘线LINE2+12 L,300008,300016,NHOU !生成网格划分辅助线LINE2+13L,300016,300015,NSHIGH !生成网格划分辅助线LINE2+14AL,LINE2+6,LINE2+8,LINE2+9,LINE2+10,LINE2+11,LINE2+12,LINE2+13,LINE2+14 !生成此区零度剖面L,300012,300016,NR/2 !生成网格划分辅助线LINE2+15L,300012,300013,NSHIGH !生成此区零度剖面的外边缘线LINE2+16 AL,LINE2+7,LINE2+14,LINE2+15,LINE2+16 !生成此区零度剖面L,300009,300012,NHOU !生成网格划分辅助线LINE2+17 LARC,300008,300009,400001,RSHELL !生成此区零度剖面的外边缘线LINE2+18 LESIZE,LINE2+18,,,NR/2 !指定弧线剖分数AL,LINE2+13,LINE2+15,LINE2+17,LINE2+18 !生成此区零度剖面LARC,300009,300010,400001,RSHELL !生成此区零度剖面的外边缘线LINE2+19 LESIZE,LINE2+19,,,NR/2 !指定弧线剖分数L,300010,300011,NHOU !生成网格划分辅助线LINE2+20L,300011,300012,NR/2 !生成此区零度剖面的外边缘线LINE2+21 AL,LINE2+17,LINE2+19,LINE2+20,LINE2+21 !生成此区零度剖面RECTANG,DSHELL/2.0D0,DSHELL/2.0D0+TSHELL,SHIGH+RSHELL,SHIGH+RSHELL+HPSH !生成管板管程侧凸肩零度剖面壳程筒体零度剖面RECTANG,DSHELL/2.0D0,DSHELL/2.0D0+TSHELL,SHIGH+RSHELL+HPSH,SHIGH+RSHELL+HPSH+HSH ELLLSEL,S,LOC,X,DSHELL/2.0D0+TSHELL-TCHANNEL/2.0D0 !选中线段LESIZE,ALL,,,NHOU !指定线段剖分数ALLSEL !全选LSEL,S,LOC,X,DSHELL/2.0D0+TSHELL/2.0D0 !选中线段LESIZE,ALL,,,NHOU !指定线段剖分数ALLSEL !全选保温层零度剖面RECTANG,DSHELL/2.0D0+TSHELL,DSHELL/2.0D0+TSHELL+TINSULAT,-RCHANNEL-HPCH-HCHANNEL ,-RCHANNEL-HPCHRECTANG,DSHELL/2.0D0+TSHELL,DSHELL/2.0D0+TSHELL+TINSULAT,-RCHANNEL-HPCH,-RCHANNE LRECTANG,DSHELL/2.0D0+TSHELL,DSHELL/2.0D0+TSHELL+TINSULAT,-RCHANNEL,TCHANNEL RECTANG,DSHELL/2.0D0+TSHELL,DSHELL/2.0D0+TSHELL+TINSULAT,TCHANNEL,SHIGH-TSHELL RECTANG,DSHELL/2.0D0+TSHELL,DSHELL/2.0D0+TSHELL+TINSULAT,SHIGH-TSHELL,SHIGH+RSHELL RECTANG,DSHELL/2.0D0+TSHELL,DSHELL/2.0D0+TSHELL+TINSULAT,SHIGH+RSHELL,SHIGH+RSHEL L+HPSHRECTANG,DSHELL/2.0D0+TSHELL,DSHELL/2.0D0+TSHELL+TINSULAT,SHIGH+RSHELL+HPSH,SHIGH+RSHELL+HPSH+HSHELLLSEL,S,LOC,Z,-(RCHANNEL+HPCH+HCHANNEL/2.0D0) !选中线段LESIZE,ALL,,,NCHANNEL !指定线段剖分数ALLSEL !全选LSEL,S,LOC,Z,-(RCHANNEL+HPCH/2.0D0) !选中线段LESIZE,ALL,,,NHPCH !指定线段剖分数ALLSEL !全选LSEL,S,LOC,Z,(-RCHANNEL+TCHANNEL)/2.0D0 !选中线段LESIZE,ALL,,,NR/2 !指定线段剖分数ALLSEL !全选LSEL,S,LOC,Z,TCHANNEL+(SHIGH-TCHANNEL-TSHELL)/2.0D0 !选中线段LESIZE,ALL,,,NHIGH !指定线段剖分数ALLSEL !全选LSEL,S,LOC,Z,SHIGH-TSHELL+(TSHELL+RSHELL)/2.0D0 !选中线段LESIZE,ALL,,,NR/2 !指定线段剖分数ALLSEL !全选LSEL,S,LOC,Z,SHIGH+RSHELL+HPCH/2.0D0 !选中线段LESIZE,ALL,,,NHPSH !指定线段剖分数ALLSEL !全选LSEL,S,LOC,Z,SHIGH+RSHELL+HPSH+HSHELL/2.0D0 !选中线段LESIZE,ALL,,,NHSHELL !指定线段剖分数ALLSEL !全选LSEL,S,LOC,X,DSHELL/2.0D0+TSHELL+TINSULAT/2.0D0 !选中线段LESIZE,ALL,,,NINSULAT !指定线段剖分数ALLSEL !全选管箱的体MSHAPE,0,2D !指定面单元形状为四边形TYPE,1 !指定面单元类型MAT,2 !设定材料属性AMESH,NAREA1+1 !对管箱零度剖面进行网格划分TYPE,2 !指定体单元类型ESIZE,,NR0 !定义单元尺寸VROTAT,NAREA1+1,,,,,,10001,10002,90,1 !产生管箱体单元管板的体MSHAPE,2,2D !指定面单元形状为四边形TYPE,1 !指定面单元类型MAT,3 !设定材料属性AMESH,NAREA1+2,NAREA1+9 !对管板零度剖面进行网格划分ESIZE,,NR0 !定义单元尺寸VROTAT,NAREA1+2,NAREA1+3,NAREA1+4,NAREA1+5,,,10001,10002,90,1 !产生管板体单元VROTAT,NAREA1+6,NAREA1+7,NAREA1+8,NAREA1+9,,,10001,10002,90,1 !产生管板体单元筒体的体MSHAPE,2,2D !指定面单元形状为四边形TYPE,1 !指定面单元类型MAT,1 !设定材料属性AMESH,NAREA1+10 !对筒体零度剖面进行网格划分ESIZE,,NR0 !定义单元尺寸VROTAT,NAREA1+10,,,,,,10001,10002,90,1 !产生筒体体单元保温层的体MSHAPE,2,2D !指定面单元形状为四边形TYPE,1 !指定面单元类型MAT,5 !设定材料属性AMESH,NAREA1+11,NAREA1+17 !对保温层零度剖面进行网格划分ESIZE,,NR0 !定义单元尺寸VROTAT,NAREA1+11,NAREA1+12,NAREA1+13,NAREA1+14,,,10001,10002,90,1 !产生保温层体单元VROTAT,NAREA1+15,NAREA1+16,NAREA1+17,,,,10001,10002,90,1 !产生保温层体单元ACLEAR,NAREA1+1,NAREA1+17,1 !清除所有零度剖面面单元NUMMRG,ALL !所有重合的点、线、面信息合并NUMCMP,ALL !所有点、线、面信息依次重排序FINISH! ***************求解***************/SOLULSEL,S,RADIUS,,DSHELL/2.0D0 !选线LSEL,A,RADIUS,,DSHELL/2.0D0-RSHELL !选线LSEL,A,RADIUS,,ROL !选线LSEL,U,LOC,Z,-RCHANNEL-HPCH-HCHANNEL,SHIGH-TSHELL !去掉多余的线ASLL,S,0 !选中与线相关的面ASEL,U,LOC,Z,SHIGH-TSHELL/2.0D0 !去掉多余的面ASEL,U,LOC,Z,SHIGH+RSHELL !去掉多余的面ASEL,U,LOC,Z,SHIGH+RSHELL+HPSH !去掉多余的面ASEL,U,LOC,Z,SHIGH+RSHELL+HPSH+HSHELL !去掉多余的面SFA,ALL,,CONV,ASHELL,TEPSHEL !对管板上表面施加温度边界条件ALLSEL !全选LSEL,S,RADIUS,,DP/2.0D0 !选线LSEL,R,LOC,Z,SHIGH+RSHELL+HPSH+HSHELL !选线ASLL,S,0 !选中与线相关的面ASEL,U,LOC,Z,SHIGH+RSHELL+HPSH+HSHELL !去掉多余的面SFA,ALL,,CONV,ASHELL,TEPSHEL !对换热管外表面施加温度边界条件ALLSEL !全选LSEL,S,RADIUS,,DP/2.0D0-TPIPE !选线ASLL,S,0 !选中与线相关的面LSEL,ALL !选中所有线段LSEL,S,RADIUS,,DP/2.0D0 !选线ASLL,U,0 !不选与线相关的面SFA,ALL,,CONV,ATUBE,TEPTUBE !对换热管内表面施加温度边界条件ALLSEL !全选LSEL,S,RADIUS,,DSHELL/2.0D0+TSHELL-TCHANNEL !选线LSEL,A,RADIUS,,DSHELL/2.0D0+TSHELL-TCHANNEL-RCHANNEL !选线LSEL,A,RADIUS,,ROL !选线LSEL,U,LOC,Z,TCHANNEL,SHIGH+RSHELL+HPSH+HSHELL !去掉多余的线ASLL,S,0 !选中与线相关的面ASEL,U,LOC,Z,TCHANNEL/2.0D0 !去掉多余的面ASEL,U,LOC,Z,-RCHANNEL !去掉多余的面ASEL,U,LOC,Z,-RCHANNEL-HPCH !去掉多余的面ASEL,U,LOC,Z,-RCHANNEL-HPCH-HCHANNEL !去掉多余的面SFA,ALL,,CONV,ACHANNEL,TEPTUBE !对管板下表面和管箱内表面施加温度边界条件ALLSEL !全选LSEL,S,RADIUS,,DP/2.0D0-TPIPE !选线LSEL,R,LOC,Z,0 !选线ASLL,S,0 !选中与线相关的面ASEL,U,LOC,Z,TCHANNEL/2.0D0 !去掉多余的面SFA,ALL,,CONV,ACHANNEL,TEPTUBE !对换热管下表面施加温度边界条件ALLSEL !全选LSEL,S,RADIUS,,DSHELL/2.0D0+TSHELL+TINSULAT !选线ASLL,S,0 !选中与线相关的面ASEL,U,LOC,Z,-RCHANNEL-HPCH-HCHANNEL !去掉多余的面ASEL,U,LOC,Z,-RCHANNEL-HPCH !去掉多余的面ASEL,U,LOC,Z,-RCHANNEL !去掉多余的面ASEL,U,LOC,Z,TCHANNEL !去掉多余的面ASEL,U,LOC,Z,SHIGH-TSHELL !去掉多余的面ASEL,U,LOC,Z,SHIGH+RSHELL !去掉多余的面ASEL,U,LOC,Z,SHIGH+RSHELL+HPSH !去掉多余的面ASEL,U,LOC,Z,SHIGH+RSHELL+HPSH+HSHELL !去掉多余的面SFA,ALL,,CONV,AENVIR,TEPENVIR !对保温层外表面施加温度边界条件ALLSEL !全选SOLVE !求解FINISH! ***************后处理***************/POST1PLNSOL,TEMP !绘出温度场分布云图! ***************前处理***************/PREP7ET,2,SOLID45 !热单元类型转换为结构单元类型MP,EX,1,2.1E5 !定义材料的弹性模量,1号材料对应筒体钢板MP,NUXY,1,0.3 !定义材料的泊松比MP,ALPX,1,1.272E-5 !定义热膨胀系数MP,EX,2,202780 !定义材料的弹性模量,2号材料对应管箱钢板MP,NUXY,2,0.3 !定义材料的泊松比MP,ALPX,2,1.75E-5 !定义热膨胀系数MP,EX,3,2E5 !定义材料的弹性模量,3号材料对应管板锻件MP,NUXY,3,0.3 !定义材料的泊松比MP,ALPX,3,1.655E-5 !定义热膨胀系数MP,EX,4,93710 !定义材料的弹性模量,4号材料对应换热管钢管MP,NUXY,4,0.3 !定义材料的泊松比MP,ALPX,4,1.618E-5 !定义热膨胀系数FINISH/SOLUVSEL,S,MAT,,1,4,1,1 !选中参加应力计算的材料NSEL,R,LOC,Z,SHIGH+RSHELL+HPSH+HSHELL !选点D,ALL,UZ,0 !约束ALLSEL !全选VSEL,S,MAT,,1,4,1,1 !选中参加应力计算的材料NSEL,R,LOC,X,0 !选点D,ALL,UX,0 !约束ALLSEL !全选VSEL,S,MAT,,1,4,1,1 !选中参加应力计算的材料NSEL,R,LOC,Y,0 !选点D,ALL,UY,0 !约束ALLSEL !全选LSEL,S,RADIUS,,DSHELL/2.0D0 !选线LSEL,A,RADIUS,,DSHELL/2.0D0-RSHELL !选线LSEL,A,RADIUS,,ROL !选线LSEL,U,LOC,Z,-RCHANNEL-HPCH-HCHANNEL,SHIGH-TSHELL !去掉多余的线ASLL,S,0 !选中与线相关的面ASEL,U,LOC,Z,SHIGH-TSHELL/2.0D0 !去掉多余的面ASEL,U,LOC,Z,SHIGH+RSHELL !去掉多余的面ASEL,U,LOC,Z,SHIGH+RSHELL+HPSH !去掉多余的面ASEL,U,LOC,Z,SHIGH+RSHELL+HPSH+HSHELL !去掉多余的面SFA,ALL,,PRES,PS !对管板上表面施加压力载荷ALLSEL !全选LSEL,S,RADIUS,,DP/2.0D0 !选线LSEL,R,LOC,Z,SHIGH+RSHELL+HPSH+HSHELL !选线ASLL,S,0 !选中与线相关的面ASEL,U,LOC,Z,SHIGH+RSHELL+HPSH+HSHELL !去掉多余的面SFA,ALL,,PRES,PS !对换热管外表面施加压力载荷ALLSEL !全选LSEL,S,RADIUS,,DP/2.0D0-TPIPE !选线ASLL,S,0 !选中与线相关的面LSEL,ALL !选中所有线段LSEL,S,RADIUS,,DP/2.0D0 !选线ASLL,U,0 !不选与线相关的面SFA,ALL,,PRES,PT !对换热管内表面施加压力载荷ALLSEL !全选LSEL,S,RADIUS,,DSHELL/2.0D0+TSHELL-TCHANNEL !选线LSEL,A,RADIUS,,DSHELL/2.0D0+TSHELL-TCHANNEL-RCHANNEL !选线LSEL,A,RADIUS,,ROL !选线LSEL,U,LOC,Z,TCHANNEL,SHIGH+RSHELL+HPSH+HSHELL !去掉多余的线ASLL,S,0 !选中与线相关的面ASEL,U,LOC,Z,TCHANNEL/2.0D0 !去掉多余的面ASEL,U,LOC,Z,-RCHANNEL !去掉多余的面ASEL,U,LOC,Z,-RCHANNEL-HPCH !去掉多余的面ASEL,U,LOC,Z,-RCHANNEL-HPCH-HCHANNEL !去掉多余的面SFA,ALL,,PRES,PT !对管板下表面和管箱内表面施加压力载荷ALLSEL !全选LSEL,S,RADIUS,,DP/2.0D0-TPIPE !选线LSEL,R,LOC,Z,0 !选线ASLL,S,0 !选中与线相关的面ASEL,U,LOC,Z,TCHANNEL/2.0D0 !去掉多余的面SFA,ALL,,PRES,PT !对换热管下表面施加压力载荷ALLSEL !全选*SET,PN,(DSHELL/2.0D0+TSHELL-TCHANNEL)**2*PT/((DSHELL/2.0D0+TSHELL)**2-(DSHELL/2.0+T SHELL-TCHANNEL)**2) !轴向平衡力VSEL,S,MAT,,1,4,1,1 !选中参加应力计算的材料ASEL,R,LOC,Z,-RCHANNEL-HPCH-HCHANNEL !选面SFA,ALL,,PRES,-PN !对面施加轴向平衡力ALLSELLDREAD,TEMP,,,,,,RTH !读入先前求解的温度场TREF,20 !定义参考温度VSEL,S,MAT,,1,4,1,1 !选中参加应力计算的材料SOLVE !求解FINISH! ***************后处理***************/POST1PLNSOL,S,INT !绘出应力强度云图! ***************应力评定路径1--1****************SET,FIRPOT,27934 !定义内壁节点号*SET,SECPOT,27932 !定义外壁节点号PATH,L1,2 !定义路径名称及路径参数设置PPATH,1,FIRPOT !定义路径起点PPATH,2,SECPOT !定义路径终点PRSECT, ,0 !列出路径线性化处理结果! ***************应力评定路径2--2****************SET,FIRPOT,27972 !定义内壁节点号*SET,SECPOT,27968 !定义外壁节点号PATH,L2,2 !定义路径名称及路径参数设置PPATH,1,FIRPOT !定义路径起点PPATH,2,SECPOT !定义路径终点PRSECT, ,0 !列出路径线性化处理结果! ***************应力评定路径3--3****************SET,FIRPOT,121 !定义内壁节点号*SET,SECPOT,7471 !定义外壁节点号PATH,L3,2 !定义路径名称及路径参数设置PPATH,1,FIRPOT !定义路径起点PPATH,2,SECPOT !定义路径终点PRSECT, ,0 !列出路径线性化处理结果FINISH。
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命令得到更详细的解释;后继载荷步在后继载荷步中,用户可以随意杀死或重新激活单元;象上面提到的,要正确的施加和删除约束和结点载荷;用下列命令杀死单元：Command:EKILLGUI: Main Menu>Solution>-Load Step Opts-Other>Kill Elements用下列命令重新激活单元：Command: EALIVEGUI: Main Menu>Solution>-Load Step Opts-Other>Activate Elem第二个或后继载荷步：TIME,...ESEL,...EKILL,... 杀死选择的单元ESEL,...EALIVE,... 重新激活选择的单元...FDELE,... 删除不活动自由度的结点载荷D,... 约束不活动自由度...F,... 在活动自由度上施加合适的结点载荷DDELE,... 删除重新激活的自由度上的约束SAVESOLVE四、u /grid, keykey: “0”或“off”无网络“1”或“on” xy网络“2”或“x”只有x线“3”或“y”只有y线u xvar, nn: “0”或“1”将x轴作为时间轴“n”将x轴表示变量“n”“-1”u /axlab, axis, lab 定义轴线的标志axis: “x”或“y”lab: 标志,可长达30个字符u plvar, nvar, nvar2, ……,nvar10 画出要显示的变量作为纵坐标五、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”仅用绝对值六、VDELE, NV1, NV2, NINC, KSWP: 删除未分网格的体nv1:初始体号nv2:最终的体号ninc:体号之间的间隔kswp=0:只删除体kswp=1:删除体及组成关键点,线面如果nv1=all,则nv2,ninc不起作用七、VSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWPType,是选择的方式,有选择s,补选a,不选,全选all、反选inv等,其余方式不常用Item, Comp 是选取的原则以及下面的子项如 volu 就是根据实体编号选择,loc 就是根据坐标选取,它的comp就可以是实体的某方向坐标其余还有材料类型、实常数等MIN, VMAX, VINC,这个就不必说了吧,例：vsel,s,volu,,14vsel,a,volu,,17,23,2上面的命令选中了实体编号为 14,17,19,21,23的五个实体u rforce, nvar, node, item, comp, name 指定待存储的节点力数据nvar: 变量号node: 节点号item compF x,M x, y,zname: 给此变量一个名称,8个字符u add, ir, ia,ib,ic,name,--,--,facta, factb, factc 将ia,ib,ic变量相加赋给ir变量ir, ia,ib,ic：变量号name: 变量的名称Fini退出四大模块,回到BEGIN层/cle 清空内存,开始新的计算1．定义参数、数组,并赋值.2． /prep7进入前处理定义几何图形：关键点、线、面、体定义几个所关心的节点,以备后处理时调用节点号;设材料线弹性、非线性特性设置单元类型及相应KEYOPT设置实常数设置网格划分,划分网格根据需要耦合某些节点自由度定义单元表存盘3．/solu加边界条件设置求解选项定义载荷步求解载荷步4./post1通用后处理5./post26 时间历程后处理菜单命令7.参数化设计语言8.理论手册Fini退出四大模块,回到BEGIN层/cle 清空内存,开始新的计算1 定义参数、数组,并赋值.u 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进入前处理定义几何图形：关键点、线、面、体u csys,kcnkcn , 0 迪卡尔zuobiaosi1 柱坐标2 球4 工作平面5 柱坐标系以Y轴为轴心n 已定义的局部坐标系u numstr, label, value 设置以下项目编号的开始nodeelemkplineareavolu注意：vclear, aclear, lclear, kclear 将自动设置节点、单元开始号为最高号,这时如需要自定义起始号,重发numstru K, npt, x,y,z, 定义关键点Npt：关键点号,如果赋0,则分配给最小号u Kgen,itime,Np1,Np2,Ninc,Dx,Dy,Dz,kinc,noelem,imoveItime：拷贝份数Np1,Np2,Ninc：所选关键点Dx,Dy,Dz：偏移坐标Kinc：每份之间节点号增量noelem: “0”如果附有节点及单元,则一起拷贝;“1”不拷贝节点和单元imove：“0”生成拷贝“1”移动原关键点至新位置,并保持号码,此时itime,kinc,noelem被忽略注意：MAT,REAL,TYPE 将一起拷贝,不是当前的MAT,REAL,TYPEu A, P1, P2, ……… P18 由关键点生成面u AL, L1,L2, ……,L10 由线生成面面的法向由L1按右手法则决定,如果L1为负号,则反向;线需在某一平面内坐标值固定的面内u vsba, nv, na, sep0,keep1,keep2 用面分体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：拷贝对象编号dx, dy, dz ：位移增量kinc: 对应关键点号增量noelem,：0：同时拷贝节点及单元1：不拷贝节点及单元imove： 0：拷贝体1：移动体u cm, cname, entity 定义组元,将几何元素分组形成组元cname: 由字母数字组成的组元名entity: 组元的类型volu, area, line, kp, elem, nodeu cmgrp, aname, cname1, ……,cname8 将组元分组形成组元集合aname: 组元集名称cname1……cname8: 已定义的组元或组元集名称u cmlist,nameu cmdele,nameu cmplot, label1定义几个所关心的节点,以备后处理时调用节点号;u n,node,x,y,z,thxy, thyz, thzx 根据坐标定义节点号如果已有此节点,则原节点被重新定义,一般为最大节点号;设材料线弹性、非线性特性u 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: 材料编号缺省为当前材料号c 材料特性值,或材料之特性,温度曲线中的常数项c1-c4: 材料的特性-温度曲线中1次项,2次项,3次项,4次项的系数u Tb, lab, mat, ntemp,npts,tbopt,eosopt 定义非线性材料特性表Lab: 材料特性表之种类Bkin: 双线性随动强化Bis 双线性等向强化Mkin: 多线性随动强化最多5个点Mis 多线性等向强化最多100个点Dp: dp模型Mat: 材料号Ntemp: 数据的温度数对于bkin: ntemp缺省为6mis ntemp缺省为1,最多20bis ntemp缺省为6,最多为6dp: ntemp, npts, tbopt 全用不上Npts: 对某一给定温度数据的点数u TBTEMP,temp,kmod 为材料表定义温度值temp: 温度值kmod: 缺省为定义一个新温度值如果是某一整数,则重新定义材料表中的温度值注意：此命令一发生,则后面的TBDATA和TBPT均指此温度,应该按升序若Kmod为crit, 且temp为空,则其后的tbdata数据为solid46,shell99,solid191中所述破坏准则如果kmod为strain,且temp为空,则其后tbdata数据为mkin中特性;u TBDATA, stloc, c1,c2,c3,c4,c5,c6给当前数据表定义数据配合tbtemp,及tb使用stloc: 所要输入数据在数据表中的初始位置,缺省为上一次的位置加1 每重新发生一次tb或tbtemp命令上一次位置重设为1,发生tb后第一次用空闲此项,则c1赋给第一个常数u tbpt, oper, x,y 在应力-应变曲线上定义一个点oper: defi 定义一个点dele 删除一个点x,y：坐标设置单元类型及相应KEYOPTu ET, itype, ename, kop1……kop6, inopr 设定当前单元类型Itype：单元号Ename：单元名设置实常数u Keyopt, itype, knum, valueitype: 已定义的单元类型号knum: 单元的关键字号value: 数值注意：如果 ,则必须使用keyopt命令,否则也可在ET命令中输入设置网格划分,划分网格映射网格划分1.面映射网格划分条件：a. 3或4条边b.面的对边必须划分为相同的单元或其划分与一个过渡形网格的划分相匹配c. 该面如有3条边,则划分的单元不必须为偶数,并且各边单元数相等d. mahkeye. mshpattern如果多于四条边,可将线合并成Lcomb可用amap命令,先选面,再选4个关键点即可指定面的对边的分割数,以生成过渡映射四边形网格,只适用于有四条边的面2. 体映射网格划分1若将体划分为六面体单元,必须满足以下条件a. 该体的外形为块状六面体、楔形或棱形五面体、四面体b. 对边必须划分为相同的单元数,或分割符合过渡网格形式c. 如果体是棱形或四面体,三角形面上的单元分割数必须是偶数2 当需要减少围成体的面数以进行映射网格划分时,可以对面相加或连接;如果连接而有边界线,线也必须连接在一起;3体扫掠生成网格步骤：a. 确定体的拓扑是否能够进行扫掠;侧面不能有孔；体内不能有封闭腔;源面与目标面必须相对b. 定义合适的单元类型c. 确定扫掠操作中如何控制生成单元层的数目 lesized. 确定体的哪一个边界面作为源面、目标面e. 有选择地对源面、目标面和边界面划分网格3. 关于连接线和面的一些说明连接仅是映射网格划分的辅助工具4. 用desize定义单元尺寸时单元划分应遵守的级别高：lesizekesizeesizedesize用smartzing定义单元尺寸时单元划分应遵守的级别高：lesizekesizesmartsizeu LESIZE,NL1,Size, Angsiz,ndiv,space,kforc,layer1,layer2,kyndiv 为线指定网格尺寸NL1: 线号,如果为all,则指定所有选中线的网格;Size: 单元边长,程序据size计算分割份数,自动取整到下一个整数Angsiz: 弧线时每单元跨过的度数Ndiv: 分割份数Space: “+”: 最后尺寸比最先尺寸“-“: 中间尺寸比两端尺寸free: 由其他项控制尺寸kforc 0: 仅设置未定义的线,1：设置所有选定线,2：仅改设置份数少的,3：仅改设置份数多的kyndiv: 0,No,off 表示不可改变指定尺寸1,yes,on 表示可改变u ESIZE,size,ndiv 指定线的缺省划分份数已直接定义的线,关键点网格划分设置不受影响u desize, minl, minh,……控制缺省的单元尺寸minl: n 每根线上低阶单元数缺省为3defa 缺省值stat 列出当前设置off 关闭缺省单元尺寸minh: n 每根线上高阶单元数缺省为2u mshape, key, dimension 指定单元形状key: 0 四边形2D,六面体3D1 三角形 2D, 四面体3DDimension: 2D 二维3D 三维u smart,off 关闭智能网格u mshkey, key 指定自由或映射网格方式key: 0 自由网格划分1 映射网格划分2 如果可能的话使用映射,否则自由即使自由smartsizing也不管用了u Amesh, nA1,nA2,ninc 划分面单元网格nA1,nA2,ninc 待划分的面号,nA1如果是All,则对所有选中面划分u SECTYPE, ID, TYPE, SUBTYPE, NAME, REFINEKEY定义一个截面号,并初步定义截面类型ID: 截面号TYPE: BEAM:定义此截面用于梁SUBTYPE: RECT 矩形CSOLID:圆形实心截面CTUBE: 圆管I: 工字形HREC: 矩形空管ASEC: 任意截面MESH: 用户定义的划分网格NAME: 8字符的截面名称字母和数字组成REFINEKEY: 网格细化程度：0~5对于薄壁构件用此控制,对于实心截面用SECDATA控制u SECDATA, VAL1, VAL2, …….VAL10 描述梁截面说明：对于SUBTYPE=MESH, 所需数据由SECWRITE产生,SECREAD读入u SECNUM,SECID 设定随后梁单元划分将要使用的截面编号u LATT, MAT, REAL, TYPE, --, KB, KE, SECNUM为准备划分的线定义一系列特性MAT: 材料号REAL: 实常数号TYPE: 线单元类型号KB、KE: 待划分线的定向关键点起始、终止号SECNUM: 截面类型号u SECPLOT,SECID,MESHKEY 画梁截面的几何形状及网格划分SECID:由SECTYPE命令分配的截面编号MESHKEY：0：不显示网格划分1：显示网格划分u /ESHAPE, SCALE 按看似固体化分的形式显示线、面单元SCALE: 0:简单显示线、面单元1：使用实常数显示单元形状u esurf, xnode, tlab, shape 在已存在的选中单元的自由表面覆盖产生单元xnode: 仅为产生surf151 或surf152单元时使用tlab: 仅用来生成接触元或目标元top 产生单元且法线方向与所覆盖的单元相同,仅对梁或壳有效,对实体单元无效Bottom产生单元且法线方向与所覆盖的单元相反,仅对梁或壳有效,对实体单元无效Reverse 将已产生单元反向Shape: 空与所覆盖单元形状相同Tri 产生三角形表面的目标元注意：选中的单元是由所选节点决定的,而不是选单元,如同将压力加在节点上而不是单元上u Nummrg,label,toler, Gtoler,actiontch 合并相同位置的itemlabel: 要合并的项目node: 节点, Elem,单元,kp: 关键点也合并线,面及点mat: 材料,type: 单元类型,Real: 实常数cp：耦合项,CE：约束项,CE: 约束方程,All：所有项toler: 公差Gtoler：实体公差Action: sele 仅选择不合并空合并注意：可以先选择一部分项目,再执行合并;如果多次发生合并命令,一定要先合并节点,再合并关键点;合并节点后,实体荷载不能转化到单元,此时可合并关键点解决问题;u Lsel, type, item, comp, vmin, vmax, vinc, kswp 选择线type: s 从全部线中选一组线r 从当前选中线中选一组线a 再选一部线附加给当前选中组aunoneuunselectinve: 反向选择item: line 线号loc 坐标length 线长comp: x,y,zkswp: 0 只选线1 选择线及相关关键点、节点和单元u Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点为下一步做准备Type: S: 选择一组新节点缺省R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0”使用正负号“1”仅用绝对值u NSLL,type, nkey 选择与所选线相联系的节点u nsla, type, nkey: 选择与选中面相关的节点type：s 选一套新节点r 从已选节点中再选a 附加一部分节点到已选节点u 从已选节点中去除一部分nkey: 0 仅选面内的节点1 选所有和面相联系的节点如面内线,关键点处的节点u esel, type, item, comp, vmin, vmax, vinc, kabs 选择一组单元Type: S: 选择一组单元缺省R: 在当前组中再选一部分作为一组A: 为当前组附加单元U: 在当前组中不选一部分单元All: 选所有单元None: 全不选Inve: 反向选择当前组Stat: 显示当前选择状态Item： Elem: 单元号Type: 单元类型号Mat: 材料号Real: 实常数号Esys: 单元坐标系号u ALLSEL, LABT, ENTITY 选中所有项目LABT: ALL: 选所有项目及其低级项目BELOW: 选指定项目的直接下属及更低级项目ENTITY: ALL: 所有项目缺省VOLU:体高级AREA:面LINE :线KP:关键点ELEM:单元NODE:节点低级u Tshap,shape 定义接触目标面为2D、3D的简单图形Shape: line:直线Arc:顺时针弧Tria:3点三角形Quad:4点四边形………….根据需要耦合某些节点自由度u cp, nset, lab,,node1,node2,……node17nset: 耦合组编号lab: ux,uy,uz,rotx,roty,rotznode1-node17: 待耦合的节点号;如果某一节点号为负,则此节点从该耦合组中删去;如果node1=all,则所有选中节点加入该耦合组;注意：1,不同自由度类型将生成不同编号2,不可将同一自由度用于多套耦合组u CPINTF, LAB, TOLER 将相邻节点的指定自由度定义为耦合自由度LAB：UX,UY,UZ,ROTX,ROTY,ROTZ,ALLTOLER: 公差,缺省为说明：先选中欲耦合节点,再执行此命令定义单元表说明：1,单元表仅对选中单元起作用,使用单元表之前务必选择一种类型的单元2,单元表各行为选中各单元,各列为每单元的不同数据u ETABLE, LAB, ITEM, COMP 定义单元表,添加、删除单元表某列LAB:用户指定的列名REFL, STAT, ERAS 为预定名称ITEM: 数据标志查各单元可输出项目COMP: 数据分量标志存盘u save, fname, ext,dir, slab 存盘fname : 文件名最多32个字符缺省为工作名ext: 扩展名最多32个字符缺省为dbdir: 目录名最多64个字符缺省为当前slab: “all”存所有信息“model”存模型信息“solv”存模型信息和求解信息3 /soluu /solu 进入求解器加边界条件u D, node, lab, value, value2, nend, ninc, lab2, lab3, ……lab6 定义节点位移约束Node : 预加位移约束的节点号,如果为all,则所有选中节点全加约束,此时忽略nend和ninc.Lab: ux,uy,uz,rotx,roty,rotz,allValue,value2: 自由度的数值缺省为0Nend, ninc: 节点范围为：node-nend,编号间隔为nincLab2-lab6: 将lab2-lab6以同样数值施加给所选节点;注意：在节点坐标系中讨论设置求解选项u antype, status, ldstep, substep, actionantype: static or 1 静力分析buckle or 2 屈曲分析modal or 3 模态分析trans or 4 瞬态分析status: new 重新分析缺省,以后各项将忽略rest 再分析,仅对static,full transion 有效ldstep: 指定从哪个荷载步开始继续分析,缺省为最大的,runn数指分析点的最后一步substep: 指定从哪个子步开始继续分析;缺省为本目录中,runn文件中最高的子步数action, continue: 继续分析指定的ldstep,substep说明：继续以前的分析因某种原因中断有两种类型singleframe restart: 从停止点继续需要文件：必须在初始求解后马上存盘单元矩阵或 .osav : 如果.esav坏了,将.osav改为.esavresults file: 不必要,但如果有,后继分析的结果也将很好地附加到它后面注意：如果初始分析生成了.rdb, .ldhi, 或rnnn 文件;必须删除再做后继分析步骤： 1进入anasys 以同样工作名2进入求解器,并恢复数据库3antype, rest4指定附加的荷载5指定是否使用现有的矩阵缺省重新生成kuse: 1 用现有矩阵6求解multiframe restart:从以有结果的任一步继续用不着u pred,sskey, --,lskey….. 在非线性分析中是否打开预测器sskey: off 不作预测当有旋转自由度时或使用solid65时缺省为offon 第一个子步后作预测除非有旋转自由度时或使用solid65时缺省为on -- ：未使用变量区lskey: off 跨越荷载步时不作预测缺省on 跨越荷载步时作预测此时sskey必须同时on注意：此命令的缺省值假定solcontrol为onu autots, key 是否使用自动时间步长key:on: 当solcontrol为on时缺省为onoff: 当solcontrol为off时缺省为off1: 由程序选择当solcontrol为on且不发生autots命令时在 .log文件中纪录“1”注意：当使用自动时间步长时,也会使用步长预测器和二分步长u NROPT, option,--,adptky 指定牛顿拉夫逊法求解的选项OPTION: AUT程序选择FULL:完全牛顿拉夫逊法MODI:修正的牛顿拉夫逊法INIT：使用初始刚阵UNSYM：完全牛顿拉夫逊法,且允许非对称刚阵ADPTKY:ON: 使用自适应下降因子OFF：不使用自适应下降因子u NLGEOM,KEYKEY: OFF:不包括几何非线性缺省ON：包括几何非线性u ncnv, kstop, dlim, itlim, etlim, cplim 终止分析选项kstop: 0 如果求解不收敛,也不终止分析1 如果求解不收敛,终止分析和程序缺省2如果求解不收敛,终止分析,但不终止程序dlim：最大位移限制,缺省为itlim: 累积迭代次数限制,缺省为无穷多etlim：程序执行时间秒限制,缺省为无穷cplim：cpu时间秒限制,缺省为无穷u solcontrol ,key1, key2,key3,vtol 指定是否使用一些非线性求解缺省值key1: on 激活一些优化缺省值缺省CNVTOL Toler=%Minref=对力和弯矩NEQIT 最大迭代次数根据模型设定在15~26之间ARCLEN 如用弧长法则用较更先进的方法PRED 除非有rotx,y,z或solid65,否则打开LNSRCH 当有接触时自动打开CUTCONTROL Plslimit=15%, npoint=13SSTIF 当NLGEOM,on时则打开NROPT,adaptkey 关闭除非：摩擦接触存在；单元12,26,48,49,52存在；当塑性存在且有单元20,23,24,60存在AUTOS 由程序选择off 不使用这些缺省值key2: on 检查接触状态此时key1为on此时时间步会以单元的接触状态据keyopt7的假定为基础当keyopt2=on 时,保证时间步足够小key3: 应力荷载刚化控制,尽量使用缺省值空：缺省,对某些单元包括应力荷载刚化,对某些不包括查nopl:对任何单元不包括应力刚化incp:对某些单元包括应力荷载刚化查vtol：u outres, item, freq, cname 规定写入数据库的求解信息item: all 所有求解项basic 只写nsol, rsol, nload, strsnsol 节点自由度rsol 节点作用荷载nload 节点荷载和输入的应变荷载strs 节点应力freq: 如果为n,则每n步包括最后一步写入一次none: 则在此荷载步中不写次项all: 每一步都写last: 只写最后一步静力或瞬态时为缺省定义载荷步u nsubst, nsbstp, nsbmx, nsbmn, carry 指定此荷载步的子步数nsbstp: 此荷载步的子步数如果自动时间步长使用autots,则此数定义第一子步的长度；如果solcontrol打开,且3D面-面接触单元使用,则缺省为1-20步；如果solcontrol打开,并无3D接触单元,则缺省为1子步；如果solcontrol关闭,则缺省为以前指定值；如以前未指定,则缺省为1nsbmx, nsbmn：最多,最少子步数如果自动时间步长打开u time, time 指定荷载步结束时间注意：第一步结束时间不可为“0”u f, node, lab, value, value2, nend, ninc 在指定节点加集中荷载node:节点号lab: Fx,Fy,Fz,Mx,My,Mzvalue: 力大小value2: 力的第二个大小如果有复数荷载nend,ninc：在从node到nend的节点增量为ninc上施加同样的力注意：1节点力在节点坐标系中定义,其正负与节点坐标轴正向一致u sfa, area, lkey, lab, value, value2 在指定面上加荷载area: n 面号all 所有选中号lkey: 如果是体的面,忽略此项lab: presvalue: 压力值u SFBEAM, ELEM, LKEY, LAB, VALI, VALJ, VAL2I, VAL2J, IOFFST, JOFFST 对梁单元施加线荷载ELEM: 单元号,可以为ALL,即选中单元LKEY: 面载类型号,见单元介绍;对于BEAM188,1为竖向；2为横向；3为切向VALI,VALJ: I, J节点处压力值VAL2I,VAL2J: 暂时无用IOFFST, JOFFST: 线载距离I, J 节点距离u lswrite, lsnum 将荷载与荷载选项写入荷载文件中lsnum ：荷载步文件名的后缀,即荷载步数当 stat 列示当前步数init 重设为“1”缺省为当前步数加“1”注意1. 尽量加面载,不加集中力,以免奇异点2. 面的切向荷载必须借助面单元求解载荷步u lssolve, lsmin, lsmax, lsinc 读入并求解多个荷载步lsmin, lsmax, lsinc ：荷载步文件范围4 /post1通用后处理u set, lstep, sbstep, fact, king, time, angle, nset 设定从结果文件读入的数据lstep :荷载步数sbstep：子步数,缺省为最后一步time：时间点如果弧长法则不用nset： data set numberu dscale, wn, dmult 显示变形比例wn: 窗口号或all,缺省为1dmult, 0或auto : 自动将最大变形图画为构件长的5%u pldisp, kund 显示变形的结构kund： 0 仅显示变形后的结构1 显示变形前和变形后的结构2 显示变形结构和未变形结构的边缘u get, par, node, n, u, xy,z 获得节点n的xy,z位移给参数par等价于函数 ux,uy,uznodex,y,z: 获得x,y,z节点号arnodex,y,z：获得和节点n相连的面注意：此命令也可用于/solu模块u fsum, lab, item 对单元之节点力和力矩求和lab: 空在整体迪卡尔坐标系下求和rsys 在当前激活的rsys坐标系下求和item: 空对所有选中单元不包括接触元求和cont: 仅对接触节点求和u PRSSOL, ITEM, COMP 打印BEAM188、BEAM189截面结果说明：只有刚计算完还未退出ANSYS时可用,重新进入ANSYS时不可用item comp 截面数据及分量标志S COMP X,XZ,YZ应力分量PRIN S1,S2,S3主应力SINT应力强度,SEQV等效应力EPTO COMP 总应变PRIN 总主应变,应变强度,等效应变EPPL COMP 塑性应变分量PRIN 主塑性应变,塑性应变强度,等效塑性应变u plnsol, item, comp, kund, fact 画节点结果为连续的轮廓线item: 项目见下表comp: 分量kund: 0 不显示未变形的结构1 变形和未变形重叠2 变形轮廓和未变形边缘fact: 对于接触的2D显示的比例系数,缺省为1item comp discriptionu x,y,z,sum 位移rot x,y,z,sum 转角s x,y,z,xy,yz,xz 应力分量1,2,3 主应力Int,eqv 应力intensity,等效应力epeo x,y,z,xy,yz,xz 总位移分量1,2,3 主应变Int,eqv 应变intensity,等效应变epel x,y,z,xy,yz,xz 弹性应变分量1,2,3 弹性主应变Int,eqv 弹性intensity,弹性等效应变eppl x,y,z,xy,yz,xz 塑性应变分量u PRNSOL, item, comp 打印选中节点结果item: 项目见上表comp: 分量u PRETAB, LAB1, LAB2, ……LAB9 沿线单元长度方向绘单元表数据LABn : 空：所有ETABLE命令指定的列名列名：任何ETABLE命令指定的列名u PLLS, LABI, LABJ, FACT, KUND 沿线单元长度方向绘单元表数据LABI:节点I的单元表列名LABJ:节点J的单元表列名FACT: 显示比例,缺省为1kund: 0 不显示未变形的结构1 变形和未变形重叠2 变形轮廓和未变形边缘5 /post26 时间历程后处理u nsol, nvar, node, item, comp,name在时间历程后处理器中定义节点变量的序号nvar：变量号从2到nv根据numvar定义node: 节点号item compu x, y,zrot x, y,zu ESOL, NVAR, ELEM, NODE, ITEM, COMP, NAME 将结果存入变量NVAR: 变量号,2以上ELEM: 单元号NODE: 该单元的节点号,决定存储该单元的哪个量,如果空,则给出平均值ITEM:COMP:NAME: 8字符的变量名, 缺省为ITEM加COMPu rforce, nvar, node, item, comp, name 指定待存储的节点力数据nvar: 变量号node: 节点号item compF x,M x, y,zname: 给此变量一个名称,8个字符u add, ir, ia,ib,ic,name,--,--,facta, factb, factc将ia,ib,ic变量相加赋给ir变量ir, ia,ib,ic：变量号name: 变量的名称u /grid, keykey: “0”或“off”无网络“1”或“on” xy网络“2”或“x”只有x线“3”或“y”只有y线u xvar, nn: “0”或“1”将x轴作为时间轴“n”将x轴表示变量“n”“-1”u /axlab, axis, lab 定义轴线的标志axis: “x”或“y”lab: 标志,可长达30个字符u plvar, nvar, nvar2, ……,nvar10 画出要显示的变量作为纵坐标u prvar, nvar1, ……,nvar6 列出要显示的变量6 PLOTCONTROL菜单命令u pbc, ilem, ……,key, min, max, abs 在显示屏上显示符号及数值item: u 所加的位移约束rot 所加的转角约束key: 0 不显示符号1 显示符号2 显示符号及数值u /SHOW, FNAME, EXT, VECT, NCPL 确定图形显示的设备及其他参数FNAME: X11:屏幕文件名：各图形将生成一系列图形文件JPEG: 各图形将生成一系列JPEG图形文件说明：没必要用此命令,需要的图形文件可计算后再输出7 参数化设计语言u do, par, ival, fval, inc 定义一个do循环的开始par: 循环控制变量ival, fval, inc：起始值,终值,步长正,负u enddo 定义一个do循环的结束u if,val1, oper, val2, base: 条件语句val1, val2: 待比较的值也可是字符,用引号括起来oper: 逻辑操作当实数比较时,误差为1e-10eq, ne, lt, gt, le, ge, ablt, abgtbase: 当oper结果为逻辑真时的行为lable: 用户定义的行标志stop: 将跳出anasysexit: 跳出当前的do循环cycle: 跳至当前do循环的末尾then: 构成if-then-else结构。
ANSYS热分析指南——ANSYS稳态热分析word精品文档59页
ANSYS热分析指南(第三章)第三章稳态热分析3.1稳态传热的定义ANSYS/Multiphysics，ANSYS/Mechanical，ANSYS/FLOTRAN和ANSYS/Professional这些产品支持稳态热分析。
稳态传热用于分析稳定的热载荷对系统或部件的影响。
通常在进行瞬态热分析以前，进行稳态热分析用于确定初始温度分布。
也可以在所有瞬态效应消失后，将稳态热分析作为瞬态热分析的最后一步进行分析。
稳态热分析可以计算确定由于不随时间变化的热载荷引起的温度、热梯度、热流率、热流密度等参数。
这些热载荷包括：对流辐射热流率热流密度（单位面积热流）热生成率（单位体积热流）固定温度的边界条件稳态热分析可用于材料属性固定不变的线性问题和材料性质随温度变化的非线性问题。
事实上，大多数材料的热性能都随温度变化，因此在通常情况下，热分析都是非线性的。
当然，如果在分析中考虑辐射，则分析也是非线性的。
3.2热分析的单元ANSYS和ANSYS/Professional中大约有40种单元有助于进行稳态分析。
有关单元的详细描述请参考《ANSYS Element Reference》，该手册以单元编号来讲述单元，第一个单元是LINK1。
单元名采用大写，所有的单元都可用于稳态和瞬态热分析。
其中SOLID70单元还具有补偿在恒定速度场下由于传质导致的热流的功能。
这些热分析单元如下：表3-1二维实体单元表3-2三维实体单元表3-3辐射连接单元表3-4传导杆单元表3-5对流连接单元表3-6壳单元表3-7耦合场单元表3-8特殊单元3.3热分析的基本过程ANSYS热分析包含如下三个主要步骤：前处理：建模求解：施加荷载并求解后处理：查看结果以下的内容将讲述如何执行上面的步骤。
首先，对每一步的任务进行总体的介绍，然后通过一个管接处的稳态热分析的实例来引导读者如何按照GUI路径逐步完成一个稳态热分析。
最后，本章提供了该实例等效的命令流文件。
8-2传热及温度应力分析ANSYS算例
(13) 定义材料参数 Main Menu → Preprocessor → Material Props → Material Models → Material Models Available: Structural（双击打开子菜单） → Linear（双击） → Elastic （双击）→ Isotropic（双击） → EX： 2.0e5 (弹性模量) ，PRXY：0.3 (泊松比)→ OK →转到Material Models Available: Thermal Expansion（双击） →Secant Coefficient （双击）→ Isotropic（双击） →ALPX：1.2E-5(平均线膨胀系数) → OK →关闭材料 定义菜单（点击菜单的右上角X）
(3) 设置计算类型 Main Menu: Preferences… → select Thermal, steady → OK
(4) 选择单元类型 Main Menu: Preprocessor → Element Type → Add/Edit/Delete → Add → Thermal Solid，
【ANSYS 应用实例 2.1】 焊接接头稳态传热过程的数值模拟
如图 2-1 所示，圆形的冷凝管通过法兰接头进行对接。接头的制作方法如下：先把法兰 移动到圆管接头位置，然后沿圆周焊接两道次，把法兰连接到圆管上。用螺栓把两个法兰接 头拉紧，法兰之间压上一块垫片。圆管内的液体温度为 0℃，蒸汽冷凝在圆管的外表面上， 蒸汽温度为 100℃。圆管内表面换热系数为 5000W/m2K，外表面换热系数为 20000 W/m2K。
UNIT2-2
TH-FEA(应用实例-UNIT2)
清华大学 曾攀
(6) 生成几何模型 Main Menu: Preprocessor → Modeling → Create → Keypoints → In Active CS → NPT
一个玻璃钢外壳ANSYS热应力分析的命令流
一个玻璃钢外壳ANSYS热应力分析的命令流/filname,test1 ！定义分析名称/title,this is a plane thermal analysis abort test1 ！定义标题/prep7 ！进入前处理et,1,plane55 ！设置单元类型mp,kxx,1,0.4 ！定义玻璃钢的导热系数k,1,0,0.021,0 ！创建关键点k,2,0.0793,0.021,0k,3,0.0793,0.027,0k,4,0.0553,0.027,0k,5,0.0553,0.0317,0k,6,0.0647,0.041,0k,7,0.0793,0.041,0k,8,0.0793,0.045,0k,9,0,0.045,0k,10,0,0,0k,11,0.1,0,0a,1,2,3,4,5,6,7,8,9, ！由关键点生成面aesize,1,0.002 ！设置网络划分数量amesh,all ！划分网络/solu ！进入求解器antype,static ！设置为稳态分析dl,8,1,temp,175,0 ！给编号为8的线上施加175度的栽荷lsel,s,loc,y,0,0.044 ！选择除8号线外其它的所有的线（通过坐标）dl,all,1,temp,30,0 ！施加30度的边界条件allsel,all ！选择所有outpr,basic,all ！控制输出文件类型solve ！求解finish/prep7etchg,tts ！转换单无，由热——结构mp,ex,1,2.1e10 ！定义弹性模量mp,prxy,1,0.3 ！定义泊松比mp,alpx,1,1.3e-5 ！定义玻璃钢的线膨胀系数/solu ！进入求解器antype,staticlsel,s,loc,x,0.02,0.079lsel,r,loc,y,0,0.044dl,all,,ux,0dl,all,,uy,0allsel,alldl,2,1,symmdl,7,1,symmdl,9,1,symmsfl,8,pres,140,0 tunif,30 ldread,temp,,,,,,rth solve finish。
Workbench教程_流体热单元和热应力分析
Workbench-Simulation Heat Transfer 11.0Workshop 10流体热单元和热应力分析分析Workshop Supplement •一个有内部热源的物块(4m x 10m x 1m)，通过三个水流通道(直径0.25m)进行冷却。
•计算流体的最大温度和全部的热流•通道由一维流体热单元“Pipe” 进行建模水流通道以线框方式查看分析参数Workshop Supplement•材料属性ρhVD =–使用钢材料Re启动几何文件Workshop Supplement•双击几何文件passages.agdb•查看实体模型。
附加三条线，用来表示水流。
–Solid–Center,Right,Left•查看参数管理器中的内容查看参数管中的内容•返回到Project页，然后启动一个新的仿真分析。
高亮仿真Workshop Supplement使用米作为长度单位指定水的材料属性Workshop Supplement •在Geometry下，高亮“Center” ，找到默认设置—Structural Steel。
•在水的细节栏中，左击Structural Steel然后选择新材料…•右击New Material，输入值并重命名为Water。
N M t i l输入值并重命名为W t设置材料属性Workshop Supplement •类似的，设置其他两条水流通道的材料属性…–在Right和Left下，将Structural Steel替代为S l S lWater。
钢的材料属性Workshop Supplement •接受Simulation定义的钢材料属性。
料属性分析类型Workshop Supplement •New Analysis –Steady State Thermal水流的几何–命令行Workshop Supplement•单元类型:默认情况下Si l ti将线定义为梁我们使用命令行来覆盖这个默认–默认情况下，Simulation将线定义为梁。
(完整版)史上最全的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分析实例与工程应用命令流学习笔记1大纲静力分析：2杆、3梁、5薄膜和板壳、4实体单元梁单元：简化计算，结构总体受力情况实体单元：较复杂的结构，局部细节的受力情况稳定性分析：6振动、模态分析：7简单振动和梁的振动、8膜板和实体振动2杆系结构的静力分析2.1铰接杆在外力作用下的变形二维杆单元LINK1*AFUN,DEG：三角函数默认为弧度，改为角度后处理：结构变形图、显示节点位移和杆件应力2.2人字形屋架的静力分析后处理：杆单元的轴力、轴向应力、轴向应变2.3超静定拉压杆的反力计算后处理：节点反力2.4平行杆件与刚性梁连接的热应力问题定义3点的UY为耦合自由度，即三者的UY位移相等温度（增量）后处理：寻找特定位置的节点和单元，并从单元表中提取它们的内力2.5端部有间隙的杆的热膨胀二维带厚度的平面应力单元PLANE42、二维接触单元CONTACT26温度（始、末）后处理：定义水平应力和铅直应力单元表，并提取3号单元的应力结果*Status,ParmFINISH定义数组变量，将计算结果通过数组变量输出到文件3梁的弯曲静力分析3.1单跨等截面超静定梁的平面弯曲二维弹性梁单元BEAM3后处理：定义以两端弯矩和剪力的单元表，并列出单元表数据并用单元表数据绘制剪力图和弯矩图更细的节点划分方案，更精细3.2四跨连续梁的内力计算体素建模：keypoint, line, area, volume便于细分单元3.3七层框架结构计算3.4工字形截面外伸梁的平面弯曲3.5矩形截面梁的纵横弯曲分析考虑应力强化效应后处理：迭代过程3.6空间刚架静力分析三维梁单元BEAM43.7悬臂梁的双向弯曲三维8节点耦合场实体单元SOLID5三维20节点固体单元SOLID92三维10节点耦合场实体单元SOLID98三维结构实体自适应单元SOLID147定义宏程序，对应四种工况，各种结果差别不大3.8圆形截面悬臂杆的弯扭组合变形三维直管单元PIPE16（只定义外直径，不定义内直径）3.9悬臂等强度梁的弯曲四边形壳单元SHELL63（这里用退化的三角形单元，并使用节点耦合自由度保证模型的对称变形）三维非对称锥形梁单元BEAM44（定义横截面主轴，单元宽度线性变化）计算结果都很好，但壳体单元更能模拟出等强度梁的实际几何形状，更直观，截面定义更简单。
ANSYS命令流学习笔记6-热应力分析
!ANSYS命令流学习笔记6!热应力分析!学习重点：!1、理解热力耦合的直接法和间接法!间接法：先进行热分析，然后将求得的节点温度作为载荷施加到结构应力分析中。
!直接法：直接采用具有温度和位移自由度的耦合单元，同时得到热分析和结构应力分析的结果。
直接法又分弱耦合和强耦合选择强耦合时，形成不对称矩阵，线性系统可以直接求解。
选择弱耦合时，对称矩阵，还是把热和结构分别进行求解，并将热结果施加在结构上，是间接法的变形，至少经过两次迭代。
弱耦合可以保证精度。
!2、如何利用坐标值来选择单元或几何。
熟练应用nsel，lsel，asel命令。
选择不同的单元，指定不同单元类型，或者材料属性!3、后处理强度理论的理解。
不同的材料可以发生不同形式的失效。
而且同一种材料在不同的受力状态下，也可以发生不同的失效模式。
如碳钢单向拉伸，以屈服模式失效。
但制成螺钉时，其根部应力集中引起三向拉伸，会出现断裂。
铸铁单向拉伸断裂失效，但是钢球挤压铸铁板时，接触点三向受压状态，铸铁出现屈服。
无论脆性还是塑性材料，在三向拉应力相近时应用第一强度理论（最大拉应力），以断裂失效判定。
在三向压应力相近时，都会引起塑性变形，采用第三或第四强度理论。
！第三强度理论，最大切应力理论。
各向同性的材料，最大剪应力校核，适用于塑性材料，屈服失效。
偏保守。
σ1-σ3≤[σ]。
莫尔强度理论可以看做第三强度理论的推广，但是实际上莫尔强度理论以试验资料为基础，经过逻辑综合得到的。
！第四强度理论，最大形状改变比能理论，适用于塑性材料的屈服失效，比第三理论适用范围广。
Squa{1/2*[ (σ1-σ2)^2 + (σ2-σ3)^2 +(σ3-σ1)^2 ] } ≤[σ]!案例如下：!某液体管路内部通有液体，外部包有保温层，保温层与空气接触，结构如图2.1所示。
已知管路由铸铁制造，其导热系数为70W/(m·℃)，弹性模量为200GPa，泊松比为0.3，热膨胀系数为1.2×10-5/℃；保温层的导热系数为0.02W/(m·℃)，弹性模量为20GPa，泊松比为0.4，热膨胀系数为1.2×10-5/℃；管路内液体压力0.3MPa，温度为70℃，对流换热系数为1W/(m2·℃)；空气温度为-40℃，对流换热系数为0.5W/(m2·℃)。
ANSYS应力场命令流
ANSYS应力场命令流FINISH/SOLUANTYPE,,REST解除死单元的约束/////////////////////////////////////////////////ESEL,S,LIVE !选择所有活动单元NSLE,S !选择所有活动结点NSEL,INVEDDELE,ALL,ALL 解除所有节点的自由度约束ALLSEL,ALL第一条焊缝开始分析0.4s结束第一条焊缝开始分析0.8s/prep7ALLSEL,ALL!激活所有的单元EALIVE,all!杀死所有的焊缝上的单元VSEL,S,,,1,2,1 选择第1-2条焊缝ESLV,S 选择第1-2条焊缝上的单元Ekill,all 杀死第1-2条焊缝上的单元ALLSEL,ALL!激活焊点后面的焊缝单元，现在主要是对第1条焊缝上的单元VSEL,S,,,1 选择第1条焊缝NSLV,R,1 再选择第1条焊缝上所有的节点NSEL,R,LOC,X,0,0.001ESLN,R 再选择0-0.4s间的单元，不包括下一个焊缝的单元EALIVE,all !激活焊点后面的焊缝单元ALLSEL,ALL!读入温度场结果，杀死所有的超过800摄氏度的单元/////////////////////////ETCHG,STT/POST1VSEL,S,,,1, 选择第1条焊缝ESLV,S 选择第1条焊缝上的单元SET,3,LAST,1NSEL,S,TEMP,,800,3000,,0ESLN,R,1/prep7Ekill,all 杀死800-3000摄氏度的单元ALLSEL,ALLETCHG,TTS 转换成结构单元再死掉表面效应单元ESEL,S,TYPE,,2Ekill,all 杀死第1-2条焊缝上的表面效应单元ALLSEL,ALLTREF,26LDREAD,TEMP,,,0.8,,hanjiejianmo,rth,!/SOLU约束所有不活动的结点自由度（可选ESEL,S,LIVE !选择所有活动单元NSLE,S !选择所有活动结点NSEL,INVE !选择所有非活动结点（不与活动单元相连的结点）D,all,,,,,,ALL,,,,, !约束所有不活动的结点自由度（可选）ALLSEL,ALL////////////////约束焊件的端部DL,1,,UX,DL,1,,UZ,DL,3,,UZ,D,7426,,,,,,UY,,,,,D,7930,,,,,,UY,,,,,分析类型选项trnopt,full 指定瞬态分析选项timint,on 时间积分效应nlgeom,on 大变形选项SSTIF, 1 硬化线性分析，在大变形选项打开时使用nropt,full,,on !在完全瞬态分析中,指定Newton－raphson选项EQSLV,,,0,///////////////////////////////////网上的该法CUTCONTROL,PLSLIMIT,2.5,CUTCONTROL,CRPLIMITexp,0.2,0CUTCONTROL,CRPLIMITimp,0.25,1CUTCONTROL,NPOINT,15////////////////////////////////////载荷步选项TIME,0.8NSUBST,20,KBC,0autots,on !指定使用自动步长跟踪非线性选项LNSRCH,OFFLNSRCH,1pred,on,,onCNVTOL,F,,0.005,2,, /////////////////////////////增大收敛值NEQIT,40 /////////////////////////////增大迭代次数TOFFST,273outres,all,last 规定写入数据库的求解信息/PBF,TEMP,,1SOLVE FINISH。
ANSYS热应力分析实例
A N S Y S热应力分析实例-CAL-FENGHAI-(2020YEAR-YICAI)_JINGBIAN热流体在代有冷却栅的管道里流动，如图为其轴对称截面图。
管道及冷却栅的材料均为不锈钢，导热系数为1.25Btu/hr-in-oF，弹性模量为28E6lb/in2泊松比为0.3。
管内压力为1000 lb/in2，管内流体温度为450 oF，对流系数为1 Btu/hr-in2-oF,外界流体温度为70 oF，对流系数为0.25 Btu/hr-in2-oF。
求温度及应力分布。
7.3.2菜单操作过程7.3.2.1设置分析标题1、选择“Utility Menu>File>Change Title”，输入Indirect thermal-stress Analysis of a cooling fin。
2、选择“Utility Menu>File>Change Filename”，输入PIPE_FIN。
7.3.2.2进入热分析，定义热单元和热材料属性1、选择“Main Menu>Preprocessor>Element Type>Add/Edit/Delete”，选择PLANE55，设定单元选项为轴对称。
2、设定导热系数：选择“Main Menu>Preprocessor>MaterialPorps>Material Models”，点击Thermal，Conductivity，Isotropic，输入1.25。
7.3.2.3创建模型1、创建八个关键点，选择“MainMenu>Preprocessor>Creat>Keypoints>On Active CS”，关键点的坐标如下：编号 1 2 3 4 5 6 7 8X 5 6 12 12 6 6 5 5Y 0 0 0 0.25 0.25 1 1 0.252、组成三个面：选择“MainMenu>Preprocessor>Creat>Area>Arbitrary>Throuth Kps”，由1,2,5,8组成面1；由2,3,4,5组成面2；由8,5,6,7组成面3。
ANSYS热应力分析命令流
/FILNAME,Double,1 ！定义工作文件名。
/TITLE,Temperature Analysis ！定义工作标题。
!*/PREP7!定义单元。
ET,1,SOLID70!*!定义材料属性。
MPTEMP,,,,,,,, !定义材料1。
MPTEMP,1,0MPDATA,KXX,1,,238*3.6 !定义材料1的传热系数KXX1。
MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,C,1,,500 !定义材料1的比热C1。
MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,1,,3e-6 !定义材料1密度DENS1。
!*MPTEMP,,,,,,,, !定义材料2。
MPTEMP,1,0MPDATA,KXX,2,,15*3.6 !定义材料2的传热系数KXX2。
MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,C,2,,100 !定义材料2的比热C2。
MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,2,,2.2e-6 !定义材料2密度DENS2。
!*!建立几何模型。
BLC4,-80,-10,160,20,700VOFFST,3,20, ,!*!网格划分。
FLST,5,20,4,ORDE,2FITEM,5,1FITEM,5,-20CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINECMSEL,,_YLESIZE,_Y1,10, , , , , , ,1 !定义网格大小。
!*TYPE, 1MAT, 1REAL,ESYS, 0SECNUM,CM,_Y,VOLUVSEL, , , , 1CM,_Y1,VOLUCHKMSH,'VOLU'CMSEL,S,_YVSWEEP,_Y1CMDELE,_YCMDELE,_Y1CMDELE,_Y2TYPE, 1MAT, 2REAL,ESYS, 0SECNUM,CM,_Y,VOLUVSEL, , , , 2CM,_Y1,VOLUCHKMSH,'VOLU'CMSEL,S,_YVSWEEP,_Y1CMDELE,_YCMDELE,_Y1CMDELE,_Y2!定义网格大小完成。
基于ANSYS的某刚性梁的热应力分析
Ｎ， ，１ 一１ ６ ，
图 １外 力 合 温 度 作 用 下 的 平 行 杆 刚 架模 型
ＦＩ工 Ｉ Ｅ，１ ４ ， １ ． 元 类 型 、 材 料 属 性 及 单 元 实 常 数 我 们 选 用 二 维 杆 单 元 ２单 ＬＮ ｔ 可 以 用 于 许 多 工 程 领 域 ， 以 将 其 看 作 是 桁 架 、 杆 或 者 弹 Ｉ Ｋ１ 它 ｎ 可 连 簧 等构 件 。二 维 单 元 是 带 有 两 个 节 点 ， 个 节点 有 ２个 沿 ｘ 和 Ｙ方 向 每
【 摘 要 】 本文用 Ａ ＹＳ软件对一平行杆件 与刚性梁连接 的热应 力问题进行 了有限元分析 ， 算了其 热应力 ， ＮＳ 计 并与理论解进行 了对比 ， 其 Ａ ＹＳ解与 理 论 解 完全 一 致 。从 而 论 证 了用 有 限 元进 行 结 构 分 析 的优 越 性 。 ＮＳ 【 关键词 】 ＹＳ 有 限元 ； ＡＮＳ ； 刚性梁 ； 热应力 １ 刚性 梁 热 应 力 的 理 论 解 和 有 限 元解 。 １１问 题 概 述 长 度 ｆ１ 的 两 根 铜 杆 和一 根 铁 杆 ＿ 平 行 地 铰 结 ． ＝ｍ ｌ １ ， 在 刚 性 壁 和 刚 性 梁 之 间 （ 图 １所 示 ） 它 们 的 横 截 面 都 等 于 Ａ＝ 如 ， ６ ｍ ， 和 铁 的 弹 性 模 量 分 别 等 于 丘＝ ０ Ｇ ａ和 Ｅ＝ １ＧＰ ， 膨 胀 ５ ｍ 刚 １０ Ｐ 。２ ０ ａ热 系数 分 别 等 于 吼＝ ．ｘ ℃和 ＝ ．ｘ Ｏ １ ｌ ６ １ ｌ ℃。刚 性壁 上 铰 节 点 之 间 的 ３ 距离 为 １ 试计算在 刚性梁 中部 的集中力为 Ｑ ４ Ｏ Ｏ ｍ， ＝ Ｏ Ｏ Ｎ作用下 ， 温度 升 高 Ｖ ＝ ０Ｃ 后 ， 杆 件 中 的应 力 。 Ｔ １ ｏ以 ３根 在 两 个 相 距 ／ｌ 的 刚 性 面 之 间 有 一 根 等 截 面 杆 ＿ 杆 件 材 料 的 弹 ＝ｍ ｌ １ ， 性 模 量 Ｅ ２ ０ Ｐ ， 在距 离 上端 ａ ０３ 和 ０６ 截 面 位 置 处 分 别 为受 ＝ １Ｇ ａ ＝． ｍ ．ｍ 到沿 杆 件 轴 向 的 两 个 集 中力 Ｆ＝ ０ ０ 。１０ Ｎ和 Ｆ＝ ０ Ｎ （ 图 １所 示 ）本 ２５０ 如 ， 文利 用 二 维 杆 单 元 ＵＮ 确 定 两 个 刚性 面 对 杆 件 的支 反 力 Ｒ 和 Ｒ 。 Ｋ１ 。 ２
Ansys 第36例 热应力分析(间接法)实例—液体管路doc资料
A n s y s第36例热应力分析(间接法)实例—液体管路第36例热应力分析（间接法）实例—液体管路本例介绍了利用间接法进行热应力计算的方法和步骤：首先进行热分析得到结构节点温度分布，然后把温度作为载荷施加到结构上并进行结构分析。
36.1概述利用间接法计算热应力，首先进行热分析，然后进行结构分析。
热分析可以是瞬态的，也可以是稳态的，需要将热分析求得的节点温度作为体载荷施加到结构上。
当热分析是瞬态的时，需要找到温度梯度最大的时间点，并将该时间点的结构温度场作为体载荷施加到结构上。
由于间接法可以使用所有热分析和结构分析的功能，所以对于大多数情况都推荐使用该方法。
间接法进行热应力计算的主要步骤如下。
36.1.1热分析瞬态热分析的过程在前例已经介绍过，下面介绍稳态热分析。
稳态热分析用于研究稳定的热载荷对结构的影响，有时还用于瞬态热分析时计算初始温度场。
稳态热分析主要步骤如下。
1．建模稳态热分析的建模过程与其他分析相似，包括定义单元类型、定义单元实常数、定义材料特性、建立几何模型和划分网格等。
但需注意的是：稳态热分析必须定义材料的导热系数。
2．施加载荷和求解(1)指定分析类型。
Main Menu→Solution→Analysis Type→New Analysis,选择 Static.(2)施加载荷。
可以施加的载荷有恒定的温度、热流率、对流、热流密度、生热率，Main Menu→Solution→Define Loads→Apply→Thermal.(3)设置载荷步选项。
普通选项包括时间（用于定义载荷步和子步）、每一载荷步的子步数，以及阶跃选项等, Main Menu→Solution→Load Step Opts→T ime/Frequenc→Time→Time Step.非线性选项包括：—迭代次数（默认25），Main Menu→Solution→Load Step Opts→Nonlinear→Equilibrium Iter;打开自动时间步长，Main Menu→Solution→Load Step Opts→Time/Frequenc→Time→Time Step等.输出选项包括：控制打印的输出，Main Menu→Solution→Load Step Opts→Output Ctrls→Solu Printout;控制结果文件的输出，Main Menu→Solution→Load Step Opts→Output Ctrls→DB/Results File o(4)设置分析选项。
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 温度荷载实例分析命令流
!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。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
/FILNAME,Double,1 ！定义工作文件名。
/TITLE,Temperature Analysis ！定义工作标题。
!*
/PREP7
!定义单元。
ET,1,SOLID70
!*
!定义材料属性。
MPTEMP,,,,,,,, !定义材料1。
MPTEMP,1,0
MPDATA,KXX,1,,238*3.6 !定义材料1的传热系数KXX1。
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,1,,500 !定义材料1的比热C1。
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,3e-6 !定义材料1密度DENS1。
!*
MPTEMP,,,,,,,, !定义材料2。
MPTEMP,1,0
MPDATA,KXX,2,,15*3.6 !定义材料2的传热系数KXX2。
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,2,,100 !定义材料2的比热C2。
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,2,,2.2e-6 !定义材料2密度DENS2。
!*
!建立几何模型。
BLC4,-80,-10,160,20,700
VOFFST,3,20, ,
!*
!网格划分。
FLST,5,20,4,ORDE,2
FITEM,5,1
FITEM,5,-20
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1,10, , , , , , ,1 !定义网格大小。
!*
TYPE, 1
MAT, 1
REAL,
ESYS, 0
SECNUM,
CM,_Y,VOLU
VSEL, , , , 1
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VSWEEP,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
TYPE, 1
MAT, 2
REAL,
ESYS, 0
SECNUM,
CM,_Y,VOLU
VSEL, , , , 2
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VSWEEP,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!定义网格大小完成。
!*
FINISH
/SOL
ANTYPE,4 !定义瞬态分析类型。
TRNOPT,FULL
LUMPM,0
OUTRES,ERASE
OUTRES,ALL,ALL !定义输出每一步结果。
TIME,1 !定义时间及步长。
AUTOTS,-1
DELTIM,0.05, , ,1
KBC,1
TSRES,ERASE
!*
FLST,2,6035,1,ORDE,2
FITEM,2,1
FITEM,2,-6035
IC,P51X,TEMP,20, !定义初始温度。
FLST,2,2,5,ORDE,2
FITEM,2,1
FITEM,2,9
/GO
SFA,P51X,1,HFLUX,100 !定义热流。
!求解。
FINISH
/SOL
/STATUS,SOLU
SOLVE
!重新进入前处理。
/PREP7
ETCHG,TTS
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2.06e5
MPDATA,PRXY,1,,0.3
MPTEMP,,,,,,,,
MPTEMP,1,0
UIMP,1,REFT,,,20
MPDATA,ALPX,1,,0.5e-4 !定义材料1属性。
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,1.5e5
MPDATA,PRXY,2,,0.3
MPTEMP,,,,,,,,
MPTEMP,1,0
UIMP,2,REFT,,,20
MPDATA,ALPX,2,,5e-4 !定义材料2属性。
!重新进入求解模块。
FINISH
/SOL
!*
ANTYPE,0
NSUBST,1,0,0
FLST,2,2,5,ORDE,2
FITEM,2,2
FITEM,2,11
/GO
DA,P51X,ALL,
LDREAD,TEMP,,,0.3, ,'Double','rth',' ' ! 读入热分析的节点温度。
TREF,20, ! 设置参考温度。
!*
/STATUS,SOLU ! 求解计算。
SOLVE。