;new ^Cj3\G4,
8v=t -GJW
;restore therm_model.sav m @)Ya*=<
Jyp 7+M]
;model mohr ;本构模型 = GiN~$d
;model th_iso p[;@9!t
Z2-"N B
def get_znum ;获取单元总数 8~O0P=
aY DM)b}
znt=zone_head j %':M
loop while znt # null =4OV
}z= I
znum=znum+1 x1" 8K
znt=z_next(znt) \FnR'ne
end_loop N(O* "1b
oxJAI4{y
4
end r[#*..Y
get_znum J<&?Hb*|
?KE:KV[Y
tJm1Q#||
def weibull_dis ;weibull分布 @ 0/EKWF
):n'B` f}z
znt=zone_head &z>iqm"Ww
loop while znt # null Dv4 H^
weibull_young=((-ln(urand))^(1.0/emod_weib_m))*ex_young ;弹性模量weibull分布 eQMa9_
weibull_thexp=((-ln(urand))^(1.0/exth_weib_m))*ex_thexp ;热膨胀系数weibull分布 zhY] !
R~;8v 1>K
;----------某个量不变化时候时候-------------- f=Oj01Ut*
weibull_young=18.0000e9 ;弹性模量不变化 7&(h_} Z
;weibull_thexp=4.6200e-5 ;膨胀系数不变化 NS "1zR+
;------------------------------------- tq L2' (=
<S12=<c?'
bulk_mod=weibull_young/3.0/(1.0-2*p_rat) `trcYmR=k
shear_mod=weibull_young/2.0/(1.0+p_rat) A-h[vP!v|
6 LqF*$+$`
z_prop(znt,'bulk')=bulk_mod ;给每个单元赋参数 .}E@ 7^X
z_prop(znt,'shear')=shear_mod o)'06FF\$
:W+%jn
z_prop(znt,'thexp')=weibull_thexp D4?cnwU
>D_)z/v?"
z_gnum=z_gnum+1 JM53sx4&
$2a_!/
table(1,z_gnum)=weibull_young V@\u<LO0G
table(2,z_gnum)=weibull_thexp 6zGeGW
c<{~j~+
znt=z_next(znt) UlPGB2B
end_loop cs[nFfM
end 3PkU>+.6
v|@EuN14<
;-------------------均匀系数、尺度参数、变异系数------------------ F'@ 9kdp
;------------------弹性模量-------------------------- /4upw`35]
;1.5,19.94018,0.6790---2.0,20.31144,0.5228---4.0,19.85878,0.2805---8.0,19.11437,0.1882 0%yPuY>
;15, 18.63933,0.0810---20 18.489980.0625---40, 18.25188,0.0306---100,18.10319,0.0130 sWi4+PAM0
w BoP&l
set emod_weib_m=1.5 ;弹摸均匀系数 Sae*VvT6
set ex_young=19.94018e9 ;尺度参数 Cn_r?1{W
set p_rat=0.3 ;泊松比 N,*'")k9
M}
+s_h9
;--------------------膨胀系数----------------------- :w7?]y6~S
;1.5,5.11798,0.6790---2.0,5.21327,0.5228---4.0,5.09709,0.2805---8.0,4.90602,0.1882 &egP3
;15, 4.78409,0.0810---20 4.74576,0.0625---40, 4.68465,0.0306---100,4.64648,0.0130 J?1Eh14KZ
<X ?xr f
set exth_weib_m=1.5 ;膨胀系数均匀系数 *|g l1S
set ex_thexp=5.11798e-5 ;尺度参数 rmdg~
<<-BQ
l~
(^]3l%Ed
weibull_dis
