三次样条Hermite插值基构造四边形薄板单元.pdfVIP

中国科学: 数学 2015 年 第45 卷 第9 期: 1523 1536 三次样条 插值基构造四边形薄板单元 献给徐利治教授 华诞 陈娟 李崇君 东北财经大学数学与数量经济学院, 大连 116025; 大连理工大学数学科学学院, 大连 116024 E-mail: chenjdufe@163.com, chongjun@ 收稿日期: 2015-01-21; 接受日期: 2015-03-18; * 通信作者 国家自然科学基金(批准号:和 资助项目 摘要 薄板弯曲单元被广泛地应用到工程问题的有限元计算中. 然而, 由于协调的薄板弯曲位移型单 元要求挠度和转角 (即位移的函数值和导数值) 都是连续的, 导致很难直接构造协调的位移型薄板单 元. 在数学上, 样条是满足一定协调性的分片光滑的多项式, 有限元的形函数可以视为样条函数. 本文 基于三角形面积坐标和B 网方法, 利用三次样条Hermite 插值基重构了两个协调的薄板弯曲单元. 由 于单元形函数是基于四边形构造的, 避免了等参变换, 可以有效地降低网格畸变对计算精度的影响. 关键词 样条插值基 四边形单元 薄板单元 B 网方法 主题分类 65D07, 65N30 引言 在有限元的工程应用中, 板单元的构造一直是研究热点. 协调的薄板弯曲位移型单元要求挠度和 转角, 也即位移的函数值和导数值都是连续的, 因此, 直接构造协调的位移型薄板弯曲单元是非常困 难的 (参见文献 [1, 2]). 例如, 基于三角形的Argyris 元需要21 个参数, 形函数空间是5 次多项式. 为 了减少节点参数和降低形函数空间多项式的次数, 通常采用非协调元, 如6 个参数的Morley 元, 形函 数空间是2 次多项式. 相比三角形单元, 构造协调的四边形板元则更加困难. Sander [3] 和 Fraeijs de Veubeke [4] 曾给出 了一个构造四边形协调单元的方法, 由连接四边形的对角线形成四个三角区域, 分别构造三次位移场 函数. 采用对角线作为斜坐标系来构造基函数, 但它的表达式很复杂, 造成应用上的不便. 对于非协调 元, Batoz 和 Tahar [5] 通过考虑挠度 和转角 各自独立插值, 仅在角点满足离散的 Kirchhoff 假设, 构造了四边形薄板弯曲单元 DKQ (discrete Krichhoff quadrilateral). 它是基于四边形等参单 元Q8 构造的, 优点是形式简单、精度较好, 得到了广泛的应用. DKQ 单元的不足是由于采用等参变换 导致其对网格的抗畸变能力不够高. 许多学者都希望提高DKQ 的精度, 例如, Jeyachandrabose 等人[6] 结合了离散的Kirchhoff 假设和最小二乘法, 导出了一种改进的四边形薄板弯曲单元IDKQ (improved discrete Krichhoff quadrilateral); Chen 和 Cheung [7] 在应变矩阵中加入精化应变矩阵, 得到了精化不 英文引用格式 Chen J, Li C J. The cubic spline Hermite interpolation bases for thin plate bending quadrilateral elements (in Chinese). Sci Sin Math, 2015, 45: 1523–1536, doi: 10.1360/N012015-00035 陈娟等: 三次样条 Hermite 插值基构造四边形薄板单元 协调元 RDKQ (refined discrete Krichhoff quadrilateral). 文献 [8–10] 利用四边形面积坐标方法, 或广 义协调方法, 构造了一些薄板弯曲单元. 等参元是将矩阵单元的插值基函数通过等参变换, 映射为四边形单元的插值基函数, 在有限元计 算时, 需要计算等参变换Jacobi 矩阵的逆, 因此, 当网格畸变时, 会导致等参单元的精度下降 (参见文 [12]