我觉得下面这本书不错,有料,写得用心。好东西和大家分享一下:
h,"4SSL |r bWYl.b 《计算土力学》应该是本严谨的书,我向张教授问书中附录程序代码,张教授回复给我一个下载网址:
http://www.cm.nitech.ac.jp/cho/models/models.html。下面是其中一段代码,大家觉得还行的话,可以到上面网址去下载。
>3&O::]3 yv)ux:P&+ Appendix IV 软岩的弹粘塑性模型的计算程序
@O[5M2|r 为了区分数字Ô和字母O,所有的数字0均用Ô表示。
2U-3Q]/I} A-IV-1 主程序
qoAj]
") IMPLICIT REAL*8(A-H,O-Z)
|\n_OS7 C
n[(Qr9 PARAMETER (IMZ=1)
B-l'vVx C
|.wEm;Bz COMMON NN
b:w?PC~O COMMON /TN1/TNONE(IMZ),VP(IMZ),VPÔ(IMZ)
MeUaTJFEB COMMON /CRITICAL/S1S3(IMZ)
qpXWi
&g COMMON /TNSAND/DEPSON(6),E(3)
U)3DQ6T99 COMMON /TNSAND1/LTYPE
RVeEkv[qp COMMON /TNSAND2/RATEÔ
"9n3VX) DIMENSION DEP(6,6),SIG(6),SIGT(6),DE(6),DEPP(6),DSIG(6)
vv=VRhwF C
+o9":dl OPEN(1,FILE='in.dat')
85GKymz$P OPEN(2,FILE='out.txt')
gynh#&r C
2+X\}s1vN READ(1,*)LTYPE1,NNN
?@V[#. READ(1,*)POI,EPSÔ,ZMF,ZRDK,YOUNG,BETA,AN,ALFA,CNN,OCR
K#"O
a
h READ(1,*)SIGMAÔ,STRAIN,RATE,SIGMQR
x.q+uU$^ C
6&/T@LQYrh C
`hb%+-lj+ C
AFAAuFE" C POI :ν Poisson's ratio
A,4fEmWM C EPSÔ :eÔ when σ'm=98kPa Reference void ratio
[]D@Q+1 C ZMF :Rf=σ1/σ3 Stress ratio at critical state
9:M`
j C ZRDK :λ-κ Difference between compression and swelling index
Q1+dCCY#F C YOUNG :E Young's modulus
[Jv@J\ C BETA :β Parameter of tij rock subloading time-dependent model
s-v C ANN :a Parameter involved in Cam-Clay subloading model &
#f+$Ddg* C modified Cam-Clay model & tij subloading model &
f&
CBU C tij rock subloading time-dependent model
Jkek-m C ALFA : 二次固结系数
7a~X:# C CNN :C 二次固结系数
cb}[S:&| C SIGMAÔ :σ'mÔ Initial mean effective stress
,2W8=ON C OCR :OCR Overconsolidation ratio
[ 1u-Q%?# C STRAIN :应变制限値 单向加载时最大的应变;交变加载时应变的振幅
,I,\ml
C STRESS :应力制限值 交变加载时应力的振幅
|vw"[7_aS C SIGMQR :蠕变应力
^U!0-y C
6AhM=C C LTYPE1=1: 中空扭转剪切(排水试验:应力一定;非排水试验:垂直位移一定)
[W7CXZDd C LTYPE1=2: 最小主应力一定的三轴压缩
\o^2y.q:> C LTYPE1=3: 平均主应力一定的三轴压缩
&c,kQo+pA C LTYPE1=4: 平面应变状态侧向围压一定的压缩((SIG(1)-SIG(3)))
~Rr~1I&mR, C LTYPE1=5: 等向压缩(SIG(1)=SIG(2)=SIG(3))
,KFapz! C
R~6$oeWAw C
{yExQbN C NNN:单向加载时达到最大应变时计算的步数;交变加载时单向振幅的计算步数
V=*wKuB C NTOTAL:全部计算的步数
RVQh2'w C NSPACE: 输出时的步数间隔
j{Q9{}<e \7'+h5a 。。。。。。