我觉得下面这本书不错，有料，写得用心。好东西和大家分享一下：
lo:~a J8 j$P I,` 《计算土力学》应该是本严谨的书，我向张教授问书中附录程序代码，张教授回复给我一个下载网址：
http://www.cm.nitech.ac.jp/cho/models/models.html。下面是其中一段代码，大家觉得还行的话，可以到上面网址去下载。
T VkcDS *\q8BZ Appendix IV 软岩的弹粘塑性模型的计算程序
Slx2 z%'> 为了区分数字Ô和字母O,所有的数字0均用Ô表示。
SMO*({/ A-IV-1 主程序
"iek,Y}j7 IMPLICIT REAL*8(A-H,O-Z)
3M*Y= ?pI C
H@W0gK(cS; PARAMETER (IMZ=1)
rO^xz7K^ C
b(&] >z COMMON NN
^, wnp@ COMMON /TN1/TNONE(IMZ),VP(IMZ),VPÔ(IMZ)
l,bZG3,6 COMMON /CRITICAL/S1S3(IMZ)
Mq+< mX7 COMMON /TNSAND/DEPSON(6),E(3)
l%PnB
)F COMMON /TNSAND1/LTYPE
OmNn,PCl8 COMMON /TNSAND2/RATEÔ
M&e8zS DIMENSION DEP(6,6),SIG(6),SIGT(6),DE(6),DEPP(6),DSIG(6)
,+-h 7^{` C
1n3$V:00 OPEN(1,FILE='in.dat')
e^frVEV OPEN(2,FILE='out.txt')
lU$X4JBzS C
/0\QL+^! READ(1,*)LTYPE1,NNN
s=/^lOOO READ(1,*)POI,EPSÔ,ZMF,ZRDK,YOUNG,BETA,AN,ALFA,CNN,OCR
dq1TRFu READ(1,*)SIGMAÔ,STRAIN,RATE,SIGMQR
B$7[8h C
_3&/(B%H C
ZR
mPP C
J|?[.h7tO C POI :ν Poisson's ratio
;Jo*|pju C EPSÔ :eÔ when σ'm=98kPa Reference void ratio
*`.LA@bHU C ZMF :Rf=σ1/σ3 Stress ratio at critical state
T{Hf P C ZRDK :λ-κ Difference between compression and swelling index
)@$
&FFIu C YOUNG :E Young's modulus
#M:B3C!ouY C BETA :β Parameter of tij rock subloading time-dependent model
2}!R
T C ANN :a Parameter involved in Cam-Clay subloading model &
3 _U\VGm C modified Cam-Clay model & tij subloading model &
bL
MkPty C tij rock subloading time-dependent model
Nh :4ys!P C ALFA : 二次固结系数
\E&th p C CNN :C 二次固结系数
B$`d&7I;D C SIGMAÔ :σ'mÔ Initial mean effective stress
wRi~Yb? C OCR :OCR Overconsolidation ratio
kPedX C STRAIN :应变制限値 单向加载时最大的应变；交变加载时应变的振幅
`axQd%:AC C STRESS :应力制限值 交变加载时应力的振幅
`&,_xUA C SIGMQR :蠕变应力
1:5P%$? b C
|rFJ*.nD C LTYPE1=1: 中空扭转剪切(排水试验:应力一定;非排水试验:垂直位移一定)
X&,N}9>B C LTYPE1=2: 最小主应力一定的三轴压缩
v?{vg?vI C LTYPE1=3: 平均主应力一定的三轴压缩
7YD\ !2b C LTYPE1=4: 平面应变状态侧向围压一定的压缩((SIG(1)-SIG(3)))
l]>!`'sJL C LTYPE1=5: 等向压缩(SIG(1)=SIG(2)=SIG(3))
8E=vR 8 C
~=Fk/ C
o[imNy~ ~ C NNN:单向加载时达到最大应变时计算的步数；交变加载时单向振幅的计算步数
*V<)p%l. C NTOTAL：全部计算的步数
#4l HaFq C NSPACE: 输出时的步数间隔
=x oBC&u \`&xprqAw 。。。。。。
