我觉得下面这本书不错,有料,写得用心。好东西和大家分享一下:
ZD`0(CkXb >`Zw0S 《计算土力学》应该是本严谨的书,我向张教授问书中附录程序代码,张教授回复给我一个下载网址:
http://www.cm.nitech.ac.jp/cho/models/models.html。下面是其中一段代码,大家觉得还行的话,可以到上面网址去下载。
'645Fr[lg DzG$\%G2R} Appendix IV 软岩的弹粘塑性模型的计算程序
)W$@phY(I 为了区分数字Ô和字母O,所有的数字0均用Ô表示。
./E<v A-IV-1 主程序
5jK9cF$> IMPLICIT REAL*8(A-H,O-Z)
1A;f[Rze C
C!S(!Z, PARAMETER (IMZ=1)
5vqh09-FB C
[[VB'Rs COMMON NN
y,vrMWDy COMMON /TN1/TNONE(IMZ),VP(IMZ),VPÔ(IMZ)
.
I#dR* COMMON /CRITICAL/S1S3(IMZ)
PitDk
1T COMMON /TNSAND/DEPSON(6),E(3)
hYU4%"X COMMON /TNSAND1/LTYPE
jij-pDQnv COMMON /TNSAND2/RATEÔ
Vh5Z'4N DIMENSION DEP(6,6),SIG(6),SIGT(6),DE(6),DEPP(6),DSIG(6)
}C9P-- C
iP+3) OPEN(1,FILE='in.dat')
_t;Mi/\P OPEN(2,FILE='out.txt')
PvqG5-L~W C
-4nSiI READ(1,*)LTYPE1,NNN
137:T: READ(1,*)POI,EPSÔ,ZMF,ZRDK,YOUNG,BETA,AN,ALFA,CNN,OCR
R}Z2rbt READ(1,*)SIGMAÔ,STRAIN,RATE,SIGMQR
y?yWM8 C
,lr\XhO C
tAH0o\1; C
1JSKK.LuJV C POI :ν Poisson's ratio
ok<!/"RX$ C EPSÔ :eÔ when σ'm=98kPa Reference void ratio
!O*uQB C ZMF :Rf=σ1/σ3 Stress ratio at critical state
/;a b"b C ZRDK :λ-κ Difference between compression and swelling index
\N"=qw^ t C YOUNG :E Young's modulus
o<nkK+=Afm C BETA :β Parameter of tij rock subloading time-dependent model
t{RdqAF C ANN :a Parameter involved in Cam-Clay subloading model &
m #}%l3$ C modified Cam-Clay model & tij subloading model &
= s$UU15 C tij rock subloading time-dependent model
x#,nR]C C ALFA : 二次固结系数
-|nHwSrCZ/ C CNN :C 二次固结系数
/DN!" C SIGMAÔ :σ'mÔ Initial mean effective stress
x4nmDEpa C OCR :OCR Overconsolidation ratio
-/ h'uG C STRAIN :应变制限値 单向加载时最大的应变;交变加载时应变的振幅
'r_NA!R C STRESS :应力制限值 交变加载时应力的振幅
bO\E)%zp C SIGMQR :蠕变应力
e!JC5Al7 C
.0f6b C LTYPE1=1: 中空扭转剪切(排水试验:应力一定;非排水试验:垂直位移一定)
|}^BF%8V: C LTYPE1=2: 最小主应力一定的三轴压缩
"fX9bh^ C LTYPE1=3: 平均主应力一定的三轴压缩
$ @Fvl-lK C LTYPE1=4: 平面应变状态侧向围压一定的压缩((SIG(1)-SIG(3)))
Y; OqdO C LTYPE1=5: 等向压缩(SIG(1)=SIG(2)=SIG(3))
B gB]M3Il C
x;p7n2_ C
>lIQM3 C NNN:单向加载时达到最大应变时计算的步数;交变加载时单向振幅的计算步数
S,%BhQ[ C NTOTAL:全部计算的步数
ohj(1jt C NSPACE: 输出时的步数间隔
RbGq$vYol/ &$+nuUA 。。。。。。