C1连续型广义有限元格式

• 作者: 田荣
• 来源期刊: 力学学报
• 学科: 工学
• 关键词: 单位分解插值 广义有限元 无额外自由度 C1连续
• 年，卷（期）: 2019,（1）
• 所属期刊栏目: 生物、工程及交叉力学
• 页码范围: 263-277
• 页数: 15页
• 分类号: TB115|O346.1
• 字数: 5721字
• 语种: 中文
• DOI: 10.6052/0459-1879-18-188

摘要

C1连续,即一阶导数连续.C1连续型插值格式具有同时适用于离散PDE的弱形式与强形式的优点 ——即一种插值格式可以在使用PDE弱形式还是强形式之间做出选择,从而构造出更加高效的数值方法.由于单位分解广义有限元方法(PUFEM,Babuˇska and Melenk(1997)),允许用户根据局部解的特征自定义任意高阶局部近似,具有精度高、程序实现与传统有限元相容性好的特点而受到广泛关注.但是,其总体近似函数的光滑性是由其所采用的单位分解函数 —— 一般为标准有限元形函数 —— 的光滑性所决定,因此多为C0连续.如何在C0连续标准有限元形函数的基础上,构造出满足C1连续的总体近似函数,是一个仍未解决的问题.本文在作者前期研究的无额外自由度的单位分解插值格式的基础上,仅基于C0标准有限元形函数,构造出至少C1连续的无额外自由度单位分解格式.针对Poisson方程,讨论了该格式对PDE弱形式与强形式的离散.测试结果表明,方法可以同时用于弱形式与强形式的数值求解,而且可以在不改变网格和自由度数的前提下,获得高阶收敛.使用该插值格式的条件是:网格须是直角坐标网格(不要求均匀).该插值格式可以同时用于流体力学问题和使用欧拉背景网格求解动量方程的固体力学方法,如材料物质点法(material point method).对于强形式的欧拉网格求解,该插值格式与"差分"不同之处在于,它具有有限元一样的在任意点处进行"插值"的特点.对于弱形式的积分求解,由于该插值格式具有导数连续性,可以允许积分网格独立于插值网格.这一特点将使得弱形式的数值积分的实施更加灵活方便.

作者信息

序号	姓名	单位	发文数	被引次数	H指数	G指数
1	田荣	中国工程物理研究院高性能数值模拟软件中心	2	4	1.0	2.0