基于SBFEM的多裂纹问题断裂分析.pdfVIP

第48卷第3期 大连理工大学学报 V01．48．、No．3 2 00 2 00 8 8年5月 of of May JournalDalian UniversityTechnology 基于SBFEM的多裂纹问题断裂分析 刘钧玉， 林 皋。， 杜建国 (大连理工大学土木水利学院，辽宁大连 116024) 摘要：多裂纹问题中的应力强度因子是断裂力学中需要计算的重要参数．在子结构法思想 的基础上利用比例边界有限元法计算了弹性多裂纹问题的I型裂纹应力强度因子．对于多裂 纹的弹性问题根据裂纹的数目确定相似中心的数日。在每一个子块内保持比例边界有限元法 的特点．利用该数值技巧可以求解任意多裂纹问题的应力强度因子．数值算例表明该方法是 有效且精确的．最后给出了正交各向异性材料双边非对称裂纹问题的计算结果，进而推广了 比例边界有限元法的应用范围． 关键词：断裂力学≯应力强度因子；多裂纹问题；比例边界有限元法；正交各向异 性材料 中图分类号：0346．1；0241文献标志码：A 0 引 言 应力强度因子‘妇及裂纹和多材料角的T一应力和 平面奇异应力场的高阶奇异条件[51．本文基于子 具有嵌入几何奇异性(如裂纹、不同材料接 结构法思想利用比例边界有限元法求解弹性多裂 点、缺口等)弹性体的应力分析是工程应用中一个 纹问题的应力强度因子，并与解析解和文献中可 十分重要的问题；许多复杂的工程要求精确地确 获得的计算结果进行对比． 定奇异参数如弹性断裂力学中的应力强度因子． 由于解析方法只能求解规则形状、简单荷载作用的 1 比例边界有限元法的基本方程 问题，大部分实际问题都需要利用数值方法进行求 比例边界有限元法将在频率域内以位移表达 解．最普遍的求解应力强度因子的方法是有限元 的弹性动力偏微分方程(1)转换为常微分方程形 法、边界元法、比例边界有限元法、边界配置法等， 式比例边界有限元方程(2)，可以解析地求解． 此外还可用有限元法和无限元法联合求解[1】． (1) LT口+p+cE，2pU=0 比例边界有限元法集合了有限元法和边界元 法的特点，该方法只需用有限元离散部分边界，从 (2) 而实现了将问题降低一维的目的，在计算应力强 上述式(1)中部分参数如图1所示．式中和图中y 度因子时，对于裂纹面通过相似中心的裂纹不需 为研究区域≯ID为质量密度；P为体荷载；口为应力； 要离散裂纹面，减少了计算量．此外，比例边界有 L为微分算子fs为维数，S。为位移边界；S。为荷载 限