时滞Lienard非线性系统的Hopf分岔

本文研究了Ｌｉｅｎａｒｄ非线性时滞系统的线性稳定性和Ｈｏｐｆ分岔，考究了特征方程随两参数变化时根的分布，应用中心流形和范式分析失去线性稳定性出现的Ｈｏｐｆ分岔及其稳定性。     

基于有限变形理论的数值流形方法研究简

原有数值流形方法通过积累每一时步的小变形而得到结构最终的大变形,然而,当结构发生大变形、大转动时往往产生较大计算误差. 针对该问题,从动量守恒方程以及应力边界条件的积分弱形式出发,引入流形方法的插值函数,建立了基于有限变形理论的数值流形方法. 通过对比改进前后流形方法的计算迭代格式,指出了原有流形方法计算大变形问题时的误差来源. 最后,通过大变形悬臂梁和旋转块体算例对有限变形流形方法进行了验证. 数值结果表明,改进后的流形方法能够很好地处理大变形大转动问题,消除了转动所带来的计算误差,其计算结果与解析解及ABAQUS 软件求得的数值解相吻合.     

基于有限变形理论的数值流形方法研究

原有数值流形方法通过积累每一时步的小变形而得到结构最终的大变形,然而,当结构发生大变形、大转动时往往产生较大计算误差. 针对该问题,从动量守恒方程以及应力边界条件的积分弱形式出发,引入流形方法的插值函数,建立了基于有限变形理论的数值流形方法. 通过对比改进前后流形方法的计算迭代格式,指出了原有流形方法计算大变形问题时的误差来源. 最后,通过大变形悬臂梁和旋转块体算例对有限变形流形方法进行了验证. 数值结果表明,改进后的流形方法能够很好地处理大变形大转动问题,消除了转动所带来的计算误差,其计算结果与解析解及ABAQUS 软件求得的数值解相吻合.      

二阶流形元法与结构变形分析

在原有一阶流形元法的基础上开发了二阶流形元法数值仿真模型和相应的计算程序，并给出了计算实例。结果表明，二阶流形元法可以以较高的精度分析一般结构的变形和接触应力问题，对大变形问题独有优势。并能很好地模拟不连续介质的破坏过程及块体破坏后的运动。     

基于单位分解法的无网格数值流形方法

在数值流形方法和单位分解法的基础上，提出了无网格数值流形方法. 无网格数值流形 方法在分析时采用了双重覆盖系统，即数学覆盖和物理覆盖. 数学覆盖提供的节点形成求解 域的有限覆盖和单位分解函数；而物理覆盖描述问题的几何区域及其域内不连续性. 与原有 的数值流形方法相比，无网格数值流形方法的数学覆盖形状更加灵活，可以用一系列节点的 影响域来建立数学覆盖和单位分解函数，具有无网格方法的特性，从而摆脱了传统的数值流 形方法中网格所带来的困难. 与无网格方法相比，由于采用了有限覆盖技术，试函数的构造 不受域内不连续的影响，克服了原有的无网格方法在处理不连续问题时所遇到的困难. 详细推导了无网格数值流形方法的试函数和求解方程，最后给出了算例，验证了该方法的正 确性.     

弹性力学的复变量数值流形方法

数值流形方法通过引入数学和物理双重网格,将插值域和积分域分别定义在两个不同的覆盖上来完成系统能量泛函积分运算. 当采用高阶函数构造位移函数时,广义节点自由度将大大增加. 在求解系统的平衡方程中,运算量是与自由度的三次方成正比的,因此数值流形方法的计算量是较大的. 为此,在复变量理论的基础上,采用一维基函数建立二维问题的逼近试函数,然后将其应用于弹性力学的数值流形方法,提出了复变量数值流形方法,推导了弹性力学的复变量数值流形方法的公式. 与传统的数值流形方法相比,复变量数值流形方法具有计算量小、精度高的优点.      

数值流形方法中线性相关性问题的研究

数值流形方法的形函数由覆盖函数和局部近似函数组成,形函数之间往往存在线性相关性。在现有研究成果的基础上对形函数线性相关性进行了分析,指出线性相关性的根源在于覆盖函数具有单位分解特性,并与单元形状有关。研究了线性相关性与整体刚度矩阵奇异性以及求解收敛性之间的关系,指出形函数线性相关不一定导致整体刚度矩阵奇异。对8结点六面体高阶流形单元的局部近似函数及单元形状与线性相关性之间的关系进行了分析,构造出一种完全线性独立的流形单元。通过算例分析了8结点六面体流形单元局部近似函数中一次完全多项式对求解精度和收敛性的影响,发现采用一次完全多项式局部近似函数的形函数虽然线性相关,但求解仍然收敛,且精度高于线性无关的单元。     

三节点二次协调元

基于数值流形方法的思想，提出三节点二次协调单元，有效地改善了三节点线性单元的计算精度，数值算例表明该方法非常有效。     

弹性力学中的一种非协调数值流形方法

通过引入数学和物理双重网格，将插值域与 积分域分别定义在不同的覆盖上，即在数学网格上进行插值函数的构造，物理网格上完成 系统能量泛函积分运算，最后通过覆盖权函数将二者联结在一起. 它的优点是单元网格划 分随意，不受复杂边界形状和二相材料界面的限制，单元可以是任意形状，是较之于有限 元方法更一般的数值模拟方法. 在4节点四边形数值流形方法中，由于单元总体位移函数 包含的完全多项式不完全，使得计算精度不够精确，为此，在单元总体位移函数上附 加非协调位移基本项，使之趋于完全，提出了弹性力学问题的一种改进的数值流形 方法------非协调数值流形方法. 通过内部自由度静力凝聚处理，导出了消除内参后的单元应变矩阵 和单元刚度矩阵，使得在不增加广义节点自由度的前提下，大大提高了数值流形方法的计 算精度和计算效率. 同时对非协调项进行了显式处理，可以对工程实践起到更切实的帮助. 数值试验表明，它们能够保证收敛，有较高的精度，对畸变不敏感，从而证明了该方法的 可行性.     

粘弹性桩的混沌运动分析

研究了在轴向载荷和周期性横向载荷共同作用下非线性粘弹性嵌岩桩的混沌运动情况。假定桩和土体分别满足Leaderman非线性粘弹性和线性粘弹性本构关系，得到的运动方程为非线性偏微分．积分方程；利用Galerkin方法将方程简化为非线性常微分一积分方程，同时利用非线性动力系统中的数值方法，进行了数值计算，得到了不同载荷参数、几何参数、材料参数时粘弹性桩发生周期运动、多周期运动及混沌运动的时程曲线、相图、功率谱、Poincare截面图，同时得到了挠度-载荷、挠度-几何参数、挠度-材料参数等分叉图，考察了各种参数的影响。数值结果表明非线性粘弹性桩在一定的条件下可以通过倍周期分叉的方式进入混沌运动状态，且桩的载荷参数、几何参数、材料参数对其运动状态有较大的影响。     

热传导问题的非协调数值流形方法

数值流形方法通过引入数学与物理双重网格,将插值域与积分域分别定义在两个不同的覆盖上,其优点是网格划分随意,不受复杂边界形状和材料界面的限制,是较之于有限元方法更一般化的数值模拟方法。在计算精度方面,数值流形方法远远高于有限元法。但它的精度还是不够理想。为此本文在单元总体位移场上附加非协调位移基本项,使单元位移函数趋于完全,构造了非协调流形单元来改善流形单元的计算精度和计算效率,并将其应用于热传导问题,推导了势问题的非协调数值流形方法。     

某类非线性系统积分流形的研究——多频共振稳定性准则的积分流形方法证明

本文遵循Bogoliubov积分流形的思想,研究系统φ=w(x) φ(t, x, φ),=εX(t, x, φ)在x空间某点或某高维超曲面之ε领域内存在稳定积分流形之充分条件。当x空间某点为系统之孤立共振点或某高维超曲面为系统之孤立共振超曲面时,所得定理可证明文献[4]中有关多频共振稳定性判别准则之正确性。     

弹性力学的复变量无网格流形方法

为克服无网格流形方法配点过多、计算速度慢、容易形成病态方程组等缺点,将复变量移动最小二乘法与无网格流形方法相结合,提出了弹性力学的复变量无网格流形方法。分别采用线性基本与二次基进行计算,并与无网格流形方法相比。研究表明该方法计算量小、精度高。     

多孔介质非线性渗流问题的摄动解

考虑变形多孔介质渗透参数(渗透率和孔隙度)与孔隙压力呈负指数变化的特点，建立了多孔介质渗流问题的数学模型，采用积分变换方法求出了一维非线性渗流问题的摄动解，并对常数渗透参数和指数渗透参数的渗流问题进行对比分析，计算结果表明：两者之间的差别较大，且渗透参数的变化对于流体渗流中后期过程有着重要的影响，但对渗流早期影响不大，这对于定量研究工程中非线性渗流问题模型参数的相对重要性提供了可靠的理论依据。     

时滞对于参数激励系统周期运动的影响

本文以结构物承受周期激励地震波和具有时滞的弹性地基为背景，研究时滞对于结构物振动响应的影响规律。问胚的数学模型是一个非线性参数激励时滞系统，采用中心流形和平均法得到Hopf分岔方程，研究结果表明时滞量对结构物的地震响应有重要影响，它可使响应振幅增大，意味着结构物遭破坏的危险性增大，同时，结果也表明可以通过阻尼控制器克服时滞引起的危险性。最后，对于不同的时滞量，比较理论分析与数值模拟结果基本吻合，说明了本文结果的可靠性。     

计算结构动力响应的分段精细时程积分方法

利用钟万勰等发展的精细时程积分方法，提出了解线性定常结构动力系统响应的分段精细进程积分方法，它能适用于各种激励作用下系统的动力响应。对于载荷项采用线性和两次多项式进行拟合，采用精细时程积分方法和叠代方法对动力响应进行计算，与传统的离散积分方法如纽马克方法和威尔逊方法等及状态方程直接积分方法进行数值比较，本方法具有很高的精度和计算效率。     

多输入—多输出实验模态参数识别（1）：统一的数学模型

本文利用状态空间方程给出一般粘性阻尼线性振动系统的线性高阶微分多项式模型,及其所对应的离散差分模型.这个模型描述了系统多输入激励信号和多输出响应信号之间的关系.近十年来发展起来的多输入—多输出实验模态参数识别都是建立在这个模型之上的,在这个意义上,我们称之为多输入—多输出实验模态参数识别技术的统一数学模型.当考虑到测量频段外模态和测量数据中噪声的影响时,本文也给出了相应的统一数学模型.     

三维Laplace方程边界元中线性单元的精确积分法

边界元中的边界积分计算影响计算精度和计算速度。非奇异积分一般采用数值积分，当配置点接近积分单元时，计算精度降低。未知函数线性插值得到的解是连续解，但计算难度增大。本文采用积分区域变换，将三维Laplace问题的二维积分化为一维积分，这样奇异积分和非奇异积分能采用精确积分的方法计算，使求解精度，计算速度都得到提高。     

瑞典圆弧法的积分形式及其广义数学模型

为提高边坡稳定系数的精度、降低其计算复杂度并扩大其应用范围,基于应用较为广泛的极限平衡理论——瑞典圆弧法,转化传统数学模型的连加为积分,得到了积分数学模型;通过对同一边坡在强度折减模型、传统数学模型、积分数学模型这三种不同模型下的稳定系数进行对比,证明了积分公式相对传统公式在精度和计算复杂度上的优越性;并推导出了边坡稳定系数计算的广义数学模型,进而给出了成层土复杂边坡的简化计算方法。结果表明:积分数学模型不用考虑土条划分和预估滑动面形式,其计算值相对强度折减模型仅低4.8%,相对传统数学模型高9.17%;计算时间相对传统数学模型低40%,相对强度折减法仅高20%;广义数学模型在三个假定的前提下,解决了"异形坡面"边坡的稳定系数计算;成层土边坡的简化方法,一定程度上进一步解决了非连续成层土体边坡的稳定系数计算问题。与传统数学模型相比,积分数学模型提高了计算精度并减少了工作量,弥补了传统数学模型因土条宽度选择不当和提出的各种假设所带来的不足,验证了广义数学模型和成层土边坡的简化计算方法对于求解复杂边坡的稳定系数的可行性。     

结构阻尼时域本构模型及其应用

阻尼合金作为一种新型结构功能材料在不少领域已获应用,由于其阻尼值较大且随频率呈复杂交化关系,传统的线性粘性阻尼理论或经典的非频变结构阻尼理论难以精确地描述其耗能行为。本文应用粘弹性阻尼理论,根据阻尼合金储能模量和损耗因子在频域的实测数据．应用最优化方法拟合出标准线性体模型中的本构参数;根据积分形式的三参量本构关系和变形体虚功原理,推导出了有限元形式的动力学方程;讨论了三参数初值的选取;对包含卷积积分的有限元动力学方程通过数学推导将其化为三阶线性微分方程组,再转化为标准状态变量方程,应用数值求解。数值计算实例证明了所提方法的正确和有效性。