| 查看: 1853 | 回復(fù): 13 | |||
chiko75銅蟲 (小有名氣)
|
[求助]
ANSYS中,do循環(huán)中使用生死單元再單點重啟動,如何接著最后一個子步重啟動? 已有1人參與
|
||
|
請問一下,我在ANSYS中試用重啟動功能時,在do循環(huán)里激活單元再單點重啟動,重啟動后發(fā)現(xiàn)所有單元都激活了。沒有重啟動的命令時計算結(jié)果是正常的。請問應(yīng)該怎么改?命令如下: !邊界條件 time,1e-5 solve esel,s,mat,,1 ekill,all allsel *do,i,1,10 locy=i*step_inc local,13,0,0,locy,0 esel,s,loc,x,-width/2,width/2 esel,r,loc,y,-step_inc,step_inc esel,r,loc,z,0,height ealive,all allsel !載荷 allsel solve finish /post1 !后處理中沒有任何命令 /solu resume antype,,rest,,,continue *enddo |
|
你是在知乎上那一篇文件中看到的 "要用solve初始化"?這肯定是錯誤的。如果是 "用 /solu 初始化" 到還勉強說得通。不過,初始化意味著什么你了解嗎?每一次初始化都會使后續(xù)的求解從頭開始,而不管之前計算到了哪一步,這顯然不是你需要的 (從已計算的最后一步的結(jié)構(gòu)繼續(xù)計算)。 進一步的檢查分析問題則比較復(fù)雜了。你可以先執(zhí)行一次循環(huán),并在循環(huán)的開始和重啟動命令之后,以及執(zhí)行 parasove 命令和 pararest 命令前后,分別將死單元放到不同的 component 中并輸出其中的死單元編號到硬盤上,然后比較一下幾組死單元是否相同。如不相同則需要進一步檢查,看看死單元是什么情況下改變的。 至于循環(huán)尾部的 solve 被注解掉,那么就要回到循環(huán)的開始處,經(jīng)過一些操作后才會進行求解,這其中是否有什么問題也需要仔細檢查。 另外,本想試一下你的命令流,那么還需要你提供命令流中兩個 /input 命令所使用的輸入文件,否則無法試算和查找問題。。 |
|
有兩個地方改一下試試: 1 循環(huán)中,第一個 solve 之后增加一個 save 命令。 2 循環(huán)中,后處理之后,回到 /solu,一個 resume 命令是多余的,且會導(dǎo)致其它問題,需要刪除掉。 不過,上面增加 save 命令后,這個 resume 命令也就失去了作用,不刪除也許是可以的。只是為了保險,還是刪除較好。 3 此外,在循環(huán)中第一個 solve 之后,save 之前,還可以增加一個命令“ parsave, all 然后,在后處理之后,回到求解模塊時,再增加一個命令: paeres, new 以保證重新回到求解模塊時模型的狀態(tài)與之前相同。 你先改改試試,如仍有問題,需要你提供完整的命令流,以便分析。 |
|
昨天試算和修改你的命令流,主要是添加了很多輸出中間結(jié)果的命令,進行檢查分析。最終確定重啟動確實會導(dǎo)致生死單元的定義丟失,通過重新定義生死單元可以解決這個問題。 看你的命令流,只是殺死了材料 1 的所有單元,所以定義生死單元還是很方便的。不過,有一個問題我沒有修改,那就是你定義的約束條件好像涉及了少數(shù)材料 1 的單元,而你激活了與約束條件有關(guān)的單元,因而計算過程中死單元的數(shù)量不斷在變化。不知道這是否你的本意。 除了為了檢查中間結(jié)果的部分外,修改的重點有一以幾個: 1 在循環(huán)中,solve 命令之后增加了 save 命令,以保證盤上存放的是最后一個計算成功的子步的模型。 2 在循環(huán)中,后處理之后,回到求解模塊時,恢復(fù)了 resume 命令,以保證模型與上一次求解后的狀態(tài)一致。 3 由于上兩步不能保證死單元不變,所以在 Antype Rest 命令之后重新定義了死單元。 修改后的命令流如下。其中,用 *cfopen 命令打開了一個用于存放中間結(jié)果的文件,計算結(jié)束后可以看一下該文件,可以看到死計算過程單及元數(shù)量的變化。 !*** chiko75 finish /clear /filn, chiko75-03,1 /title, chiko75 -03 /units,si /prep7 ET,1,SOLID70!responsive structural element for solid185 ET,2,SOLID90!responsive structural element for solid186 !et,2,suff_h152 !keyopt,2,4,0 !keyopt,2,5,1 !keyopt,2,8,3 !keyopt,2,9,1 mat3_top=1 mat_h13=2 !*** !*** /input,FeCrBSi_top_ther,dat,E:\acdE\ANSYSdb\tesueship\matpara !*** /input,H13_ther,dat,E:\acdE\ANSYSdb\tesueship\matpara /input, mat3_top,txt !*** 讀入兩個材料屬性文件 /input, mat_h13,txt !*** mpdele,enth,1 mpdele,enth,2 w_matrix=0.05 l_matrix=0.08 h_matrix=0.010 h_clad=0.002 esize_ref=0.001 D=0.002 pi=acos(-1) dists_inc=esize_ref v_scan=10e-3 temp_preh=298.15 count=4 !l_matrix/v_scan /COM,create volume block,-D/2,D/2,0,l_matrix,0,h_clad cm,clad,volu block,-1.5*D,1.5*D,0,l_matrix,0,-h_matrix block,-(1.5*D+2*esize_ref),-1.5*D,0,l_matrix,0,-h_matrix block,1.5*D,(1.5*D+2*esize_ref),0,l_matrix,0,-h_matrix block,-6*D,-(1.5*D+2*esize_ref),0,l_matrix,0,-h_matrix block,(1.5*D+2*esize_ref),6*D,0,l_matrix,0,-h_matrix block,-w_matrix/2,-6*D,0,l_matrix,0,-h_matrix block,6*D,w_matrix/2,0,l_matrix,0,-h_matrix allsel vglue,all vsel,s,loc,x,-w_matrix/2,-6*D vsel,a,loc,x,6*D,w_matrix/2 cm,matrix_hex,volu vsel,s,loc,x,-6*D,6*D vsel,r,loc,z,0,-l_matrix cm,matrix_tet,volu allsel cmgrp,matrix,matrix_hex,matrix_tet allsel lsel,s,loc,x,-1.5*D,1.5*D lesize,all,esize_ref lsel,s,loc,x,-(1.5*D+esize_ref) lsel,a,loc,x,(1.5*D+esize_ref) lsel,u,tan1,z,1 lsel,u,tan1,z,-1 lesize,all,esize_ref lsel,s,loc,x,-4*D,-w_matrix lsel,a,loc,x,4*D,w_matrix lsel,r,loc,z,-h_matrix/2 lesize,all,4*esize_ref lsel,s,loc,x,-4*D,-w_matrix lsel,a,loc,x,4*D,w_matrix lsel,r,tan1,z,0 lesize,all,4*esize_ref mat,1 type,1 mshkey,1 mshape,0 vmesh,clad mat,2 type,1 mshkey,1 mshape,0 vmesh,matrix_hex !vmesh,7 type,2 mshkey,0 mshape,1 !vmesh,5 !vmesh,6 vmesh,matrix_tet numcmp,area vimp,matrix_tet,0,3 numcmp,all !*** 壓縮實體編號,去掉空白的編號 save /solu *cfopen, chiko75-rl3,txt ! 定義一個輸出文件,記錄一些中間結(jié)果 antype,4 timint,1,therm solcontrol,on,, tintp,,,,1,,,,,, thopt,full,0.05,64 neqit,100 trnopt,full nropt,full,,on pred,on,,on cutcontrol,noiterpredict,0 cutcontrol,cutbackfactor,0.5 cutcontrol,plslimit,10, cnvtol,temp,,0.01 cnvtol,heat,,0.01,,1e-6 kbc,0 lumpm,0 lnsrch,1 autots,1 ncnv,2 outres,all,all!basic,last tunif,temp_preh allsel nsel,s,loc,z,-h_matrix sf,all,conv,5,temp_preh allsel time,1e-4 deltim,5e-5,5e-5,5e-5 *vwrite, (/' First solu & solve time=1e-4') !*** 寫中間結(jié)果文件,下同 solve /solu !*** 這個 /solu 是多余的,可以刪除試試,是否影響計算結(jié)果 *vwrite, (/' Second solu - undeleted') esel,s,mat,,1 ekill,all cm, emat_1,elem !*** 為材料 1 死單元建組 !*** 將當(dāng)前死單元數(shù)寫入中間結(jié)果文件,下同 *get,ekilled, elem, 0, count cm, ekilled, elem *vwrite, ekilled (' Ekill all elements for mat 1, count =',f6.0 ) esel,s,live eplot !rescontrol,define,all,all i=1 *do,i,0,count ! count = 4 disy=i*dists_inc t=(disy+dists_inc)/v_scan !!!!! *vwrite, i,t (//' In *do, i=',f5.0, ' time=',f10.8) esel,s,mat,,1 nsle,s,1 nsel,r,loc,y,disy-dists_inc,disy+dists_inc,0.00001 esln,s,1 !*** 可能會包含部分 mat 1 的單元 ealive,all !*** 導(dǎo)致部分死單元被激活 bf,all,hgen,4e10 esel,s,live eplot csys,0 allsel esel,s,live nsle,s,1, nsel,r,ext nsel,u,loc,z,-h_matrix !------改 sf,all,conv,50,temp_preh !------改 time,t nsubst,5,1000,5,1 allsel !求解 SOLVE parsav,all,parareco-i1,txt !** 保存參數(shù) save !*** 求解后將模型存盤 !*** 檢查死單元 cm, e_befor,elem !*** 保存當(dāng)前單元 esel,s,live esel,inve cm, ekilled_2,elem !*** 為死單元建組 *get,ekilled_2, elem, 0, count !!!!! *vwrite, i, ekilled_2 (' After solve, i =',f5.0,' killed elem count = ',f7.0) cmsel, s, e_befor !*** 恢復(fù)原來單元組 /solu !*** 這個 /solu 是多余的,可以刪除試試,是否影響計算結(jié)果 !!!!! *vwrite,i,t (/' 3-th solu i = ',f5.0,' time=',f10.8) parsav,all,parareco-i2,txt !** 再次保存參數(shù),與上一次比較 BFDELE,all,HGEN allsel nsel,u,loc,z,-h_matrix sfdele,all,conv allsel save finish /post1 !etable,temper,temp !esel,s,etab,temper,1160 !cm,temper1,elem finish !save !/filename, !*** 檢查死單元 cm, e_befor,elem !*** 保存當(dāng)前單元 esel,s,live esel,inve cm, ekilled_3,elem !*** 為死單元建組 *get,ekilled_3, elem, 0, count *vwrite, i, ekilled_3 (' After solve, i =',f5.0,' killed elem count = ',f6.0) cmsel, s, e_befor !*** 恢復(fù)原來單元組 /solu resume !**** 恢復(fù)為上面存儲的模型 parres,new,parareco-i1,txt !!!!! *vwrite (/' 4-th solu, before rest ') parsav,all,parareco-i3,txt !** 再次保存參數(shù),與第一次比較 !*** 檢查死單元 cm, e_befor,elem !*** 保存當(dāng)前單元 esel,s,live esel,inve cm, ekilled_4,elem !*** 為死單元建組 *get,ekilled_4, elem, 0, count *vwrite, i, ekilled_4 (' Before rest, i =',f5.0,' killed elem count = ',f6.0) cmsel, s,e_befor !*** 恢復(fù)原來單元組 antype,,rest !*** 檢查死單元 cm, e_befor,elem !*** 保存當(dāng)前單元 esel,s,live esel,inve cm, ekilled_5,elem !*** 為死單元建組 *get,ekilled_5, elem, 0, count *vwrite, i, ekilled_5 (' After rest, i =',f5.0,' killed elem count = ',f6.0) cmsel, s, e_befor !*** 恢復(fù)原來單元組 !*** 如果沒有死單元,則重新定義材料 1 的定義為死單元 *if, ekilled_5,eq,0,then ESEL,S,MAT,,1 ekill,all *get,ekilled_5, elem, 0, count *vwrite, i, ekilled_5 (' Modify: After rest, i =',f5.0,' killed elem count = ',f6.0) cmsel, s, e_befor !*** 恢復(fù)原來單元組 *endif parres,new,parareco-i3,txt !*** 檢查死單元 cm, e_befor,elem !*** 保存當(dāng)前單元 esel,s,live esel,inve cm, ekilled_6,elem !*** 為死單元建組 *get,ekilled_6, elem, 0, count *vwrite, i, ekilled_6 (' After rest¶res, i =',f5.0,' killed elem count = ',f6.0) cmsel, s, e_befor !*** 恢復(fù)原來單元組 !parres,new,parareco,txt !esel,s,,,temper1 !ekill,all ! allsel !kuse,1 !solve allsel *ENDDO *cfclose ! * |
|
你好! 周末摸索了一下你的問題,發(fā)現(xiàn)確實是重啟動導(dǎo)致死單元的定義被刪除,不過,在重啟動命令之后重新定義死單元就可以解決這個問題。 用于摸索你的問題和在重啟動后重新添加死單元定義的命令流如下 ,你可以試試看: !*** chiko75 finish /clear /filn, chiko75-03,1 /title, chiko75 -03 /units,si /prep7 ET,1,SOLID70!responsive structural element for solid185 ET,2,SOLID90!responsive structural element for solid186 !et,2,suff_h152 !keyopt,2,4,0 !keyopt,2,5,1 !keyopt,2,8,3 !keyopt,2,9,1 mat3_top=1 mat_h13=2 !*** !*** /input,FeCrBSi_top_ther,dat,E:\acdE\ANSYSdb\tesueship\matpara !*** /input,H13_ther,dat,E:\acdE\ANSYSdb\tesueship\matpara /input, mat3_top,txt !*** 讀入兩個材料屬性文件 /input, mat_h13,txt !*** mpdele,enth,1 mpdele,enth,2 w_matrix=0.05 l_matrix=0.08 h_matrix=0.010 h_clad=0.002 esize_ref=0.001 D=0.002 pi=acos(-1) dists_inc=esize_ref v_scan=10e-3 temp_preh=298.15 count=4 !l_matrix/v_scan /COM,create volume block,-D/2,D/2,0,l_matrix,0,h_clad cm,clad,volu block,-1.5*D,1.5*D,0,l_matrix,0,-h_matrix block,-(1.5*D+2*esize_ref),-1.5*D,0,l_matrix,0,-h_matrix block,1.5*D,(1.5*D+2*esize_ref),0,l_matrix,0,-h_matrix block,-6*D,-(1.5*D+2*esize_ref),0,l_matrix,0,-h_matrix block,(1.5*D+2*esize_ref),6*D,0,l_matrix,0,-h_matrix block,-w_matrix/2,-6*D,0,l_matrix,0,-h_matrix block,6*D,w_matrix/2,0,l_matrix,0,-h_matrix allsel vglue,all vsel,s,loc,x,-w_matrix/2,-6*D vsel,a,loc,x,6*D,w_matrix/2 cm,matrix_hex,volu vsel,s,loc,x,-6*D,6*D vsel,r,loc,z,0,-l_matrix cm,matrix_tet,volu allsel cmgrp,matrix,matrix_hex,matrix_tet allsel lsel,s,loc,x,-1.5*D,1.5*D lesize,all,esize_ref lsel,s,loc,x,-(1.5*D+esize_ref) lsel,a,loc,x,(1.5*D+esize_ref) lsel,u,tan1,z,1 lsel,u,tan1,z,-1 lesize,all,esize_ref lsel,s,loc,x,-4*D,-w_matrix lsel,a,loc,x,4*D,w_matrix lsel,r,loc,z,-h_matrix/2 lesize,all,4*esize_ref lsel,s,loc,x,-4*D,-w_matrix lsel,a,loc,x,4*D,w_matrix lsel,r,tan1,z,0 lesize,all,4*esize_ref mat,1 type,1 mshkey,1 mshape,0 vmesh,clad mat,2 type,1 mshkey,1 mshape,0 vmesh,matrix_hex !vmesh,7 type,2 mshkey,0 mshape,1 !vmesh,5 !vmesh,6 vmesh,matrix_tet numcmp,area vimp,matrix_tet,0,3 numcmp,all !*** 壓縮實體編號,去掉空白的編號 save /solu *cfopen, chiko75-rl3,txt ! 定義一個輸出文件,記錄一些中間結(jié)果 antype,4 timint,1,therm solcontrol,on,, tintp,,,,1,,,,,, thopt,full,0.05,64 neqit,100 trnopt,full nropt,full,,on pred,on,,on cutcontrol,noiterpredict,0 cutcontrol,cutbackfactor,0.5 cutcontrol,plslimit,10, cnvtol,temp,,0.01 cnvtol,heat,,0.01,,1e-6 kbc,0 lumpm,0 lnsrch,1 autots,1 ncnv,2 outres,all,all!basic,last tunif,temp_preh allsel nsel,s,loc,z,-h_matrix sf,all,conv,5,temp_preh allsel time,1e-4 deltim,5e-5,5e-5,5e-5 *vwrite, (/' First solu & solve time=1e-4') !*** 寫中間結(jié)果文件,下同 solve /solu !*** 這個 /solu 是多余的,可以刪除試試,是否影響計算結(jié)果 *vwrite, (/' Second solu - undeleted') esel,s,mat,,1 ekill,all cm, emat_1,elem !*** 為材料 1 死單元建組 !*** 將當(dāng)前死單元數(shù)寫入中間結(jié)果文件,下同 *get,ekilled, elem, 0, count cm, ekilled, elem *vwrite, ekilled (' Ekill all elements for mat 1, count =',f6.0 ) esel,s,live eplot !rescontrol,define,all,all i=1 *do,i,0,count ! count = 4 disy=i*dists_inc t=(disy+dists_inc)/v_scan !!!!! *vwrite, i,t (//' In *do, i=',f5.0, ' time=',f10.8) esel,s,mat,,1 nsle,s,1 nsel,r,loc,y,disy-dists_inc,disy+dists_inc,0.00001 esln,s,1 !*** 可能會包含部分 mat 1 的單元 ealive,all !*** 導(dǎo)致部分死單元被激活 bf,all,hgen,4e10 esel,s,live eplot csys,0 allsel esel,s,live nsle,s,1, nsel,r,ext nsel,u,loc,z,-h_matrix !------改 sf,all,conv,50,temp_preh !------改 time,t nsubst,5,1000,5,1 allsel !求解 SOLVE parsav,all,parareco-i1,txt !** 保存參數(shù) save !*** 求解后將模型存盤 !*** 檢查死單元 cm, e_befor,elem !*** 保存當(dāng)前單元 esel,s,live esel,inve cm, ekilled_2,elem !*** 為死單元建組 *get,ekilled_2, elem, 0, count !!!!! *vwrite, i, ekilled_2 (' After solve, i =',f5.0,' killed elem count = ',f7.0) cmsel, s, e_befor !*** 恢復(fù)原來單元組 /solu !*** 這個 /solu 是多余的,可以刪除試試,是否影響計算結(jié)果 !!!!! *vwrite,i,t (/' 3-th solu i = ',f5.0,' time=',f10.8) parsav,all,parareco-i2,txt !** 再次保存參數(shù),與上一次比較 BFDELE,all,HGEN allsel nsel,u,loc,z,-h_matrix sfdele,all,conv allsel save finish /post1 !etable,temper,temp !esel,s,etab,temper,1160 !cm,temper1,elem finish !save !/filename, !*** 檢查死單元 cm, e_befor,elem !*** 保存當(dāng)前單元 esel,s,live esel,inve cm, ekilled_3,elem !*** 為死單元建組 *get,ekilled_3, elem, 0, count *vwrite, i, ekilled_3 (' After solve, i =',f5.0,' killed elem count = ',f6.0) cmsel, s, e_befor !*** 恢復(fù)原來單元組 /solu resume !**** 恢復(fù)為上面存儲的模型 parres,new,parareco-i1,txt !!!!! *vwrite (/' 4-th solu, before rest ') parsav,all,parareco-i3,txt !** 再次保存參數(shù),與第一次比較 !*** 檢查死單元 cm, e_befor,elem !*** 保存當(dāng)前單元 esel,s,live esel,inve cm, ekilled_4,elem !*** 為死單元建組 *get,ekilled_4, elem, 0, count *vwrite, i, ekilled_4 (' Before rest, i =',f5.0,' killed elem count = ',f6.0) cmsel, s,e_befor !*** 恢復(fù)原來單元組 antype,,rest !*** 檢查死單元 cm, e_befor,elem !*** 保存當(dāng)前單元 esel,s,live esel,inve cm, ekilled_5,elem !*** 為死單元建組 *get,ekilled_5, elem, 0, count *vwrite, i, ekilled_5 (' After rest, i =',f5.0,' killed elem count = ',f6.0) cmsel, s, e_befor !*** 恢復(fù)原來單元組 !*** 如果沒有死單元,則重新定義材料 1 的定義為死單元 *if, ekilled_5,eq,0,then ESEL,S,MAT,,1 ekill,all *get,ekilled_5, elem, 0, count *vwrite, i, ekilled_5 (' Modify: After rest, i =',f5.0,' killed elem count = ',f6.0) cmsel, s, e_befor !*** 恢復(fù)原來單元組 *endif parres,new,parareco-i3,txt !*** 檢查死單元 cm, e_befor,elem !*** 保存當(dāng)前單元 esel,s,live esel,inve cm, ekilled_6,elem !*** 為死單元建組 *get,ekilled_6, elem, 0, count *vwrite, i, ekilled_6 (' After rest¶res, i =',f5.0,' killed elem count = ',f6.0) cmsel, s, e_befor !*** 恢復(fù)原來單元組 !parres,new,parareco,txt !esel,s,,,temper1 !ekill,all ! allsel !kuse,1 !solve allsel *ENDDO *cfclose |
銅蟲 (小有名氣)
|
謝謝,不過還是不行。具體命令如下 finish /clear /units,si /prep7 ET,1,SOLID70!responsive structural element for solid185 ET,2,SOLID90!responsive structural element for solid186 !et,2,suff_h152 !keyopt,2,4,0 !keyopt,2,5,1 !keyopt,2,8,3 !keyopt,2,9,1 mat3_top=1 mat_h13=2 /input,FeCrBSi_top_ther,dat,E:\acdE\ANSYSdb\tesueship\matpara /input,H13_ther,dat,E:\acdE\ANSYSdb\tesueship\matpara mpdele,enth,1 mpdele,enth,2 w_matrix=0.05 l_matrix=0.08 h_matrix=0.010 h_clad=0.002 esize_ref=0.001 D=0.002 pi=acos(-1) dists_inc=esize_ref v_scan=10e-3 temp_preh=298.15 count=4 !l_matrix/v_scan /COM,create volume block,-D/2,D/2,0,l_matrix,0,h_clad cm,clad,volu block,-1.5*D,1.5*D,0,l_matrix,0,-h_matrix block,-(1.5*D+2*esize_ref),-1.5*D,0,l_matrix,0,-h_matrix block,1.5*D,(1.5*D+2*esize_ref),0,l_matrix,0,-h_matrix block,-6*D,-(1.5*D+2*esize_ref),0,l_matrix,0,-h_matrix block,(1.5*D+2*esize_ref),6*D,0,l_matrix,0,-h_matrix block,-w_matrix/2,-6*D,0,l_matrix,0,-h_matrix block,6*D,w_matrix/2,0,l_matrix,0,-h_matrix allsel vglue,all vsel,s,loc,x,-w_matrix/2,-6*D vsel,a,loc,x,6*D,w_matrix/2 cm,matrix_hex,volu vsel,s,loc,x,-6*D,6*D vsel,r,loc,z,0,-l_matrix cm,matrix_tet,volu allsel cmgrp,matrix,matrix_hex,matrix_tet allsel lsel,s,loc,x,-1.5*D,1.5*D lesize,all,esize_ref lsel,s,loc,x,-(1.5*D+esize_ref) lsel,a,loc,x,(1.5*D+esize_ref) lsel,u,tan1,z,1 lsel,u,tan1,z,-1 lesize,all,esize_ref lsel,s,loc,x,-4*D,-w_matrix lsel,a,loc,x,4*D,w_matrix lsel,r,loc,z,-h_matrix/2 lesize,all,4*esize_ref lsel,s,loc,x,-4*D,-w_matrix lsel,a,loc,x,4*D,w_matrix lsel,r,tan1,z,0 lesize,all,4*esize_ref mat,1 type,1 mshkey,1 mshape,0 vmesh,clad mat,2 type,1 mshkey,1 mshape,0 vmesh,matrix_hex !vmesh,7 type,2 mshkey,0 mshape,1 !vmesh,5 !vmesh,6 vmesh,matrix_tet numcmp,area vimp,matrix_tet,0,3 save /solu antype,4 timint,1,therm solcontrol,on,, tintp,,,,1,,,,,, thopt,full,0.05,64 neqit,100 trnopt,full nropt,full,,on pred,on,,on cutcontrol,noiterpredict,0 cutcontrol,cutbackfactor,0.5 cutcontrol,plslimit,10, cnvtol,temp,,0.01 cnvtol,heat,,0.01,,1e-6 kbc,0 lumpm,0 lnsrch,1 autots,1 ncnv,2 outres,all,all!basic,last tunif,temp_preh allsel nsel,s,loc,z,-h_matrix sf,all,conv,5,temp_preh allsel time,1e-4 deltim,5e-5,5e-5,5e-5 solve /solu esel,s,mat,,1 ekill,all esel,s,live eplot !rescontrol,define,all,all i=1 *do,i,0,count disy=i*dists_inc t=(disy+dists_inc)/v_scan esel,s,mat,,1 nsle,s,1 nsel,r,loc,y,disy-dists_inc,disy+dists_inc,0.00001 esln,s,1 ealive,all bf,all,hgen,4e10 esel,s,live eplot csys,0 allsel esel,s,live nsle,s,1, nsel,r,ext nsel,u,loc,z,-h_matrix !------改 sf,all,conv,50,temp_preh !------改 time,t nsubst,5,1000,5,1 allsel !求解 SOLVE parsav,all,parareco,txt save /solu BFDELE,all,HGEN allsel nsel,u,loc,z,-h_matrix sfdele,all,conv allsel save finish /post1 !etable,temper,temp !esel,s,etab,temper,1160 !cm,temper1,elem finish !save !/filename, /solu !resume parres,new,parareco,txt antype,,rest !parres,new,parareco,txt !esel,s,,,temper1 !ekill,all allsel !kuse,1 !solve allsel *ENDDO 材料物性參數(shù)如下: !mat3_top=1 MPTEMP MPTEMP,1,298.15 MPTEMP,2,373.15 MPTEMP,3,473.15 MPTEMP,4,573.15 MPTEMP,5,673.15 MPTEMP,6,773.15 MPTEMP,7,873.15 MPTEMP,8,973.15 MPTEMP,9,975.33443 MPTEMP,10,998.17961 MPTEMP,11,999.8139 MPTEMP,12,1015.91817 MPTEMP,13,1073.15 MPTEMP,14,1173.15 MPTEMP,15,1373.15 MPTEMP,16,1411.33015 MPTEMP,17,1456.42042 MPTEMP,18,1473.15 MPTEMP,19,1573.15 MPTEMP,20,1618.42234 MPTEMP,21,1658.14359 MPTEMP,22,1673.15 MPTEMP,23,1773.15 MPTEMP,24,1873.15 MPTEMP,25,1973.15 MPTEMP,26,2073.15 MPTEMP,27,2173.15 MPTEMP,28,2273.15 MPTEMP,29,2373.15 MPTEMP,30,2473.15 MPTEMP,31,2573.15 MPTEMP,32,2673.15 MPTEMP,33,2773.15 MPTEMP,34,2873.15 MPTEMP,35,2973.15 MPTEMP,36,3073.15 MPTEMP,37,3173.15 MPTEMP,38,3273.15 FACT_DENS=1000.0 MPDATA,DENS,mat3_top,1,7.02202213*FACT_DENS MPDATA,DENS,mat3_top,2,6.99741469*FACT_DENS MPDATA,DENS,mat3_top,3,6.964275398*FACT_DENS MPDATA,DENS,mat3_top,4,6.930778305*FACT_DENS MPDATA,DENS,mat3_top,5,6.896943469*FACT_DENS MPDATA,DENS,mat3_top,6,6.862790649*FACT_DENS MPDATA,DENS,mat3_top,7,6.828339298*FACT_DENS MPDATA,DENS,mat3_top,8,6.557242353*FACT_DENS MPDATA,DENS,mat3_top,9,6.548547164*FACT_DENS MPDATA,DENS,mat3_top,10,7.178593816*FACT_DENS MPDATA,DENS,mat3_top,11,7.215224808*FACT_DENS MPDATA,DENS,mat3_top,12,7.214457885*FACT_DENS MPDATA,DENS,mat3_top,13,7.1932084*FACT_DENS MPDATA,DENS,mat3_top,14,7.153846107*FACT_DENS MPDATA,DENS,mat3_top,15,7.182897395*FACT_DENS MPDATA,DENS,mat3_top,16,7.130974654*FACT_DENS MPDATA,DENS,mat3_top,17,7.097697794*FACT_DENS MPDATA,DENS,mat3_top,18,7.084619144*FACT_DENS MPDATA,DENS,mat3_top,19,6.981560988*FACT_DENS MPDATA,DENS,mat3_top,20,6.905731826*FACT_DENS MPDATA,DENS,mat3_top,21,6.869300497*FACT_DENS MPDATA,DENS,mat3_top,22,6.857610676*FACT_DENS MPDATA,DENS,mat3_top,23,6.77826579*FACT_DENS MPDATA,DENS,mat3_top,24,6.696569749*FACT_DENS MPDATA,DENS,mat3_top,25,6.612749342*FACT_DENS MPDATA,DENS,mat3_top,26,6.52702854*FACT_DENS MPDATA,DENS,mat3_top,27,6.43962735*FACT_DENS MPDATA,DENS,mat3_top,28,6.350760781*FACT_DENS MPDATA,DENS,mat3_top,29,6.260637928*FACT_DENS MPDATA,DENS,mat3_top,30,6.169461176*FACT_DENS MPDATA,DENS,mat3_top,31,6.077425517*FACT_DENS MPDATA,DENS,mat3_top,32,5.984717994*FACT_DENS MPDATA,DENS,mat3_top,33,5.891517235*FACT_DENS MPDATA,DENS,mat3_top,34,5.797993114*FACT_DENS MPDATA,DENS,mat3_top,35,5.704306492*FACT_DENS MPDATA,DENS,mat3_top,36,5.610609061*FACT_DENS MPDATA,DENS,mat3_top,37,5.517043266*FACT_DENS MPDATA,DENS,mat3_top,38,5.423742303*FACT_DENS MPTEMP MPTEMP,1,298.15 MPTEMP,2,373.15 MPTEMP,3,473.15 MPTEMP,4,573.15 MPTEMP,5,673.15 MPTEMP,6,773.15 MPTEMP,7,873.15 MPTEMP,8,973.15 MPTEMP,9,1073.15 MPTEMP,10,1173.15 MPTEMP,11,1373.15 MPTEMP,12,1473.15 MPTEMP,13,1573.15 MPTEMP,14,1673.15 MPTEMP,15,1773.15 MPTEMP,16,1873.15 MPTEMP,17,1973.15 MPTEMP,18,2073.15 MPTEMP,19,2173.15 MPTEMP,20,2273.15 MPTEMP,21,2373.15 MPTEMP,22,2473.15 MPTEMP,23,2573.15 MPTEMP,24,2673.15 MPTEMP,25,2773.15 MPTEMP,26,2873.15 MPTEMP,27,2973.15 MPTEMP,28,3073.15 MPTEMP,29,3173.15 MPTEMP,30,3273.15 FACT_C=1000.0 MPDATA,C,mat3_top,1,0.52752*FACT_C MPDATA,C,mat3_top,2,0.55072*FACT_C MPDATA,C,mat3_top,3,0.57544*FACT_C MPDATA,C,mat3_top,4,0.59716*FACT_C MPDATA,C,mat3_top,5,0.61792*FACT_C MPDATA,C,mat3_top,6,0.63911*FACT_C MPDATA,C,mat3_top,7,0.7672*FACT_C MPDATA,C,mat3_top,8,0.8071*FACT_C MPDATA,C,mat3_top,9,0.74449*FACT_C MPDATA,C,mat3_top,10,0.76564*FACT_C MPDATA,C,mat3_top,11,1.72865*FACT_C MPDATA,C,mat3_top,12,1.35303*FACT_C MPDATA,C,mat3_top,13,2.05203*FACT_C MPDATA,C,mat3_top,14,0.84526*FACT_C MPDATA,C,mat3_top,15,0.88219*FACT_C MPDATA,C,mat3_top,16,0.90248*FACT_C MPDATA,C,mat3_top,17,0.9105*FACT_C MPDATA,C,mat3_top,18,0.91712*FACT_C MPDATA,C,mat3_top,19,0.92188*FACT_C MPDATA,C,mat3_top,20,0.92292*FACT_C MPDATA,C,mat3_top,21,0.92367*FACT_C MPDATA,C,mat3_top,22,0.92438*FACT_C MPDATA,C,mat3_top,23,0.92507*FACT_C MPDATA,C,mat3_top,24,0.92574*FACT_C MPDATA,C,mat3_top,25,0.9264*FACT_C MPDATA,C,mat3_top,26,0.92704*FACT_C MPDATA,C,mat3_top,27,0.92762*FACT_C MPDATA,C,mat3_top,28,0.9282*FACT_C MPDATA,C,mat3_top,29,0.92876*FACT_C MPDATA,C,mat3_top,30,0.92927*FACT_C MPTEMP MPTEMP,1,298.15 MPTEMP,2,373.15 MPTEMP,3,473.15 MPTEMP,4,573.15 MPTEMP,5,673.15 MPTEMP,6,773.15 MPTEMP,7,873.15 MPTEMP,8,973.15 MPTEMP,9,975.33443 MPTEMP,10,998.17961 MPTEMP,11,999.8139 MPTEMP,12,1015.91817 MPTEMP,13,1073.15 MPTEMP,14,1173.15 MPTEMP,15,1373.15 MPTEMP,16,1411.33015 MPTEMP,17,1456.42042 MPTEMP,18,1473.15 MPTEMP,19,1573.15 MPTEMP,20,1618.42234 MPTEMP,21,1658.14359 MPTEMP,22,1673.15 MPTEMP,23,1773.15 MPTEMP,24,1873.15 MPTEMP,25,1973.15 MPTEMP,26,2073.15 MPTEMP,27,2173.15 MPTEMP,28,2273.15 MPTEMP,29,2373.15 MPTEMP,30,2473.15 MPTEMP,31,2573.15 MPTEMP,32,2673.15 MPTEMP,33,2773.15 MPTEMP,34,2873.15 MPTEMP,35,2973.15 MPTEMP,36,3073.15 MPTEMP,37,3173.15 MPTEMP,38,3273.15 FACT_KXX=1.0 MPDATA,KXX,mat3_top,1,10.78639865*FACT_KXX MPDATA,KXX,mat3_top,2,11.37679024*FACT_KXX MPDATA,KXX,mat3_top,3,12.15276873*FACT_KXX MPDATA,KXX,mat3_top,4,12.91691904*FACT_KXX MPDATA,KXX,mat3_top,5,13.67015615*FACT_KXX MPDATA,KXX,mat3_top,6,14.41326944*FACT_KXX MPDATA,KXX,mat3_top,7,15.14695328*FACT_KXX MPDATA,KXX,mat3_top,8,15.82635113*FACT_KXX MPDATA,KXX,mat3_top,9,15.8357737*FACT_KXX MPDATA,KXX,mat3_top,10,17.27579801*FACT_KXX MPDATA,KXX,mat3_top,11,17.36736784*FACT_KXX MPDATA,KXX,mat3_top,12,17.63439611*FACT_KXX MPDATA,KXX,mat3_top,13,18.15612414*FACT_KXX MPDATA,KXX,mat3_top,14,19.05002667*FACT_KXX MPDATA,KXX,mat3_top,15,28.46718607*FACT_KXX MPDATA,KXX,mat3_top,16,29.65452812*FACT_KXX MPDATA,KXX,mat3_top,17,30.72842726*FACT_KXX MPDATA,KXX,mat3_top,18,31.07635294*FACT_KXX MPDATA,KXX,mat3_top,19,33.5424953*FACT_KXX MPDATA,KXX,mat3_top,20,35.06679548*FACT_KXX MPDATA,KXX,mat3_top,21,35.46549703*FACT_KXX MPDATA,KXX,mat3_top,22,35.68068478*FACT_KXX MPDATA,KXX,mat3_top,23,37.11465705*FACT_KXX MPDATA,KXX,mat3_top,24,38.54862933*FACT_KXX MPDATA,KXX,mat3_top,25,39.9826016*FACT_KXX MPDATA,KXX,mat3_top,26,41.41657388*FACT_KXX MPDATA,KXX,mat3_top,27,42.85054615*FACT_KXX MPDATA,KXX,mat3_top,28,44.28451843*FACT_KXX MPDATA,KXX,mat3_top,29,45.71849071*FACT_KXX MPDATA,KXX,mat3_top,30,47.15246298*FACT_KXX MPDATA,KXX,mat3_top,31,48.58643526*FACT_KXX MPDATA,KXX,mat3_top,32,50.02040753*FACT_KXX MPDATA,KXX,mat3_top,33,51.45437981*FACT_KXX MPDATA,KXX,mat3_top,34,52.88835208*FACT_KXX MPDATA,KXX,mat3_top,35,54.32232436*FACT_KXX MPDATA,KXX,mat3_top,36,55.75629664*FACT_KXX MPDATA,KXX,mat3_top,37,57.19026892*FACT_KXX MPDATA,KXX,mat3_top,38,58.62424112*FACT_KXX MPTEMP MPTEMP MPTEMP,1,298.15 MPTEMP,2,373.15 MPTEMP,3,473.15 MPTEMP,4,573.15 MPTEMP,5,673.15 MPTEMP,6,773.15 MPTEMP,7,873.15 MPTEMP,8,973.15 MPTEMP,9,975.33443 MPTEMP,10,998.17961 MPTEMP,11,999.8139 MPTEMP,12,1015.91817 MPTEMP,13,1073.15 MPTEMP,14,1173.15 MPTEMP,15,1373.15 MPTEMP,16,1411.33015 MPTEMP,17,1456.42042 MPTEMP,18,1473.15 MPTEMP,19,1573.15 MPTEMP,20,1618.42234 MPTEMP,21,1658.14359 MPTEMP,22,1673.15 MPTEMP,23,1773.15 MPTEMP,24,1873.15 MPTEMP,25,1973.15 MPTEMP,26,2073.15 MPTEMP,27,2173.15 MPTEMP,28,2273.15 MPTEMP,29,2373.15 MPTEMP,30,2473.15 MPTEMP,31,2573.15 MPTEMP,32,2673.15 MPTEMP,33,2773.15 MPTEMP,34,2873.15 MPTEMP,35,2973.15 MPTEMP,36,3073.15 MPTEMP,37,3173.15 MPTEMP,38,3273.15 FACT_ENTH=1000.0 MPDATA,ENTH,mat3_top,1,100 MPDATA,ENTH,mat3_top,2,300 MPDATA,ENTH,mat3_top,3,500 MPDATA,ENTH,mat3_top,4,1000 MPDATA,ENTH,mat3_top,5,9.90639*FACT_ENTH MPDATA,ENTH,mat3_top,6,72.74899*FACT_ENTH MPDATA,ENTH,mat3_top,7,137.74765*FACT_ENTH MPDATA,ENTH,mat3_top,8,215.45355*FACT_ENTH MPDATA,ENTH,mat3_top,9,217.21913*FACT_ENTH MPDATA,ENTH,mat3_top,10,253.35072*FACT_ENTH MPDATA,ENTH,mat3_top,11,255.04297*FACT_ENTH MPDATA,ENTH,mat3_top,12,268.8865*FACT_ENTH MPDATA,ENTH,mat3_top,13,311.14525*FACT_ENTH MPDATA,ENTH,mat3_top,14,386.64614*FACT_ENTH MPDATA,ENTH,mat3_top,15,678.20587*FACT_ENTH MPDATA,ENTH,mat3_top,16,755.79953*FACT_ENTH MPDATA,ENTH,mat3_top,17,814.53698*FACT_ENTH MPDATA,ENTH,mat3_top,18,836.66302*FACT_ENTH MPDATA,ENTH,mat3_top,19,1000.295*FACT_ENTH MPDATA,ENTH,mat3_top,20,1108.57351*FACT_ENTH MPDATA,ENTH,mat3_top,21,1159.07857*FACT_ENTH MPDATA,ENTH,mat3_top,22,1171.72341*FACT_ENTH MPDATA,ENTH,mat3_top,23,1258.08278*FACT_ENTH MPDATA,ENTH,mat3_top,24,1347.6612*FACT_ENTH MPDATA,ENTH,mat3_top,25,1438.32037*FACT_ENTH MPDATA,ENTH,mat3_top,26,1529.715*FACT_ENTH MPDATA,ENTH,mat3_top,27,1621.68255*FACT_ENTH MPDATA,ENTH,mat3_top,28,1713.93227*FACT_ENTH MPDATA,ENTH,mat3_top,29,1806.26271*FACT_ENTH MPDATA,ENTH,mat3_top,30,1898.66532*FACT_ENTH MPDATA,ENTH,mat3_top,31,1991.13773*FACT_ENTH MPDATA,ENTH,mat3_top,32,2083.67828*FACT_ENTH MPDATA,ENTH,mat3_top,33,2176.28536*FACT_ENTH MPDATA,ENTH,mat3_top,34,2268.95727*FACT_ENTH MPDATA,ENTH,mat3_top,35,2361.69078*FACT_ENTH MPDATA,ENTH,mat3_top,36,2454.48176*FACT_ENTH MPDATA,ENTH,mat3_top,37,2547.33011*FACT_ENTH MPDATA,ENTH,mat3_top,38,2640.23601*FACT_ENTH mptemp *SET,FACT_C *SET,FACT_DENS *SET,FACT_ENTH *SET,FACT_KXX !mat_h13=2 MPTEMP MPTEMP,1,298.15 MPTEMP,2,373.15 MPTEMP,3,473.15 MPTEMP,4,573.15 MPTEMP,5,673.15 MPTEMP,6,773.15 MPTEMP,7,873.15 MPTEMP,8,875.01567 MPTEMP,9,973.15 MPTEMP,10,1073.15 MPTEMP,11,1092.32776 MPTEMP,12,1117.50602 MPTEMP,13,1155.61939 MPTEMP,14,1173.15 MPTEMP,15,1273.15 MPTEMP,16,1368.65072 MPTEMP,17,1373.15 MPTEMP,18,1428.82356 MPTEMP,19,1437.90645 MPTEMP,20,1473.15 MPTEMP,21,1573.15 MPTEMP,22,1655.59955 MPTEMP,23,1655.80609 MPTEMP,24,1673.15 MPTEMP,25,1712.75856 MPTEMP,26,1723.71087 MPTEMP,27,1761.69286 MPTEMP,28,1773.15 MPTEMP,29,1873.15 MPTEMP,30,1973.15 MPTEMP,31,2073.15 MPTEMP,32,2173.15 MPTEMP,33,2273.15 MPTEMP,34,2373.15 MPTEMP,35,2473.15 MPTEMP,36,2573.15 MPTEMP,37,2673.15 MPTEMP,38,2773.15 MPTEMP,39,2873.15 MPTEMP,40,2973.15 MPTEMP,41,3073.15 MPTEMP,42,3173.15 MPTEMP,43,3273.15 FACT_DENS=1000.0 MPDATA,DENS,mat_h13,1,7.805597106*FACT_DENS MPDATA,DENS,mat_h13,2,7.784283775*FACT_DENS MPDATA,DENS,mat_h13,3,7.754637969*FACT_DENS MPDATA,DENS,mat_h13,4,7.723617555*FACT_DENS MPDATA,DENS,mat_h13,5,7.691255191*FACT_DENS MPDATA,DENS,mat_h13,6,7.657584709*FACT_DENS MPDATA,DENS,mat_h13,7,7.622641022*FACT_DENS MPDATA,DENS,mat_h13,8,7.621953051*FACT_DENS MPDATA,DENS,mat_h13,9,7.586668528*FACT_DENS MPDATA,DENS,mat_h13,10,7.549354465*FACT_DENS MPDATA,DENS,mat_h13,11,7.542039185*FACT_DENS MPDATA,DENS,mat_h13,12,7.618098094*FACT_DENS MPDATA,DENS,mat_h13,13,7.595766337*FACT_DENS MPDATA,DENS,mat_h13,14,7.586050348*FACT_DENS MPDATA,DENS,mat_h13,15,7.530676272*FACT_DENS MPDATA,DENS,mat_h13,16,7.477851745*FACT_DENS MPDATA,DENS,mat_h13,17,7.475536509*FACT_DENS MPDATA,DENS,mat_h13,18,7.446875484*FACT_DENS MPDATA,DENS,mat_h13,19,7.442316523*FACT_DENS MPDATA,DENS,mat_h13,20,7.424166336*FACT_DENS MPDATA,DENS,mat_h13,21,7.372437013*FACT_DENS MPDATA,DENS,mat_h13,22,7.329281374*FACT_DENS MPDATA,DENS,mat_h13,23,7.328732933*FACT_DENS MPDATA,DENS,mat_h13,24,7.31737591*FACT_DENS MPDATA,DENS,mat_h13,25,7.255845915*FACT_DENS MPDATA,DENS,mat_h13,26,7.178741338*FACT_DENS MPDATA,DENS,mat_h13,27,7.037752863*FACT_DENS MPDATA,DENS,mat_h13,28,7.028857964*FACT_DENS MPDATA,DENS,mat_h13,29,6.949761253*FACT_DENS MPDATA,DENS,mat_h13,30,6.868208899*FACT_DENS MPDATA,DENS,mat_h13,31,6.784433655*FACT_DENS MPDATA,DENS,mat_h13,32,6.69866536*FACT_DENS MPDATA,DENS,mat_h13,33,6.611129766*FACT_DENS MPDATA,DENS,mat_h13,34,6.522047477*FACT_DENS MPDATA,DENS,mat_h13,35,6.431632985*FACT_DENS MPDATA,DENS,mat_h13,36,6.340093896*FACT_DENS MPDATA,DENS,mat_h13,37,6.247630215*FACT_DENS MPDATA,DENS,mat_h13,38,6.154433763*FACT_DENS MPDATA,DENS,mat_h13,39,6.060687713*FACT_DENS MPDATA,DENS,mat_h13,40,5.966566279*FACT_DENS MPDATA,DENS,mat_h13,41,5.872234406*FACT_DENS MPDATA,DENS,mat_h13,42,5.777847649*FACT_DENS MPDATA,DENS,mat_h13,43,5.683552083*FACT_DENS MPTEMP MPTEMP,1,298.15 MPTEMP,2,373.15 MPTEMP,3,473.15 MPTEMP,4,573.15 MPTEMP,5,673.15 MPTEMP,6,773.15 MPTEMP,7,873.15 MPTEMP,8,973.15 MPTEMP,9,1073.15 MPTEMP,10,1173.15 MPTEMP,11,1273.15 MPTEMP,12,1373.15 MPTEMP,13,1473.15 MPTEMP,14,1573.15 MPTEMP,15,1673.15 MPTEMP,16,1773.15 MPTEMP,17,1873.15 MPTEMP,18,1973.15 MPTEMP,19,2073.15 MPTEMP,20,2173.15 MPTEMP,21,2273.15 MPTEMP,22,2373.15 MPTEMP,23,2473.15 MPTEMP,24,2573.15 MPTEMP,25,2673.15 MPTEMP,26,2773.15 MPTEMP,27,2873.15 MPTEMP,28,2973.15 MPTEMP,29,3073.15 MPTEMP,30,3173.15 MPTEMP,31,3273.15 FACT_C=1000.0 MPDATA,C,mat_h13,1,0.4446*FACT_C MPDATA,C,mat_h13,2,0.47537*FACT_C MPDATA,C,mat_h13,3,0.51395*FACT_C MPDATA,C,mat_h13,4,0.55669*FACT_C MPDATA,C,mat_h13,5,0.60819*FACT_C MPDATA,C,mat_h13,6,0.67373*FACT_C MPDATA,C,mat_h13,7,0.81549*FACT_C MPDATA,C,mat_h13,8,0.95838*FACT_C MPDATA,C,mat_h13,9,0.86826*FACT_C MPDATA,C,mat_h13,10,0.67803*FACT_C MPDATA,C,mat_h13,11,0.72558*FACT_C MPDATA,C,mat_h13,12,0.64371*FACT_C MPDATA,C,mat_h13,13,0.66086*FACT_C MPDATA,C,mat_h13,14,0.67942*FACT_C MPDATA,C,mat_h13,15,0.93276*FACT_C MPDATA,C,mat_h13,16,0.80655*FACT_C MPDATA,C,mat_h13,17,0.82573*FACT_C MPDATA,C,mat_h13,18,0.82777*FACT_C MPDATA,C,mat_h13,19,0.82957*FACT_C MPDATA,C,mat_h13,20,0.83103*FACT_C MPDATA,C,mat_h13,21,0.83132*FACT_C MPDATA,C,mat_h13,22,0.83153*FACT_C MPDATA,C,mat_h13,23,0.83173*FACT_C MPDATA,C,mat_h13,24,0.83192*FACT_C MPDATA,C,mat_h13,25,0.8321*FACT_C MPDATA,C,mat_h13,26,0.83226*FACT_C MPDATA,C,mat_h13,27,0.83241*FACT_C MPDATA,C,mat_h13,28,0.83248*FACT_C MPDATA,C,mat_h13,29,0.83253*FACT_C MPDATA,C,mat_h13,30,0.83257*FACT_C MPDATA,C,mat_h13,31,0.83262*FACT_C MPTEMP MPTEMP,1,298.15 MPTEMP,2,373.15 MPTEMP,3,473.15 MPTEMP,4,573.15 MPTEMP,5,673.15 MPTEMP,6,773.15 MPTEMP,7,873.15 MPTEMP,8,875.01567 MPTEMP,9,973.15 MPTEMP,10,1073.15 MPTEMP,11,1092.32776 MPTEMP,12,1117.50602 MPTEMP,13,1155.61939 MPTEMP,14,1173.15 MPTEMP,15,1273.15 MPTEMP,16,1368.65072 MPTEMP,17,1373.15 MPTEMP,18,1428.82356 MPTEMP,19,1437.90645 MPTEMP,20,1473.15 MPTEMP,21,1573.15 MPTEMP,22,1655.59955 MPTEMP,23,1655.80609 MPTEMP,24,1673.15 MPTEMP,25,1712.75856 MPTEMP,26,1723.71087 MPTEMP,27,1761.69286 MPTEMP,28,1773.15 MPTEMP,29,1873.15 MPTEMP,30,1973.15 MPTEMP,31,2073.15 MPTEMP,32,2173.15 MPTEMP,33,2273.15 MPTEMP,34,2373.15 MPTEMP,35,2473.15 MPTEMP,36,2573.15 MPTEMP,37,2673.15 MPTEMP,38,2773.15 MPTEMP,39,2873.15 MPTEMP,40,2973.15 MPTEMP,41,3073.15 MPTEMP,42,3173.15 MPTEMP,43,3273.15 FACT_ENTH=1000.0 MPDATA,ENTH,mat_h13,1,1.22414*FACT_ENTH MPDATA,ENTH,mat_h13,2,19.30644*FACT_ENTH MPDATA,ENTH,mat_h13,3,68.76324*FACT_ENTH MPDATA,ENTH,mat_h13,4,122.24202*FACT_ENTH MPDATA,ENTH,mat_h13,5,180.39682*FACT_ENTH MPDATA,ENTH,mat_h13,6,244.34854*FACT_ENTH MPDATA,ENTH,mat_h13,7,315.86824*FACT_ENTH MPDATA,ENTH,mat_h13,8,317.39159*FACT_ENTH MPDATA,ENTH,mat_h13,9,403.58612*FACT_ENTH MPDATA,ENTH,mat_h13,10,500.98604*FACT_ENTH MPDATA,ENTH,mat_h13,11,517.40312*FACT_ENTH MPDATA,ENTH,mat_h13,12,571.92769*FACT_ENTH MPDATA,ENTH,mat_h13,13,600.77262*FACT_ENTH MPDATA,ENTH,mat_h13,14,612.59991*FACT_ENTH MPDATA,ENTH,mat_h13,15,682.70577*FACT_ENTH MPDATA,ENTH,mat_h13,16,754.59285*FACT_ENTH MPDATA,ENTH,mat_h13,17,757.48782*FACT_ENTH MPDATA,ENTH,mat_h13,18,793.60412*FACT_ENTH MPDATA,ENTH,mat_h13,19,799.64554*FACT_ENTH MPDATA,ENTH,mat_h13,20,822.84589*FACT_ENTH MPDATA,ENTH,mat_h13,21,889.90497*FACT_ENTH MPDATA,ENTH,mat_h13,22,946.62651*FACT_ENTH MPDATA,ENTH,mat_h13,23,947.00369*FACT_ENTH MPDATA,ENTH,mat_h13,24,961.23459*FACT_ENTH MPDATA,ENTH,mat_h13,25,1026.53923*FACT_ENTH MPDATA,ENTH,mat_h13,26,1098.55322*FACT_ENTH MPDATA,ENTH,mat_h13,27,1247.8653*FACT_ENTH MPDATA,ENTH,mat_h13,28,1257.07667*FACT_ENTH MPDATA,ENTH,mat_h13,29,1339.20152*FACT_ENTH MPDATA,ENTH,mat_h13,30,1421.87806*FACT_ENTH MPDATA,ENTH,mat_h13,31,1504.74751*FACT_ENTH MPDATA,ENTH,mat_h13,32,1587.78077*FACT_ENTH MPDATA,ENTH,mat_h13,33,1670.90225*FACT_ENTH MPDATA,ENTH,mat_h13,34,1754.0446*FACT_ENTH MPDATA,ENTH,mat_h13,35,1837.20707*FACT_ENTH MPDATA,ENTH,mat_h13,36,1920.38899*FACT_ENTH MPDATA,ENTH,mat_h13,37,2003.58942*FACT_ENTH MPDATA,ENTH,mat_h13,38,2086.80705*FACT_ENTH MPDATA,ENTH,mat_h13,39,2170.04013*FACT_ENTH MPDATA,ENTH,mat_h13,40,2253.28445*FACT_ENTH MPDATA,ENTH,mat_h13,41,2336.53424*FACT_ENTH MPDATA,ENTH,mat_h13,42,2419.78918*FACT_ENTH MPDATA,ENTH,mat_h13,43,2503.04912*FACT_ENTH MPTEMP MPTEMP,1,298.15 MPTEMP,2,373.15 MPTEMP,3,473.15 MPTEMP,4,573.15 MPTEMP,5,673.15 MPTEMP,6,773.15 MPTEMP,7,873.15 MPTEMP,8,875.01567 MPTEMP,9,973.15 MPTEMP,10,1073.15 MPTEMP,11,1092.32776 MPTEMP,12,1117.50602 MPTEMP,13,1155.61939 MPTEMP,14,1173.15 MPTEMP,15,1273.15 MPTEMP,16,1368.65072 MPTEMP,17,1373.15 MPTEMP,18,1428.82356 MPTEMP,19,1437.90645 MPTEMP,20,1473.15 MPTEMP,21,1573.15 MPTEMP,22,1655.59955 MPTEMP,23,1655.80609 MPTEMP,24,1673.15 MPTEMP,25,1712.75856 MPTEMP,26,1723.71087 MPTEMP,27,1761.69286 MPTEMP,28,1773.15 MPTEMP,29,1873.15 MPTEMP,30,1973.15 MPTEMP,31,2073.15 MPTEMP,32,2173.15 MPTEMP,33,2273.15 MPTEMP,34,2373.15 MPTEMP,35,2473.15 MPTEMP,36,2573.15 MPTEMP,37,2673.15 MPTEMP,38,2773.15 MPTEMP,39,2873.15 MPTEMP,40,2973.15 MPTEMP,41,3073.15 MPTEMP,42,3173.15 MPTEMP,43,3273.15 FACT_KXX=1.0 MPDATA,KXX,mat_h13,1,26.81964092*FACT_KXX MPDATA,KXX,mat_h13,2,29.09009252*FACT_KXX MPDATA,KXX,mat_h13,3,31.50733865*FACT_KXX MPDATA,KXX,mat_h13,4,32.93937262*FACT_KXX MPDATA,KXX,mat_h13,5,33.28016053*FACT_KXX MPDATA,KXX,mat_h13,6,32.69629538*FACT_KXX MPDATA,KXX,mat_h13,7,31.57973044*FACT_KXX MPDATA,KXX,mat_h13,8,31.5404691*FACT_KXX MPDATA,KXX,mat_h13,9,29.84011757*FACT_KXX MPDATA,KXX,mat_h13,10,28.854242*FACT_KXX MPDATA,KXX,mat_h13,11,28.76518791*FACT_KXX MPDATA,KXX,mat_h13,12,27.07334962*FACT_KXX MPDATA,KXX,mat_h13,13,27.49975399*FACT_KXX MPDATA,KXX,mat_h13,14,27.6958099*FACT_KXX MPDATA,KXX,mat_h13,15,28.83582751*FACT_KXX MPDATA,KXX,mat_h13,16,29.95437684*FACT_KXX MPDATA,KXX,mat_h13,17,30.00814562*FACT_KXX MPDATA,KXX,mat_h13,18,30.67351431*FACT_KXX MPDATA,KXX,mat_h13,19,30.78679377*FACT_KXX MPDATA,KXX,mat_h13,20,31.20782738*FACT_KXX MPDATA,KXX,mat_h13,21,32.40185987*FACT_KXX MPDATA,KXX,mat_h13,22,33.38530275*FACT_KXX MPDATA,KXX,mat_h13,23,33.38649655*FACT_KXX MPDATA,KXX,mat_h13,24,33.58165321*FACT_KXX MPDATA,KXX,mat_h13,25,33.80010711*FACT_KXX MPDATA,KXX,mat_h13,26,33.74368979*FACT_KXX MPDATA,KXX,mat_h13,27,32.97180222*FACT_KXX MPDATA,KXX,mat_h13,28,33.17149632*FACT_KXX MPDATA,KXX,mat_h13,29,34.91446362*FACT_KXX MPDATA,KXX,mat_h13,30,36.65743091*FACT_KXX MPDATA,KXX,mat_h13,31,38.4003982*FACT_KXX MPDATA,KXX,mat_h13,32,40.1433655*FACT_KXX MPDATA,KXX,mat_h13,33,41.8863328*FACT_KXX MPDATA,KXX,mat_h13,34,43.62930006*FACT_KXX MPDATA,KXX,mat_h13,35,45.37226737*FACT_KXX MPDATA,KXX,mat_h13,36,47.11523464*FACT_KXX MPDATA,KXX,mat_h13,37,48.85820192*FACT_KXX MPDATA,KXX,mat_h13,38,50.60116912*FACT_KXX MPDATA,KXX,mat_h13,39,52.34413634*FACT_KXX MPDATA,KXX,mat_h13,40,54.0871035*FACT_KXX MPDATA,KXX,mat_h13,41,55.8300706*FACT_KXX MPDATA,KXX,mat_h13,42,57.57303761*FACT_KXX MPDATA,KXX,mat_h13,43,59.31600449*FACT_KXX MPTEMP *SET,FACT_C *SET,FACT_DENS *SET,FACT_ENTH *SET,FACT_KXX |
|
怎么個不行,請你具體說明一下,有什么警告和/或錯誤信息。 簡單看了一下你的命令流,發(fā)現(xiàn)幾個問題如下: 在循環(huán)中,/post1...finishi 之后,持續(xù)近日求解模塊,執(zhí)行了重啟動,但是接下來的 solve 被你注解掉了,即沒有再求解。然后回到循環(huán)開始處有一次 solve,然后有有 /solu 。這樣的安排不知道軟件是如何執(zhí)行的。 總值,你的命令流中有多處 /solu 命令在 solve 之后的情況,而按照 ANSYS 的規(guī)定,每一次進入 /solu 都會被認為是開始一個新的分析,將從頭開始進行分析,對你的情況就有可能出現(xiàn)一線不到的問題。所以你需要整理一下你的 /solu 命令,不要隨處執(zhí)行這一命令,要去掉多余的 /solu 命令。 其次,你的材料屬性定義放在命令流的最后,這也是不對的。你需要把它們移到合適的位置。 |
銅蟲 (小有名氣)
|
我是按照您的建議修改后發(fā)現(xiàn)還是出現(xiàn)了重啟動后所有單元都激活了的情況,至于注釋掉的命令可以不管,因為我在知乎上看到要用solve初始化,但失敗了。其次,我用/input的命令導(dǎo)入材料物性,為了方便還是把物性寫出來了。 發(fā)自小木蟲Android客戶端 |
銅蟲 (小有名氣)
送紅花一朵 |
謝謝您的耐心解答。 我在重啟動前將死單元放到不同的 component 中并輸出其中的死單元編號到硬盤上,重啟動命令之后再執(zhí)行pararest 命令,死單元的component消失了而死單元的編號數(shù)組還在且不變。 我也懷疑重啟動是否能用于生死單元和*do循環(huán)中了,ANSYS Help 19.2 5.8.1.1.1 Multiframe Restart中提到的limitation中有“The .Rnnn file does not save EKILL and EALIVE commands. If EKILL or EALIVE are required in the restarted session, issue them again. Restarting with tabular input (EALIVE,%table%) is not supported.”和“Multiframe restart does not support the arc-length method (ARCLEN), reading and solving multiple load steps (LSSOLVE), or nested *DO loops.”后面一個限制應(yīng)該是不能在restart后嵌套*do循環(huán)吧? Help的5.8.1.2.也提到用solve來初始化,不過我這是單點重啟動,和多點重啟動應(yīng)該不同吧? 我/input的文件內(nèi)容就是我的命令流里*ENDDO后面的物性參數(shù),分別以“!mat_top=1”和“!mat_h13=2”開頭 |
銅蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
|
15779376950 2026-03-01 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 281求調(diào)劑 +3 | Y?l?h 2026-03-04 | 3/150 |
|
|
[考研] 331求調(diào)劑 +3 | zzZ&zZ 2026-03-03 | 3/150 |
|
|
[考研] 264求調(diào)劑 +7 | 26調(diào)劑 2026-03-03 | 7/350 |
|
|
[考研] 291求調(diào)劑 +4 | Afy123456 2026-03-03 | 7/350 |
|
|
[考研] 298求調(diào)劑 +3 | 人間唯你是清歡 2026-03-03 | 4/200 |
|
|
[考研] 一志愿中科大能動297求調(diào)劑,本科川大 +3 | 邵11 2026-03-03 | 3/150 |
|
|
[考研] 284求調(diào)劑 +6 | 天下熯 2026-03-02 | 6/300 |
|
|
[考研] 299求調(diào)劑 +5 | kkcoco25 2026-03-02 | 9/450 |
|
|
[考研] 調(diào)劑 +5 | 13853210211 2026-03-02 | 7/350 |
|
|
[考研] 0854復(fù)試調(diào)劑 276 +5 | wmm9 2026-03-01 | 7/350 |
|
|
[考研] 化學(xué),材料,環(huán)境類求調(diào)劑 +7 | 考研版棒棒 2026-03-02 | 7/350 |
|
|
[考研] 材料工程274求調(diào)劑 +5 | Lilithan 2026-03-01 | 5/250 |
|
|
[考研] 275求調(diào)劑 +7 | 明遠求學(xué) 2026-03-01 | 7/350 |
|
|
[考研] 一志愿東北大學(xué)材料專碩328,求調(diào)劑 +3 | shs1083 2026-03-02 | 3/150 |
|
|
[考研] 一志愿鄭大材料學(xué)碩298分,求調(diào)劑 +6 | wsl111 2026-03-01 | 6/300 |
|
|
[考研] 272求調(diào)劑 +6 | 田智友 2026-02-28 | 6/300 |
|
|
[考研] 313求調(diào)劑 +3 | 水流年lc 2026-02-28 | 3/150 |
|
|
[論文投稿]
求助coordination chemistry reviews 的寫作模板
10+3
|
ljplijiapeng 2026-02-27 | 4/200 |
|
|
[考研] 304求調(diào)劑 +3 | 52hz~~ 2026-02-28 | 5/250 |
|