/BATCH  
!/input,menust,tmp,''
! /GRA,POWER  
! /GST,ON 
! /PLO,INFO,3 
! /GRO,CURL,ON
! /CPLANE,1   
! /REPLOT,RESIZE  
WPSTYLE,,,,,,,,0

!所有的參數  
*set,UW,2317.25				!單位重
*set,CF,1089.7515				!新拌比熱	(硬固比熱*1.15)
*set,CH,947.61				!硬固比熱
*set,KUC,132804.964				!硬固熱傳導係數
*set,KUC_F,167998.2795			!新拌熱傳導係數=新拌熱擴散係數*新拌比熱*單位重*時間,	(新拌熱擴散係數為硬固熱擴散係數*1.1)
*set,Ea,39149.10464				!表徵活化能
*set,Alphau,0.8			!極限水化度
*set,Tau,2.917				!水化度參數tau
*set,Beta,0.9				!水化度參數beta
*set,HT,167551992.0				!總單位體積水化熱
*set,A1_CONV,253550.0			!面一散熱係數
*set,A2_CONV,253550.0			!面二散熱係數
*set,A3_CONV,0			!面三散熱係數
*set,A4_CONV,0			!面四散熱係數
*set,A5_CONV,253550.0			!面五散熱係數
*set,A6_CONV,253550.0			!面六散熱係數

!*  
/NOPR   
KEYW,PR_SET,1   
KEYW,PR_STRUC,0 
KEYW,PR_THERM,1 
KEYW,PR_FLUID,0 
KEYW,PR_ELMAG,0 
KEYW,MAGNOD,0   
KEYW,MAGEDG,0   
KEYW,MAGHFE,0   
KEYW,MAGELC,0   
KEYW,PR_MULTI,0 
/GO 

/PREP7  
!*  
ET,1,SOLID70

BLOCK,-0.9,0.9,-0.2,0.2,-1.0,1.0,           !分析尺寸(m)
ESIZE,0.2,0,                  		!元素大小(m)
MSHAPE,0,3D 
MSHKEY,1 
!*  
CM,_Y,VOLU  
VSEL, , , ,1 
CM,_Y1,VOLU 
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2 
ANTYPE,4
!*  
TRNOPT,FULL 
LUMPM,0 
!*  
TUNIF,25,		!初始溫度
/GO 

*get,ele_num,elem,0,count

*DO,x,1,ele_num*48*7+ele_num,1    !element的數量(由整體尺寸及mesh尺寸決定)*一天的time step數量(本研究一階time step定為1/48天)*分析齡期
MPTEMP,1,,,,,,
MPDATA,DENS,x,,UW       !單位重
*ENDDO

*DO,x,1,ele_num,1
MPTEMP,1,,,,,,
MPDATA,C,x,,CF      	!設定新拌比熱
MPTEMP,1,,,,,,
MPDATA,KXX,x,,KUC_F  	!設定新拌熱擴散係數*新拌比熱(1048.0065)*單位重(2359.8)*時間(24)
MPDATA,KYY,x,,KUC_F
MPDATA,KZZ,x,,KUC_F
*ENDDO

*DO,x,1,ele_num,1      
MPCHG,x,x
*ENDDO

!OutDoorTEMP
*DIM,outdoorTEMP,,169,1  !定義氣溫矩陣
*set,   outdoorTEMP(1,      1)      ,             25
*set,   outdoorTEMP(2,      1)      ,             25
*set,   outdoorTEMP(3,      1)      ,             25
*set,   outdoorTEMP(4,      1)      ,             25
*set,   outdoorTEMP(5,      1)      ,             25
*set,   outdoorTEMP(6,      1)      ,             25
*set,   outdoorTEMP(7,      1)      ,             25
*set,   outdoorTEMP(8,      1)      ,             25
*set,   outdoorTEMP(9,      1)      ,             25
*set,   outdoorTEMP(10,     1)      ,             25
*set,   outdoorTEMP(11,     1)      ,             25
*set,   outdoorTEMP(12,     1)      ,             25
*set,   outdoorTEMP(13,     1)      ,             25
*set,   outdoorTEMP(14,     1)      ,             25
*set,   outdoorTEMP(15,     1)      ,             25
*set,   outdoorTEMP(16,     1)      ,             25
*set,   outdoorTEMP(17,     1)      ,             25
*set,   outdoorTEMP(18,     1)      ,             25
*set,   outdoorTEMP(19,     1)      ,             25
*set,   outdoorTEMP(20,     1)      ,             25
*set,   outdoorTEMP(21,     1)      ,             25
*set,   outdoorTEMP(22,     1)      ,             25
*set,   outdoorTEMP(23,     1)      ,             25
*set,   outdoorTEMP(24,     1)      ,             25
*set,   outdoorTEMP(25,     1)      ,             25
*set,   outdoorTEMP(26,     1)      ,             25
*set,   outdoorTEMP(27,     1)      ,             25
*set,   outdoorTEMP(28,     1)      ,             25
*set,   outdoorTEMP(29,     1)      ,             25
*set,   outdoorTEMP(30,     1)      ,             25
*set,   outdoorTEMP(31,     1)      ,             25
*set,   outdoorTEMP(32,     1)      ,             25
*set,   outdoorTEMP(33,     1)      ,             25
*set,   outdoorTEMP(34,     1)      ,             25
*set,   outdoorTEMP(35,     1)      ,             25
*set,   outdoorTEMP(36,     1)      ,             25
*set,   outdoorTEMP(37,     1)      ,             25
*set,   outdoorTEMP(38,     1)      ,             25
*set,   outdoorTEMP(39,     1)      ,             25
*set,   outdoorTEMP(40,     1)      ,             25
*set,   outdoorTEMP(41,     1)      ,             25
*set,   outdoorTEMP(42,     1)      ,             25
*set,   outdoorTEMP(43,     1)      ,             25
*set,   outdoorTEMP(44,     1)      ,             25
*set,   outdoorTEMP(45,     1)      ,             25
*set,   outdoorTEMP(46,     1)      ,             25
*set,   outdoorTEMP(47,     1)      ,             25
*set,   outdoorTEMP(48,     1)      ,             25
*set,   outdoorTEMP(49,     1)      ,             25
*set,   outdoorTEMP(50,     1)      ,             25
*set,   outdoorTEMP(51,     1)      ,             25
*set,   outdoorTEMP(52,     1)      ,             25
*set,   outdoorTEMP(53,     1)      ,             25
*set,   outdoorTEMP(54,     1)      ,             25
*set,   outdoorTEMP(55,     1)      ,             25
*set,   outdoorTEMP(56,     1)      ,             25
*set,   outdoorTEMP(57,     1)      ,             25
*set,   outdoorTEMP(58,     1)      ,             25
*set,   outdoorTEMP(59,     1)      ,             25
*set,   outdoorTEMP(60,     1)      ,             25
*set,   outdoorTEMP(61,     1)      ,             25
*set,   outdoorTEMP(62,     1)      ,             25
*set,   outdoorTEMP(63,     1)      ,             25
*set,   outdoorTEMP(64,     1)      ,             25
*set,   outdoorTEMP(65,     1)      ,             25
*set,   outdoorTEMP(66,     1)      ,             25
*set,   outdoorTEMP(67,     1)      ,             25
*set,   outdoorTEMP(68,     1)      ,             25
*set,   outdoorTEMP(69,     1)      ,             25
*set,   outdoorTEMP(70,     1)      ,             25
*set,   outdoorTEMP(71,     1)      ,             25
*set,   outdoorTEMP(72,     1)      ,             25
*set,   outdoorTEMP(73,     1)      ,             25
*set,   outdoorTEMP(74,     1)      ,             25
*set,   outdoorTEMP(75,     1)      ,             25
*set,   outdoorTEMP(76,     1)      ,             25
*set,   outdoorTEMP(77,     1)      ,             25
*set,   outdoorTEMP(78,     1)      ,             25
*set,   outdoorTEMP(79,     1)      ,             25
*set,   outdoorTEMP(80,     1)      ,             25
*set,   outdoorTEMP(81,     1)      ,             25
*set,   outdoorTEMP(82,     1)      ,             25
*set,   outdoorTEMP(83,     1)      ,             25
*set,   outdoorTEMP(84,     1)      ,             25
*set,   outdoorTEMP(85,     1)      ,             25
*set,   outdoorTEMP(86,     1)      ,             25
*set,   outdoorTEMP(87,     1)      ,             25
*set,   outdoorTEMP(88,     1)      ,             25
*set,   outdoorTEMP(89,     1)      ,             25
*set,   outdoorTEMP(90,     1)      ,             25
*set,   outdoorTEMP(91,     1)      ,             25
*set,   outdoorTEMP(92,     1)      ,             25
*set,   outdoorTEMP(93,     1)      ,             25
*set,   outdoorTEMP(94,     1)      ,             25
*set,   outdoorTEMP(95,     1)      ,             25
*set,   outdoorTEMP(96,     1)      ,             25
*set,   outdoorTEMP(97,     1)      ,             25
*set,   outdoorTEMP(98,     1)      ,             25
*set,   outdoorTEMP(99,     1)      ,             25
*set,   outdoorTEMP(100,    1)      ,             25
*set,   outdoorTEMP(101,    1)      ,             25
*set,   outdoorTEMP(102,    1)      ,             25
*set,   outdoorTEMP(103,    1)      ,             25
*set,   outdoorTEMP(104,    1)      ,             25
*set,   outdoorTEMP(105,    1)      ,             25
*set,   outdoorTEMP(106,    1)      ,             25
*set,   outdoorTEMP(107,    1)      ,             25
*set,   outdoorTEMP(108,    1)      ,             25
*set,   outdoorTEMP(109,    1)      ,             25
*set,   outdoorTEMP(110,    1)      ,             25
*set,   outdoorTEMP(111,    1)      ,             25
*set,   outdoorTEMP(112,    1)      ,             25
*set,   outdoorTEMP(113,    1)      ,             25
*set,   outdoorTEMP(114,    1)      ,             25
*set,   outdoorTEMP(115,    1)      ,             25
*set,   outdoorTEMP(116,    1)      ,             25
*set,   outdoorTEMP(117,    1)      ,             25
*set,   outdoorTEMP(118,    1)      ,             25
*set,   outdoorTEMP(119,    1)      ,             25
*set,   outdoorTEMP(120,    1)      ,             25
*set,   outdoorTEMP(121,    1)      ,             25
*set,   outdoorTEMP(122,    1)      ,             25
*set,   outdoorTEMP(123,    1)      ,             25
*set,   outdoorTEMP(124,    1)      ,             25
*set,   outdoorTEMP(125,    1)      ,             25
*set,   outdoorTEMP(126,    1)      ,             25
*set,   outdoorTEMP(127,    1)      ,             25
*set,   outdoorTEMP(128,    1)      ,             25
*set,   outdoorTEMP(129,    1)      ,             25
*set,   outdoorTEMP(130,    1)      ,             25
*set,   outdoorTEMP(131,    1)      ,             25
*set,   outdoorTEMP(132,    1)      ,             25
*set,   outdoorTEMP(133,    1)      ,             25
*set,   outdoorTEMP(134,    1)      ,             25
*set,   outdoorTEMP(135,    1)      ,             25
*set,   outdoorTEMP(136,    1)      ,             25
*set,   outdoorTEMP(137,    1)      ,             25
*set,   outdoorTEMP(138,    1)      ,             25
*set,   outdoorTEMP(139,    1)      ,             25
*set,   outdoorTEMP(140,    1)      ,             25
*set,   outdoorTEMP(141,    1)      ,             25
*set,   outdoorTEMP(142,    1)      ,             25
*set,   outdoorTEMP(143,    1)      ,             25
*set,   outdoorTEMP(144,    1)      ,             25
*set,   outdoorTEMP(145,    1)      ,             25
*set,   outdoorTEMP(146,    1)      ,             25
*set,   outdoorTEMP(147,    1)      ,             25
*set,   outdoorTEMP(148,    1)      ,             25
*set,   outdoorTEMP(149,    1)      ,             25
*set,   outdoorTEMP(150,    1)      ,             25
*set,   outdoorTEMP(151,    1)      ,             25
*set,   outdoorTEMP(152,    1)      ,             25
*set,   outdoorTEMP(153,    1)      ,             25
*set,   outdoorTEMP(154,    1)      ,             25
*set,   outdoorTEMP(155,    1)      ,             25
*set,   outdoorTEMP(156,    1)      ,             25
*set,   outdoorTEMP(157,    1)      ,             25
*set,   outdoorTEMP(158,    1)      ,             25
*set,   outdoorTEMP(159,    1)      ,             25
*set,   outdoorTEMP(160,    1)      ,             25
*set,   outdoorTEMP(161,    1)      ,             25
*set,   outdoorTEMP(162,    1)      ,             25
*set,   outdoorTEMP(163,    1)      ,             25
*set,   outdoorTEMP(164,    1)      ,             25
*set,   outdoorTEMP(165,    1)      ,             25
*set,   outdoorTEMP(166,    1)      ,             25
*set,   outdoorTEMP(167,    1)      ,             25
*set,   outdoorTEMP(168,    1)      ,             25
*set,   outdoorTEMP(169,    1)      ,             25


SFA,1,1,CONV,A1_CONV,outdoorTEMP(1,1)		
SFA,2,1,CONV,A2_CONV,outdoorTEMP(1,1)		
SFA,3,1,CONV,A3_CONV,outdoorTEMP(1,1)		
SFA,4,1,CONV,A4_CONV,outdoorTEMP(1,1)		
SFA,5,1,CONV,A5_CONV,outdoorTEMP(1,1)		
SFA,6,1,CONV,A6_CONV,outdoorTEMP(1,1)		
 
FINISH  

/SOL
*get,nnode,node,0,count   	!求出模型節點總數
*DIM,t,,35000,nnode       	!定義節點溫度矩陣
*DIM,te,,35000,nnode      	!定義節點等效齡期矩陣
*DIM,a,,35000,nnode       	!定義水化度矩陣
*DIM,avgDOH,,35000,ele_num  !定義element水化度矩陣
nsel,all
*do,k,1,nnode,1
*set,t(1,k),26.3                                                                                                            !給予節點溫度矩陣之初始溫升
!*set,t(1,3654),25.1
!*set,t(1,3641),23.7
*set,te(1,k),1/48*exp((-Ea/8.314)*(1/(t(1,k)+273)-1/296))                                                            !給予節點等效齡期矩陣之初始等效齡期(1/48：time step長度；40312.687：表徵活化能；8.314：理想氣體常數)
*set,a(1,k),Alphau*exp(-(Tau/te(1,k))**Beta)                                                                            !給予節點水化度矩陣初始水化度(0.7716：極限水化度；1.711：τ；0.8588：β) 
BF,k,HGEN,HT*a(1,k)*Beta*(Tau**Beta)*(exp((-Ea/8.314)*((1/((t(1,k)+273)))-1/296)))/((te(1,k))**(1+Beta))  !施加每節點之熱體荷載(143280160：HT；1.8588:β+1)
*enddo

OUTRES,ALL,ALL, 
!*  
OUTRES,ALL,ALL, 
!*  
TIME,1/48 
AUTOTS,0
DELTIM,0.001, , ,1
KBC,0   
!*  
TSRES,ERASE 
!*  
TIME,1/48  
AUTOTS,0
NSUBST,10, , ,1  
KBC,0   
!*  
TSRES,ERASE 
TIMINT,1
TINTP, ,0.25250625,0.505, , , , 
!*  
! /STATUS,SOLU
SOLVE   

*do,i,1/48,7,1/48                               

LDREAD,TEMP,i/(1/48),10, ,2,'ttcL','rth',' '    !設前一步階分析之節點溫度為下一步階分析的初始溫度
! APLOT   
FLST,2,4,5,ORDE,2   
FITEM,2,3   
FITEM,2,-6  
/GO  

nsel,all
*do,j,1,nnode,1
*get,t((i*48)+1,j),node,j,temp                                                                                                                                !取每一步階分析結果所有節點之溫度並輸入溫度矩陣中
*set,te((i*48)+1,j),te((i*48),j)+1/48*exp((-Ea/8.314)*(1/(0.5*(t((i*48)+1,j)+t((i*48),j))+273)-1/296))                                                 !給予每一步階節點等效齡期矩陣之等效齡期
*set,a((i*48)+1,j),Alphau*exp(-(Tau/te((i*48)+1,j))**Beta)                                                                                                !水化度曲線
BF,j,HGEN,HT*a((i*48)+1,j)*Beta*(Tau**Beta)*exp((-Ea/8.314)*(1/(0.5*(t((i*48)+1,j)+t((i*48),j))+273)-1/296))/((te((i*48)+1,j))**(1+Beta))   !施加每節點之熱體荷載
*enddo 

!將每個element中的8個節點水化度平均，作為element的水化度(需先取得每個element對應的node編號)

*get,n1,elem,1,node,1
*get,n2,elem,1,node,1
*get,n3,elem,1,node,1
*get,n4,elem,1,node,1
*get,n5,elem,1,node,1
*get,n6,elem,1,node,1
*get,n7,elem,1,node,1
*get,n8,elem,1,node,1

*do,x,1,ele_num,1

*del,n1,
*del,n2,
*del,n3,
*del,n4,
*del,n5,
*del,n6,
*del,n7,
*del,n8,
*get,n1,elem,x,node,1
*get,n2,elem,x,node,2
*get,n3,elem,x,node,3
*get,n4,elem,x,node,4
*get,n5,elem,x,node,5
*get,n6,elem,x,node,6
*get,n7,elem,x,node,7
*get,n8,elem,x,node,8
*set,avgDOH(i*48,   x),1/8*(a((i*48)+1,   n1)+a((i*48)+1, n2)+a((i*48)+1, n3)+a((i*48)+1,  n4)+a((i*48)+1,n5)+a((i*48)+1,n6)+a((i*48)+1,n7)+a((i*48)+1,n8))

*enddo

*do,m,1,169,1
*if,i,ge,(m-1)/24,and,i,lt,m/24,then

SFA,1,1,CONV,A1_CONV,outdoorTEMP(m,1)		
SFA,2,1,CONV,A2_CONV,outdoorTEMP(m,1)		
SFA,3,1,CONV,A3_CONV,outdoorTEMP(m,1)		
SFA,4,1,CONV,A4_CONV,outdoorTEMP(m,1)		
SFA,5,1,CONV,A5_CONV,outdoorTEMP(m,1)		
SFA,6,1,CONV,A6_CONV,outdoorTEMP(m,1)	

*endif

*enddo

*DO,x,1,ele_num,1
MPCHG,x+ele_num*i*48,x                                                                                                                   !element的material properity在每階time step結束後換成另一個新的material properity
MPDATA,C,x+ele_num*i*48,,CH*(1.15-0.15*avgDOH((i*48),1)/Alphau)                                                                      !重新計算新的(material properity)比熱(911.31：完全硬固(水化)的比熱)
MPDATA,KXX,x+ele_num*i*48,,KUC*(1.265-0.28*avgDOH((i*48),1)/Alphau+0.015*avgDOH((i*48),1)*avgDOH((i*48),1)/(Alphau*Alphau))      !重新計算新的(material properity)熱傳導係數(117856.5138：完全硬固的熱傳導係數=完全硬固比熱*完全硬固熱擴散係數*單位重)
MPDATA,KYY,x+ele_num*i*48,,KUC*(1.265-0.28*avgDOH((i*48),1)/Alphau+0.015*avgDOH((i*48),1)*avgDOH((i*48),1)/(Alphau*Alphau))
MPDATA,KZZ,x+ele_num*i*48,,KUC*(1.265-0.28*avgDOH((i*48),1)/Alphau+0.015*avgDOH((i*48),1)*avgDOH((i*48),1)/(Alphau*Alphau))
*ENDDO

TIME,1/48+i 
AUTOTS,0
DELTIM,0.001, , ,1
KBC,0   
!*  
TSRES,ERASE 
!*  
TIME,1/48+i  
AUTOTS,0
NSUBST,10, , ,1  
KBC,0   
!*  
TSRES,ERASE 
TIMINT,1
TINTP, ,0.25250625,0.505, , , , 
!*  
! /STATUS,SOLU
SOLVE   
ALLSEL,ALL 

*enddo
 

/POST26 
FILE,'ttcL','rth','.'
/UI,COLL,1  
NUMVAR,200  
SOLU,191,NCMIT  
STORE,MERGE 
FILLDATA,191,,,,1,1 
REALVAR,191,191 
!*  

!設定儲存Node

NSOL,2,240,TEMP,, e,
STORE,MERGE
!*  
NSOL,3,198,TEMP,, w,
STORE,MERGE
!*  
NSOL,4,216,TEMP,, s,
STORE,MERGE
!*  
NSOL,5,222,TEMP,, n,
STORE,MERGE
!*  
NSOL,6,219,TEMP,, c,
STORE,MERGE
!*   

! Name: e   
! ID:  23   
! Function: (nsol(326 ,TEMP,)+nsol(326 ,TEMP,))/2
NSOL,199,326,TEMP, 
NSOL,200,326,TEMP, 
ADD,198,200,199
FILLDATA,200,,,,2,0 
REALVAR,200,200 
QUOT,23,198,200,,e  
!   
STORE,MERGE 

! Name: w   
! ID:  24   
! Function: (nsol(263 ,TEMP,)+nsol(263 ,TEMP,))/2
NSOL,199,263,TEMP, 
NSOL,200,263,TEMP, 
ADD,198,200,199
FILLDATA,200,,,,2,0 
REALVAR,200,200 
QUOT,24,198,200,,w  
!   
STORE,MERGE 

! Name: s   
! ID:  25   
! Function: (nsol(294 ,TEMP,)+nsol(303 ,TEMP,))/2
NSOL,199,294,TEMP, 
NSOL,200,303,TEMP, 
ADD,198,200,199
FILLDATA,200,,,,2,0 
REALVAR,200,200 
QUOT,25,198,200,,s  
!   
STORE,MERGE 

! Name: n   
! ID:  26   
! Function: (nsol(286 ,TEMP,)+nsol(295 ,TEMP,))/2
NSOL,199,286,TEMP, 
NSOL,200,295,TEMP, 
ADD,198,200,199
FILLDATA,200,,,,2,0 
REALVAR,200,200 
QUOT,26,198,200,,n
!   
STORE,MERGE 

! Name: c   
! ID:  27   
! Function:  (nsol(290,TEMP,)+nsol(299,TEMP,))/2
NSOL,200,290,TEMP,  
NSOL,199,299,TEMP,  
ADD,198,200,199
FILLDATA,200,,,,2,0 
REALVAR,200,200 
QUOT,27,198,200,,c  
!   
STORE,MERGE 

! Save time history variables to file avg_test.csv   
*DIM,_P26_EXPORT,TABLE,3370,5   
VGET,_P26_EXPORT(1,0),1 
VGET,_P26_EXPORT(1,1),23
VGET,_P26_EXPORT(1,2),24
VGET,_P26_EXPORT(1,3),25
VGET,_P26_EXPORT(1,4),26
VGET,_P26_EXPORT(1,5),27
/OUTPUT,'0.8_2.917_0.9_1.8x2.0','csv','.'
*VWRITE,'TIME','e','w','s','n','c'  
%C, %C, %C, %C, %C, %C  
*VWRITE,_P26_EXPORT(1,0),_P26_EXPORT(1,1),_P26_EXPORT(1,2),_P26_EXPORT(1,3),_P26_EXPORT(1,4),_P26_EXPORT(1,5)  
%G, %G, %G, %G, %G, %G  
/OUTPUT,TERM
!/INPUT,scratch,gui 
save 
/QUIT
/EXIT,NOSAVE