我模擬的溫度應力不對，有人可以看看嗎。

可以微信紅包100，后續還有問題

!-test1

!有重力加速度

!對稱 X=0

!長度10/2

!時間7天

!提取溫度差，加載到結構模型，求解溫度應力，參考溫度為0

!48小時拆模

!澆筑溫度為20

!有重力

/clear

*set, day,7           !計算天數，先調試1天

*set, hour, day*24       !計算小時數

*set, BEG_TEMP, 24       !初始溫度

*set, air_max, 28        !最高氣溫 

*set, air_min, 22        !最低氣溫

*set, l,4             !澆注墻體長度的1/2

*set, h, 3.8             !澆注墻體的高度

*set, t, 0.5             !澆注墻體的厚度 

*set, plank, 0.8         !底板厚

*set, board, 2           !底板寬計算長度 

*set, initial, 20        !底板澆筑溫度

*set, houremove, 480     !拆除模版時間 

!!!!!!!!!!!!!!!!!!!!!

/prep7

et,1,solid70

mp,dens,1,2400          !素混凝土容重

mp,kxx,1,8.51           !混凝土的導熱系數

mp,c,1,0.944            !混凝土的比熱容

!墻體、底板、模板建模

block,,l,0.4,h,,-t             !墻體X長L-0，Y高H-0.4，Z厚-T

block,,l,0.4,0,,-t             !高0.4墻體X長L-0，Y高-0.4， Z厚-T

block,,l,-plank,0,0,-t         !墻體下的底板X長L-0，Y高PLANK，Z厚-T

block,,l,-plank,0,-t,-board    !X長L-0，Y高PLANK，Z厚T-BOARD

!把已建立的體連接起來，同時創建模板

vglue, 1,2

vglue, 3,5

vglue, 3,4

!墻體及底板分網,寬度方向分5份;長度方向分成35*3/2份;高度分16份

lsel,s,length,,t,t           !選擇長度為t的線

lesize,all,,,4               !劃分為4等分 

lsel,s,length,,l,l           !選擇長度為l的線

lesize,all,0.3,             !單元長度0.3m 

lsel,s,length,,h-0.4,h-0.4   !選擇長度為h-0.4的線

lesize,all,,,16              !劃分為16等分 

lsel,s,length,,0.4,0.4       !選擇長度為0.4的線

lesize,all,,,2               !劃分為2等分 

lsel,s,length,,plank,plank   !選擇長度為plank的線

lesize,all,,,2               !劃分為2等分

lsel,s,length,,board-t,board-t   !選擇長度為board-t的線

lesize,all,,,5               !劃分為5等分

vmesh,1

vmesh,2

vmesh,3

vmesh,5

!選中混凝土表面,建立conc1,conc2組合面

allsel

asel,s,,,1

asel,a,,,2

asel,a,,,4

asel,a,,,6

asel,a,,,8

asel,a,,,11

asel,a,,,7

asel,a,,,14

asel,a,,,18

asel,a,,,25

asel,a,,,12

asel,a,,,19

cm,conc1,area    !選中的面組成1組,編號為conc1

!建立對稱面約束

allsel

asel,s,loc,x,0   !建立X坐標為0的面

da,all,symm

DTRAN !將實體模型中的約束轉換到有限元模型中

*get,n_num,node,,count !得到節點總數

*dim,he,array,hour,n_num !存放節點溫度

*dim,t_temp,array,hour,n_num !存放節點溫差

*dim,hgen,array,hour,1 !存放單位時間產熱

!計算單位時間產熱hgen(i)

*do, i, 1, hour

HE00=330*536*(1-exp(-0.69*(i/24)**0.56))

HE01=330*536*(1-exp(-0.69*(((i+1)/24))**0.56))

*set,hgen(i),HE01-HE00

*enddo

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!求解

/SOLU

ANTYPE,TRANSIENT,NEW   !定義分析類型

nropt, full    

allsel, all      !

nsel,s,,,all     !選所有節點

ic,all,temp,initial   !初始溫度20

*do, i, 1, hour,1

time,i

timint, on

vsel,s,,,1

eslv,s

BFe,ALL,HGEN,,hgen(i)         !加水化熱

cmsel,s,conc1

!在混凝土表面施加對流系數和空氣溫度

*set,ta,20

!sfa,all,,CONV,53,ta

sfa,all,,CONV,53,air_max/2+air_min/2+(air_max/2-air_min/2)*cos(15*(i-7))

allsel

TRNOPT,FULL

deltim,1,,,1  !

AUTOTS,ON

KBC,0

outres,all,1

solve

!得到節點溫差

*do,ii,1,n_num

*get,te,node,ii,temp

*SET,he(i,ii),te  

*if,i,eq,1,then

t_temp(i,ii)=he(i,ii)-initial

*else

t_temp(i,ii)=he(i,ii)-he(i-1,ii)

*endif

*enddo

allsel

*enddo

save

/PREP7 

allsel,all  

sfdele,all,all   !刪除節點上的面荷載   

bfdele,all,all   !刪除面上的體載荷  

SFEDELE,ALL,ALL,CONV !刪對流邊界條件 

SFDELE,ALL,CONV 

ETCHG,TTS

mp,dens,999,2400  

mp,prxy,999,0.167 

mp,alpx,999,0.8e-5

mp,ex,999,0.0022696e10  

esel,s,mat,,1   

mpchg,999,all 

 

*SET,Ec,3.25e10 

*dim,exx,array,hour

*do,i,1,hour,1   

exx(i)=(1-exp(-0.28*(i/24)**0.52))*Ec 

*enddo  

*do,i,1,hour,1   

mp,dens,i,2400  

mp,prxy,i,0.167 

mp,alpx,i,0.8e-5

mp,ex,i,exx(i)  

*enddo  

esel,s,mat,,999   

mpchg,1,all 

allsel  

!施加約束

asel,s,loc,x,0

da,all,symm

!da,all,ux,0

nsel,s,loc,x,l

d,all,ux,0

nsel,s,loc,y,-plank

d,all,all

nsel,s,loc,z,0

d,all,uz,0

nsel,s,loc,y,0,10

nsel,r,loc,z,-t

d,all,uz,0

nsel,s,loc,z,-board

d,all,uz,0

FINISH  

/SOL

antype,static 

ALLSEL,ALL 

ACEL,0,9.8,0

allsel,all

outres,all,all  

tref,20

  

*do,HOU,1,hour

time,HOU

*do,j,1,n_num

bf,j,temp,t_temp(hou,j)

!LDREAD,TEMP,,,HOU,,test8,rth,

*enddo

allsel  

vsel,s,,,1

eslv,s,1

!隨著時間變化，材料屬性變化

mpchg,HOU,all

allsel

solve  

*enddo  

finish