我觉得下面这本书不错，有料，写得用心。好东西和大家分享一下：
+ # GQ, $}lbT15a 《计算土力学》应该是本严谨的书，我向张教授问书中附录程序代码，张教授回复给我一个下载网址：
http://www.cm.nitech.ac.jp/cho/models/models.html。下面是其中一段代码，大家觉得还行的话，可以到上面网址去下载。
q2. XoCf /E
Bo3` Appendix IV 软岩的弹粘塑性模型的计算程序
@.G[s)x 为了区分数字Ô和字母O,所有的数字0均用Ô表示。
cU ? 0(z7 A-IV-1 主程序
<.pU,T/ IMPLICIT REAL*8(A-H,O-Z)
M (jgd C
eAX
)^q PARAMETER (IMZ=1)
s >e=?W C
)p
T?/ J COMMON NN
Wi[ ~fI8^! COMMON /TN1/TNONE(IMZ),VP(IMZ),VPÔ(IMZ)
rrQQZ5fh b COMMON /CRITICAL/S1S3(IMZ)
b|cyjDMAA COMMON /TNSAND/DEPSON(6),E(3)
hfnN@Kg?B} COMMON /TNSAND1/LTYPE
K#*reJ}K COMMON /TNSAND2/RATEÔ
5nv<^>[J DIMENSION DEP(6,6),SIG(6),SIGT(6),DE(6),DEPP(6),DSIG(6)
!lEY=1nHOJ C
|_o=^?z' OPEN(1,FILE='in.dat')
(:._"jp] OPEN(2,FILE='out.txt')
qP{/[uj[K C
0dhF&*h|L READ(1,*)LTYPE1,NNN
.% .9n\b READ(1,*)POI,EPSÔ,ZMF,ZRDK,YOUNG,BETA,AN,ALFA,CNN,OCR
C K:y? READ(1,*)SIGMAÔ,STRAIN,RATE,SIGMQR
,stN C
Yiry[ "[]Q C
)ap_Z6 C
T _sTC) &a C POI :ν Poisson's ratio
+
` s@ C EPSÔ :eÔ when σ'm=98kPa Reference void ratio
B+[A]dgS C ZMF :Rf=σ1/σ3 Stress ratio at critical state
# ?q&r_@ @ C ZRDK :λ-κ Difference between compression and swelling index
/GIxR6i C YOUNG :E Young's modulus
V2$h8\a C BETA :β Parameter of tij rock subloading time-dependent model
':gUOra|I C ANN :a Parameter involved in Cam-Clay subloading model &
RRmL d/( C modified Cam-Clay model & tij subloading model &
fQ/
0R C tij rock subloading time-dependent model
T?:glp[4I C ALFA : 二次固结系数
f ` :i.Sr C CNN :C 二次固结系数
D@.tkzU@E C SIGMAÔ :σ'mÔ Initial mean effective stress
L !=4N!j C OCR :OCR Overconsolidation ratio
{LX.iH9}l C STRAIN :应变制限値 单向加载时最大的应变；交变加载时应变的振幅
_7IKzUn9g[ C STRESS :应力制限值 交变加载时应力的振幅
[Q M u2 C SIGMQR :蠕变应力
A/s>PhxV C
Sl-v W C LTYPE1=1: 中空扭转剪切(排水试验:应力一定;非排水试验:垂直位移一定)
M7+nW ; e% C LTYPE1=2: 最小主应力一定的三轴压缩
zXv3:uRp. C LTYPE1=3: 平均主应力一定的三轴压缩
Ul 2R'"FB C LTYPE1=4: 平面应变状态侧向围压一定的压缩((SIG(1)-SIG(3)))
e_s &L,ze C LTYPE1=5: 等向压缩(SIG(1)=SIG(2)=SIG(3))
._8K suJG C
?47@ o1 C
A]Y V s C NNN:单向加载时达到最大应变时计算的步数；交变加载时单向振幅的计算步数
qtiz a~u C NTOTAL：全部计算的步数
YGv<VOWG2 C NSPACE: 输出时的步数间隔
4 !+pc-}- &0 7]LF$] 。。。。。。
