我觉得下面这本书不错，有料，写得用心。好东西和大家分享一下：
%o-*~GQ@B 8eNGPuo L) 《计算土力学》应该是本严谨的书，我向张教授问书中附录程序代码，张教授回复给我一个下载网址：
http://www.cm.nitech.ac.jp/cho/models/models.html。下面是其中一段代码，大家觉得还行的话，可以到上面网址去下载。
%zGPF h!.(7qdd Appendix IV 软岩的弹粘塑性模型的计算程序
Rp#SqRy` 为了区分数字Ô和字母O,所有的数字0均用Ô表示。
{|cA[#j# A-IV-1 主程序
Dqki }k~{ IMPLICIT REAL*8(A-H,O-Z)
Tn|re Xc0e C
p\AS f PARAMETER (IMZ=1)
v|e>zm < C
-Ac^#/[ 0 COMMON NN
!**q20-aP COMMON /TN1/TNONE(IMZ),VP(IMZ),VPÔ(IMZ)
pR $c<p COMMON /CRITICAL/S1S3(IMZ)
Y+E@afsKs COMMON /TNSAND/DEPSON(6),E(3)
\hz)oC COMMON /TNSAND1/LTYPE
$[d}g COMMON /TNSAND2/RATEÔ
q:( K^ DIMENSION DEP(6,6),SIG(6),SIGT(6),DE(6),DEPP(6),DSIG(6)
eUl[ gHP C
lWR OPEN(1,FILE='in.dat')
+ x1sV *S OPEN(2,FILE='out.txt')
S}<(9@]z C
S $Wd}2> READ(1,*)LTYPE1,NNN
Q]\x O/ READ(1,*)POI,EPSÔ,ZMF,ZRDK,YOUNG,BETA,AN,ALFA,CNN,OCR
.s+e
hZ READ(1,*)SIGMAÔ,STRAIN,RATE,SIGMQR
GC~::m~ C
K vgZx(. C
h W-[omr0 C
u[% #/ C POI :ν Poisson's ratio
P VPwYmte C EPSÔ :eÔ when σ'm=98kPa Reference void ratio
j2z$kw% C ZMF :Rf=σ1/σ3 Stress ratio at critical state
m~v
Ie c C ZRDK :λ-κ Difference between compression and swelling index
wBf
bpoE7 C YOUNG :E Young's modulus
EpiagCS C BETA :β Parameter of tij rock subloading time-dependent model
|67UN U C ANN :a Parameter involved in Cam-Clay subloading model &
xnA rYm C modified Cam-Clay model & tij subloading model &
*m 7e>]- C tij rock subloading time-dependent model
U!h!z`RU54 C ALFA : 二次固结系数
ZISR]xay C CNN :C 二次固结系数
5g=" # C SIGMAÔ :σ'mÔ Initial mean effective stress
,F1$Of/'@\ C OCR :OCR Overconsolidation ratio
uQDu<@5^[ C STRAIN :应变制限値 单向加载时最大的应变；交变加载时应变的振幅
(L\tp>
E- C STRESS :应力制限值 交变加载时应力的振幅
NJ~'`{3v C SIGMQR :蠕变应力
D4G{= Y}G C
uBM1;9h C LTYPE1=1: 中空扭转剪切(排水试验:应力一定;非排水试验:垂直位移一定)
$ 8s&=OW C LTYPE1=2: 最小主应力一定的三轴压缩
wG B'c's* C LTYPE1=3: 平均主应力一定的三轴压缩
oq|K:<l C LTYPE1=4: 平面应变状态侧向围压一定的压缩((SIG(1)-SIG(3)))
+CACs7tV C LTYPE1=5: 等向压缩(SIG(1)=SIG(2)=SIG(3))
@[^H*^1|g C
,i}"e(f C
W{%M+a[#l C NNN:单向加载时达到最大应变时计算的步数；交变加载时单向振幅的计算步数
\ %K6T)9 C NTOTAL：全部计算的步数
0
[s1!Cm!i C NSPACE: 输出时的步数间隔
9X-DR (HEjmQjE 。。。。。。
