ansys1__热分析
- 格式:ppt
- 大小:438.00 KB
- 文档页数:27


ansys项目热分析部分FINISH/CLEAR/FILNAME,THERMOANALYSIS,1/TITLE,The Analysis of THERMALKEYW,THERMAL,1/PREP7K1=50ET,1,SOLID70 !定义工件的单元类型EMUNIT,MKS !采用国际单位制mptemp,1,0,100,200,300,400,500 !定义材料属性温度范围mptemp,,1200,1250,1350,1400,1450,1500mptemp,,1505,1705,1905,2105,2305,2500mpdata,dens,1,1,7820,7800,7800,7800,7800,7800 !定义材料密度(kg/m3)mpdata,dens,1,,7800,7800,7800,7800,7800,7800mpdata,dens,1,,7800,7800,7800,7800,7800,7800mpdata,kxx,1,1,52,50.7,48.6,46.1,42.3,38.9 !热传导系数(W/(m*K))mpdata,kxx,1,,30,30,30,30,30,30mpdata,kxx,1,,30,30,30,30,30,30mpdata,c,1,1,450,469,481,508.5,536,569 !考虑相变潜热的比热容(J/kg*K)mpdata,c,1,,700,2172,5116,6589,8061,9533mpdata,c,1,,9533,7757,5982,4206,2431,700L=0.3 !加热线长度LSIZE=0.01 !单元尺寸V=0.3/30 !速度TINC=LSIZE/V !单步时间长度,每个单元加热时间A=2.175E6 !整个丙烯燃烧过程的总热值(J/mol)eita=0.2575 !此条件下的热效率q_c=0.000182518 !丙烯流量(m^3/s)r0=0.004245 !热源有效加热半径(m)qmax=3*A*eita*q_c/3.1415926/r0/r0 !加热斑点中心最大热流密度(W*m/mol)tt=L/V !总时间建模,划分网格!创建几何模型X1=500e-3 !设置工件尺寸及位置(m)Y1=500e-3t=14e-3 !设置工件厚度BLOCK,0,X1,0,Y1,0,t !工件模型ESIZE,0.01VSEL,S,VOLU,,1 !划分工件网格MSHAPE,0,3DMSHKEY,1VMESH,ALLFINISH!划分网格LSEL,S,LOC,Z,t/2LSEL,R,LOC,X,0LESIZE,ALL,0.003,,,2.5LSEL,S,LOC,Z,t/2LSEL,R,LOC,X,X1/2,X1LESIZE,ALL,0.003,,,0.4LSEL,S,LOC,Z,tLSEL,R,LOC,X,0LESIZE,ALL,0.01LSEL,S,LOC,Z,tLSEL,R,LOC,Y,0LESIZE,ALL,0.01/SOLUANTYPE,TRANS,NEWTUNIF,20 !初始温度20度tref,25nropt,auto !打开牛顿莱普生自适应法autos,onkbc,0pred,on !打开非线性分析的预测lnsrch,on !打开先搜索MAX_TIME=1+tt/TINCMAX_X=1+L/LSIZEMAX_Y=1+L/LSIZE*DIM,HEATFLUX,TABLE,MAX_X,MAX_Y,MAX_TIME,X,Y,TIME *DO,K,1,MAX_TIME,1*DO,I,1,MAX_X,1HEATFLUX(I,0,K)=(I-1)*LSIZE*ENDDO*DO,J,1,MAX_Y,1HEATFLUX(0,J,K)=(J-1)*LSIZE*ENDDO*DO,I,1,MAX_X,1*DO,J,1,MAX_Y,1XCENTER=X1/2YCENTER=V*(K-1)*TINCr=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE-YCENTER)**2) !离加热斑点中心的距离(m)*IF,r,LE,r0,THENHEATFLUX(I,J,K)=qmax*exp(-3*r*r/r0/r0) !离加热斑点中心r处的热流密度(W*m/mol)*ELSEHEATFLUX(I,J,K)=0*ENDIF*ENDDO*ENDDOHEATFLUX(0,0,K)=(K-1)*TINC*ENDDO!开始计算TIMINT,ONKBC,1*DO,TM,1E-6,tt,TINCASEL,S,LOC,Z,tSFA,ALL,1,HFLUX,%HEATFLUX%NSUBST,1TIME,TMAUTOTS,0ALLSEL,ALLOUTRES,ALL,ALLSOLVEWPCSYS,-1/VIEW,1,1,1,1/PSF,HFLUX,,2 !显示表面负载/REPLOT*ENDDOFINISHASEL,S,LOC,Z,tSFADELE,ALL,HFLUXALLSEL,ALL/SOLUANTYPE,TRANStrnopt,full !瞬态分析选项,完全分析nropt,full,,on !定义完全牛顿--拉普森方法,激活自适应下降(默认)pred,on !打开预测校正tintp,0.005,,,1,0.5,0.2 !定义瞬态综合参数tref,20 !参考温度20摄氏度t=0 !求解时间初始值dt0=1e-6 !建立初始条件的小时间段dt1=0.5 !起始加热点对应的时间LSIZE=0.01l=0.3 !加热线长度(m)lnum=l/LSIZE !加热线网格划分份数V=l/145 !热源移动速度TINC=LSIZE/V !单步时间长度TIME,dt0TIMINT,OFFKBC,1NSUBST,5ASEL,S,LOC,Z,tNSLA,S,1ESLN,S,0*GET,NEMAX,ELEM,,NUM,MAX*GET,NEMIN,ELEM,,NUM,MINESEL,ALLIC,ALL,TEMP,20ALLSEL,ALLOUTRES,ALL,ALLSOLVEA=2.175E6 !整个丙烯燃烧过程的总热值(J/mol)eita=0.2575 !此条件下的热效率q_c=0.000182518 !丙烯流量(m^3/s)r0=0.04245e-3 !热源有效加热半径(m)qmax=3*A*eita*q_c/3.1415926/r0/r0 !加热斑点中心最大热流密度(W*m/mol)*do,im,1,1,1b=im*LSIZEASEL,S,LOC,Z,tNSLA,S,1ESLN,S,0*do,i,NEMIN,NEMAX,1*if,esel(i),eq,1,then !如果该单元在上面的集合里,则xe=centrx(i) !读取该单元的中心坐标ye=centry(i)ze=centrz(i)*if,abs(xe-X1/2),le,r0,thenr=abs(sqrt((xe-X1/2)*(xe-X1/2)+(ye-b)*(ye-b)))!该单元中心距离热源中心*if,r,le,r0,then !该单元中心在加热半径范围内qr=qmax*exp(-3*r*r/r0/r0) !该单元中心处的热流大小sfe,i,nmface(i),hflux,,qr !在每个单元所指定的面上施加热源载荷*endif*elseqr=0*endif*endif*enddotm=dt1+im*TINCtime,tmtimint,onautos,onKBC,1deltim,0.01allsel,alloutres,all,allsolveASEL,S,LOC,Z,tNSLA,S,1ESLN,S,0*do,i,1,6sfedele,all,i,hflux !删除每个单元六个面上的热载荷*enddo*enddo热流密度热流密度也是一种面荷载。
第一章简介一、热分析的目的热分析用于计算一个系统或部件的温度分布及其它热物理参数,如热量的获取或损失、热梯度、热流密度(热通量〕等。
热分析在许多工程应用中扮演重要角色,如内燃机、涡轮机、换热器、管路系统、电子元件等。
二、ANSYS的热分析∙在ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Thermal、ANSYS/FLOTRAN、ANSYS/ED五种产品中包含热分析功能,其中ANSYS/FLOTRAN不含相变热分析。
∙ANSYS热分析基于能量守恒原理的热平衡方程,用有限元法计算各节点的温度,并导出其它热物理参数。
∙ANSYS热分析包括热传导、热对流及热辐射三种热传递方式。
此外,还可以分析相变、有内热源、接触热阻等问题。
三、ANSYS 热分析分类∙稳态传热:系统的温度场不随时间变化∙瞬态传热:系统的温度场随时间明显变化四、耦合分析∙热-结构耦合∙热-流体耦合∙热-电耦合∙热-磁耦合∙热-电-磁-结构耦合等第二章 基础知识一、符号与单位2-℃ 二、传热学经典理论回顾热分析遵循热力学第一定律,即能量守恒定律:●对于一个封闭的系统(没有质量的流入或流出〕式中: Q —— 热量;W —— 作功;∆U ——系统内能; ∆KE ——系统动能; ∆PE ——系统势能;● 对于大多数工程传热问题:0==PE KE ∆∆; ● 通常考虑没有做功:0=W , 则:U Q ∆=;●对于稳态热分析:0=∆=U Q ,即流入系统的热量等于流出的热量; ●对于瞬态热分析:dtdUq =,即流入或流出的热传递速率q 等于系统内能的变化。
三、热传递的方式1、热传导热传导可以定义为完全接触的两个物体之间或一个物体的不同部分之间由于温度梯度而引起的内能的交换。
热传导遵循付里叶定律:dxdTkq -='',式中''q 为热流密度(W/m 2),k 为导热系数(W/m-℃),“-”表示热量流向温度降低的方向。