一种直降式滑橇减震器性能研究

为验证无人机伞降回收系统的直降式滑橇着陆装置的动力学性能,对其着陆缓冲过程进行数学模拟构造了该装置的动力学模型,通过数值求解得出了该装置缓冲过程的动态特性。对其着陆压缩过程的动力学模型进行了简化,计算结果可用于减震器参数的确定。落震试验结果表明:动力学模型能准确模拟着陆装置的缓冲过程;以5m/s着陆时过载仅为6.8g;本文所采用动力学模型的数值求解结果与实验结果一致,误差小于5%。     

一种微机电非线性耦合系统奇点稳定性研究

工程中许多微机电系统都采用电容驱动原理，这类结构实际上存在着强烈的静电和机械两个物理场的非线性耦合，因此系统的动态特性比较复杂。本文基于一种扭转微镜系统，通过数值分析方法，研究其非线性动态特性，根据理论分析和数值计算证明该系统在相平面中存在两个奇点，一个是稳定中心，一个是鞍点，且两个奇点位置均随施加电压的变化而逐渐靠近，从而得出系统的静态分叉点；同时分析了有阻尼和无阻尼时电压的变化对相轨道的影响，以及阻尼对吸合电压的影响，吸合电压随阻尼的增大而提高，这些研究结果不仅对扭转微镜的设计和应用提供了理论和方法，而且对用电容驱动的微机电系统的设计亦有参考价值。     

温度场下岛-桥结构屈曲薄膜的动力学分析

基于岛-桥结构的柔性电子器件已被用于健康监测和皮肤电子等领域。但是柔性电子器件在工作中极易受工作温度变化等激励产生振动,进而影响器件的灵敏度与可靠性。因此本文研究在温度场作用下岛-桥结构屈曲薄膜的动力学问题。首先,基于Euler-Bernoulli梁理论,建立温度场作用下岛-桥结构屈曲薄膜的动力学控制方程。其次,通过引入新变量,将原动力学方程引入Hamilton体系中,得到相应的Hamilton正则方程。随后,采用辛Runge-Kutta方法求解该Hamilton正则方程,并与经典Runge-Kutta方法对比,数值结果显示了辛算法在求解非线性动力学方程时高精度、高数值稳定性的优势,进一步讨论了温度变化量、预应变、阻尼系数等对屈曲薄膜动力学响应的影响。本文研究为柔性电子器件动力学设计提供了理论参考。     

用机电耦合模型研究转子系统的非平稳过程

建立恰当的转子系统模型,对非平稳过程进行正确完备的描述是进行转子系统设计、实施状态监测和故障诊断的关键。通过考虑驱动电机的影响,建立了异步电机-受驱转子系统的机电耦合数学模型,并利用该模型对机电转子系统由串联补偿电容和负载波动引起的两种机电耦合非平稳过程进行了数值研究。研究结果表明：建立的机电耦合模型可以多方位地研究转子系统在非平衡过程中的机电耦合动力学规律。     

水动力学问题的弱可压缩模型求解方法研究

从弱可压缩水动力学方程出发，采用坐标变换的方法处理自由表面，建立了能够模拟有自由表面流动问题的定常、非定常的三维水动力学模型和对流扩散模型，模型采用浮湍流模型进行封闭，并对模型求解的数值方法进行了研究。     

分数阶导数粘弹性模型的有限元法

在分析分数阶导数三元件模型理论的基础上,把分数阶导数三元件模型引入有限元模型中,推导出具有分数阶导数三元件本构关系的粘弹性结构动力学有限元格式。同时,应用分数阶导数型粘弹性结构动力学方程的数值算法求解了该有限元格式的数值解。并以二维沥青路面结构为例进行了路面动态粘弹性响应分析。算例分析表明,该方法能够正确有效地进行路面动态粘弹性分析。     

含压电层的功能梯度柱壳的自由振动分析

基于三维弹性理论，导出了带有压电层的圆柱形梯度壳的动力学方程以及相应的边界条件，用幂级数展开法得到了求解该圆柱形梯度壳自由振动的三维精确公式．通过实例模型求解了该壳体的自由振动的固有频率；分析了不同电学边界条件对壳体的振动频率的影响。结果可评估各种近似理论解和数值解的正确性。     

有初始缺隐的矩形板非线性屈曲分析

推导出几何非线性和物理非线性且存在初始缺隐的矩形薄板在纵向荷载作用下增量形式的基本微方程，并用伽辽金配点法求解该方程。数值计算结果与已有文献的理论和实验结果相当一致。     

多体系统动力学求解与动画输出的并行处理

在介绍了多体系统动力学的建模方法及数值仿真的原理之后，着重介绍如何利用局域网的技术，在多台计算机之间实现数值求解和动画计算与输出的并行处理。用一个柔性索动力学仿真算例，论证了该技术的可行性，为多体系统动力学仿真的并行处理打下了基础。     

完整约束柔体系统动力学方程的数值求解方法

本文从Ｂａｕｍｇａｒｔｅ的约束久修正法出发，推出了非生动力学微分／代数混合方程约束韧约的自动修正算法，给出了完整约束柔体系统动力学方程的数值求解方法，数值算例说明了该精度和有效性。     

Prediction of buckling characteristics of carbon nanotubes

In this paper, to investigate the buckling characteristics of carbon nanotubes, an equivalent beam model is first constructed. The molecular mechanics potentials in a C–C covalent bond are transformed into the form of equivalent strain energy stored in a three dimensional (3D) virtual beam element connecting two carbon atoms. Then, the equivalent stiffness parameters of the beam element can be estimated from the force field constants of the molecular mechanics theory. To evaluate the buckling loads of multi-walled carbon nanotubes, the effects of van-der Waals forces are further modeled using a newly proposed rod element. Then, the buckling characteristics of nanotubes can be easily obtained using a 3D beam and rod model of the traditional finite element method (FEM). The results of this numerical model are in good agreement with some previous results, such as those obtained from molecular dynamics computations. This method, designated as molecular structural mechanics approach, is thus proved to be an efficient means to predict the buckling characteristics of carbon nanotubes. Moreover, in the case of nanotubes with large length/diameter, the validity of Euler’s beam buckling theory and a shell model with the proper material properties defined from the results of present 3D FEM beam model is investigated to reduce the computational cost. The results of these simple theoretical models are found to agree well with the existing experimental results.     

Compressive buckling analysis of plates in unilateral contact

The unilateral contact buckling behaviour of delaminated plates in a composite member is studied in this paper, where the two-dimensional mechanical problem is simplified to a one-dimensional mathematical model following the assumption of a buckling mode function in terms of the lateral coordinates. After the governing differential equations for the plate parts in the contact and non-contact regions have been solved, the problem is reduced to just two nonlinear algebraic equations allowing the buckling coefficient to be obtained through an iteration method. A simplified method is also deduced based on an assumed zero length of contact. The numerical results of the buckling coefficients based on both methods show good agreement.     

DYNAMIC BUCKLING OF DOUBLE-WALLED CARBON NANOTUBES UNDER STEP AXIAL LOAD

An approximate method is presented in this paper for studying the dynamic buckling of double-walled carbon nanotubes （DWNTs） under step axial load. The analysis is based on the continuum mechanics model, which takes into account the van der Waals interaction between the outer and inner nanotubes. A buckling condition is derived, from which the critical buckling load and associated buckling mode can be determined. As examples, numerical results are worked out for DWNTs under fixed boundary conditions. It is shown that, due to the effect of van der Waals forces, the critical buckling load of a DWNT is enhanced when inserting an inner tube into a single-walled one. The paper indicates that the critical buckling load of DWNTs for dynamic buckling is higher than that for static buckling. The effect of the radii is also examined. In addition, some of the results are compared with the previous ones.     

Interactive buckling in long thin-walled rectangular hollow section struts

An analytical model describing the nonlinear interaction between global and local buckling modes in long thin-walled rectangular hollow section struts under pure compression founded on variational principles is presented. A system of nonlinear differential and integral equations subject to boundary conditions is formulated and solved using numerical continuation techniques. For the first time, the equilibrium behaviour of such struts with different cross-section joint rigidities is highlighted with characteristically unstable interactive buckling paths and a progressive change in the local buckling wavelength. With increasing joint rigidity within the cross-section, the severity of the unstable post-buckling behaviour is shown to be mollified. The results from the analytical model are validated using a nonlinear finite element model developed within the commercial package Abaqus and show excellent comparisons. A simplified method to calculate the local buckling load of the more compressed web undergoing global buckling and the corresponding global mode amplitude at the secondary bifurcation is also developed. Parametric studies on the effect of varying the length and cross-section aspect ratio are also presented that demonstrate the effectiveness of the currently developed models.     

圆柱壳撞水时的弹塑性动力屈曲研究

本文研究圆柱壳在流固中冲击载荷下的弹塑性动力曲问题。建立了液体－气体－固体三相的数学模型。其中结构部分控制方程由弹塑性力学中关于加速度的最小原理获得，本构关系采用增量理论，液体假设为不可压缩，空气层认为最等熵压缩。分别讨论了不同冲击高度时压力变化规律，屈曲对载荷的影响。屈曲沿壳长的分布及发展规律等。     

Flexural-torsional bifurcations of a cantilever beam under potential and circulatory forces II. Post-critical analysis

The post-critical behavior of a cantilever beam with rectangular cross-section, under simultaneous action of conservative and non-conservative loads, is analyzed. An internally constrained Cosserat rod model is adopted to describe the dynamics of the beam in finite displacement regime. The bifurcation equations for simple buckling (divergence), simple flutter (Hopf) and double-zero (Takens-Bogdanova-Arnold) bifurcations are derived by means of the multiple time scales method. Due to the nilpotent eigenvalue at the double-zero critical point, the evaluation of the generalized Keldysh's eigenfunctions is required. Finally, some numerical results are shown and the bifurcation scenario of the beam is discussed.     

SIMULATION OF TRANSVERSE VIBRATION OF SLENDER ROD STRING IN THE VERTICAL CYLINDER UNDER BUCKLING DISPLACEMENT EXCITATION1)

对于铅直圆筒内受交变拉压轴向载荷作用的细长杆柱,当杆柱底端所受到的轴向压力大于杆柱屈曲的临界载荷时,细长杆柱在圆筒内将产生螺旋屈曲,杆柱的屈曲变形将激励杆柱在圆筒内产生横向振动。以细长杆在圆筒内的瞬时屈曲构型作为杆柱横向振动的位移激励,建立了屈曲位移激励下的细长杆在圆筒内横向振动与杆管碰撞规律的仿真模型。采用有限差分法对井深进行离散,Newmark 法对时间进行离散,以恢复系数法建立了细长杆和圆筒的碰撞条件,对细长杆在圆筒内的横向振动数学模型进行了数值仿真。仿真结果表明,细长杆和圆筒内壁的碰撞现象主要发生在细长杆底端受压发生屈曲后,且几乎沿整个圆筒都有发生,从圆筒顶部至底部的碰撞力峰值逐渐增大;而在杆柱底端受拉时碰撞现象很少,碰撞力也较小。     

Plate/shell structure topology optimization of orthotropic material for buckling problem based on independent continuous topological variables

The purpose of the present work is to study the buckling problem with plate/shell topology optimiza-tion of orthotropic material.A model of buckling topology optimization is established based on the independent,con-tinuous, and mapping method, which considers structural mass as objective and buckling critical loads as constraints. Firstly, composite exponential function (CEF) and power function(PF)as filter functions are introduced to recognize the element mass,the element stiffness matrix,and the ele-ment geometric stiffness matrix.The filter functions of the orthotropic material stiffness are deduced. Then these fil-ter functions are put into buckling topology optimization of a differential equation to analyze the design sensitiv-ity.Furthermore,the buckling constraints are approximately expressed as explicit functions with respect to the design vari-ables based on the first-order Taylor expansion.The objective function is standardized based on the second-order Taylor expansion. Therefore,the optimization model is translated into a quadratic program.Finally,the dual sequence quadratic programming(DSQP)algorithm and the global convergence method of moving asymptotes algorithm with two different filter functions(CEF and PF)are applied to solve the opti-mal model.Three numerical results show that DSQP&CEF has the best performance in the view of structural mass and discretion.     

Analysis of dynamical buckling and post buckling for beams by finite segment method

IntroductionThe stability or dynamical buckling of structures is the parametrically exacted vibration ofstructures,and itis a nonlinear vibration problem[1].Dynamical stabilityof beams was studiedusing Galerkin variation method by Russian scholar,and B.B.…     

Buckling of yeast modeled as viscoelastic shells with transverse shearing

Yeast cells can be regarded as micron-sized and liquid-filled cylindrical shells. Owing to the rigid cell walls, yeast cells can bear compressive forces produced during the biotechnological process chain. However, when the compressive forces applied on the yeast go beyond a critical value, mechanical buckling will occur. Since the buckling of the yeast can change the networks in its cellular control, the experimental research of the buckling of the yeast has received considerable attention recently. In this paper, we apply a viscoelastic shell model to study the buckling of the yeast. Meanwhile, the turgor pressure in the yeast due to the internal liquid is taken into account as well. The governing equations are based on the first-order shear deformation theory. The critical axial compressive force in the phase space is obtained by the Laplace transformation, and the Bellman numerical inversion method is then applied to the analytical result to obtain the corresponding numerical results in the physical phase. The concepts of instantaneous critical buckling force, durable critical buckling force, and delay buckling are set up in this paper. And the effects of the transverse shear deformation and the turgor pressure on the buckling phenomena are also given. The numerical results show that the transverse shearing effect will decrease the instantaneous critical buckling force and the durable critical buckling force, while the turgor pressure will increase both of them.