岩土工程数值分析学习笔记.docx
《岩土工程数值分析学习笔记.docx》由会员分享，可在线阅读，更多相关《岩土工程数值分析学习笔记.docx（26页珍藏版）》请在冰豆网上搜索。
岩土工程数值分析学习笔记
岩土工程数值分析读书笔记
理论部分
0岩土工程数值分析简介
岩土工程问题解析分析是以弹塑性力学理论和结构力学作为理论依据，适用于解决连续介质、各向同性材料、未知量少、边界条件简单的工程问题，存在很大的局限性。
岩土工程问题数值分析是借助于计算机的计算能力，适用于解决材料复杂、边界条件复杂、任意荷载、任意几何形状，适用范围广。
岩土工程数值分析发展过程：
20世纪40年代，使用差分法解决了土工中的渗流及固结问题，如土坝渗流及浸润线的求法、土坝及地基的固结等。
20世纪60年代，使用有限元法成解决了土石坝的静力问题的求解。
20世纪70年代，使用有限元法解决了土石坝及高楼（包括地基）的抗震分析。
20世纪80年代，边界元法异军突起，解决了半无限域的边界问题；地基的静力及动力问题都使用边界元法得到了有效地解决。
岩土工程数值分析的方法有两类，一类方法是将土视为连续介质，随后又将其离散化，如有限单元法、有限差分法、边界单元法、有限元线法、无单元法以及各种方法的耦合。
另一类计算方法是考虑岩土材料本身的不连续性，如裂缝及不同材料间界面的界面模型和界面单元的使用，离散元法，不连续变形分析，流形元法，颗粒流等数值计算方法。
1数值分析过程中存在的问题及解决措施
问题：
（1）对岩土工程数值分析方法缺乏系统的知识和深入的理解，出现问题时不知道在什么情况下属于理论问题或数学模型问题；在什么情况下是属于计算方法问题或本构模型问题；在什么情况下是参数的确定问题或计算本身的问题等。
（2）各种本构模型固有的局限性。
具有多相性土的物理力学性质太复杂，难以准确地用数学模型和本构模型描述。
例如邓肯一张模型不能反映剪胀性，不能反映压缩与剪切的交叉影响；
（3）现有的试验手段和设备不能提供适当、合理和精确的参数。
靠少数样本点所获得的参数难以准确地描述整个空间场地的物理力学性能；土的参数因土样扰动难以高质量的获取，其精度很差。
（4）数学模型还会给人造成一种错觉，让人觉得其计算结果也一定会更好、更可靠。
这样可能
使人们忽略了精确的数学公式也照样会有出错的可能性。
只有当输入参数的质量和精度很高,
并能与数学模型的精度相匹配时，才有可能得到较为准确的计算结果。
措施：
（1）加强对土的本构模型的教学与培训，了解和掌握各种土的本构模型的优点和局限性以及模型参数的离散性。
（2）在使用数值分析方法的同时，不断地积累使用经验，包括他人的经验。
（3）在丰富的使用经验基础上，建立相应的使用规范。
2有限元法软件plaxis理论学习
（1）模型的选用和分析:
Mohr-Coulomb模型（MC）:
这种模型被推荐用于问题的初步分析，对于每个土层，可以估计出一个平均刚度常数，由于这个刚度是常数，计算往往会相对较快，可以得到变形的一个初步印象。
五个输入参数，即：
表示土体弹性的E和V表示土体塑性的？
和C,以及剪
胀角％
节理岩石模型（JR）:
节理模型是一种各向异性的弹塑性模型，特别适用于模拟包括层理尤其是断层方向在内的岩层行为等。
两个参数？
和Co
Hardening-Soil模型（HS）:
是一种改进了的模拟岩土行为的模型，适合所有类型的土；不能用来模拟滞后或者反复循环加载情形。
三个参数三轴加载刚度E50、三轴卸载刚度Eur和固
结仪加载刚度Eoedo我们一般取Eur=3E50和Eoed=E50作为不同土体类型的平均值。
软土蠕变模型（SSC）:
软土蠕变模型是一个新近开发的应用于地基和路基等的沉陷问题的模型，通常会过高地预计弹性岩土的行为范围。
软土模型（SS）:
软土模型是一种Cam-Clay类型的模型，特别适用于接近正常固结的粘性土的主压缩，可以被Hardening-Soil模型取代。
改进的Cam-Clay模型（MCC）:
主要用于模拟接近正常固结的粘性土，在实际应用中是不被推荐的。
如果要对所考虑的问题进行一个简单迅速的初步分析，我们建议使用Mohr-Coulomb模
型。
当缺乏好的土工数据时，进一步的高级分析是没有用的。
在许多情况下，当你拥有主导土层的好的数据时，可以利用Hardening-Soil模型
来进行一个额外的分析。
毫无疑问，同时拥有三轴试验和固结仪试验结果的可能性是很小的。
但是，原位实验数据的修正值对高质量实验数据来说是一个有益的补充最后，软土蠕变模型可以用于分析蠕变（即：
极软土的次压缩）。
用不同的土工模型来分析同一个岩土问题显得代价过高，但是它们往往是值得的。
首先，用Mohr-Coulomb模型来分析是相对
较快而且简单的；其次，这一过程通常会减小计算结果的误差。
（2）各模型的参数：
MOHR-COULOMB模型的基本参数：
E:
'杨氏模量；v:
'白桑比；？
：
内摩擦角；c：
内聚力；収剪胀角，这些参数可以从土样的基本试验得到。
下一畀他]alAuu]ih/g
杨氏模量E:
在土力学中，初始斜率用E0表示，50%强度处的割线模量由E50表示。
对于具有大范围线弹性行为的材料来说，使Eo是符合实际的，但是对于土体加载问题一般使用E50。
如果考虑隧道和开挖问题中的卸载问题要用Eur替换E50.
泊桑比V：
在许多情况下得到的V值是介于0.3和0.4之间的。
一般地说，除了一维压缩，这个范围的值还可以用在加载条件下，
在卸载条件下，使用0.15和0.25之间的值更为普遍。
内摩擦角:
PLAXIS可以处理无粘性砂土c=0,但是有一些选项的执行不太好。
为了避免复杂性，我们建议不熟练的用户至少输入一个较小值（使用c>0.2kPa）。
ijj的小的负值仅仅对极
剪胀角：
以度的方式指定的。
除了严重的超固结土层以外，粘性土通常没有什么剪胀性（书=0。
砂土的剪胀性依赖于密度和摩擦角。
对于石英砂土来说，书=?
-30，书的值比？
的值小30度，然而，剪胀角在多数情况下为零。
松的砂土是实际的。
HARDENING-SOIL模型的参数:
软土模型参数:
/T：
修缩指标/：
修1E的膨胀指标
C:
内聚力
问：
（P：
内摩擦角妙：
剪胀角禹级参数（使用缺省哎世）：
%一:
卸载/晅新加载比
Kf：
If：
常固结中侧应力系数
下一凱皿
（3）单位和符号规定
为方便起见.以卜-列出经常便用的味套单位£
模力角移n喪掠
标卅常他
[m]fkNl宙心]
[m]
[e]fkPa]=[kNfkPal
「dg氏]
fklN
『sec]
曲]
fkipsl^fklb/sqft]fkipsl
fdeM-l_
2T
[gj
[kNill-J
[klKcunj
渗选系散
（tiVdayJ
[fl'secJ
[Z
Lkibl
绒荷磁
[kNmJ
[klb/nj
分衍荷戎
[kPaJ
L口戶］
应丿J
[kPii]
[Ip叩
az指向平面外。
轴对称分析里，x代表径向坐标，y代表示径向应力，az表示环向应力。
在所有输出数据里，
符号规定：
在平面应变分析里，表轴向坐标，Z代表切向。
此时，压应力（包括孔隙压力）和压力设为负值，而拉应力和拉力设为正值。
CXx
I十#一口;卜01屯血王上打母a・・门L*丄證…,W"r
rate【血」
<1
警
■家建模过程：
绘制几何轮廓线f设定边界条件+添加荷载f添加材料特性f划分网格+设置初始条件+计算工序分层f计算
警告：
如果修改已经建好的模型，要意识到必须重新生成网格，有时还包括重新生成初始条件,以便和修改后的模型一致。
（4）荷载问题：
荷载子菜单包含的选项，用于在几何模型里引入分布荷载、线荷载或集中荷载以及指定位移。
荷载和指定位移既可用于模型边界，也可用于模型内部。
指定位移是为了控制某些点的位移而强加于模型的特殊条件。
大小为零的指定位移即为约束。
注意：
如果在一条几何线上同时施加指定位移和荷载，那么在计算过程中指定位移优先于荷
载，在全约束线上施加指定位移，约束则优先于指定位移。
因而，在全约束线段上没有必要施加指定位移、集中荷载、分布荷载。
如果只有一个位移方向是预加的，而其他方向是自由的，那么就可以在自由位移的方向上施加集中荷载、分布荷载。
（5）材料性状种类一材料类型
排水性状：
该设置不产生超静水压。
很明显，它适用于干土，以及由于高渗透性（砂土）和/或缓慢加载的完全排水条件。
该选项也可以用来模拟长期岩土性质，其中不需要模拟不排水加载和固结的精确历史。
不排水性状：
该设置用于研究超静水压的完全发展过程。
由于低透水性（粘土）和/或快速加载，孔隙水渗流有时可以忽略。
注意，要输入有效的模型参数。
非多孔性状：
使用该设置的类组，不论是初始孔压还是超静水压，都不予考虑。
有关应用可以在模拟混凝土或结构性状的例子里找到。
非多孔性状通常和线弹性模型组合应用。
输入的饱和容重和渗透性参数对非多孔材料是无效的。
（6）网格生成
网格的基本单元类型为15节点的三角形单元或6节点的三角形单元。
全局疏密度水平：
很粗疏、
为粗疏。
很粗疏：
/J,=25
大约50个中”元
/?
.~50
大约100个单元
n.=100
大约250个单兀
细新
/?
-200
人约500个单元
/7,-100
人约L000个单-元
局部疏密度：
在可能出现强烈的应力集中或大变形梯度的区域，最好是用更精确（细密）的有限元元网格来模拟，这通常是在几何模型里出现边缘、角或结构对象的时候用到。
这时候的局部单元尺寸因子就应当等于0.5。
关于网格生成的几点建议：
为了优化有限元计算，可以先用相对粗疏的网格作初始分析。
这样可以检查几何模型的大小是否合适，还可以预计发生应力集中和大变形梯度的位置。
所得成果主要用来加密有限元网格，生成加密的有限元模型。
为了能够顺利地生成详细的有限元网格，首先应当从网格子菜单里选择需要的全局疏密度。
此外，如果还要进行局部加密，那么应当首先加密类组，然后加密几何线，最后再加密几何点。
需要时，可以直接给几何点设局部单元尺寸因子。
（7）初始条件
在几何模型和有限元网格都建好之后，必须明确初始应力状态和初始构造。
初始条件由两个不同的模式组成：
一个模式是生成初始水压（水力条件模式），另一个模式是定义初始
几何构造和生成初始有效应力场（几何构造模式）。
注意：
从初始条件返回到几何图形创建模式是允许的，但是一般不这么做，因为这样做可
能会丢失一些有关初始条件的信息。
（8）计算
计算分类：
分三种基本计算类型，塑性计算、固结分析和Phi-c折减（安全分析）。
塑性计算：
在弹-塑性变形分析当中，如果不需要考虑超静水压随时间的减小，那么应当选塑性计算。
一般性塑性计算根据未变形的几何图形得出刚度矩阵。
这类计算适用于大多数岩土工程的实际应用情况。
固结分析：
要分析饱和粘性土的超静水压随时间的变化与消散，应该选择固结计算。
一般情况下，执行完不排水塑性计算之后，再执行不添加另外荷载的固结分析。
也可以在固结分析过程中施加荷载。
但是，在接近破坏时要引起注意，因为此时迭代过程可能会不收敛。
Phi-c折减（安全分析）：
可以用减小土的强度参数的方法来执行安全分析。
这个过程就叫做Phi-c折减，是一个单独的计算类型。
如果要计算全局安全系数，就应当选择Phi-c折减。
建议在计算工序列表的末尾再定义所有安全分析，并且用起始工序号参数注明是在给哪个计算工序求算安全系数。
plaxis例题模拟
算例一砂土层上圆形基础的沉降（选用柔性基础）
詈QE
loacJ
E：
I:
i11:
1-11t-L-ZZ-Z二！
;Z-iZ!
：
--J：
.二ai■二i±»zI上■二-i-z-ri-z-t■r和二；■二二'土二・二z-:
二:
r
■■■
--.■--.----~~■~=~J-Ir
■■■■*■=-・■
--I--II--I1-1a-J-I■■-I■■-I-■J-■J-,.™-I--I-■:
-F;■-_E■--.■
rm『门rmnTrnTTmTn
團1砂土丘上阖畑：
《«11曲？
5^睁
图1表示放置在4m厚砂土层上半径为1m的一个圆形基础。
砂土层下是深厚的坚硬岩石
层。
计算土体在上部荷载作用下产生的位移和应力。
选用轴对称模型，竖直方向取4m（到岩
石层顶部），水平方向取5m（5倍的基础半径）。
计算过程如下：
—
CpPTTfllwrings.
K£>«1Hrtn
图5新建项目参数设置
"PtAiiiaJC^Lubi^-l^ficLdtVitu~fduiLik~~Het
sag3»也丽
iIErm（PTI1F'rm|
Jw™
ftrtw.'D.：
[t~[拂a%画
I：
•曲
■p■八
hn
|pfcrtc
4吐n14ccb±..1
ZdorricKUi1nu;cna.Il-LdoAtnn
1LoafkuFVXt
1TtK11F
1tU[^wQD呼
OjO]DC
1DT
JL^eccDTCtudtn
ajM=.D
「F.
1
Pl^tepnOfWrti«
*letsfialset
[denhfkaton:
|gy
Materialtype:
血嵐
to™013
图6几何模型
图9计算参数设置
-Rrcpffhes
0（£
Oncel
ysb
图7基础模型参数设置
\777
4-■+
+*+
++.乂4放大
图10定义计算位置
f
八\/\j
//
EtieCllVC5IICS5C5
Esdi飾eePI阴応[“ml詡靶「仇«kh」ni2
DelcrfTicdkkrih
EFlf的t讪eIc即ecerrent129.J5*iD%
Ic:
xcEcrrtnlsceJl口Z.dDIvt亡si
图11变形网格
图8地基模型的初始应力场
AA
WflEEsEAfflaB
II;;:
'
II
II'
TH
IL
II'.
»——e——e——e——e——e——e——e——«
Totaldiipla^Einenli（Utol'）
EsliErreUlDt129.1.5*10-3m
图12变形总位移场
多PlflxisS.2Output'[Lrsson•Plarteforccsl
音AbEdrtGmEriryDefoTmatjonsForcesWrndowHelp
矍&弗昌QQ0歹a普［
Flale
Element
Node
X
[ml
¥
[m]
n
fkh抑]
Q
[kr/m]
M
[kWiT/m]
[kh/m]
1
1
433
l.DDD
4.000
0.979
1.239
D.DOQ
16.561
ItESCrt1'
291
（1.917
4,0（10
4.5ie
'16,251
D.254
2B753
型
4.000
a.S53
-aojfio
-0.693
3DJ59
293
（1.75fl
1,000
14.106
诩"?
29
”2.511
33120
294
0耐7
4.000
19.741
-担朝3
-4.903
35A53
2
294
IQQQ
19,563
-43,7{]9
-4,903
35161
WEEon1-
295
4.000
24.097
-41.063
-7.497
36.979
286
fl.jflO
4,000
29,412
-4ZJ10
<0,086
3B.6ia
297
0,417
4.QQ0
33.506
-39.470
-13.5SS
4D』巧
231
（1.3J3
1,000
37,573
”旺55B
<1,7^6
41,169
3
291
fljn
4.QQ0
36.993
-33.984
-iijq6
牡血
佬EE知1-
267
4,000
39-006
-25,915
-16,576
4ZrO65
26a
d.lfi?
4.QQ0
41J12
-17J74
-17.9^1
4232
269
（I.on
4,000
42,801
-9-005
-18,782
43,045
270
4.QOO
43.179
-0.065
-19.125
^3170
tJBPB
m~abpuTCfihCunef
d.m
图13
内力数据
™fileEditViewFormatWindowHelp
】Hp4CtkOjip5
SurMJjt^QE
Q&W日S購忆区
图14基础的荷载位移曲线
计算结果分析:
基础的沉降值随着距离荷载位置的增大而逐渐衰减；在350KN的均布荷载下，基础的最大沉降量是0.129m；最大弯矩是43,178KN/m;最大剪力是44.063KN。
算例二：
不排水条件下河堤的受力分析
20
201
ultimate
Ah
sand
一
6
=3
initialporepre^wra
ev&Tflporeprewurein肄ndr«yet
图J受水位变化影响的河堤
如图1所示的堤身高5m,由渗透性差的粘土构成。
其地基上部6m由软土构成，上层3m为一粘土（clay）层，下层3m为一泥炭（peat）层。
软土层基本上是不透水的，所以河流水位的短期变化并不影响这一部分的孔隙压力分布。
软土层下较大深度范围内为渗透性很好的砂土（sand）层，其上部4m包括在有限元模型之中。
这里假设砂土层中的水与河流的水相通，这意味着砂土层的水头紧密跟随河流的水位变化。
rSetasdefeUt
NSHt
图15新建项目设置
岩土工程数值分析读书笔记
图17划分网格
Plixi-s-i._SCalcfJI*EdieViAw匚*hculaceH*lp
Kp博
Inp-ucOwpuiCup*»0
Defet±..-]
JdHfi11AuCjD1=l
1P9-iad±!
rici
1£3rTFtam
1ll=iHlL
1Tmw
1WT
1F
IfpHax
□
o
Zg
Zg
OlOO-■■
c
c
■C'rtlHWIM
3
Q
pIabhir
1*4十mJmjl-t^^irrv
口8…
Q
图19初始应力的设置
riJ*E-dltvlsnES
■•A—Q—.wiI
aif.h—HK
hri/A
IdL。
■■匕JTil~fl
图20第二计算工序一般设置
1
图21第二计算工序定义孔隙水压力
HI
++
£应力篮£4血翹注謎宜
fl
#
#
fl
fl
fl
#
fl
图22孔隙水压场
It
jB/j"；
IB
I、\W'J^
“•疣恢%Will
打打
「■.………於"""
（/f;hnI{I
////?
?
rrF
ft
#
I
2贰
一~一
1—^
1
Nntt
—11—
\-
—-
PPL
\VI
4-
"1\
有效应力几乎为零区域
应
Un
4
图24水位上升后堤坝的有效应力
应图25堤址后有效应力横断面
I■■^1I
■B^a—
TKnT
卞0
ttt#
计算结果分析:
图24可以发现模型右边砂土层顶部的有效应力几乎为零，这是由于砂土层内孔隙压力的
增加造成的，从这一应力图，也可以发现土堤的变形在其后的粘土层中产生了一个被动应力区。
粘土和泥炭层的不排水性能导致了超静水压的发展。
图23可以看出堤坝在水位变化情况
下，土体的总变形，堤顶和堤址后土体隆起，故在施工过程中需要严格处理堤址后土体。
算例三水下土体开挖施工
横向支撐
10
开挖部分
A圭地下连续墙
10
鉴于几何模型的对称性，分析时只考虑其中一半（左边）。
开挖过程模拟为三个不同的阶段。
地下墙使用板单元来模拟，两侧地下墙和土的相互作用用界面单元来模拟，界面单元可以模拟墙体与土之间产生的滑动摩擦，此摩擦较土之间的摩擦要小些。
地下墙之间的横向支撑用弹簧单元来模拟。
lit
11
€>
■^sL
■w
材料定义:
砂土层、粘土层材料定义如上例,
不再叙述，地下连续墙和横向锚杆定义如下:
Materialsets
Plateproperties
图28定义地下连续墙的材料
Material
inch口「properties
翌IsadI
图29定义锚杆的材料
全局网格选用中等级别，在地下连续墙处网格加密。
图30划分网格
14/20
'八**■'_'L*:
:
■'***-n\
*（诗諾十「卄「
矚i蹦iS養管酸宦：
齣1讓SI攣壽辭iiSBiiffi匡Mffll曲
I「.J.,r-『MSI\■'"!
\
\■==:
三i：
="iS-ir'-'-'■
；；=.iHlmif.:
:
:
:
:
.・-"-IELI]21■■■■
!
m;:
：
:
一二二i：
itr：
^
：
ill1■:
lEI
'■■■J!
；；；
严卄耳IIIIh.■:
l»MFIIIII邺I■■二T・p・二"z:
ii盒待ziaiiIjj||||O1IXHniii-
III■■IIIEZ：
=：
==Z:
:
:
==-_
■M■.1■■■■-
l=:
ir:
:
:
==：
=：
==__■-.亠vi■iiiirii*>liliii<■*3IVLfllilllIIViiTiixii11aIIIIlBIB■■H■J匚二：
二：
：
二：
二-:
:
:
;；l*f>>>1-1■■«ILI
岂載萱盟阴議劭!
姿更：
：
IHiiiiisiSHiiiipiikiiz二:
;；iEEii：
iv=石1|.里r■—d；*:
■■HI"IV^ii|Vi-rv-H«i-i<--—：
■：
V■—M■■HHM■B■XI
j..l叩.■■II.*■■■bll■■中I
i盘iiii:
It!
■!
升级会员 