我觉得下面这本书不错,有料,写得用心。好东西和大家分享一下:
dljE.peL Q Be6\oq 《计算土力学》应该是本严谨的书,我向张教授问书中附录程序代码,张教授回复给我一个下载网址:
http://www.cm.nitech.ac.jp/cho/models/models.html。下面是其中一段代码,大家觉得还行的话,可以到上面网址去下载。
380` >"D HBlk~eZ Appendix IV 软岩的弹粘塑性模型的计算程序
50,'z?-_ 为了区分数字Ô和字母O,所有的数字0均用Ô表示。
D|-]<r1" A-IV-1 主程序
n6/Ous IMPLICIT REAL*8(A-H,O-Z)
(Ou%0
KW C
;Shu PARAMETER (IMZ=1)
@- U\!Tf C
$%bSRvA COMMON NN
F-&tSU, COMMON /TN1/TNONE(IMZ),VP(IMZ),VPÔ(IMZ)
T[oC='I+O COMMON /CRITICAL/S1S3(IMZ)
pSs*Z6c)@ COMMON /TNSAND/DEPSON(6),E(3)
r/Y J, 2! COMMON /TNSAND1/LTYPE
USg"wJY COMMON /TNSAND2/RATEÔ
C/kf?:j DIMENSION DEP(6,6),SIG(6),SIGT(6),DE(6),DEPP(6),DSIG(6)
3BFOZV+ C
e;]tO-Nu OPEN(1,FILE='in.dat')
[9m3@Yd' OPEN(2,FILE='out.txt')
AGlBvRX7e C
G@]3EP READ(1,*)LTYPE1,NNN
^HKXm#vAB READ(1,*)POI,EPSÔ,ZMF,ZRDK,YOUNG,BETA,AN,ALFA,CNN,OCR
.wfydu)3 READ(1,*)SIGMAÔ,STRAIN,RATE,SIGMQR
CMt<oT6.? C
$O"ss>8Se C
%yR XOt2( C
8q2a8I9g C POI :ν Poisson's ratio
mQ"~x] C EPSÔ :eÔ when σ'm=98kPa Reference void ratio
HW@wia C ZMF :Rf=σ1/σ3 Stress ratio at critical state
eg0_ < C ZRDK :λ-κ Difference between compression and swelling index
Iy<>-e"| C YOUNG :E Young's modulus
NR4+&d C BETA :β Parameter of tij rock subloading time-dependent model
w,UE0i9I C ANN :a Parameter involved in Cam-Clay subloading model &
J4Gzp~{ C modified Cam-Clay model & tij subloading model &
Q6h+. C tij rock subloading time-dependent model
HBp??.r C ALFA : 二次固结系数
x9F* $G C CNN :C 二次固结系数
n}Z%-w$K# C SIGMAÔ :σ'mÔ Initial mean effective stress
Vb?wwx7= C OCR :OCR Overconsolidation ratio
dXDyY C STRAIN :应变制限値 单向加载时最大的应变;交变加载时应变的振幅
K_+M?ap_ C STRESS :应力制限值 交变加载时应力的振幅
9 5!xJdq C SIGMQR :蠕变应力
f=91
Z_M C
OF*E1BM C LTYPE1=1: 中空扭转剪切(排水试验:应力一定;非排水试验:垂直位移一定)
P6%qNR/ x C LTYPE1=2: 最小主应力一定的三轴压缩
)F hbN@3 C LTYPE1=3: 平均主应力一定的三轴压缩
pImq<Z C LTYPE1=4: 平面应变状态侧向围压一定的压缩((SIG(1)-SIG(3)))
<&[`
+ C LTYPE1=5: 等向压缩(SIG(1)=SIG(2)=SIG(3))
Lf9s'o}.R C
IsT}T}p,t C
.~I:Hcf/ C NNN:单向加载时达到最大应变时计算的步数;交变加载时单向振幅的计算步数
=eY C NTOTAL:全部计算的步数
r } Wdj C NSPACE: 输出时的步数间隔
`}t5` :#k z@UH[>^gj 。。。。。。