1、力学专业毕业论文 精品论文 L-稳定格式求解结构动力学方程和多体系统动力学方程关键词：结构动力学方程 多体系统动力学方程 数值求解 Runge-Kutta 方法 L-稳定格式 微分-代数方程 常微分方程摘要：航天工程实际问题中包含大量结构和机构系统，如卫星包含天线和帆板等柔性结构，天线和帆板展开过程又是机构运动。研究结构的动力学行为，通常是对结构动力学方程进行研究；研究机构的运动，通常利用多体系统动力学的建模理论，建立多体系统动力学方程。空间离散的结构动力学方程为常微分方程的初值问题，多体系统动力学方程为微分-代数方程的初值问题。本文的主要目的是试图研究这两类方程的数值求解。 采用 Runge
2、-Kutta 方法求解结构动力学方程。针对增量形式的结构动力学方程，推导了 Runge-Kutta 方法的迭代形式，研究了降低计算量的两种办法，发现了一些传统格式与 Runge-Kutta方法之间的关系。基于单自由度自由振动方程，推导了 Runge-Kutta 方法求解结构动力学方程的逼近算子，应用逼近算子的谱半径可以分析 Runge-Kutta 方法求解二阶常微分方程的稳定性。通过谱半径图分析出各种格式的数值阻尼的强弱，发现高精度的 L-稳定格式具有理想的数值阻尼，在精确求解低频振动的同时很好的抑制高频震荡。通过算例研究，发现较小的物理阻尼在精确解上能快速的耗散高频震荡，对数值计算却影响很小
3、。从而得出结论：在数值计算中通过增加物理阻尼抑制高频震荡是不可行，而采用有数值阻尼的算法才是可行的。 应用适合求解高指标微分.代数方程的求解器，求解高指标形式的多体系统动力学方程。发现采用高指标形式的多体系统动力学方程，使用隐式形式微分方程需要引入新的求解变量少，而且能保持原系统矩阵的稀疏特性。但是存在两个缺点：(1)加速度变量和拉格朗日乘子的求解精度低于位移变量、速度变量的求解精度；(2)对协调一致的初始条件很依赖。针对这两个缺点，提出了对加速度变量和拉格朗日乘子修正的办法。采用修正能提高加速度变量和拉格朗日乘子的求解精度，提高整个计算过程中的收敛性。数值算例表明，采用修正能降低求解的计算量
4、，特别对于容许误差很小的情况，计算量降低更显著。 基于微分方程精确解的 Taylor 展式，提出了一种新的构造块格式方法；并采用指数函数的 Pade 逼近作为块格式的稳定性函数，构造出 L-稳定的块格式。实际上，块格式也是一类 Runge-Kutta 方法，从而通过 Runge-Kutta 方法的精度条件，分析了 L-稳定块格式的精度。针对隐式形式微分方程，推导了块格式的并行迭代和局部误差估计方法，采用变步长实现了 6 阶精度 L-稳定格式，编制了求解器 PBIDE(ParallelBlocksolverforImplicitDifferentialEquations)。 针对 Runge-K
5、utta 方法在求解微分.代数方程中出现精度降低的情况，提出了多步块格式。分析了多步块格式稳定性和精度条件，证明其为刚性精确的，且级精度与格式精度相同。数值算例表明在求解指标-2 的微分-代数方程不会出现精度降低现象。构造出了 2 级 3 阶精度的多步块格式，此格式的稳定性函数为指数函数的 4 阶精度广义 Pade 量近，且是 L-稳定的；构造出了 3 级 4 阶精度的多步块格式，此格式的稳定性函数与 Runge-Kutta 的稳定性函数一样，为指数函数的有理分式逼近，也是 L-稳定格式。推导了多步块格式的 Nordsieek 表达。基于 2 级 3 阶精度的多步块格式，对隐式形式微分方程，采
6、用变步长进行了实现，编制了求解器MBIDE(MultistepBlocksolverforImplicitDifferentialEquations)。 对求解器 PBIDE 和 MBIDE 进行测试，与现在使用广泛的求解器进行比较，比较了它们的运行精度和计算量。选用结构动力学的算例、多体系统动力学的算例和几个其它类型的算例对本文的两个求解器进行了测试，测试的结果表明这两个求解器是效率很高的。正文内容航天工程实际问题中包含大量结构和机构系统，如卫星包含天线和帆板等柔性结构，天线和帆板展开过程又是机构运动。研究结构的动力学行为，通常是对结构动力学方程进行研究；研究机构的运动，通常利用多体系统动力
7、学的建模理论，建立多体系统动力学方程。空间离散的结构动力学方程为常微分方程的初值问题，多体系统动力学方程为微分-代数方程的初值问题。本文的主要目的是试图研究这两类方程的数值求解。 采用 Runge-Kutta 方法求解结构动力学方程。针对增量形式的结构动力学方程，推导了 Runge-Kutta 方法的迭代形式，研究了降低计算量的两种办法，发现了一些传统格式与 Runge-Kutta 方法之间的关系。基于单自由度自由振动方程，推导了 Runge-Kutta 方法求解结构动力学方程的逼近算子，应用逼近算子的谱半径可以分析 Runge-Kutta 方法求解二阶常微分方程的稳定性。通过谱半径图分析出各
8、种格式的数值阻尼的强弱，发现高精度的 L-稳定格式具有理想的数值阻尼，在精确求解低频振动的同时很好的抑制高频震荡。通过算例研究，发现较小的物理阻尼在精确解上能快速的耗散高频震荡，对数值计算却影响很小。从而得出结论：在数值计算中通过增加物理阻尼抑制高频震荡是不可行，而采用有数值阻尼的算法才是可行的。应用适合求解高指标微分.代数方程的求解器，求解高指标形式的多体系统动力学方程。发现采用高指标形式的多体系统动力学方程，使用隐式形式微分方程需要引入新的求解变量少，而且能保持原系统矩阵的稀疏特性。但是存在两个缺点：(1)加速度变量和拉格朗日乘子的求解精度低于位移变量、速度变量的求解精度；(2)对协调一致
9、的初始条件很依赖。针对这两个缺点，提出了对加速度变量和拉格朗日乘子修正的办法。采用修正能提高加速度变量和拉格朗日乘子的求解精度，提高整个计算过程中的收敛性。数值算例表明，采用修正能降低求解的计算量，特别对于容许误差很小的情况，计算量降低更显著。 基于微分方程精确解的 Taylor 展式，提出了一种新的构造块格式方法；并采用指数函数的 Pade 逼近作为块格式的稳定性函数，构造出 L-稳定的块格式。实际上，块格式也是一类 Runge-Kutta 方法，从而通过 Runge-Kutta 方法的精度条件，分析了 L-稳定块格式的精度。针对隐式形式微分方程，推导了块格式的并行迭代和局部误差估计方法，采
10、用变步长实现了 6 阶精度 L-稳定格式，编制了求解器 PBIDE(ParallelBlocksolverforImplicitDifferentialEquations)。 针对 Runge-Kutta 方法在求解微分.代数方程中出现精度降低的情况，提出了多步块格式。分析了多步块格式稳定性和精度条件，证明其为刚性精确的，且级精度与格式精度相同。数值算例表明在求解指标-2 的微分-代数方程不会出现精度降低现象。构造出了 2 级 3 阶精度的多步块格式，此格式的稳定性函数为指数函数的 4 阶精度广义 Pade 量近，且是 L-稳定的；构造出了 3 级 4 阶精度的多步块格式，此格式的稳定性函数与
11、 Runge-Kutta 的稳定性函数一样，为指数函数的有理分式逼近，也是 L-稳定格式。推导了多步块格式的 Nordsieek 表达。基于 2 级 3 阶精度的多步块格式，对隐式形式微分方程，采用变步长进行了实现，编制了求解器MBIDE(MultistepBlocksolverforImplicitDifferentialEquations)。 对求解器 PBIDE 和 MBIDE 进行测试，与现在使用广泛的求解器进行比较，比较了它们的运行精度和计算量。选用结构动力学的算例、多体系统动力学的算例和几个其它类型的算例对本文的两个求解器进行了测试，测试的结果表明这两个求解器是效率很高的。航天工程
12、实际问题中包含大量结构和机构系统，如卫星包含天线和帆板等柔性结构，天线和帆板展开过程又是机构运动。研究结构的动力学行为，通常是对结构动力学方程进行研究；研究机构的运动，通常利用多体系统动力学的建模理论，建立多体系统动力学方程。空间离散的结构动力学方程为常微分方程的初值问题，多体系统动力学方程为微分-代数方程的初值问题。本文的主要目的是试图研究这两类方程的数值求解。 采用 Runge-Kutta 方法求解结构动力学方程。针对增量形式的结构动力学方程，推导了 Runge-Kutta 方法的迭代形式，研究了降低计算量的两种办法，发现了一些传统格式与 Runge-Kutta 方法之间的关系。基于单自由
13、度自由振动方程，推导了 Runge-Kutta 方法求解结构动力学方程的逼近算子，应用逼近算子的谱半径可以分析 Runge-Kutta 方法求解二阶常微分方程的稳定性。通过谱半径图分析出各种格式的数值阻尼的强弱，发现高精度的 L-稳定格式具有理想的数值阻尼，在精确求解低频振动的同时很好的抑制高频震荡。通过算例研究，发现较小的物理阻尼在精确解上能快速的耗散高频震荡，对数值计算却影响很小。从而得出结论：在数值计算中通过增加物理阻尼抑制高频震荡是不可行，而采用有数值阻尼的算法才是可行的。 应用适合求解高指标微分.代数方程的求解器，求解高指标形式的多体系统动力学方程。发现采用高指标形式的多体系统动力学
14、方程，使用隐式形式微分方程需要引入新的求解变量少，而且能保持原系统矩阵的稀疏特性。但是存在两个缺点：(1)加速度变量和拉格朗日乘子的求解精度低于位移变量、速度变量的求解精度；(2)对协调一致的初始条件很依赖。针对这两个缺点，提出了对加速度变量和拉格朗日乘子修正的办法。采用修正能提高加速度变量和拉格朗日乘子的求解精度，提高整个计算过程中的收敛性。数值算例表明，采用修正能降低求解的计算量，特别对于容许误差很小的情况，计算量降低更显著。 基于微分方程精确解的 Taylor 展式，提出了一种新的构造块格式方法；并采用指数函数的 Pade 逼近作为块格式的稳定性函数，构造出 L-稳定的块格式。实际上，块
15、格式也是一类 Runge-Kutta 方法，从而通过 Runge-Kutta 方法的精度条件，分析了 L-稳定块格式的精度。针对隐式形式微分方程，推导了块格式的并行迭代和局部误差估计方法，采用变步长实现了 6 阶精度 L-稳定格式，编制了求解器PBIDE(ParallelBlocksolverforImplicitDifferentialEquations)。 针对Runge-Kutta 方法在求解微分.代数方程中出现精度降低的情况，提出了多步块格式。分析了多步块格式稳定性和精度条件，证明其为刚性精确的，且级精度与格式精度相同。数值算例表明在求解指标-2 的微分-代数方程不会出现精度降低现象。
16、构造出了 2 级 3 阶精度的多步块格式，此格式的稳定性函数为指数函数的 4 阶精度广义 Pade 量近，且是 L-稳定的；构造出了 3 级 4 阶精度的多步块格式，此格式的稳定性函数与 Runge-Kutta 的稳定性函数一样，为指数函数的有理分式逼近，也是 L-稳定格式。推导了多步块格式的 Nordsieek 表达。基于 2 级 3 阶精度的多步块格式，对隐式形式微分方程，采用变步长进行了实现，编制了求解器MBIDE(MultistepBlocksolverforImplicitDifferentialEquations)。 对求解器 PBIDE 和 MBIDE 进行测试，与现在使用广泛的
17、求解器进行比较，比较了它们的运行精度和计算量。选用结构动力学的算例、多体系统动力学的算例和几个其它类型的算例对本文的两个求解器进行了测试，测试的结果表明这两个求解器是效率很高的。航天工程实际问题中包含大量结构和机构系统，如卫星包含天线和帆板等柔性结构，天线和帆板展开过程又是机构运动。研究结构的动力学行为，通常是对结构动力学方程进行研究；研究机构的运动，通常利用多体系统动力学的建模理论，建立多体系统动力学方程。空间离散的结构动力学方程为常微分方程的初值问题，多体系统动力学方程为微分-代数方程的初值问题。本文的主要目的是试图研究这两类方程的数值求解。 采用 Runge-Kutta 方法求解结构动力
18、学方程。针对增量形式的结构动力学方程，推导了 Runge-Kutta 方法的迭代形式，研究了降低计算量的两种办法，发现了一些传统格式与 Runge-Kutta 方法之间的关系。基于单自由度自由振动方程，推导了 Runge-Kutta 方法求解结构动力学方程的逼近算子，应用逼近算子的谱半径可以分析 Runge-Kutta 方法求解二阶常微分方程的稳定性。通过谱半径图分析出各种格式的数值阻尼的强弱，发现高精度的 L-稳定格式具有理想的数值阻尼，在精确求解低频振动的同时很好的抑制高频震荡。通过算例研究，发现较小的物理阻尼在精确解上能快速的耗散高频震荡，对数值计算却影响很小。从而得出结论：在数值计算中
19、通过增加物理阻尼抑制高频震荡是不可行，而采用有数值阻尼的算法才是可行的。 应用适合求解高指标微分.代数方程的求解器，求解高指标形式的多体系统动力学方程。发现采用高指标形式的多体系统动力学方程，使用隐式形式微分方程需要引入新的求解变量少，而且能保持原系统矩阵的稀疏特性。但是存在两个缺点：(1)加速度变量和拉格朗日乘子的求解精度低于位移变量、速度变量的求解精度；(2)对协调一致的初始条件很依赖。针对这两个缺点，提出了对加速度变量和拉格朗日乘子修正的办法。采用修正能提高加速度变量和拉格朗日乘子的求解精度，提高整个计算过程中的收敛性。数值算例表明，采用修正能降低求解的计算量，特别对于容许误差很小的情况
20、，计算量降低更显著。 基于微分方程精确解的 Taylor 展式，提出了一种新的构造块格式方法；并采用指数函数的 Pade 逼近作为块格式的稳定性函数，构造出 L-稳定的块格式。实际上，块格式也是一类 Runge-Kutta 方法，从而通过 Runge-Kutta 方法的精度条件，分析了 L-稳定块格式的精度。针对隐式形式微分方程，推导了块格式的并行迭代和局部误差估计方法，采用变步长实现了 6 阶精度 L-稳定格式，编制了求解器PBIDE(ParallelBlocksolverforImplicitDifferentialEquations)。 针对Runge-Kutta 方法在求解微分.代数方
21、程中出现精度降低的情况，提出了多步块格式。分析了多步块格式稳定性和精度条件，证明其为刚性精确的，且级精度与格式精度相同。数值算例表明在求解指标-2 的微分-代数方程不会出现精度降低现象。构造出了 2 级 3 阶精度的多步块格式，此格式的稳定性函数为指数函数的 4 阶精度广义 Pade 量近，且是 L-稳定的；构造出了 3 级 4 阶精度的多步块格式，此格式的稳定性函数与 Runge-Kutta 的稳定性函数一样，为指数函数的有理分式逼近，也是 L-稳定格式。推导了多步块格式的 Nordsieek 表达。基于 2 级 3 阶精度的多步块格式，对隐式形式微分方程，采用变步长进行了实现，编制了求解器
22、MBIDE(MultistepBlocksolverforImplicitDifferentialEquations)。 对求解器 PBIDE 和 MBIDE 进行测试，与现在使用广泛的求解器进行比较，比较了它们的运行精度和计算量。选用结构动力学的算例、多体系统动力学的算例和几个其它类型的算例对本文的两个求解器进行了测试，测试的结果表明这两个求解器是效率很高的。航天工程实际问题中包含大量结构和机构系统，如卫星包含天线和帆板等柔性结构，天线和帆板展开过程又是机构运动。研究结构的动力学行为，通常是对结构动力学方程进行研究；研究机构的运动，通常利用多体系统动力学的建模理论，建立多体系统动力学方程。空
23、间离散的结构动力学方程为常微分方程的初值问题，多体系统动力学方程为微分-代数方程的初值问题。本文的主要目的是试图研究这两类方程的数值求解。 采用 Runge-Kutta 方法求解结构动力学方程。针对增量形式的结构动力学方程，推导了 Runge-Kutta 方法的迭代形式，研究了降低计算量的两种办法，发现了一些传统格式与 Runge-Kutta 方法之间的关系。基于单自由度自由振动方程，推导了 Runge-Kutta 方法求解结构动力学方程的逼近算子，应用逼近算子的谱半径可以分析 Runge-Kutta 方法求解二阶常微分方程的稳定性。通过谱半径图分析出各种格式的数值阻尼的强弱，发现高精度的 L
24、-稳定格式具有理想的数值阻尼，在精确求解低频振动的同时很好的抑制高频震荡。通过算例研究，发现较小的物理阻尼在精确解上能快速的耗散高频震荡，对数值计算却影响很小。从而得出结论：在数值计算中通过增加物理阻尼抑制高频震荡是不可行，而采用有数值阻尼的算法才是可行的。 应用适合求解高指标微分.代数方程的求解器，求解高指标形式的多体系统动力学方程。发现采用高指标形式的多体系统动力学方程，使用隐式形式微分方程需要引入新的求解变量少，而且能保持原系统矩阵的稀疏特性。但是存在两个缺点：(1)加速度变量和拉格朗日乘子的求解精度低于位移变量、速度变量的求解精度；(2)对协调一致的初始条件很依赖。针对这两个缺点，提出
25、了对加速度变量和拉格朗日乘子修正的办法。采用修正能提高加速度变量和拉格朗日乘子的求解精度，提高整个计算过程中的收敛性。数值算例表明，采用修正能降低求解的计算量，特别对于容许误差很小的情况，计算量降低更显著。 基于微分方程精确解的 Taylor 展式，提出了一种新的构造块格式方法；并采用指数函数的 Pade 逼近作为块格式的稳定性函数，构造出 L-稳定的块格式。实际上，块格式也是一类 Runge-Kutta 方法，从而通过 Runge-Kutta 方法的精度条件，分析了 L-稳定块格式的精度。针对隐式形式微分方程，推导了块格式的并行迭代和局部误差估计方法，采用变步长实现了 6 阶精度 L-稳定格
26、式，编制了求解器PBIDE(ParallelBlocksolverforImplicitDifferentialEquations)。 针对Runge-Kutta 方法在求解微分.代数方程中出现精度降低的情况，提出了多步块格式。分析了多步块格式稳定性和精度条件，证明其为刚性精确的，且级精度与格式精度相同。数值算例表明在求解指标-2 的微分-代数方程不会出现精度降低现象。构造出了 2 级 3 阶精度的多步块格式，此格式的稳定性函数为指数函数的 4 阶精度广义 Pade 量近，且是 L-稳定的；构造出了 3 级 4 阶精度的多步块格式，此格式的稳定性函数与 Runge-Kutta 的稳定性函数一样
27、，为指数函数的有理分式逼近，也是 L-稳定格式。推导了多步块格式的 Nordsieek 表达。基于 2 级 3 阶精度的多步块格式，对隐式形式微分方程，采用变步长进行了实现，编制了求解器MBIDE(MultistepBlocksolverforImplicitDifferentialEquations)。 对求解器 PBIDE 和 MBIDE 进行测试，与现在使用广泛的求解器进行比较，比较了它们的运行精度和计算量。选用结构动力学的算例、多体系统动力学的算例和几个其它类型的算例对本文的两个求解器进行了测试，测试的结果表明这两个求解器是效率很高的。航天工程实际问题中包含大量结构和机构系统，如卫星包
28、含天线和帆板等柔性结构，天线和帆板展开过程又是机构运动。研究结构的动力学行为，通常是对结构动力学方程进行研究；研究机构的运动，通常利用多体系统动力学的建模理论，建立多体系统动力学方程。空间离散的结构动力学方程为常微分方程的初值问题，多体系统动力学方程为微分-代数方程的初值问题。本文的主要目的是试图研究这两类方程的数值求解。 采用 Runge-Kutta 方法求解结构动力学方程。针对增量形式的结构动力学方程，推导了 Runge-Kutta 方法的迭代形式，研究了降低计算量的两种办法，发现了一些传统格式与 Runge-Kutta 方法之间的关系。基于单自由度自由振动方程，推导了 Runge-Kut
29、ta 方法求解结构动力学方程的逼近算子，应用逼近算子的谱半径可以分析 Runge-Kutta 方法求解二阶常微分方程的稳定性。通过谱半径图分析出各种格式的数值阻尼的强弱，发现高精度的 L-稳定格式具有理想的数值阻尼，在精确求解低频振动的同时很好的抑制高频震荡。通过算例研究，发现较小的物理阻尼在精确解上能快速的耗散高频震荡，对数值计算却影响很小。从而得出结论：在数值计算中通过增加物理阻尼抑制高频震荡是不可行，而采用有数值阻尼的算法才是可行的。 应用适合求解高指标微分.代数方程的求解器，求解高指标形式的多体系统动力学方程。发现采用高指标形式的多体系统动力学方程，使用隐式形式微分方程需要引入新的求解
30、变量少，而且能保持原系统矩阵的稀疏特性。但是存在两个缺点：(1)加速度变量和拉格朗日乘子的求解精度低于位移变量、速度变量的求解精度；(2)对协调一致的初始条件很依赖。针对这两个缺点，提出了对加速度变量和拉格朗日乘子修正的办法。采用修正能提高加速度变量和拉格朗日乘子的求解精度，提高整个计算过程中的收敛性。数值算例表明，采用修正能降低求解的计算量，特别对于容许误差很小的情况，计算量降低更显著。 基于微分方程精确解的 Taylor 展式，提出了一种新的构造块格式方法；并采用指数函数的 Pade 逼近作为块格式的稳定性函数，构造出 L-稳定的块格式。实际上，块格式也是一类 Runge-Kutta 方法
31、，从而通过 Runge-Kutta 方法的精度条件，分析了 L-稳定块格式的精度。针对隐式形式微分方程，推导了块格式的并行迭代和局部误差估计方法，采用变步长实现了 6 阶精度 L-稳定格式，编制了求解器PBIDE(ParallelBlocksolverforImplicitDifferentialEquations)。 针对Runge-Kutta 方法在求解微分.代数方程中出现精度降低的情况，提出了多步块格式。分析了多步块格式稳定性和精度条件，证明其为刚性精确的，且级精度与格式精度相同。数值算例表明在求解指标-2 的微分-代数方程不会出现精度降低现象。构造出了 2 级 3 阶精度的多步块格式，
32、此格式的稳定性函数为指数函数的 4 阶精度广义 Pade 量近，且是 L-稳定的；构造出了 3 级 4 阶精度的多步块格式，此格式的稳定性函数与 Runge-Kutta 的稳定性函数一样，为指数函数的有理分式逼近，也是 L-稳定格式。推导了多步块格式的 Nordsieek 表达。基于 2 级 3 阶精度的多步块格式，对隐式形式微分方程，采用变步长进行了实现，编制了求解器MBIDE(MultistepBlocksolverforImplicitDifferentialEquations)。 对求解器 PBIDE 和 MBIDE 进行测试，与现在使用广泛的求解器进行比较，比较了它们的运行精度和计算
33、量。选用结构动力学的算例、多体系统动力学的算例和几个其它类型的算例对本文的两个求解器进行了测试，测试的结果表明这两个求解器是效率很高的。航天工程实际问题中包含大量结构和机构系统，如卫星包含天线和帆板等柔性结构，天线和帆板展开过程又是机构运动。研究结构的动力学行为，通常是对结构动力学方程进行研究；研究机构的运动，通常利用多体系统动力学的建模理论，建立多体系统动力学方程。空间离散的结构动力学方程为常微分方程的初值问题，多体系统动力学方程为微分-代数方程的初值问题。本文的主要目的是试图研究这两类方程的数值求解。 采用 Runge-Kutta 方法求解结构动力学方程。针对增量形式的结构动力学方程，推导
34、了 Runge-Kutta 方法的迭代形式，研究了降低计算量的两种办法，发现了一些传统格式与 Runge-Kutta 方法之间的关系。基于单自由度自由振动方程，推导了 Runge-Kutta 方法求解结构动力学方程的逼近算子，应用逼近算子的谱半径可以分析 Runge-Kutta 方法求解二阶常微分方程的稳定性。通过谱半径图分析出各种格式的数值阻尼的强弱，发现高精度的 L-稳定格式具有理想的数值阻尼，在精确求解低频振动的同时很好的抑制高频震荡。通过算例研究，发现较小的物理阻尼在精确解上能快速的耗散高频震荡，对数值计算却影响很小。从而得出结论：在数值计算中通过增加物理阻尼抑制高频震荡是不可行，而采
35、用有数值阻尼的算法才是可行的。 应用适合求解高指标微分.代数方程的求解器，求解高指标形式的多体系统动力学方程。发现采用高指标形式的多体系统动力学方程，使用隐式形式微分方程需要引入新的求解变量少，而且能保持原系统矩阵的稀疏特性。但是存在两个缺点：(1)加速度变量和拉格朗日乘子的求解精度低于位移变量、速度变量的求解精度；(2)对协调一致的初始条件很依赖。针对这两个缺点，提出了对加速度变量和拉格朗日乘子修正的办法。采用修正能提高加速度变量和拉格朗日乘子的求解精度，提高整个计算过程中的收敛性。数值算例表明，采用修正能降低求解的计算量，特别对于容许误差很小的情况，计算量降低更显著。 基于微分方程精确解的
36、 Taylor 展式，提出了一种新的构造块格式方法；并采用指数函数的 Pade 逼近作为块格式的稳定性函数，构造出 L-稳定的块格式。实际上，块格式也是一类 Runge-Kutta 方法，从而通过 Runge-Kutta 方法的精度条件，分析了 L-稳定块格式的精度。针对隐式形式微分方程，推导了块格式的并行迭代和局部误差估计方法，采用变步长实现了 6 阶精度 L-稳定格式，编制了求解器PBIDE(ParallelBlocksolverforImplicitDifferentialEquations)。 针对Runge-Kutta 方法在求解微分.代数方程中出现精度降低的情况，提出了多步块格式。
37、分析了多步块格式稳定性和精度条件，证明其为刚性精确的，且级精度与格式精度相同。数值算例表明在求解指标-2 的微分-代数方程不会出现精度降低现象。构造出了 2 级 3 阶精度的多步块格式，此格式的稳定性函数为指数函数的 4 阶精度广义 Pade 量近，且是 L-稳定的；构造出了 3 级 4 阶精度的多步块格式，此格式的稳定性函数与 Runge-Kutta 的稳定性函数一样，为指数函数的有理分式逼近，也是 L-稳定格式。推导了多步块格式的 Nordsieek 表达。基于 2 级 3 阶精度的多步块格式，对隐式形式微分方程，采用变步长进行了实现，编制了求解器MBIDE(MultistepBlocks
38、olverforImplicitDifferentialEquations)。 对求解器 PBIDE 和 MBIDE 进行测试，与现在使用广泛的求解器进行比较，比较了它们的运行精度和计算量。选用结构动力学的算例、多体系统动力学的算例和几个其它类型的算例对本文的两个求解器进行了测试，测试的结果表明这两个求解器是效率很高的。航天工程实际问题中包含大量结构和机构系统，如卫星包含天线和帆板等柔性结构，天线和帆板展开过程又是机构运动。研究结构的动力学行为，通常是对结构动力学方程进行研究；研究机构的运动，通常利用多体系统动力学的建模理论，建立多体系统动力学方程。空间离散的结构动力学方程为常微分方程的初值问
39、题，多体系统动力学方程为微分-代数方程的初值问题。本文的主要目的是试图研究这两类方程的数值求解。 采用 Runge-Kutta 方法求解结构动力学方程。针对增量形式的结构动力学方程，推导了 Runge-Kutta 方法的迭代形式，研究了降低计算量的两种办法，发现了一些传统格式与 Runge-Kutta 方法之间的关系。基于单自由度自由振动方程，推导了 Runge-Kutta 方法求解结构动力学方程的逼近算子，应用逼近算子的谱半径可以分析 Runge-Kutta 方法求解二阶常微分方程的稳定性。通过谱半径图分析出各种格式的数值阻尼的强弱，发现高精度的 L-稳定格式具有理想的数值阻尼，在精确求解低
40、频振动的同时很好的抑制高频震荡。通过算例研究，发现较小的物理阻尼在精确解上能快速的耗散高频震荡，对数值计算却影响很小。从而得出结论：在数值计算中通过增加物理阻尼抑制高频震荡是不可行，而采用有数值阻尼的算法才是可行的。 应用适合求解高指标微分.代数方程的求解器，求解高指标形式的多体系统动力学方程。发现采用高指标形式的多体系统动力学方程，使用隐式形式微分方程需要引入新的求解变量少，而且能保持原系统矩阵的稀疏特性。但是存在两个缺点：(1)加速度变量和拉格朗日乘子的求解精度低于位移变量、速度变量的求解精度；(2)对协调一致的初始条件很依赖。针对这两个缺点，提出了对加速度变量和拉格朗日乘子修正的办法。采
41、用修正能提高加速度变量和拉格朗日乘子的求解精度，提高整个计算过程中的收敛性。数值算例表明，采用修正能降低求解的计算量，特别对于容许误差很小的情况，计算量降低更显著。 基于微分方程精确解的 Taylor 展式，提出了一种新的构造块格式方法；并采用指数函数的 Pade 逼近作为块格式的稳定性函数，构造出 L-稳定的块格式。实际上，块格式也是一类 Runge-Kutta 方法，从而通过 Runge-Kutta 方法的精度条件，分析了 L-稳定块格式的精度。针对隐式形式微分方程，推导了块格式的并行迭代和局部误差估计方法，采用变步长实现了 6 阶精度 L-稳定格式，编制了求解器PBIDE(Paralle
42、lBlocksolverforImplicitDifferentialEquations)。 针对Runge-Kutta 方法在求解微分.代数方程中出现精度降低的情况，提出了多步块格式。分析了多步块格式稳定性和精度条件，证明其为刚性精确的，且级精度与格式精度相同。数值算例表明在求解指标-2 的微分-代数方程不会出现精度降低现象。构造出了 2 级 3 阶精度的多步块格式，此格式的稳定性函数为指数函数的 4 阶精度广义 Pade 量近，且是 L-稳定的；构造出了 3 级 4 阶精度的多步块格式，此格式的稳定性函数与 Runge-Kutta 的稳定性函数一样，为指数函数的有理分式逼近，也是 L-稳定
43、格式。推导了多步块格式的 Nordsieek 表达。基于 2 级 3 阶精度的多步块格式，对隐式形式微分方程，采用变步长进行了实现，编制了求解器MBIDE(MultistepBlocksolverforImplicitDifferentialEquations)。 对求解器 PBIDE 和 MBIDE 进行测试，与现在使用广泛的求解器进行比较，比较了它们的运行精度和计算量。选用结构动力学的算例、多体系统动力学的算例和几个其它类型的算例对本文的两个求解器进行了测试，测试的结果表明这两个求解器是效率很高的。航天工程实际问题中包含大量结构和机构系统，如卫星包含天线和帆板等柔性结构，天线和帆板展开过程
44、又是机构运动。研究结构的动力学行为，通常是对结构动力学方程进行研究；研究机构的运动，通常利用多体系统动力学的建模理论，建立多体系统动力学方程。空间离散的结构动力学方程为常微分方程的初值问题，多体系统动力学方程为微分-代数方程的初值问题。本文的主要目的是试图研究这两类方程的数值求解。 采用 Runge-Kutta 方法求解结构动力学方程。针对增量形式的结构动力学方程，推导了 Runge-Kutta 方法的迭代形式，研究了降低计算量的两种办法，发现了一些传统格式与 Runge-Kutta 方法之间的关系。基于单自由度自由振动方程，推导了 Runge-Kutta 方法求解结构动力学方程的逼近算子，应
45、用逼近算子的谱半径可以分析 Runge-Kutta 方法求解二阶常微分方程的稳定性。通过谱半径图分析出各种格式的数值阻尼的强弱，发现高精度的 L-稳定格式具有理想的数值阻尼，在精确求解低频振动的同时很好的抑制高频震荡。通过算例研究，发现较小的物理阻尼在精确解上能快速的耗散高频震荡，对数值计算却影响很小。从而得出结论：在数值计算中通过增加物理阻尼抑制高频震荡是不可行，而采用有数值阻尼的算法才是可行的。 应用适合求解高指标微分.代数方程的求解器，求解高指标形式的多体系统动力学方程。发现采用高指标形式的多体系统动力学方程，使用隐式形式微分方程需要引入新的求解变量少，而且能保持原系统矩阵的稀疏特性。但
46、是存在两个缺点：(1)加速度变量和拉格朗日乘子的求解精度低于位移变量、速度变量的求解精度；(2)对协调一致的初始条件很依赖。针对这两个缺点，提出了对加速度变量和拉格朗日乘子修正的办法。采用修正能提高加速度变量和拉格朗日乘子的求解精度，提高整个计算过程中的收敛性。数值算例表明，采用修正能降低求解的计算量，特别对于容许误差很小的情况，计算量降低更显著。 基于微分方程精确解的 Taylor 展式，提出了一种新的构造块格式方法；并采用指数函数的 Pade 逼近作为块格式的稳定性函数，构造出 L-稳定的块格式。实际上，块格式也是一类 Runge-Kutta 方法，从而通过 Runge-Kutta 方法的
47、精度条件，分析了 L-稳定块格式的精度。针对隐式形式微分方程，推导了块格式的并行迭代和局部误差估计方法，采用变步长实现了 6 阶精度 L-稳定格式，编制了求解器PBIDE(ParallelBlocksolverforImplicitDifferentialEquations)。 针对Runge-Kutta 方法在求解微分.代数方程中出现精度降低的情况，提出了多步块格式。分析了多步块格式稳定性和精度条件，证明其为刚性精确的，且级精度与格式精度相同。数值算例表明在求解指标-2 的微分-代数方程不会出现精度降低现象。构造出了 2 级 3 阶精度的多步块格式，此格式的稳定性函数为指数函数的 4 阶精度
48、广义 Pade 量近，且是 L-稳定的；构造出了 3 级 4 阶精度的多步块格式，此格式的稳定性函数与 Runge-Kutta 的稳定性函数一样，为指数函数的有理分式逼近，也是 L-稳定格式。推导了多步块格式的 Nordsieek 表达。基于 2 级 3 阶精度的多步块格式，对隐式形式微分方程，采用变步长进行了实现，编制了求解器MBIDE(MultistepBlocksolverforImplicitDifferentialEquations)。 对求解器 PBIDE 和 MBIDE 进行测试，与现在使用广泛的求解器进行比较，比较了它们的运行精度和计算量。选用结构动力学的算例、多体系统动力学的
49、算例和几个其它类型的算例对本文的两个求解器进行了测试，测试的结果表明这两个求解器是效率很高的。航天工程实际问题中包含大量结构和机构系统，如卫星包含天线和帆板等柔性结构，天线和帆板展开过程又是机构运动。研究结构的动力学行为，通常是对结构动力学方程进行研究；研究机构的运动，通常利用多体系统动力学的建模理论，建立多体系统动力学方程。空间离散的结构动力学方程为常微分方程的初值问题，多体系统动力学方程为微分-代数方程的初值问题。本文的主要目的是试图研究这两类方程的数值求解。 采用 Runge-Kutta 方法求解结构动力学方程。针对增量形式的结构动力学方程，推导了 Runge-Kutta 方法的迭代形式，研究了降低计算量的两种办法，发现了一些传统格式与 Runge-Kutta 方法之间的关系。基于单自由度自由振动方程，推导了 Runge-Kutta 方法求解结构动力学方程的逼近算子，应用逼近算子的谱半径可以分析 Runge-Kutta 方法求解二阶常微分方程的稳定性。通过谱半径图分析出各种格式的数值阻尼的强弱，发现高精度的 L-稳定格式具有理想的数值阻尼，在精确求解低频振动的同时很好的抑制高频震荡。通过算例研究，发现较小的物理阻尼在精确解上能快速的耗散高频震荡，对数值计算却影响很小。从而得出结论：在
