| analyze |
2008-07-23 17:39 |
张峰教授《计算土力学》书中附录程序代码下载地址
我觉得下面这本书不错，有料，写得用心。好东西和大家分享一下： =+-Yxh|*
GP&vLt51
《计算土力学》应该是本严谨的书，我向张教授问书中附录程序代码，张教授回复给我一个下载网址：http://www.cm.nitech.ac.jp/cho/models/models.html。下面是其中一段代码，大家觉得还行的话，可以到上面网址去下载。 AtF3%Z v2
DyeQJ7p
Appendix IV 软岩的弹粘塑性模型的计算程序 *uo'VJI7_,
为了区分数字Ô和字母O,所有的数字0均用Ô表示。 <$ F\Nk|x
A-IV-1 主程序 K'X2dG*
IMPLICIT REAL*8(A-H,O-Z) |N }P(GF
C 5B&;uY
PARAMETER (IMZ=1) *x ON W
C % ]I ZLJ
COMMON NN yaG= j
COMMON /TN1/TNONE(IMZ),VP(IMZ),VPÔ(IMZ) G:pEE:W[
COMMON /CRITICAL/S1S3(IMZ) ^5A
t?I8
COMMON /TNSAND/DEPSON(6),E(3) q EP
4
COMMON /TNSAND1/LTYPE abJ"
[
COMMON /TNSAND2/RATEÔ qf=1?=l291
DIMENSION DEP(6,6),SIG(6),SIGT(6),DE(6),DEPP(6),DSIG(6) ^| /] (
C &g.@u~SI1
OPEN(1,FILE='in.dat') d'/TdVM
OPEN(2,FILE='out.txt') \0mb
3Q'
C vfl5Mx4
READ(1,*)LTYPE1,NNN H"C[&r
READ(1,*)POI,EPSÔ,ZMF,ZRDK,YOUNG,BETA,AN,ALFA,CNN,OCR jS+AGE?5e
READ(1,*)SIGMAÔ,STRAIN,RATE,SIGMQR N nk@h
C Ux#x#N
C a)S+8uU
C `2`\]X_A{
C POI :ν Poisson's ratio ^2$ lJ
C EPSÔ :eÔ when σ'm=98kPa Reference void ratio 7<jZ`qd q_
C ZMF :Rf=σ1/σ3 Stress ratio at critical state OM|Fwr$
C ZRDK :λ-κ Difference between compression and swelling index Pt&(npjN,
C YOUNG :E Young's modulus I
j$lDJS
C BETA :β Parameter of tij rock subloading time-dependent model $uap8nN
C ANN :a Parameter involved in Cam-Clay subloading model & sf]s",t~J
C modified Cam-Clay model & tij subloading model & )r[&RGz6
C tij rock subloading time-dependent model c-g)eV|)S
C ALFA : 二次固结系数 $&"V ^@
C CNN :C 二次固结系数 , SQZD,3v4
C SIGMAÔ :σ'mÔ Initial mean effective stress `E!N9qI?t$
C OCR :OCR Overconsolidation ratio 7C$
5
C STRAIN :应变制限値 单向加载时最大的应变；交变加载时应变的振幅 ov wQ2TuK
C STRESS :应力制限值 交变加载时应力的振幅 # [ e
C SIGMQR :蠕变应力 <L{(M j%Z
C {0;3W7
C LTYPE1=1: 中空扭转剪切(排水试验:应力一定;非排水试验:垂直位移一定) ? W( 6
C LTYPE1=2: 最小主应力一定的三轴压缩 u5~Ns&o&N
C LTYPE1=3: 平均主应力一定的三轴压缩 Qb! PRCHQ
C LTYPE1=4: 平面应变状态侧向围压一定的压缩((SIG(1)-SIG(3))) WHAQu]{
C LTYPE1=5: 等向压缩(SIG(1)=SIG(2)=SIG(3)) A<6%r7&B'
C f>s?4
C yA) (*PFz
C NNN:单向加载时达到最大应变时计算的步数；交变加载时单向振幅的计算步数 v^ /Q 8Q
C NTOTAL：全部计算的步数 P7
PB t
C NSPACE: 输出时的步数间隔 F q!fWl
MwfOy@|N
。。。。。。 |
|
