ansys溫度荷載實(shí)例分析命令流

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!設(shè)定碎的導(dǎo)熱系數(shù)!定義碎的密度

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*畫網(wǎng)格

5、aatt,1,1,AESIZE,all,0.02MSHAPE,0,2DMSHKEY,2AMESH,ALL!將所有面網(wǎng)格劃分時(shí)的單元尺寸設(shè)置為0.02m!利用四邊形單元進(jìn)行網(wǎng)格劃分!采用映射網(wǎng)格劃分單元!對(duì)所有網(wǎng)格進(jìn)行劃分單元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*計(jì)算日照時(shí)間內(nèi)，即14個(gè)小時(shí)內(nèi)的空氣的溫度值*DIM,t_outdoor,13*DO,t,6,18*日照強(qiáng)度!定義室外空氣溫度數(shù)組為t_outdoor!給室外空氣溫度數(shù)組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!*綜合換熱系數(shù)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!設(shè)置求解類型瞬態(tài)熱分析!指定瞬態(tài)分析的求解方法為完全法!打開時(shí)間積分效應(yīng)!設(shè)置構(gòu)件的初始溫度為20攝氏度!將除SVAR和LOCI以外的所有計(jì)算數(shù)據(jù)寫入數(shù)據(jù)庫(kù)和文件中l(wèi)=0*do,m,1,13,1*do,r,3600,3600,3600time,(m-1)*3600+rl=l+1nsubst,1,100,1!設(shè)置每個(gè)荷載的子步數(shù)數(shù)為1,最大值為100,最小值為1autots,on!打開自動(dòng)時(shí)間步長(zhǎng)跟蹤eqslv,JCG!指定方程求解器為JCGkbc,0!使用遞增方式加載線性荷載!施加對(duì)流荷載!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)!給面溫度!施加太陽(yáng)輻射荷載!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定義當(dāng)前荷載步的對(duì)流換熱系數(shù)及周圍環(huán)境21定義當(dāng)前荷載步的對(duì)流換熱系數(shù)及周圍環(huán)境18定義當(dāng)前荷載步的對(duì)流換熱系數(shù)及周圍環(huán)境10定義當(dāng)前荷載步的對(duì)流換熱系數(shù)及周圍環(huán)境給體1施加當(dāng)前荷載步的等效生熱率給體4施加當(dāng)前荷載步的等效生熱率給體3施加當(dāng)前荷載步的等效生熱率給體2施加當(dāng)前荷載步的等效生熱率allssolve*enddo*enddo進(jìn)入熱應(yīng)力求解*/PREP7ETCHG,TT

19、S!熱至U結(jié)構(gòu)分析轉(zhuǎn)換/SOLantype,4!瞬態(tài)timint,1,struct!Turnsontransienteffectstimint,0,thermtimint,0,magtimint,0,elecautots,on!打開自動(dòng)時(shí)間步長(zhǎng)，有利于非線性收斂tref,20trnopt,full!完全瞬態(tài)nropt,full!設(shè)置牛頓-拉普森選項(xiàng)kbc,0!漸變荷載nlgeom,on!打開大變形效應(yīng)LUMPM,0PSTRES,ONNROPT,FULL,on!定義熱應(yīng)力計(jì)算參考溫度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',''!讀入熱分析的計(jì)算結(jié)果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讀入熱分析的計(jì)算結(jié)果*ENDDO*DO,I,1,300LDREAD,TEMP,635

21、5+50*i,'model_BD1',Hh',''!OUTRES,ALL,ALL,TIME,6355+50*i讀入熱分析的計(jì)算結(jié)果DELTIM,50,50,50SOLVE*ENDDO*DO,I,1,300LDREAD,TEMP,21355+100*i,'model_BD1','rth',''!OUTRES,ALL,ALL,TIME,21355+100*i讀入熱分析的計(jì)算結(jié)果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!選擇編號(hào)為1的體元素eslv,r!選擇當(dāng)前所選體元素上的所有單元nsle,r!選擇當(dāng)前所選單元上的所有節(jié)點(diǎn)*get,node_total,node,0,count!*get,node_num1,node,0,num,min最小編號(hào)，并賦值給變量node_num1(1)*do,n,2,node_total,1*get,node_num1(n),node,node_num1(n-1),nxth值給數(shù)組node_nu

24、m1*enddo*do,f,1,node_total,1*get,node_temp1(f),node,node_num1(f),temp值，并賦值給數(shù)組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)輻射凈強(qiáng)度*elsevalue_4=value_3*(t_earth(m)*4-tarea*4)*(1-cos_ph(1)/2)射凈強(qiáng)度*endif*set,long_wave(1,l),value_4/0.03的等效生熱率提取當(dāng)前所選節(jié)點(diǎn)的個(gè)數(shù)!提取當(dāng)前所選節(jié)點(diǎn)中節(jié)點(diǎn)的將當(dāng)前所選擇的節(jié)點(diǎn)的編號(hào)賦提取當(dāng)前所選節(jié)點(diǎn)的節(jié)點(diǎn)溫度!計(jì)算當(dāng)前所選節(jié)點(diǎn)的溫度值計(jì)算當(dāng)前所選節(jié)點(diǎn)溫度值的平!計(jì)算當(dāng)前時(shí)刻的天空溫度值，!計(jì)算與天空之間值長(zhǎng)波!計(jì)算與地面之間值長(zhǎng)波輻!計(jì)算長(zhǎng)波輻射凈強(qiáng)度allselvsel,s,4!選擇編號(hào)為4的體元素提取當(dāng)前所選節(jié)點(diǎn)的個(gè)數(shù)!提取當(dāng)前所選節(jié)點(diǎn)中節(jié)!將當(dāng)前所選擇的節(jié)點(diǎn)的!

26、提取當(dāng)前所選節(jié)點(diǎn)的節(jié)計(jì)算當(dāng)前所選節(jié)點(diǎn)的溫!計(jì)算當(dāng)前所選節(jié)點(diǎn)溫度!計(jì)算當(dāng)前時(shí)刻的天空溫!計(jì)算與天空之間值提取當(dāng)前所選節(jié)點(diǎn)的個(gè)數(shù)!提取當(dāng)前所選節(jié)點(diǎn)中節(jié)將當(dāng)前所選擇的節(jié)點(diǎn)的編eslv,r!選擇當(dāng)前所選體元素上的所有單元nsle,r!選擇當(dāng)前所選單元上的所有節(jié)點(diǎn)*get,node_total,node,0,count*get,node_num1,node,0,num,min點(diǎn)的最小編號(hào)，并賦值給變量node_num1(1)*do,n,2,node_total,1*get,node_num1(n),node,node_num1(n-1),nxth編號(hào)賦值給數(shù)組node_num1*enddo*do,f,1

27、,node_total,1*get,node_temp1(f),node,node_num1(f),temp點(diǎn)溫度值，并賦值給數(shù)組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)長(zhǎng)波輻射凈強(qiáng)度*elsevalue_4=value_3*(t_earth(m)*4-tarea*4)*(1-cos_ph(2)/2)波輻射凈強(qiáng)度*endif*set,long_wave(2,l),value_4/0.03強(qiáng)度的等效生熱率allselvsel,s,3!選擇編號(hào)為3的體元素eslv,r!選擇當(dāng)前所選體元素上的所有單元nsle,r!選擇當(dāng)前所選單元上的所有節(jié)點(diǎn)*get,node_total,node,0,count*get,node_num1(1),node,0,num,min點(diǎn)的最小編號(hào)，并賦值給變量node_num1(1)*do,n,2,node_total,1*get,no

29、de_num1(n),node,node_num1(n-1),nxth號(hào)賦值Z數(shù)組node_num1*enddo*do,f,1,node_total,1!提取當(dāng)前所選節(jié)點(diǎn)的節(jié)點(diǎn)!計(jì)算當(dāng)前所選節(jié)點(diǎn)的溫!計(jì)算當(dāng)前所選節(jié)點(diǎn)溫度值!計(jì)算當(dāng)前時(shí)刻的天空溫!計(jì)算與天空之間值長(zhǎng)波!計(jì)算與地面之間值長(zhǎng)波輻!計(jì)算長(zhǎng)波輻射凈強(qiáng)度提取當(dāng)前所選節(jié)點(diǎn)的個(gè)!提取當(dāng)前所選節(jié)點(diǎn)中將當(dāng)前所選擇的節(jié)點(diǎn)的提取當(dāng)前所選節(jié)點(diǎn)的節(jié)計(jì)算當(dāng)前所選節(jié)點(diǎn)的*get,node_temp1(f),node,node_num1,temp溫度值，并賦值給數(shù)組node_temp1*enddototal=0*do,f,1,node_total,1total=total+node_temp1(f)度值之和，并賦值給變量total*enddotarea=total/node