我觉得下面这本书不错，有料，写得用心。好东西和大家分享一下：
k)\Yl`4au 1fRYXqx 《计算土力学》应该是本严谨的书，我向张教授问书中附录程序代码，张教授回复给我一个下载网址：
http://www.cm.nitech.ac.jp/cho/models/models.html。下面是其中一段代码，大家觉得还行的话，可以到上面网址去下载。
]D&$k P( yBO88rfh> Appendix IV 软岩的弹粘塑性模型的计算程序
+s&+G![ 为了区分数字Ô和字母O,所有的数字0均用Ô表示。
b\~rL,7( A-IV-1 主程序
$5nOi aQL IMPLICIT REAL*8(A-H,O-Z)
YBQ O ]3f C
vW6Pf^yJ PARAMETER (IMZ=1)
!!f)w!wW C
-Zfzl`r COMMON NN
_5#f9,m1 COMMON /TN1/TNONE(IMZ),VP(IMZ),VPÔ(IMZ)
(_-<3)q4 COMMON /CRITICAL/S1S3(IMZ)
A" ` 6 2 COMMON /TNSAND/DEPSON(6),E(3)
0}$",M !p COMMON /TNSAND1/LTYPE
[ >#?C*s COMMON /TNSAND2/RATEÔ
Z[KXDQn8 DIMENSION DEP(6,6),SIG(6),SIGT(6),DE(6),DEPP(6),DSIG(6)
\> nPg5OT C
VR5$[-E3 OPEN(1,FILE='in.dat')
{ /12.y=)~ OPEN(2,FILE='out.txt')
#
4`* `)% C
Lu}oC2 READ(1,*)LTYPE1,NNN
7~T E=t READ(1,*)POI,EPSÔ,ZMF,ZRDK,YOUNG,BETA,AN,ALFA,CNN,OCR
cl{x5>.'# READ(1,*)SIGMAÔ,STRAIN,RATE,SIGMQR
2J rr;"r C
&!J X
C
"XMTj <D C
Hx2.2 A^ C POI :ν Poisson's ratio
[>&Nhn0 iY C EPSÔ :eÔ when σ'm=98kPa Reference void ratio
":?>6'*1 C ZMF :Rf=σ1/σ3 Stress ratio at critical state
a_'W1ek-@ C ZRDK :λ-κ Difference between compression and swelling index
wG2lCv`d C YOUNG :E Young's modulus
Y>Q9? >}Q C BETA :β Parameter of tij rock subloading time-dependent model
-wlob `3 C ANN :a Parameter involved in Cam-Clay subloading model &
3toY #!1Ch C modified Cam-Clay model & tij subloading model &
h;f5@#F C tij rock subloading time-dependent model
.z7X Ymv C ALFA : 二次固结系数
& PrV+Lv C CNN :C 二次固结系数
\\
M2_mT C SIGMAÔ :σ'mÔ Initial mean effective stress
)zI<C=])" C OCR :OCR Overconsolidation ratio
( :]. ?o C STRAIN :应变制限値 单向加载时最大的应变；交变加载时应变的振幅
|332G64K C STRESS :应力制限值 交变加载时应力的振幅
+SkD/"5ng C SIGMQR :蠕变应力
\ntmD?kA C
QsI#Ae,O#; C LTYPE1=1: 中空扭转剪切(排水试验:应力一定;非排水试验:垂直位移一定)
j2deb`GD C LTYPE1=2: 最小主应力一定的三轴压缩
'(Uyju= C LTYPE1=3: 平均主应力一定的三轴压缩
g1muT.W]S C LTYPE1=4: 平面应变状态侧向围压一定的压缩((SIG(1)-SIG(3)))
_WtX8 C LTYPE1=5: 等向压缩(SIG(1)=SIG(2)=SIG(3))
4JFi|oK0H C
cEi{+rfZd| C
d
/&aC#'B C NNN:单向加载时达到最大应变时计算的步数；交变加载时单向振幅的计算步数
U 2/H,D C NTOTAL：全部计算的步数
3SVI|A5(d C NSPACE: 输出时的步数间隔
=A9>Ej/ Y }'C'PR 。。。。。。
