我觉得下面这本书不错，有料，写得用心。好东西和大家分享一下：
Y=D\ 5c$\DZ( 《计算土力学》应该是本严谨的书，我向张教授问书中附录程序代码，张教授回复给我一个下载网址：
http://www.cm.nitech.ac.jp/cho/models/models.html。下面是其中一段代码，大家觉得还行的话，可以到上面网址去下载。
1D1qOg"LE 8cGoo u6 Appendix IV 软岩的弹粘塑性模型的计算程序
Sy m}#F\s 为了区分数字Ô和字母O,所有的数字0均用Ô表示。
~\+Bb8+hpJ A-IV-1 主程序
$|t={s34 IMPLICIT REAL*8(A-H,O-Z)
OMAvJzK . C
&qP0-x) PARAMETER (IMZ=1)
%+oqAY m+s C
\. a 7F4h COMMON NN
[1e]_9)p COMMON /TN1/TNONE(IMZ),VP(IMZ),VPÔ(IMZ)
x/ix%!8J COMMON /CRITICAL/S1S3(IMZ)
> D% COMMON /TNSAND/DEPSON(6),E(3)
|nBs(>b COMMON /TNSAND1/LTYPE
iKu4s COMMON /TNSAND2/RATEÔ
w+$~ ds DIMENSION DEP(6,6),SIG(6),SIGT(6),DE(6),DEPP(6),DSIG(6)
LfN,aW C
R>B6@|}? OPEN(1,FILE='in.dat')
PYTwyqS OPEN(2,FILE='out.txt')
Py# TXzEcC C
MzH'<`;BP READ(1,*)LTYPE1,NNN
}#S1!TU READ(1,*)POI,EPSÔ,ZMF,ZRDK,YOUNG,BETA,AN,ALFA,CNN,OCR
J,J6bfR/ READ(1,*)SIGMAÔ,STRAIN,RATE,SIGMQR
0i>p1/kv C
P !I Lji! C
O l; DJV C
iK!F VKi} C POI :ν Poisson's ratio
mbns%%GJU C EPSÔ :eÔ when σ'm=98kPa Reference void ratio
K>TEt5 C ZMF :Rf=σ1/σ3 Stress ratio at critical state
l{*Ko~g C ZRDK :λ-κ Difference between compression and swelling index
1024L; C YOUNG :E Young's modulus
-us:!p1T C BETA :β Parameter of tij rock subloading time-dependent model
QH_Ds,oH= C ANN :a Parameter involved in Cam-Clay subloading model &
2"c 5 < C modified Cam-Clay model & tij subloading model &
@w;$M]o1 C tij rock subloading time-dependent model
R_b4S%jhx C ALFA : 二次固结系数
M 5$JB nN C CNN :C 二次固结系数
q:=jv6T# C SIGMAÔ :σ'mÔ Initial mean effective stress
4s@Tn>%SP C OCR :OCR Overconsolidation ratio
5?[hr5E.E C STRAIN :应变制限値 单向加载时最大的应变；交变加载时应变的振幅
jVh:Bw C STRESS :应力制限值 交变加载时应力的振幅
GK;IY =8W C SIGMQR :蠕变应力
7A h C
a#l ytp C LTYPE1=1: 中空扭转剪切(排水试验:应力一定;非排水试验:垂直位移一定)
{yDQncq'^ C LTYPE1=2: 最小主应力一定的三轴压缩
X#HH7V> C LTYPE1=3: 平均主应力一定的三轴压缩
$%4< q0 - C LTYPE1=4: 平面应变状态侧向围压一定的压缩((SIG(1)-SIG(3)))
#8 ~ygEa} C LTYPE1=5: 等向压缩(SIG(1)=SIG(2)=SIG(3))
7$x%A&] C
Yr>0Qg], C
iC Ce8nK C NNN:单向加载时达到最大应变时计算的步数；交变加载时单向振幅的计算步数
YRMe<upo C NTOTAL：全部计算的步数
w[ $oH^7 C NSPACE: 输出时的步数间隔
DP;:%L} P-<1vfThH 。。。。。。
