ansys温度荷载实例分析命令流

1、!XXXXXXXXXXXXXXXXXX!定义单元类型!XXXXXXXXXXXXXXXXXXFINI/CLE/TITLE,QSHangPylonTest/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,kzz,1,3.3!mp,dens,1,2400mp,c,1,925mp,prxy,1,0.2mp,ex,1,3.45e10mp,alpx,1,1e-5mp,kyy,1,3.3!设定碎的导热系数!定义碎的密度

2、!定义碎的比热!泊松比C50!混凝土的弹性模量C50!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,5loca

3、l,11,1,21.2,0,0,4/2.5,1l,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,1l,10,11CSYS,0l,11,6al,all*plane2asel,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,1

4、7,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*画网格

5、aatt,1,1,AESIZE,all,0.02MSHAPE,0,2DMSHKEY,2AMESH,ALL!将所有面网格划分时的单元尺寸设置为0.02m!利用四边形单元进行网格划分!采用映射网格划分单元!对所有网格进行划分单元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.32

6、lsel,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,l

7、oc,x,-100,0cm,Nligne,linensll,scm,Nnode,nodeALLSTYPE,2ESURF*计算日照时间内，即14个小时内的空气的温度值*DIM,t_outdoor,13*DO,t,6,18*日照强度!定义室外空气温度数组为t_outdoor!给室外空气温度数组t_outdoor赋值*SET,t_outdoor(t-5),31.5+5.5*sin(3.1415926535898*(t-9)/12)*ENDDODIM,W_insolation,13DIM,N_insolation,13!北边日照28E_insolation(1)=240$S_insolation(1)

8、=8933W_insolation(1)E_insolation(2)=521$S_insolation(2)=80$W_insolation(2)N_insolation(2)=160E_insolation(3)=629$S_insolation(3)=115$W_insolation(3)=N_insolation(3)=166E_insolation(4)=640$S_insolation(4)=148$W_insolation(4)N_insolation(4)=152E_insolation(5)=536$S_insolation(5)=184$W_insolation(5)N_i

9、nsolation(5)=142E_insolation(6)=374$S_insolation(6)=229$W_insolation(6)N_insolation(6)=153E_insolation(7)=178$S_insolation(7)=240$W_insolation(7)N_insolation(7)=153E_insolation(8)=178$S_insolation(8)=229$W_insolation(8)N_insolation(8)=153E_insolation(9)=165$S_insolation(9)=184$W_insolation(9)N_insol

10、ation(9)=142E_insolation(10)=148$S_insolation(10)=148$W_insolation(10)N_insolation(10)=152E_insolation(11)=115$S_insolation(11)=115$W_insolation(11)N_insolation(11)=166E_insolation(12)=80$S_insolation(12)=80$W_insolation(12)N_insolation(12)=160E_insolation(13)=33$S_insolation(13)=33$W_insolation(13)

11、N_insolation(13)=89$N_insolation(1)$99$550=54169127142153153322461448206!*综合换热系数h=12.47+3.33*1.5*综合大气温度DIM,Et_synthetiser,13DIM,St_synthetiser,13DIM,Wt_synthetiser,13DIM,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)+

12、(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,transtrnopt,fulltimint,ontunif,27.6!tref,20outres,all!antype,4!autots,on!trnopt,full!lnsrch,on!outres,basic,last!outpr,

13、basic,last!设置求解类型瞬态热分析!指定瞬态分析的求解方法为完全法!打开时间积分效应!设置构件的初始温度为20摄氏度!将除SVAR和LOCI以外的所有计算数据写入数据库和文件中l=0*do,m,1,13,1*do,r,3600,3600,3600time,(m-1)*3600+rl=l+1nsubst,1,100,1!设置每个荷载的子步数数为1,最大值为100,最小值为1autots,on!打开自动时间步长跟踪eqslv,JCG!指定方程求解器为JCGkbc,0!使用递增方式加载线性荷载!施加对流荷载!alls!nsel,s,Enode!d,all,temp,Et_synthetis

14、er(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_synt

15、hetiser(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,

16、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,

17、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)!给面温度!sfa,21,conv,10,t_outdoor(m)!给面温度!sfa,18,conv,10,t_outdoor(m)!给面温度!sfa,10,conv,10,t_outdoor(m)!给面温度!施加太阳辐射荷载!bfv,1,HGEN,heat_eq(1,m)+long_wave(1)!bfv,4,HGEN,heat_eq(2

18、,m)+long_wave(2)!bfv,3,HGEN,heat_eq(3,m)+long_wave(3)!bfv,2,HGEN,heat_eq(4,m)+long_wave(4)!bfe,all,hgen,1,1e105定义当前荷载步的对流换热系数及周围环境21定义当前荷载步的对流换热系数及周围环境18定义当前荷载步的对流换热系数及周围环境10定义当前荷载步的对流换热系数及周围环境给体1施加当前荷载步的等效生热率给体4施加当前荷载步的等效生热率给体3施加当前荷载步的等效生热率给体2施加当前荷载步的等效生热率allssolve*enddo*enddo进入热应力求解*/PREP7ETCHG,TT

19、S!热至U结构分析转换/SOLantype,4!瞬态timint,1,struct!Turnsontransienteffectstimint,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

20、,UY,0DL,ALL,UZ,0*DO,I,1,1355LDREAD,TEMP,1*i,'model_BD1'Eh',''!读入热分析的计算结果OUTRES,ALL,ALL,TIME,1*iDELTIM,1,1,1SOLVE*ENDDO*DO,I,1,500LDREAD,TEMP,1355+10*i,'model_BD1'Eh',''!OUTRES,ALL,ALL,TIME,1355+10*iDELTIM,10,10,10SOLVE读入热分析的计算结果*ENDDO*DO,I,1,300LDREAD,TEMP,635

21、5+50*i,'model_BD1',Hh',''!OUTRES,ALL,ALL,TIME,6355+50*i读入热分析的计算结果DELTIM,50,50,50SOLVE*ENDDO*DO,I,1,300LDREAD,TEMP,21355+100*i,'model_BD1','rth',''!OUTRES,ALL,ALL,TIME,21355+100*i读入热分析的计算结果DELTIM,100,100,100SOLVE*ENDDOSET,LASTCSYS,0PATH,zdirec,2,50PPATH,1,0

22、,0,0.17PPATH,2,0,0.7,0.17PDEF,S,xPLPATH,Sx/IMAGE,SAVE,zpath1,BMPSET,LASTCSYS,0PATH,zdirec,2,50PPATH,1,-1.72,0.4,0.17PPATH,2,1,0.4,0.17PDEF,S,xPLPATH,Sx/IMAGE,SAVE,zpath2,BMPSET,LASTCSYS,0PATH,zdirec,2,50PPATH,1,0,0.7,0PPATH,2,0,0.7,0.35PDEF,S,zPLPATH,Sz/IMAGE,SAVE,zpath3,BMPSET,LASTCSYS,0PATH,zdirec

23、,2,50PPATH,1,-1.72,0.7,0.17PPATH,2,1,0.7,0.17PDEF,S,zPLPATH,Sz/IMAGE,SAVE,zpath4,BMPvsel,s,1!选择编号为1的体元素eslv,r!选择当前所选体元素上的所有单元nsle,r!选择当前所选单元上的所有节点*get,node_total,node,0,count!*get,node_num1,node,0,num,min最小编号，并赋值给变量node_num1(1)*do,n,2,node_total,1*get,node_num1(n),node,node_num1(n-1),nxth值给数组node_nu

24、m1*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

25、(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的体元素提取当前所选节点的个数!提取当前所选节点中节!将当前所选择的节点的!

26、提取当前所选节点的节计算当前所选节点的温!计算当前所选节点温度!计算当前时刻的天空温!计算与天空之间值提取当前所选节点的个数!提取当前所选节点中节将当前所选择的节点的编eslv,r!选择当前所选体元素上的所有单元nsle,r!选择当前所选单元上的所有节点*get,node_total,node,0,count*get,node_num1,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

27、,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)/

28、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,no

29、de_num1(n),node,node_num1(n-1),nxth号赋值Z数组node_num1*enddo*do,f,1,node_total,1!提取当前所选节点的节点!计算当前所选节点的温!计算当前所选节点温度值!计算当前时刻的天空温!计算与天空之间值长波!计算与地面之间值长波辐!计算长波辐射净强度提取当前所选节点的个!提取当前所选节点中将当前所选择的节点的提取当前所选节点的节计算当前所选节点的*get,node_temp1(f),node,node_num1,temp温度值，并赋值给数组node_temp1*enddototal=0*do,f,1,node_total,1total=total+node_temp1(f)度值之和，并赋值给变量total*enddotarea=total/node