首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 93 毫秒
1.
改进型扩展比例边界有限元法   总被引:3,自引:3,他引:0  
江守燕  李云  杜成斌 《力学学报》2019,51(1):278-288
结合了扩展有限元法(extended finite elementmethods,XFEM)和比例边界有限元法(scaled boundary finite elementmethods,SBFEM)的主要优点,提出了一种改进型扩展比例边界有限元法(improvedextended scaled boundary finite elementmethods,$i$XSBFEM),为断裂问题模拟提供了一条新的途径.类似XFEM,采用两个正交的水平集函数表征材料内部裂纹面,并基于水平集函数判断单元切割类型;将被裂纹切割的单元作为SBFE的子域处理,采用SBFEM求解单元刚度矩阵,从而避免了XFEM中求解不连续单元刚度矩阵需要进一步进行单元子划分的缺陷;同时,借助XFEM的主要思想,将裂纹与单元边界交点的真实位移作为单元结点的附加自由度考虑,赋予了单元结点附加自由度明确的物理意义,可以直接根据位移求解结果得出裂纹与单元边界交点的位移;对于含有裂尖的单元,选取围绕裂尖单元一圈的若干层单元作为超级单元,并将此超级单元作为SBFE的一个子域求解刚度矩阵,超级单元内部的结点位移可通过SBFE的位移模式求解得到,应力强度因子可基于裂尖处的奇异位移(应力)直接获得,无需借助其他的数值方法.最后,通过若干数值算例验证了建议的$i$XSBFEM的有效性,相比于常规XFEM,$i$XSBFEM的基于位移范数的相对误差收敛性较好;采用$i$XSBFEM通过应力法和位移法直接计算得到的裂尖应力强度因子均与解析解吻合\较好.   相似文献   

2.
在比例边界等几何分析的建模剖分过程中会出现交界面网格非匹配现象,给计算分析带来困难。为了能够处理此类问题,提出了基于非重叠Mortar方法的比例边界等几何分析。该方法能在将全域分解为若干子域时,针对每个子域分别建模和剖分网格,交界面网格无需逐点匹配;采用交叉点修正的非均匀有理B-样条基函数构造Lagrange乘子空间,子域交界面连续条件可通过Mortar条件满足;并根据交界面连续条件进行自由度凝聚,得到对称正定的系数矩阵,可直接求解。分片试验、U形结构和半无限空间上的柔性基础等数值算例验证了本文方法的有效性,计算精度满足要求。  相似文献   

3.
广义扩展有限元法及其在裂纹扩展分析中的应用   总被引:1,自引:0,他引:1  
结合广义有限元法(GFEM)和扩展有限元法(XFEM)的特点,提出了一种新的数值方法——广义扩展有限元法(GXFEM)。阐述了广义扩展有限元法的基本原理,对相关公式进行推导,探讨数值实施中需注意的重要问题,给出利用广义扩展有限元法进行断裂分析时应力强度因子的计算方法,编写了广义扩展有限元法程序。通过算例进行了应力强度因子的计算,模拟了结构裂纹的扩展过程。算例结果表明,利用广义扩展有限元法计算裂纹扩展问题,不需要进行过密的网格划分,且网格在裂纹扩展后无需重新剖分,具有相当高的计算精度。  相似文献   

4.
基于有限断裂法和比例边界有限元法提出了一种裂缝开裂过程模拟的数值模型。采用基于有限断裂法的混合断裂准则作为起裂及扩展的判断标准,当最大环向应力和能量释放率同时达到其临界值时,裂缝扩展。结合多边形比例边界有限元法,可以半解析地求解裂尖区域附近的应力场和位移场,在裂尖附近无需富集即可获得高精度的解。计算能量释放率时,只需将裂尖多边形内的裂尖位置局部调整,无需改变整体网格的分布,网格重剖分的工作量降至最少。裂缝扩展步长通过混合断裂准则确定,避免了人为假设的随意性,并可以实现裂缝变步长扩展的模拟,更符合实际情况。通过对四点剪切梁的复合型裂缝扩展过程的模拟,对本文模型进行了验证,并应用于重力坝模型的裂缝扩展模拟,计算结果表明,本文提出的模型简单易行且精度较高。  相似文献   

5.
混凝土断裂的连续-非连续方法   总被引:1,自引:0,他引:1  
采用有限元形函数作为单位分解函数,位移间断用富集节点的附加自由度表示,建立了允许在单元内部位移非连续的局部富集公式以表征混凝土的开裂区域.富集基函数由节点形函数和节点形函数与间断函数的乘积的并集构成.非连续位移的扩展路径完全与网格结构无关.不同于以非协调应变为基础的嵌入非连续模型,对单元的类型没有限制而且间断位移可以贯穿单元边界.局部富集思想与扩展有限元类似,但富集点自由度保持节点位移的物理意义不变,使相邻单元无需进行富集运算.在变分公式中引入混凝土粘结本构定律,推导了考虑断裂过程区非线性影响的基本方程.对混凝土粘结裂纹扩展的数值模拟说明了该计算方法的有效性.  相似文献   

6.
扩展有限元法 (extended finite element method, XFEM) 因具有裂纹几何独立于模拟网格、裂纹扩展时无需网格重分重映、计算精度高等优点,成为裂纹分析的主流数值方法之一. 但该方法在工程实践中存在单元被裂纹分割的几何困难 —— 现有精确几何分割方法实现复杂、计算量大、鲁棒性差. 为克服这一困难, 本文提出一种基于单元水平集的模板分割方法, 用于非连续单元子剖分和数值积分. 首先, 遍历单元水平集值所有形态并建立标准单元分割模板库; 然后, 根据单元水平集值, 对非标准单元进行形态查询和模板插值; 最后, 套用标准单元分割模板实现单元高效分割和子剖分. 将该方法与常规XFEM、改进型XFEM进行结合,从而应用于孔洞、夹杂、裂纹等非连续问题分析中. 算例分析表明, 本文提出的模板分割方法具有较高计算精度. 由于不引入复杂几何操作, 该模板分割方法同时具有较高计算效率和鲁棒性, 故可为XFEM类方法在实际工程应用中提供有效支撑.   相似文献   

7.
一致性高阶无网格法能高效精确地求解连续体问题,尤其是能得到高精度的应力场。本文将该方法拓展到应力解析精度至关重要的裂纹问题(即非连续体问题)的数值分析。采用背景积分网格描述裂纹几何,基于无需增加节点额外自由度的虚拟节点法描述裂纹处位移场的间断,提出了虚拟节点的引入算法和断裂单元的数值积分方法。为进一步模拟裂纹扩展,采用相互作用积分方法计算应力强度因子,裂纹的扩展方向由最大周向应力准则确定。数值结果表明,本文发展方法能够精确地通过间断分片试验;相较于标准的高阶无网格法和低阶一致性无网格法,本文的一致性高阶无网格法显著改善了应力强度因子的计算精度,能够准确预测裂纹扩展路径。  相似文献   

8.
裂纹问题的一致性高阶无网格法   总被引:2,自引:0,他引:2  
一致性高阶无网格法能高效精确地求解连续体问题,尤其是能得到高精度的应力场。本文将该方法拓展到应力解析精度至关重要的裂纹问题(即非连续体问题)的数值分析。采用背景积分网格描述裂纹几何,基于无需增加节点额外自由度的虚拟节点法描述裂纹处位移场的间断,提出了虚拟节点的引入算法和断裂单元的数值积分方法。为进一步模拟裂纹扩展,采用相互作用积分方法计算应力强度因子,裂纹的扩展方向由最大周向应力准则确定。数值结果表明,本文发展方法能够精确地通过间断分片试验;相较于标准的高阶无网格法和低阶一致性无网格法,本文的一致性高阶无网格法显著改善了应力强度因子的计算精度,能够准确预测裂纹扩展路径。  相似文献   

9.
杜成斌  黄文仓  江守燕 《力学学报》2022,54(4):1026-1039
混凝土是一种被广泛应用于土木和水利工程中的准脆性材料, 在各种内外部因素的作用下, 开裂是混凝土结构最为普遍的破坏形式, 准确模拟结构的开裂过程, 对于结构的安全评估至关重要. 将比例边界有限元与非局部宏微观损伤模型相结合提出一种准脆性材料开裂模拟新方法. 以比例边界有限元子域的比例中心作为物质点, 通过两比例中心(物质点对)之间的物质键的正伸长率来定义微细观损伤, 将某点影响域内物质键的微细观损伤加权平均得到该点的宏观拓扑损伤. 再引入能量退化函数, 将宏观拓扑损伤嵌入到比例边界有限元的基本框架中. 充分利用比例边界有限元网格允许存在悬挂节点的优势, 采用四叉树网格离散技术进行快速、高质量的多级网格划分与过渡. 通过一个I型开裂与一个混合型开裂的两个典型算例, 验证了该方法可捕获结构裂纹扩展路径与荷载变形曲线. 与现有的方法相比, 本文的损伤模型可得到更准确的局部开裂损伤带, 结果更为合理, 且具有更高的计算精度和计算效率. 当损伤过程区网格尺寸小于影响域半径的1/5时, 计算结果不存在网格敏感性问题.   相似文献   

10.
基于扩展有限元法的混凝土细观断裂破坏过程模拟   总被引:1,自引:0,他引:1  
扩展有限元法(XFEM)是分析不连续力学问题(特别是断裂问题)的一种有效的数值方法。在常规的有限元位移模式中,基于单位分解的思想加入一个跳跃函数和渐进缝尖位移场来对不连续体附近的节点自由度进行局部加强,从而反映了位移的不连续性。介绍了扩展有限元的基本原理,给出了扩展有限元进行混凝土开裂及裂纹扩展的分析方法,最后采用扩展有限元法模拟了湿筛混凝土单轴拉伸作用下及WinklerL-型混凝土板的细观断裂破坏过程。分析了混凝土裂纹萌生、扩展的过程及破坏形态,数值结果与实验结果吻合良好。研究表明:扩展有限元法通过特定的位移模式,使裂纹两侧不连续位移场的表达独立于网格划分,能有效地模拟混凝土材料细观断裂破坏过程。  相似文献   

11.
The scaled boundary finite element method (SBFEM) is a recently developed numerical method combining advantages of both finite element methods (FEM) and boundary element methods (BEM) and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors (SIFs) directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.The project supported by the National Natural Science Foundation of China (50579081) and the Australian Research Council (DP0452681)The English text was polished by Keren Wang.  相似文献   

12.
殷德胜  尹栓  周宜红 《计算力学学报》2014,31(6):735-741,748
比例边界有限元法SBFEM(Scaled Boundary Finite Element Method)是一种半解析数值方法,在裂缝分析特别是强度因子计算上具有相当高的精度。本文提出了一种用于裂缝分析的基于虚拟结构面的SBFEM与常规FEM的耦合分析方法。首先选取裂缝周边一定范围的计算域,并将结构分成不含裂缝区域和含裂缝区域两部分。然后,对不含裂缝区域,采用FEM进行网格离散;对含裂缝区域,采用SBFEM进行网格离散;两者相互独立,在这两个域内,分别采用各自相应的位移模式。最后通过在SBFEM网格的外边界设置虚拟耦合结构面的模式,实现有限元网格和比例边界有限元网格的耦合。通过两个经典的含裂缝平板的算例研究,探讨了本文方法在I型开裂和混合型开裂分析中,影响应力强度因子精度的因素。算例表明,SBFEM具有的降维和半解析性质,使本文方法在裂缝分析中的前处理简单易行,且计算结果具有相当高的计算精度。  相似文献   

13.
数值流形方法是一种非常灵活的数值计算方法,连续体的有限单元方法和块体系统的非连续变形分析方法只是这一数值方法的特例.数值流形方法中高阶位移函数的构造可通过提高权函数的阶次来实现,这种方法往往需要沿单元边界配置适当的边内节点,这些结点的出现增加了前处理的复杂性,特别是对于大型复杂的空间问题.另一方面,在数值流形方法中可通过缩小单元尺寸(h加密)来提高求解精度.当模拟裂纹扩展时,这种细化策略可用来克服裂纹尖端的奇异性.一个传统的解决方案是细化整个网格,但这会导致计算效率的显著降低.将适合分析的T样条(analysis-suitable T-spline,AST)引入数值流形方法中来建立高阶数值流形方法的分析格式,有效的避免了该问题的出现.AST样条基函数具有线性无关,单位分解,局部加密等许多重要性质,使得其非常适合用于工程设计及分析.在引入AST样条后,可通过改变数学覆盖的构造形式建立不同阶次的数值流形方法分析格式;AST样条自身的局部加密性质也使得数值流形方法中的数学网格局部加密更容易实现.算例结果表明:随着AST样条基函数阶次的提高,数值流形方法的计算结果有了明显的改善;基于AST样条基函数的数值流形方法在保持计算精度的前提下降低了自由度的数量.  相似文献   

14.
从虚功方程出发,结合扩展有限元离散技术与接触条件的非线性互补表述,建立了摩擦接触裂纹问题的扩展有限元非线性互补模型,将不等式接触条件转化为非线性互补类的非光滑方程组,并采用基于广义导数的非光滑阻尼牛顿法求解方程组,无需引入任何额外人工变量以及迭代求解。以含中心倾斜裂纹平板和边裂纹平板为例,运用相互作用积分法计算摩擦接触裂纹的应力强度因子,将其结果与理论解进行对比分析,该方法都能给出精确的计算结果;基于扩展有限元方法对单轴压缩作用下倾斜裂纹扩展过程进行了数值模拟,计算结果表明,受压裂纹数值结果与实验结果比较吻合,从而验证了本文方法的有效性与正确性。  相似文献   

15.
The extended finite element method (X-FEM) is a numerical method for modeling strong (displacement) as well as weak (strain) discontinuities within a standard finite element framework. In the X-FEM, special functions are added to the finite element approximation using the framework of partition of unity. For crack modeling in isotropic linear elasticity, a discontinuous function and the two-dimensional asymptotic crack-tip displacement fields are used to account for the crack. This enables the domain to be modeled by finite elements without explicitly meshing the crack surfaces, and hence quasi-static crack propagation simulations can be carried out without remeshing. In this paper, we discuss some of the key issues in the X-FEM and describe its implementation within a general-purpose finite element code. The finite element program Dynaflow™ is considered in this study and the implementation for modeling 2-d cracks in isotropic and bimaterial media is described. In particular, the array-allocation for enriched degrees of freedom, use of geometric-based queries for carrying out nodal enrichment and mesh partitioning, and the assembly procedure for the discrete equations are presented. We place particular emphasis on the design of a computer code to enable the modeling of discontinuous phenomena within a finite element framework.  相似文献   

16.
基于单元破裂的岩石裂纹扩展模拟方法   总被引:3,自引:0,他引:3  
传统离散元方法在处理破裂问题时, 采用界面上的准则进行判断, 裂纹只能沿着单元边界扩展. 当物理问题存在宏观或微观裂隙时, 在界面上应用准则具有其合理性; 而裂纹沿着单元边界扩展, 使得裂纹路径受网格影响较大, 扩展方向受到限制. 针对上述情况, 可以基于单元破裂的方式, 构建连续- 非连续单元法, 并应用于岩石裂纹扩展问题的模拟. 该方法在连续计算时, 将单元离散为具有物理意义的弹簧系统, 在局部坐标系下由弹簧特征长度、面积求解单元变形和应力, 通过更新局部坐标系和弹簧特征量, 可进一步计算块体大位移、大转动, 连续问题计算结果与有限元一致, 同时提高了计算效率. 在此基础上, 引入最大拉应力与莫尔-库伦的复合准则, 判断单元破裂状态和破裂方向, 并采用局部块体切割的方式, 在单元内形成初始裂纹. 裂纹两侧相应增加新的计算节点, 同时引入内聚力模型描述裂纹两侧的法向、切向作用与张开度及滑移变形之间的关系. 按此方式, 裂纹尖端处的扩展路径可穿过单元内部和单元边界, 在扩展方向的选取上更为准确. 最后, 通过三点弯曲梁、单切口平板拉伸、双切口试样等典型数值试验, 模拟裂纹在拉伸、压剪等各种应力状态下的扩展问题, 并对岩石单轴压缩试验的破坏过程进行模拟, 分析裂纹形成与应力-应变曲线各阶段之间的对应关系. 结果表明: 连续-非连续单元法通过单元内部破裂的方式, 可以显示模拟裂纹萌生、扩展、贯通直至形成宏观裂缝的过程.  相似文献   

17.
In this paper, the extended finite element method (X-FEM) is considered for the analysis of fretting fatigue problems. A two-dimensional implementation of the X-FEM is carried out within the finite element software ABAQUS? by means of user subroutines, and crack propagation in fretting fatigue problems is investigated. On utilizing the non-linear contact capabilities of this code, the numerical technique is applied to a specimen-indenter model. The use of the X-FEM facilitates very accurate stress intensity factor computations on relatively coarse meshes, and furthermore, no remeshing is required for crack growth simulations. The implementation is applied to complete and incomplete contact fretting problems. A study of crack growth is conducted for several bulk loads applied to the specimen, and the influence of the initial crack angle is ascertained. The numerical simulations reveal the merits of applying the X-FEM to fretting fatigue problems.  相似文献   

18.
数值流形方法及其在岩石力学中的应用   总被引:9,自引:0,他引:9  
李树忱  程玉民 《力学进展》2004,34(4):446-454
数值流形方法是目前岩石力学分析的主要方法之一.该方法起源于不连续变形分析,主要用于统一求解连续和非连续问题,其核心技术是在分析时采用了双重网格:数学网格提供的节点形成求解域的有限覆盖和权函数;而物理网格为求解的积分域.数学网格被用来建立数学覆盖,数学覆盖与物理网格的交集定义为物理覆盖,由物理覆盖的交集形成流形单元.流形方法的优点在于它使用了独立的数学和物理网格,具有和有限元明显不同的定义形式,且数学网格对于同一问题不同的求解精度的需求可以很方便地细化.由于该方法考虑了块体运动学,可以模拟节理岩体裂隙的开裂和闭合过程,因而在岩石力学中得到了广泛应用,近年来许多学者对该方法进行了研究.本文简要叙述了节理岩体的数值方法从连续到非连续的发展过程,详细地介绍了数值流形方法的组成和数值流形方法在岩石力学及其相关领域的研究和发展概况,最后就作者所关心的一些问题,如三维问题的数值流形方法、数值流形方法在物理非线性问题和裂纹扩展问题中的应用、相关的耦合方法等进行了探讨.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号