用三维边界元法分析复合材料平板表面裂纹问题.pdfVIP

用三维边界元法分析复合材料平板表面裂纹问题 刘云忠 (中国航天工业总公司四院四十三所北京 100076) 摘要本文建立了一种三维边界元法模式，用于分析复合材料平板表面裂纹闸题。计 算了带小韧带中心裂纹和具有穿透或部分穿透半椭圆型裂纹复合材料平饭的应力强度园子． 探讨了表面裂纹的前、后穿透行为。 关键调三维边界元法表面裂纹复台材料应力强度因子 1引言 为了把线性断裂力学理论成功地运用到实际问题当中去，计算嵌入在复合材料结构 原件中表面裂纹的应力强度因子值是很有实际意义的。在以前，人们大都采用有限元法 来计算嵌入在复合材料结构原件中表面裂纹的应力强度因子。用有限元法分析三维裂纹 问题，在准备输入数据方面需要付出大量的劳动。然而，边界元法在此方面却要优于有 限元法，因为它只需要准备边界数据。但是、边界元法在计算时却要比有限元法赞时，因 为它包含很多矩阵代数，并且它的线性联立方程组的系数矩阵是满矩阵。为了节省CPU 时间，本文建立了一种边界元法模式，称之为三维边界元法，用于表面裂纹问题的分析。 通过分析，本文探讨了表面裂纹的前、后穿透行为。 2建立的三维边界元法模式的特点 在边界元法中，物体的表面被处理成一系列的单元．把节点的位移和表面应力当作 未知数，节点的位移和表面应力用列向量{“)和{9)来定义： …一≯心’……、’ (1) ¨q}一{吼，口2，…’qN} ()和{q}的关系由式(2)给出： [H]{“}一[G]{q) (2) 边界条件引入式(2)可以得到线性联立方程组。 本文建立的三维边界元法模式，有如下特征： (1)假如q，一。的边界条件得到满足，那么，相应地[G]矩阵第i排的元素可以不必 计算。这将节省计算[G]矩阵的时间，因为在实际问题中有许多应力自由边界。 (2)在三维边界元法中所使用的元素类型是六节点、八节点等参单元和九节点拉格 朗日单元(如图l所示)。三角形单元的数值积分通过[-[ammar七节点积分公式来实现， 559 四边形单元的数值积分通过Gauss六 ” 节点积分公式来实现。对于四边形单元， 在单元表面用36个Gauss节点来计算 积分。 △伯琢缚 (3)在裂纹尖端附近，位移和应力 图1在兰维边界元法中的单元 具有／7和1／／7阶奇异性(r是到裂 fa)三角形辱参单元 纹前缘的距离)。这种奇异性可以通过把 (b)四边形等参单元 (c)四边形拉恪朗日单元 中间节点转换到距裂纹前缘1／4节点处 来表示。更连一步，把下面的修正函数叠加到应力奇异单元的形状函数中去： 撕)一去[1一L(I—CZ-)] (3) 这里，L为单元的长度。1／v厂F阶应力奇异性可以通过式(3)来表示。 3带小韧带中心裂纹板的分析 为了研究表面裂纹的前、后穿透行为，计算了带有小韧带中心裂纹复合材料平板的 应力强度因子。文献[1]推荐使用Mindin解作为基本解来计算带有小韧带中心裂绞复合 材料平板的应力强度因子，本文采用Kelvin解进行计算。计算过程中，在裂纹尖端使用 如下应力奇异单元形状…： 女(r)一何 (4) 应力强度因子用式(5)来计算，并且通过最小二乘法进行外推。 K，=苷与√擎(平面应，变) ，(5)