23-04-23-材料参数-泊松比
1.3轴向拉压杆横截面上的变形

目录引言 (2)一杆件受拉压的内力、应力、变形 (2)1.1轴向拉压的内力、轴力图 (2)1.2 轴向拉压杆横截面上的应力 (5)1.3 轴向拉压杆横截面上的变形 (7)1.4 圣维南原理 (9)1.5 工程结构实例分析 (11)二圆轴扭转 (15)2.1、扭转的力学模型及ANSYS建模 (15)2.2、圆轴扭转时,横截面上的内力偶矩------扭矩 (15)2.3、圆轴扭转时,横截面上的应力、强度条件 (15)(1) 横截面上的切应力 (15)(2) 极惯性矩与抗扭截面系数 (15)三、梁弯曲的内力、变形、应力 (20)3.1 梁的弯曲内力、变形 (20)3.2 弯曲应力 (27)3.3 工程实例: (31)四、压杆稳定 (35)4.1、压杆稳定的概念 (35)4.2、临界压力 (35)4.3、三类压杆的临界载荷 (36)4.4、压杆稳定性计算 (36)4.5 工程实例4 (38)引 言《材料力学》是机械、土木类工科学生重要的技术基础课,其计算方法和思想在工程计算中应用非常广泛。
为了使学生对课内知识体系有一个比较清晰的感性认识,锻炼学生的求真精神和实践动手能力,进一步培养学生的综合创造力,兴趣小组的学生们在教师的指导下基于ANSYS 有限元分析软件对《材料力学》的某些知识点进行数值计算与模拟,得到相关的数据、云图或动画,从而对理论公式进行形象验证,更开阔了学生的视野,提高了学生的CAE 水平。
本研究内容包括三部分:(1)对《材料力学》课程中的基本内容,包括拉压、剪切、扭转、弯曲的内力、应力、变形、压杆稳定、动载荷、疲劳强度、圣维南原理等重要理论知识点情况通过ANSYS 进行分析,得到内力、变形、应力、应变相关的数据、云图或动画;(2)对重要知识点的典型例题通过ANSYS 进行计算,并与理论计算结果进行对比验证。
(3)对《材料力学》理论知识能够解决的典型工程实际问题进行建模、分析与计算。
一 杆件受拉压的内力、应力、变形1.1轴向拉压的内力、轴力图在工程结构和机械中,发生轴向拉伸或压缩的构件是很常见的。
例题4——材料非线性与几何非线性分析

研究生上机习题课4-1本例题属于平面应力问题。
采用 mm ,ms,kN,kg单位制。
1、运行Ansys 9.0 ——Ansys Product Launcher——License——Ansys Multiphysics——点击FileManagement ,点击Working directory最右边的… ,改变工作目录,d:\teaching,在Initial Jobname 中输入Example41,点击Run。
2、点击Preference——选择Structural,OK。
3、创建几何模型。
4、做矩形。
Preprocessor——Modeling——Create——Areas——Rectangle——By Dimensions——X输入0,330;Y输入0,60。
5、作半圆。
移动当前坐标系,功能菜单Work Plane——Offset Wp to Keypoint,点击矩形右边的两个关键点,则当前坐标系移动到矩形右边的中点。
Preprocessor——Modeling——Create——Areas——circle——By Dimensions——外部半径输入30,内部不输;角度分别输入-90,90。
6、作1/4圆。
移动当前坐标系,功能菜单Work Plane——Offset Wp to Global Origin,则当前坐标系移动到原点。
Preprocessor——Modeling——Create——Areas——circle——By Dimensions——外部半径输入120,内部不输;角度分别输入-90,0。
7、将3个面积加起来。
Preprocessor——Modeling——Operate——Booleans——Add——Areas,全部选择,Pick All。
8、做r=15的倒角。
Preprocessor——Modeling——Create——Line——Line Fillet,点击1/4圆的边和矩形的下边。
一种极性聚合物负泊松比泡沫材料及其制备方法_CN109734955A

5 .根据权利要求1所述的极性聚合物负泊松比泡沫材料的制备方法,其特征在于,所述 物理发泡剂为二氧化碳、氮气、丁烷、戊烷中的一种或多种的混合物。
备方法 ( 57 )摘要
本发明 涉 及一 种极性聚 合 物负泊松比 泡 沫 材料及其制备方法,属于功能性聚合物泡沫材料 技术领域。该方法解决了现有技术中负泊松比聚 合物泡沫材料的制备方法制备过程繁琐、制备效 率 低 、对聚 合物有限 制、不 利于大规模生 产的 技 术问 题。本发明的 制备方法 ,将极性聚合物直接 成型或 与助剂复合后成型 ,制成待发泡 制件 ,置 于发泡设备内 ,引入物理发泡剂和水 ,在80~180 ℃的饱和温度以及3~25MPa的饱和压力下,饱和 时间0 .5h以上,快速泄压,得到泡沫材料。该方法 可以一步制备极性聚合物负泊松比泡沫材料,制 备工艺简单 ,制备效率高 ,使 用水辅助完成制备 , 安全性高 ,绿色环保,可制备尺寸较大的负泊松 比泡沫材料,适合大规模生产。
9 .根据权利要求1所述的极性聚合物负泊松比泡沫材料的制备方法,其特征在于,所述 快速泄压的泄压速度为2~50MPa/s。
10 .权利要求1~9任何一项所述的极性聚合物负泊松比泡沫材料的制备方法制备的极 性聚合物负泊734955 A
说 明 书
1/5 页
一种极性聚合物负泊松比泡沫材料及其制备方法
(74)专利代理机构 长春众邦菁华知识产权代理 有限公司 22214
代理人 王莹
(51)Int .Cl . C08J 9/12(2006 .01) C08L 23/28(2006 .01) C08L 77/00(2006 .01) C08L 23/08(2006 .01)
HXD2机车牵引装置结构及其强度研究

关键词:HXD2 机车;牵引装置;结构;强度校核;应力测试;强度计算准则
设计工作。
中图分类号:U260.331+7;U264.2+23 文献标识码 :A 文章编号:1000-128X(2010)04-0005-02
Study on Structure and Strength of HXD2 Locomotive Traction Device
2 牵引装置强度校核
按UIC 615 的要求:在机车冲击情况下,牵引装置 可承受 3g 以下纵向加速度,而不会产生变形或损伤; 在机车受异常冲击在 3g ~5g 纵向加速度时,牵引支座 损伤以保护转向架。由于两端牵引支座的受力部位的 结构相同,所以仅以车体牵引支座为例计算。
牵引装置计算用参数见表 1。按 3 种工况进行校核 (见表 2)。
表 2 牵引装置计算工况表
项目
工况 1 工况 2 工况 3
描述
受启动牵引力 受 3 g 纵向加速度 受 5 g 纵向加速度
表 3 各种工况下计算过程数据及结果
项目
工况 1
数据及结果
工况 2
工况 3
牵引杆应力图 牵引支座应力图 部件受纵向力 / N 沿牵引杆向力 / N 牵引杆最大应力 / N·m- 2 牵引杆安全系数 校核结果 牵引支座最大应力 / N·m- 2 牵引支座安全系数 校核结果
Key words: HXD2 locomotive; traction device; structure; strength check; stress measurement; calculation criterion of strength
0 引言
HXD2 型电力机车转向架牵引装置采用中间推挽式 牵引杆与车体连接,牵引杆倾斜布置,与轨道成 11 ° 角。该牵引方式在我国大功率交流传动电力机车上首 次使用,为此有必要对牵引装置进行深入系统的研究。 本文对HXD2 机车转向架牵引装置结构进行介绍,提出 其强度校核准则并进行强度设计校核,通过机车的线 路试验验证其安全性。
迈达斯-有限元模态分析

2 3
Modal Analysis 39
Step
09 网格 >> 属性>>材料
操作步骤
1. 创建 >> 各向同性. 2. 输入材料属性.
号 名称 弹性模量 泊松比 质量密度
2 Matl 2e11 (N/m2) 0.3 8.245 (kg/m3)
3. 点击[确认]. 4. 点击[关闭].
2 2
1
4 3
模态分析选项
• 模态数量: 输入需要计算的自然频率数量.
• 最低/最高: 设置频率范围. 仅计算该范围内的频率.
• Sturm 序列检查: 通过迭代过程计算丢失特征值.
Modal Analysis 4
Step
00 概述
概述
线性静态分析 - 单位 : N, mm - 各向同性弹性材料 - 几何模型: Hanger.x_t
>> 变形+未变形(特征线) .
3. 分析和结果>> 模态(必需)>>
MODE1 >> 总位移.
1 2
Modal Analysis 16
Step
11 分析和结果>> Hanger >> 特征值分析结果表
操作步骤
1. 双击[特征值分析结果表]. 1
Modal Analysis 17
Step
00 概述
操作步骤
1. 坐标 : 输入“(0.125), <0, 0.5>” . 2.在工作窗口中单击鼠标右键. 3. 坐标 :输入 “(0.25), <0, 0.5>” . 4.在工作窗口中单击鼠标右键. 5. 坐标 :输入 “(0.375), <0, 0.5>” . 6.在工作窗口中单击鼠标右键. 7. 点击[取消] .
桥头跳车力学分析及台背刚柔过渡设计参数

重心位置 Λ
0. 3 0. 3 0. 3 0. 3
过渡半径 R m 667 333 200 100
最大加速度 a
m s2 ±1. 158 ±2. 32 ±3. 86 ±7. 72
后轴荷载增大系数 Κ
1. 115 1. 232 1. 386 1. 772
注: + 号表示加速度方向向上, - 号表示向下, 分别为车辆从桥头→搭板和搭板→路面时的离心加速度, 二者大小相等方向相反。
m
g
F2 (1-
Λ) = 1+
a g
计算中取搭板长度为 5 m , 车速为 100 km h,
计算结果列于表 1。从式 (1) 和表 1 的荷载作用力
桥头设置搭板后, 桥台不发生沉降, 一旦路堤发 分析结果可以看出, 由于桥台和路堤的强度和刚度
生沉降后, 桥台→搭板→路堤三者路线纵坡呈折线, 存在很大的差异, 路堤发生一定沉降后, 在设置桥
模量平均达 172 M Pa, 计算得到复合地基压缩模量 达 179 M Pa, 为加固前的 3 倍以上, 加固效果明显, 达到了设计意图, 实现了刚柔过渡。
参照文献[ 2 ]的研究, 取车辆后轴弹簧刚度为 k = 2 060 kN m。 令弹簧受激励后的最大压缩量为 △x , 取沉降差 h 为 15 mm , 在弹簧的弹性范围内有 △x = h, 路面板受力 F 由 2 部分组成: 后轴实际作 用力 G 和附加冲击作用力 d F。由模型可知, 附加冲 击力 d F 可由下式计算得到:
并结合湖南省益常高速公路的工程实际情况, 提出 了桥头刚柔过渡设计的部分设计参数, 供大家参考。
2 汽车荷载分析模型和受力分析 根据笔者多年的工程实践和研究, 桥头台背回
Abaqus切削仿真材料库

Abaqus切削仿真材料库:在abaqus中定义材料的参数共有4栏:generalmechanicalthermalothergeneral下面的子选项卡有:◆Density 密度◆Depvar◆Regularization◆User defined field◆User output variablesMechanical下面的子选项卡有◆Elasticity◆Plasticity◆Damage for ductile metals◆Damage for traction separated laws◆Damage for fiber reinforced composites ◆Damage for elastomers◆Deformation plasticity◆Damping◆Expansion◆Brittle cracking◆Eos◆ViscosityThermal 下面的子选项卡有◆Conductivity◆Heat generation◆Inelastic heat fraction◆Joule heat fraction(焦耳热比例)◆Latent heat(汽化潜热)◆Specific heatOther下面的子选项卡有◆Acoustic medium◆Electrical◆Mass diffusion◆Pore fluid◆GasketElastiityc下面的子选项卡有◆Elastic◆Hyperelastic◆Hyper foam◆Low density foam◆Hyper elastic◆Porous elastic◆ViscoelasticPlasticity下面的子选项卡有◆PlasticPlastic又分为以下几类IsotropicKinematicJohnson-CookUserCombinedPlastic下面又有suboptionsRate dependentPotentialCyclic hardeningOrnlCycled plasticAnneal temperature◆Cap plasticity◆Cap iron plasticity◆Clay plasticity◆Concrete damaged plasticity◆Concrete smeared cracking◆Crushable foam◆Drucker prager◆Mohr coulomb plasticity◆Porous metal plasticity◆Creep◆Swelling◆ViscousDamage for ductile metals下面又有以下选项卡 Ductile damageJohnson-cook damageShear damageFLD damageFLSD damageM-K damageMSFLD damage目前接触过下列部分材料的参数AISI434045号钢AISI D2铝合金7055T-7451Ti6Al4V以下单位制为mm、吨、mw、N、s、Mpa⏹AISI 4340密度density:7830e-12吨/mm3弹性模量E: 207000Mpa泊松比μ:0.3比热容:480J/Kg.℃=480e6mJ/ton.℃导热系数conductivity:44W/m.℃=44mw/mm.℃延展系数expansion:1.35e-5m/m.℃=1.35e-5mm/mm.℃熔化温度:1520℃室温:25摄氏度J-C塑性参数:A=792Mpa;B=510Mpa;n=0.26;C=0.014;m=1.03 J-C本构方程:J-C损伤参数:D1=0.05; D2=3.44; D3=-2.12; D4=0.002; D5=0.61J-C 损伤准则(damage initiation)p-静水压力;q-等效应力;p/q应力三轴度stress-triaxialityAISI 4340相当于中国的40CrNiMo脆性断裂韧度75Mpa.m1/2⏹45号钢密度density:7800e-12弹性模量E:209000Mpa泊松比μ:0.3延展系数:1.17e-5比热容:456J/Kg.℃=456e6mJ/ton.℃J-C塑性参数:A=496Mpa;B=434Mpa;n=0.307;m=0.804;C=?0.0084参考应变率2e-4⏹7050-T7451静力拉伸性能准静态压缩实验与SHPB高速冲击实验J-C塑性参数:A=490;B=206.9;n=0.344;C=0.005;m=1.80,参考应变率1e-3熔点:600℃,室温:25℃注:在abaqus中必须用真实应力-真实应变定义塑性。
基于MatDEM的黄泥扁滑坡支护效果分析

第42卷 第5期2023年 9月 地质科技通报B u l l e t i n o f G e o l o g i c a l S c i e n c e a n d T e c h n o l o g yV o l .42 N o .5S e p .2023王飞,徐楚,黎伟.基于M a t D E M 的黄泥扁滑坡支护效果分析[J ].地质科技通报,2023,42(5):52-60.W a n g F e i ,X u C h u ,L i W e i .S u p p o r t i n g e f f e c t a n a l y s i s o f t h e H u a n g n i b i a n l a n d s l i d e b a s e d o n M a t D E M [J ].B u l l e t i n o f G e o l o gi c a l S c i e n c e a n d T e c h n o l o g y,2023,42(5):52-60.基金项目:国家自然科学基金项目(42102327)作者简介:王 飞(1999 ),男,现正攻读地质工程专业硕士学位,主要从事岩土工程数值模拟的研究工作㊂E -m a i l :764731655@q q.c o m 通信作者:黎 伟(1989 ),男,副教授,主要从事特殊土工程力学特性及其加固改良技术工作㊂E -m a i l :l i w e i 890508@126.c o m©E d i t o r i a l O f f i c e o f B u l l e t i n o f G e o l o g i c a l S c i e n c e a n d T e c h n o l o g y .T h i s i s a n o pe n a c c e s s a r t i c l e u n d e r t h e C C B Y -N C -N D l i c e n s e .基于M a t D E M 的黄泥扁滑坡支护效果分析王 飞,徐 楚,黎 伟(中国地质大学(武汉)工程学院,武汉430074)摘 要:黄泥扁滑坡是一个小型中层土质滑坡,已出现滑移的迹象,急需对其进行支护治理,分析该滑坡的破坏机理㊁验证现有的支护结构是否满足滑坡的稳定性要求以及探索合理且有效的支护部位和支护形式具有十分重要的意义㊂基于M a t D E M 建立了黄泥扁滑坡在20年一遇洪水位下的三维模型,模拟了滑坡在无支护㊁挡土墙支护㊁挡土墙和抗滑桩共同支护3种情况下的变形,并根据位移场分析了不同支护方案的支护效果㊂模拟结果显示:①该滑坡在无支护的情况下处于不稳定状态;②在无支护的情况下,滑体从滨河路的坡脚处剪出,在有挡土墙支护的情况下,滑体从快速通道的坡脚处剪出,表明现有的挡土墙支护结构并不满足滑坡的稳定性需求,需要对快速通道的坡脚进行重点支护;③采用抗滑桩对快速通道的坡脚进行支护,并在抗滑桩之间加设挡土板,此方案不仅在很大程度上提高了滑坡的整体稳定性,而且有效地解决了部分土体从桩间滑落的问题㊂研究成果可以为后续的工程设计㊁施工等提供一定的参考,同时也验证了M a t D E M 在滑坡分析中的适用性㊂关键词:M a t D E M ;离散元模拟;挡土墙;抗滑桩;支护效果;黄泥扁滑坡2022-01-11收稿;2022-03-24修回;2022-04-18接受中图分类号:P 642.22 文章编号:2096-8523(2023)05-0052-09d o i :10.19509/j .c n k i .d z k q.t b 20220015 开放科学(资源服务)标识码(O S I D ):S u p p o r t i n g e f f e c t a n a l y s i s o f t h e H u a n gn i b i a n l a n d s l i d e b a s e d o n M a t D E MW a n g Fe i ,X u C h u ,L i W e i (F a c u l t y of E ng i n e e r i n g ,Chi n a U n i v e r s i t y of G e o s c i e n c e s (W u h a n ),W u h a n 430074,C h i n a )A b s t r a c t :[O b je c t i v e ]T h e H u a n g n i b i a n l a n d s l i d e i s a s m a l l -t o m e d i u m -s i z e d c o h e s i v e s o i l l a n d s l i d e t h a t h a s s h o w n s i g n s of m o v e m e n t .I t i s u rg e n t t o i m p l e m e n t s u p po r t a n d r e m e d i a t i o n m e a s u r e s f o r t h i s l a n d -s l i d e .T h i s s t u d y a i m s a t a n a l y s i n g t h e f a i l u r e m e c h a n i s m o f t h e l a n d s l i d e ,v e r i f y i n g t h e s t a b i l i t y r e qu i r e -m e n t s o f e x i s t i n g s u p p o r t s t r u c t u r e s ,a n d e x p l o r i n g r e a s o n a b l e a n d e f f e c t i v e s u p po r t l o c a t i o n s a n d f o r m s .[M e t h o d s ]T h i s p a p e r e s t a b l i s h e s a t h r e e -d i m e n s i o n a l m o d e l o f t h e H u a n gn i b i a n l a n d s l i d e u n d e r t h e c o n d i -t i o n o f a o n c e -i n -20-ye a rf l o o d l e v e l b a s e d o n M a t D E M.T h e d e f o r m a t i o n o f t h e l a n d s l i d e u n d e r t h r e e c o n -d i t i o n s o f n o s u p p o r t ,r e t a i n i ng w a l l s u p p o r t ,a n d th e c o m bi n a t i o n o f r e t a i n i n g w a l l a n d a n t i -s l i d e p i l e s u p-p o r t w e r e s i m u l a t e d ,a n d t h e n t h e s u p p o r t i n g e f f e c t o f d i f f e r e n t s u p p o r t i n g s c h e m e s a c c o r d i n g to t h e d i s -p l a c e m e n t f i e l d w a s a n a l ys e d .[R e s u l t s ]T h e s i m u l a t i o n r e s u l t s s h o w t h a t :①T h e l a n d s l i d e i s u n s t a b l e u n -d e r t h e o n c e -i n -20-y e a r f l o o d l e v e l ;②T h e s l i d i n g b o d y i s c u t f r o m t h e s l o pe f o o t o f B i n h e R o a d i n t h e c a s e o f n o s u p p o r t a n d f r o m t h e s l o p e f o o t o f t h e E x p r e s s w a y i n t h e c a s e o f r e t a i n i n g w a l l s u p p o r t ,i n d i c a t i n gCopyright ©博看网. All Rights Reserved.第5期王飞等:基于M a t D E M的黄泥扁滑坡支护效果分析t h a t t h e e x i s t i n g r e t a i n i n g w a l l s u p p o r t s t r u c t u r e d o e s n o t m e e t t h e s t a b i l i t y r e q u i r e m e n t s o f t h e l a n d s l i d e, a n d t h e r e f o r e t h e s l o p e f o o t o f t h e E x p r e s s w a y n e e d s t o b e m a i n l y s u p p o r t e d;a n d③T h e s c h e m e t h a t a r-r a n g e s t h e a n t i-s l i d e p i l e s a t t h e s l o p e f o o t o f t h e E x p r e s s w a y a n d a d d s a r e t a i n i n g p l a t e b e t w e e n t h e a n t i-s l i d e p i l e s i s e f f e c t i v e t o s u p p o r t t h e l a n d s l i d e a n d p r e v e n t s o i l b e t w e e n p i l e s f r o m s l i p p i n g.[C o n c l u s i o n] T h e a n a l y s i s r e s u l t s c a n p r o v i d e r e f e r e n c e s f o r s u b s e q u e n t e n g i n e e r i n g d e s i g n a n d c o n s t r u c t i o n a n d v e r i f y t h e a p p l i c a b i l i t y o f M a t D E M i n l a n d s l i d e a n a l y s i s.K e y w o r d s:M a t D E M;d i s c r e t e e l e m e n t s i m u l a t i o n;r e t a i n i n g w a l l;a n t i-s l i d e p i l e;l a n d s l i d e s u p p o r t e f f e c t;H u a n g n i b i a n l a n d s l i d eR e c e i v e d:2022-01-11;R e v i s e d:2022-03-24;A c c e p t e d:2022-04-18滑坡的支护结构比较多样,有抗滑桩㊁挡土墙㊁锚索㊁柔性结构㊁生态结构等[1],在实际工程中,往往需要结合滑坡的诱发因素㊁变形破坏机制㊁几何边界条件㊁现状等情况[2],以及便捷性㊁安全性㊁经济性㊁技术性等原则[3]来综合考虑滑坡的支护形式㊂数值模拟技术可以很好地结合多种因素对某一地质问题进行具体分析,是一种解决地质问题的有力工具[4]㊂随着计算机性能的飞速发展,数值模拟技术在地质灾害分析中的应用越来越广泛,运用数值模拟技术分析滑坡在不同支护形式下的支护效果,可以为实际工程的设计㊁施工提供一定参考,国内外许多学者为此开展了大量研究㊂宁宇等[5]利用有限差分法和自主研发的程序分析了 联合h型桩 在某滑坡中的支护效果,并与其他支护形式进行了对比,验证了 联合h型桩 更有利于提高滑坡的稳定性;陶连金等[6]运用离散元法和U D E C软件来模拟边坡的开挖和治理,通过对比不同支护情况下的边坡位移云图,证明了锚索抗滑桩的支护效果更好;T r o j n a r[7]运用有限元法对比分析了石笼墙和抗滑桩对某滑坡的支护效果,验证了抗滑桩的支护比石笼墙更经济㊁有效;K a h y a o g l u等[8]运用有限元法模拟了双排桩对某失稳滑坡的支护效果,并分析了桩身实测位移与模拟位移之间的差异,验证了双排桩支护的有效性㊂岩土体在微观上的本质是由一系列颗粒㊁胶结物及孔隙所组成的结构体[9],许多学者认为运用离散元法能够有效地模拟岩土体的非连续性㊁不均匀性和变形破坏较大等特点[10-15]㊂基于此,本研究运用离散元软件M a t D E M进行数值模拟,该软件由南京大学刘春等[16]研发,是一款默认采用线弹性理论为接触模型的离散元模拟软件,其结合了矩阵运算和高性能G P U计算,具有强大的运算能力,目前许多国内学者运用该软件对滑坡进行了一系列研究㊂X u e等[17]基于M a t D E M分析了大型土石混合体边坡的稳定性;朱晨光等[18]基于M a t D E M分析了高速远程滑坡在滑动过程中滑带的生热规律和温度分布特征;L u o等[19]基于M a t D E M模拟了 贵州纳雍岩崩 的破坏过程;G i a n v i t o等[20]基于M a t D E M分析了已失稳的滑坡将可能产生的新破坏;汪志林等[21]基于M a t D E M分析了含双裂隙组灰岩边坡的渐进失稳过程;栗晓松等[22]基于M a t D E M分析了烟家沟滑坡的演化过程㊂前人运用M a t D E M对滑坡进行的一系列研究,均验证了M a t D E M在滑坡分析中的可靠性和适用性,但将其运用于滑坡治理方面的研究相对较少,因此,笔者在前人的基础上,运用该软件对黄泥扁滑坡进行数值模拟,对比不同支护方案在20年一遇洪水工况下的支护效果,并提出能够满足滑坡稳定性要求的支护方案㊂1滑坡概况根据勘察成果可知:黄泥扁滑坡位于四川省达州市通川区滨水公园中部地带,具体位于滨河路与快速通道之间,为州河的右岸边坡(图1-a);研究区域钻孔揭露的地层有:第四系全新统人工填土层(Q m l)㊁第四系全新统残坡积粉质黏土层(Q e l+d l4),中侏罗统上沙溪庙组(J2s)(强风化和中风化)砂质泥岩层,其中,滑体由人工填土层和粉质黏土层组成,滑带为强风化砂质泥岩层(图1-b)㊂该地区的构造简单,区内无断层,影响该滑坡稳定性的因素主要是水文地质条件(疏松多孔的滑土体,在雨水的渗流作用下,易处于饱和状态,其抗剪强度大大降低)和人类工程活动(快速通道㊁滨河路的建设)㊂勘察结果显示:该段滑坡主要危及对象为坡体上的快速通道和滨河路,快速通道安全等级为一级,根据现行的‘公路滑坡防治设计规范“可确定该滑坡的防治安全等级为一级㊁安全系数(非正常工况)为1.20㊂根据实地调查发现:该滑坡已出现明显滑移迹象,滑坡上所建的两栋建筑物(已拆,见图1-a)出现开裂现象,裂缝宽5~20mm(图2)㊂由此判断该滑坡处于强变形阶段,主滑方向172ʎ,东西宽约80 m,南北长约100m,滑体最厚处达20m,总体积约3.2万m3,为小型中层土质滑坡㊂35Copyright©博看网. All Rights Reserved.h t t p s ://d z k j q b .c u g.e d u .c n 地质科技通报2023年图1 黄泥扁滑坡工程地质平面图(a )和剖面图(b)F i g .1 E n g i n e e r i n g g e o l o g i c a l p l a n e (a )a n d p r o f i l e (b )o f t h e H u a n gn i b i a n l a n d s l i de 图2 建筑物墙体开裂F i g .2 B u i l d i n g wa l l c r a c k s 2 模型建立及参数确定2.1模型的建立利用M a t D E M 构建滑坡的基本模型时需要分两步进行:①构造初始堆积模型㊂构建一定尺寸的空模型箱,往模型箱内生成一定半径的颗粒,通过压力板压实颗粒来模拟重力沉积;本研究依据勘察得到的滑坡边界,确定了模型箱的长㊁宽㊁高分别为115,85,70m ,并根据计算机的计算能力以及模型的精确程度,设置了单元的平均半径为0.7m ㊂②切割地层并设置材料参数㊂导入真实的三维地层数据去切割上一步生成的初始堆积模型,并将材料参数赋予各层单元㊂本研究依据勘察报告中给出的地形平面图和地层剖面图,提取出每个层面的高程坐标数据集,经切割后的模型见图3,材料参数的设置见第2.2节㊂图3中显示了4种颜色的材料,从上往下分别代表:人工填土㊁粉质黏土㊁强风化砂质泥岩㊁中风化砂质泥岩,坡脚处为州河,坡体上有2条重要的道路(滨河路和快速通道),水位线为20年一遇的洪水位界线㊂2.2参数的确定在模拟过程中需要先确定岩土体的杨氏模量㊁泊松比㊁抗拉强度㊁抗压强度㊁内摩擦系数以及密度45Copyright ©博看网. All Rights Reserved.第5期王 飞等:基于M a t D E M 的黄泥扁滑坡支护效果分析①为人工填土;②为粉质黏土;③为强风化砂质泥岩;④为中风化砂质泥岩图3 滑坡的基本模型F i g.3 B a s i c m o d e l o f t h e l a n d s l i d e 等宏观力学参数㊂本研究通过参考勘察资料,结合同类土体的力学参数取值来确定所选用的参数(表1)㊂在颗粒随机堆积的模型中,其宏观力学参数通常比设定值小,为了获得力学性质更准确的材料,需要对材料进行一定的训练,训练的原理及过程如下:①输入材料属性来建立一个颗粒随机堆积的块体模型;②对这个块体模型进行单轴压缩㊁剪切等一系列模拟实验,得到模型各参数的模拟值;③利用模拟值和设定值之间的比率去调整设定值;④重复以上步骤,直至各力学参数收敛于设定值㊂由于颗粒在堆积㊁胶结后,模型的力学性质主要受单元的接触模型及其参数控制,所以需要将经过训练的宏观参数转换为适用于紧密规则堆积的离散元模型,宏㊁微观力学参数之间的解析解如下[23]:K n =2E d 4(1-2ν)(1)K s =2(1-5ν)E d4(1+ν)(1-2ν)(2)X b =3K n +K s62K n (K n +K s )㊃T u ㊃d 2(3)F s 0=1-2μp 6㊃C u ㊃d 2(4)μp =-22+2[(1+μ2i )+μi ]22+2[(1+μ2i )+μi ]2(5)上述公式中:法向刚度(K n )㊁切向刚度(K s )㊁断裂位移(X b )㊁初始抗剪力(F s 0)㊁摩擦系数(μp )均为材料的5种微观力学参数;杨氏模量(E )㊁泊松比(ν)㊁抗压强度(C u )㊁抗拉强度(T u )㊁内摩擦系数(μi )均为材料的5种宏观力学参数;d 为颗粒直径㊂经过转换后的各材料微观参数见表2㊂由公式(2)可知:只有当泊松比ɤ0.2时,才能保证切向刚度为正数㊂这是因为:上述转换公式适用于三维紧密规则堆积的离散元模型,其最大可能表1 各材料的宏观力学参数取值T a b l e 1 M a c r o s c o p i c m e c h a n i c a l pa r a m e t e r s o f t h e m a t e r i a l s 状态岩性杨氏模量E /M P a 泊松比ν抗拉强度T u/M P a 抗压强度C u /M P a 内摩擦系数μi密度ρ/(k g㊃m -3)天然饱和①人工填土50.170.0020.20.462110②粉质黏土100.180.0040.40.222060③1强风化砂质泥岩2000.102.5006.20.502490③2中风化砂质泥岩4000.155.00015.40.602590挡土墙和抗滑桩200000.1510.00050.00.602850①人工填土30.180.0010.10.362130②粉质黏土80.190.0020.20.172090③1强风化砂质泥岩1000.121.2003.10.302510③2中风化砂质泥岩3000.162.5007.40.482600表2 各材料的微观力学参数取值T a b l e 2 M i c r o s c o p i c m e c h a n i c a l pa r a m e t e r s o f t h e m a t e r i a l s 状态岩性法向刚度K n /(MN ㊃m -1)切向刚度K s /(MN ㊃m -1)断裂位移Xb /mm 初始抗剪力F s 0/MN摩擦系数μp 天然饱和①人工填土1.90.20.300.20.46②粉质黏土3.50.30.400.60.22③1强风化砂质泥岩75.037.30.503.10.50③2中风化砂质泥岩147.559.60.805.00.60挡土墙和抗滑桩4630.0195.00.082.80.60①人工填土1.30.20.300.10.36②粉质黏土3.10.30.400.40.17③1强风化砂质泥岩56.315.50.402.40.30③2中风化砂质泥岩134.123.20.704.50.4855Copyright ©博看网. All Rights Reserved.h t t p s ://d z k j q b .c u g.e d u .c n 地质科技通报2023年泊松比为0.2[24],对于泊松比较大的材料,其泊松比通常来自于其特定的结构,需要通过特定的堆积来形成较大的泊松比㊂比如:当直接设置不规则堆积体的材料泊松比为0.2时(本文的模型就是不规则堆积体),实际得到的泊松比会大于该值,具体的值需要对材料进行训练才能确定㊂3 数值模拟3.1无任何支护本研究采用等效模拟的方法来模拟水对土体以及整个滑坡的影响㊂将处于20年一遇洪水位以下的颗粒筛选出来,假设这些土体已达到饱和状态,将表2中的参数赋给它们,并赋予它们一个竖直向上的浮力来减弱它们自身的体力㊂整个模拟从开始破坏到趋于稳定,共历时约39s ,图4中分别截取了第1.5,3,6,9,18,39s 时的位移云图,可以看出:在第1.5s 时(图4-a),滑坡的中上部出现一条明显的裂缝,并向下延伸成滑面,此裂缝出现的位置距坡顶约40m ,与勘察过程中发现的地裂缝(图1)所在位置基本吻合;到第3s 时(图4-b),裂缝和滑面沿着图中红色箭头所指的方向进一步延伸;到第6s 时(图4-c ),坡体上的裂缝进一步扩大,滑面贯通,整个滑体从滨河路的坡脚处剪出;到第9s 时(图4-d),坡体已经完全破坏,裂缝以上的单元也开始出现滑落;到第18s 时(图4-e),破坏几乎曼延至坡顶;到第39s 时(图4-f),整个系统才趋于稳定㊂基于此,可分析出该滑坡在无支护工况下的破坏机理为:水位的上升降低了滑带土和滑体土的强度,同时水的浮力减小了主滑体的抗滑力,继而引发了拉张裂缝的形成,随后裂缝逐渐扩展至滑面,直至滑面贯通,最终造成滑坡的整体破坏㊂图4 无支护时的滑坡位移云图F i g .4 L a n d s l i d e d i s p l a c e m e n t d i a g r a m o f t h e c a s e w i t h o u t s u p po r t 此外,本研究采用G e o s t u d i o 软件计算了图3所示剖面的稳定性系数(图5),计算结果从以下两个方面验证了M a t D E M 模拟结果的可靠性:①G e o s t u d i o 计算得到的稳定性系数为0.86,表明该滑坡在无支护工况下是极其不稳定的;②图5所示的滑体和滑面与图4在第6s 显示的滑体和滑面大致相同㊂3.2挡土墙支护实际的工程中,在滨河路的坡脚处已经建有一座挡土墙,本研究根据实际的情况,在模型同样的位置也构建了一座挡土墙(图6-c),以分析现有的支护形式是否满足滑坡的稳定性要求㊂挡土墙的截面形图5 无支护时的G e o s t u d i o 计算结果F i g.5 G e o s t u d i o c a l c u l a t i o n r e s u l t s o f t h e c a s e w i t h -o u t s u p po r t 65Copyright ©博看网. All Rights Reserved.第5期王 飞等:基于M a t D E M 的黄泥扁滑坡支护效果分析状及尺寸见图6-a ,整段挡土墙长85m ,高13m (图6-b)㊂整个模拟从开始破坏到趋于稳定,共历时约42s ,图7中分别截取了第1.2,3,6,9,18,42s 时的位移云图,可以看出:在第1.2s 时(图7-a),滑体的中上部产生了轻微的变形;在第3s 时(图7-b ),变形已经十分明显;在第6s 时(图7-c),滑面已经基本贯通,但是由于挡土墙的存在,阻挡了整个滑体从滨河路的坡脚处剪出,所以滑体只能从快速通道的坡脚处剪出;在第9s 时(图7-d),剪出的土体几乎掩埋了整条滨河路,上部的破坏继续向坡顶蔓延;在第18s 时(图7-e),已经有大量的土颗粒越过了挡土墙;在第42s 时(图7-f),更多的土颗粒越过挡土墙堆积在河道中,此时整个系统趋于稳定㊂基于此,可分析出该滑坡在挡土墙支护工况下的破坏机理与无支护工况下的破坏机理是类似的,都是由于土体的强度和滑体的抗滑力降低而导致的完全失稳,而挡土墙的存在仅仅是改变了剪出口的位置㊂所以,在无支护工况和挡土墙支护工况下的相同位置都出现了地裂缝㊂a .挡土墙截面尺寸(单位:m );b .挡土墙模型(单位:m );c .挡土墙支护的滑坡模型图6 挡土墙结构示意图F i g .6 S t r u c t u r e d i a g r a m o f t h e r e t a i n i n g wa ll 图7 挡土墙支护时的滑坡位移云图F i g .7 L a n d s l i d e d i s p l a c e m e n t d i a g r a m o f t h e c a s e w i t h r e t a i n i n g w a l l s u p po r t 对比无支护时的位移云图(图8-a ,b ),可以发现,在无支护时,滑体从滨河路的坡脚处剪出,而有挡土墙支护时,滑体从快速通道的坡脚处剪出,所以,为了提高滑坡的整体稳定性,还需要对快速通道的坡脚进行支护㊂从图8-c 可以看出:挡土墙在一定程度上消减了土颗粒的移动速度,并阻挡了一部分滑体滑入州河中㊂同种工况下,G e o s t u d i o 计算得到挡土墙支护75Copyright ©博看网. All Rights Reserved.h t t p s ://d z k j q b .c u g.e d u .c n 地质科技通报2023年a .无支护时的剪出口位置;b .挡土墙支护时的剪出口位置;c .滑体平均速度和平均位移曲线图8 无支护和挡土墙支护下的模拟结果对比图F i g .8 C o m p a r i s o n o f s i m u l a t i o n r e s u l t s b e t w e e n t h e c a s e u n d e r n o s u p p o r t a n d t h a t w i t h r e t a i n i n g w a l l s u p po r t 时的稳定性系数为1.03(图9),根据‘公路滑坡防治设计规范“,该滑坡在挡土墙支护的情况下仍处于欠稳定状态,图9还显示了滑体从快速通道的坡脚处剪出,进一步说明了对快速通道坡脚处进行支护的必要性㊂图9 挡土墙支护时的G e o s t u d i o 计算结果F i g.9 G e o s t u d i o c a l c u l a t i o n r e s u l t s o f t h e c a s e w i t h r e t a i n i n g w a l l s u p po r t 3.3挡土墙和抗滑桩共同支护快速通道坡脚处的滑体厚度约16m ,厚度较大,不宜用挡土墙进行支护㊂抗滑桩适用于滑动面埋深较大㊁滑坡推力较大的滑坡体,因此,选用抗滑桩对该滑坡进行支护是比较合理的㊂抗滑桩的设置参照‘公路滑坡防治设计规范“,形式选用矩形人工挖孔桩,桩身截面长度3m ,截面宽度2m ,桩长24m ,嵌固段长度8m (图10),桩间距为5倍,共9根,桩体的材料参数和挡土墙一致㊂整个模拟过程共历时9s ,从位移云图(图11)可以看出:在0.6s 时,坡体发生了轻微的变形(图11-a);在6s 时,大量桩间土从桩间挤出,主滑体的位移明显增大(图11-b );在9s 时,桩间土继续向外挤出,但主滑体的位移相对于第6s 并没有增大很多(图11-c)㊂模拟结果表明:抗滑桩对整个滑坡起到了很好的支护效果,但对快速通道和滨河路的支护并不理想,大量桩间土的滑落一方面会造成快速通道结构的破坏,另一方面会掩图10 抗滑桩位置示意图F i g .10 D i s t r i b u t i o n m a p of t h e a n t i -s l i d e p i l e 埋滨河路㊂所以,将桩间距设置为5倍,或者仅采用抗滑桩对快速通道坡脚进行支护的方案还是不够的㊂为了有效解决土体从桩间滑落这一问题,可以适当缩小桩间距,但鉴于以上的模拟结果来看,9根桩已经在很大程度上提高了滑坡的整体稳定性,此时再缩小桩间距来增加抗滑桩,会大大增加工程造价,此方案在经济上是不合理的㊂另外可以在抗滑桩之间加设一道挡土板,此方案不仅工程造价低,而且施工工艺也较埋置抗滑桩简单㊂挡土板高6m ,厚0.5m ,材料参数和挡土墙一致㊂该方案的模拟结果如图12所示,对比图11可以发现:在仅有抗滑桩支护时,上部滑体的平均位移在第6s 达到约2m ,而在抗滑桩 挡土板共同支护时,同时刻滑体的平均位移不到0.5m ,表明挡土板的存在,不仅有效地解决了土体从桩间滑落的问题,也在一定程度上改善了上部主滑体的稳定性㊂同种工况下,G e o s t u d i o 计算得到的稳定性系数为1.27(图13),根据‘公路滑坡防治设计规范“,该滑坡在挡土墙和抗滑桩共同支护的情况下处于稳定状态,进一步说明了M a t D E M 模拟结果的可靠性㊂值得一提的是:根据规范,在计算滑坡的稳定性系数时,需要考虑地震工况,但鉴于本研究用G e o s -t u d i o 计算稳定性系数的主要目的是验证M a t D E M模拟结果的合理性,所以G e o s t u d i o 的计算没有考虑85Copyright ©博看网. All Rights Reserved.第5期王 飞等:基于M a t D E M的黄泥扁滑坡支护效果分析图11 抗滑桩支护时的滑坡位移云图F i g .11 L a n d s l i d e d i s p l a c e m e n t d i a g r a m o f t h e c a s e w i t h a n t i -s l i d e p i l e s s u p po r t图12 抗滑桩-挡土板支护时的滑坡位移云图F i g .12 L a n d s l i d e d i s p l a c e m e n t d i a g r a m o f t h e c a s e w i t h a n t i -s l i d e p i l e s a n d r e t a i n i n g p l a t e s u p po rt 图13 抗滑桩支护时的G e o s t u d i o 计算结果F i g.13 G e o s t u d i o c a l c u l a t i o n r e s u l t s o f t h e c a s e w i t h t h e a n t i -s l i d e p i l e s u p po r t 地震工况㊂4 结 论(1)在20年一遇洪水且无任何支护的情况下,黄泥扁滑坡处于不稳定状态,其破坏在最开始形成一条拉张裂缝,随后这条裂缝逐渐扩展至滑面,直至滑面贯通,最后造成滑坡的整体破坏㊂(2)在20年一遇的洪水工况下,实际工程中修建于滨河路坡脚处的挡土墙并不能阻挡滑坡的整体破坏,对提高滑坡整体稳定性的作用并不大㊂在无支护和挡土墙支护的模拟中,裂缝出现的位置与实际勘察发现的地裂缝位置大致相同,证明了模拟结果的可靠性㊂(3)在无任何支护的情况下,失稳的滑坡从滨河路的坡脚处剪出;在挡土墙支护的情况下,失稳的滑坡从快速通道的坡脚处剪出,所以,对快速通道的坡脚进行支护是必要的㊂(4)将抗滑桩布置于快速通道的坡脚处,并采用5倍桩间距,此方案可以阻止滑坡的整体失稳,在很大程度上提高了滑坡的整体稳定性;但是此方案存在土体从桩间滑落的隐患,该隐患会威胁到两条道路的正常运行㊂为了解决此隐患,可以在抗滑桩之间加设一道挡土板㊂(所有作者声明不存在利益冲突)参考文献:[1] 罗丽娟,赵法锁.滑坡防治工程措施研究现状与应用综述[J ].自然灾害学报,2009,18(4):158-164.L u o L J ,Z h a o F S .S t a t u s o f r e s e a r c h a n d a p p l i c a t i o n o f e n gi -n e e r i n g m e a s u r e s f o r p r e v e n t i n g a n d c o n t r o l l i n g la n d s l i d e [J ].J o u r n a l o f N a t u r a l D i s a s t e r s ,2009,18(4):158-164(i n C h i n e s e w i t h E n gl i s h a b s t r a c t ).[2] 张倬元.滑坡防治工程的现状与发展展望[J ].地质灾害与环境保护,2000,11(2):89-97,181.Z h a n g Z Y.T h e p r e s e n t s t a t u s ,t e c h n i c a l a d v a n c e a n d d e v e l o p-m e n t t r e n d s o f l a n d s l i d e r e m e d i a l m e a s u r e s [J ].J o u r n a l o f G e o -l o gi c a l H a z a r d s a n d E n v i r o n m e n t P r e s e r v a t i o n ,2000,11(2):89-97,181(i n C h i n e s e w i t h E n gl i s h a b s t r a c t ).[3] 魏金花.地质灾害治理工程施工中边坡稳定问题及滑坡治理方法[J ].西部资源,2020(3):28-30.95Copyright ©博看网. All Rights Reserved.h t t p s://d z k j q b.c u g.e d u.c n地质科技通报2023年W e i J H.S l o p e s t a b i l i t y p r o b l e m a n d l a n d s l i d e t r e a t m e n t m e t h-o d i n c o n s t r u c t i o n o f g e o l o g i c a l h a z a r d c o n t r o l p r o j e c t[J].W e s t e r n R e s o u r c e,2020(3):28-30(i n C h i n e s e w i t h E n g l i s h a b-s t r a c t).[4]黄润秋,张倬元,王士天.工程地质学中的数值模拟与拟合[J].水文地质工程地质,1991(5):7-10.H u a n g R Q,Z h a n g Z Y,W a n g S T.N u m e r i c a l s i m u l a t i o n a n df i t t i ng i n e n g i n e e r i n g g e o l o g y[J].H y d r o g e o l o g y a n d E n g i n e e r-i n g G e o l o g y,1991(5):7-10(i n C h i n e s e w i t h E n g l i s h a b s t r a c t).[5]宁宇,黄青富,郝李坤,等.联合h型桩在滑坡体阻滑中应用数值模拟研究[J].科学技术与工程,2021,21(23):10004-10012.N i n g Y,H u a n g Q F,H e i L K,e t a l.N u m e r i c a l s i m u l a t i o n r e-s e a r c h o n a p p l i c a t i o n o f c o m b i n e d h-p i l e i n s l i d i n g r e s i s t a n c e o f l a n d s l i d e[J].S c i e n c e T e c h n o l o g y a n d E n g i n e e r i n g,2021,21(23):10004-10012(i n C h i n e s e w i t h E n g l i s h a b s t r a c t). [6]陶连金,沈小辉,王开源,等.某大型高速公路滑坡稳定性分析及锚桩加固的模拟研究[J].工程地质学报,2012,20(2):259-265.T a o L J,S h e n X H,W a n g K Y,e t a l.S t a b i l i t y a n a l y s i s a n d s u p p o r t i n g e f f e c t m o d e l i n g o f a l a r g e s c a l e l a n d s l i d e o n h i g h-w a y[J].J o u r n a l o f E n g i n e e r i n g G e o l o g y,2012,20(2):259-265(i n C h i n e s e w i t h E n g l i s h a b s t r a c t).[7] T r o j n a r K.N u m e r i c a l a n a l y s i s o f t h e l a n d s l i d e g e o h a z a r d s:C a s e s t u d y w i t h g a b i o n s a n d p i l e s s o l u t i o n s[C]ʊA n o n.P r o-c e ed i n g s o f C E E2019.C h a m:S p r i n ge r I n t e r n a t i o n a l P u b l i s h-i n g,2020:474-479.[8] K a h y a o g l u M R,I m a n c h G,O z d e n G.N u m e r i c a l s i m u l a t i o n s o fl a n d s l i d e-s t a b i l i z i n g p i l e s:A r e m e d i a t i o n p r o j e c t i n S o k e,T u r-k e y[J].E n v i r o n m e n t a l E a r t h S c i e n c e s,2017,76(19):656.[9]吴恒,张信贵,易念平,等.水土作用与土体细观结构研究[J].岩石力学与工程学报,2000,19(2):199-204.W u H,Z h a n g X G,Y i N P,e t a l.R e s e a r c h o n w a t e r-s o i l i n t e r-a c t i o n a n d m e s o s t r u c t u r e o f s o i l m a s s[J].C h i n e s e J o u r n a l o fR o c k M e c h a n i c s a n d E n g i n e e r i n g,2000,19(2):199-204(i n C h i-n e s e w i t h E n g l i s h a b s t r a c t).[10]C u n d a l l P A,S t r a c k O D L.A d i s c r e t e n u m e r i c a l m o d e l f o rg r a n u l a r a s s e m b l i e s[J].G e o t e c h n i q u e,1980,30(3):331-336.[11]王泳嘉,邢纪波.离散单元法同拉格朗日元法及其在岩土力学中的应用[J].岩土力学,1995,16(2):1-14.W a n g Y J,X i n g J B.D i s c r e t e e l e m e n t m e t h o d a n d l a g r a n g i a ne l e m e n t m e t h o d a n d t h e i r a p p l i c a t i o n s i n g e o m e c h a n i c s[J].R o c k a n d S o i l M e c h a n i c s,1995,16(2):1-14(i n C h i n e s e w i t hE n g l i s h a b s t r a c t).[12]张泽林,吴树仁,唐辉明,等.反倾岩质边坡的时效变形破坏研究[J].地质科技情报,2014,33(5):181-187.Z h a n g Z L,W u S R,T a n g H M,e t a l.T i m e d e p e n d e n t d e f o r m-a t i o n o f a n t i t h e t i c d i p r o c k s l o p e[J].G e o l o g i c a l S c i e n c e a n dT e c h n o l o g y I n f o r m a t i o n,2014,33(5):181-187(i n C h i n e s e w i t hE n g l i s h a b s t r a c t).[13]张江晖,徐守余,蒋静,等.含夹层碳酸盐岩储层裂缝发育规律研究[J].地质科技情报,2019,38(2):75-80.Z h a n g J H,X u S Y,J i a n g J,e t a l.F r a c t u r e d e v e l o p m e n t i n t h ec a r b o n a t e r e s e r v o i r w i t h i n t e r l a y e r s[J].G e o l o g i c a l S c i e n c e a n dT e c h n o l o g y I n f o r m a t i o n,2019,38(2):75-80(i n C h i n e s e w i t hE n g l i s h a b s t r a c t).[14]井旭,谢婉丽,单帅.原状及重塑黄土双轴试验微观力学特征离散元模拟[J].地质科技通报,2021,40(3):184-193.J i n g X,X i e W L,S h a n S.D i s c r e t e e l e m e n t s i m u l a t i o n s t u d y o n m i c r o m e c h a n i c a l c h a r a c t e r i s t i c s o f u n d i s t u r b e d a n d r e m o l d e d l o e s s i n b i a x i a l t e s t[J].B u l l e t i n o f G e o l o g i c a l S c i e n c e a n d T e c h n o l o g y,2021,40(3):184-193(i n C h i n e s e w i t h E n g l i s h a b-s t r a c t).[15]陆敏凤,唐朝晖,柴波,等.矿渣类颗粒介质结构对力链发展规律的影响[J].地质科技通报,2022,41(4):274-281.L u M F,T a n g Z H,C h a i B,e t a l.G r a n u l e s t r u c t u r e i n f l u e n c e s o n f o r c e c h a i n d e v e l o p m e n t i n t h e s l a g m e d i u m[J].B u l l e t i n o fG e o l o g i c a l S c i e n c e a n d T e c h n o l o g y,2022,41(4):274-281.[16]刘春,施斌,顾凯,等.岩土体大型三维离散元模拟系统的研发与应用[C]ʊ佚名.2014年全国工程地质学术大会论文集.[出版地不详]:[出版社不详],2014:560-566.L i u C,S h i B,G u K,e t a l.D e v e l o p m e n t a n d a p p l i c a t i o n o f l a r g e-s c a l e d i s c r e t e e l e m e n t s i m u l a t i o n s y s t e m f o r r o c k a n d s o i l[C]ʊA n o n.P a p e r s o f2014N a t i o n a l E n g i n e e r i n g G e o l o g y A c a d e m i cC o n f e r e n c e.[S.l.]:[s.n.],2014:560-566(i n C h i n e s e w i t hE n g l i s h a b s t r a c t).[17]X u e Y D,Z h o u J,H u a n g H W,e t a l.A n a l y s i s o f l a r g e s o i l r o c km i x t u r e s l o p e b a s e d o n D E M[J].I n t e r n a t i o n a l S o c i e t y f o r R o c k M e c h a n i c s a n d R o c k E n g i n e e r i n g,2019,27(3):651-658.[18]朱晨光,刘春,许强,等.滑坡滑带摩擦热离散元数值模拟研究[J].工程地质学报,2019,27(3):651-658.Z h u C G,L i u C,X u Q,e t a l.D i s c r e t e e l e m e n t n u m e r i c a l s i m u-l a t i o n r e s e a r c h o n f r i c t i o n h e a t i n s l i d i n g z o n e o f t h e l a n d s l i d e [J].J o u r n a l o f E n g i n e e r i n g G e o l o g y,2019,27(3):651-658(i nC h i n e s e w i t h E n g l i s h a b s t r a c t).[19]L u o H,X i n g A G,J i n K P,e t a l.D i s c r e t e e l e m e n t m o d e l i n g o ft h e N a y o n g r o c k a v a l a n c h e,G u i z h o u,C h i n a c o n s t r a i n e d b y d y-n a m i c p a r a m e t e r s f r o m s e i s m i c s i g n a l i n v e r s i o n[J].R o c k M e-c h a n i c s a nd R o c k E n g i ne e r i n g,2021,15(7):1359-1375.[20]G i a n v i t o S,F a n X M,X u Q,e t a l.S o m e c o n s i d e r a t i o n s o n t h eu s e o f n u m e r i c a l m e t h o d s t o s i m u l a t e p a s t l a n d s l i d e s a n d p o s s i-b l e n e w f a i l u r e s:T h ec a s e o f t h e r e c e n t X i n m o l a nd s l i d e(S i-c h u a n,C h i n a)[J].L a nd s l i de s,2018,15(7):1359-1375.[21]汪志林,叶海旺,李子旋,等.含双裂隙组灰岩边坡渐进失稳过程分析[J].矿冶工程,2021,41(2):20-23,32.W a n g Z L,Y e H W,L i Z X,e t a l.A n a l y s i s o f p r o g r e s s i v e i n s t a-b i l i t y p r oc e s s o f l i m e s t o n e s l o p e w i t h t w o g r o u p s o f f i s s u r e s[J].M i n i n g a n d M e t a l l u r g i c a l E n g i n e e r i n g,2021,41(2):20-23,32(i n C h i n e s e w i t h E n g l i s h a b s t r a c t).[22]栗晓松,范文,曹琰波,等.基于M a t D E M的烟家沟滑坡演化过程数值模拟分析[J].地质与资源,2021,30(2):199-206.L i X S,F a n W,C a o Y B,e t a l.M a t d e m-b a s e d n u m e r i c a l s i m u-l a t i o n a n a l y s i s o f Y a n j i a g o u l a n d s l i d e e v o l u t i o n p r o c e s s[J].G e-o l o g y a n d R e s o u r c e s,2021,30(2):199-206(i n C h i n e s e w i t hE n g l i s h a b s t r a c t).[23]L i u C,X u Q,S h i B,e t a l.M e c h a n i c a l p r o p e r t i e s a n d e n e r g yc o n v e r s i o n o f3D c l o s e-p a c k ed l a t t i ce m o d e lf o r b r i t t l e r o c k s[J].C o m p u t e r s&G e o s c i e n c e s,2017,103:12-20. [24]刘春.地质与岩土工程矩阵离散元分析[M].北京:科学出版社,2019.L i u C.M a t r i x d i s c r e t e e l e m e n t a n a l y s i s o f g e o l o g y a n d g e o t e c h-n i c a l e n g i n e e r i n g[M].B e i j i n g:S c i e n c e P u b l i s h i n g,2019(i nC h i n e s e).06Copyright©博看网. All Rights Reserved.。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
泊松比的定义
法国力学家西莫恩·德尼·泊松( Simeon Denis Poisson,1781-1840)于1829年发表《弹性体平衡和运动研究报告》,提出泊松比(Poisson's ratio)概念。
泊松比是指材料在单向受拉或受压时,横向正应变与轴向正应变的比值,也叫横向变形系数,它是反映材料横向变形的弹性常数。
材料沿载荷方向产生伸长(或缩短)变形的同时,在垂直于载荷的方向会产生缩短(或伸长)变形。
垂直方向上的应变ε
y 与载荷方向上的应变ε
x
之比的负值称为材料的泊松
比。
以μ表示泊松比,则μ=-εy/εx。
在材料弹性变形阶段内,μ是一个常数。
理论上,各向同性材料的三个弹性常数E、G、μ中,只有两个是独立的,
因为它们之间存在如下关系:
G=E/[2(1+μ)]
材料常数G是剪切模量(modulus of rigidity),是剪切应力与应变的比值,又称切变模量或刚性模量,材料的力学性能指标之一。
是材料在剪切应力作用下,在弹性变形比例极限范围内,切应力与切应变的比值。
它表征材料抵抗切应变的能力。
剪切模量越大,则表示材料的刚性越强。
剪切模量的倒数称为剪切柔量,是单位剪切力作用下发生切应变的量度,可表示材料剪切变形的难易程度。
泊松比试验测定方法
(1)机械法
包括引伸计法和电阻应变法。
测试原理:对试样施加轴向力,在其弹性范围内测定相应的轴向变形和横向变形,然后计算其泊松比。
试验时通过横向和纵向引伸计自动记录方法绘制横向—纵向应变曲线,得到材料的泊松比。
(2)声学方法
声学测试泊松比的方法通常是根据弹性波理论,通过测定纵横波速来推算材料的泊松比。
以超声脉冲回波法为例,采用超声波脉冲回波法通过测量弹性波在固体样品中得传播速度来获取材料的弹性参数。
弹性波在各向同性的弹性介质中传播时,根据弹性波的在固体中的传播理论,不同模式的声波在固体中得传播速度与材料的相应的弹性模量和密度相关。
该方法先根据性波的在固体中的传播理论公式测得玻璃的杨氏模量和剪切模量,然后通过公式计算得到泊松比。
除了超声波脉冲回波法以外,用于泊松比测试的声学方法还包括布里渊散射(SBS)、表面声波(AM)、声学显微(AM)等,由于非金属材料,通常其声阻和内阻尼都较大的,要测得准确的声速和振动测试都还比较困难,目前这几种方法都还处于研究阶段,主要用于金属材料的泊松比测试。
(3)光学方法
用于材料泊松比测试的光学方法以弹性力学中板或梁的纯弯曲理论为基础,从干涉图样中推算泊松比,现在最主要的测试方法有光干涉测量法、光导热塑全息照相法、数字散斑内相关法(DSCM)等。
光学方法是非接触式的测量方法,理论上来讲其测量精度优于传统的机械法,但要求照相、摄像器材及振动台具,且对光电测试法的分辨率要求较高,虽然这些方法国内外都已做出大量研究,但并没有形成统一的测试标准,实验所能达到的测量误差亦很难评估。
常用材料泊松比
泊松比为什么不大于0.5?
如下图所示的微元体,变形前长宽高分别为dx、dy、dz,体积为V=dxdydz。
受力变形后,三个边的长度分别变成:
d x+ε1d x=(1+ε1)d x
d y+ε2d y=(1+ε2)d y
d z+ε3d z=(1+ε3)d z
由上式得到变形后的体积V1:
V1=(1+ε1)(1−ε2)(1+ε3)d x d y d z
略去高阶量后简化为:
V1=(1+ε1+ε2+ε3)d x d y d z
微元体受力后的体积变化,又称为单位体积的体积改变量,或体应变,用θ表示:
θ=V1−V
V
=ε1+ε2+ε3
由广义胡克定律,在考虑泊松效应时,xyz三个方向都受力的情况下,ε1、ε2、ε3的计算公式如下:
ε1=1
E
[σ1−μ(σ2−σ3)]
ε2=1
E
[σ2−μ(σ1+σ3)]
ε3=1
E
[σ3−μ(σ1+σ3)]
代入到体积变化计算式中有:
θ=V1−V
V
=ε1+ε2+ε3=
1−2μ
E
(σ1+σ2+σ3)=
σm
K
K为体积模量,σm为静水压力。
K=
E
3(1−2μ)
σm=σ1+σ2+σ3
3
当材料为完全不可压缩时,体积应变θ=0,即:
σm/K=σm (1/k) =σm [3(1-2μ)/E]=0
此时μ=0.5。
泊松比最小值为-1
由杨氏模量、剪切模量和泊松比之间的关系式:G=E/[2(1+μ)]
因杨氏模量和剪切模量均为正值,可推导泊松比最小值为-1。
负泊松比是指受拉伸时, 材料在弹性范围内横向发生膨胀;而受压缩时, 材料的横向反而发生收缩。
这种现象在热力学上是可能的,但通常材料中并没有普遍观察到负泊松比效应的存在。
近年来发现的一些特殊结构的材料具有负泊松比效应,由于其奇特的性能而倍受材料科学家和物理学家们的重视。
负泊松比材料可以分为多孔状负泊松比材料(包括泡沫(Foam)材料和蜂巢状(Honeycomb)结构材料)、负泊松比复合材料及分子负泊松比材料等。
负泊松比材料由于具有不同于普通材料的独特性质,在很多方面具备了其他材料所不能比拟的优势,尤其是材料的物理机械性能有了很大的提高,如提高了材料的剪切模量、材料的抗缺口性能、抗断裂性能以及材料的回弹韧性。
另外,由于材料的泊松比影响到应力波的传输和反射,应力的消除和在裂纹附近的应力分布,所以负泊松比材料适合制造紧固件或安全带,在受外力时材料的横向膨胀可以抵消外力的作用,从而提高这些部件的抗负荷能力。