我觉得下面这本书不错,有料,写得用心。好东西和大家分享一下:
*km!<L7Y rUW/d3y 《计算土力学》应该是本严谨的书,我向张教授问书中附录程序代码,张教授回复给我一个下载网址:
http://www.cm.nitech.ac.jp/cho/models/models.html。下面是其中一段代码,大家觉得还行的话,可以到上面网址去下载。
"MPr'3 S] R.:T_% Appendix IV 软岩的弹粘塑性模型的计算程序
@#>YU 为了区分数字Ô和字母O,所有的数字0均用Ô表示。
[KW9J}] A-IV-1 主程序
r]B`\XWz IMPLICIT REAL*8(A-H,O-Z)
i$g6C C
c?%}J\<n PARAMETER (IMZ=1)
VD90JU]X< C
68()2v4X COMMON NN
A5IW[Gu! COMMON /TN1/TNONE(IMZ),VP(IMZ),VPÔ(IMZ)
_x]q`[Dih COMMON /CRITICAL/S1S3(IMZ)
1uJpn COMMON /TNSAND/DEPSON(6),E(3)
q71V]! COMMON /TNSAND1/LTYPE
WxJV
zHtR COMMON /TNSAND2/RATEÔ
A2%RcKY7 DIMENSION DEP(6,6),SIG(6),SIGT(6),DE(6),DEPP(6),DSIG(6)
9!,f4&G` C
Z&6*8#wn OPEN(1,FILE='in.dat')
XNr8,[c OPEN(2,FILE='out.txt')
_X~xfmU C
{O_`eS READ(1,*)LTYPE1,NNN
A(1WQUu j READ(1,*)POI,EPSÔ,ZMF,ZRDK,YOUNG,BETA,AN,ALFA,CNN,OCR
e8U6D+jY READ(1,*)SIGMAÔ,STRAIN,RATE,SIGMQR
THcX.%ToT C
Kwo0%2Onkd C
`OHdo$Y9 C
Dq<!wtFG[ C POI :ν Poisson's ratio
R,8460e7 C EPSÔ :eÔ when σ'm=98kPa Reference void ratio
3Lm7{s?=Z- C ZMF :Rf=σ1/σ3 Stress ratio at critical state
|o#pd\ C ZRDK :λ-κ Difference between compression and swelling index
~mvD|$1z C YOUNG :E Young's modulus
{q/D,Rh8 C BETA :β Parameter of tij rock subloading time-dependent model
i{}Q5iy C ANN :a Parameter involved in Cam-Clay subloading model &
?#|in} C modified Cam-Clay model & tij subloading model &
OP98 sd&T C tij rock subloading time-dependent model
,H@ x. C ALFA : 二次固结系数
*d}{7UMy# C CNN :C 二次固结系数
la_FZ C SIGMAÔ :σ'mÔ Initial mean effective stress
r}|a*dh'R C OCR :OCR Overconsolidation ratio
@ 5V3I^ C STRAIN :应变制限値 单向加载时最大的应变;交变加载时应变的振幅
o>y@1%aU C STRESS :应力制限值 交变加载时应力的振幅
W744hq@P% C SIGMQR :蠕变应力
n7,LfO# C
yC9~X='D C LTYPE1=1: 中空扭转剪切(排水试验:应力一定;非排水试验:垂直位移一定)
%5Zhq> C LTYPE1=2: 最小主应力一定的三轴压缩
R[!%d6jDE C LTYPE1=3: 平均主应力一定的三轴压缩
a9p6[qOcd C LTYPE1=4: 平面应变状态侧向围压一定的压缩((SIG(1)-SIG(3)))
Df}3^J~JX C LTYPE1=5: 等向压缩(SIG(1)=SIG(2)=SIG(3))
:Y[?@/m4 C
0:zDt~Ju C
,H5o/qNU`{ C NNN:单向加载时达到最大应变时计算的步数;交变加载时单向振幅的计算步数
(2'q~Z+>' C NTOTAL:全部计算的步数
_MzdbUb5, C NSPACE: 输出时的步数间隔
wQrD(Dv(yA `m\l#r2C 。。。。。。