共查询到18条相似文献,搜索用时 187 毫秒
1.
非线性动力方程的增维精细积分法 总被引:30,自引:0,他引:30
对线性定常结构的动力系统提出的精细积分法,能得到在数值上逼近于精确解的结果。但是对于非齐次动力方程却涉及到矩阵求逆的困难,而且通常与时间有关的非齐次项不能进入精细积分的细化过程。采用增维的方法,将非齐次动力方程化为齐次方程,在实施精细积分的过程中不必进行矩阵求逆。这种处理方法对于程序实现和提高数值计算的稳定性十分有利,而且在大型问题中可明显提高计算效率,数值算例显示本文方法是有效的。 相似文献
2.
3.
提出一种针对非线性动力方程的改进精细积分方法。该方法是在时间步长内采用分段的三次样条函数拟合非齐次项,保持高精度拟合的同时避免了求导运算和高次多项式插值带来的Runge现象。通过引入4×2个变量将动力方程增加四维转化为齐次方程,并建立相应的通解格式,避免了状态空间下系统矩阵求逆。将指数矩阵分为四个子模块,利用各模块的特点分别进行理论推导及基于精细积分法进行分步、分块计算得到相应的理论解和高精度数值解,无需反复计算整个指数矩阵,提高了解算效率。针对含未知状态量的非齐次项,引入预测-校正的方法进行迭代求解。数值计算结果表明了本文方法的有效性。 相似文献
4.
5.
6.
7.
针对u-p形式的饱和两相介质波动方程,采用精细时程积分方法计算固相位移u,采用向后差分算法求解流体压力p,建立了饱和两相介质动力固结问题时域求解的精细时程积分方法。针对标准算例,对该方法的计算精度进行了校核。开展了该方法相关算法特性的研究,对采用不同数值积分方法计算非齐次波动方程特解项计算精度的差异进行了对比研究,并对采用不同积分点数目的高斯积分法计算特解项条件下计算精度的差异进行了对比研究。研究结果表明,(1)该方法具有良好的计算精度。(2)计算非齐次波动方程特解项的数值积分方法中,梯形积分法的计算精度最差,高斯积分法、辛普生积分法和科茨积分法都具有较好的计算精度。(3)增加高斯积分点数目对于提高计算精度的作用并不显著。 相似文献
8.
9.
采用双重互易边界元法结合精细积分法求解二维含热源的瞬态热传导问题。针对边界积分方程中热源项和温度关于时间导数项引起的域积分,采用双重互易法处理,将域积分转换为边界积分。采用边界元法将边界积分方程离散后,得到关于时间的微分方程组,并利用精细积分法处理其中的指数型矩阵;对于微分方程组中由边界条件和热源项引起的非齐次项,采用解析的方法计算。为了比较精细积分-双重互易边界元法的计算效果,同时使用有限差分法计算温度对时间的导数项。通过数值算例验证了本文方法的有效性和精确性。计算结果表明:时间步长对于精细积分-双重互易边界元法的结果影响较小,而有限差分法对时间步长比较敏感且只在时间步长选取较小时有效;当选取较大时间步长时,精细积分-双重互易边界元法依然具有良好的计算精度。 相似文献
10.
常规位移有限元的结构振动方程是n个二阶常微分方程组.采用一般交分原理推导,将结构振动问题引入Hamiltoil体系,将得到2n个一阶常微分方程组.精细积分法宜于处理一阶方程,应用于线性定常结构动力问题求解,可以得到在数值上逼近精确解的结果.对于非齐次动力方程,当结构具有刚体位移时,系统矩阵将出现奇异.本文借鉴全元选大元高斯-约当法求解线性方程组的经验,提出全元选大元法求奇异矩阵零本征解的方法,该方法可以简便快速地寻求奇异矩阵零本征值对应的子空间.利用Hamiltoil体系已有研究成果及Hamilton系统的共轭辛正交归一关系,迅速将零本征值对应的子空间分离出来,通过投影排除奇异部分,然后用精细积分法求得问题的解.数值算例表明,该方法对Hamilton系统奇异问题,处理方便,计算量小,易于实现,同时保持了精细算法的优点. 相似文献
11.
12.
悬索桥颤振稳定性分析的精细时程积分法 总被引:3,自引:0,他引:3
研究精细时程积分法在悬索桥颤振稳定性分析中的应用,首先,将 气流整个作为一个系统,组集系统关于模态广义坐标的状态空间方程,然后,应用精细时程积分法计算状态的向量的时程响应,根据状态向量时程响应的对数衰减率判断系统的颤振稳定性,最后,以英国塞文悬索桥为数值算例,验证了本文方法的正确性。 相似文献
13.
A PREDICT-CORRECT NUMERICAL INTEGRATION SCHEME FOR SOLVING NONLINEAR DYNAMIC EQUATIONS 总被引:1,自引:0,他引:1
A new numerical integration scheme incorporating a predict-correct algorithm forsolving the nonlinear dynamic systems was proposed in this paper. A nonlinear dynamic systemgoverned by the equation v=F(v,t) was transformed into the form as v=Hv f(v,t). Thenonlinear part f(v,t) was then expanded by Taylor series and only the first-order term retained inthe polynomial. Utilizing the theory of linear differential equation and the precise time-integrationmethod, an exact solution for linearizing equation was obtained. In order to find the solution of theoriginal system, a third-order interpolation polynomial of v was used and an equivalent nonlinearordinary differential equation was regenerated. With a predicted solution as an initial value andan iteration scheme, a corrected result was achieved. Since the error caused by linearization couldbe eliminated in the correction process, the accuracy of calculation was improved greatly. Threeengineering scenarios were used to assess the accuracy and reliability of the proposed method andthe results were satisfactory. 相似文献
14.
基于模态分析法的载荷识别方法利用模态矩阵获得系统的非耦合形式,推导单自由度系统的载荷识别公式,但要求系统为比例阻尼。在模态模型基础上,对一般系统建立基于二阶系统解耦的动载荷时域识别模型。首先,利用基于Lancaster结构的二阶系统解耦方法推导出系统的非耦合形式;然后,采用精细逐步积分方法,在载荷为阶跃力的假设下,推导出载荷识别数学公式;最后,由系统实时响应反求结构载荷的时间历程。数值算例验证了本文方法针对比例阻尼系统较模态模型具有更高的精度,而且对非比例阻尼系统也有效可行。 相似文献
15.
16.
对线性定常结构动力系统提出的精细积分方法,在数值精度等方面表现出极大优越性,但是当矩阵尺度很大时在数值计算与存储中将产生困难,对此,本文对瞬态热传导方程,根据结构的概念,将结构分为若干个子结构,对各子结构分别进行指数矩阵运算并通过了结构间界面的物理量相联系,从而提高精细积分方法的计算效率。 相似文献
17.
非线性最优控制系统的时程精细计算研究 总被引:2,自引:1,他引:2
针对非线性最优控制问题 ,通过一阶 Taylor级数展开 ,得到线性化的动力学方程 ,进而在方程原变量的基础上 ,引入对偶向量 (Lagrange乘子向量 ) ,将动力学方程从 Lagrange体系引入到了 Hamilton体系 ,在全状态下 ,从一个新的角度对非线性最优控制问题进行了描述 ,进一步基于时程精细积分理论 ,对其方程进行了有效的精细求解 ,并通过算例说明了文中方法的有效性 相似文献
18.
相邻结构的碰撞时程分析是计算地震碰撞力反应谱的基础。结构碰撞时程分析要求采用稳定性好、精度高及计算效率高的数值分析方法。精细积分法将二阶动力微分方程通过增元降阶的方式转换成Hamilton对偶变量体系,得到了动力微分方程的精确解。基于此,本文将精细积分法引入结构碰撞时程分析及地震碰撞力反应谱计算中。在推导精细积分法公式的基础上,在MATLAB环境下编制了结构碰撞时程分析程序和碰撞力反应谱计算程序,并实现了碰撞力反应谱程序的并行化。经算例验证,精细积分法应用于结构碰撞时程分析及地震碰撞力反应谱计算是可行的,程序计算结果准确。 相似文献