STABILITY ANALYSIS OF REUSABLE LAUNCH VEHICLE LANDING STRUCTURE$^{\bf 1)}$

近年来,包括中国在内的诸多国家相继开展垂直起降重复使用火箭的研究,运载火箭在平台上垂直着陆时的着陆稳定性为实现运载火箭重复使用的关键问题. 由于在运载火箭设计初期结构设计尚未完成,不具有供着陆稳定性分析的详细的动力学模型,难以开展着陆过程动力学仿真,故对运载火箭着陆稳定性评估方法的研究尤为必要. 本文基于广义碰撞定律,对二维运动模式下运载火箭与着陆平台的多点碰撞过程进行了分析,切向采用库伦摩擦模型给出了切向运动学恢复系数的表达式. 本文首先通过机械能约束和接触碰撞中的单边约束给出了一般运动形式下广义运动学恢复系数的值域,再对两种典型运动模式,给出了该两种典型运动模式下广义运动学恢复系数的值域. 然后考虑着陆腿中缓冲器的作用,将运载火箭与平台的碰撞近似为完全非弹性碰撞,得到了其广义运动学恢复系数,并结合运动学分析和能量法提出了一种基于碰撞后速度的着陆稳定性的判别方法. 最后以某型运载火箭着陆样机的参数为例,分析了碰撞前速度、着陆腿跨距、摩擦系数对着陆稳定性的影响,结果表明,本文提出的稳定性判别方法较能量法更为精确,可以考虑触地速度、角速度、摩擦系数等参数间的耦合关系.     

UNCONSTRAINED MODAL ANALYSIS OF DYNAMIC MODEL OF FLEXIBLE SPACECRAFT$^{\bf 1)}$

挠性航天器动力学建模中的挠性耦合影响系数是动力学建模中的重要力学概念,它反映了航天器姿态和轨道运动与挠性附件的弹性振动效应. 挠性耦合影响系数间的恒等式关系,即惯性完备性准则,是挠性航天器动力学模型降阶和模态截断的重要依据. 以中心刚体带挠性附件航天器为研究对象,采用约束模态和非约束模态法描述挠性附件结构变形,利用欧拉-拉格朗日方程建立挠性航天器的动力学模型. 基于 Hughes 的研究成果,对挠性航天器的非约束模态恒等式及其用于动力学模型降阶的惯性完备性准则进行了证明和应用研究. 探讨了两种动力学模型惯量间的关系,并利用约束模态惯性完备性准则,推导了非约束模态惯性完备性准则. 最后,对中心刚体带双侧太阳帆板和带单侧太阳帆板构成的挠性航天器模型进行数值仿真计算,求出挠性附件非约束模态平动耦合系数,分析了非约束模态特征值和平动耦合系数随着刚柔质量比的变化情况,并尝试用非约束模态惯性完备性准则的质量特征恒等式对挠性航天器模型进行了检验.     

DYNAMIC MODELING AND ANALYSIS OF THE LOWER LIMB PROSTHESIS WITH FOUR-BAR LINKAGE PROSTHETIC KNEE$^{\bf 1)}$

相比于单轴式膝关节,四连杆膝关节具有更好的仿生特性和运动安全性,因而在下肢假肢研究中得到广泛关注. 本研究以一款四连杆膝关节被动假肢为研究对象,主要关注足-地交互作用力以及膝关节单边接触力等强非线性因素对下肢假肢步态的影响. 为此,采用 Kelvin-Voigt 模型和库伦模型描述足-地接触力和摩擦力,并采用 Kelvin-Voigt 模型描述膝关节单边接触力,从而基于第一类拉格朗日方程建立假肢动力学模型. 本研究以步态实验测得的髋关节运动数据为模型的驱动信号,针对假肢的步态特征进行了数值分析. 计算结果显示,当膝关节液压阻尼器的刚度较小时,强非线性作用力会使假肢产生显著的亚谐波响应,进而导致步态周期失谐. 进一步研究发现,提胯行为能够避免步态周期失谐,这也为残疾人行走时的提胯等代偿行为提供了一种新的力学解释. 为了评价假肢步态与健康人实测步态的一致性,本研究进一步定义了步态相关系数并分析了膝关节液压阻尼器刚度、阻尼参数对相关系数的影响. 结果表明,通过合理的刚度、阻尼参数设计,两者步态的相关系数可达到 0.9 以上,这为四连杆膝关节被动假肢进一步优化提供了理论支撑.     

DYNAMIC MODELING OF NONIDEAL SYSTEM BASED ON GAUSS'S PRINCIPLE$^{\bf 1)}$

高斯原理给出了通过求函数极值、从可能运动中鉴别出真实运动的规则, 它可以使得多体系统动力学问题不需通过求解微分(代数)方程, 而是采用求解最小值的优化方法来解决, 从而提供了一种适用于优化算法的建模思路, 因此, 如何定义恰当的高斯拘束函数是动力学优化方法得以实现的前提. 对于理想系统而言, 约束对系统的作用可以通过约束方程来体现, 故高斯拘束可表达为系统质点加速度的函数, 系统的动力学问题因此可以描述为目标函数为高斯拘束函数、优化变量为质点加速度的约束最优化问题; 当系统中需要考虑干摩擦等非理想因素时, 部分相互作用不能被所定义的约束方程所涵盖而需要采用额外的物理规律来描述, 这种相互作用破坏了原有的针对理想系统的高斯拘束函数的极值特性. 基于变分类的高斯原理, 推导并证明了目标函数以理想约束力所表达的非理想系统的极值原理, 针对目前文献中用于非理想系统的高斯原理进行了讨论, 指出其实际为文中的极值原理在非理想约束力与理想约束力无明显关联时的一种特殊表达形式, 当非理想约束力与理想约束力有明显的函数关系(如库仑摩擦定律中滑动摩擦力与法向约束力间的线性关系)时, 该形式失效; 同时根据文中的极值原理, 得到了考虑库仑摩擦时非理想的多体系统动力学问题的优化模型. 例子中分析了优化模型及相应的线性互补性模型的关系, 分析发现在满足刚体滑动问题的唯一性条件下二者互为充分必要条件, 从而证明了文中优化模型的可靠性; 并采用优化计算方法进行了动力学模拟, 模拟结果显示了将高斯原理与优化算法相结合的可行性及有效性.     

DYNAMICAL CONTROL OF STABILITY FOR BIRKHOFFIAN SYSTEM$^{\bf 1)}$

本文研究 Birkhoff 系统和广义 Birkhoff 系统平衡稳定性的动力学控制. 首先建立系统的运动方程和平衡方程. 其次,研究 Birkhoff 系统中控制参数出现在 Birkhoff 函数中平衡稳 定性的动力学控制. 方法是通过选取控制参数使得 Birkhoff 函数 $B$ 成为定号函数,而其时间导数 $\dot {B}$ 为与 $B$ 反号的常号函数. 再次,研究广义 Birkhoff 系统平衡稳定性的动力学控制,通过选取 Birkhoff 函数或附加项中包含控制参数的方法,使得 Birkhoff 函数是定号函数,而其时间导数为反号的常号函数,从而控制系统的平衡稳定性. 最后举例说明结果的应用.     

自行车动力学建模及稳定性分析研究综述

自行车发明于两个多世纪前. 这一看似古老的交通工具在为人们提供出行便利的同时,其独特的运动特性及动力学性质 也吸引了来自数学、物理及力学等多个学科相关学者的兴趣. 大体上,自行车可以描述为具有 7 个自由度和 4 个非完整约束的多刚体系统. 但由于前后车轮之间复杂的运动耦合关系,使得自行车的约束方程和动力学模型变得异常复杂, 导致对自行车的稳定性存在一些模糊认识. 本文针对经典的 Carvallo-Whipple 自行车构型,系统回顾了历史上自行车动力学研究中的相关问题,这些问题包括:(1) 自行车在复杂曲面上的几何约束和非完整约束的数学描述;(2) 自行车系统内在的对称性及守恒量; (3) 自行车动力学的各类建模方法; (4) 自行车运动的相对平衡点及稳定性分析,包括水平面上的匀速直线运动及旋转对称曲面上的匀速圆周运动;(5) 影响自行车自稳定性的结构参数等. 本文最后对自行车动力学实验和控制方面的研究工作进行了回顾,并对自行车今后的研究给出了展望.      

REVIEW OF THE THIRD NATIONAL SYMPOSIUM ON EXPLOSION AND IMPACT DYNAMICS FOR YOUNG SCHOLARS$^{\bf 1)}$

简要介绍了第三届全国爆炸与冲击动力学青年学者学术研讨会的情况,概况总结了与会学者所做的报告内容. 报告包括 4 个特邀报告、19 个专题邀请报告和 17 个主题邀请报告,其中主题邀请报告分为爆轰和爆炸动力学、结构动力学和多尺度高性能计算、材料动力学和实验测试技术、复合材料结构的动态力学行为、轻质结构的吸能特性和优化设计等 5 个主题. 研讨会为从事爆炸与冲击动力学的青年学者开展学术交流、建立学术友谊、增强学术合作提供了优秀平台,对学科的建设和发展起到了积极推进作用.     

DYNAMIC RESPONSE OF SATURATED POROUS ELASTIC FOUNDATION UNDER POROSITY ANISOTROPY$^{\bf 1)}$

由于不同的沉积条件和应力状态, 天然土体通常表现出一定的各向异性特征. 文章研究地基上表面受温度载荷和机械载荷时, 孔隙率各向异性参数变化对饱和多孔弹性地基热-水-力耦合动力响应问题的影响. 基于Lord-Shulman广义热弹性理论, 结合孔隙率各向异性基本假设, 建立了孔隙率各向异性饱和多孔弹性地基热-水-力耦合动力响应模型, 利用正则模态法推导出无量纲竖向位移、超孔隙水压力、竖向应力和温度分布的解析表达式并加以图示. 正则模态法是一种利用加权残差求得解析解的方法, 相较于其他方法能快速求解偏微分方程. 当孔隙率各向异性参数为1时, 可将该各向异性耦合动力响应模型退化为热-水-力耦合动力响应模型验证该地基模型的合理性. 着重分析了孔隙率各向异性参数变化对不同物理量的影响. 结果表明: 孔隙率各向异性参数变化对物理量均有一定影响. 在地基上表面受温度载荷作用时, 对超孔隙水压力和竖向应力影响最为明显; 在地基上表面受机械载荷作用时, 对超孔隙水压力和温度影响明显. 整体而言, 无论地基上表面受何种载荷, 随着各向异性参数增大, 峰值逐渐减小, 在地基深度增加方向峰值所在位置向靠近地基上表面方向移动.     

ANALYSIS OF UNSTEADY CAVITATION FLOW OVER HYDROFOIL BASED ON DYNAMIC MODE DECOMPOSITION$^{\bf 1)}$

空化是船舶和水下航行体推进器中经常发生的一种特殊流动现象,它具有强烈的非定常性,空化的发生往往会影响推进器的水动力性能和效率. 为探究绕水翼非定常空化流场结构,本文基于 Schnerr-Sauer 空化模型和 SST $k$-$\omega $ 湍流模型,开展绕二维水翼非定常空化流动数值预报与流场结构分析. 通过将数值预报的空泡形态演变和压力数据与试验结果对比,验证了建立的数值方法的有效性. 并基于动力学模态分解方法对空化流场的速度场进行模态分解,分析了各个模态的流场特征. 结果表明,第一阶模态对应频率为 0,代表平均流场;第二阶模态对应频率约为空泡脱落频率,揭示了空泡在水翼前缘周期性地生长与脱落,第三阶模态对应频率约为第二阶模态的 2 倍,揭示了两个大尺度旋涡在水翼后方存在融合行为. 第四阶模态对应频率约为第二阶模态的 3 倍,具有更高的频率,表征流场中存在一些小尺度旋涡的融合行为. 最后对不同空化数下的空化流场进行了模态分解分析,发现脱落空泡的旋涡结构随着空化数的减小而增大,第二阶模态频率随着空化数的减小而减小.     

STUDY ON A DYNAMICS MODEL OF TAPPING MODE AFM AND ENERGY DISSIPATION MECHANISM$^{\bf 1)}$

轻敲模式下探针从远离到间歇性接触样品表面,是一个连续的能量耗散过程.针对该连续过程的能量耗散机理研究仅零星存在于各个文献之中,对于连续过程中各个阶段的能量耗散机理也没有一个系统的解释和实验验证.本文提出了新的位移激励下原子力显微镜探针-样品系统简化模型并得到了一维振子系统等效阻尼的计算方法,并通过该方法计算了探针在远离样品表面时的空气黏性阻尼和靠近样品时的空气压膜阻尼,分析了探针从远离样品到间歇性接触样品表面这一过程中的环境耗散机理变化,得到了原子力显微镜系统理论品质因数与探针工作位置的关系曲线;在此基础上设计了轻敲模式下的微悬臂梁扫频实验,得到了系统实验品质因数与探针工作位置的关系曲线,进而验证了理论模型的准确性. 本文通过对轻敲模式下AFM环境耗散机理进行理论分析和实验验证,希望可以对轻敲模式下AFM动力学特性及其阻尼作用机理有更近一步的认识,同时对微纳米机电系统 (MEMS/NEMS) 能量耗散机理的研究提供理论参考和实验方法.     

DYNAMIC MODELING AND STABILITY CONTROL OF FOLDING WING AIRCRAFT 1)

折叠翼飞行器在变形过程中,其动力学模型呈现多刚体、多自由度和强非线性特点,同时气动力/力矩、压心、质心和转动惯量等参数也会大幅度变化,严重影响飞行稳定性. 由此,本论文将对飞行器的多刚体动力学建模与变形稳定控制进行研究.基于凯恩方法建立了折叠翼飞行器的多刚体动力学模型,并从中得到了变形所产生的附加力和力矩表达式.通过气动计算拟合出气动参数与折叠角之间的函数关系,由此分析了不同折叠角速度下飞行器的纵向动态特性, 结果表明,折叠翼飞行器变形过程中速度、高度和俯仰角均会发生变化,飞行器无法保持稳定飞行.为此提出了一种基于自抗扰理论的飞行器变形过程中的稳定控制方法.将折叠翼飞行器纵向非线性动力学模型中存在的非线性项、耦合项以及参数时变项都视为系统内外总扰动,利用扩张状态观测器对总扰动进行实时估计和补偿, 针对补偿后的系统设计PD控制器,实现了速度通道和高度通道的解耦控制.通过Lyapunov稳定性原理证明了系统的稳定性, 并进行数学仿真验证. 仿真结果表明,基于自抗扰理论设计的稳定控制器能够解决飞行器变形所带来的强非线性和参数时变等问题,保证飞行器的高精度稳定控制.     

MULTI-MODE OSCILLATIONS AND HAMILTON ENERGY FEEDBACK CONTROL OF A CLASS OF MEMRISTOR NEURON$^{\bf 1)}$

根据法拉第电磁感应定律,在离子穿越细胞膜或者在外界电磁辐射下,细胞内外的电生理环境会产生电磁感应效应,继而会影响神经元的电活动行为. 基于此,本文考虑电磁感应影响下的 Hindmarsh-Rose (HR) 神经元模型,研究了其混合模式振荡放电特征,并设计一个 Hamilton 能量反馈控制器,将其控制到不同的周期簇放电状态. 首先,通过理论分析发现磁通 HR 神经元系统的 Hopf 分岔使其平衡点的稳定性发生了改变,并产生极限环,进而研究了 Hopf 分岔点附近膜电压的放电特征. 基于双参数数值仿真发现该系统具有丰富的分岔结构,在不同的参数平面上存在倍周期分岔、伴有混沌的加周期分岔、无混沌的加周期分岔以及共存的混合模式振荡. 最后,为了有效控制膜电压的混合模式振荡,利用亥姆霍兹理论计算出磁通 HR 神经元系统的 Hamilton 能量函数并设计 Hamilton 能量反馈控制器,通过数值仿真分析了膜电压在不同反馈增益下的簇放电状态,发现该控制器能够有效地控制膜电压到不同的周期簇放电模式. 本文的研究结果为探究电磁感应下神经元的分岔结构及其能量控制领域提供了有用的理论支撑.     

CONSTITUTIVE MODELS AND HYPERGRAVITY PHYSICAL SIMULATION OF SOILS$^{\bf 1)}$

数值模拟和物理模拟是分析土体沉降和稳定性的主要手段. 本构模型作为描述土体应力应变关系的数学表达式, 是数值模拟的基础. 土体具有碎散性, 这一基本物理特性导致了其具有压硬性、摩擦性和剪胀性, 这是土的力学特性区别于金属的主要特征, 在土体的本构模型中必须反映这3个基本特性. 传统土力学将土体的变形和强度分离考虑, 分别采用弹性理论和基于刚塑性模型的极限平衡理论分析, 虽然应用广泛, 但由于不能全面地反映土的基本力学特性, 计算结果的精度常常难以满足定量分析的需要. 剑桥模型作为第一个全面反映压硬性、摩擦性和剪胀性的弹塑性本构模型, 实现了变形和强度的统一, 能较好地描述饱和正常固结黏土的应力应变关系, 被视为是现代土力学的开端; 统一硬化模型通过引入一个独特的硬化参数进一步发展了剑桥模型, 将适用范围扩大到超固结黏土. 作者认为, 未来岩土体本构模型研究的挑战是: 如何考虑岩土体在受力过程中土骨架相变与多场耦合, 以解决目前本构模型尚无法定量分析的能源、交通、环境和水利相关的重大岩土工程问题. 超重力物理模拟具有缩尺效应和缩时效应, 克服了常重力物理模拟中模型的应力水平低于原型的缺点, 特别适用于大尺度、长历时问题的模拟. 相较数值模拟, 超重力物理模拟的优势在于能够检验本构模型的合理性, 揭示本构模型无法描述的未知特性. 最后, 介绍了采用数值模拟和物理模拟联合分析大直径钢管桩水平受荷特性的工程案例.     

DYNAMICS MODELING,SIMULATION, AND CONTROL OF ROBOTS WITH FLEXIBLE JOINTS AND FLEXIBLE LINKS$^{\bf 1)}$

本文探究了铰柔性对机器人动力学响应和动力学控制的影响. 首先, 建立由$n$个柔性铰和$n$个柔性杆组成的空间机器人模型, 运用递推拉格朗日动力学方法, 得到柔性机器人系统的刚柔耦合动力学方程. 在动力学建模过程中, 除了考虑杆件的拉伸变形、弯曲变形、扭转变形以及非线性耦合变形对机器人系统动力学行为的影响, 还考虑了铰的柔性对机器人动力学响应和控制的影响. 其中, 柔性铰模型是基于Spong的柔性关节简化模型, 将柔性铰看成线性扭转弹簧, 不仅考虑了铰阻尼的存在, 还考虑了柔性铰的质量效应. 其次, 编写了空间柔性铰柔性杆机器人仿真程序, 研究铰的刚度系数和阻尼系数对系统动力学响应的影响. 研究表明: 随着柔性铰刚度系数的增大, 柔性机器人的动态响应幅值减小, 振动频率变大. 随着柔性铰阻尼系数的增大, 柔性机器人的动态响应幅值减小, 振动幅值的衰减速度变快. 可通过调节柔性铰的刚度和阻尼来减小柔性铰柔性杆机器人的振动, 因此铰阻尼的研究具有重要工程意义. 最后, 研究了铰柔性在机器人系统动力学控制中的影响. 在刚性铰机械臂和柔性铰机械臂完成相同圆周运动时, 通过逆动力学方法求解得到两种情况下的关节驱动力矩. 研究表明: 引入柔性铰会使控制所需的驱动力矩变小, 对机器人控制的影响显著.     

BUFFER AND COMPLIANT DYNAMIC SURFACE CONTROL OF SPACE ROBOT CAPTURING SATELLITE BASED ON COMPLIANT MECHANISM$^{\bf 1)}$

研究了空间机器人在轨捕获非合作卫星过程避免关节受碰撞冲击破坏的缓冲从顺控制问题, 为此在机械臂与关节电机之间配置了一种柔性机构, 其作用在于: (1)在接触、碰撞阶段可通过其内置弹簧的变形来吸收被捕获卫星对空间机器人关节产生的冲击力矩; (2)在镇定运动阶段, 结合与之配合的缓冲从顺控制策略来适时开、关关节电机, 以保证关节受到的冲击力矩受限在安全范围. 首先, 利用多刚体系统理论获得配置柔性机构空间机器人及目标卫星分体系统动力学方程; 之后, 结合整个系统动量守恒关系, 捕获操作后系统运动几何关系及力的传递规律, 建立了两者形成联合体系统的动力学方程, 并计算了碰撞过程的冲击效应与冲击力. 为了实现失稳联合体系统的镇定控制, 提出了一种基于动态面的缓冲从顺控制方案. 上述控制方案可在实现吸收捕获操作产生的冲击力矩的同时, 还能在冲击力矩过大时适时开启、关闭关节电机, 以避免关节电机发生破坏; 此外, 动态面的引入避免了反演法存在的计算膨胀问题, 有效减少了计算量. 基于Lyapunov函数法证明了系统的稳定性, 并通过系统数值仿真结果验证了上述缓冲从顺控制策略的正确性.     

蒸汽管网的动力学建模与迭代方法1)

热电联产过程中的蒸汽输送过程难以监测,流量和压力变化比较复杂,管网运维依赖人员经验,容易造成能源浪费。本文建立了蒸汽供热管网的动态仿真模型,模型将蒸汽管网假设为单相可压缩流体网络,由节点和管路组成。使用管网质量方程、动量方程和能量方程建立管网物理模型,并采用隐式欧拉算法、稀疏矩阵算法求解由物理模型所得的非线性微分方程组。此外,本文建立了简单管网,并结合Fluent仿真对算法进行了验证。最后,建立了实际管网模型,与某企业蒸汽管网实际运行数据进行对比,结果表明所提出的模型不仅具有计算准确性,而且有实际工程应用价值。     

COORDINATED PATH PLANNING BY INTEGRATING IMPROVED RRT$^{\small*}$ AND QUARTIC SPLINE$^{\bf 1)}$

针对空间机器人抓捕空间非合作目标的在轨服务任务,同时考虑机器人运动学约束和动力学约束,提出一种分层式的自由漂浮双臂空间机器人协调路径规划方法. 首先,在路径规划层面上基于 RRT* 算法分别规划双臂末端执行器在笛卡尔空间下的初始可行路径,为双臂设置独立的采样空间,保证路径规划过程中双臂系统不发生自身碰撞. 然后,在轨迹规划层面上利用四次样条曲线平滑 RRT* 算法生成的初始路径,设计满足样条曲线的一阶、二阶及三阶微分连续约束,同时考虑机械臂末端执行器的初末速度约束条件、初始加速度约束条件,得到适合于空间机器人执行的动力学可行的平滑 轨迹.最后,计算所规划路径的最大速度、最大加速度与机械臂末端执行器物理极限值的比值,取最小上限,即为最少路径规划时间. 所提路径规划方法能够设计出满足特定路径点约束的协调路径,且所设计的路径考虑了机械臂的物理限制条件,通过对自由漂浮双臂空间机器人进行仿真试验,验证了所提路径规划算法的有效性.     

作旋转运动功能梯度材料矩形Mindlin板的刚柔耦合动力学建模1)

研究了作大范围运动功能梯度材料矩形板的刚柔耦合动力学问题。从连续介质力学理论出发,基于Mindlin板理论,采用无网格径向基点插值法(radial point interpolation method,RPIM)对矩形板变形场进行离散,考虑柔性板变形位移中二阶非线性耦合变形量,即横向弯曲引起的纵向缩短量,运用第二类拉格朗日方程推导了大范围运动功能梯度材料板的刚柔耦合动力学方程。分别采用一次近似耦合模型和传统零次近似耦合模型对不同转速下的功能梯度悬臂板进行了仿真,结果说明随着转速的提高,传统零次模型将发散,而一次近似模型依然收敛。将仿真结果与假设模态法和有限元法对比,验证了本文基于RPIM建立的动力学模型的正确性及在同等计算条件下的优越性。研究了功能梯度指数对作旋转运动的功能梯度矩形板动力学特性的影响。结果表明,随着功能梯度指数的增大,板的横向变形增大,且固有频率随之减小,说明功能梯度指数的增大会使系统柔性增加。同时,对中心刚体-功能梯度悬臂板的自由振动频率转向现象进行了讨论。     

NON-AXISYMMETRIC MIXED BOUNDARY VALUE PROBLEM FOR DYNAMIC INTERACTION BETWEEN SATURATED POROUS FOUNDATION AND RECTANGULAR PLATE$^{\bf 1)}$

在同一界面的不同区域具有多种边界条件, 称之为混合边界, 这是一个熟知的力学问题. 对这类问题进行精确分析时, 必须要进行混合边值问题的求解. 而对于一般的三维非轴对称情形, 混合边值问题的求解往往存在数学困难. 本文利用Hilbert定理和双重Fourier变换, 给出了一种求解三维非轴对称混合边值问题的解析方法, 利用该方法对具有混合透水边界的饱和多孔地基上矩形板的振动弯曲进行了解析研究(板与地基接触面为不透水边界, 其余为透水边界). 首先, 基于Kirchhoff理论和Biot多孔介质理论建立矩形板与饱和多孔地基的动力控制方程, 进行耦合求解. 针对板土接触面和非接触面的混合边值问题, 采用双重Fourier变换构造出两对二维对偶积分方程, 以接触应力和接触面孔隙压力为基本未知量, 用Jacobi正交多项式将未知量展开, 再利用Schmidt法对二维对偶积分方程完成求解, 最终推导出板土系统在动力作用下的位移和应力解析式. 通过将本文计算模型退化为单一弹性地基, 与已有研究结果进行对比, 验证了本文方法的正确性和有效性. 最后, 通过数值算例, 对饱和多孔地基上矩形板的动力响应及参数影响做出分析和讨论. 此外, 本文提出的解析法具有一般性, 可广泛应用于复杂接触问题和多场耦合问题的求解.     

NONLINEAR DYNAMIC RESPONSE OF PRE-DEFORMED BLADE WITH VARIABLE ROTATIONAL SPEED UNDER 2:1 INTERNAL RESONANCE$^{\bf 1)}$

在工程实际中,涡轮机叶片的转速在很多应用场景下不是一个定常值,比如发动机在启动、变速、停机等工况下,转子输入与输出功率失衡,伴随产生扭振,产生速度脉冲. 另外,由于服役环境、安装误差等因素会引起叶片在所难免的预变形. 本文主要研究预变形叶片,在变转速条件下的非线性动力学行为. 考虑叶片转速由一定常转速和一简谐变化的微小扰动叠加而成. 应用拉格朗日原理得到变转速叶片的动力学控制方程,并采用假设模态法将偏微分方程转为常微分方程,通过引入无量纲,使方程更具有一般性. 运用多尺度方法求解了该参激振动系统,得到了在 2:1 内共振情形下的平均方程,进而获得系统的稳态响应. 详细研究温度梯度、阻尼以及转速扰动幅值等系统参数对叶片动力学响应的影响规律,同时考察了立方项在 2:1 内共振下对方程的影响. 对原动力方程进行正向、反向扫频积分来观察其跳跃现象,并对解析解进行验证. 结果发现参数的变化对叶片均有不同程度影响,在 2:1 内共振下立方项对系统响应的影响很小,解析解与数值解吻合很好.