基于对角隐式Runge-Kutta公式的无约束优化方法

通过把一个无约束优化问题转化为一个等价的常微分方程,利用二阶半对角隐式Runge Kutta公式构造了求解无约束优化问题的LRKOPT算法。LRKOPT算法具有与IMPBOT方法相似的数值特性,但LRKOPT算法可以看成是最速下降方向与牛顿法方向的非线性组合,而IMPBOT方法为它们两者之间的线性组合。在目标函数为一致凸函数的假设条件下,证明了LRKOPT方法的具有全局收敛和局部超线性收敛性。数值结果表明LRKOPT方法具有很好的数值稳定性并且LRKOPT方法的计算效率优于IMPBOT方法。

四级四阶对角隐式辛Runge-Kutta方法参数计算

1.引 言设有Hamilton系统这坐 H（p1,…,pn,q1,…,qn）是 Hamilton函数,它和t无关.记z=（p1,…;Pn,q1,…,qn）T和Hamilton系统（1）的右端项为f（z）,则（1）可表示为dz/dt=f（z）. 冯康用辛几何的观点提出了计算Hamilton系统的辛差分格式[1].Runge-Kutta方法是求非线性常微分方程（组）数值解的重要单步方法.若能找到具有辛性的Runge-Kutta方法,对于求解非线性 Hamilton系统数值解将具有非常重要的意义.J.M.

一些具有较少迭代次数和较好稳定性的新型 PDIRK方法(英文)

为求解常微分方程刚性问题 ,本文构造了一些新型的并行对角隐式迭代 Runge-Kutta方法 .与已有的一些方法相比 ,这些方法具有更少的迭代次数和更好的稳定性 (强 A-稳定或 L-稳定 )且更适于并行计算 .

钽的层裂实验数值模拟

采用连续介质力学基唯象模型模拟分析了钽的平板撞击层裂行为。该模型包括了材料的非线性弹性(状态方程)、率相关塑性和孔洞的形核及生长等多种效应,并且采用一种对角隐式Runge-Kutta方法来求解本构率方程组,提高了热粘塑性本构关系计算的稳定性及精度。将数值模拟结果和相关实验数据进行了对比分析,结果表明,对于样品中的拉应力峰值明显高于材料层裂强度的实验(中、高速平板撞击实验),理论模型具有较好的预估能力,但对于临界层裂问题(低速平板撞击实验),该模型对材料损伤与失效过程的描述可能不够准确,需要进一步改进。