我觉得下面这本书不错，有料，写得用心。好东西和大家分享一下：
=_`4H Dr T/" 6iv\1 《计算土力学》应该是本严谨的书，我向张教授问书中附录程序代码，张教授回复给我一个下载网址：
http://www.cm.nitech.ac.jp/cho/models/models.html。下面是其中一段代码，大家觉得还行的话，可以到上面网址去下载。
%jK-}0Tu wPDA_ns~ Appendix IV 软岩的弹粘塑性模型的计算程序
"mnW qRpX 为了区分数字Ô和字母O,所有的数字0均用Ô表示。
PEPBnBA&1 A-IV-1 主程序
CqUK[#kW( IMPLICIT REAL*8(A-H,O-Z)
{iqH 27\E C
#L{QnV.3 PARAMETER (IMZ=1)
JBCcR,\kM* C
38dXfl COMMON NN
i !jZZj-{ COMMON /TN1/TNONE(IMZ),VP(IMZ),VPÔ(IMZ)
&WIi w$@ COMMON /CRITICAL/S1S3(IMZ)
6=kA COMMON /TNSAND/DEPSON(6),E(3)
xS` %3+| COMMON /TNSAND1/LTYPE
Z,Z34:- COMMON /TNSAND2/RATEÔ
Y?2I
/ DIMENSION DEP(6,6),SIG(6),SIGT(6),DE(6),DEPP(6),DSIG(6)
XQhBnam%
C
) =<,$ |g OPEN(1,FILE='in.dat')
&j,rq?eh$ OPEN(2,FILE='out.txt')
FrO)3 1z C
nJ#uz: (w, READ(1,*)LTYPE1,NNN
VpJ/M(UD- READ(1,*)POI,EPSÔ,ZMF,ZRDK,YOUNG,BETA,AN,ALFA,CNN,OCR
: luVsQ READ(1,*)SIGMAÔ,STRAIN,RATE,SIGMQR
D L_{q6ZK C
s#ijpc>h C
?k~(E`ZE3 C
maLKUSgo C POI :ν Poisson's ratio
2*-s3 >VK C EPSÔ :eÔ when σ'm=98kPa Reference void ratio
A('=P} I^ C ZMF :Rf=σ1/σ3 Stress ratio at critical state
k+2~=# C ZRDK :λ-κ Difference between compression and swelling index
.N8AkQ (Ok C YOUNG :E Young's modulus
D*Y4B ?, C BETA :β Parameter of tij rock subloading time-dependent model
{Q[{H'Oa C ANN :a Parameter involved in Cam-Clay subloading model &
p6]4YGw*^ C modified Cam-Clay model & tij subloading model &
o~,dkV C tij rock subloading time-dependent model
_*1/4^ C ALFA : 二次固结系数
{[PoLOCI C CNN :C 二次固结系数
w1/p wzn C SIGMAÔ :σ'mÔ Initial mean effective stress
U(DK~#} C OCR :OCR Overconsolidation ratio
@gQ{*dN C STRAIN :应变制限値 单向加载时最大的应变；交变加载时应变的振幅
DPJ#Y -0 C STRESS :应力制限值 交变加载时应力的振幅
q4<3 O"c1 C SIGMQR :蠕变应力
U[=VW0 C
Fd-PjW/E8 C LTYPE1=1: 中空扭转剪切(排水试验:应力一定;非排水试验:垂直位移一定)
NyI0 []z C LTYPE1=2: 最小主应力一定的三轴压缩
@.T
'>;izr C LTYPE1=3: 平均主应力一定的三轴压缩
[k.|iCD C LTYPE1=4: 平面应变状态侧向围压一定的压缩((SIG(1)-SIG(3)))
&a% |L=FY C LTYPE1=5: 等向压缩(SIG(1)=SIG(2)=SIG(3))
Ln0rm9FV- C
|ul25/B
B C
* VH!<k[n C NNN:单向加载时达到最大应变时计算的步数；交变加载时单向振幅的计算步数
W[[3'J TF C NTOTAL：全部计算的步数
n1 `D:XrE C NSPACE: 输出时的步数间隔
>6Y\CixN 6+(g4MW 。。。。。。
