有限元线法用于3D物体的断裂力学性能分析研究.pdfVIP

有限元线法用于3D物体的断裂力学性能分析 胡少伟1王红霞2 (1．南京水利科学研究院，2．上海交通大学) 摘 要：有限元线法(FEMOL)是继有限元法(离散法)之后最近发展起来的半解析方法． 本文提供了三维有限元线法(3DFEMOL)，并且首次详尽地给出了针对断裂力学问题的3D FEMOL分析过程．把这个方法应用于含裂缝三维物体实例分析，标明用较粗糙的网格就能够使结 果得到好的精度．特别是用于受拉试件的应力强度因子计算显示比其它方法更能够得到好的精 度．这标明该方法在含裂缝物体的断裂分析中有巨大的应用前景，并且为FEMOL分析方法用于非 线性材料断裂力学行为分析提供了可能性． 关键词：三维有限元线法(3DFEMOL)；含裂缝3D物体；应力强度因子(SIF) 一、 引言 随着高性能计算机的发展与应用， 在相关文献中涌现了各种各样的数值计算方法。 然而， 绝大多数数值方法仅能提供(局部)近似解，如有限差分法、有限元法、直接能量法、特征函数 扩展法和边界元法。 由于绝大多数开裂计算都涉及几何奇异性问题，一般简单的数值方法都是 无解的。 这篇论文的目的就是把有限元线法(3DFEMOL)用于3D物体的断裂力学性能分析。介于完全 解析法与离散法之间的3DFEMOL在这篇论文中首先被提供。该方法以常微分方程求解器为支撑 软件，利用有限元技术将控制的偏微分方程半离散化为用结线函数表示的常微分方程组，然后调 用高质量的常微分方程求解器，得到用户预先指定误差限的解答。 FEMOL固有的优势是用较为粗糙的网格获得好的精度，而且在较为粗糙的网格上用连续函数 是可以获得部分解答的。至于几何奇异问题，应用一个位移公式能够取得满意精度并使得结果不 再奇异。 二、3DFEMOL基本过程 取代表性问题方程(1)作为例子 V -V2“=f(x，Y，z)in On (1) ”=u(x，Y，z) QD D疗 _2譬=q(’x，Y，z)，，zJ D疗1rⅣrⅣ 其中，三维区域V由平滑边界厶和．R组成，即加=尼u兀，Rn表示外法向方向。 8 (一)区域划分 在FEMOL中，把区域Q划分为一系列四面体单元。单元的形状依赖于单元映射函数。一 般地，四面体单元的全部12条边界允许是曲边或者直边，并且允许在任何区域方便且灵活的建 立映射函数。 厶…表记为整体网格线元数和五Z另…表示整体结点数。 L4 L 13 图1．一个FEMOL网格 (二)单元映射 样的方式可以建立映射单元的次数pp=J，3分别表示线性和立方体单元)。 ·3￡ ／J搿／1 z≈z I± 7 2l 1 312 ．12： ； 气／了 叼 ： 驻； 匕———● 0 l，¨j，一7 ：《 一一