邓肯张模型

邓肯张模型
邓肯张模型

以土的常三轴实验学习Duncan-Chang本构关系模型

一、实验过程

1、试样制备

试验土样取自于南水北调焦作段一处工程,取回后,人工制成含水量15%的土体。在实验制样过程中,由于含水量较高,所以在通过制样器后,土柱未能成型,于是在原来土样基础上,添加了较干的土,再在制样器侧壁涂抹凡士林。最后制成高度7厘米,直径3.5厘米的土柱实验样品

2、不固结不排水(UU)剪切试验

试验是在土木工程学院深部矿井重点实验室进行的,试验装置如图1所示。

图1 常三轴实验仪

主要试验步骤为

(1)记录体变管的初始读数;

(2)对试样加周围压力,并在周围压力下固结。当孔隙水压力的读数接近零时,说明固结完成,记下排水管的读数;

(3)开动马达,合上离合器,按0.0065%/min的剪切应变速率对试样加载。按百分表读数为0,30,60,90,120,150,180,210,240,300,360,420,480,540,600,660,?,的间隙记读排水管读数和量力环量表读数,直到试样破坏为止。

二、邓肯张双曲线模型

到目前为止,国内外学者提出的土体本构模型不计其数,但是真正广泛用于工程实际的模型却为数不多,邓肯-张模型为其中之一。该模型是一种建立在增量广义虎克定律基础上的非线性弹性模型,可经反映应力~应变关系的非线性,模型参数只有8个,且物理意义明确,易于掌握,并可通过静三轴试验全部确定,便于在数值计算中运用,因而,得到了广泛地应用。

1、邓肯-张双曲线模型的本质

邓肯-张双曲线模型的本质在于假定土的应力应变之间的关系具有双曲线性质,见图2(a)。

图2(a ) 12()~a σσε- 双曲线

图2(b) 1131/()~εσσε-关系

图2 三轴试验的应力应变典型关系理论图

1963年,康纳(Kondner )根据大量土的三轴试验的应力应变关系曲线,提出可以用双曲线拟合出一般土的三轴试验1

3()~a σσε-曲线,即

a

a

b a εεσσ+=

-31 (1)

其中,,a b 为试验常数。对于常规三轴压缩试验,轴应变1a εε=。邓肯等

人根据这一双曲线应力应变关系提出了一种目前被广泛应用的增量弹性模型,一般被称为邓肯——张(Duncan-Chang )模型。

2.切线变形模量t E 邓肯—张计算公式

在常三轴压缩试验中,式(1)也可以写成

1

112

a b εεσσ=+- (2)

将常规三轴压缩试验的结果按

1113

~εεσσ-的关系进行调整,则二者近似成线

性关系,如图2(b)。其中,

a 为直线的截距;

b 为直线的斜率。

在常规三轴压缩试验中,由于2

30d d σσ==,所以切线模量为

132

11()()t d a

E d a b σσεε-==

+ (3)

在试验的起始点,10ε=,t

i E E =,则:

1

i E a

=

(4) 这表明

a 代表的是在这个试验中的起始变形模量(初始切线模量)的倒数。在

式(1)中,如果1

ε→∞,则:

131

()ult

b

σσ-= (5) 或者 131

()ult

b σσ=

- (6)

由此可看出b 代表的是双曲线的渐进线所对应的极限偏差应力13()ult σσ-的倒数。

在土的试样中,如果应力应变曲线近似于双曲线关系,则往往是根据一定的应变值(如115%ε=来确定土的强度13()f σσ-,而不可能在试验中使1ε无限大,求取13()ult σσ-;对于有峰值点的情况,取1313()()f σσσσ-=-峰,这样

1313()()f ult σσσσ-<-。定义破坏比f R 为:

破坏比:

1313()()f f ult

R σσσσ-=

- (7)

f R 值一般在0.75~1.0之间

13121

()()f ult f

R b σσσσ==

-- (8)

将式(8)、式(4)代入式(3)中,得

2

113111

()t f i i f E R E E εσσ??

????=

??+

??-????

(9)

式(9)中t E 表示为应变1ε的函数,使用时不够方便,可将t E 表示为应力的函数形式。从式(1)可以得到

13113()1()a b σσεσσ-=

-- (10)

将式(10)代入式(3),得

2

2

2

131********

1

()()1

11()1()1()t a

E ab b a a a b b b σσσσσσσσσσ=

=

=

??

??

??--++??

????------?

??

??

? (11) 将式(8)、式(4)代入式(11),得

2

13

131()t i f

f E E R σσσσ??

-=-??-????

(12)

三、实验数据处理

1、选取本实验中相关数据,绘制成如下表格

ε 1 σ1-σ3(kPa)

ε1/100 σ1-σ3/100(kPa)

ε1/σ1-σ3 0.68 31.5 0.0068 0.315 0.021587302 1.32 46.94 0.0132 0.4694 0.028121006 2 69.9 0.02 0.699 0.028612303 2.65 84.89 0.0265 0.8489 0.031216869 3.27 99.59 0.0327 0.9959 0.032834622 3.95 121.84 0.0395 1.2184 0.032419567 4.97 151.32 0.0497 1.5132 0.032844303 5.65 157.09 0.0565 1.5709 0.035966643 6.3 163.44 0.063 1.6344 0.038546256 6.95 169.63 0.0695 1.6963 0.040971526 7.58 175.81 0.0758 1.7581 0.043114726 8.25

189.29

0.0825

1.8929

0.043583919

8.9 195.03 0.089 1.9503 0.045634005

9.57 208.05 0.0957 2.0805 0.045998558

10.2 206.48 0.102 2.0648 0.049399458

10.9 219.24 0.109 2.1924 0.049717205

11.53 217.45 0.1153 2.1745 0.053023684

12.2 223.01 0.122 2.2301 0.054706067

12.88 235.04 0.1288 2.3504 0.054799183

13.5 240.23 0.135 2.4023 0.056196145

14.2 245.28 0.142 2.4528 0.05789302

14.9 249.89 0.149 2.4989 0.059626236

15.6 254.68 0.156 2.5468 0.061253338

16.23 252.53 0.1623 2.5253 0.064269592

16.85 257.33 0.1685 2.5733 0.065480123

17.5 255.37 0.175 2.5537 0.068528018

18.15 253.26 0.1815 2.5326 0.071665482

18.77 257.84 0.1877 2.5784 0.072797083 2、利用Origin9.0软件绘制图像如下图所示

围压

图3(a )

131~-εσσ)( 双曲线 0.00

0.020.040.060.080.100.120.140.160.180.20

0.00

0.010.020.030.040.050.060.070.08ε1/σ1-σ3

ε1/100

( ε1/σ1-σ3)— ε1实验数据关系图—— ( ε1/σ1-σ3)— ε1 理论直线图

图3(b ) 1131/()~εσσε-关系

图3 三轴试验的应力应变典型关系实验图

将以上两图与图2(a)、(b)对比发现,实验数据能够较好的反映出Duncan-Chang模型典型的双曲线图。

四、实验与课程结语

通过查阅相关文献和书籍得知,Duncan-Chang双曲线模型能较好地反映土体的非线性性态,概念清楚,易于理解,在岩土工程和地下工程的数值分析中得到广泛应用。在学习《高等岩土力学》这门课之后,安排常三轴实验对于Duncan-Chang双曲线模型有了更深刻的理解。通过实验发现,常三轴实验所得数据能够很好的拟合出公式所给的Duncan-Chang典型双曲线的图形,说明此次实验是比较成功的。但是由于时间和材料有限,本次的实验次数太少,故而仅能得到如上的结论。

另外,在阅读相关文献时候,发现在不少文中讲到实际工程应用中,作具体数值计算时,要用到具体的参数值,而一些文献、书籍和资料中给出了一般性的土的Duncan-Chang模型参数,但是考虑到具体工程需要具体分析,所以对于今后的学习和工作需要注意。

五、参考文献

[1] 张海霞,陈勇,秦帮民.三轴不固结不排水的实验方法.水利水电进展[J].2007.

[2] 冯卫星,常绍东,胡万毅. 北京细砂土邓肯张模型参数实验研究. 岩石力学与工程学报. 1999.

[3] 钱家欢. 土力学[M] . 南京: 河海大学出版社, 1995: 123O132.

[4] 王保田, 张福海. 土力学与地基处理[M] . 南京: 河海大学出版社, 2005: 116O127.

塑性力学读书报告 邓肯-张模型研究认识 学院:建设工程 姓名:王吉亮 学号:2006631011 专业:地质工程

教师:金英玉

邓肯-张模型研究认识 王吉亮(83分) 摘 要:从邓肯-张模型的本源开始,分析研究了邓肯-张模型与E-B 模型的建立过程和模型中参数如何确定的问题,结合对该模型的认识,提出该模型具有的缺点与不足。 关键词:邓肯-张模型;E-B 模型;参数确定 CONGNITION ON THE STUDY OF DUNCAN-CHANG MODEL Wang Jiliang Abstract: rom the parent of Duncan-Chang model, studing the establish procedure of Duncan-Chang model and E-B model, introducing the problem of how to define the indexes in the model. Associate the congnition on this model, present the shortcomings. Keywords: Duncan-Chang model; E-B model; indexes define 1 引言 邓肯-张模型是一个非线性本构模型,既然是一个本构模型,可想而之他反应的是应力与应变之间的关系。说它是非线性的,那么反映应力应变关系的模量就不是一个常数E那么简单。在介绍该模型之前,先要介绍一个概念,就是反映非线性关系的增量广义胡克定律: 1123()t t t v d d d d E E σεσσ= -+ (1) 1963年,康纳(Kondner )根据大量土的三 轴试验的应力应变关系曲线,提出可以用双曲线拟合出一般土的三轴试验13()~a σσε-曲线,即: 13a a a b εσσε-= + (2) 其中,a 、b 为试验常数。对于常规三轴压缩试验,1a εε=。邓肯等人根据这一双曲线应力应变关系提出了一种目前被广泛的增量弹性模型, 一般被称为邓肯-张(Duncan-Chang )模型。 在常规三轴压缩试验中,13a a a b εσσε-=+可以写成: 1113 a b εεσσ=+- (3) 将常规三轴压缩试验的结果按 11 13 ~εεσσ-的关系进行整理,则二者近似成线性关系(见图1)。其中,a 为直线的截距;b 为直线的斜率。 在常规三轴压缩试验中,由于 230d d σσ==,所以切线模量为 ε1/(σ1 -σ3 ) -σ3 )ult 图1 1113 ~εεσσ-线性关系图 132 11()() t d a E d a b σσεε-= =+ (4) 在试验的起始点,10ε=,t i E E =,则: 1 i E a = ,这表明a 表示的是在这个试验中的起始变形模量E i 的倒数。如果1ε→∞,则: 131 ()ult b σσ-= (5)

5.370569 应力差(б1-б3/100kPa轴向变形ε1 体应变εv ε3ε1/(б1-б3 0.5080.002250.00074-0.0007554.42913E-051.0020.004490.0013- 0.0015954.48104E-051.4630.006740.00176-0.002494.60697E-051.8490.008980.00223-0.0033754.85668E-052.1490.011230.00223-0.00455.22569E-052.3310.013480.00214-0.005675.78293E-052.4770.015720.00176-0.006986.34639E-052.5870.017970.00158-0.0081956.94627E-052.6650.020210.00111-0.009557.58349E-052.730.022460.00056-0.010958.22711E-052.7770.024710.00009-0.012318.89809E-052.8120.026950.00003-0.013469.58393E-052.8450.03032-0.00012-0.015220.0001065732.8780.03369-0.00195-0.017820.000117062.8890.03706-0.00288-0.019970.000128282.8940.04043-0.00381-0.022120.0001397032.8890.04492-0.00474-0.024830.0001554862.8790.04941-0.00567-0.027540.0001716222.8690.0539-0.0065-0.03020.000187872.8530.05839-0.00752- 0.0329550.0002046622.8360.06289-0.00827-0.035580.0002217562.809 0.06738 -0.00891 -0.038145 0.000239872 应力差(б1-б3/100kPa轴向变形ε1 体应变εv ε3ε1/(б1-б3 0.9090.001250.00075-0.000250.1375137512.2860.00350.00151- 0.0009950.1531058623.5330.005750.00236-0.0016950.1627512034.3410.0080.00302-

邓肯-张模型是一个非线性本构模型,既然是一个本构模型,可想而之他反应的是应力与应变之间的关系。说它是非线性的,那么反映应力应变关系的模量就不是一个常数E那么简单。在介绍该模型之前,先要介绍一个概念,就是反映非线性关系的增量广义胡克定律: 1123()t t t v d d d d E E σεσσ= -+ (1) 1963年,康纳(Kondner )根据大量土的三 轴试验的应力应变关系曲线,提出可以用双曲线拟合出一般土的三轴试验13()~a σσε-曲线,即: 13a a a b εσσε-= + (2) 其中,a 、b 为试验常数。对于常规三轴压缩试验,1a εε=。邓肯等人根据这一双曲线应力应变关系提出了一种目前被广泛的增量弹性模型,一般被称为邓肯-张(Duncan-Chang )模型。 在常规三轴压缩试验中,13a a a b εσσε-=+可以写成: 1113 a b εεσσ=+- (3) 将常规三轴压缩试验的结果按 11 13 ~εεσσ-的关系进行整理,则二者近似成线性关系(见图 1)。其中,a 为直线的截距;b 为直线的斜率。 在常规三轴压缩试验中,由于 230d d σσ==,所以切线模量为 ε1 /(σ1 -σ3 ) 1 b=1/(σ1 -σ3 )ult   a =1/E i 图1 1113 ~εεσσ-线性关系图 132 11()() t d a E d a b σσεε-= =+ (4) 在试验的起始点,10ε=,t i E E =,则: 1i E a = ,这表明a 表示的是在这个试验中的起始变形模量E i 的倒数。如果1ε→∞,则: 131 ()ult b σσ-= (5) 由此可以看出b 代表的是双曲线的渐近线所对应的极限偏差应力13()ult σσ-的倒数。 在土的试样中,如果应力应变曲线近似于双曲线关系,则往往是根据一定的应变值(如 115%ε=)来确定土的强度13()f σσ-,而不可 能在试验中使1ε无限大,求取13()ult σσ-;对于有峰值点的情况,取1313()()f σσσσ-=-峰, 这样1313()()f σσσσ--ult <。定义破坏比R f 为: 1313()()f f R σσσσ-=-ult (6) 而 13131 ()()f f R b σσσσ== --ult (7) 将上式与1 i E a = 代入 132 11()() t d a E d a b σσεε-==+ (8) 得到:

邓肯张模型FORTRAN子程序源代码 SUBROUTINE UMA T(STRESS,STA TEV,DDSDDE,SSE,SPD,SCD, 1 RPL,DDSDDT,DRPLDE,DRPLDT,STRAN,DSTRAN, 2 TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,MA TERL,NDI,NSHR,NTENS, 3 NSTA TV,PROPS,NPROPS,COORDS,DROT,PNEWDT,CELENT, 4 DFGRD0,DFGRD1,NOEL,NPT,KSLAY,KSPT,KSTEP,KINC) C INCLUDE 'ABA_PARAM.INC' C CHARACTER*80 MA TERL DIMENSION STRESS(NTENS),STA TEV(NSTA TV), 1 DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS), 2 STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1), 3 PROPS(NPROPS),COORDS(3),DROT(3,3), 4 DFGRD0(3,3),DFGRD1(3,3) C DIMENSION PS(3),DSTRESS(NTENS),TDSTRESS(NTENS),TSTRESS(NTENS) PARAMETER (ONE=1.0D0,TWO=2.0D0,THREE=3.0D0,SIX=6.0D0) K=PROPS(1) N=PROPS(2) RF=PROPS(3) C=PROPS(4) FAI=PROPS(5)/180.0*3.1415926 G=PROPS(6) D=PROPS(7) F=PROPS(8) KUR=PROPS(9) PA=PROPS(10) DFAI=PROPS(11)/180.0*3.1415926 S1S3O=STA TEV(1) S3O=STA TEV(2) SSS=STA TEV(3) CALL GETPS(STRESS,PS,NTENS) FAI=FAI-DFAI*LOG10(S3O/PA) CALL GETEMOD(PS,K,N,RF,C,FAI,ENU,PA,KUR,EMOD,S,S3O,G,D,F 1 ,SSS,S1S3O) EBULK3=EMOD/(ONE-TWO*ENU) EG2=EMOD/(ONE+ENU) EG=EG2/TWO EG3=THREE*EG ELAM=(EBULK3-EG2)/THREE CALL GETDDSDDE(DDSDDE,NTENS,NDI,ELAM,EG2,EG)

研究生课程作业邓肯张模型参数计算 学生姓名李俊 学科专业岩土工程 学号201420105614 任课教师周小文教授 作业提交日期2014年12月

1.计算轴向应变 c h h ?∑= 1ε 式中 1ε-轴向应变; h ?∑-固结下沉量,由轴向位移计测得 0h -土样初始高度 c h —按实测固结下沉的试样高度 c h ?—试样固结下沉量 2.计算按实测固结下沉的试样高度,面积: 式中 Ac -按实测固结下沉的试样面积 0V -土样初始体积 3.计算剪切过程中试样的平均面积: 式中 a A -剪切过程中平均断面积 c V -按实测固结下沉的试样的体积 i V ?-排水剪中剪切时的试样体积变化 按体变管或排水管读数求得 1h ?-固结下沉量,由轴向位移计测得 3. 计算主应力差 c i c h V V A ?-= 01 h h V V A c i c a ?-?-= C c c A h V ?=

103 1?=-a A CR σσ 式中 31σσ- - 主应力差 1σ―大主应力 3σ-小主应力 C -测力计率定系数 R -测力计读数 2 数据处理 2.1 3σ=100kPa 数据初步计算 当3σ=100kPa 时,各数据初步计算如表1所示。 围压100kPa 数据初步计算表 表1

2.1.1 由切线模量计算数据 对公式 ) (311 σσε-=a +b 1ε进行直线拟合,如图1所示。 图1 1131 /()~εσσε-拟合曲线 a =0.0002,1 i E a = =5000kPa b ==0.0028,()131 ult b σσ-= =263.16kPa ()13f σσ-=204.26kPa ,()()1313f f ult R σσσσ-= -=0.7762 2.1.2 由泊松比计算数据 对公式()313/f D εεε-=+-进行直线拟合,如图2所示。

邓肯-张模型的关键点是材料的弹性模量随大小主应力差及小主应力(围压)的变化而变化,用APDL实现之的基本思路是:给每个单元定义一个材料号,分级施加荷载,在每个荷载步结束时提取出各单元的大小主应力,据此计算出下个荷载步的弹性模量Et,修改各单元之MP,用于下一步计算。 以下是一个简单算例,copy出去可直接运行。 !!!常规三轴试验模拟 !!!by taomingxing,NWPU !!!2003.7.16 FINISH /CLEAR /TITLE,Numerical Simulation of three axes testing of soils /PREP7 *dim,SUy,array,50 !Settlement records *dim,MaxPs,array,120 !Max history p1-p3 *dim,MaxDs,array,120 !Max history Ds !*dim,EEt,array,50 !Et of elememt !!!Duncan-Chang Model !!!Symbols:c-粘滞力,Fai-内摩擦角,Sf-破坏强度(p1-p3)f,Ds-应力水平,Pa-大气压,P3-围压 *CREATE,Duncan-Chang !Creat Macro file *afun,deg !Unit of angle *set,Pa,1e5 *set,P1,-ArrS3(i) !注意:岩土工程中应力为拉负压正

*set,P3,-ArrS1(i) *if,P3,LT,0.1*Pa,then P3=0.1*Pa !围压最小取值 *endif Sf=2*(c*cos(Fai)+P3*sin(Fai))/(1-sin(Fai)) !Mohr-Coulomb破坏强度(p1-p3)f Ds=(P1-P3)/Sf !应力水平, *if,Ds,GT,0.95,then Ds=0.95 !应力水平最大取值 *endif !判断加卸荷,如果(P1-P3)小于历史最大值视为卸荷-再加荷过程 *if,MaxPs(i),LT,P1-P3,then Ei=k*Pa*(P3/Pa)**n Et=Ei*(1-Rf*Ds)**2 !加荷情况的切线模量 MaxPs(i)=P1-P3 !保存历史最大应力 *elseif,MaxPs(i),GE,P1-P3 Et=Kur*Pa*(P3/Pa)**n !卸荷模量 *endif mp,ex,i,Et !修改单元i的Et mp,nuxy,i,Mu *END !!!单元类型 et,1,42 !平面四节点单元

以土的常三轴实验学习Duncan-Chang本构关系模型 一、实验过程 1、试样制备 试验土样取自于南水北调焦作段一处工程,取回后,人工制成含水量15%的土体。在实验制样过程中,由于含水量较高,所以在通过制样器后,土柱未能成型,于是在原来土样基础上,添加了较干的土,再在制样器侧壁涂抹凡士林。最后制成高度7厘米,直径3.5厘米的土柱实验样品 2、不固结不排水(UU)剪切试验 试验是在土木工程学院深部矿井重点实验室进行的,试验装置如图1所示。 图1 常三轴实验仪

主要试验步骤为 (1)记录体变管的初始读数; (2)对试样加周围压力,并在周围压力下固结。当孔隙水压力的读数接近零时,说明固结完成,记下排水管的读数; (3)开动马达,合上离合器,按0.0065%/min的剪切应变速率对试样加载。按百分表读数为0,30,60,90,120,150,180,210,240,300,360,420,480,540,600,660,?,的间隙记读排水管读数和量力环量表读数,直到试样破坏为止。 二、邓肯张双曲线模型 到目前为止,国内外学者提出的土体本构模型不计其数,但是真正广泛用于工程实际的模型却为数不多,邓肯-张模型为其中之一。该模型是一种建立在增量广义虎克定律基础上的非线性弹性模型,可经反映应力~应变关系的非线性,模型参数只有8个,且物理意义明确,易于掌握,并可通过静三轴试验全部确定,便于在数值计算中运用,因而,得到了广泛地应用。 1、邓肯-张双曲线模型的本质 邓肯-张双曲线模型的本质在于假定土的应力应变之间的关系具有双曲线性质,见图2(a)。

图2(a ) 12()~a σσε- 双曲线 图2(b) 1131/()~εσσε-关系 图2 三轴试验的应力应变典型关系理论图

ANSYS邓肯-张材料模型 楼主给的在ANSYS上实现邓肯-张模型的方法很有用, 但其中还有几点需要修正的,这也是楼上的兄弟们有疑问的原因。我把楼主的代码运行了一下,然后对照作了修改,现在上传一下,有问题的兄弟可以仔细对照一下,在这里我对其中几个比较明显的问题说明一下: 1.MP命令不能直接给单元加材料,这是对的。在这里,楼主遗漏了一下命令:MPCHG,具体见下面的修改过的代码。 2.关于密度的问题。这些要在宏中定义,每修改一种材料(即调用一次邓肯-张子程序)就要修改一次材料的密度,其他有关材料的问题可以类推。 3.关于施加重力的问题。要在调用宏后,在同一个循环中重新定义一下重力。 以下是我修改过的楼主的代码,希望对兄弟们有所帮助。 !用APDL得到初步成果,贴于此供感兴趣的朋友参考,不当之处敬请指正,!欢迎加以完善。 !基本思路: !邓肯-张模型的关键点是材料的弹性模量随大小主应力差 !及小主应力(围压)的变化而变化,用APDL实现之的基本思路是: !给每个单元定义一个材料号,分级施加荷载,在每个荷载步结束时提取出各 !单元的大小主应力,据此计算出下个荷载步的弹性模量Et,修改各单元之MP,!用于下一步计算。 !以下是一个简单算例,copy出去可直接运行。 !!!常规三轴试验模拟

!********************************************************** FINISH /CLEAR /TITLE,Numerical Simulation of three axes testing of soils /PREP7 *dim,SUy,array,50!Settlement records *dim,MaxPs,array,120!Max history p1-p3 *dim,MaxDs,array,120!Max history Ds !*dim,EEt,array,50!Et of elememt !!!Duncan-Chang Model !!!Symbols:c-粘滞力,Fai-内摩擦角,Sf-破坏强度(p1-p3)f, !Ds-应力水平,Pa-大气压,P3-围压 !********************************************************************** *CREATE,Duncan-Chang!Creat Macro file *afun,deg!Unit of angle *set,Pa,1e5 *set,P1,-ArrS3(i)!注意:岩土工程中应力为拉负压正 *set,P3,-ArrS1(i) *if,P3,LT,0.1*Pa,then P3=0.1*Pa!围压最小取值 *endif Sf0=2*(c0*cos(Fai)+P3*sin(Fai))/(1-sin(Fai))!Mohr-Coulomb破坏强度

相关文档
最新文档