我觉得下面这本书不错,有料,写得用心。好东西和大家分享一下:
Y\\nJuJo r0bPaAKw 《计算土力学》应该是本严谨的书,我向张教授问书中附录程序代码,张教授回复给我一个下载网址:
http://www.cm.nitech.ac.jp/cho/models/models.html。下面是其中一段代码,大家觉得还行的话,可以到上面网址去下载。
Rx-i.Et Z 7E)7sd Appendix IV 软岩的弹粘塑性模型的计算程序
aC 0Jfo 为了区分数字Ô和字母O,所有的数字0均用Ô表示。
>AFpO*q" A-IV-1 主程序
X6 cb#s0| IMPLICIT REAL*8(A-H,O-Z)
f`rz)C03 C
H;&t"Ql. PARAMETER (IMZ=1)
U#
B C
.w)t<7 y COMMON NN
kn|l 3+ COMMON /TN1/TNONE(IMZ),VP(IMZ),VPÔ(IMZ)
TvwIro COMMON /CRITICAL/S1S3(IMZ)
1R9/AP COMMON /TNSAND/DEPSON(6),E(3)
:!hH`l}p COMMON /TNSAND1/LTYPE
1 to<at-NN COMMON /TNSAND2/RATEÔ
>f8,YisH DIMENSION DEP(6,6),SIG(6),SIGT(6),DE(6),DEPP(6),DSIG(6)
5~0;R`D C
]zY'w,?D\F OPEN(1,FILE='in.dat')
1`\kXaG OPEN(2,FILE='out.txt')
>L4$DKO C
Mp=+*I[ READ(1,*)LTYPE1,NNN
/MtacR READ(1,*)POI,EPSÔ,ZMF,ZRDK,YOUNG,BETA,AN,ALFA,CNN,OCR
RtL'fd READ(1,*)SIGMAÔ,STRAIN,RATE,SIGMQR
[nD4\x+ C
B(dL`]@Xm C
XePBA
J C
nJg2O@mRJ C POI :ν Poisson's ratio
Tyl"N{ _ C EPSÔ :eÔ when σ'm=98kPa Reference void ratio
rM |RGe C ZMF :Rf=σ1/σ3 Stress ratio at critical state
KVy5/A/8c C ZRDK :λ-κ Difference between compression and swelling index
!,"G/}'^; C YOUNG :E Young's modulus
6<nO2 GW C BETA :β Parameter of tij rock subloading time-dependent model
axOy~%%c C ANN :a Parameter involved in Cam-Clay subloading model &
XS/TYdXB8 C modified Cam-Clay model & tij subloading model &
ir#^5e@ C tij rock subloading time-dependent model
s$6#3%h C ALFA : 二次固结系数
4gENV{L C CNN :C 二次固结系数
|_m;@.44?U C SIGMAÔ :σ'mÔ Initial mean effective stress
x0GZ2*vfsb C OCR :OCR Overconsolidation ratio
h_!"CF<n C STRAIN :应变制限値 单向加载时最大的应变;交变加载时应变的振幅
v}D! C STRESS :应力制限值 交变加载时应力的振幅
fTgN2U C SIGMQR :蠕变应力
e[!>ezaIY C
'Y Zs6rcJ C LTYPE1=1: 中空扭转剪切(排水试验:应力一定;非排水试验:垂直位移一定)
eO G%6C%a C LTYPE1=2: 最小主应力一定的三轴压缩
B35f5m7r C LTYPE1=3: 平均主应力一定的三轴压缩
AFED YRX C LTYPE1=4: 平面应变状态侧向围压一定的压缩((SIG(1)-SIG(3)))
h]P$L> C LTYPE1=5: 等向压缩(SIG(1)=SIG(2)=SIG(3))
RfRaWbn C
mX_`rvYII C
yG4MqR)J C NNN:单向加载时达到最大应变时计算的步数;交变加载时单向振幅的计算步数
c Zr4 C NTOTAL:全部计算的步数
JqZ5DjI: C NSPACE: 输出时的步数间隔
Z.JTq~`I "Fiv
]^ 。。。。。。