参外联合激励下的簇发振荡及其机理

当周期激励频率远小于系统固有频率时，会存在快慢耦合效应，与单项激励不同，参外联合激励不仅会导致快子系统平衡曲线和分岔行为的复杂化，也会产生一些特殊的非线性现象，为此，本文以两耦合Hodgkin-Huxley细胞模型为例，引入周期参外联合激励，探讨在频域不同尺度耦合时该系统的簇发振荡的特点及其分岔机制．通过建立相应的快慢子系统，得到慢变参数变化下的快子系统的各种分岔模式以及相应的分岔行为，结合转换相图，揭示耦合系统随激励幅值变化时的动力学行为及其机理．研究表明，在激励幅值较小时，系统表现为概周期振荡，两频率分别近似于快子系统平衡曲线由Hopf分岔引起的两稳定极限环的振荡频率．概周期解随激励幅值的增加进入簇发振荡，导致这些簇发振荡的主要原因是在慢变参数变化的部分区间内，存在唯一稳定的平衡曲线，使得系统的轨迹逐渐趋向该平衡曲线，产生沉寂态，并随着慢变参数的变化，由分岔进入激发态．同时，快子系统中参与簇发振荡的稳定吸引子随激励幅值的变化也会不同，导致不同形式的簇发振荡．另外，与单项激励下的情形不同，联合激励时快子系统的部分稳定吸引子掩埋在其它稳定吸引子内，从而失去对簇发振荡的影响．     

HyperFLOW软件非结构网格亚跨声速湍流模拟的验证与确认

计算流体力学(computational fluid dynamics,CFD)数值模拟在航空航天等领域发挥越来越重要的作用,然而CFD数值模拟结果的可信度仍然需要通过不断地验证与确认来提高.本文给出了从制造解精度测试、简单到复杂外形湍流模拟网格收敛性研究等三个方面开展CFD软件验证与确认的方法,并对自主研发的CFD软件平台HyperFLOW在非结构网格上模拟亚跨声速湍流问题的能力进行了验证与确认.首先通过基于Euler方程和标量扩散方程的制造解精度测试,分别验证了HyperFLOW在非结构网格上对Euler方程和黏性项的求解精度,结果表明其能够在任意非结构网格上达到设计的二阶精度. 其次,通过NASATurbulence Modeling Resource中的湍流平板、二维翼型近尾迹流动、二维Bump等几个典型的亚声速湍流算例的网格收敛性研究,量化考察了数值结果的观测精度阶和网格收敛性指数,并与国外知名CFD解算器CFL3D,FUN3D的计算结果进行了对比,验证了HyperFLOW对简单湍流问题的模拟能力,且具有良好的网格收敛性和计算精度(阶). 最后,通过NASA CommonResearchModel标模定升力系数的网格收敛性研究和升阻极曲线预测,验证了软件在复杂外形亚跨声速湍流流动数值模拟中也具有良好的可信度.      

基于多尺度方法的平动圆柱贮箱航天器刚——液耦合动力学研究

以充液航天器为工程背景,借助多尺度方法研究刚--液耦合动力学系统非线性动力学特性.利用多维模态方法,将描述横向外激励下圆柱贮箱中液体非线性晃动的自由边界问题转换为液体模态系数相互耦合的有限维非线性常微分方程组.推导液体晃动产生的作用于贮箱壁的晃动力和晃动力矩的解析表达式,进而建立航天器刚体部分平动和液体晃动耦合的非线性动力学方程组.应用多尺度方法对刚--液耦合系统的动力学特性进行解析分析,通过固有频率的特征方程求解耦合系统固有频率,推导外激励频率接近耦合系统第一阶固有频率时液体晃动稳态解的幅值频率响应方程.结合数值方法,研究了液体晃动稳态解的幅值频率响应曲线和激励--幅值响应曲线.结果表明, 随充液比变化,液体晃动稳态解的幅值频率响应曲线会发生软、硬弹簧特性转换现象和"跳跃"现象;幅值频率响应曲线的软、硬弹簧特性转换点受重力加速度和弹簧刚度系数影响;以上所得研究结果表明,考虑非线性效应时的刚--液耦合系统动力学特性与传统的线性系统模型所显示的动力学特性具有本质区别.本文的研究工作对进一步分析充液航天器刚--液耦合非线性动力学特性具有重要参考价值.      

船舶与海洋平台结构冰载荷的高性能扩展多面体离散元方法

船舶与海洋平台结构的冰载荷是寒区海洋工程结构物设计中的关键参数,而离散元方法是有效计算结构冰载荷的重要手段. 本文采用基于闵可夫斯基和原理的扩展多面体离散元方法模拟船舶与海洋平台结构的相互作用过程. 其中,构造扩展多面体的近似包络函数并建立了基于优化模型的快速接触搜索算法;考虑单元间粘结作用的刚度软化过程建标识码元间的粘结-破碎模型. 同时,发展了 CPU-GPU 协同异构环境下的高性能并行算法. 为分析海冰与海洋结构作用中的冰载荷,采用ISO标准验证了扩展多面体离散元分析结构冰载荷的准确性. 采用离散元方法计算了船舶结构的冰载荷,研究了船舶结构表明的线载荷分布特点,并采用船舶结构冰阻力经验公式验证了计算结果的合理性. 采用离散元方法计算了平整冰区与多桩腿平台结构的相互作用,分析各桩腿上的冰载荷特点. 针对碎冰区的海冰管理过程,采用离散元方法分析了船舶结构绕行过程中的船舶和海洋平台结构冰载荷. 本文方法可有效应用于海洋结构冰载荷分析,能为极地船舶与海洋平台结构的设计和安全运行提供科学的分析手段.      

具有恒定冲击载荷的梯度泡沫金属材料设计

多胞材料可通过大变形大量地吸收冲击能量，引入密度梯度可进一步提高其耐撞性。梯度多胞材料的宏观力学响应对材料密度分布极为敏感，不同类型的细观构型的影响也极为不同。已有的研究工作主要局限在对给定的密度梯度分析其动态响应，较少对耐撞性设计方法进行研究。本文针对梯度闭孔泡沫金属材料，基于非线性塑性冲击波模型发展了耐撞性反向设计方法，以维持冲击物受载恒定为目标，运用级数法获得了简化模型和渐近解。利用变胞元尺寸法构建了连续梯度变化的三维Voronoi细观有限元模型，并利用ABAQUS/Explicit有限元软件对理论设计进行数值验证。结果表明，反向设计理论简化模型的渐近解对于梯度闭孔泡沫金属材料的耐撞性设计是有效的，所提出的耐撞性设计方法在控制冲击吸能过程和冲击物受载方面具有指导意义。     

不同泄压条件对方管内爆炸压力特性的影响

为研究泄压膜约束条件对甲烷/空气预混气体爆炸压力特性的影响，在方形火焰燃烧传播测试管道中布置压力传感器，开展不同泄压膜材料、泄压膜层数及泄压口位置实验。结果表明:牛皮纸和聚丙烯薄膜约束泄爆过程中，每增加一层泄压膜，管道内最大泄爆压力平均上升11.2%和12.3%。各强度泄压膜约束条件下，管道内最大泄爆压力随着泄压口位置接近点火端，均呈现Z形规律，当泄压口设置在距尾部端面0.25 m时，各曲线达到最小值，当泄压口设置在距尾部端面0.50 m时，各曲线出现最大值。     

非贯通裂隙岩体损伤演化率相关性及变形特征

含非贯通裂隙岩体是自然界中岩体的主要赋存形式，其裂隙几何特征对岩体的强度及变形均产生显著影响。应变率对岩体的损伤演化及黏滞效应也具有显著的率相关性。首先，运用模型元件的方法，将非贯通裂隙岩体动态破坏过程视为具复合损伤、静态弹性特性、动态黏滞特性的非均质点组成，对黏弹性响应的Maxwell体进行改进，将细观损伤体与裂隙损伤演化的宏观损伤体根据等效应变假设并联组成宏细观复合损伤体，构建综合考虑岩体宏细观缺陷的动态损伤模型；其次，基于断裂力学及应变能理论，对岩体宏观裂隙动态扩展的能量机制进行分析，综合考虑初始裂隙应变能、裂隙动态损伤演化过程应变能、裂隙闭合应变能，得到裂隙岩体宏观动态损伤变量计算公式；最后，将模型计算结果与实验结果进行比较，模型计算结果与实验结果吻合较好，证明了模型的合理性，同时利用模型讨论了裂隙倾角、应变率、岩石性质对岩体变形特征的影响规律。     

基于自适应参数估计的三轴磁传感器实时校正方法

针对传统基于椭球模型三轴磁传感器无法实时校正的问题,设计了基于自适应参数估计的三轴磁传感器实时校正方法,通过对参数的实时估计与处理,实现了三轴磁传感器误差的校正。首先,对三轴磁传感器全误差模型进行分析与处理,建立了基于模值的参数估计方程;其次,详细分析了参数模型的噪声特性,针对非高斯状态相关噪声,提出了采用自适应参数估计方法;最后,设计仿真与实验分析,完成三轴磁传感器误差校正。实验结果表明,所提出的算法可以实现误差参数的实时校正,提升了磁传感器误差校正的自主性。同时,基于自适应参数估计方法,使得校正结果模值标准差相较于传统方法减小了3倍。     

陀螺电机用H型动压气体轴承启动摩擦特性分析

为减小陀螺电机用H型动压气体轴承启动过程中的摩擦与磨损,建立了考虑表面粗糙接触与润滑的H型动压气体轴承模型,对其启动摩擦特性进行研究。首先,通过求解轴粗糙表面在转子重力作用下的弹性变形,确定启动前转子启动位置和启动力矩。然后,对表面粗糙接触与润滑方程、转子五自由度运动方程进行联立求解,得到启动过程中轴承承载力、运动轨迹、接触面积等特性,并通过接触面积减为0来判断转子的浮起时间和浮起转速。最后,研究轴承表面粗糙度和腔型结构等参数对轴承启动摩擦特性的影响。分析结果表明:减小粗糙表面峰顶标准差可降低浮起转速、减小磨损范围;随着腔宽比的增大,浮起转速增大、磨损范围减小;而随着腔深的增大,浮起转速和磨损范围均先减小后增大,腔深为1.0μm和2.0μm时浮起转速和磨损范围分别最小。     

一种考虑界面相物性影响的界面单元法

界面是由复杂的界面相简化而成的,界面破坏实际是界面相材料的破坏。数值计算为了方便,如经典模型和内聚力模型等,都把很薄的界面相作无厚度化处理。导致只能考虑界面的面力,而无法考虑界面相内的应力(平行于界面方向的应力)。使界面失效准则先天性地排除了界面相内部应力的影响,从界面相材料失效机理的角度来看这是不够严谨的。本文将界面相材料等效为一种弹性连续体,由界面本构关系推导得到了一种新的界面单元。该单元具有界面参数易确定、对界面相物性可以进行等效描述等优点。通过商用有限元软件ABAQUS和用户子程序UEL实现了数值分析,并与直接物理模型的数值模拟结果进行对比,证明了本方法的简便及准确性。通过对不同界面相厚度结构的进一步分析,探讨了本文方法的可行范围。