我觉得下面这本书不错,有料,写得用心。好东西和大家分享一下:
=DwH*U/YR Q 5TyS8 《计算土力学》应该是本严谨的书,我向张教授问书中附录程序代码,张教授回复给我一个下载网址:
http://www.cm.nitech.ac.jp/cho/models/models.html。下面是其中一段代码,大家觉得还行的话,可以到上面网址去下载。
0LuY"(LR NZ9,9 Appendix IV 软岩的弹粘塑性模型的计算程序
[c )\?MWW 为了区分数字Ô和字母O,所有的数字0均用Ô表示。
R$`T"C" A-IV-1 主程序
NH6!|T IMPLICIT REAL*8(A-H,O-Z)
')~HOCBSE C
8#-}3~l[ PARAMETER (IMZ=1)
~,1X>N" C
YP97D n COMMON NN
3e1"5~?'< COMMON /TN1/TNONE(IMZ),VP(IMZ),VPÔ(IMZ)
;%9ZL[- COMMON /CRITICAL/S1S3(IMZ)
@},k\Is COMMON /TNSAND/DEPSON(6),E(3)
6
VDF@V$E COMMON /TNSAND1/LTYPE
Z6pDQ^Ii COMMON /TNSAND2/RATEÔ
40/[uW" DIMENSION DEP(6,6),SIG(6),SIGT(6),DE(6),DEPP(6),DSIG(6)
FQ^uX]<3j C
^T>.04";x OPEN(1,FILE='in.dat')
6Z:|"AwC2 OPEN(2,FILE='out.txt')
^r$5];n
C
'iK*#b8l READ(1,*)LTYPE1,NNN
l).Ijl}AH; READ(1,*)POI,EPSÔ,ZMF,ZRDK,YOUNG,BETA,AN,ALFA,CNN,OCR
izA3 INT READ(1,*)SIGMAÔ,STRAIN,RATE,SIGMQR
N?:S?p9R@ C
1-<Xi-=^{t C
AlV2tffY^ C
tJ3s#q6 C POI :ν Poisson's ratio
]P^ +~ C EPSÔ :eÔ when σ'm=98kPa Reference void ratio
"5bk82." C ZMF :Rf=σ1/σ3 Stress ratio at critical state
$R4\jIewV C ZRDK :λ-κ Difference between compression and swelling index
#xB%v C YOUNG :E Young's modulus
r&;AG@N/ C BETA :β Parameter of tij rock subloading time-dependent model
O[5ti=W C ANN :a Parameter involved in Cam-Clay subloading model &
)I_I?e C modified Cam-Clay model & tij subloading model &
A+&Va\|x C tij rock subloading time-dependent model
" OtLJ C ALFA : 二次固结系数
G<DUy^$i C CNN :C 二次固结系数
*z~Y *Q0
C SIGMAÔ :σ'mÔ Initial mean effective stress
/@e\I0P^ C OCR :OCR Overconsolidation ratio
i"@?eq#h C STRAIN :应变制限値 单向加载时最大的应变;交变加载时应变的振幅
#'oKkrl C STRESS :应力制限值 交变加载时应力的振幅
^&%?Q_] C SIGMQR :蠕变应力
J4; ".Y= C
F9" K C LTYPE1=1: 中空扭转剪切(排水试验:应力一定;非排水试验:垂直位移一定)
~XRr }z_Lq C LTYPE1=2: 最小主应力一定的三轴压缩
|xOOdy6 )~ C LTYPE1=3: 平均主应力一定的三轴压缩
@VAhmYz C LTYPE1=4: 平面应变状态侧向围压一定的压缩((SIG(1)-SIG(3)))
8*ysuL# C LTYPE1=5: 等向压缩(SIG(1)=SIG(2)=SIG(3))
R_}(p2 C
:OM>z4mQ C
3z<t# C NNN:单向加载时达到最大应变时计算的步数;交变加载时单向振幅的计算步数
h2]Od(^[ C NTOTAL:全部计算的步数
M
5h U.3.L C NSPACE: 输出时的步数间隔
=&'j;j "%Ak[04' 。。。。。。