;new U%0Ty|$Y
hsu{ey p
;restore therm_model.sav o yo(1 >
Jk,; JQ
;model mohr ;本构模型 O-)-YVU
;model th_iso tT d\|
vn/.}GkpU
def get_znum ;获取单元总数 8L w B
B
% `Z! 4L
znt=zone_head 26CS6( sn
loop while znt # null s+tGFjq
znum=znum+1 q -8t'7
znt=z_next(znt) pW4 cX
end_loop g^zs,4pPU<
}i!J/tJ)b
end F;^F+H
get_znum @k&qb!Qah
QeF3qXI
vloF::1
def weibull_dis ;weibull分布 2w F8 P)
TX
znt=zone_head \d::l{VB
loop while znt # null Sna4wkbS
weibull_young=((-ln(urand))^(1.0/emod_weib_m))*ex_young ;弹性模量weibull分布 YJ$1N!rG
weibull_thexp=((-ln(urand))^(1.0/exth_weib_m))*ex_thexp ;热膨胀系数weibull分布 [9: 9Ql_h
Jmx Ko+-
;----------某个量不变化时候时候-------------- XrZ*1V
weibull_young=18.0000e9 ;弹性模量不变化 arf8xqR-U]
;weibull_thexp=4.6200e-5 ;膨胀系数不变化 #m={yck *
;------------------------------------- sVex
(X
-_$$ Te
bulk_mod=weibull_young/3.0/(1.0-2*p_rat) F2X0%te
shear_mod=weibull_young/2.0/(1.0+p_rat) m.,U:>
K _&4D'
z_prop(znt,'bulk')=bulk_mod ;给每个单元赋参数 3G kv4,w<
z_prop(znt,'shear')=shear_mod 1'? 4m0W1
`p+Zz"/
z_prop(znt,'thexp')=weibull_thexp \l~*PG2
l& ?ii68/
z_gnum=z_gnum+1 7`u$
X.:]=,aGW
table(1,z_gnum)=weibull_young ->N8#XH2=
table(2,z_gnum)=weibull_thexp k1Q ? '<`
rx}r~0 i
znt=z_next(znt) 5atYOe p
end_loop 0wS+++n$5
end H /}]FmjN
r(9~$_ (vK
;-------------------均匀系数、尺度参数、变异系数------------------ "Wx]RN:
;------------------弹性模量-------------------------- <_Q1k>
;1.5,19.94018,0.6790---2.0,20.31144,0.5228---4.0,19.85878,0.2805---8.0,19.11437,0.1882 |fo0
;15, 18.63933,0.0810---20 18.489980.0625---40, 18.25188,0.0306---100,18.10319,0.0130 !W?gR.0$=
f5 ` g
set emod_weib_m=1.5 ;弹摸均匀系数 g:EVhuK
set ex_young=19.94018e9 ;尺度参数 hJPlq0C
set p_rat=0.3 ;泊松比 ZRYEqSm
^Ox3XC
;--------------------膨胀系数----------------------- Z@r.pRr'
;1.5,5.11798,0.6790---2.0,5.21327,0.5228---4.0,5.09709,0.2805---8.0,4.90602,0.1882 9?k_y ZV
;15, 4.78409,0.0810---20 4.74576,0.0625---40, 4.68465,0.0306---100,4.64648,0.0130 j)q\9#sI/(
)vxUT{; sH
set exth_weib_m=1.5 ;膨胀系数均匀系数 \QB;Ja _
set ex_thexp=5.11798e-5 ;尺度参数 ]kboG%Dl?9
|ZQ@fmvL/p
c/ uNM
weibull_dis
