共查询到18条相似文献,搜索用时 109 毫秒
1.
提出了一种有限元模拟裂纹扩展的单元子划分结合子结构的方法。本方法中,裂纹可以进入或穿过一个单元,或沿单元的边界扩展,因此裂纹可以沿任意路径扩展而不受初始网格的限制。对上述几类包含裂纹的单元按照裂纹的路径进行子划分,覆盖一条裂纹的所有子划分单元就组成了一个子结构,子结构规模随裂纹的扩展而增大。子结构中因单元子划分而新增的结点自由度,通过自由度的凝聚用初始网格结点的自由度表示,因此结构整体分析的总自由度不变。以上述方法为基础建立了裂纹萌生和扩展的准则。用本文的方法分析了单(双)材料无限大平面中心(界面)裂纹的裂尖场,验证了本文方法的精度,并模拟了颗粒复合材料中微裂纹在颗粒、基体和界面中逐步扩展的过程,考核了本文方法对复杂裂纹扩展问题模拟的适用性。 相似文献
2.
扩展有限元法 (extended finite element method, XFEM) 因具有裂纹几何独立于模拟网格、裂纹扩展时无需网格重分重映、计算精度高等优点,成为裂纹分析的主流数值方法之一. 但该方法在工程实践中存在单元被裂纹分割的几何困难 —— 现有精确几何分割方法实现复杂、计算量大、鲁棒性差. 为克服这一困难, 本文提出一种基于单元水平集的模板分割方法, 用于非连续单元子剖分和数值积分. 首先, 遍历单元水平集值所有形态并建立标准单元分割模板库; 然后, 根据单元水平集值, 对非标准单元进行形态查询和模板插值; 最后, 套用标准单元分割模板实现单元高效分割和子剖分. 将该方法与常规XFEM、改进型XFEM进行结合,从而应用于孔洞、夹杂、裂纹等非连续问题分析中. 算例分析表明, 本文提出的模板分割方法具有较高计算精度. 由于不引入复杂几何操作, 该模板分割方法同时具有较高计算效率和鲁棒性, 故可为XFEM类方法在实际工程应用中提供有效支撑. 相似文献
3.
表面裂纹疲劳扩展的数值模拟 总被引:2,自引:1,他引:1
建立了一种无形状约束的模拟表面裂纹在线弹性断裂力学条件下疲劳扩展的数值方法,并研究了表面疲劳裂纹形状演化和裂纹尖端应力强度因子(SIF)的分布特征。该方法以三维有限单元技术和Paris疲劳裂纹扩展规律为基础,并在裂纹扩展增量计算中考虑了裂纹闭合影响。本文第一部分主要介绍模拟三维疲劳裂纹扩展的数值方法的理论背景和相关的技术细节。着重分析和讨论基于三维有限单元法计算裂纹SIF所涉及的几个主要问题:裂纹尖端单元网格密度对估算精度的影响;自由表面的影响及其修正方法;裂纹尖端非正交单元网格的影响及修正方法。 相似文献
4.
广义扩展有限元法及其在裂纹扩展分析中的应用 总被引:1,自引:0,他引:1
结合广义有限元法(GFEM)和扩展有限元法(XFEM)的特点,提出了一种新的数值方法——广义扩展有限元法(GXFEM)。阐述了广义扩展有限元法的基本原理,对相关公式进行推导,探讨数值实施中需注意的重要问题,给出利用广义扩展有限元法进行断裂分析时应力强度因子的计算方法,编写了广义扩展有限元法程序。通过算例进行了应力强度因子的计算,模拟了结构裂纹的扩展过程。算例结果表明,利用广义扩展有限元法计算裂纹扩展问题,不需要进行过密的网格划分,且网格在裂纹扩展后无需重新剖分,具有相当高的计算精度。 相似文献
5.
对于已含初始裂纹平台管节点的寿命预测很大程度上依靠应力强度因子的精确值,而复杂载荷条件下的节点应力强度因子的计算尚无参数方程直接确定。本文提出了一种含表面裂纹的K节点的有限元网格产生方法,即把整个K节点划分为几个子区域,每个子区域的网格具有不同类型的单元和不同的密度。这种方法在控制网格密度,尤其是控制沿着裂纹边缘单元的边长比方面有其独特的优越性,当所有子区域的网格自动产生后,容易得到整个结构的有限元模型。同时用J积分和位移外推插值法分别计算了一个K型节点沿着裂纹前缘的应力强度因子值,发现:试验得到的应力强度因子值和提出的模型计算结果非常吻合,证明了所提有限元模型的准确性。 相似文献
6.
对于已含初始裂纹平台管节点的寿命预测很大程度上依靠应力强度因子的精确值,而复杂载荷条件下的节点应力强度因子的计算尚无参数方程直接确定.本文提出了一种含表面裂纹的K节点的有限元网格产生方法,即把整个K节点划分为几个子区域,每个子区域的网格具有不同类型的单元和不同的密度.这种方法在控制网格密度,尤其是控制沿着裂纹边缘单元的边长比方面有其独特的优越性,当所有子区域的网格自动产生后,容易得到整个结构的有限元模型.同时用J积分和位移外推插值法分别计算了一个K型节点沿着裂纹前缘的应力强度因子值,发现:试验得到的应力强度因子值和提出的模型计算结果非常吻合,证明了所提有限元模型的准确性. 相似文献
7.
非网格重剖分模拟宏观裂纹体的扩展有限单元法(Ⅰ:基础理论) 总被引:1,自引:1,他引:0
基于有限元计算网格,扩展有限单元法通过建立特殊的广义节点插值形式来描述含裂缝体的不连续位移场,避免了有限元法模拟裂缝时需要的网格重划分。进而,本文从虚功原理出发,在有限元法框架内完整地推导了能模拟宏观裂纹力学场的扩展有限元法实现公式,在理论上更全面地考虑了内部裂纹面上分布外载荷及缝内粘连材料刚度的影响,并提出了构建统一的扩展有限单元刚度阵形成模式,保证了与传统有限单元方式的协调一致。文中对方法的实现过程也做了详细阐述,给出了通用的计算公式,确保了算法的可行性。 相似文献
8.
9.
10.
数学网格和物理网格分离的有限单元法(II):粘聚裂纹扩展问题中的应用 总被引:3,自引:3,他引:0
强化有限单元法将物理网格与数学网格分离开来,可以方便地描述非连续变形;粘聚区域模型是模拟断裂过程区作用最简单有效的方法,且可以避免裂纹尖端的应力奇异性.本文以平面问题为例,将强化有限单元法与粘聚区域模型相结合,利用富集数学节点描述任意粘聚裂纹扩展过程中的非连续变形问题,提出了裂纹扩展过程中数学节点富集和数学单元定义的方法.本文还导出了与平面4~8节点平面等参单元对应的8~16节点粘聚裂纹单元列武.最后,通过三点弯梁的裂纹扩展过程模拟验证了本文提出的粘聚裂纹扩展模拟方法的有效性. 相似文献
11.
基于单元破裂的岩石裂纹扩展模拟方法 总被引:3,自引:0,他引:3
传统离散元方法在处理破裂问题时, 采用界面上的准则进行判断, 裂纹只能沿着单元边界扩展. 当物理问题存在宏观或微观裂隙时, 在界面上应用准则具有其合理性; 而裂纹沿着单元边界扩展, 使得裂纹路径受网格影响较大, 扩展方向受到限制. 针对上述情况, 可以基于单元破裂的方式, 构建连续- 非连续单元法, 并应用于岩石裂纹扩展问题的模拟. 该方法在连续计算时, 将单元离散为具有物理意义的弹簧系统, 在局部坐标系下由弹簧特征长度、面积求解单元变形和应力, 通过更新局部坐标系和弹簧特征量, 可进一步计算块体大位移、大转动, 连续问题计算结果与有限元一致, 同时提高了计算效率. 在此基础上, 引入最大拉应力与莫尔-库伦的复合准则, 判断单元破裂状态和破裂方向, 并采用局部块体切割的方式, 在单元内形成初始裂纹. 裂纹两侧相应增加新的计算节点, 同时引入内聚力模型描述裂纹两侧的法向、切向作用与张开度及滑移变形之间的关系. 按此方式, 裂纹尖端处的扩展路径可穿过单元内部和单元边界, 在扩展方向的选取上更为准确. 最后, 通过三点弯曲梁、单切口平板拉伸、双切口试样等典型数值试验, 模拟裂纹在拉伸、压剪等各种应力状态下的扩展问题, 并对岩石单轴压缩试验的破坏过程进行模拟, 分析裂纹形成与应力-应变曲线各阶段之间的对应关系. 结果表明: 连续-非连续单元法通过单元内部破裂的方式, 可以显示模拟裂纹萌生、扩展、贯通直至形成宏观裂缝的过程. 相似文献
12.
基于扩展有限单元法的短纤维/橡胶复合材料裂纹扩展分析 总被引:1,自引:0,他引:1
为了研究短纤维/橡胶复合材料裂纹在拉伸载荷作用下的扩展演化规律,应用扩展有限元法对预制裂纹的短纤维/橡胶复合材料的裂纹扩展进行了数值模拟。采用随机顺序吸附算法在ABAQUS软件中生成短纤维分布模型,分析了失效准则参数(最大许用主应力和裂纹表面能)、短纤维细观参数(体积含量、长度和取向角)对裂纹扩展行为的影响规律,探究了短纤维的增韧机理,数值结果与试验结果吻合。研究结果表明:降低基体材料的损伤阈值可以实现裂纹的扩展,而不改变裂纹扩展特征;短纤维增强相的存在会改变裂纹的扩展方向,使裂纹的扩展路径呈现复杂性从而提高材料抗撕裂性能;裂纹扩展路径曲折度随短纤维含量、长度和取向角的增加而增大,试样拉断时所需要消耗的功增大,抗撕裂性提高。 相似文献
13.
Yah Xiangqiao 《Acta Mechanica Solida Sinica》2005,18(3):234-241
This paper presents an extension of a displacement discontinuity method with cracktip elements (a boundary element method) proposed by the author for fatigue crack growth analysis in plane elastic media under mixed-mode conditions. The boundary element method consists of the non-singular displacement discontinuity elements presented by Crouch and Starfield and the crack-tip displacement discontinuity elements due to the author. In the boundary element implementation the left or right crack-tip element is placed locally at the corresponding left or right crack tip on top of the non-singular displacement discontinuity elements that cover the entire crack surface and the other boundaries. Crack growth is simulated with an incremental crack extension analysis based on the maximum circumferential stress criterion. In the numerical simulation, for each increment of crack extension, remeshing of existing boundaries is not required because of an intrinsic feature of the numerical approach. Crack growth is modeled by adding new boundary elements on the incremental crack extension to the previous crack boundaries. At the same time, the element characteristics of some related elements are adjusted according to the manner in which the boundary element method is implemented. As an example, the fatigue growth process of cracks emanating from a circular hole in a plane elastic plate is simulated using the numerical simulation approach. 相似文献
14.
Mechanism of quasi-static crack branching in brittle solids has been analyzed by a modified displacement discontinuity method. It has been assumed that the pre-existing cracks in brittle solids may propagate at the crack tips due to the initiation and propagation of the kink (or wing) cracks. The originated wing cracks will act as new cracks and can be further propagated from their tips according to the linear elastic fracture mechanics (LEFM) theory. The kink displacement discontinuity formulations (considering the linear and quadratic interpolation functions) are specially developed to calculate the displacement discontinuities for the left and right sides of a kink point so that the first and second mode kink stress intensity factors can be estimated. The crack tips are also treated by boundary displacement collocation technique considering the singularity variation of the displacements and stresses near the crack tip. The propagating direction of the secondary cracks can be predicted by using the maximum tangential stress criterion. An iterative algorithm is used to predict the crack propagating path assuming an incremental increase of the crack length in the predicted direction (straight and curved cracks have been treated). The same approach has been used for estimating the crack propagating direction and path of the original and wing cracks considering the special crack tip elements. Some example problems are numerically solved assuming quasi-static conditions. These results are compared with the corresponding experimental and numerical results given in the literature. This comparison validates the accuracy and applicability of the proposed method. 相似文献
15.
This paper presents a numerical approach for modeling the interaction between multiple cracks in a rectangular plate under cyclic loads. It involves the formulation of fatigue growth of multiple crack tips under ruixed-mode loading and an extension of a hybrid displacement discontinuity method (a boundary element method) to fatigue crack growth analyses. Because of an intrinsic feature of the boundary element method, a general growth problem of multiple cracks can be solved in a single-region formulation. In the numerical simulation, remeshing of existing boundaries is not necessary for each increment of crack extension. Crack extension is conveniently modeled by adding new boundary elements on the incremental crack extension to the previous crack boundaries. As an example, the numerical approach is used to analyze the fatigue growth of three parallel cracks in a rectangular plate. The numerical results illustrate the validation of the numerical approach and can reveal the effect of the geometry of the cracked plate on the fatigue growth. 相似文献
16.
This paper presents a numerical approach of fatigue growth analysis of cracks emanating from a hole in infinite elastic plate
subjected to remote loads. It involves a generation of Bueckner’s principle and a hybrid displacement discontinuity method
(a boundary element method) proposed recently by the senior author of the paper. Because of an intrinsic feature of the boundary
element method, a general crack growth problem can be solved in a single region formulation. In the numerical simulation,
for each increment of crack extension, remeshing of existing boundaries is not necessary. Crack extension is modeled conveniently
by adding new boundary elements on the incremental crack extension to the previous crack boundaries. As an example, fatigue
growth process of an inclined crack in an infinite plate under uniaxial cycle load is modeled to illustrate the effectiveness
of the numerical approach. In addition, fatigue growth of cracks emanating from a circular hole in infinite elastic plate
subjected to remote loads is investigated by using the numerical approach. Many numerical results are given 相似文献
17.
从传统的混凝土分离裂缝模型和弥散裂缝模型出发,使用不同的模型模拟不同发展阶段的混凝土裂缝,并利用无网格方法可以灵活布置节点和边界的优点,使分离裂缝模型的使用得到很大的简化。通过在宏观裂缝表面布置基于试验结果的混凝土与混凝土界面单元,正确地模拟了裂缝表面抗剪能力的变化。算例表明,使用本文方法可以准确模拟斜拉破坏混凝土梁的破坏过程,且精度和数值稳定性要高于传统的有限元弥散裂缝方法。 相似文献
18.
节理岩体爆生裂纹扩展动态焦散线模型实验研究 总被引:5,自引:0,他引:5
应用动态焦散线测试系统,模拟含节理岩体断裂爆破过程,进行了PMMA模型透射式动态焦散线实验,着重研究了爆炸初始裂纹与节理面夹角不同的情况下,裂纹尖端动态强度因子的变化规律、裂纹穿过节理面的扩展规律、以及炮孔与节理距离不同时裂纹穿过节理扩展的规律。实验结果分析表明,爆生裂纹穿过节理面时,裂纹尖端的动态强度因子和裂纹扩展速度显著下降,穿过节理面后,强度因子又增强;裂纹穿过节理面时,裂纹会沿节理面偏离一段距离后沿原方向继续扩展。实验结果还表明,炮孔与节理间距适宜的情况下,裂纹才会穿过节理继续扩展,间距太小和过大都不利于裂纹的扩展。 相似文献