我觉得下面这本书不错，有料，写得用心。好东西和大家分享一下：
'IX1WS&\" DzGUKJh6 《计算土力学》应该是本严谨的书，我向张教授问书中附录程序代码，张教授回复给我一个下载网址：
http://www.cm.nitech.ac.jp/cho/models/models.html。下面是其中一段代码，大家觉得还行的话，可以到上面网址去下载。
/xk7Z
q l1 Kv`v\ Appendix IV 软岩的弹粘塑性模型的计算程序
'(5GR I< 为了区分数字Ô和字母O,所有的数字0均用Ô表示。
49MEGl;K0\ A-IV-1 主程序
1Z}5ykM3 IMPLICIT REAL*8(A-H,O-Z)
.rtA sbp.! C
45c?0tj PARAMETER (IMZ=1)
j9 NF| C
b)I-do+ COMMON NN
n)kbQ] COMMON /TN1/TNONE(IMZ),VP(IMZ),VPÔ(IMZ)
xG/Q% A COMMON /CRITICAL/S1S3(IMZ)
mU #F> COMMON /TNSAND/DEPSON(6),E(3)
",b:rgpRp COMMON /TNSAND1/LTYPE
Xfc$M(a
K{ COMMON /TNSAND2/RATEÔ
r"dIB@ DIMENSION DEP(6,6),SIG(6),SIGT(6),DE(6),DEPP(6),DSIG(6)
N_c44[ z 1 C
zTg\\z; OPEN(1,FILE='in.dat')
>$Z G=& OPEN(2,FILE='out.txt')
mw
28E\U C
>B{NxL3-> READ(1,*)LTYPE1,NNN
uxf,95<g) READ(1,*)POI,EPSÔ,ZMF,ZRDK,YOUNG,BETA,AN,ALFA,CNN,OCR
.H)H9cmf READ(1,*)SIGMAÔ,STRAIN,RATE,SIGMQR
P 1mg;!tq C
\1<'XVS C
: .x((
FU C
efQ8jO C POI :ν Poisson's ratio
Rml'{S C EPSÔ :eÔ when σ'm=98kPa Reference void ratio
tY"eo Pm e C ZMF :Rf=σ1/σ3 Stress ratio at critical state
V9oBSP'kt C ZRDK :λ-κ Difference between compression and swelling index
cT'Bp)a C YOUNG :E Young's modulus
^ #:;6^Su C BETA :β Parameter of tij rock subloading time-dependent model
d)AYY}pw C ANN :a Parameter involved in Cam-Clay subloading model &
}emUpju<C C modified Cam-Clay model & tij subloading model &
wm`<+K C tij rock subloading time-dependent model
)s
$]+HQs C ALFA : 二次固结系数
(T T=i C CNN :C 二次固结系数
<< >+z5D+ C SIGMAÔ :σ'mÔ Initial mean effective stress
/w?e(v< C OCR :OCR Overconsolidation ratio
\(\a= C STRAIN :应变制限値 单向加载时最大的应变；交变加载时应变的振幅
i|^Q{3?o# C STRESS :应力制限值 交变加载时应力的振幅
.C^P6S2oJ C SIGMQR :蠕变应力
/L^g. ~ C
aoN\n]g C LTYPE1=1: 中空扭转剪切(排水试验:应力一定;非排水试验:垂直位移一定)
B#V""[Y9 C LTYPE1=2: 最小主应力一定的三轴压缩
dOeM0_o C LTYPE1=3: 平均主应力一定的三轴压缩
~71U s C LTYPE1=4: 平面应变状态侧向围压一定的压缩((SIG(1)-SIG(3)))
,sL'T[tuiU C LTYPE1=5: 等向压缩(SIG(1)=SIG(2)=SIG(3))
RAO+<m C
c< $<n C
bB_L L C NNN:单向加载时达到最大应变时计算的步数；交变加载时单向振幅的计算步数
T)mQ+&| C NTOTAL：全部计算的步数
U?0|2hR~ C NSPACE: 输出时的步数间隔
ftYJ 3/ WH u} mj )Nk 。。。。。。
