应用逐步加载法求解圆板弯曲非线性问题

本文采用逐步加载法将圆板弯曲的非线性微分方程组线性化,再用变分方法求解线性化方程.文中推得各次加载时的载荷与挠度,应力与挠度关系的递推公式.圆形薄板在轴对称弯曲情况下的非线性问题可用卡门方程表示     

金属材料的非线性疲劳累积损伤模型及强度衰减分析

用连续介质损伤力学的方法推导得到了一个非线性的损伤发展方程，该方程计及了加载过程应力与损伤的完全耦合效应；又由损伤发展方程得到了一个计及加载循环周内和循环周间损伤积累非线性效应的完全非线性疲劳累积损伤模型，已有的部分非线性疲劳累积损伤模型成为它的特例。由损伤发展方程确定了材料剩余强度随损伤发展的衰减规律和临界损伤值。     

基于弹性力学第一性原理的数据驱动力学建模

提出了一种基于弹性力学第一性原理的数据驱动力学建模方法,其能够从基于弹性力学方程的数值计算结果建立简洁且能准确捕捉变形机制的力学模型。基于有限元计算得到的高精度数据和无监督数据驱动控制方程识别方法Seq-SVF,从梁的载荷和位移数据中自动识别出了Timoshenko梁形式的弯曲控制微分方程,得到了三种不同加载条件下剪切影响系数关于结构尺寸和力学参数的函数表达式。揭示了经典模型适用的加载条件,同时还给出了一种未发现的新模型。通过将基于弹性力学的第一性原理计算与数据驱动范式相结合,克服了传统建模方法的局限性和对人类经验的强依赖性,为建立简洁的力学模型提供了一种新途径。     

电磁内爆套筒动力屈曲的实验和数值模拟

用实验及数值模拟方法研究了电磁内爆套筒的屈曲响应规律。用电容器组脉冲发生器装置作为驱动源,设计4种不同材料、不同尺寸的金属套筒,通过调节电容器组的充电电压得到不同的加载电流,研究材料及几何参数、加载脉冲特性对套筒屈曲的影响;采用瞬态非线性有限元方法对实验结果进行了数值模拟。研究表明,在同样加载电流作用下相同材料套筒形成的屈曲波数随半径/厚度比增大而增大,不同材料套筒的屈曲波数随屈服强度/塑性强化模量比值增大而增大,而相同材料的屈曲波数不随加载电流的大小而改变;模拟计算结果与实验结果基本吻合。     

气流作用下同轴带电射流的不稳定性研究

通过对气体驱动同轴电流动聚焦的实验模型进行简化,开展了电场力和惯性力共同作用下同轴带电射流的不稳定性理论研究.在流动为无黏、不可压缩、无旋的假设下,建立了三层流体带电射流物理模型并得到了扰动在时间域内发展演化的解析形式色散关系,利用正则模方法求解色散方程发现了流动的不稳定模态,进而分析了主要控制参数对不稳定模态的影响.结果表明,在参考状态下轴对称模态的最不稳定增长率最大,因此轴对称扰动控制整个流场.外层气流速度越高,气体惯性力越大,射流的界面越容易失稳.内外层液-液同轴射流之间的速度差越大,射流越不稳定.表面张力对射流不稳定性起到促进作用.轴向电场对射流不稳定性具有双重影响:当加载电场强度较小时,射流不稳定性被抑制;当施加电压大于某一临界值时,轴向电场会促进射流失稳.临界电压的大小与界面上自由电荷密度和射流表面扰动发展关系密切.这些结果与已有的实验现象吻合,能够对实验的过程控制提供理论指导.     

不确定结构动力区间分析方法研究

将结构中的不确定性参数用区间数来表示,建立了系统的区间动力特征方程和动力控制方程.通过将Monte Carlo方法引入到区间分析中来,并改进了一般Monte Carlo方法在区间模拟中存在的缺陷,对区间动力特征方程和动力控制方程进行了求解,得到了结构的固有频率、位移响应和应力响应的区间计算结果.算例表明改进的Monte Carlo方法求解简单,具有普适性,能够给出较高精度的响应区间.     

具有可变均压槽的气体静压推力轴承性能研究

具有可变均压槽的气体静压推力轴承,是作者提出的一种新结构轴承.本文建立了这种新结构轴承的性能计算方法,用有限差分法耦合求解气体静压推力轴承的雷诺方程和弹性薄片的变形控制方程,得到轴承的各项性能参数,并通过试验加以对比验证,计算结果与试验具有较好的一致性.结果表明该新型气体静压推力轴承的设计方案能够提高轴承的刚度.     

带有末端集中质量的双连杆柔性机械臂主动控制

对带有末端集中质量的双连杆柔性机械臂的主动控制进行了研究,给出系统的动力学方程,采用非线性解耦反馈控制方法分别得出系统大范围运动方程和柔性臂的动力学方程,采用机械臂逆动力学方法和LQR方法分别设计大范围运动控制律和压电作动器控制律.仿真结果显示,本文控制方法能够有效地进行机械臂的轨迹跟踪,柔性臂的弹性振动可以得到有效抑制.     

检验人工听骨力学性质的解析方法

传导性耳聋多由鼓膜或听骨链病变所致,人工中耳已在临床上广泛应用于传导性耳聋治疗.根据人工中耳结构试验和力学假设,用变分原理建立在任意时刻$t$点,人工鼓膜受迫振动位移与激发声压的关系,推导控制方程. 采用分离变量法及贝塞尔函数对方程求解,得到鼓膜和听骨位移的解析解. 通过与文献试验数据比较,感觉曲线形状与走势基本一致,证明了该方程是正确的,能够反映声音的传导机制. 利用方程计算得到:质量小的材料传音效果好于质量大的;顶盘面积一定时,质量为某一特定值时传音效果更佳; 该结果与相关文献的试验结果相符合. 研究表明: 该方程对于人工中耳数据的获得非常方便,并且能够优化人工听骨材料,特别是对质量的控制能达到科学的效果,避免测量因素及其他因素的干扰.该方法为耳生物力学的研究及应用提供了新的视点.      

压电智能环形板的主动控制

对在不同位置粘有任意多组压电传感器和压电执行器的轴对称弹性环形薄板的振动控制进行了研究．根据压电执行元件的等效作用量得到了压电智能环板的振动控制方程和传感方程,再利用分离变量法以及由传感器测得的电量和作用在执行器上电压之间的控制模式得到振动方程的全解．实行了对整体结构的主动控制．对不同的压电片布置进行了数值计算．结果表明：当离散分布压电元件布置越密,振动衰减的效果越佳     

随机荷载作用下随机结构线性反应的概率密度演化分析

提出了随机荷载作用下随机结构线性静力反应的概率密度演化方法。基于力学平衡方程,导出了随机荷载作用下随机结构反应的状态方程,进而引入扩展状态向量,建立了随机荷载作用下的随机结构静力反应的概率密度演化方程,讨论了其差分数值求解技术,进行了八层框架结构在随机荷载作用下的反应的算例分析,在单一随机参数结构的情况下,与随机结构反应的精确解答进行了对比;对于多个随机参数结构随机反应,则与MonteCarlo分析结果进行了比较,研究表明,本文提出的方法具有很高的精度及良好的实用性。     

Analytical modeling of static behavior of electrostatically actuated nano/micromirrors considering van der Waals forces

In this paper,the effect of van der Waals(vdW)force on the pull-in behavior of electrostatically actuatednano/micromirrors is investigated.First,the minimum potential energy principle is utilized to find the equation governing the static behavior of nano/micromirror under electrostatic and vdW forces.Then,the stability of static equilibrium points is analyzed using the energy method.It is foundthat when there exist two equilibrium points,the smaller oneis stable and the larger one is unstable.The effects of different design parameters on the mirror’s pull-in angle andpull-in voltage are studied and it is found that vdW forcecan considerably reduce the stability limit of the mirror.Atthe end,the nonlinear equilibrium equation is solved numerically and analytically using homotopy perturbation method(HPM).It is observed that a sixth order perturbation approximation can precisely model the mirror’s behavior.The results of this paper can be used for stable operation design andsafe fabrication of torsional nano/micro actuators.     

Exact three-dimensional analysis for static torsion of piezoelectric rods

Based on the theory of elasticity, exact analytical and numerical solutions of piezoelectric rods under static torsion are studied. In this paper, direct solution method is used. The main scope is to check the extension of validity of assumptions in previous papers that had been made based on linear distribution of electric potential through the cross section and their influences on deflection and the angle of rotation. Stress and electric induction functions are employed to obtain the exact solution of the static and electrostatic equilibrium equations under torsional loading. It is shown that previous assumptions are valid only in some types of piezoelectric materials, while in other types these assumptions lead to considerable deviations from accurate modeling. The present analytical solutions are compared with three-dimensional finite element analysis results and absolute agreements are found. At the end of this article, torsional rigidity, shape-effects on induced piezoelectric deformation and the range of valid region for linear distribution of electric potential assumption have been studied.     

基于高阶位移模式的压电层合板动力有限元模型

建立了含压电片层合板的有限元动力学模型。以位于压电层上下表面处的电场强度和层间电压为未知量,给出了三次函数的电势分布模式,采用Reddy的高阶剪切理论描述板的位移场,假设板厚度方向的正应力为零给出了减缩的本构方程,采用有限元方法,基于Hamilton原理导出结构的动力学方程,然后用静态缩聚的方法压缩掉电场自由度和次要的位移自由度。最后用四边形矩形单元求解了一对称铺层和非对称铺层悬臂板的固有频率,并与ANSYS结果对比验证了本文模型的精确性。     

Static and dynamic responses of a microcantilever with a T-shaped tip mass to an electrostatic actuation

In this study, nonlinear static and dynamic responses of a microcantilever with a T-shaped tip mass excited by electrostatic actuations are investigated. The electrostatic force is generated by applying an electric voltage between the horizontal part of T-shaped tip mass and an opposite electrode plate. The cantilever microbeam is modeled as an Euler–Bernoulli beam. The T-shaped tip mass is assumed to be a rigid body and the nonlinear effect of electrostatic force is considered. An equation of motion and its associated boundary conditions are derived by the aid of combining the Hamilton principle and Newton's method.An exact solution is obtained for static deflection and mode shape of vibration around the static position. The differential equation of nonlinear vibration around the static position is discretized using the Galerkin method. The system mode shapes are used as its related comparison functions. The discretized equations are solved by the perturbation theory in the neighborhood of primary and subharmonic resonances.In addition, effects of mass inertia, mass moment of inertia as well as rotation of the T-shaped mass, which were ignored in previous works, are considered in the analysis. It is shown that by increasing the length of the horizontal part of the T-shaped mass, the amount of static deflection increases,natural frequency decreases and nonlinear shift of the resonance frequency increases. It is concluded that attaching an electrode plate with a T-shaped configuration to the end of the cantilever microbeam results in a configuration with larger pull-in voltage and smaller nonlinear shift of the reso-nance frequency compared to the configuration in which the electrode plate is directly attached to it.     

The extended Durbin method and its application for piezoelectric wave propagation problems

Based on the well-known Durbin method, an efficient numerical method was developed for the inversion of the two-sided Laplace transform. The accuracy of the method was verified using examples. As an application of the method, transient elastic waves propagating in a two-layered piezoelectric medium subjected to anti-plane concentrated loading and in-plane electric displacement loading were investigated. One-sided and two-sided Laplace transforms were applied to determine the shear stresses and electric displacements in the double Laplace transform domain. Subsequently, the Durbin method for one-sided Laplace transform inversion and the extended Durbin method for two-sided Laplace transform inversion were used to implement the numerical inversions. Additionally, the numerical results of the transient stresses and electric displacements were evaluated and discussed. It showed that the arrival time of transient waves satisfies physical phenomena, and the transient solution oscillates near the static solution and rapidly approximates the static solution.     

有限到无限:基于离散静力平衡推导悬链线方程和梁位移方程

悬链线问题是一类非常经典的力学问题,推导出悬链线方程的方法从静力平衡到变分原理,非常之多。本文基于离散的静力平衡模型,采用从有限到无限的求解方法,建立并求解了悬链线方程。在此基础上,本文进一步建立了含有扭簧作用的离散静力平衡模型,推导并求解梁位移方程。最后在不同参数下求解梁位移方程得到其渐变解,并分析了多种平衡状态的存在情况。     

在中支座可滑动的两跨连续索的静力分析

基于非线性弹性理论,从单索的平衡出发,在索元的无应力长度已知情况下,提出了两跨连续索在中间节点滑移的初始平衡位置确定及均布荷载作用下的静力响应分析方法,建立了其非缄性计算的基本理论、公式和求解方法,并进行了实例计算。     

Equilibria of axially moving beams in the supercritical regime

Equilibria of axially moving beams are computationally investigated in the supercritical transport speed ranges. In the supercritical regime, the pattern of equilibria consists of the straight configuration and of non-trivial solutions that bifurcate with transport speed. The governing equations of coupled planar is reduced to a partial-differential equation and an integro-partial-differential equation of transverse vibration. The numerical schemes are respectively presented for the governing equations and the corresponding static equilibrium equation of coupled planar and the two governing equations of transverse motion for non-trivial equilibrium solutions via the finite difference method and differential quadrature method under the simple support boundary. A steel beam is treated as example to demonstrate the non-trivial equilibrium solutions of three nonlinear equations. Numerical results indicate that the three models predict qualitatively the same tendencies of the equilibrium with the changing parameters and the integro-partial-differential equation yields results quantitatively closer to those of the coupled equations.     

微重时卧圆柱内静液面方程与数值计算

由力学平衡方程导出了微重时卧圆柱内静液面方程，推出了边界接触角条件，使用Ｒｕｎｇｅ Ｋｕｔａ法在计算机上得出了数值结果，绘出了静液面形状，并对结果进行了分析．