我觉得下面这本书不错，有料，写得用心。好东西和大家分享一下：
A4uDuB;;ZQ 0NWtu]9QC 《计算土力学》应该是本严谨的书，我向张教授问书中附录程序代码，张教授回复给我一个下载网址：
http://www.cm.nitech.ac.jp/cho/models/models.html。下面是其中一段代码，大家觉得还行的话，可以到上面网址去下载。
9ApGn!` dd%-bI^ Appendix IV 软岩的弹粘塑性模型的计算程序
` m Tc 为了区分数字Ô和字母O,所有的数字0均用Ô表示。
#T8$NZA A-IV-1 主程序
r=ds'n" IMPLICIT REAL*8(A-H,O-Z)
4$!iw3N ( C
M.R] hI PARAMETER (IMZ=1)
ec` $2 u C
N%&D(_ COMMON NN
aF\?X &| COMMON /TN1/TNONE(IMZ),VP(IMZ),VPÔ(IMZ)
)C C rO COMMON /CRITICAL/S1S3(IMZ)
W e*)RXm % COMMON /TNSAND/DEPSON(6),E(3)
HR83{B21 COMMON /TNSAND1/LTYPE
x2g P, p- COMMON /TNSAND2/RATEÔ
y:+s*x6Vg DIMENSION DEP(6,6),SIG(6),SIGT(6),DE(6),DEPP(6),DSIG(6)
ZM%z"hO9R C
s%R'c_cGZ OPEN(1,FILE='in.dat')
,0Y5O?pu\ OPEN(2,FILE='out.txt')
Y4|g^>{<ni C
nob^
I 5? READ(1,*)LTYPE1,NNN
qP0_#l& READ(1,*)POI,EPSÔ,ZMF,ZRDK,YOUNG,BETA,AN,ALFA,CNN,OCR
[ ,fd Nxc8 READ(1,*)SIGMAÔ,STRAIN,RATE,SIGMQR
*! y.!v* C
&$</|F)y C
lhA<wV1-9G C
ku}I; k | C POI :ν Poisson's ratio
\3@A C7 C EPSÔ :eÔ when σ'm=98kPa Reference void ratio
l6Q75i)eF C ZMF :Rf=σ1/σ3 Stress ratio at critical state
|+MV%QG; C ZRDK :λ-κ Difference between compression and swelling index
\aO.LwYm;: C YOUNG :E Young's modulus
P>t[35 /1 C BETA :β Parameter of tij rock subloading time-dependent model
{1ic* cZS C ANN :a Parameter involved in Cam-Clay subloading model &
U )N_ / C modified Cam-Clay model & tij subloading model &
+v tI1LC;_ C tij rock subloading time-dependent model
|E)IJj
3 C ALFA : 二次固结系数
A{Q :,S) C CNN :C 二次固结系数
2 <@27 C5 C SIGMAÔ :σ'mÔ Initial mean effective stress
+t XOP| X C OCR :OCR Overconsolidation ratio
r?= 7#/] C STRAIN :应变制限値 单向加载时最大的应变；交变加载时应变的振幅
/RmHG
H! C STRESS :应力制限值 交变加载时应力的振幅
U'acVcD C SIGMQR :蠕变应力
?I`ru:iG C
rpT{0 >5 C LTYPE1=1: 中空扭转剪切(排水试验:应力一定;非排水试验:垂直位移一定)
_( 'KNA~ C LTYPE1=2: 最小主应力一定的三轴压缩
UMJ>6 Ko8 C LTYPE1=3: 平均主应力一定的三轴压缩
9 v<Sng C LTYPE1=4: 平面应变状态侧向围压一定的压缩((SIG(1)-SIG(3)))
j\ SW~}d9 C LTYPE1=5: 等向压缩(SIG(1)=SIG(2)=SIG(3))
| <ZkJR3B C
cAE.I$T( C
!*gTC1bv B C NNN:单向加载时达到最大应变时计算的步数；交变加载时单向振幅的计算步数
:J}@*>c C NTOTAL：全部计算的步数
e
r;3TG~ C NSPACE: 输出时的步数间隔
8HLcDS# \o j#*aL^ 。。。。。。
