奇怪吸引子分片机理的讨论

用已经提出的二维平面映射的奇怪吸引子结构讨论了这类平面映射的奇怪吸引子分片的机理，给出了有关判据。以Ｈｅｎｏｎ吸引子为例说明了这种机理的合理性．     

Lorenz吸引子的微观结构

轨迹保稳降维是一种分析高维非线性系统稳定性的方法。其要点是先在高维空问中求取轨迹；再将F轨迹映射为n-1个R^2映像，并在变换中严格保持感兴趣的稳定特性：分析各映像轨迹的稳定性：最后聚合为原轨迹特性的描述。本文按此分析Lorenz吸引子的结构稳定性。例如，将其中的z变量处理为时变参量后，(x,y）子系统成为时变的线性2维系统，可得分岔集{zcr}及奇点特性沿z轴的变化规律。故对于特定的轨迹(x,y,z），可将其在各坐标平面上的投影轨迹分成短线段的有序队列，各相邻线段对应于特性不同的奇点，从而揭示Lorenz吸引子全局分岔的精细结构及其通往高维混沌的道路。     

基于多边界切割插值的改进子模型分析方法

应用传统子模型技术进行结构的有限元分析,可以提高计算精度,得到局部的精确解,但仅适用于切割边界为体单元和壳单元的模型,当有其他类型单元通过切割边界时将产生很大的误差。本文对传统子模型技术及存在的问题进行了系统分析,提出了多边界插值方法和改进多边界插值方法,并应用于深水张力腿平台的整体模型和子模型分析,探讨了传统边界插值方法、多边界插值方法和改进多边界插值方法的计算精度。计算结果表明,相比传统边界插值方法,多边界插值方法具有更高的计算精度,随着内外切割边界距离的增大,误差逐渐减小并趋于稳定;改进多边界插值方法能进一步提高内外切割边界距离较小时的计算精度,可以满足梁单元通过切割边界的情况,拓宽了子模型技术的应用范围。     

广义van der Pol非线性振子的多吸引子共存和跳跃现象

高余维 经分岔和跳威风昌非线性动力学中的极为重要的研究内容。在本文里我们研究了广义ｖａｎ ｄｅｒ Ｐｏｌ非线性振子中的余维３经分岔、多吸引子共存和极限环振动的跳跃问题，用数值模拟方法研究了这个系统，数值结果证明了理论结果的正确性。     

基于单元子分法的结构多尺度边界单元法

建立在基于单元子分法的一种有效自适应格式以及多区域边界元三步求解技术基础上提出了一种计算结构多尺度问题的多区域边界元法。首先,通过高斯积分误差分析公式确定边界单元在满足精度要求下所需要的高斯点数,当所需高斯点数超过规定数目时该单元就被自动划分成一定数量的子单元,从而消除结构多尺度所引起的近奇异性。在单元子分技术的基础上采用多区域边界元三步求解技术来处理材料非均质问题:第一步消除各子域的内部未知量,第二步消除各子域独自拥有的边界未知量,第三步根据位移相容性条件和面力平衡条件建立系统方程组并求解公共界面节点位移以及每个子域的其他未知量。数值算例结果表明本方法可以用较少的计算时间得到满意的结果,是处理结构多尺度问题的一种有效方法。     

准周期激励Duffing系统中奇异非混沌吸引子的判别方法

奇异非混沌动力学是非线性动力学领域中的新课题。本文以准周期激励Duffing振子为例,对其产生的奇异非混沌吸引子(strange nonchaotic attractors, SNAs)进行分析。通过三维庞加莱截面和定量方法如傅里叶变换、李雅普诺夫指数、李雅普诺夫维数、关联维数和盒维数检测SNAs是否存在。研究结果表明,傅里叶变换无法判断混沌与奇异非混沌行为。而李雅普诺夫指数、李雅普诺夫维数可以作为检测系统混沌与非混沌指标。关联维数和盒维数显著表明系统奇异与非奇异性,从而阐明适用于准周期驱动Duffing振子中存在SNAs的判别方法,并为其他类似系统检测SNAs提供指导。     

约束多体系统动力学的局部子空间分解法

由约束矩阵及其正交补的列矢量张成约束多体系统的局部约束子空间和局部切向子空间．将系统的广义加速度和广义力沿这两个子空间分解,推导描述系统运动的纯微分方程和求约束反力的公式,并给出了违约修正的一种有效方法     

基于Newmark格式的多柔体系统子循环算法研究

提出一种基于Newmark积分格式的FMD(Flexible Multi-body Dynamics)子循环算法,可用于求解缩并型的FMD方程.推导了FMD方程中快变分量和慢变分量的同步更新及子步更新公式;在积分过程中,采用能量平衡计算校验算法的稳定性;结果表明该算法可以在保持合适的精度要求下,有效地提高计算效率.     

黏弹性多层介质中SH波动的一种吸收边界条件

提出一种高精度时域吸收边界条件,与有限元法结合用于模拟瞬态标量SH波在达朗贝尔黏弹性多层介质中传播问题.建立时域吸收边界条件的过程是:首先将半无限域沿着竖向半离散得到半离散的位移方程以及人工边界处的力-位移关系,再通过引入模态分解,将物理空间下的量转化到模态空间,从而获得半无限域模态空间下的频域动力刚度;其次采用一种在全频范围内收敛的连分式精确逼近单层介质模态空间下标量形式的频域动力刚度,将标量连分式扩展为矩阵形式用来表示多层介质的频域动力刚度;最后通过引入辅助变量技术,将模态空间下基于连分式的频域动力刚度关系转化为时域吸收边界条件,进一步转换到物理空间后得到物理空间下的时域吸收边界条件.单层介质和五层介质的数值算例表明,建立的高精度时域吸收边界条件对于达朗贝尔黏弹性单层介质是精确且稳定的;对于达朗贝尔黏弹性多层介质,为了保证其高精度特性,需要将人工边界放置在距离感兴趣区域约为0.5倍无限域高度的位置处.     

黏弹性多层介质中SH波动的一种吸收边界条件

提出一种高精度时域吸收边界条件,与有限元法结合用于模拟瞬态标量SH波在达朗贝尔黏弹性多层介质中传播问题.建立时域吸收边界条件的过程是:首先将半无限域沿着竖向半离散得到半离散的位移方程以及人工边界处的力-位移关系,再通过引入模态分解,将物理空间下的量转化到模态空间,从而获得半无限域模态空间下的频域动力刚度;其次采用一种在全频范围内收敛的连分式精确逼近单层介质模态空间下标量形式的频域动力刚度,将标量连分式扩展为矩阵形式用来表示多层介质的频域动力刚度;最后通过引入辅助变量技术,将模态空间下基于连分式的频域动力刚度关系转化为时域吸收边界条件,进一步转换到物理空间后得到物理空间下的时域吸收边界条件.单层介质和五层介质的数值算例表明,建立的高精度时域吸收边界条件对于达朗贝尔黏弹性单层介质是精确且稳定的;对于达朗贝尔黏弹性多层介质,为了保证其高精度特性,需要将人工边界放置在距离感兴趣区域约为0.5倍无限域高度的位置处.     

多柔体系统响应计算的子循环计算方法研究

过去近30年中,柔性多体系统动力学研究取得了巨大的进展,人们的兴趣集中在柔性多体系统建模、计算及实验研究等3个方面. Belytschko等于1979年提出的子循环算法已经成功地应用于结构动力响应的有限元计算中,然而有关柔性多体动力学的子循环算法研究尚未见报道. 该文提出了一种适合于柔性多体系统响应计算的中心差分类子循环算法,在将非线性微分-代数混合方程组(DAEs)缩并为纯微分方程组(ODE)的基础上,推导出快、慢变分量的同步更新公式和子步更新公式;在变量的数值积分过程中,采用能量平衡计算检查算法的稳定性;算例结果表明该算法可以在保持合适的精度要求下,有效地提高响应的计算效率;对积分步长进行摄动修正可以保持算法的稳定性.      

考虑可制造性约束的声子晶体多目标拓扑优化

对声子晶体进行拓扑优化可得到具有目标带隙特性的声子晶体结构, 在减振、隔声等领域具有潜在应用价值. 然而, 优化结果常会出现孤立材料单元而导致的制造困难问题. 针对该问题, 本文提出协同考虑带隙性能和可制造性约束的二维多相声子晶体多目标拓扑优化方法. 以在特定频率段带隙最宽和结构质量最小为优化目标, 在对微结构进行连通性分析的基础上, 引入考虑可制造性因素的附加约束, 并利用有限元法和具有精英选择策略的非支配排序遗传算法对该优化模型进行数值求解. 通过数值算例验证了本文模型及方法的合理性和有效性. 结果表明, 附加可制造性约束的拓扑优化模型可有效避免二维声子晶体构型中出现孤立材料单元的情况, 优化结果在满足带隙性能预期指标的同时也可兼顾到实际可制造性要求. 与仅仅考虑带隙性能的单目标优化结果相比, 本文提出的同时兼顾带隙性能和可制造因素的多目标优化模型可以针对实际应用和制造条件, 实现不同目标间的平衡, 具有显著优势和良好的应用前景.