我觉得下面这本书不错,有料,写得用心。好东西和大家分享一下:
sM%l:Fv +8"P*z, 《计算土力学》应该是本严谨的书,我向张教授问书中附录程序代码,张教授回复给我一个下载网址:
http://www.cm.nitech.ac.jp/cho/models/models.html。下面是其中一段代码,大家觉得还行的话,可以到上面网址去下载。
|CQ0{1R1 0kw) -)= Appendix IV 软岩的弹粘塑性模型的计算程序
qv4r!x 为了区分数字Ô和字母O,所有的数字0均用Ô表示。
6$zd2N? A-IV-1 主程序
<AP.m4N) _ IMPLICIT REAL*8(A-H,O-Z)
,DEcCHr, C
\'|t>|zhp PARAMETER (IMZ=1)
563ExibH C
n-,mC/4 COMMON NN
KuL+~ COMMON /TN1/TNONE(IMZ),VP(IMZ),VPÔ(IMZ)
&qIdT;^=I COMMON /CRITICAL/S1S3(IMZ)
"|R75m,Id COMMON /TNSAND/DEPSON(6),E(3)
mX?t|:[b COMMON /TNSAND1/LTYPE
OI3j!L2f COMMON /TNSAND2/RATEÔ
XN{zl* ` DIMENSION DEP(6,6),SIG(6),SIGT(6),DE(6),DEPP(6),DSIG(6)
@=Dc(5`[ C
a:4!z;2
| OPEN(1,FILE='in.dat')
?ef7%0 OPEN(2,FILE='out.txt')
!DHfw-1K C
yf-2E_yB READ(1,*)LTYPE1,NNN
P^U.VXY} READ(1,*)POI,EPSÔ,ZMF,ZRDK,YOUNG,BETA,AN,ALFA,CNN,OCR
rj?c READ(1,*)SIGMAÔ,STRAIN,RATE,SIGMQR
Vock19P C
}([}A`@ C
oiD{Z C
BWB}bq C POI :ν Poisson's ratio
ml!c0< C EPSÔ :eÔ when σ'm=98kPa Reference void ratio
g`tV^b") C ZMF :Rf=σ1/σ3 Stress ratio at critical state
pRUQMPn ( C ZRDK :λ-κ Difference between compression and swelling index
s?4nR:ZC} C YOUNG :E Young's modulus
6z:/ma^
C BETA :β Parameter of tij rock subloading time-dependent model
r`RLDN!` C ANN :a Parameter involved in Cam-Clay subloading model &
->7zVAX C modified Cam-Clay model & tij subloading model &
$@L2zl1 C tij rock subloading time-dependent model
0F%?<:
& C ALFA : 二次固结系数
WMWUP ZsGS C CNN :C 二次固结系数
<Q`3;ca^ C SIGMAÔ :σ'mÔ Initial mean effective stress
fvV"H{V, C OCR :OCR Overconsolidation ratio
nKI?Sc C STRAIN :应变制限値 单向加载时最大的应变;交变加载时应变的振幅
[$:L|V!{ C STRESS :应力制限值 交变加载时应力的振幅
VZtFgN$J C SIGMQR :蠕变应力
8U7dd[ C
m'k>U4 C LTYPE1=1: 中空扭转剪切(排水试验:应力一定;非排水试验:垂直位移一定)
Lr=^0 C LTYPE1=2: 最小主应力一定的三轴压缩
;>F1?5P{ C LTYPE1=3: 平均主应力一定的三轴压缩
,}9
tJY@E C LTYPE1=4: 平面应变状态侧向围压一定的压缩((SIG(1)-SIG(3)))
Y0m?ZVt C LTYPE1=5: 等向压缩(SIG(1)=SIG(2)=SIG(3))
#B#xSmak C
yJ6g{#X4K< C
2uV5hSHYe C NNN:单向加载时达到最大应变时计算的步数;交变加载时单向振幅的计算步数
hF`<I.z} C NTOTAL:全部计算的步数
]v?jfy C NSPACE: 输出时的步数间隔
'tU \~3k :JZV=@<T 。。。。。。