;new u !fZ>kS
e$mVA}>Ybp
;restore therm_model.sav B =DV!oUg
cvx"XxE,
;model mohr ;本构模型 CI+li H
;model th_iso 0t*PQ%
,V&E"D{u
def get_znum ;获取单元总数 H D,6
49Y:}<Yd
znt=zone_head e"Z,!Q^-L
loop while znt # null |@?='E?h
znum=znum+1 <Zh\6*3:ab
znt=z_next(znt) LOgB_$9_3
end_loop +RK/u
Eh"Y<]$
end 0
eOdE+
get_znum # X6=`Xe#
P# ;pQC
}Fu2%L>
def weibull_dis ;weibull分布 N Cm=l
Yt r*"-
znt=zone_head P)ne^ _
loop while znt # null 3[~LmA
weibull_young=((-ln(urand))^(1.0/emod_weib_m))*ex_young ;弹性模量weibull分布 P+3
]g{2w
weibull_thexp=((-ln(urand))^(1.0/exth_weib_m))*ex_thexp ;热膨胀系数weibull分布 !=+;9Ry$z
s GrI%3[e"
;----------某个量不变化时候时候-------------- h^J :k
weibull_young=18.0000e9 ;弹性模量不变化 |6`7kb;p
;weibull_thexp=4.6200e-5 ;膨胀系数不变化 HS1{4/
;------------------------------------- q@4Cw&AI+
~fO#En
bulk_mod=weibull_young/3.0/(1.0-2*p_rat) `"$9L[>
shear_mod=weibull_young/2.0/(1.0+p_rat) >q&e.- qL
*LvdrPxU=
z_prop(znt,'bulk')=bulk_mod ;给每个单元赋参数 TNcMr bWA
z_prop(znt,'shear')=shear_mod 0+A#k7c6p
O [=W%2I!i
z_prop(znt,'thexp')=weibull_thexp /4+*!X
T$V8 n_;
z_gnum=z_gnum+1 rPHM_fW(O@
)P.,h&h/
table(1,z_gnum)=weibull_young 2[KHmdgtB
table(2,z_gnum)=weibull_thexp !Z; Nv
zS?DXE
znt=z_next(znt) ^SwU]e
end_loop E/E|*6R
end #;"lBqxY`
BJsN~` =r
;-------------------均匀系数、尺度参数、变异系数------------------ . ;D'
;------------------弹性模量-------------------------- zaE !=-U
;1.5,19.94018,0.6790---2.0,20.31144,0.5228---4.0,19.85878,0.2805---8.0,19.11437,0.1882 t "6u
;15, 18.63933,0.0810---20 18.489980.0625---40, 18.25188,0.0306---100,18.10319,0.0130 zXd#kw;
I*(7(>zgyv
set emod_weib_m=1.5 ;弹摸均匀系数 XK(aH~7xme
set ex_young=19.94018e9 ;尺度参数 1DF8-|+
set p_rat=0.3 ;泊松比 vLCm,Bb2L
>DW%i\k1V~
;--------------------膨胀系数----------------------- *#'&a(h B!
;1.5,5.11798,0.6790---2.0,5.21327,0.5228---4.0,5.09709,0.2805---8.0,4.90602,0.1882 tTJ$tx
;15, 4.78409,0.0810---20 4.74576,0.0625---40, 4.68465,0.0306---100,4.64648,0.0130 BkDq9>
L-e6^%eU
set exth_weib_m=1.5 ;膨胀系数均匀系数 Ge76/T%{Q
set ex_thexp=5.11798e-5 ;尺度参数 'qidorT>N
{_4zm&
hllb\Y)XL
weibull_dis
