我觉得下面这本书不错，有料，写得用心。好东西和大家分享一下：
Wv||9[Rd Eh*(N (` 《计算土力学》应该是本严谨的书，我向张教授问书中附录程序代码，张教授回复给我一个下载网址：
http://www.cm.nitech.ac.jp/cho/models/models.html。下面是其中一段代码，大家觉得还行的话，可以到上面网址去下载。
sa TS8p z c 2gi 3 Appendix IV 软岩的弹粘塑性模型的计算程序
"hk {"0E 为了区分数字Ô和字母O,所有的数字0均用Ô表示。
r{ K Q3j9O A-IV-1 主程序
=~>g--^U IMPLICIT REAL*8(A-H,O-Z)
!/e*v>3u& C
[1nUq!uTm PARAMETER (IMZ=1)
= Z ql6D C
H?rC I S0 COMMON NN
;!A=YXB COMMON /TN1/TNONE(IMZ),VP(IMZ),VPÔ(IMZ)
h[mJ=LIrg COMMON /CRITICAL/S1S3(IMZ)
(K_{a+$[ COMMON /TNSAND/DEPSON(6),E(3)
@N ]]Cf>x COMMON /TNSAND1/LTYPE
7,zE?KG / COMMON /TNSAND2/RATEÔ
&["e1k i DIMENSION DEP(6,6),SIG(6),SIGT(6),DE(6),DEPP(6),DSIG(6)
qw
03] a C
J0Gjo9L OPEN(1,FILE='in.dat')
pI_dV44W OPEN(2,FILE='out.txt')
YmCu\+ u C
*M/ :W =,t READ(1,*)LTYPE1,NNN
m(D+!I9 READ(1,*)POI,EPSÔ,ZMF,ZRDK,YOUNG,BETA,AN,ALFA,CNN,OCR
]!j%Ad READ(1,*)SIGMAÔ,STRAIN,RATE,SIGMQR
hd]ts. C
;W]9DBAB C
n,d)Wwe_`y C
\ef:H& r C POI :ν Poisson's ratio
,q;?zcC7 C EPSÔ :eÔ when σ'm=98kPa Reference void ratio
<_c8F!K)T C ZMF :Rf=σ1/σ3 Stress ratio at critical state
md,KRE C ZRDK :λ-κ Difference between compression and swelling index
>g m C YOUNG :E Young's modulus
jweX"G54R C BETA :β Parameter of tij rock subloading time-dependent model
J;=T"C& C ANN :a Parameter involved in Cam-Clay subloading model &
wh)F&@6 R! C modified Cam-Clay model & tij subloading model &
[r! f&R C tij rock subloading time-dependent model
je5[.VT M C ALFA : 二次固结系数
7MGvw-Tpb7 C CNN :C 二次固结系数
&isKU 8n
C SIGMAÔ :σ'mÔ Initial mean effective stress
B$n 1 k 45 C OCR :OCR Overconsolidation ratio
F0~<p[9Nx C STRAIN :应变制限値 单向加载时最大的应变；交变加载时应变的振幅
f !#+cM C STRESS :应力制限值 交变加载时应力的振幅
}m'n1tm;
C SIGMQR :蠕变应力
i\O^s ] C
D *j^f7ab C LTYPE1=1: 中空扭转剪切(排水试验:应力一定;非排水试验:垂直位移一定)
\2]_NU5. C LTYPE1=2: 最小主应力一定的三轴压缩
Tw7] C LTYPE1=3: 平均主应力一定的三轴压缩
~GcWG4 C LTYPE1=4: 平面应变状态侧向围压一定的压缩((SIG(1)-SIG(3)))
lh[?`+A C LTYPE1=5: 等向压缩(SIG(1)=SIG(2)=SIG(3))
>ZW|wpO C
B%tWi C
6']HmM C NNN:单向加载时达到最大应变时计算的步数；交变加载时单向振幅的计算步数
[(gXjt- C NTOTAL：全部计算的步数
LM+d3|gSV C NSPACE: 输出时的步数间隔
NJ]3qH xK_oV+ 。。。。。。
