基于扩展有限元法重力坝动态裂纹开裂模拟.docVIP

基于扩展有限元法重力坝动态裂纹开裂模拟.doc

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
基于扩展有限元法重力坝动态裂纹开裂模拟

基于扩展有限元法重力坝动态裂纹开裂模拟   摘要:基于ABAQUS平台的XFEM功能模拟了Koyna坝在地震荷载下的裂纹开裂扩展。ABAQUS平台的XFEM功能可以模拟动力作用下裂纹的开裂扩展,具有广阔的工程应用前景,能够正确捕捉裂纹开裂的位置以及扩展的路径,但其网格依赖性较强,也不能模拟裂纹的分叉,这方面的研究仍然需要进一步完善。   关键词:XFEM;重力坝;裂纹      引言    扩展有限元法(eXtended Finite Element Method,XFEM)是1999年由美国西北大学Belytschko教授[1, 2]首先提出, 适用于求解不连续问题的一种新型数值方法,它克服了常规有限元方法在求解不连续问题时需要个按照不连续界面划分网格的缺陷,仅需按照材料的几何外形划分网格,无需考虑不连续面的存在。经过近十几年的发展,XFEM已广泛渗透到各个领域[3-5],其理论也不断地得到完善。然而,运用XFEM求解动力作用下的裂纹扩展问题甚少,在国内仅有方修君 [7] 等基于ABAQUS平台用户子程序UEL开发了XFEM程序模拟Koyna坝的地震开裂过程。目前,XFEM已被嵌入进商业有限元软件ABAQUS中,可以模拟静、动作用下的裂纹扩展问题。本文主要基于ABAQUS平台的扩展有限元功能模拟Koyna重力坝在地震作用下的裂纹扩展过程,展示了XFEM法的工程应用前景。   动力XFEM的支配方程    对动力作用下含裂纹体进行受力分析,有    , (1)   其中,表示材料的密度;为柯西应力张量;为体积力向量;表示某点的位置向量;为求解域。    在边界域上,需满足以下条件:    ,,    , (2)   其中,为力边界;为位移边界;为裂纹边界(不连续面力边界);为上指定的面力,为上指定的位移;为的外法向;为不连续边界上的面力。    不考虑裂尖场的奇异性,XFEM的位移模式可表示为    (3)   式中,为求解域中所有结点的集合;为常规有限元的结点处的形函数;为常规有限元部分在结点处的未知量;为Heaviside函数;为Heaviside改进的结点集;为与Heaviside改进相关的结点未知量。    由于和只是位置的函数,因此速度和加速度可以离散为   (4)   (5)    由式(3)~(5),基于虚功原理原理可得原问题的XFEM形式的支配方程为   (6)   其中,                  重力坝的地震响应分析    Koyna混凝土重力坝位于印度的Koyna河上, 坝高103 m , 坝底宽约70 m,坝顶宽度14. 8 m其几何外形见图1。1967年, 该坝坝址区域遭受一次6. 5级强烈地震作用, 地震期间水位约为91.75 m,地震后坝体发生开裂。作为少数几个在强震中破坏且有比较完整记录的重力坝之一,Koyna 坝一直是大坝抗震分析中的经典研究对象。    本文基于ABAQUS平台的XFEM功能,模拟了该重力坝在遭受地震时裂缝的开裂扩展过程,实际记录的Koyna水平向地震加速度为0.49 g,竖向地震加速度为0.34 g,见图2。施加的荷载包括坝体自重、静水压力、地震荷载,未计及动水压力的影响。计算时采用了两种计算网格,见图3。      图1 Koyna坝几何尺寸          (a) 水平向分量(b) 竖向分量   图2 Koyna地震时程         (a) 网格1 (b) 网格2   图3 计算网格    图4给出了两种计算网格下的裂纹形态,可以看出ABAQUS平台的XFEM计算结果的网格依赖性较强。图4(c)同时给出了其他研究者的基于有限元法的计算结果,虽然基于ABAQUS平台的XFEM的网格依赖性较强,但可以捕捉裂纹的位置,它不能模拟裂纹的分叉。图5给出了两种计算网格下坝顶顺河向位移时程曲线,两种计算网格下坝顶位移时程曲线的变换规律一致,网格2时,由于地震未结束时裂缝已贯穿,导致计算终止,位移时程曲线没有遍历整个地震过程(20 s)。      (a) 网格1(b) 网格2         图4不同计算网格下裂纹形态及与文献结果比较      图5 坝顶顺河向位移时程曲线   结论    本文基于ABAQUS平台的XFEM功能模拟了Koyna坝在地震荷载下的裂纹开裂扩展,ABAQUS平台的XFEM功能可以模拟动力作用下裂纹的开裂扩展,具有广阔的工程应用前景,能够正确捕捉裂纹开裂的位置以及扩展的路径,但其网格依赖性较强,也不能模拟裂纹的分叉,这方面的研究仍然需要进一步完善。   参考文献    [1] Belytschko T, Black T. Elastic c

文档评论(0)

fangsheke66 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档