一种三维非定常线弹性问题的自适应并行算法
发布时间:2020-03-22 16:15
【摘要】:三维非定常线弹性模型广泛存在于固体力学、结构力学和流固耦合等众多领域.此类问题的快速求解是科学与工程计算的研究热点.本文针对一种三维非定常线弹性模型,在拉氏(Lagrangian)网格下给出了时空方向分别采用有限差分和线性有限元离散的全离散隐格式及相应的求解算法,数值实验表明该格式的有限元误差函数在离散L_2范数下具有O(τ+h_2)的饱和误差阶.接着,给出时空步长对上述离散格式的系数矩阵条件数的影响分析,进而设计了全离散隐格式的自适应求解算法,数值实验表明该方法的求解效率优于常用的CG法和AMG-CG法.最后,在MPI并行环境中,为上述自适应算法设计出相应的并行求解算法;进一步结合不同解法器的并行可扩展性,设计了一种基于时空尺度关系和进程数的自适应并行算法及相应的求解模块,数值实验表明改进的自适应并行算法具有更好求解效率,如对于初始时刻的空间剖分步长为1/64和进程数为64的情况下,新的自适应并行算法的求解效率相对改进前提升了30.66%.
【图文】:
t°邋tl邋t2逦tMA邋tM逡逑图3.2时间剖分示意图.逡逑引入速度函数逡逑v(t)邋:=邋v(x(t),t)邋:=邋Dtu(t),邋Vie邋(0,T].逦(3.9)逡逑对任给的时刻im+1邋(m邋=邋0,…,,M邋-邋1),设己知严时刻的速度和位移逡逑u(f^,下面给出Atu⑷在tm+1时刻的有限差分近似公式.逡逑为此,首先给出近似公式逡逑Dttn(tm+1)邋=邋Dtv(tm+1)逡逑?邋—逦^ ̄V(",(向后差商,局部截断误差为0(Af))逡逑(3.10)逡逑接着导出v(f"+1)的近似公式,为此给出逡逑v(^+^)邋?逦y(tm+n逦(平均速度)逦(3.11)逡逑V(严+!)邋?(中心差商)逦(3.12)逡逑由(3.11),(3.12)近似可得逡逑v(r^)邋=邋2u(mm)邋_邋v(n
图3.4邋1个六面体单元分裂为6个四面体单元的剖分规则示意图逡逑
【学位授予单位】:湘潭大学
【学位级别】:硕士
【学位授予年份】:2018
【分类号】:O241.8
本文编号:2595294
【图文】:
t°邋tl邋t2逦tMA邋tM逡逑图3.2时间剖分示意图.逡逑引入速度函数逡逑v(t)邋:=邋v(x(t),t)邋:=邋Dtu(t),邋Vie邋(0,T].逦(3.9)逡逑对任给的时刻im+1邋(m邋=邋0,…,,M邋-邋1),设己知严时刻的速度和位移逡逑u(f^,下面给出Atu⑷在tm+1时刻的有限差分近似公式.逡逑为此,首先给出近似公式逡逑Dttn(tm+1)邋=邋Dtv(tm+1)逡逑?邋—逦^ ̄V(",(向后差商,局部截断误差为0(Af))逡逑(3.10)逡逑接着导出v(f"+1)的近似公式,为此给出逡逑v(^+^)邋?逦y(tm+n逦(平均速度)逦(3.11)逡逑V(严+!)邋?(中心差商)逦(3.12)逡逑由(3.11),(3.12)近似可得逡逑v(r^)邋=邋2u(mm)邋_邋v(n
图3.4邋1个六面体单元分裂为6个四面体单元的剖分规则示意图逡逑
【学位授予单位】:湘潭大学
【学位级别】:硕士
【学位授予年份】:2018
【分类号】:O241.8
【参考文献】
相关期刊论文 前3条
1 肖映雄;周志阳;舒适;;几类典型网格下三维弹性问题的代数多层网格法[J];工程力学;2011年06期
2 白象忠;郝亚娟;;非线性流体弹性力学研究进展[J];力学进展;2008年05期
3 张平;肖映雄;舒适;;代数多层网格法及其在固体力学计算中的应用研究[J];湘潭大学自然科学学报;2008年03期
相关博士学位论文 前1条
1 李政;精细油藏数值模拟中的高效求解器研究[D];昆明理工大学;2017年
相关硕士学位论文 前7条
1 李奇;一种二维非定常线弹性问题的自适应并行算法[D];湘潭大学;2017年
2 龚能;智能电网设备PMU的优化配置研究[D];湘潭大学;2016年
3 王锦;一种求解三维Neumann边界条件线弹性问题线性有限元方程的高效预条件子[D];湘潭大学;2016年
4 岳孝强;几种并行AMG法及其在辐射扩散问题中的应用[D];湘潭大学;2012年
5 梁文涛;一种求解三维弹性问题有限元方程的并行DDM预条件子[D];湘潭大学;2010年
6 曹放;一种基于HYPRE的高次Lagrange有限元方程的并行代数多层网格法[D];湘潭大学;2008年
7 聂玮芳;基于集群的并行Splatting算法关键技术的研究[D];天津大学;2006年
本文编号:2595294
本文链接:https://www.wllwen.com/kejilunwen/yysx/2595294.html