;new pIhy3@bY
L#~z#
;restore therm_model.sav +?%LX4Y
[dIlt"2fV
;model mohr ;本构模型 @s%!R
;model th_iso DQ`\HY
7i##g,
def get_znum ;获取单元总数 0K:3?Ik
y{#9&ct&
znt=zone_head K)&oDwk
loop while znt # null x@Ze%$'
znum=znum+1 tPz!C&.=
znt=z_next(znt) ',l}$]y5
end_loop nl/~7({
A>B_~=
end 9k*1_
get_znum 2S@aG%-)
grZ?F~P8
=8W'4MC
def weibull_dis ;weibull分布 +6;1.5Tc
LQ{z}Ay
znt=zone_head SLjSNuOP
loop while znt # null ;a9`z+ K
weibull_young=((-ln(urand))^(1.0/emod_weib_m))*ex_young ;弹性模量weibull分布 Z` zyEP A
weibull_thexp=((-ln(urand))^(1.0/exth_weib_m))*ex_thexp ;热膨胀系数weibull分布 (&@,Z I;
ud$*/ )/
;----------某个量不变化时候时候-------------- )Rj,PF-9Z[
weibull_young=18.0000e9 ;弹性模量不变化 nJr:U2d
;weibull_thexp=4.6200e-5 ;膨胀系数不变化 \wb0%>
0
;------------------------------------- 3cB=9Y{<
j= vlsW
bulk_mod=weibull_young/3.0/(1.0-2*p_rat) kNC]q,ljt5
shear_mod=weibull_young/2.0/(1.0+p_rat) I%'6IpR"d
=67tQx58
z_prop(znt,'bulk')=bulk_mod ;给每个单元赋参数 .[:2M9Rx
z_prop(znt,'shear')=shear_mod @q++eGm\Q
SUaXm#9
z_prop(znt,'thexp')=weibull_thexp S&z8-D=8k
3aERfIJyE
z_gnum=z_gnum+1 Q$:>yveR*
U%4s@{7
table(1,z_gnum)=weibull_young cqZuG}VR
table(2,z_gnum)=weibull_thexp q6AL}9]9
]s>y se
znt=z_next(znt) T(q/$p&q
end_loop .]y"04@]
end =&x
u"V
@lu`oyM
;-------------------均匀系数、尺度参数、变异系数------------------ .3:s4=(f
;------------------弹性模量-------------------------- CZ|h` ";P2
;1.5,19.94018,0.6790---2.0,20.31144,0.5228---4.0,19.85878,0.2805---8.0,19.11437,0.1882 lCxPR'C|
;15, 18.63933,0.0810---20 18.489980.0625---40, 18.25188,0.0306---100,18.10319,0.0130 IRY/0v
~`!{5:v
set emod_weib_m=1.5 ;弹摸均匀系数 %O02xr=
set ex_young=19.94018e9 ;尺度参数 q7aH=dhw
set p_rat=0.3 ;泊松比 0oNy
K['Gp>l
;--------------------膨胀系数----------------------- +I[Hxf ~
;1.5,5.11798,0.6790---2.0,5.21327,0.5228---4.0,5.09709,0.2805---8.0,4.90602,0.1882 >>R,P
Ow-
;15, 4.78409,0.0810---20 4.74576,0.0625---40, 4.68465,0.0306---100,4.64648,0.0130 :#nv:~2]
v7#`b}'W
set exth_weib_m=1.5 ;膨胀系数均匀系数 |N+uEiJ
set ex_thexp=5.11798e-5 ;尺度参数 WP ~]pduT
^w+)A;?W
r~Vb*~U"
weibull_dis