我觉得下面这本书不错,有料,写得用心。好东西和大家分享一下:
{r&mNbz 2{rWAPHgz 《计算土力学》应该是本严谨的书,我向张教授问书中附录程序代码,张教授回复给我一个下载网址:
http://www.cm.nitech.ac.jp/cho/models/models.html。下面是其中一段代码,大家觉得还行的话,可以到上面网址去下载。
vPNbV Fvl\. Appendix IV 软岩的弹粘塑性模型的计算程序
Cdz&'en^ 为了区分数字Ô和字母O,所有的数字0均用Ô表示。
0!\C@wnH A-IV-1 主程序
_-_iw&F IMPLICIT REAL*8(A-H,O-Z)
9&eY<'MgP C
Xxg|01 PARAMETER (IMZ=1)
/grTOf& C
3U@jw,K!{A COMMON NN
j~-N2b6z COMMON /TN1/TNONE(IMZ),VP(IMZ),VPÔ(IMZ)
L3, /7 COMMON /CRITICAL/S1S3(IMZ)
#Q8_:dPY COMMON /TNSAND/DEPSON(6),E(3)
-qIi.]/f"9 COMMON /TNSAND1/LTYPE
`MOw\Z).. COMMON /TNSAND2/RATEÔ
_`udd)Y2 DIMENSION DEP(6,6),SIG(6),SIGT(6),DE(6),DEPP(6),DSIG(6)
fs'SCwx C
!cyrt< OPEN(1,FILE='in.dat')
1!v{#w{u7 OPEN(2,FILE='out.txt')
0Qt!w( C
odhS0+d^ READ(1,*)LTYPE1,NNN
PEf yHf7` READ(1,*)POI,EPSÔ,ZMF,ZRDK,YOUNG,BETA,AN,ALFA,CNN,OCR
, _e[P READ(1,*)SIGMAÔ,STRAIN,RATE,SIGMQR
YhN<vZ}U!~ C
|H?t+Dyn)q C
t eY@)F C
f*o C POI :ν Poisson's ratio
"JB4Uaa C EPSÔ :eÔ when σ'm=98kPa Reference void ratio
QRQZ{m C ZMF :Rf=σ1/σ3 Stress ratio at critical state
*:n~j9V- C ZRDK :λ-κ Difference between compression and swelling index
Z3S+")^ C YOUNG :E Young's modulus
xuK"pS C BETA :β Parameter of tij rock subloading time-dependent model
NU=ru/ C ANN :a Parameter involved in Cam-Clay subloading model &
[:y:_ECs6 C modified Cam-Clay model & tij subloading model &
[CJ<$R ! C tij rock subloading time-dependent model
JsJP%'^/R C ALFA : 二次固结系数
:0J`4 C CNN :C 二次固结系数
PK[mf\G\ C SIGMAÔ :σ'mÔ Initial mean effective stress
f% pT-# C OCR :OCR Overconsolidation ratio
|.N[NY C STRAIN :应变制限値 单向加载时最大的应变;交变加载时应变的振幅
k>($[;k|b C STRESS :应力制限值 交变加载时应力的振幅
+S4n416K C SIGMQR :蠕变应力
V]zc-gYI C
DDeU: C LTYPE1=1: 中空扭转剪切(排水试验:应力一定;非排水试验:垂直位移一定)
\d@5*q C LTYPE1=2: 最小主应力一定的三轴压缩
F)^0R%{C C LTYPE1=3: 平均主应力一定的三轴压缩
6'# 5Dqw"r C LTYPE1=4: 平面应变状态侧向围压一定的压缩((SIG(1)-SIG(3)))
;st0Ekni) C LTYPE1=5: 等向压缩(SIG(1)=SIG(2)=SIG(3))
b>uD-CSA C
B>53+GyMV C
m+(Cl#+ C NNN:单向加载时达到最大应变时计算的步数;交变加载时单向振幅的计算步数
=)Xj[NNRT C NTOTAL:全部计算的步数
xo[o^go C NSPACE: 输出时的步数间隔
7ch9Pf n28JWkK8 。。。。。。