连续时间线性约束LQ控制问题的时程精细积分方法

本文基于连续时间线性约束LQ控制问题,给出时段的消元公式。由于消元过程与消元次序无关,故可在此过程中引入2N类高精度时程积分方法,求出Riccati方程后,对状态向量进一步采用时程精细积分法,可确定系统非常精确的状态向量,该方法不仅保证了系统的计算精度,而且有很好的数值稳定性,数值例题说明了本文方法的有效性。     

精细积分法应用于地震碰撞力反应谱计算研究

相邻结构的碰撞时程分析是计算地震碰撞力反应谱的基础。结构碰撞时程分析要求采用稳定性好、精度高及计算效率高的数值分析方法。精细积分法将二阶动力微分方程通过增元降阶的方式转换成Hamilton对偶变量体系,得到了动力微分方程的精确解。基于此,本文将精细积分法引入结构碰撞时程分析及地震碰撞力反应谱计算中。在推导精细积分法公式的基础上,在MATLAB环境下编制了结构碰撞时程分析程序和碰撞力反应谱计算程序,并实现了碰撞力反应谱程序的并行化。经算例验证,精细积分法应用于结构碰撞时程分析及地震碰撞力反应谱计算是可行的,程序计算结果准确。     

结构动力响应精细时程法的一种并行算法

结构动力响应精细时程法的并行算法分为两类：基于特解的并行算法和基于直接积分法的并行算法;后者因为不需知道荷载的具体形式而更具应用价值。精细时程法的时程积分由齐次方程的通解和非齐次项的积分构成,基于直接积分法的并行算法很好地并行了非齐次项的积分,而对通解项采用串行计算。设计了一种不均衡步数的负载分配策略,能够减少处理器等待自身初值的时间,相对均衡步数的分配策略,能够获得更高的加速比,给出了相应的证明和算例验证。     

子域精细积分及偏微分方程数值解

对于偏微分方程半解析法的方程，精细时程积分虽然能求出高度准确的解，但往往面临矩阵尺度太大的困难，另一方面差分法虽然有带宽小的优点，但有稳定性及精度方面的问题，本文提出子域精细积分法，既可利用精细积分的数值优点，又有带宽小的好处，数值例题表明了子域精细积分法的效能。     

子域精细积分及偏微分方程数值解

对于偏微分方程半解析法的方程，精细时程积分虽然能求出高度准确的解，但往往面临矩阵尺度太大的困难；另一方面差分法虽然有带宽小的优点，但有稳定性及精度方面的问题．本文提出子域精细积分法，既可利用精细积分的数值优点，又有带宽小的好处．数值例题表明了子域精细积分法的效能．     

非线性度对修正双步长显式法及常用逐步积分法的影响

为了掌握非线性度对逐步积分法的影响,研究了几种积分算法在不同非线性度振动系统中的响应。通过3个典型非线性算例,对修正双步长显式法、蛙跳式中心差分法、Newmark法、广义α法和精细积分法的计算精度和稳定性能等进行了比较。结果表明:非线性度对广义α法、精细积分法和Newmark法的稳定性有影响;高非线性度对Newmark法的计算稳定性影响最大;时间步长越小,算法精度和计算量越高;相同小步长情况下,精细积分法的精度最高,而修正双步长显式法的计算量最小;在时间步长较大时,低非线性度会引起精细积分法不稳定,修正双步长显式法的精度最高,修正双步长显式法在非线性系统中具有很强的鲁棒性。     

结构动态载荷识别的精细逐步积分法

对比例阻尼系统给出了基于精细逐步积分法的动态载荷识别方法。首先将系统进行模态坐标变换得到无耦合运动方程 ,然后应用精细逐步积分法构造一种高效精确的载荷识别公式 ,再由量测到的结构动态响应求出动态力的时间历程。数值算例验证了本方法的识别精度是好的     

增维精细积分法的适用范围研究

对线性定常结构动力系统提出的增维精细积分法,能够将非齐次动力方程转化为齐次动力方程,不用对状态矩阵求逆就能方便高效地求解出结构的动力响应。本文在仔细分析增维精细积分法性质的基础上,提出了其适用条件,进一步拓宽了其应用范围,并给出了将荷载项展开成傅里叶级数时,相应增维精细积分法的表达式。同时,在一个时间步长内,通过对非齐次项作线性化假设,成功地将增维精细积分法应用到了非线性动力分析领域。本文方法计算格式统一,易于编程,具有很高的计算效率。数值算例证明了本文方法的有效性。     

二阶双曲型方程的精细时程积分法

对于二阶双曲型偏微分方程初边值问题，可以用有限差分法进行求解。通常的有限差分法在使用过程中受到精确度和稳定性的限制，本文提出求解二阶双曲型方程的精细时程积分法。由于这种方法是半解析方法，在时间域上可以精确计算，所以这种方法不仅精确度高，而且还绝对稳定。文末的数值算例进一步验证了上述结构，而且对大的时间步长（例如△t=0.5）仍然获得精度很高的数值结果。可见，精细时程积分法是一种很实用的方法。     

非参数概率系统非平稳随机响应研究

对复杂系统的测量和知识的有限性造成的不确定很难完全由随机参数模型进行描述,采用随机矩阵模型更具有一般性和合理性.本文应用随机矩阵模拟不确定线性动力系统有限元模型中质量阵、阻尼阵和刚度阵的概率不确定性.综合运用虚拟激励法和精细时程积分法建立了非参数概率系统非平稳随机响应的有效计算方法.数值模拟结果表明,对于高精度制造,模型的不确定性是不能忽略的.本文提出的算法为此类问题求解提供了一条有效途径.     

结构动力方程的增维精细积分法

对线性定常结构动力系统提出的精细积分方法,能够得到在数值上逼近于精确解的结果,但对于非齐次动力方程涉及到矩阵求逆的困难。提出采用增维的办法,将非齐次动力方程转化为齐次动力方程,在实施精细积分过程中不必进行矩阵求逆,这种方法对于程序实现和提高数值稳定性十分有利,而且在大型问题中计算效率较高,从而改进了精细积分方法的应用,数值例题显示了本文方法的有效性。     

结构动力方程的样条精细积分法

结合精细积分法和样条函数拟合技术的优点,提出了求解结构动力方程的一种有效方法.首先对非齐次项用三次正规化B样条函数进行拟合,然后利用正规化B样条函数形状相同、仅相差一个平移量的特点,构造了一个高效的特解求解方法.按此方法只需求出一个标准B样条项所对应的特解,然后通过时间坐标的平移并结合叠加原理,即可求出任意时刻的特解值.由于特解计算中采用数值积分的方法,避免了矩阵求逆,因而本方法具有较大的适用范围.算例结果证明了该方法的有效性.     

波动方程的精细逐步积分法

应用现有的波动方程求解方法解决工程实际问题尚存在一定的局限性。本文在结构动力方程精细逐步积分的基础上,提出了波动方程初边值问题的精细逐步积分法,并分别给出了不同边界条件下的精细逐步积分格式。此数值方法虽然是显式积分方法,却是无条件稳定的。分别用精细逐步积分法和其它已有的方法对两个算例进行了计算,一个是有解析解的例子,该例验证了此方法的准确性,另一个例子是求解由波动方程及初始条件和边界条件组成的有杆抽油系统预测模型,此例验证了精细逐步积分法的高效性。     

饱和两相介质动力固结问题时域求解的精细时程积分方法

针对u-p形式的饱和两相介质波动方程,采用精细时程积分方法计算固相位移u,采用向后差分算法求解流体压力p,建立了饱和两相介质动力固结问题时域求解的精细时程积分方法。针对标准算例,对该方法的计算精度进行了校核。开展了该方法相关算法特性的研究,对采用不同数值积分方法计算非齐次波动方程特解项计算精度的差异进行了对比研究,并对采用不同积分点数目的高斯积分法计算特解项条件下计算精度的差异进行了对比研究。研究结果表明,(1)该方法具有良好的计算精度。(2)计算非齐次波动方程特解项的数值积分方法中,梯形积分法的计算精度最差,高斯积分法、辛普生积分法和科茨积分法都具有较好的计算精度。(3)增加高斯积分点数目对于提高计算精度的作用并不显著。     

瞬态热传导方程的子结构精细积分方法

对线性定常结构动力系统提出的精细积分方法,在数值精度等方面表现出极大优越性,但是当矩阵尺度很大时在数值计算与存储中将产生困难,对此,本文对瞬态热传导方程,根据结构的概念,将结构分为若干个子结构,对各子结构分别进行指数矩阵运算并通过了结构间界面的物理量相联系,从而提高精细积分方法的计算效率。     

Identification of discrete-time bilinear systems through equivalent linear models

This paper presents a novel method for identification of discrete-time, time-invariant state-space models of bilinear dynamical systems using the steady-state portion of a single input/multiple output time-history measurements. These measurements are recorded by exciting the system with a linear combination of sine and cosine functions of user-selected frequencies enriched by a subtle amount of random component. The proposed method relies on conversion of the bilinear system into an equivalent linear model (ELM) by an accurate approximation of the state in the bilinear term using a set of sine and cosine basis functions whose frequencies are obtained as combinations of the input frequencies. Observer/Kalman Filter Identification (OKID), a?linear time invariant (LTI) system identification algorithm, is used to identify the aforementioned ELM from which the original bilinear model is recovered. A?numerical example is also provided.     

精细时程积分法的误差分析与精度设计

通过对精细积分法递推过程的误差分析，发现该方法能莸得高精度数值结果的根本原因是：数值计算的相对误差不随递推过程的进行而扩散。数值结果的精度仅仅取决于初始Taylor级数的计算精度和指数矩阵A的最大模特征。同时，提出了一种精度估计和精度设计的方法。     

非粘滞阻尼系统时程响应分析的精细积分方法

考虑一个具有非粘滞阻尼特性的多自由度系统响应的时程分析问题.该非粘滞阻尼模型假设阻尼力与质点速度的时间历程相关,数学表达式体现为阻尼力等于质点速度与某一核函数的卷积.在利用状态空间方法将系统运动方程转换成一阶的状态方程的基础上,采用精细积分方法对状态方程进行数值求解,得到一种求解该阻尼系统时程响应的精确、高效的计算方法.通过两个数值算例表明,采用该方法得到几乎精确的数值计算结果,而且计算效率有成数量级的提高.     

基于二阶系统解耦的精细积分格式下动载荷时域识别

基于模态分析法的载荷识别方法利用模态矩阵获得系统的非耦合形式,推导单自由度系统的载荷识别公式,但要求系统为比例阻尼。在模态模型基础上,对一般系统建立基于二阶系统解耦的动载荷时域识别模型。首先,利用基于Lancaster结构的二阶系统解耦方法推导出系统的非耦合形式;然后,采用精细逐步积分方法,在载荷为阶跃力的假设下,推导出载荷识别数学公式;最后,由系统实时响应反求结构载荷的时间历程。数值算例验证了本文方法针对比例阻尼系统较模态模型具有更高的精度,而且对非比例阻尼系统也有效可行。