我觉得下面这本书不错,有料,写得用心。好东西和大家分享一下:
yQNV@T<o _hu")os 《计算土力学》应该是本严谨的书,我向张教授问书中附录程序代码,张教授回复给我一个下载网址:
http://www.cm.nitech.ac.jp/cho/models/models.html。下面是其中一段代码,大家觉得还行的话,可以到上面网址去下载。
"N?+VkZEv ~;]W T Appendix IV 软岩的弹粘塑性模型的计算程序
9$qm>,o 为了区分数字Ô和字母O,所有的数字0均用Ô表示。
c!K]J A-IV-1 主程序
?9{~> 4@ IMPLICIT REAL*8(A-H,O-Z)
*Hz^K0:8( C
j|4C\~i PARAMETER (IMZ=1)
r!PpUwod C
E>|: D COMMON NN
iBTYY{-wF COMMON /TN1/TNONE(IMZ),VP(IMZ),VPÔ(IMZ)
'cix`l|^ COMMON /CRITICAL/S1S3(IMZ)
A/7{oB:a COMMON /TNSAND/DEPSON(6),E(3)
`Ag{) COMMON /TNSAND1/LTYPE
c0&Rg# COMMON /TNSAND2/RATEÔ
*)M49a*UD DIMENSION DEP(6,6),SIG(6),SIGT(6),DE(6),DEPP(6),DSIG(6)
Om;&_!i C
ar9]"s+' OPEN(1,FILE='in.dat')
RTh`ENCKR OPEN(2,FILE='out.txt')
_tTN G2 C
6Orum/|h READ(1,*)LTYPE1,NNN
o1cErI&q" READ(1,*)POI,EPSÔ,ZMF,ZRDK,YOUNG,BETA,AN,ALFA,CNN,OCR
phnV7D(E READ(1,*)SIGMAÔ,STRAIN,RATE,SIGMQR
!K
f#@0E.. C
xG&)1sT#-\ C
?L&'- e@ C
R_XR4)(< C POI :ν Poisson's ratio
}|u4 W?H C EPSÔ :eÔ when σ'm=98kPa Reference void ratio
="wzq+ U C ZMF :Rf=σ1/σ3 Stress ratio at critical state
C?Bl{4-P}* C ZRDK :λ-κ Difference between compression and swelling index
L>yJ C YOUNG :E Young's modulus
aE;le{|!({ C BETA :β Parameter of tij rock subloading time-dependent model
,I1RV C ANN :a Parameter involved in Cam-Clay subloading model &
mpPdG C modified Cam-Clay model & tij subloading model &
*t[. =_v C tij rock subloading time-dependent model
T&4qw(\G C ALFA : 二次固结系数
/='. 4v C CNN :C 二次固结系数
]vWKR."4 C SIGMAÔ :σ'mÔ Initial mean effective stress
.sC?7O= C OCR :OCR Overconsolidation ratio
ydRC1~f0 C STRAIN :应变制限値 单向加载时最大的应变;交变加载时应变的振幅
9l,a^@Y: C STRESS :应力制限值 交变加载时应力的振幅
hd),&qoW? C SIGMQR :蠕变应力
p$Ox'A4 C
1cS}J:0P C LTYPE1=1: 中空扭转剪切(排水试验:应力一定;非排水试验:垂直位移一定)
~cTN~<{dq C LTYPE1=2: 最小主应力一定的三轴压缩
+_XzmjnDd C LTYPE1=3: 平均主应力一定的三轴压缩
.Asv%p[W C LTYPE1=4: 平面应变状态侧向围压一定的压缩((SIG(1)-SIG(3)))
%LVm3e9 C LTYPE1=5: 等向压缩(SIG(1)=SIG(2)=SIG(3))
u/cg|]x&T C
MjlP+; ! C
/x_o!<M C NNN:单向加载时达到最大应变时计算的步数;交变加载时单向振幅的计算步数
V=:'SL*3| C NTOTAL:全部计算的步数
&e6CJ C NSPACE: 输出时的步数间隔
QUU;g 2k 'EDda 。。。。。。