ANSYS温度场例题分析

合集下载

温度ansys分析

温度ansys分析4 汽车内饰压制成型模具温度场模拟与分析温度在汽车内饰压制成型过程中是一个极其重要的参数，无论是模压料的充模流动阶段还是固化阶段，都是在一定的温度下进行的；如果在充模阶段温度控制的不当将直接影响制品的表面质量和力学性能，具体的说，若模具温度过低则会导致模压料流动性降低，难以充满模腔，若模具温度过高则会引起模压料在模具内未完全成型前就开始固化，并且有可能使制品表面的树脂发生分解，同样得不到合格的汽车内饰制品；另一个重要方面就是在压制的过程中要尽量保持模具温度的均匀分布，如果温度分布不均匀就会导致模压料局部提前固化，还会使制品固化度不均匀甚至发生局部树脂分解，同时也会使得制品脱模后产生较大的翘曲变形。

因此有必要对模具的加热过程及其温度场进行模拟，根据分析结果对模具的加热设计进行优化。

在世界计算机辅助工程领域中有许多CAE软件都具有热分析的功能，我们以目前使用最为广泛的大型通用有限元软件ANSYS来分析汽车顶篷内饰压制成型模具的加热过程及其温度场分布。

4．1 ANSYS有限元分析软件4．1．1 ANSYS简介ANSYS是一种应用广泛的大型通用有限元软件，具有完备的预处理器和后处理器(又称前处理模块和后处理模块)。

目前已广泛应用于核工业、石油化工、航空航天、机械制造、能源、汽车交通、国防、军工、电子、土木工程、造船、轻工、日用家电等工业及科学研究中.ANSYS软件含有多种分析能力。

包括简单线性静态分析和复杂非线性动态分析。

可用来求解结构、流体、电力、电磁场及碰撞等问题。

它包含了前处理器、求解器及后处理器和优化等模块，将有限元分析、计算机图形学和优化技术相组合，已成为解决现代工程学问题必不可少的有力工具。

4．1．2 ANSYS热分析模块ANSYS在处理热分析问题方面具有强大的功能，其不但具有快速的网格划分能力和强大的结果后处理功能，而且还具有非常友好的人机交互界面。

在ANSYS 软件中有五个模块可以进行热分析，如图4．1所示，包括：ANSYS／Multiphysics、ANSYS／Mechanical、ANSYS／Thermal、ANSYS／Flotran和ANSYS／ED。

ANSYS 温度场分析

加载可能有问题，注意对照书中的例子，对热分析不了解。但我看你是把材料1都选中了，然后加了温度40度，这样不对，应该把那个面选中，然后选中面上所有的节点，耦合这些节点的自由度，然后加温度。你这个D和DA的加载命令最好跟书中保持一致。追问那个40度是材料本身的原始温度。如果只加表面的话，会不会只是指的它的表面温度。回答温度可以加在单元上吗？我理解是把所有材料1的节点选中，耦合温度自由度，然后加在其中一个节点上。追问恩，昨天我也发现那个地方错了，已经改过，但基本靠自己摸索了。先不要着急把自己的结果做出来，先看看书，理解一下各种单元的含义，理评价谢谢你这么用心。
刚学习用ANSYS进行温度场分析，情况是这要的，一块铝板放在干冰上，分析铝板（初始温度定为40度）的温度杨分布。分析出来的结果与实验相差太大。主要体现在分析出来的冷冻速度太快了。1S的时候温度就已经达到了零下40多度了。分析了各种原因，参数选的都对，不知道是不是不能直接将干冰的-78度直接加载到铝板的一个面上。命令流如下，不知道有没有错误。FINISH/CLEAR/FILNAME,3D imprinting/TITLE, The Transient Thermal Analysis Of MultiScaffold/UNITS,SI/PREP7ET,1,SOLID70LB1=214 !定义铝的导热系数TH_AL=0.003 !铝板的厚度AP=100 !定义对流换热系数MP,C,1,900MP,KXX,1,LB1MP,DENS,1,2780!**************建模*****************!创建铝板CYLIND,0,0.005,-TH_AL,0,-60,60CYLIND,0.005,0.01,-TH_AL,0,-60,60VGLUE,ALLNUMCMP,ALL!材料属性VSEL,S,,,1VSEL,A,,,2VATT,1,1,1ALLSEL,ALL!划分网格ESIZE,0.0004MSHKEY,1VSWEEP,ALL!*****进行稳态分析********/SOLUANTYPE,TRANSTIMINT,OFFESEL,S,MAT,,1D,ALL,TEMP,40TIME,0.1KBC,0ALLSEL,ALLSOLVE!*****进行瞬态分析********/SOLUANTYPE,TRANSTRNOPT,FULLTIME,50KBC,1DELTIME,1,1,2AUTOTS,ONTIMINT,ONOUTRES,ALL,ALL!**施加载荷**DDELE,ALL,TEMPDA,1,TEMP,-78DA,7,TEMP,-78ASEL,S,AREA,,2ASEL,A,AREA,,8SFA,ALL,,CONV,AP,25ALLSEL,ALLSOLVE

workbench热分析案例

热通量矢量图：通过观察热通量矢量图可以发现热量的传递方向及密度分布情况。
•6
划分网格
网格剖分：采用ansys的mesh块对导入的几何体进行网格划分，网格为四面体网格，网格最大边长为5mm。
•1
定义边界条件
墙壁外表面：采用convection边界条件，设定外界空气温度10℃，换热系数为0.36W/㎡·k。
•2
定义边界条件
墙壁内表面：
裸露于空气的表面采用 convection边界条件，拟定外界空气温度20℃，换热系数为0.36W/㎡·k，与热源接触表面采用耦合边界条件。
•3
定义边界条件
热源：与墙体平行的壁面采用 temperature边界条件，定义其温度为50℃，其余壁面均为绝热边界条件。
•4
结果及分析
温度场云图：
通过显示计算得出的温度场可以看出该模型的最小温度值出现在墙体外表面顶部与底部，在该模型中温度场关于yz平面对称。
•
结果及分析

空间结构分析温度场ansys_命令流

!此APDL程序为本书第六章《太阳辐射作用下钢结构温度场分析》中例题6.2的相关程序!（1）建立工作文件名和工作标题/CLEAR !清除目前ANSYS中所有数据/COM,Thermal !过滤ANSYS的图形截面/FILENAME, TEMP_FIELD !定义工作文件名/TITLE, TEMPERATURE FIELD ANAL YSIS OF A STEEL MEMBER !定义工作标题*AFUN,DEG !设置程序中角单位为度!（2）定义单元类型/PREP7 !进入前处理模块ET,1,SOLID70 !定义热分析单元ET,2,MESH200 !定义网格划分单元KEYOPT,2,1,6 !定义单元2的第一个关键选项KEYOPT,2,2,0 !定义单元2的第二个关键选项!（3）定义材料性能参数MP,KXX,1,45 !定义钢材KXX方向的导热系数MP,KYY,1,45 !定义钢材KYY方向的导热系数MP,KZZ,1,45 !定义钢材KZZ方向的导热系数MP,DENS,1,7850 !定义钢材的密度MP,C,1,465 !定义钢材的比热容!（4）创建有限元模型LOCAL,11,0,0,0,0,0,0,30 !定义局部坐标系11CSYS,11 !激活局部坐标系11K,1,0,0,0 !定义矩形钢管的中心轴线中的关键点1K,2,0,1,0 !定义矩形钢管的中心轴线中的关键点2L,1, 2 !连接关键点1和2，生成矩形钢管的中心轴线k,3,0.5,0,0.5 !定义矩形钢管截面的关键节点3k,4,-0.5,0,0.5 !定义矩形钢管截面的关键节点4k,5,-0.5,0,-0.5 !定义矩形钢管截面的关键节点5k,6,0.5,0,-0.5 !定义矩形钢管截面的关键节点6k,7,0.48,0,0.48 !定义矩形钢管截面的关键节点7k,8,-0.48,0,0.48 !定义矩形钢管截面的关键节点8k,9,-0.48,0,-0.48 !定义矩形钢管截面的关键节点9k,10,0.48,0,-0.48 !定义矩形钢管截面的关键节点10k,11,0.5,0,0.48 !定义矩形钢管截面的关键节点11k,12,-0.5,0,0.48 !定义矩形钢管截面的关键节点12k,13,-0.5,0,-0.48 !定义矩形钢管截面的关键节点13k,14,0.5,0,-0.48 !定义矩形钢管截面的关键节点14!定义矩形钢管截面，为便于网格划分，分四部分创建A,3,4,12,11 !利用关键点3、4、12、11，生成面1A,5,6,14,13 !利用关键点5、6、14、13，生成面2A,11,7,10,14 !利用关键点11、7、10、14，生成面3A,8,12,13,9 !利用关键点8、12、13、9，生成面4!对矩形钢管截面进行网格划分AESIZE,all,0.02 !将所有面网格划分时的单元尺寸设置为0.02m MSHAPE,0,2D !利用四边形单元进行网格划分MSHKEY,1 !采用映射网格划分单元AMESH,ALL !对所有网格进行划分单元TYPE,1 ! 定义单元属性EXTOPT,ESIZE,10,1, !在体生成的方向上单元分割为10份EXTOPT,ACLEAR,1 !在体生成时清除面单元网格EXTOPT,ATTR,0,0,0 !使用有MA T命令定义的材料属性赋给单元MA T,1 !定义单元的材料属性VDRAG,1,2 ,3 , 4, , ,1 !拉伸建立体并形成SOLID70实体单元CSYS,0 ! 激活笛卡尔坐标系!（5）计算日照时间内，即14个小时内的空气的温度值*DIM,t_outdoor,,14 !定义室外空气温度数组为t_outdoor*DO,t,6,19 !给室外空气温度数组t_outdoor赋值*SET,t_outdoor(t-5),35+5*sin(15*(t-5)-45)*ENDDO!（6）计算日照时间内，即14个小时内地面的温度值*DIM,t_earth,,14 !定义地面的温度值数组为t_earth*DO,t,6,19 !给日照时间内地面温度值数组赋值*SET,t_earth(t-5),36+8*sin(15*(t-5)-45)*ENDDO!（7）定义14个时刻太阳高度角正弦数组altitudea1=23.44 !定义太阳赤纬角23.44a2=36.26 !定义地理纬度为北纬36.26*DIM,altitude,,14 !定义日照时间内太阳高度角的正弦值数组为altitude *DO,i,1,14,1aa=cos(a2)*cos(15*(i-1)-90)*cos(a1)aa1=sin(a2)*sin(a1) !根据书中公式（）计算太阳高度角正弦值*SET,altitude(i),aa+aa1*ENDDO!（8）定义14个时刻的太阳方位角余弦数组azimuth并赋值，*DIM,azimuth,,14*DO,i,1,14,1aa=sin(a1)*cos(a2)-cos(a1)*sin(a2)*cos(15*(i-1)-90)aa2=sqrt(1-altitude(i)**2)*SET,azimuth(i),aa/aa2*ENDDO!（9）定义14个时刻的太阳方位角角度数组angle并赋值*DIM,angle,,14*DO,i,1,7,1*SET,angle(i),acos(azimuth(i))*ENDDO*DO,i,8,14,1*SET,angle(i),2*180-acos(azimuth(i))*ENDDO!（11）太阳辐射强度计算的基本参数设置*SET,factorc,0.138 !太阳散射辐射系数*SET,pg,0.35 !地面或者水平面的太阳辐射反射率*SET,factor_a,1370*(1+0.034) !地外太阳直射辐射值*SET,factor_b,0.42 !定义大气消光系数!（12）定义14个时刻的地面太阳辐射直射强度数组GND并赋值*DIM,GND,,14*DO,m,1,14,1*SET,GND(m),factor_a/exp(factor_b/altitude(m))*ENDDO!（13）定义14个时刻结构表面太阳辐射荷载的几何参数*DIM,cos_surf,,4 !表面方位角余弦值*DIM,cos_angle,,4 !表面方位角角度值*DIM,cos_title,,4 !倾角余弦值*DIM,title_angle,,4 !倾角角度值!利用面5上的三个关键点，计算面5的表面方位角余弦、表面方位角弧度、倾角余弦和倾角弧度值kp_1=4 !将面5中的关键点4的编号赋给参数kp_1kp_2=3 !将面5中的关键点3的编号赋给参数kp_2kp_3=16 !将面5中的关键点16的编号赋给参数kp_3k1x=kx(kp_1) !提取关键点4的X坐标值k1y=ky(kp_1) !提取关键点4的Y坐标值k1z=kz(kp_1) !提取关键点4的Z坐标值k2x=kx(kp_2) !提取关键点3的X坐标值k2y=ky(kp_2) !提取关键点3的Y坐标值k2z=kz(kp_2) !提取关键点3的Z坐标值k3x=kx(kp_3) !提取关键点16的X坐标值k3y=ky(kp_3) !提取关键点16的Y坐标值k3z=kz(kp_3) !提取关键点16的Z坐标值!参考空间解析几何相关知识，计算面5的a1=(k2y-k1y)*(k3z-k1z)-(k3y-k1y)*(k2z-k1z)a2=(k3x-k1x)*(k2z-k1z)-(k2x-k1x)*(k3z-k1z)a3=0a4=sqrt(a1**2+a2**2+a3**2)*set,cos_surf(1),abs(a2)/a4cosx=normkx(kp_3,kp_2,kp_1)cosy=normky(kp_3,kp_2,kp_1)*if,cosx,lt,0,and,cosy,lt,0,then*set,cos_angle(1),acos(cos_surf(1))*elseif,cosx,gt,0,and,cosy,gt,0,then*set,cos_angle(1),acos(cos_surf(1))+180*elseif,cosx,lt,0,and,cosy,gt,0,then*set,cos_angle(1),180-acos(cos_surf(1))*else*set,cos_angle(1),360-acos(cos_surf(1))*endifcosz=normkz(kp_3,kp_2,kp_1)*(-1)*set,cos_title(1),cosz*set,title_angle(1),acos(cosz)!利用面21上的三个关键点，计算面21的表面方位角余弦、表面方位角弧度、倾角余弦和倾角弧度值kp_1=13 !将面21中的关键点13的编号赋给参数kp_1kp_2=12 !将面21中的关键点12的编号赋给参数kp_2kp_3=22 !将面21中的关键点22的编号赋给参数kp_3k1x=kx(kp_1) !提取关键点13的X坐标值k1y=ky(kp_1) !提取关键点13的Y坐标值k1z=kz(kp_1) !提取关键点13的Z坐标值k2x=kx(kp_2) !提取关键点12的X坐标值k2y=ky(kp_2) !提取关键点12的Y坐标值k2z=kz(kp_2) !提取关键点12的Z坐标值k3x=kx(kp_3) !提取关键点22的X坐标值k3y=ky(kp_3) !提取关键点22的Y坐标值k3z=kz(kp_3) !提取关键点22的Z坐标值a1=(k2y-k1y)*(k3z-k1z)-(k3y-k1y)*(k2z-k1z)a2=(k3x-k1x)*(k2z-k1z)-(k2x-k1x)*(k3z-k1z)a3=0a4=sqrt(a1**2+a2**2+a3**2)*set,cos_surf(1),abs(a2)/a4cosx=normkx(kp_3,kp_2,kp_1)cosy=normky(kp_3,kp_2,kp_1)*if,cosx,lt,0,and,cosy,lt,0,then*set,cos_angle(2),acos(cos_surf(2))*elseif,cosx,gt,0,and,cosy,gt,0,then*set,cos_angle(2),acos(cos_surf(2))+180*elseif,cosx,lt,0,and,cosy,gt,0,then*set,cos_angle(2),180-acos(cos_surf(2))*else*set,cos_angle(2),360-acos(cos_surf(2))*endifcosz=normkz(kp_3,kp_2,kp_1)*(-1)*set,cos_title(2),cosz*set,title_angle(2),acos(cosz)!利用面18上的三个关键点，计算面18的表面方位角余弦、表面方位角弧度、倾角余弦和倾角弧度值kp_1=11 !将面18中的关键点11的编号赋给参数kp_1kp_2=14 !将面18中的关键点14的编号赋给参数kp_2kp_3=21 !将面18中的关键点21的编号赋给参数kp_3k1x=kx(kp_1) !提取关键点11的X坐标值k1y=ky(kp_1) !提取关键点11的Y坐标值k1z=kz(kp_1) !提取关键点11的Z坐标值k2x=kx(kp_2) !提取关键点14的X坐标值k2y=ky(kp_2) !提取关键点14的Y坐标值k2z=kz(kp_2) !提取关键点14的Z坐标值k3x=kx(kp_3) !提取关键点21的X坐标值k3y=ky(kp_3) !提取关键点21的Y坐标值k3z=kz(kp_3) !提取关键点21的Z坐标值a1=(k2y-k1y)*(k3z-k1z)-(k3y-k1y)*(k2z-k1z)a2=(k3x-k1x)*(k2z-k1z)-(k2x-k1x)*(k3z-k1z)a3=0a4=sqrt(a1**2+a2**2+a3**2)*set,cos_surf(1),abs(a2)/a4cosx=normkx(kp_3,kp_2,kp_1)cosy=normky(kp_3,kp_2,kp_1)*if,cosx,lt,0,and,cosy,lt,0,then*set,cos_angle(3),acos(cos_surf(3))*elseif,cosx,gt,0,and,cosy,gt,0,then*set,cos_angle(3),acos(cos_surf(3))+180*elseif,cosx,lt,0,and,cosy,gt,0,then*set,cos_angle(3),180-acos(cos_surf(3))*else*set,cos_angle(3),360-acos(cos_surf(3))*endifcosz=normkz(kp_3,kp_2,kp_1)*(-1)*set,cos_title(3),cosz*set,title_angle(3),acos(cosz)!利用面10上的三个关键点，计算面10的表面方位角余弦、表面方位角弧度、倾角余弦和倾角弧度值kp_1=6 !将面10中的关键点6的编号赋给参数kp_1kp_2=5 !将面10中的关键点5的编号赋给参数kp_2kp_3=19 !将面10中的关键点19的编号赋给参数kp_3k1x=kx(kp_1) !提取关键点6的X坐标值k1y=ky(kp_1) !提取关键点6的Y坐标值k1z=kz(kp_1) !提取关键点6的Z坐标值k2x=kx(kp_2) !提取关键点5的X坐标值k2y=ky(kp_2) !提取关键点5的Y坐标值k2z=kz(kp_2) !提取关键点5的Z坐标值k3x=kx(kp_3) !提取关键点19的X坐标值k3y=ky(kp_3) !提取关键点19的Y坐标值k3z=kz(kp_3) !提取关键点19的Z坐标值a1=(k2y-k1y)*(k3z-k1z)-(k3y-k1y)*(k2z-k1z)a2=(k3x-k1x)*(k2z-k1z)-(k2x-k1x)*(k3z-k1z)a3=0a4=sqrt(a1**2+a2**2+a3**2)*set,cos_surf(1),abs(a2)/a4cosx=normkx(kp_3,kp_2,kp_1)cosy=normky(kp_3,kp_2,kp_1)*if,cosx,lt,0,and,cosy,lt,0,then*set,cos_angle(4),acos(cos_surf(4))*elseif,cosx,gt,0,and,cosy,gt,0,then*set,cos_angle(4),acos(cos_surf(4))+180*elseif,cosx,lt,0,and,cosy,gt,0,then*set,cos_angle(4),180-acos(cos_surf(4))*else*set,cos_angle(4),360-acos(cos_surf(4))*endifcosz=normkz(kp_3,kp_2,kp_1)*(-1)*set,cos_title(4),cosz*set,title_angle(4),acos(cosz)!（14）定义日照时间内14个时刻的构件表面太阳方位角数组并计算赋值*dim,angle_surf,,4,14 !表面太阳方位角计算*do,m,1,14*do,j,1,4,1*set,angle_surf(j,m),abs(cos_angle(j)-angle(m))*enddo*enddo!（15）定义日照时间内14个时刻的构件表面太阳入射角数组并计算赋值*dim,cos_ps,,4,14!计算入射角*do,m,1,14cosp=sqrt(1-altitude(m)**2)*do,j,1,4,1dd1=altitude(m)*cos(title_angle(j))dd2=cosp*cos(angle_surf(j,m))*sin(title_angle(j))*set,cos_ps(j,m),dd1+dd2*enddo*enddo!（16）定义日照时间内14个时刻的构件表面与水平面之间夹角的余弦值数组并计算赋值*dim,cos_ph,,4!表面与水平面之间的夹角的余弦!利用面5上的三个关键点，计算面5水平面之间的夹角的余弦值并赋值给cos_ph(1)kp_1=4 !将面21中的三个关键点13的编号分别赋给参数kp_1kp_2=3 !将面21中的三个关键点12的编号分别赋给参数kp_2kp_3=16 !将面21中的三个关键点22的编号分别赋给参数kp_3k1x=kx(kp_1) !提取关键点4的X坐标值k1y=ky(kp_1) !提取关键点4的Y坐标值k1z=kz(kp_1) !提取关键点4的Z坐标值k2x=kx(kp_2) !提取关键点3的X坐标值k2y=ky(kp_2) !提取关键点3的Y坐标值k2z=kz(kp_2) !提取关键点3的Z坐标值k3x=kx(kp_3) !提取关键点16的X坐标值k3y=ky(kp_3) !提取关键点16的Y坐标值k3z=kz(kp_3) !提取关键点16的Z坐标值a1=(k2y-k1y)*(k3z-k1z)-(k3y-k1y)*(k2z-k1z)a2=(k3x-k1x)*(k2z-k1z)-(k2x-k1x)*(k3z-k1z)a3=(k2x-k1x)*(k3y-k1y)-(k3x-k1x)*(k2y-k1y)a4=sqrt(a1**2+a2**2+a3**2)*set,cos_ph(1),abs(a3)/a4!利用面21上的三个关键点，计算面21水平面之间的夹角的余弦值并赋值给cos_ph(2) kp_1=13 !将面21中的三个关键点13的编号分别赋给参数kp_1kp_2=12 !将面21中的三个关键点12的编号分别赋给参数kp_2kp_3=22 !将面21中的三个关键点22的编号分别赋给参数kp_3k1x=kx(kp_1) !提取关键点13的X坐标值k1y=ky(kp_1) !提取关键点13的Y坐标值k1z=kz(kp_1) !提取关键点13的Z坐标值k2x=kx(kp_2) !提取关键点12的X坐标值k2y=ky(kp_2) !提取关键点12的Y坐标值k2z=kz(kp_2) !提取关键点12的Z坐标值k3x=kx(kp_3) !提取关键点22的X坐标值k3y=ky(kp_3) !提取关键点22的Y坐标值k3z=kz(kp_3) !提取关键点22的Z坐标值a1=(k2y-k1y)*(k3z-k1z)-(k3y-k1y)*(k2z-k1z)a2=(k3x-k1x)*(k2z-k1z)-(k2x-k1x)*(k3z-k1z)a3=(k2x-k1x)*(k3y-k1y)-(k3x-k1x)*(k2y-k1y)a4=sqrt(a1**2+a2**2+a3**2)*set,cos_ph(2),abs(a3)/a4!利用面18上的三个关键点，计算面18水平面之间的夹角的余弦值并赋值给cos_ph(3) kp_1=11 !将面21中的三个关键点13的编号分别赋给参数kp_1kp_2=14 !将面21中的三个关键点14的编号分别赋给参数kp_2kp_3=21 !将面21中的三个关键点21的编号分别赋给参数kp_3k1x=kx(kp_1) !提取关键点11的X坐标值k1y=ky(kp_1) !提取关键点11的Y坐标值k1z=kz(kp_1) !提取关键点11的Z坐标值k2x=kx(kp_2) !提取关键点14的X坐标值k2y=ky(kp_2) !提取关键点14的Y坐标值k2z=kz(kp_2) !提取关键点14的Z坐标值k3x=kx(kp_3) !提取关键点21的X坐标值k3y=ky(kp_3) !提取关键点21的Y坐标值k3z=kz(kp_3) !提取关键点21的Z坐标值a1=(k2y-k1y)*(k3z-k1z)-(k3y-k1y)*(k2z-k1z)a2=(k3x-k1x)*(k2z-k1z)-(k2x-k1x)*(k3z-k1z)a3=(k2x-k1x)*(k3y-k1y)-(k3x-k1x)*(k2y-k1y)a4=sqrt(a1**2+a2**2+a3**2)*set,cos_ph(3),abs(a3)/a4!利用面10上的三个关键点，计算面10水平面之间的夹角的余弦值并赋值给cos_ph(4) kp_1=6 !将面21中的三个关键点6的编号分别赋给参数kp_1kp_2=5 !将面21中的三个关键点5的编号分别赋给参数kp_2kp_3=19 !将面21中的三个关键点19的编号分别赋给参数kp_3k1x=kx(kp_1) !提取关键点6的X坐标值k1y=ky(kp_1) !提取关键点6的Y坐标值k1z=kz(kp_1) !提取关键点6的Z坐标值k2x=kx(kp_2) !提取关键点5的X坐标值k2y=ky(kp_2) !提取关键点5的Y坐标值k2z=kz(kp_2) !提取关键点5的Z坐标值k3x=kx(kp_3) !提取关键点19的X坐标值k3y=ky(kp_3) !提取关键点19的Y坐标值k3z=kz(kp_3) !提取关键点19的Z坐标值a1=(k2y-k1y)*(k3z-k1z)-(k3y-k1y)*(k2z-k1z)a2=(k3x-k1x)*(k2z-k1z)-(k2x-k1x)*(k3z-k1z)a3=(k2x-k1x)*(k3y-k1y)-(k3x-k1x)*(k2y-k1y)a4=sqrt(a1**2+a2**2+a3**2)*set,cos_ph(4),abs(a3)/a4!(17)定义日照时间内的14个时刻的4个面的太阳直射强度数组并计算赋值*dim,vertical,,4,14!定义表面太阳辐射直射强度数组*do,m,1,14,1*do,j,1,4,1*if,cos_ps(j,m),gt,0,then*set,vertical(j,m),gnd(m)*cos_ps(j,m)*else*endif*enddo*enddo!(18)定义日照时间内的14个时刻的4个面的太阳辐射散射强度数组并计算赋值*dim,diffuse,,4,14!定义表面太阳辐射散射强度数组*DIM,fuhao,,4 !定义*do,m,1,14,1!利用面5上的三个关键点，计算当前时刻的面5的太阳辐射散射强度并赋值给diffuse(1,m) kp_1=4 !将面5中的三个关键点4的编号分别赋给参数kp_1kp_2=3 !将面5中的三个关键点3的编号分别赋给参数kp_2kp_3=16 !将面5中的三个关键点16的编号分别赋给参数kp_3cos=normkz(kp_3,kp_2,kp_1)*if,cos,lt,0,then*set,fuhao(j),100*set,diffuse(1,m),gnd(m)*factorc*(1+cos_ph(1))/2*else*set,diffuse(1,m),0*set,fuhao(1),0*endif*enddo*do,m,1,14,1kp_1=13 !将面21中的三个关键点13的编号分别赋给参数kp_1kp_2=12 !将面21中的三个关键点12的编号分别赋给参数kp_2kp_3=22 !将面21中的三个关键点22的编号分别赋给参数kp_3cos=normkz(kp_3,kp_2,kp_1)*if,cos,lt,0,then*set,fuhao(j),100*set,diffuse(2,m),gnd(m)*factorc*(1+cos_ph(2))/2*else*set,diffuse(2,m),0*set,fuhao(2),0*endif*enddo*do,m,1,14,1kp_1=11 !将面18中的三个关键点11的编号分别赋给参数kp_1kp_2=14 !将面18中的三个关键点14的编号分别赋给参数kp_2kp_3=21 !将面18中的三个关键点21的编号分别赋给参数kp_3cos=normkz(kp_3,kp_2,kp_1)*if,cos,lt,0,then*set,diffuse(3,m),gnd(m)*factorc*(1+cos_ph(3))/2*else*set,diffuse(3,m),0*set,fuhao(3),0*endif*enddo*do,m,1,14,1kp_1=6 !将面10中的三个关键点6的编号分别赋给参数kp_1kp_2=5 !将面10中的三个关键点5的编号分别赋给参数kp_2kp_3=19 !将面10中的三个关键点19的编号分别赋给参数kp_3cos=normkz(kp_3,kp_2,kp_1)*if,cos,lt,0,then*set,fuhao(j),100*set,diffuse(4,m),gnd(m)*factorc*(1+cos_ph(4))/2*else*set,diffuse(4,m),0*set,fuhao(4),0*endif*enddo!(19)定义日照时间内的14个时刻的4个面的太阳辐射反射强度数组并计算赋值*dim,reflect,,4,14!表面太阳辐射反射强度*do,m,1,14,1kp_1=4 !将面5中的三个关键点4的编号分别赋给参数kp_1kp_2=3 !将面5中的三个关键点2的编号分别赋给参数kp_2kp_3=16 !将面5中的三个关键点16的编号分别赋给参数kp_3cos=normkz(kp_3,kp_2,kp_1)*if,cos,gt,0,then*set,reflect(1,m),gnd(m)*(cos(22.76*3.14/180)+factorc)*pg*(1-cos_ph(1))/2*else*set,reflect(1,m),0*endif*enddo*do,m,1,14,1kp_1=13 !将面21中的三个关键点13的编号分别赋给参数kp_1kp_2=12 !将面21中的三个关键点12的编号分别赋给参数kp_2kp_3=22 !将面21中的三个关键点22的编号分别赋给参数kp_3cos=normkz(kp_3,kp_2,kp_1)*if,cos,gt,0,then*set,reflect(2,m),gnd(m)*(cos(22.76*3.14/180)+factorc)*pg*(1-cos_ph(2))/2*else*set,reflect(2,m),0*endif*enddo*do,m,1,14,1kp_1=11 !将面18中的三个关键点11的编号分别赋给参数kp_1kp_2=14 !将面18中的三个关键点14的编号分别赋给参数kp_2kp_3=21 !将面18中的三个关键点21的编号分别赋给参数kp_3cos=normkz(kp_3,kp_2,kp_1)*if,cos,gt,0,then*set,reflect(3,m),gnd(m)*(cos(22.76*3.14/180)+factorc)*pg*(1-cos_ph(3))/2*else*set,reflect(3,m),0*endif*enddo*do,m,1,14,1kp_1=6 !将面10中的三个关键点6的编号分别赋给参数kp_1kp_2=5 !将面10中的三个关键点5的编号分别赋给参数kp_2kp_3=19 !将面10中的三个关键点19的编号分别赋给参数kp_3cos=normkz(kp_3,kp_2,kp_1)*if,cos,gt,0,then*set,reflect(4,m),gnd(m)*(cos(22.76*3.14/180)+factorc)*pg*(1-cos_ph(4))/2*else*set,reflect(4,m),0*endif*enddo!(20)定义日照时间内的14个时刻的4个面的太阳辐射总强度数组并计算赋值*dim,radi_all,,4,14!表面辐射总强度*do,m,1,14,1*do,j,1,4,1*set,radi_all(j,m),(vertical(j,m)+diffuse(j,m)+reflect(j,m))*0.55*enddo*enddo!(21)定义日照时间内的14个时刻的4个面的太阳辐射等效生热率数组并计算赋值!等效生热率*dim,heat_eq,,4,14*do,m,1,14,1*do,j,1,4,1*set,heat_eq(j,m),radi_all(j,m)/0.02*enddo*enddofinish!(21) 进入求解器，设置求解选项进行数值计算/soluantype,trans !设置求解类型瞬态热分析trnopt,full ! 指定瞬态分析的求解方法为完全法timint,on ! 打开时间积分效应tunif,32.5 !设置构件的初始温度为32.5摄氏度outres,all ! 将除SV AR和LOCI以外的所有计算数据写入数据库和文件中!求解参数设定value_1=0.8 ! 定义构件的辐射发射率参数value_2=5.67e-8 !定义斯蒂芬-玻尔慈曼常数参数value_3=value_1*value_2 !将value_1和value_2相乘并赋值给value_3*dim,long_wave,,4,840 !求解过程中，每个荷载步中给四个外表面所施加的长波辐射强度*dim,node_num1,,2000 !拱上实体包含的节点编号数组*dim,node_temp1,,2000 !拱上实体包含节点的温度值数组l=0*do,m,1,14,1*do,r,3600,3600,3600time,(m-1)*3600+rl=l+1nsubst,1,100,1 !设置每个荷载的子步数数为1，最大值为100，最小值为1autots,on ! 打开自动时间步长跟踪eqslv,JCG !指定方程求解器为JCGkbc,0 !使用递增方式加载!施加对流荷载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施加当前荷载步的等效生热率solvevsel,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在土木工程结构中，温度应力在很多情况下对结构的影响很大。

自紧密封U型管式换热器管板温度场的ANSYS分析

越大。
关键词
有限元分析
管板
温度场
表皮效应Ａ
中图法分类号
ＴＩ２４Ｋ７．；
文献标志码
管壳式换热器有许多类型，每种换热器的管而
的过程中选用的是ＳＥＬ３单元和ＳＬＤ０热分ＨＬ６ＯＩ７
板都存在较大的温差。温差越大，由此而产生的温
对流传热系数／Ｗ・（・）［ｍ ℃
腐蚀量／ｍｍ
］９８９６．
５
ｌ６２０
３
时，自紧式密封就更显示出其优越性。陔自紧密封
本文中的换热器管板是在高温高压下工作的，
且管程、壳程温差较大，管板将会产生较大的温则差应力，因此，对管板进行温度场的分析，而了解从管板的温度场分布情况是非常有ｃｅｈＥｇｇ２．
动力技术
自紧密封Ｕ型管式换热器管板温度场的ＡＹＮＳＳ分析
林玉娟孙士财冯永利李雷明。
（东北石油大学机械科学与工程学院，大庆１３１中国石油克拉玛依石化公司，６３８；克拉玛依８４０；３０３中国石油大庆炼化公司，大庆１３１）６４１
４期
林玉娟，：等自紧密封ｕ型管式换热器管板温度场的ＡＳＳ分析ＮＹ
表１几何模型的主要结构尺寸

ansys 求当结构内侧温度为0度而外侧温度为-10度时的变形

ansys 求当结构内侧温度为0度而外侧温度为-10度时的变形1.引言1.1 概述概述本文旨在研究在给定条件下的结构变形问题。

具体而言，我们关注的是当结构内侧温度为0度，而外侧温度为-10度时，结构的变形情况。

为了解决这一问题，我们采用了ANSYS模拟方法。

ANSYS作为一种常用的工程仿真软件，能够对复杂结构进行力学分析，并得出相应的结果。

通过在ANSYS中设置合适的边界条件和材料属性，我们可以模拟不同温度下的结构变形情况。

在本文中，我们首先进行了背景介绍，介绍了该问题的背景和相关研究现状。

然后，我们详细讨论了使用ANSYS进行模拟的方法。

通过分析和讨论模拟结果，我们将得出对结构变形情况的评估和结论。

通过本文的研究，我们希望能够揭示在不同温度条件下结构的变形情况，为相关工程实践提供参考和指导。

同时，我们也期待这项研究能够为进一步探索结构变形领域提供一定的理论和实践基础。

1.2文章结构文章结构部分的内容可以按照以下方式编写：1.2 文章结构本文主要分为引言、正文和结论三个部分。

引言部分主要包括对本篇文章的概述、文章结构和研究目的的介绍。

首先概述研究的背景和意义，说明研究的目的和重要性。

然后介绍文章的结构，即各个章节的内容和逻辑顺序。

最后明确研究的目的，指出本文的研究重点和亮点。

正文部分主要包括背景介绍和ANSYS模拟方法。

在背景介绍章节中，可以从结构分析的角度出发，说明为什么需要对结构内侧温度为0度而外侧温度为-10度时的变形进行研究。

介绍相关的理论知识和前人研究成果，展示研究的基础和现有问题。

在ANSYS模拟方法章节中，详细介绍使用ANSYS软件进行结构分析的步骤和方法。

包括建模、网格划分、边界条件设置等。

同时，可以说明为什么选择ANSYS作为研究工具，以及其在结构分析领域的优势和应用。

结论部分主要包括对结果的分析和结构变形评估。

在结果分析章节中，对模拟结果进行定量或定性的分析，解释温度变化对结构变形的影响。

ANSYS温度场例题分析

短圆柱体的热传导过程问题：一短圆柱体，直径和高度均为1m，现在其上端面施加大小为100℃的均匀温度载荷，圆柱体下端面及侧面的温度均为0℃，试求圆柱体内部的温度场分布（假设圆柱体不与外界发生热交换）。

圆柱体材料的热传导系数为30W/(m·℃)。

求解：第一步：建立工作文件名和工作标题在ANSYS软件中建立相应的文件夹，并选择Thermal复选框。

第二部：定义单元类型在单元类型（element type）中选择thermal solid和quad 4node 55，在单元类型选择数字（element type reference number）输入框中输入1，在单元类型选择框里选择Axisymmetric，其余默认即可。

第三步：定义材料性能参数在材料性能参数对话框中输入圆柱体的导热系数30.第四步：创建几何模型、划分网格创建数据点，输入点坐标。

在第一个输入框中输入关键点编号1，并输入第一个关键点坐标0、0、0，重复输入第二个、第三个、第四个关键点，相应的坐标分别为2（0.5,0,0）；3（0.5,1,0）；4（0,1,0）。

结果如下图1所示：在模型中创建直线，选择编号为1、2的关键点生成一条直线，在选取2、3生成一条直线，同样选择编号为3、4和编号为4、1的关键点生成另外两条直线。

结果如下图2所示：之后在plot numbering controls对话框，分别打开KP Keypoint numbers、LINE line numbers、AREA Area numbers，建立直线L1、L2、L3、L4线段。

生成几何模型，如下图所示：在L1、L3线段上划分20个单元，并将L2、L4划分成40个单元格，并在模型上选取编号为A1的平面，如下图所示：将结果进行保存。

第五步：加载求解选择分析类型Steady-State，在Select Entities对话框，第一个下拉列表框中选择Lines，在第二个下拉列表中选择By Num，第三个单选框中选择From Full。

例4受热载荷作用的厚壁圆筒的有限元建模与温度场求解(ANSYS有限元分析入门)

受热载荷作用的厚壁圆筒的热—结构耦合分析计算分析模型如图4-1 所示, 习题文件名: Thermal。

（轴对称问题）R1=0.3 R2=0.5 底面固定圆筒内壁温度:500℃，外壁温度:100℃。

两端自由且绝热，热导率0.07弹性模量2.06e11泊松比0.3热膨胀系数1e-5高为1m图4-1受热载荷作用的厚壁圆筒的计算分析模型（截面图）4.1进入ANSYS程序→AnsysED 10.0→ANSYS（启动）→file →change jobname：Thermal（更改工作名）→change title（更改标题）thermal→ok→plot （绘制）→replot（重新绘制）4.2选择单元类型ANSYS Main Menu: Preprocessor（前处理器）→Element T ype（单元类型）→Add/Edit/Delete… （添加/编辑/删除）→Add…（添加）→select Thermal Solid Quad 4node 55 →OK (back to Element T ypes window)→Options… →select K3: Axisymmetric（轴对称问题）OK→Close (the Element T ype window)4.3定义材料参数ANSYS Main Menu: Preprocessor →Material Props（材料参数）→Material Models（材料模型）→Thermal（热参数）→Conductivity（热导率）→Isotropi c→input KXX:0.07→OK→Structural（结构参数）→Linea r→Elastic→Isotropic →input EX:2.06e11，PRXY：0.3→OK→Structural（结构参数）→Thermal Expansion（热膨胀）→Secant Coefficient→Isotropic →input ALPX:1e-5→OK（由于是热——结构耦合分析，需要同时定义热参数和结构参数）4.4生成几何模型生成特征点ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints（创建关键点）→In Active CS（在当前坐标系）→依次输入四个点的坐标：input:1(0.3,0),2(0.5,0),3(0.5,1),4(0.3,1)→OK✓生成圆柱体截面ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Arbitrary →Through KPS →依次连接四个特征点，1(0.3,0),2(0.5,0),3(0.5,1),4(0.3,1) →OK4.5网格划分ANSYS Main Menu: Preprocessor →Meshing →Mesh T ool（网格划分工具）→(Size Controls) lines: Set →拾取两条水平边:OK→input NDIV: 5 →Apply→拾取两条竖直边:OK →input NDIV: 15 （以上两步为设置各个线上网格划分的密度）→OK →Mesh T ool→Mesh →Pick All (in Picking Menu) →Close( the Mesh T ool window)4.6模型施加约束✓分别给两条直边施加约束ANSYS Main Menu: Solution（求解）→Define Loads（定义载荷）→Apply（施加）→Thermal（热载荷）→T emperature（边界温度）→On Lines →拾取左边, Value: 500 →Apply →拾取右边，Value：100 →OK（对两边的温度进行设置）4.8 分析计算ANSYS Main Menu: Solution →Solve （求解）→Current LS（当前坐标系）→OK(to close the solve Current Load Step window) →OK4.9 结果显示ANSYS Main Menu: General Postproc（通用后处理器）→Plot Results（绘制结果）→Contour Plot→Nodal Solu（节点解）→Nodal T emperature(节点温度)→OK4.10 开始结构分析ANSYS Main Menu: Preprocessor（前处理器）→Element T ype（单元类型）→Switch Elem T ype（转换单元类型）4.11 施加边界位移条件ANSYS Main Menu: Solution（求解）→Define Loads（定义载荷）→Apply（施加）→Structura l→Displacement→On Lines→拾取底边→OK→ALL DOF→OK4.12 将热分析应力分布读入当前工况ANSYS Main Menu: Solution（求解）→Define Loads（定义载荷）→Apply（施加）→Thermal →T emperature →From Therm Analy（从热分析读入）→Browse→从工作目录中读取thermal.rth 文件（.rth是热分析结果文件的后缀）4.13 分析计算ANSYS Main Menu: Solution →Solve （求解）→Current LS（当前坐标系）→OK(to close the solve Current Load Step window) →OK4.14 结果显示ANSYS Main Menu: General Postproc（通用后处理器）→Plot Results（绘制结果）→Contour Plot→Nodal Solu（节点解）→Stress（应力）→1st principal stress（第一强度理论）4.15 退出系统ANSYS Utility Menu: File→Exit…→Save Everything→OK。

Ansys有限元分析温度场模拟指导书

实验名称：温度场有限元分析一、实验目的1. 掌握Ansys分析温度场方法2. 掌握温度场几何模型二、问题描述井式炉炉壁材料由三层组成，最外一层为膨胀珍珠岩，中间为硅藻土砖构成，最里层为轻质耐火黏土砖，井式炉可简化为圆筒，筒内为高温炉气，筒外为室温空气，求内外壁温度及温度分布。

井式炉炉壁体材料的各项参数见表1。

表1 井式炉炉壁材料的各项参数三、分析过程1. 启动ANSYS，定义标题。

单击Utility Menu→File→Change Title菜单，定义分析标题为“Steady-state thermal analysis of submarine”2.定义单位制。

在命令流窗口中输入“/UNITS, SI”，并按Enter 键3. 定义二维热单元。

单击Main Menu→Preprocessor→Element Type→Add/Edit/Delete 菜单，选择Quad 4node 55定义二维热单元PLANE554.定义材料参数。

单击Main Menu→Preprocessor→Material Props→Material Models菜单5. 在右侧列表框中依次单击Thermal→Conductivity→Isotropic，在KXX文本框中输入膨胀珍珠岩的导热系数0.04，单击OK。

6. 重复步骤4和5分别定义硅藻土砖和轻质耐火黏土砖的导热系数为0.159和0.08，点击Material新建Material Model菜单。

7.建立模型。

单击Main Menu→Preprocessor→Modeling→Create→Areas→Circle→By Dimensions菜单。

在RAD1文本框中输入0.86，在RAD2文本框中输入0.86-0.065，在THERA1文本框中输入-3，在THERA2文本框中输入3，单击APPL Y按钮。

8.重复第7步，输入RAD1=0.86-0.065，RAD2=0.86-0.245，单击APPL Y；输入RAD1=0.86-0.245，RAD2=0.86-0.36，单击OK。

Ansys计算温度场操作流程学习资料

精品文档Instruction of Ansys temperature field calculationQuestion 1: Con sider an infin ite (in one directio n) plate with in itial temperature T 0. One end of the plate is exposed to the environment of which the temperature is T e (III type boundary condition). Analyze the temperature distribution in the plate during the period of 2000s.问题1考虑一个方向无限长的平板，初始温度为T0，一段暴露在温度为T e的环境中，分析其在2000s内温度分布情况。

12 24Basic parameters基本物性参数Geometry 几何：a=1 m, b=0.1 mo 3 oMaterial 材料：2=54 W/m • C,尸7800 kg/m , c p=465 J/kg C Loads载荷：T0=0 °C, T e=1000 o C, h=50 W/m2 o CJobname and directory setting设置文件名、存储路径Menu | File | Change Job nameMenu | File | Change DirectoryPreprocessinc前处理⑴Define Element Type定义单元类型Preprocessor | Eleme nt Type | Add/Edit/DeleteAdd: Thermal Mass | Solid | Quad 4node 55(2) Set Material Properties 设置材料属性Preprocessor | Material Props | Material ModelsThermal: Co nductivity: Isotropic KXX=54Thermal: Den sity=7800Thermal: Specific Heat=465精品文档Modeling 建模(1) Create Node 1 建立节点1Preprocessor | Modeling | Create | Nodes | In Active CS No.: 1, (x, y, z) = (0,0,0)(2) Create Node 12建立节点12 Preprocessor | Modeling | Create | Nodes | In Active CS No.: 2, (x, y, z) = (0,1,0)(3) Fill Between Node 1 and 12在节点1，12 间填充其余节点Preprocessor | Modeling | Create | Nodes | Fill Between Nds Number of nodes to fill: 10Spacing ratio: 1(均匀网格)(4) Create Node 13~24 by copying复制生成节点13〜24Preprocessor | Modeling | Copy | Nodes | Copy Pick All 选择所有节点Total number of copies: 2复制2 份(包含原先的1 份) X-offset: 0.1 设置X 方向偏移量(5) Create Element 1 生成单元1Preprocessor | Modeling | Create | Elements | Auto Numbered | Thru Nodes Select 1, 13, 14, 2(with mouse)(6) Copy Element 2~10 by copying复制生成其余单元Preprocessor | Modeling | Copy | Elements | Auto Numbered Pick All 选择所有单元Total number of copies: 11复制11 份(包含原先的1 份)Node number increment: 1设置节点增量Applying Loads加载(设置边界条件和初始条件)(1) Apply convention loads on Node 12 and 2在4 节点12，24 上加对流载荷Preprocessor | Loads | Define Loads | Apply | Thermal | Convection | On Nodes Pick 12, 24VALI Film coefficient: 50 换热系数VAL2L Bulk temperature: 1000 环境温度(2) Set initial temperature(uniform temperature)设置初始条件(均匀温度)Solution | Define Loads | Settings | Uniform Temp Uniform Temperature: 0Solution 求解(1) Set Analysis Type as Transien t置分析类型为瞬态分析Solution | Analysis Type | New Analysis Transient精品文档(2) Set solution control设置分析控制参数Solution | Analysis Type | So'ln Controls”Transient”Label: Stepped Loading”Basic”Label: Time at end of loadstep: 10000; Automatic time stepping: on; Time increment: 10; Frequency: Write every substep⑶Solve求解Solution | Solve | Current LSViewing Results 查看结果(1) TimeHist Postpro时间历程后处理器Add data 添加数据:TimeHist Postpro | Add data: Nodal Solution: DOF Solution: Nodal Temperature:Variable Name: T1; Node 1Variable Name: T5; Node 5Variable Name: T9; Node 9Variable Name: T11; Node 11(2) List/Graph 列表/作图Select the variables选择变量：List data/ Graph data⑶Save image保存图片Menu | PlotCtrls | Capture Image： File | Save as(4) General Postpro通用后处理器Select the data of certain time选择某一时刻的数据：General Postpro | Read Results | By pickPlot the con tour map of temperature distribution 绘制该时刻温度分布等值线图：General Postpro | Plot Results | Contour Plot | Nodal Solu | Nodal Solution： DOF Solution：Temperature(5) Save image保存图片Menu | PlotCtrls | Capture Image： File | Save as(6) Animate 生成动画Menu | PlotCtrls | Animate | Over TimeNumber of animation frames： 100Question 2： Consider a solid cylinder with initial temperature T0. The top surface of the cylinder is exposed to the environment of which the temperature is T e (III type boundary condition) and the other surfaces are heat insulation (II type boundary condition). Analyze the temperature distribution in the cylinder during the period of 1000s.问题2:考虑一个初始温度为T o的圆柱体，一段暴露在温度为T e的环境中，其余界面视为绝热，分析其在1000s内温度分布情况。

ANSYS温度场分析步骤

ANSYS温度场分析步骤
基于ANSYS12.0的钢板加热过程分析
一．问题描述
2000mm*2000mm*100mm的钢板，初始温度为20℃，放入温度为1120℃的加热炉内加热,已知其换热系数125W/㎡*K，钢板的比热为460J/kg*℃，密度为7850kg/m 3，导热系数为50W/m*K，计算钢板1800s后的温度场分布。

二．问题分析
此问题属于热瞬态分析（载荷随时间变化），选用SOLID70三维六面体单元进行有限元分析。

SOLID70——三维热实体，具有8各节点，每个节点一个温度自由度。

该单元可用于三维的稳态或瞬态的热分析问题。

三．操作步骤
1.定义分析文件名
Utility Menu>File>Change Jobname，输入Example。

2.定义单元类型
Main Menu>Preprocesor>Element Type>Add/Edit/Delete,选择SOLID70三维六面体单元进行有限元分析。

3.定义材料属性
①传导系数
②材料密度
③材料比热
4.建立几何模型
5.设置单元密度
6.划分单元
7.施加对流换热载荷
8.施加初始温度
9.设置求解选项
10.温度偏移量设置
11.输出控制
12.存盘
13.求解
14.显示温度场分布云图
四．总结
本例介绍了应用ANSYS对钢板加热过程进行瞬态热分析的基本步骤，应用此方法可对各种零件加热过程的温度场分布进行分析。

基于AnsysWorkbench雅阁ISG温度场仿真分析

基于AnsysWorkbench雅阁ISG温度场仿真分析本文基于Ansys Workbench对雅阁ISG的温度场进行了仿真分析。

ISG是内燃机启动器和发电机的组合装置，也称为轴承式起动机（Starter Generator，简称SG），是目前汽车发动机的“绿色”起动技术之一。

首先，我们需要构建ISG的三维模型，并设置ISG工作时的工况条件，包括工作电流、转速等。

然后，我们将模型导入Ansys Workbench中，通过选择热传导法，建立ISG的温度场分析。

在分析过程中，我们可以将ISG的温度场分为静态和动态两种情况进行分析。

其中，静态分析主要用于分析ISG在静止状态下的温度分布情况，而动态分析则可以直观地反映ISG在工作状态下的温度场分布情况。

通过静态分析，我们可以发现ISG在不同位置的温度分布存在一定的差异。

其中，发电机部分温度分布状态相对均匀，而起动机部分温度分布则表现出较强的集中性，这主要是由于起动机部分工作时电磁场的分布差异所导致的。

而通过动态分析，我们可以得知ISG在不同工作状态下的温度分布情况也会有所不同。

例如，在高负载状态下，ISG的温度分布相对均匀而稳定，在低负载状态下则出现温度分布的不均匀性。

最后，我们可以对ISG的改进进行模拟分析，以寻找最优的改进方案。

例如，可以通过对ISG内部的散热结构进行优化设计，以提高ISG的散热效率并减少温度的集中分布。

综上所述，通过Ansys Workbench的仿真分析，我们可以深入研究ISG的温度场分布情况，并寻找最优的改进方案，以提高ISG的效率和稳定性。

此外，在ISG使用过程中，温度对于ISG的运行状态有着重要的影响。

温度过高会导致ISG内部元件的热膨胀而失去原本的机械性能，从而导致ISG的故障或损坏，进一步影响到整个发动机的运行状态。

因此，在ISG的设计过程中，需要考虑机械结构和散热系统的优化，以确保其能够承受各种环境下的温度影响而稳定运行。

混凝土箱梁日照温度场温度应力ansys分析结果

SimWe仿真论坛»C06：ANSYS--实例赏评»混凝土箱梁日照温度场、温度应力ANSYS分析结果混凝土箱梁日照温度场、温度应力ANSYS分析结果混凝土箱梁在日照和气温变化等气象因素作用下，会在截面内产生非线性温度分布，引起较大的纵向、横向温度应力，在超静定结构中还会引起温度次应力。

应力大小往往会超过列车或汽车荷载效应，特别是横向温度应力对混凝土箱梁纵向裂纹的出现有很大的贡献。

下面首先发几张混凝土箱梁日照温度场ANSYS分析结果的图片，希望对这方面感兴趣的网友在此讨论。

Ⅰ：夏季日照温度场。

由于，桥轴线走向和纬度的关系，腹板在夏季腹板几乎不受日照，因此截面温度梯度主要在竖向。

peregrine2007-7-14 15:07夏季，t=10:00的温度场peregrine2007-7-14 15:09夏季，t=14:00的温度场[[i] 本帖最后由 peregrine 于 2007-7-14 15:15 编辑 [/i]]peregrine2007-7-14 15:15回复 #3 peregrine 的帖子夏季，t=03:00，夜间负温差peregrine2007-7-14 15:19Ⅱ：冬季温度场。

本箱梁冬季腹板也会受到一定的日照。

冬季，t=16:00bridge-7-18 21:481、底板温度基本是处于均匀温度状态原来做过实桥试验，上下底板也是相差很大的，是不是所处环境不同了2、“夏季，t=03:00，夜间负温差”跟实测也是差的很远，基本上是处于均匀温度状态。

3、希望提供你的计算思路，偶们好学习一下。

peregrine2007-7-19 20:15回复 #6 bridge5209 的帖子回楼上我这是根据多年气象资料计算的最不利状况下的温度分布，与楼上在某一座桥的实测数据有出入，是正常的。

1、底板温差主要受气温变化和地面或水面对太阳辐射的反射率影响，地面太阳辐射发射率随环境变化很大，难以准确确定，计算时一般偏于不利考虑，取较小值，因此计算的底板上下温差比较小，在本算例中为℃（14:00)2、夜间负温差看起来很大，但要注意的是，最高温度出现在箱梁梗胁加厚处的内部，而最低温度出现在悬臂端部板厚最薄处，特别是在悬臂端部，在很小的范围内温度降低很多，因为这个部位不仅尺寸小，而且夜间呈三面放热的状态，温度下降自然比结构主体要大得多。

T型接头焊接温度场ANSYS仿真分析

焊缝凝固过程的温度场分析初始条件：焊接件的初始温度为25度，焊缝温度为3000;对流边界条件：表面传热系数为5e-4，比热容0.2，材料密度0。

28，空气温度为25度;求2000s后整个焊接件的温度分布1、选择网格单元类型Preprocessor〉Element Type>Add/Edit/Delete>Add〉Thermal Mass〉Solid>Brick 8 node 70图1—1 定义单元类型2、设置钢板及焊缝材料属性Preprocessor〉Material Props〉Material Models〉Material Model Number 1>Thermal a．设置焊件材料密度、热传导系数、比热容，设置焊缝材料密度、热传导系数、比热容及与温度相关的涵参数,如下图所示.b．设置左右两道焊缝的焓参数，焓参数随温度变化曲线如图2-5所示。

图2—1 钢板热导率设置图2—2 设置钢板比热容图2-3 设置钢板密度图2—4 焊缝焓参数设置图2—5 左右焊缝焓参数3、建立几何模型Preprocessor〉Modeling>Create>Volumes>Block〉By Dimensions 建立焊件几何模型。

Preprocessor〉Modeling>Create〉Volumes>Cylinder>By Dimensions 建立焊缝几何模型。

建模过程如图3-1所示。

图3—1 几何模型建模过程1图3-2 几何模型建模过程2通过Reflect建立完整的几何模型，之后运用布尔运算中glue使整个模型成为一个整体，如图3-3所示.焊接模型几何参数：横板：2*1.2＊0。

4竖板：0。

4*1.2*1焊缝：R0.2*1。

2图3-3 焊件几何模型设置焊件及左右焊缝网格属性Preprocessor〉Meshing>Mesh Attributes〉Picked 选择焊件或是焊缝，分别对其进行设置。

ansys变压器温度场分析

变压器温度场的有限元分析华北科技学院机制B091班220 kV大容量油浸式电力变压器温度场的有限元分析随着电力建设的不断发展，电力设备朝着大型化方向发展，但大型变压器在内部温升控制方面一直是近年来困扰变压器制造企业的技术难题之一。

如何开发合适的温度场计算技术，准确地计算变压器在各种运行状态下内部线圈、结构件及铁芯等部位的温度，控制内部热点温度不超过其内部绝缘材料的许用温度，从而保证变压器的热寿命，提高变压器的安全可靠性，是企业急需解决的问题。

本课题采用有限元技术对大容量变压器的温度场进行数值模拟与分析。

主要研究内容有：1)了解和掌握大型变压器的结构与工作状况；2)建立典型的220 kV变压器的有限元分析模型；3)分析与计算变压器的温度场分布；表1 该变压器主要参数定义油和所有线圈的相对磁导率为1，油箱的磁导率为300，屏蔽材料的磁导率为2000。

表2 变压器结构参数(1)变压器几何尺寸在宽度方向与x轴对称:铁心窗口上、下边距绕组高度中心线相等油箱上、下盖的内侧距绕组高度中心线相等。

(2)忽略励磁电流和环流，一、二次侧绕组安匝平衡。

(3)油箱磁屏蔽的μ→∞，可作为边界而处理。

(4)旁扼截面可看成矩形截面。

电力变压器是电力网中的主要设备，其总容量达到发电设备总容量的5～6 倍。

电力变压器的技术性能、经济指标直接影响着电力系统的安全性、可靠性和经济性。

随着科学技术的发展、生产技术的进步以及新型电工材料的开发应用，变压器的各项性能指标不断刷新，单机容量越来越大，变压器中的漏磁场也随之增大，引起了人们的关注。

在额定运行情况下，漏磁场的增强引起的变压器附加损耗的增加将直接影响变压器的行效率和产品的竞争力。

严重的是，由于漏磁场在一定范围内的金属结构件中产生的涡流损耗不均匀，有可能造成这些结构件的局部过热现象。

变压器的容量越大，漏磁场就越强，从而使稳态漏磁场引起的各种附加损耗增加，如设计不当它将造成变压器的局部过热，使变压器的热性能变坏最终导致绝缘材料的热老化与击穿。

温度ansys分析

4 汽车内饰压制成型模具温度场模拟与分析温度在汽车内饰压制成型过程中是一个极其重要的参数，无论是模压料的充模流动阶段还是固化阶段，都是在一定的温度下进行的；如果在充模阶段温度控制的不当将直接影响制品的表面质量和力学性能，具体的说，若模具温度过低则会导致模压料流动性降低，难以充满模腔，若模具温度过高则会引起模压料在模具内未完全成型前就开始固化，并且有可能使制品表面的树脂发生分解，同样得不到合格的汽车内饰制品；另一个重要方面就是在压制的过程中要尽量保持模具温度的均匀分布，如果温度分布不均匀就会导致模压料局部提前固化，还会使制品固化度不均匀甚至发生局部树脂分解，同时也会使得制品脱模后产生较大的翘曲变形。

因此有必要对模具的加热过程及其温度场进行模拟，根据分析结果对模具的加热设计进行优化。

4．1 ANSYS有限元分析软件4．1．1 ANSYS简介ANSYS是一种应用广泛的大型通用有限元软件，具有完备的预处理器和后处理器(又称前处理模块和后处理模块)。

包括简单线性静态分析和复杂非线性动态分析。

可用来求解结构、流体、电力、电磁场及碰撞等问题。

在ANSYS 软件中有五个模块可以进行热分析，如图4．1所示，包括：ANSYS／Multiphysics、ANSYS／Mechanical、ANSYS／Thermal、ANSYS／Flotran和ANSYS／ED。

混凝土箱梁日照温度场温度应力ansys分析结果

应力大小往往会超过列车或汽车荷载效应，特别是横向温度应力对混凝土箱梁纵向裂纹的出现有很大的贡献。

下面首先发几张混凝土箱梁日照温度场ANSYS分析结果的图片，希望对这方面感兴趣的网友在此讨论。

Ⅰ：夏季日照温度场。

由于，桥轴线走向和纬度的关系，腹板在夏季腹板几乎不受日照，因此截面温度梯度主要在竖向。

本箱梁冬季腹板也会受到一定的日照。

3、希望提供你的计算思路，偶们好学习一下。

ansys焊接平板温度分析分解

ANSYS 操作实例问题描述某一平板由钢板和铁板焊接而成，焊接材料为铜，平板尺寸为1*1*0.2，横截结构如图1所示，平板初始温度为800℃，将平板放置于空气中进行冷却，周围空气温度为30℃，对流系数为110W/（2m·℃）。

求10分钟后平板内部的温度场及应力分布。

图 1材料温度℃弹性模量Gpa屈服强度Gpa切变模量Gpa导热系数W/(m·℃)线膨胀系数℃1-比热容J/(kg·℃)密度Kg/3m泊松比钢30 206 1.40 20.666.6 1.06e-5 460 7800 0.3 200 192 1.33 19.8400 175 1.15 18.3600 153 0.92 15.6800 125 0.68 11.2铜30 103 0.9 10.3383 1.75e-5 390 8900 0.3 200 99 0.85 0.98400 90 0.75 0.89600 79 0.62 0.75800 58 0.45 0.52铁30 118 1.04 1.1846.55.87e-645070000.3200 109 1.01 1.02 400 93 0.910.86 600 75 0.76 0.69 800520.560.51表1 材料性能参数求解步骤第一步：建立工作文件名1、选择utility menu/file/change jobname 命令，出现 change jobname 对话框，在文本框中输入工作文件名thermal ，如图2所示单击ok 按钮关闭对话框。

图 22、选择main menu/prefercens 命令，弹出对话框，做如图3改动图 3第二步：定义单元类型1、选择main menu/preprocessor/element type/add/edit/delete 命令2、单击add按钮，在生成对话框中分别选择couple field、vector quad 13,单击ok，如图4所示图 43、单击element type对话框上的options按钮，出现plane13 element type options 对话框，在element degree of freedom k1下拉列表中选择ux uy temp az选项，其余采用默认值，单击ok，如图5所示图 54、单击element type对话框上的add，在出现的对话框中分别选择couple field、scalar brick5选项，单击ok，如图6所示图 65、单击element type对话框中的close按钮，关闭对话框第三步：定义材料性能参数1、选择main menu/ preprocessor / material props /material models 命令2、在出现的对话框中依次双击thermal、conductivity、isotropic选项，在出现的对话框中输入导热系数66.6，单击ok，如图7所示图73、依次双击structral、thermal expansion、secant coefficient、isotropic选项，在出现的ALPX文本框中输入线膨胀系数1.06e-5，单击ok如图8所示图84、依次双击structral、density选项，在dens文本框中输入密度7800，单击ok，如图9所示5、依次双击thermal、specific heat选项，在C文本框中输入比热460，单击ok，如图10所示图106、依次双击structral、linear、elastic、isotropic选项，在出现的对话框中单击add temperature按钮4次，如图11所示设置，单击ok图117、依次双击structral、nolinear、inelastic、rate independent、kinematic hardening plasticcity、mises plasticty、bilinear选项，在出现的对话框中单击add temperature 按钮4次，如图12所示设置，单击ok图128、在define material model behavior对话框中选择material/new model命令，在define material id文本框中输入2，单击ok9、依次双击thermal、conductivity、isotropic选项，在出现的对话框中输入导热系数383，单击ok10、依次双击structral、thermal expansion、secant coefficient、isotropic选项，在出现的ALPX文本框中输入线膨胀系数1.75e-5，单击ok11、依次双击structral、density选项，在dens文本框中输入密度8900，单击ok12、依次双击thermal、specific heat选项，在C文本框中输入比热390，单击ok13、依次双击structral、linear、elastic、isotropic选项，在出现的对话框中单击add temperature按钮4次，如图13所示设置，单击ok图1314、依次双击structral、nolinear、inelastic、rate independent、kinematic hardening plasticcity、mises plasticty、bilinear选项，在出现的对话框中单击add temperature 按钮4次，如图14所示设置，单击ok图1415、在define material model behavior对话框中选择material/new model命令，在define material id文本框中输入3，单击ok16、依次双击thermal、conductivity、isotropic选项，在出现的对话框中输入导热系数46.5，单击ok17、依次双击structral、thermal expansion、secant coefficient、isotropic选项，在出现的ALPX文本框中输入线膨胀系数5.87e-6，单击ok18、依次双击structral、density选项，在dens文本框中输入密度7000，单击ok19、依次双击thermal、specific heat选项，在C文本框中输入比热450，单击ok20、依次双击structral、linear、elastic、isotropic选项，在出现的对话框中单击add temperature按钮4次，如图15所示设置，单击ok图1521、依次双击structral、nolinear、inelastic、rate independent、kinematic hardening plasticcity、mises plasticty、bilinear选项，在出现的对话框中单击add temperature 按钮4次，如图16所示设置，单击ok图1622、在define material model behavior对话框中选择material/exit命令，关闭该对话框第四步：创建几何模型、划分网格1、选择main menu/preprocessor/moeling/create/keypints/in active CS命令，在出现的NPT keypoints number文本框中输入关键点编号1，在X，Y，Z location in active CS 文本框中依次输入关键点坐标0、0、02、单击apply ，依次输入一下关键点编号和坐标：2（0.5,0,0）；3（1,0,0）；4（0,0.2,0）；5（0.4,0.2,0）；6（0.6,0.2,0）；7（1,0.2,0）；3、选择main menu/preprocessor/modeling/create/areas/arbitrary/through KPs 命令，在出现的菜单文本框中输入1,2,5,4，单击apply，在文本框中输入2,3,7,6，单击ok4、选择utility menu/workplane/change active CS to/Global Cylindrical命令，将当前激活坐标系转变为柱坐标系5、选择main menu/preprocessor/modeling/create/lines/lines/in Active Coord命令，在出现的菜单文本框中输入6,5，单击ok6、选择Utility Menu/PlotCtrls/Numbering命令，在出现的菜单中选择Line numbers和Area numbers选项，使其由OFF变为ON，单击ok7、选择main menu/preprocessor/meshing/mesh/areas/arbitrary/By Lines命令，在出现的菜单文本框中输入2,8,9，单击ok8、选择Utility Menu/plot/areas命令，窗口显示生成的平面几何模型，如图17所示图179、选择main menu/preprocessor/meshing/size cntrls/manua size/global/size命令，在出现的菜单SIZE element edge length文本框中输入0.05，单击ok如图18所示图1810、选择main menu/preprocessor/meshing/mesh/areas/mapped/3or4 sided命令，在出现的菜单文本框中输入1,2,3，单击ok11、选择Utility Menu/plot/elements命令，显示窗口显示网格划分结果，如图19所示图1912、选择main menu/preprocessor/meshing/size cntrls/manua size/global/size命令，将SIZE element edge length文本框中输入0.05删掉，在NDIV NO.of element divisions 文本框中输入10，单击ok现Element Extrusion Options对话框，在【type】下拉列表中选择2 solid5，在【MAT】下拉列表中选择1，其余不做改动，单击ok，如图20所示图2014、选择main menu/preprocessor/modeling/operate/extrude/areas/along normal命令，在出现的Extrude Area by菜单文本框中输入1，单击ok，出现Extrude Area along normsl 对话框，在NAREA文本框中输入1，在DIST文本框中输入1，单击ok，如图21所示图21现Element Extrusion Options对话框，在【MAT】下拉列表中选择3，其余不做改变，单击ok，如图22所示图2216、选择main menu/preprocessor/modeling/operate/extrude/areas/along normal命令，在出现的Extrude Area by菜单文本框中输入2，单击ok，出现Extrude Area along normsl 对话框，在NAREA文本框中输入2，在DIST文本框中输入1，单击ok17、选择main menu/preprocessor/modeling/operate/extrude/elem ext opts命令，出现Element Extrusion Options对话框，在【MAT】下拉列表中选择2，其余不做改变，单击ok18、选择main menu/preprocessor/modeling/operate/extrude/areas/along normal命令，在出现的Extrude Area by菜单文本框中输入3，单击ok，出现Extrude Areaalong normsl 对话框，在NAREA文本框中输入3，在DIST文本框中输入-1，单击ok19、拖拉面生成结果如图23所示图23 拖拉生成结果第五步：加载求解1、选择main menu/solution/analysis type/new analysis 命令，出现new analysis 对话框，选择分析类型为transient，如图24所示，单击ok，出现transient analysis 对话框，采用full设置，如图25所示，单击ok图24图252、选择main menu/solution/load step opts/time-frenquenc/time integration/amplitude decay命令，出现time integeration controls 对话框，进行如图26所示设置，单击ok图263、选择main menu/solution/analysis type/sol’n controls 命令，出现solution controls对话框，选择basic选项卡，参照图27进行设置；选择transient选项卡，参照图28进行设置，单击ok图27图284、选择main menu/solution/define loads/apply/structural/temperature/uniform temp 命令，出现uniform temperature 对话框，如图29进行设置，单击ok图295、选择utility menu/select/entities命令，出现select entities对话框，在第1个下拉列表中选择Areas，第2个下拉列表中选择By num/pic 选项，点击Unselect 选项，单击ok，出现Unselect areas 菜单，在文本框中输入6,13，单击ok关闭该菜单6、选择Utility menu/select/entities 命令，出现select entities对话框，在第1个下拉列表中选择Nodes，第2个下拉列表中选择Attached to 选项,在第3个选项中点击Areas,all选项，并选中From Full 选项，单击ok关闭此菜单7、选择Main menu/solution/define loads/apply/thermal/convection/on nodes命令，在出现的菜单中单击Pick all，出现Apply CONV on nodes 对话框，如图30进行设置，单击ok关闭此对话框图308、选择Utility Menu/select/everything命令9、选择Main Menu/solution/solve/Current LS 命令，出现Solve Current Load Step对话框，单击ok，ANSYS开始求解计算第六步：查看求解结果1、选择Manu Menu/general postproc/read results/last Set命令2、选择Manu Menu/general postproc/plot results/contour plot/nodal solu 命令，出现Contour Nodal Solution Data对话框，选择Node Solution/DOF Solution/Nodal Temperature,单击ok，ANSYS显示窗口显示温度场分布等值线图，如图31所示图31 温度等值线图3、选择Manu Menu/general postproc/plot results/contour plot/nodal solu 命令，出现Contour Nodal Solution Data对话框，依次选择Node Solution/DOF Solution/X、Y、Z-Compenent of displancement 选项，单击ok，便可以在ANSYS显示窗口查看X、Y、Z方向上的位移场分布等值线图，如图32、33、34所示图32 X方向位移等值线图图33 Y方向位移等值线图图34 Z方向位移等值线图4、选择Manu Menu/general postproc/plot results/contour plot/nodal solu 命令，出现Contour Nodal Solution Data对话框，依次选择Node Solution/Stress/X、Y、Z-Compenent of stress 选项，单击ok，便可以在ANSYS显示窗口查看X、Y、Z 方向上的应力场分布等值线图，如图35、36、37所示图35 X方向应力等值线图图36 Y方向应力等值线图图37 Z方向应力等值线图5、选择Manu Menu/general postproc/plot results/contour plot/nodal solu 命令，出现Contour Nodal Solution Data对话框，选择Node Solution/Stress/von Mises stress 选项，单击ok，便可以在ANSYS显示窗口查看等效应力场分布等值线图，如图38所示图38 等效应力等值线图6、选择Utilty Menu/File/Exit命令，出现Exit from ANSYS对话框，选中Quit-save everything,单击ok，关闭ANSYS。

四个ANSYS热分析经典例子

实例1：某一潜水艇可以简化为一圆筒，它由三层组成，最外面一层为不锈钢，中间为玻纤隔热层，最里面为铝层，筒内为空气，筒外为海水，求内外壁面温度及温度分布。

几何参数：筒外径30 feet总壁厚2 inch不锈钢层壁厚0.75inch玻纤层壁厚1 inch铝层壁厚0.25inch筒长200 feet导热系数不锈钢8.27BTU/hr.ft.o F玻纤0.028 BTU/hr.ft.o F铝117.4 BTU/hr.ft.o F边界条件空气温度70 o F海水温度44.5 o F空气对流系数2.5 BTU/hr.ft2.o F海水对流系数80 BTU/hr.ft2.o F沿垂直于圆筒轴线作横截面，得到一圆环，取其中1 度进行分析，如图示。

/filename，Steady1/title，Steady-state thermal analysis of submarine/units，BFTRo=15 !外径(ft)Rss=15-(0.75/12) !不锈钢层内径ft)Rins=15-(1.75/12) !玻璃纤维层内径(ft)Ral=15-(2/12) !铝层内径(ft)Tair=70 !潜水艇内空气温度Tsea=44.5 !海水温度Kss=8.27 !不锈钢的导热系数(BTU/hr.ft.oF)Kins=0.028 !玻璃纤维的导热系数(BTU/hr.ft.oF)Kal=117.4 !铝的导热系数(BTU/hr.ft.oF)Hair=2.5 !空气的对流系数(BTU/hr.ft2.oF)Hsea=80 !海水的对流系数(BTU/hr.ft2.oF)prep7et，1，plane55 !定义二维热单元mp，kxx，1，Kss !设定不锈钢的导热系数mp，kxx，2，Kins !设定玻璃纤维的导热系数mp，kxx，3，Kal !设定铝的导热系数pcirc，Ro，Rss，-0.5，0.5 !创建几何模型pcirc，Rss，Rins，-0.5，0.5pcirc，Rins，Ral，-0.5，0.5aglue，allnumcmp，arealesize，1，，，16 !设定划分网格密度lesize，4，，，4lesize，14，，，5lesize，16，，，2Mshape，2 !设定为映射网格划分mat，1amesh，1mat，2amesh，2mat，3amesh，3/SOLUSFL，11，CONV，HAIR，，TAIR !施加空气对流边界SFL，1，CONV，HSEA，，TSEA !施加海水对流边界SOLVE/POST1PLNSOL ！输出温度彩色云图finish实例2一圆筒形的罐有一接管，罐外径为3 英尺，壁厚为0.2 英尺，接管外径为0.5 英尺，壁厚为0.1英尺，罐与接管的轴线垂直且接管远离罐的端部。