我觉得下面这本书不错，有料，写得用心。好东西和大家分享一下：
?+$EPaC2 'C;K Nc 《计算土力学》应该是本严谨的书，我向张教授问书中附录程序代码，张教授回复给我一个下载网址：
http://www.cm.nitech.ac.jp/cho/models/models.html。下面是其中一段代码，大家觉得还行的话，可以到上面网址去下载。
KqNsCT+j 6 \ %#=GG Appendix IV 软岩的弹粘塑性模型的计算程序
RER93:( 为了区分数字Ô和字母O,所有的数字0均用Ô表示。
Jq+@%#G A-IV-1 主程序
%WYveY IMPLICIT REAL*8(A-H,O-Z)
@[n% q.|VB C
$. sTb PARAMETER (IMZ=1)
q~X}&}UT C
52F3 r:Rk COMMON NN
Q qcAmp COMMON /TN1/TNONE(IMZ),VP(IMZ),VPÔ(IMZ)
q ^Gj
IP COMMON /CRITICAL/S1S3(IMZ)
M?kXzb \O COMMON /TNSAND/DEPSON(6),E(3)
>R.!Qze\G COMMON /TNSAND1/LTYPE
Yi19VU|/ COMMON /TNSAND2/RATEÔ
): r'IR DIMENSION DEP(6,6),SIG(6),SIGT(6),DE(6),DEPP(6),DSIG(6)
G B>T3l" C
W;U<,g
' OPEN(1,FILE='in.dat')
3ZvQUH/{W OPEN(2,FILE='out.txt')
N'|9rB2e C
v{8r46Y~Z) READ(1,*)LTYPE1,NNN
J`wx72/-ZW READ(1,*)POI,EPSÔ,ZMF,ZRDK,YOUNG,BETA,AN,ALFA,CNN,OCR
/)rv Ndn READ(1,*)SIGMAÔ,STRAIN,RATE,SIGMQR
U;gy4rj C
=H^~"16 C
k_Lv\'Ok C
(: mF+%( C POI :ν Poisson's ratio
s&DAO r!i C EPSÔ :eÔ when σ'm=98kPa Reference void ratio
JqEo~]E] C ZMF :Rf=σ1/σ3 Stress ratio at critical state
dQ #oY|a C ZRDK :λ-κ Difference between compression and swelling index
.tK]-f2 C YOUNG :E Young's modulus
H{_6e6`e. C BETA :β Parameter of tij rock subloading time-dependent model
SK_N|X]. C ANN :a Parameter involved in Cam-Clay subloading model &
z Cu+Oi6 C modified Cam-Clay model & tij subloading model &
0,iG9D 7 C tij rock subloading time-dependent model
eEeK ] 8@ C ALFA : 二次固结系数
[kPl7[ OL C CNN :C 二次固结系数
9A}y^=!` C SIGMAÔ :σ'mÔ Initial mean effective stress
N I C.c3 C OCR :OCR Overconsolidation ratio
Xj:\B] v] C STRAIN :应变制限値 单向加载时最大的应变；交变加载时应变的振幅
9D yy& $s C STRESS :应力制限值 交变加载时应力的振幅
cO-^#di C SIGMQR :蠕变应力
q@Zeu\T,*# C
0_t9;;y : C LTYPE1=1: 中空扭转剪切(排水试验:应力一定;非排水试验:垂直位移一定)
nzU0=w}V C LTYPE1=2: 最小主应力一定的三轴压缩
aDE}'d1qo C LTYPE1=3: 平均主应力一定的三轴压缩
0f;|0siTAm C LTYPE1=4: 平面应变状态侧向围压一定的压缩((SIG(1)-SIG(3)))
^HHT>K-m C LTYPE1=5: 等向压缩(SIG(1)=SIG(2)=SIG(3))
u0$}VO5/a C
|vo Z0U C
wqyF"^It" C NNN:单向加载时达到最大应变时计算的步数；交变加载时单向振幅的计算步数
lO}I>yo}\ C NTOTAL：全部计算的步数
ntZHO}' C NSPACE: 输出时的步数间隔
fOk(i vYy C4|79UG >s 。。。。。。
