双相材料模拟的区域分解重心插值配点法 barycentric interpolation collocation method with domain decomposition for simulating bimaterial model.pdfVIP

2013年第1期 玻璃钢／复合材料 双相材料模拟的区域分解重心插值配点法 李淑萍1，王兆清n，唐炳涛2 (1．山东警察学院治安系，山东济南 250101；2．山东建筑大学工程力学研究所，山东济南250101) 摘要：针对双相材料力学性能分析，提出一种极坐标系下的区域分解重心插值配点法。根据材料界面将分析区域划分为 两个计算区域，在每一个计算区域上建立极坐标系下重心插值配点法计算公式。组合两个区域上的计算公式，施加材料界面 条件和边界条件，求解得到双相材料的位移场和应力场。双材料计算模型的数值算例表明了所提方法的计算精度。 关键词：极坐标系；重心插值配点法；区域分解法；弹性力学；双材料模型；无网格方法 文献标识码：A 中图分类号：TB332；TBl2 文章编号：1003—0999(2013)0l一0025—05 1 引 言 性能。在复合材料细观分析方法中，绝大部分采用 复合材料以其优良的力学性能在工程领域得到 直角坐标系计算，我们这里采用极坐标系计算，以更 广泛的应用。复合材料的力学性能分析对复合材料 好地模拟圆形夹杂(孔洞)等材料的不均匀性。根 设计、制造和应用具有重要意义。由于材料的非均 据材料界面将求解区域划分为两个计算区域，在每 匀性和各向异性，复合材料力学性能分析一般采用 一个区域上建立重心插值配点法计算格式，组装两 数值分析方法‘1’2。，有限元方法是最常用的数值分 个区域上的计算矩阵，施加边界条件和界面条件求 析方法旧’4o。为得到问题的高精度解答，有限元方 解得到位移场，进而利用微分矩阵计算得到应力场。 法需要划分稠密的计算网格，增加了计算成本。王 兆清等提出重心有限元方法，采用多边形作为计算 考虑定义在区间[a，b]上的函数“(石)，函数 单元模拟非均质材料∞～7【，极大地简化了网格生成 的工作量。扩展有限元方法(XFEM)采用水平积函 数捕捉材料界面，可以有效地分析非均质材料¨’9』。 9|： Lagrange插值为‘18’1 配点型数值计算方法不需要划分计算单元和数 (1) 值积分，是一种真正意义的无网格方法。根据近似 巾)=荟忐叶／蚤忐 函数的不同可分为谱配点法(拟谱方法)¨0I、径向基 函数(RBF)配点法[11|、重心插值配点法¨1—43和最 重心插值权。记重心Lagrange插值的插值基函 小二乘配点法¨5’161等。拟谱方法采用Chebyshev多 数为： 项式等谱函数插值作为近似函数ⅢJ，是一种高精度 数值计算方法，具有指数收敛精度。径向基函数配 纵小毫／荟盏扣12，…n(2) 点法采用径向基函数插值作为近似函数，其计算精 函数“(算)的重心Lagrange插值可表示为： 度依赖于径向基函数参数的选取。最小二乘配点法 (3) 采用最小二乘拟合作为近似函数，不能直接施加本 u(z)=∑l(z)ui 质边界条件。重心插值配点法采用数值稳定性好、 则函数u(戈)的m阶导数可表示为： 插值精度高的重心Lagrange插值¨8’坶1作为近似函 (