一类五阶非线性波方程的新精确解

通过利用新的G展开法,并借助Mathematica计算软件,研究了一类五阶非线性波方程的精确解,获得了方程的含有多个任意参数的新的显式行波解,分别为三角函数解、双曲函数解、指数函数解,扩大了该类方程的解的范围.     

一种压电驱动的三足爬行机器人

机器人领域涉及到力学、机械、材料、控制、电子和计算机等多个学科. 其中, 爬行机器人可在极端环境下工作, 进而可有效降低人工作业的危险性并提高工作效率. 因此, 爬行机器人一直是机器人领域的重点研究对象. 压电陶瓷是一种能够将机械能和电能互相转换的新型功能陶瓷材料. 逆压电效应是指当在电介质的极化方向施加电场, 这些电介质就在一定方向上产生机械变形或机械压力, 当外加电场撤去时, 这些变形或应力也随之消失. 本文基于压电陶瓷的逆压电效应设计了一种由3条弯曲变截面梁支撑的一体化三足爬行机器人. 利用理论力学方法对该三足爬行机器人建立整体受力分析方程, 再用哈密顿原理对变截面、变角度梁建立动力学方程, 最终得到了可求解该三足爬行机器人的压电驱动腿固有频率的方程. 设计并制作了三足爬行机器人实物, 通过实验测试了不同弯折角度、不同驱动频率、不同负载、不同电压波形对运动方向及运动速度的影响. 最后利用不对称的驱动电压使三足爬行机器人实现了左转、右转以及不加导轨的近似直线运动, 实现了设计的3个方向的运动, 最后分析了该机器人的能耗问题. 该研究可为微型爬行机器人设计和实验提供参考依据.      

三维位势问题的梯度边界积分方程的新解法

三维位势问题的边界元分析中,关于坐标变量的边界位势梯度的计算是一个困难的问题. 已有一些方法着手解决这个问题,然而,这些方法需要复杂的理论推导和大量的数值计算. 本文提出求解一般边界位势梯度边界积分方程的辅助边值问题法. 该方法构造了与原边界值问题具有相同解域的辅助边值问题,该辅助边值问题具有已知解,因此通过求解此辅助边值问题,可获得梯度边界积分方程对应的系统矩阵,然后将此系统矩阵应用于求解原边值问题,求解过程非常简单,只需求解一个线性系统即可获得原边值问题的解. 值得注意的是,在求解原边值问题时,不再需要重新计算系统矩阵,因此辅助边值问题法的效率并不很差. 辅助边值问题法避免了强奇异积分的计算,具有数学理论简单、程序设计容易、计算精度高等优点,为坐标变量梯度边界积分方程的求解提供了一个新的途径. 3个标准的数值算例验证了方法的有效性.      

基于局部加密纯无网格法非线性Cahn-Hilliard方程的模拟

为数值求解描述不同物质间相位分离现象的高阶非线性Cahn-Hilliard(C-H)方程,发展了一种基于局部加密纯无网格有限点集法(local refinement finite pointset method,LR-FPM).其构造过程为:1)将C-H方程中四阶导数降阶为两个二阶导数,连续应用基于Taylor展开和加权最小二乘法的FPM离散空间导数;2)对区域进行局部加密和采用五次样条核函数以提高数值精度;3)局部线性方程组求解中准确施加含高阶导数Neumann边值条件.随后,运用LR-FPM求解有解析解的一维/二维C-H方程,分析粒子均匀分布/非均匀分布以及局部粒子加密情况的误差和收敛阶,展示了LR-FPM较网格类算法在非均匀布点情况下的优点.最后,采用LR-FPM对无解析解的一维/二维C-H方程进行了数值预测,并与有限差分结果相比较.数值结果表明,LR-FPM方法具有较高的数值精度和收敛阶,比有限差分法更易数值实现,能够准确展现不同类型材料间相位分离非线性扩散现象随时间的演化过程.     

大型偏心结构吊装欠约束问题研究

吊装施工过程中被吊模块的水平度是作业要求的重要指标,通常需要增加配重调平。传统有限元方法需要补充约束以消除单元刚体位移,且需要重复计算平衡方程来求解调平载荷,效率不高。将模块的运动分解为随动坐标系的整体运动以及相对该坐标系的弹性变形,可将欠约束问题化为多体系统的静平衡问题。基于虚功率原理推导了吊装平顺时刻的节点力平衡方程以及相应的切线刚度矩阵,并将配重表示为基础配重与载荷系数相乘的形式。通过对节点力平衡方程求导,得到一组以载荷系数为自变量的微分方程,通过求解微分方程并结合水平度判据,可快速搜寻满足水平度要求的载荷系数。数值算例表明,该方法在解决偏心模块吊装欠约束问题方面具有明显的优势,在确定配重载荷方面具有较快的速度和合理的精度。     

自相似Euler方程的数值方法

Euler方程某些问题的解具有自相似特点，可以使用更准确的方法求解.提出了两种数值方法，分别称为自相似和准自相似方法，新方法可以使用现有守恒律方程的数值格式，无须设计特殊方法.对一维激波管问题、二维Riemann问题、激波反射以及激波折射问题进行了数值计算.对自相似Euler方程，一维计算结果显示数值解基本等同于精确解，二维结果也比现有文献计算的结果有更高的分辨率.对准自相似Euler方程，新方法可以求解不具有自相似性但接近自相似的问题，并在计算时间足够长时可以取得自相似Euler方程的效果.数值求解自相似Euler方程对自相似问题的研究，高分辨率、高精度格式的设计乃至Euler方程的精确解都有重要启示.      

有耗介质层上多导体传输线的电磁耦合时域分析方法

目前,针对空间电磁场作用有耗介质层上传输线的电磁耦合,仍缺乏有效的数值分析方法.因此,本文提出一种高效的时域混合算法,很好地解决了有耗介质层上传输线电磁耦合建模难的问题.首先,对经典传输线方程进行改进,推导了适用于有耗介质层上多导体传输线电磁耦合分析的修正传输线方程.然后,结合时域有限差分方法和相应插值技术,求解修正传输线方程,获得多导线及其端接负载上的电压和电流响应,并实现空间电磁场辐射与多导线瞬态响应的同步计算.最后,通过相应计算实例的数值模拟,与CST软件的仿真结果进行对比,验证了时域混合算法的正确性和高效性.     

受限亚音速气流中倒置悬臂壁板静气弹稳定性的理论及实验研究

板壳结构在航空航天、高速列车、能量采集等诸多工程领域已经得到了广泛应用. 将悬臂壁板倒置于轴向气流中并在壁板周围流场中设置刚性壁面可有效地调控壁板的失稳速度, 是俘能器优化设计的重要措施之一. 但针对刚性壁面作用下亚音速气流中倒置悬臂壁板的失稳机制仍需要开展深入研究. 本文以受限亚音速气流中倒置的二维悬臂壁板为对象, 以理论分析及风洞实验为手段, 研究了单侧刚性壁面效应对倒置悬臂壁板静态失稳特性的影响规律. 在理论分析中, 首先应用镜像函数法来处理壁面约束条件, 基于算子理论研究获得了以Possio积分方程为表征的壁板气动力, 壁面效应实际表征为一包含移位Tricomi算子的复合算子; 然后将壁板失稳方程的求解问题转化为定区间上的函数逼近问题; 最后, 依据Wererstrass定理并利用最小二乘法求解该最优函数, 以获得系统的失稳临界参数. 在试验研究中依据压杆稳定原理设计了壁板静态失稳的测试方法并完成了风洞实验. 理论分析结果表明, 壁板会发生发散(静气动弹性)失稳, 临界动压随壁板与壁面间距的增加而增大并最终趋于稳定(无壁面情况); 通过理论与风洞实验结果的对比分析, 验证了本文气动力及理论分析的适用性及准确性. 针对倒置悬臂壁板结构的气动弹性失稳问题, 本文提出的方法不涉及系统方程的离散及特征值求解问题, 而是将其转化为了定区间上的函数逼近问题进行求解, 这为弹性结构静气动弹性失稳问题的研究提供了一个可行的新思路.      

轴流压气机势流和尾迹相互干涉叠加效应分析

尾迹与势流的相互作用对压气机叶片表面的压力波动强度和出口的轴向速度分布及波动会产生很大影响。针对两级跨声速轴流压气机,采用数值方法研究两排转子叶片周向相对位置不同情况下,叶片表面压力波动抵消区域位置发生的变化。研究结果表明:转子叶片周向位置的变化会对静子叶片尾迹形态产生影响,会导致叶片表面压力波动抵消区域发生变化。伴随转子叶片周向相对位置的变化,存在上游转子叶片尾迹减缓下游静子叶片尾迹衰减的现象,从而改变了叶片排之间的相互干涉强度,并使叶片表面压力波动强度发生变化。     

基于耦合Ginzburg-Landau方程的石墨烯锁模光纤激光器的暗暗孤子解及其稳定性

以耦合Ginzburg-Landau方程为模型,对石墨烯锁模光纤激光器中脉冲的动力学行为进行了详细地研究.首先,采用拟解法推导了该方程组的啁啾类暗暗孤子对解;其次,对啁啾类暗暗孤子对的传输特性及稳定性进行了数值研究.结果表明:合理选择系统参数,不仅可以得到啁啾类暗暗孤子对的精确解,而且该精确解还可以在一定的扰动下稳定传输.相关结果将为设计和改善石墨烯锁模光纤激光器的性能提供一定的理论依据.