首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 93 毫秒
1.
提出了哈密顿动力系统的一个新变分原理,并基于此变分原理构造了四类保辛算法。通过新的变分原理定义修正作用量,然后将位移和动量采用拉格朗日多项式近似,并采用高斯积分对时间近似积分得到近似的修正作用量。在修正作用量的基础上,通过选择时间步两端不同的位移或动量作为独立变量,可构造四种不同类型的保辛算法。  相似文献   

2.
给出了非传统哈密顿变分原理的一种简化形式,并在此基础上利用拉格朗日多项式近似位移和动量,采用高斯积分法对时间积分,建立了针对动力学初值问题的一类高阶辛算法。在建立高阶辛算法的过程中,本文方法与基于传统哈密顿变分原理的辛算法不同,无需由端值问题向初值问题转换,因此更加简捷有效。此外,给出了线性动力问题中本文算法保辛性的证明。当位移、动量的插值次数和高斯积分点个数均为m时,本文算法是具有2m阶精度的辛算法,且是线性无条件稳定的。通过数值算例结果表明,本文算法与辛算法性质吻合,并且计算效率比同阶辛龙格库塔法提高了约50%。  相似文献   

3.
文献[1,2]给出了四种不同类型的求解哈密顿动力系统的数值方法,并证明了它们的保辛特性。本文将讨论这四类算法的具体数值性能,包括算法的线性稳定性,精度和效率等。  相似文献   

4.
在非增量算法的基础上,提出了用最优控制变分原理形成过程最优控制迭代求解的基本思路,并给出求解的基本控制方程。这一工作为有限变形力学问题的数值求解提供了一个新的处理方法。  相似文献   

5.
基于对偶变量变分原理,选择积分区间两端位移为独立变量,构造了求解完整约束哈密顿动力系统的高阶保辛算法。首先,利用拉格朗日多项式对作用量中的位移、动量及拉格朗日乘子进行近似;然后,对作用量中不包含约束的积分项采用Gauss积分近似,对作用量中包含约束的积分项采用Lobatto积分近似,从而得到近似作用量;最后,在此近似作用量的基础上,利用对偶变量变分原理,将求解完整约束哈密顿动力系统问题转化为一组非线性方程组的求解。算法具有保辛性和高阶收敛性,能够在位移的插值点处高精度地满足完整约束。算法的收敛阶数及数值性质通过数值算例验证。  相似文献   

6.
非线性轨迹优化问题的保辛自适应求解方法   总被引:1,自引:0,他引:1  
非线性轨迹优化问题一般是一个非线性最优控制问题。将非线性系统的最优控制问题导入到哈密顿体系的辛几何空间当中,基于对偶变量变分原理提出了求解非线性最优控制问题的一种保辛自适应方法。以时间区段两端协态作为独立变量,在时间区段内采用拉格朗日插值近似状态和协态变量,并利用对偶变量变分原理将非线性最优控制问题转化为非线性方程组的求解,保持了哈密顿系统的辛几何结构。并进一步,提出了基于多层次迭代的自适应算法,提高了非线性最优控制问题的求解效率。数值实验验证了该算法在求解非线性轨迹优化问题中的有效性。  相似文献   

7.
本文讨论了第二类抛物型变分不等式中的MRM(多重互易方法)方法。首先采用时间项半离散和隐格式方法将抛物型变分不等式化解为一个椭圆变分不等式,然后利用MRM-边界积分方程,将其化解为MRM-边界混合变分不等式,并给出了MRM-边界混合变分不等式解的存在唯一性。说明了该MRM-边界混合变分不等式与常规边界积分方程得到的边界混合变分不等式是一致的,并且具有更容易编程实现。这为使用MRM边界元方法数值求解抛物型变分不等式提供了方法和理论依据。文末给出了数值算例。  相似文献   

8.
关于线粘弹性动力学中各种变分原理   总被引:7,自引:0,他引:7  
罗恩 《力学学报》1990,22(4):484-489
本文提出一条简单而统一的新途径,系统地建立了线粘弹性动力学中各种简化Gurtin型变分原理,文中首先给出一个很有用的以卷积表示的积分关系式,然后从该式出发,系统地导出成互补关系的五类变量、四类变量、三类变量、二类变量及一类变量简化Gurtin型变分原理,并清楚地阐明它们之间的内在联系,而且,还发现当前在国际上有广泛影响的力学变分原理方面的名著[1]及文[2]中,所给出的四个变分原理的泛函式均有误.本文除给出这四个正确的泛函式外,还建立了一些新的更一般广义变分原理。  相似文献   

9.
用变分法估计热爆炸临界状态参数   总被引:4,自引:4,他引:0  
本文利用作者给出的关于爆炸临界状态的变分原理,并用变分的直接方法,得到了六种不同形状(平板、无限长圆柱,球,无限长方形杆,有限长圆柱,正方体)系统的临界P-K参数和临界温度的近似公式,与精确的数值计算结果比较表明,对点火临界参数,本文公式均能给出较好的估计值。  相似文献   

10.
不可压流体饱和多孔弹性梁的变分原理及有限元方法   总被引:3,自引:1,他引:2  
基于不可压饱和多孔弹性梁动力弯曲的数学模型,建立了以多孔弹性梁挠度和孔隙流体压力等效力偶为宗量的Gurtin型变分原理,并给出了特殊边界条件下解耦时的仅以挠度为宗量的变分原理.同时,作为动力响应的退化情形,讨论了拟静态情形下的相应变分原理.根据所建立的变分原理,导出了一个有限元离散公式.由于Gurtin型变分原理是关于时间的卷积型的泛函,空间的有限元离散导致一个关于时间的对称微分一积分方程组,此方程组可进一步转化为常微分方程组.利用隐式Euler法,给出了时间区域的计算格式.作为一个数值例子,分析了饱和多孔弹性悬臂梁在自由端简谐载荷作用下的动力响应,分析了流相与固相相互作用对饱和多孔弹性悬臂梁动力响应的影响.  相似文献   

11.
用有限元广义混合法分析不可压缩或几乎不可压缩弹性体   总被引:2,自引:0,他引:2  
不可压缩或几乎不可压缩问题在数学上表现为最小 势能原理中的某些项趋于无穷大,使得有限元方程产生病态。本文给出了不可压缩或几乎不可压缩弹性分析的广义混合变分原理,以此为基础建立了该类问题的有限元广义混合法。该变分原理的泛函中不含有上面这种奇异项,故其有限元方程不会产生病态。算例表明该有限元法可以同时进行可压缩、不可压缩或几乎不可压缩弹性分析,且精度良好;有限元常规位移法及Hermann法是该法的特例。  相似文献   

12.
参变量变分原理及其参数二次规划算法是由钟万勰院士1985年针对弹性接触边界非线性问题首次提出来的,经过将近40年的不断发展,目前参变量变分原理已经成功应用于各个领域,其中包括弹塑性分析、接触问题、润滑力学、岩土力学、变刚度杆系结构、先进材料性能分析、材料的蠕变与损伤、柔性结构力学和LQ最优控制等各个工程领域。本文首先回顾了参变量变分原理的起源,介绍了参变量变分原理的基本概念,然后以弹塑性分析问题为例,阐明建立参变量变分原理的理论模型以及实现数值参数二次规划求解原理,最后详细回顾了参变量变分原理的基本理论与相应数值算法在各个领域的发展及其工程应用,展示了参变量变分原理在求解各类非线性问题的特色与优势。  相似文献   

13.
Gurtin变分原理在矩形板动力初值问题中的应用   总被引:3,自引:0,他引:3  
结构动力分析是工程设计中的重要组成部分,传统动力分析方法并不能全面反映动力初值特征,而Gurtin变分原理则被认为是目前唯一能全面反映动力初值特征的变分原理。本文基于位移型Gurtin变分原理,对空间和时间同时离散,建立了一种求解板的动力初值问题的时空有限元法,并对两种边界情况板的振动问题进行了编程计算,计算结果表明时空元法精确度很高且稳定收敛。  相似文献   

14.
We consider a problem on shock wave localization in the numerical solution of one-dimensional unsteady problems of gas dynamics in Eulerian variables obtained on the basis of finite difference shock-capturing schemes. An optimization method for strong discontinuity localization proposed previously by Miranker and Pironneau is investigated by means of methods of classical variational calculus. This method may be difficult to implement when the entropy condition is included in the formulation of Miranker and Pironneau's optimization problem as an active constraint. In this connection we suggest an alternative optimization problem using artificial viscosity in the variational principle. It is shown theoretically that the application of such a variational principle yields a trajectory which coincides with the true discontinuity trajectory in the case of a shock wave moving at a constant speed. On the basis of this modification one more algorithm is proposed which reduces the shock localization problem to a problem of minimization of a univariate function. Numerical tests corroborate completely the theoretical conclusions. In particular, a higher shock localization accuracy is obtained on the basis of the proposed algorithms as compared to the original Miranker-Pironneau method.  相似文献   

15.
近些年弹性力学中出现一种新型的变分原理,称为广义混合变分原理.特点是其泛函中包含某些可以任意选择的附加函数,称为分裂因子.新原理将弹性理论中现有的各主要变分原理都统一在一个框架中,并揭示出它们之间更深一层的相互关系.在应用方面,它提供了一个新的数学手段以建立有限元分析中的新模式.这些新模式已经显示出它们的优点:适当调节分裂因子,它们给出更好的数值解答,特别是,可用它们来处理有限元方法中棘手的病态问题.本文综述了线性及非线性弹性理论中的这种新型变分原理并就其在有限元中的应用作了讨论.  相似文献   

16.
广义协调元的变分基础及几何不变性   总被引:6,自引:0,他引:6  
本文根据修正的势能原理推出了广义协调元列式。从而为广义协调元提供了一种变分依据。同时讨论了广义协调元的几何不变性,即单元与其节点号编序的相关性问题。在此基础上推导了一个任意四边形薄板弯曲单元,数值结果表明该单元能保证收敛,具有与单元节点号编序无关性,与现有的同类单元相比,具有较高的精度。  相似文献   

17.
孙雁  高强  钟万勰 《计算力学学报》2015,32(5):595-600,607
首先将非线性Schrdinger方程化为Hamilton正则方程形式,而后建立Hamilton体系下的变分原理。再用有限元法离散空间坐标,同时对时间坐标进行精细积分,最后运用混合能变分原理,提出非线性Schrdinger方程保辛数值解法。这种解法在保辛的同时,可以让能量和质量在积分格点上亦全部达到守恒。数值算例验证了该方法的有效性。  相似文献   

18.
将参变量变分原理引入铁电问题。对一类借用了经典弹塑性理论中的概念和方法的多轴铁电模型建立基于Helmholtz自由能的参变量变分原理,可以有效处理传统变分原理中由非关联流动法则或屈服面不考虑材料系数变化所引起的切线模量非对称困难。相应于参变量变分原理,引入参数二次规划算法,可获得具有可靠数值稳定性的一套铁电算法。将该算法应用于一个具体的铁电模型,数值计算结果表明本文方法的有效性。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号