![点击分享此内容可以赚币 分享](/master/images/share_but.png)
软黏土不排水循环应力应变关系的数值模拟.pdf
《软黏土不排水循环应力应变关系的数值模拟.pdf》由会员分享,可在线阅读,更多相关《软黏土不排水循环应力应变关系的数值模拟.pdf(7页珍藏版)》请在咨信网上搜索。
1、第 3期 2 0 1 5年 6月 水利水运 工程学报 AND ENGD旺 RI NG No 3 J u n 2 0 1 5 D O I :1 0 1 6 1 9 8 j c n k i 1 0 0 9 6 4 0 X 2 0 1 5 0 3 0 1 2 王哲学 , 王建化,程星磊 软黏土不排水循环应力应变关系的数值模拟 J 水利水运工程学报, 2 0 1 5 ( 3 ) :8 1 8 7 ( WA N G Z h e x u e , WA N G J i a n h u a ,C H E N G X i n g l e i N u me r i c a l s i mu l a t i o n
2、 o f u n d r a i n e d c y c l i c s t r e s s s t r a i n r e s p o n s e o f s o f t c l a y J H y d r o S c i e n c e a n d E n g i n e e r i n g ,2 0 1 5 ( 3 ) : 8 1 8 7 ) 软黏 土不排水循环应 力应 变关 系 的数值模拟 王哲学 一,王建华 一,程星磊 , ( 1 天津大学 水利工程仿真与安全国家重点实验室, 天津3 0 0 0 7 2 ; 2 天津大学 岩土工程研究所,天津 3 0 0 0 7 2 ) 摘 要 :对
3、程星磊等提出的总应力形式增量弹塑性本构模型进行二次开发, 以模拟复杂应力状态下软黏土的响 应。通过 N e w t o n R a p h s o n算法, 对材料非线性问题进行迭代求解; 针对本构模型中应力反向等一系列关键性问 题, 应用欧拉切线算法编写了有限元程序, 并结合子增量方法提高了计算精度。预测了软黏土在轴对称应力状 态下 的响应 , 得到了应力应变关系 曲线 , 将其与单元预测结果进 行 比较 , 二者 趋于一致 , 从而 验证 了有 限元程序 编写的合理性。利用该程序模拟三轴不固结不排水试验, 模拟结果与试验吻合良好, 表明该本构模型有限元程 序可以反映轴对称应力状态下软黏土的
4、不排水应力应变特性 , 可应用于更加复杂边值问题的模拟计算。 关 键 词 : 循环荷载; 饱和软黏土; 应力应变; 本构模型; 欧拉切线算法; 三轴试验 中图分类号: T U 4 3 文献标志码 : A 文章编号 : 1 0 0 9 6 4 0 X( 2 0 1 5 ) 0 3 - 0 0 8 1 - 0 7 嵌固在海洋软黏土中的锚 固基础在风 、 浪等荷载作用下 , 为维持上部结构工作性能, 必须保持 自身稳定 , 因此需要合理分析软黏土中锚 固基础在循环荷载作用下的变形失稳过程。锚固基础的变形不仅与基础本身 有关 , 海洋中软土地基对其影响更为重要。因此研究循环荷载下软黏土的不 固结不排水
5、应力应变关系具有 重要意义。J H P r 6 v o s t 1 - 2 J 提 出了基于运动硬化塑性理论 的多面模型 , 但模型的变形分析很复杂 , 需要记忆 所有假设的多个次加载面形状 , 并分析它们的移动 , 导致在数值应用 中计算过程复杂 。Y F D a f a l i a s 等 _ 3 j 建立了包括 1 个正常屈服面和 1 个次加载面的两面模型 , 大大简化 了多面模型 ; 后来他们又将两面模型中的 次加载面简化为一点 , 提出了单面模型 。两面和单面模型都 叫边界面模型_ 6 J 。由于改进初期 的边界面 模型映射中心 固定在原点 , 无法合理描述土单元 的应力应变滞 回曲
6、线 J , 且硬化模 量场 的演化过程表述较 为复杂, 不便 于应用到数值分析中。程星磊等 。 。 采用总应力形式在偏应力空问中构建硬化模量场 , 硬化模 量插值函数数学表达简洁 , 引入模型参数较少 , 硬化模量场 的演化规律表述简单 , 跟踪循环应力路径所需记 忆 的参数较少 , 便于计算应用 。通过上述分析 , 针对程星磊等建立 的软黏土循环 弹塑性本构模型进行二次开 发 , 并对循环三轴试验进行数值模拟 , 通过计算结果与试验结果的比较 , 验证该本构模型数值分 析方法的可 行 性 1 本构模型简述 采用 M i s e s 屈服准则构建边界面约束方程 : F = 1 5 i ; 一A
7、 =0 收稿 日期 : 2 0 1 4 - 0 9 3 0 基金项 目:国家 自然科学基金资助项 目( 5 1 1 7 9 1 2 0 ) ; 博 士点基金项 目( 2 0 1 3 0 0 3 2 1 1 o o 4 5 ) 作者简介:王哲学( 1 9 8 9 一) ,男, 河北承德人 , 硕士研究生 , 主要从事土动力学与海洋岩土工程研究。 E - m a i l : w a n g z h e x u e 0 1 1 7 1 2 6 e o m通信作者 : 王建华 ( t d w j h t j u e d u c n ) 学兔兔 w w w .x u e t u t u .c o m 8
8、 2 水 利 水 运 工 程 学 报 2 0 1 5年 6月 式 中: ; 为边界面上偏应力张量 ; A 。为边界面半径。 依据经典弹塑性理论 , 有总偏应变增量计算式 : d d s + 3 S ij- z ( 2 ) 式中: d e 为总偏应变增量 ; G为弹性剪切模量 ; H 为塑性硬化模量 ; d s 为当前应力点的偏应力增量 。 引入弹塑性硬化模量 , 令 : 1 H =1 ( 2 G)+1 H ( 3 ) 在初始加载或应力反 向时, 塑性变形模量为无穷大 : H =H =2 G ( 4 ) 图 1给出了应力空 间中加载和卸载时确定弹塑性 模量的映射规则 图解 。这里边界面针对九维应
9、力空 间 中的超 曲面。图中 s 为映射 中心点 , s 为 当前应力点 , S 为通过映射 中心点和当前应力点 的射线 与边界 面的 交点 , 此点称为像应力点 , s 为应力反向点。当发生应 力反向时 , 如图 1 ( b ) , s 则变为映射中心 , 其他各点定 义不变 , 此时会出现新的应力反 向点 s 。定义 和 分别为当前应 力点 和初始加 载点 到像 应力点 的距离 。 弹塑性模量插值规则为 : ( a ) 加载时 ( h ) 卸载时 图 l 弹塑性模量 的映射规则 图解 F i g 1 El a s t o - p l a s t i c ma pp i n g r u l
10、e f o r l o a d i ng a n d u nl o a di n g H=( t3 6 o ) H ( 5 ) 式中: 日 为弹塑性模量最大值 , 取 2 G; 为反映剪应变累积速率和大小的模型参数 。 在不固结不排水条件下 , 边界面半径 A 。 和弹性模量 G可通过不 固结不排水循环三轴试验得到。在三轴 应力状态下 , 式( 1 ) 简化为下式 : F=( : 一 ) 一A =0 ( 6 ) 式中 : 为试样轴向应力 ; 为试样径向应力。 由于初始应力应变响应均遵循同一初始加荷 曲线 , 且可用式( 7 ) 拟合 , 则拟合曲线 的渐近线所对应 的强 度就是边界面半径 A
11、。 , 其值等于 1 b; 曲线 的初始切线模量记为 E, 其值等于 1 a。另外 , 弹性 剪切模量 G 可 由式 ( 8 ) 求得 : = ( 口+6 ) ( 7 ) G=E 2 ( 1+ ) ( 8 ) 式 中: 口和 b 为拟合 函数 中的待定系数 ; 为泊松比, 不排水时取 0 5 。 参数 可以用下式表示 : = x 0 d y 8 。 ( 9 ) 式 中: 为反映累积偏应变大小的待定系数 ; d y 。 为八面体累积剪应变增量 , 以此描述偏应变累积变化 的速 率和趋势。对试验数据采用式( 1 0 ) 拟合 : 8 。=c ( 1 0 ) 式中: c , d为待定系数 ; 为循环
12、次数 。则有 d 7 。=c d N d ( 1 1 ) 下面结合某一具体三轴试验 , 给出待定系数 c, d, I x 的确定方法。图 2给出了该三轴试验得到的静应 力 比1 8 = 和 5时八面体累积剪应变随振次的变化关系。用式 ( 1 0 ) 拟合试验数据 , 如 图中虚线所 ,aT 0 3 0 示 , 得到不同静应力 比T S ,a T 不同循环应力 比 8 c 下 的待定系数 c 和 d值 。通过式 ( 1 1 ) 可得应力循 环一周八面体累积剪应变增量 。 。 学兔兔 w w w .x u e t u t u .c o m 第 3 期 王哲学 , 等 : 软黏土不排水循环应力应变关
13、系 的数值模拟 8 3 1 00 2 0 0 3 0 0 4 O 0 5 0 0 6 0 0 7 0 0 8 0 0 9 00l 0 0 0 _7v ( a ) 5=0 3 1 0 0 2 0 0 3 0 0 4 0 0 5 0 0 6 0 0 7 0 0 8 0 0 9 0 0 1 0 0 0 _7v ( h ) 5 = 0 5 图 2 三轴压缩时八面体累积剪应变随循环次数的变化关系曲线 F i g 2 R e l a t i o n c u r v e s o f o c t a h e d r a l c u mu l a t i v e s h e a fi n g s t r a i
14、 n a n d c y c l e t i me s i n t r i a x i a l c o mp r e s s i o n 进一步 , 可通过试算方法确定 。 , 即利用式 ( 9 ) 将 8 D 代人式 ( 5 ) , 给定一个 。 , 通过式( 2 ) 计算得到破 坏振次下的八面体累积剪应变。这里规定按照静偏应力作用下试样的八面体静剪应变与静偏应力与循环应 力共 同作用下八面体循环累积剪应变之和达到 1 4 的标准确定循环破坏次数 。若模型计算所得破坏振 次 下的八面体累积剪应变与试验结果一致 , 那么此时的 值 即为该试验应力条件下待定系数 的取值。 这样 , 不 同静应力
15、 比 T S 8 , , 不同循环应力比 , 下 的参数 c , 。 都 已给出( 见表) 。则对任意 循环应力 比下的参数值 , 可通, a 过 T 线性插值方法确定 。 8 ,cyT d 1 表 1 不同静应力比、 循环应力比下拟合函数中的待定系数 C , d及p。 取值 T a b 1 Un d e t e r mi n e d c o e f f i c i e n t s C a n d d o f fi t t i n g f u n c t i o n a n d 0 w i t h d i f f e r e n t s t a t i c s t r e s s r a t
16、i o s a n d c y c l i c s t r e s s r a t i o s f8 ,d T s J T 8 ,c y s d C d I z o T 8 r s f T 8 8 j C d 0 0 5 3 0 0 2 4 7 4 3 1 2 O1 5 3 2 7 2 0 0 3 9 2 0 5 9 4 9 4 0 7 02 0 3 4 8 2 0 0 5 2 3 0 1 5 8 8 4 1 1 5 2 O 6 3 7 4 0 0 3 6 3 1 0 2 3 3 6 0 48 7 3 3 6 8 3 0 0 3 0 5 0 4 2 3 0 O 0 4 0 4 1 5 9 4
17、0 2 7 8 6 0 0 3 3 4 0 3 5 9 1 0 0 6 03 3 4 9 8 4 0 0 3 0 0 0 7 03 0 1 0 3 61 3 8 2 5 8 4 0 0 2 4 2 1 1 8 8 9 5 0 3 3 3 21 l 9 0 4 0 至此 , 软黏土不 固结不排水条件下循环本构关系及模型参数的确定方法 已完全给出。 2 本构模型二次开发 2 1 平衡迭代过程 增量形式 的有 限元格式 中需要根据本步或迭 代计 算得到的位移增量 , 决定 s, A o - , 进而确定本步或迭代 结束时的弹塑性状态。 非线性平衡 方程采用 N e w t o n R a p h s
18、 o n迭代算 法 求解 。在迭代过程中, 把整个加载过程分成若干分段 , 荷 载以增量形式施加在结构上, 在增量 内部进行若干迭代 , 直至达到所需精度为止。整个流程融合 了增量与迭代 的 过程 。在增量内的迭代过程( 见图 3 ) 作如下说明 : 根据当前的应力状态 , 确定弹塑性矩阵D , 进而确 定刚度矩阵 。 由 A u=K A R , 得 到位移增量 。其 R 司 l R 图 3 增量 内迭代过程 Fi g 3 I t e r a t i v e p r o c e s s i n i nc r eme n t 中, A R为外加荷载增量 , Au为位移增量。由 A e=B Au,
19、 得到应变增量。采用应力更新算法 , 通过当前 2 1 O 9 8 7 6 5 4 3 2 l O 2 1 O 9 8 7 6 5 4 3 2 l 0 学兔兔 w w w .x u e t u t u .c o m 8 4 水 利 水 运 工 程 学 报 2 0 1 5 年 6月 应力状态 以及上一步得到的应变增量计算应力增量 A o, 进而得到新的应力状态 ” 。判断 内力 、 外力是 否平衡 , 位移是否满足收敛标准 : 若不满足, 返 回 , 继续进行迭代 ; 若满足, 进入下一增量计算。 2 2 应力更新算法 对第步的应力更新算法做出详细阐述, 由应力应变关系式( 2 ) , 可得 :
20、 d s =2 G d e 一1 5 ( 2 GH)s S k l d e 由上述关 系式可得到偏应力增量。静水压力增量可由弹性虎克定律得到 : d o - =3 K d 6 由此得到应力增量大小 : d o - =d o- +d 将式 ( 8 ) 和( 9 ) 代入式( 1 0 ) , 有 : ( 1 2 ) ( 1 3 ) ( 1 4 ) d o - =3 K 8 6 +2 G d e 一1 5 ( 2 G一日) s A s d e w ( 1 5 ) 得到增量形式的弹塑性应力应变关系 , 整理成矩阵形式 : d o =D印 d e ( 1 6 ) 式 中: d 为应力张量 ; d e为应
21、变张量 ; D 为弹塑性矩 阵, 由于在应力应变关系 曲线上 D叩 为切线方 向, 因此 也 叫切线矩阵。 因此, 本构关系的积分式为: A o=I D 印 d e ( 1 7 ) J 在计算过程中, 采用式( 1 6 ) 将式 ( 1 7 ) 线性化 。但是这种计算方法必然带来偏差 。可采取减小增量大小 的方法减小这种偏差。 王勖成等 提出子增量方法将总的应变增量分成若干个子增量。每一个子增量结束时的弹塑性状态 作为下一个子增量的初始状态 , 所有子增量计算结束 时的弹塑性状态作 为本步增量 的最终状态 。虽然子增 量法提高了计算精度 , 加快了后继迭代 的收敛速度 , 但是子增量数 目的增
22、加使得计算量也相应增加。因此确 定适当子增量数 目很重要。本程序根据增量步开始 时的应力偏量 s 和增量步结束时的弹性预测应力偏量 s 之间的夹角 确定子增量数 : N=1+ , 其 中: k可由计算精度确定 , 这里 k=0 O 1 。 当应力状态趋于边界面附近时 , 由式 ( 5 ) 得硬化模量趋于零 , 由式 ( 1 7 ) 得偏应力增量趋于零 , 因此应力 状态点将始终保持在边界面以内。但实际计算 中, 由于采用式( 1 6 ) 所示的欧拉切线方法对式 ( 1 7 ) 进行线性 化 , 而欧拉方法只有在偏应变增量趋于无限小的时候才能足够准确 , 因此在边界 面附近时, 应力增量 d o
23、 - 有 可能大于零 , 从而使应力状态超 出边界面。但边界面是破坏面 , 土体 的应力状态不允许超过边界面。这将导 致计算出错。采用子增量方法虽然大大减小了偏应变增量 的大小 , 但是也不足以保证应力状态不超过边 界 面。为此 , 将超出边界面的点沿边界面的半径方向径向拉回至边界面 , 并认为这一点的应力状态即为当前的 应力状态 。具体做法是 , 计算应力空间原点与当前应力点之间距离 , 按照边界面半径与该距离 比例对当前应 力点各个应力分量进行线性缩放 , 则缩放后的应力点必然落在边界面上。进一步利用修正后 的应力状态建 立 与之 相适 应 的弹塑 性矩 阵 D 。这 样 就保 证 了任
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 黏土 排水 循环 应力 应变 关系 数值 模拟
![提示](https://www.zixin.com.cn/images/bang_tan.gif)
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【mom****oy】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【mom****oy】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。