我觉得下面这本书不错,有料,写得用心。好东西和大家分享一下:
= U[$i"+ O]61guxro 《计算土力学》应该是本严谨的书,我向张教授问书中附录程序代码,张教授回复给我一个下载网址:
http://www.cm.nitech.ac.jp/cho/models/models.html。下面是其中一段代码,大家觉得还行的话,可以到上面网址去下载。
' hdLQ\J z$M-UxY Appendix IV 软岩的弹粘塑性模型的计算程序
VAe[x
` 为了区分数字Ô和字母O,所有的数字0均用Ô表示。
ZQ-6n1O A-IV-1 主程序
1Cw$^jd IMPLICIT REAL*8(A-H,O-Z)
)d|s$l$?7 C
=CjWPZShV PARAMETER (IMZ=1)
k$kxw_N5d C
oE0~F|(\1 COMMON NN
O]@s`w COMMON /TN1/TNONE(IMZ),VP(IMZ),VPÔ(IMZ)
,v?FR
}v COMMON /CRITICAL/S1S3(IMZ)
SQMtR2 COMMON /TNSAND/DEPSON(6),E(3)
G6xNR COMMON /TNSAND1/LTYPE
+Z]}ce
u" COMMON /TNSAND2/RATEÔ
KJs`[,;< DIMENSION DEP(6,6),SIG(6),SIGT(6),DE(6),DEPP(6),DSIG(6)
QcjsQTAbk C
!HeQMz OPEN(1,FILE='in.dat')
UVW4KUxR OPEN(2,FILE='out.txt')
j'lC]}kH C
wbg?IvY[ READ(1,*)LTYPE1,NNN
sN) xNz READ(1,*)POI,EPSÔ,ZMF,ZRDK,YOUNG,BETA,AN,ALFA,CNN,OCR
[lu+"V,<LJ READ(1,*)SIGMAÔ,STRAIN,RATE,SIGMQR
aadw#90 C
qt4%=E;[ C
@CoUFdbz C
*<}R=X. C POI :ν Poisson's ratio
aEun *V^, C EPSÔ :eÔ when σ'm=98kPa Reference void ratio
]K%d C ZMF :Rf=σ1/σ3 Stress ratio at critical state
&MB1'~Q,hq C ZRDK :λ-κ Difference between compression and swelling index
Z&7Yl(| C YOUNG :E Young's modulus
VA %lJ!$ C BETA :β Parameter of tij rock subloading time-dependent model
5'<mfY'B C ANN :a Parameter involved in Cam-Clay subloading model &
5'<a,,RKu C modified Cam-Clay model & tij subloading model &
mXyg\5 C tij rock subloading time-dependent model
j9-.bGtm?. C ALFA : 二次固结系数
u<+"#.[2v~ C CNN :C 二次固结系数
IIF <Zkpb C SIGMAÔ :σ'mÔ Initial mean effective stress
h<n 2pz} C OCR :OCR Overconsolidation ratio
y'z9Ya C STRAIN :应变制限値 单向加载时最大的应变;交变加载时应变的振幅
EiyHZ C STRESS :应力制限值 交变加载时应力的振幅
90ov[|MkM C SIGMQR :蠕变应力
zy|h1.gd C
hZ')<@hNP C LTYPE1=1: 中空扭转剪切(排水试验:应力一定;非排水试验:垂直位移一定)
x1"8K C LTYPE1=2: 最小主应力一定的三轴压缩
nj-LG!"a C LTYPE1=3: 平均主应力一定的三轴压缩
jLw|F-v-l< C LTYPE1=4: 平面应变状态侧向围压一定的压缩((SIG(1)-SIG(3)))
G'wyH[ d/ C LTYPE1=5: 等向压缩(SIG(1)=SIG(2)=SIG(3))
jHV)
TBr C
"s@q(J C
b=SCyGxlZ5 C NNN:单向加载时达到最大应变时计算的步数;交变加载时单向振幅的计算步数
^S(QvoaQ C NTOTAL:全部计算的步数
mApl;D X C NSPACE: 输出时的步数间隔
AXv-%k}; R?Ch8mW.! 。。。。。。