×

直升机仿真中指数-1 DAE的半显式指数Runge-Kutta方法。 (英语) Zbl 1474.65260号

摘要:在本文中,我们建议将指数积分器和半显式Runge-Kutta方法相结合,用于求解微分方程中含有刚性线性部分的指数-1 DAE系统。我们讨论了所得的半显式指数龙格-库塔(HEERK)方法的行为,用于一个简单的数值例子和耦合转子模拟。耦合转子仿真基于模块化软件设计,其中所有子系统都由状态空间形式的ODE建模。通过连接子系统的输入和输出,我们获得了一个指数为1的DAE系统。系统中的大项可以表示为刚性线性部分,其中包括强阻尼或振荡项以及转子叶片离散化的系数(三维梁方程)。我们表明,所提出的HEERK方法可以在合理的时间步长下有效地求解结果系统。

MSC公司:

65升80 微分代数方程的数值方法
65升04 刚性方程的数值方法
65英镑 常微分方程初值问题的数值方法
65升06 常微分方程的多步、Runge-Kutta和外推方法
76-10 流体力学问题的数学建模或模拟
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Arnold,M.,Strehmel,K.,Weiner,R.:索引为1的半显式微分代数方程的半显式Runge-Kutta方法。数字。数学。64(1), 409-431 (1993) ·Zbl 0791.65055号
[2] Beylkin,G.,Keiser,J.M.,Vozovoi,L.:求解非线性偏微分方程的一类新的时间离散格式。J.计算。物理学。147(2), 362-387 (1998) ·Zbl 0924.65089号
[3] Butcher,J.C.:Runge-Kutta方法的历史。申请。数字。数学。20(3), 247-260 (1996) ·Zbl 0858.65073号
[4] Butcher,J.C.,Wanner,G.:Runge-Kutta方法:一些历史注释。申请。数字。数学。22(1-3), 113-151 (1996) ·Zbl 0867.65038号
[5] Condon,M.,Gao,J.,Iserles,A.:关于高振荡半显式DAE的渐近展开解。离散连续。动态。系统。A 36(9),4813-4837(2016)·Zbl 1351.65055号
[6] Cox,S.M.,Matthews,P.C.:刚性系统的指数时间差分。J.计算。物理学。176(2), 430-455 (2002) ·Zbl 1005.65069号
[7] Hairer,E.,Lubich,C.,Roche,M.:用Runge-Kutta方法求解微分代数系统。数学课堂笔记,第1409卷。施普林格,柏林(1989)·Zbl 0683.65050号
[8] Hairer,E.,Lubich,C.,Wanner,G.:几何-数值积分,计算数学中的Springer系列,第31卷。柏林施普林格出版社(2006)·Zbl 1094.65125号
[9] Hairer,E.,Wanner,G.:《求解常微分方程II》,第14卷。施普林格,柏林(1996)·Zbl 0859.65067号
[10] Hochbruck,M.,Ostermann,A.:半线性抛物问题的显式指数Runge-Kutta方法。SIAM J.数字。分析。43(3), 1069-1090 (2005) ·Zbl 1093.65052号
[11] Hochbruck,M.,Ostermann,A.:抛物线问题的指数Runge-Kutta方法。申请。数字。数学。53(2-4), 323-339 (2005) ·Zbl 1070.65099号
[12] Hochbruck,M.,Ostermann,A.:指数积分器。Acta Numer公司。19, 209-286 (2010) ·Zbl 1242.65109号
[13] Houbolt,J.C.,Brooks,G.W.,美国国家航空咨询委员会:扭曲非均匀转子叶片的襟翼弯曲、弦向弯曲和扭转组合运动微分方程。NACA技术报告。国家航空咨询委员会(1957年)
[14] Johnson,W.:旋翼飞机空气力学。剑桥航空航天系列。剑桥大学出版社,剑桥(2013)
[15] Kohlwey,E.:使用指数1微分代数方程系统效率时间积分方法进行直升机仿真。科隆大学硕士论文(2017)
[16] Leis,U.,Wagner,S.:以非定常和粘性效应为重点,实现转子叶片空气载荷的统一表示。第13届欧洲旋翼机论坛,阿尔勒(1987)。网址:http://hdl.handle.net/20.500.11881/2876
[17] Minchev,B.V.:半线性问题的指数积分。卑尔根大学,卑尔根(2004)
[18] Moler,C.,Van Loan,C.:二十五年后,计算矩阵指数的十九种可疑方法。SIAM版本45(1),3-49(2003)·Zbl 1030.65029号
[19] Sanches,L.、Michon,G.、Berlioz,A.、Alazard,D.:各向同性和各向异性多叶转子配置的不稳定区。机械。机器。理论46(8),1054-1065(2011)·Zbl 1372.70064号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。