考虑剪切效应的旋转FGM楔形梁刚柔耦合动力学建模与仿真

本文对一类中心刚体-柔性梁系统在大范围转动下的刚柔耦合动力学问题进行了研究. 柔性梁为功能梯度材料(functionally graded materials, FGM)楔形变截面梁,材料体积分数在梁轴向呈幂律分布变化. 以弧长坐标来描述柔性FGM梁的几何位移关系,分别使用倾角和拉伸应变变量描述柔性梁的横向弯曲和纵向拉伸变形,并计及剪切效应. 采用假设模态法离散变形场,运用第二类拉格朗日方程进行方程推导,得到系统考虑剪切效应的刚柔耦合动力学模型. 基于全新的刚柔耦合动力学建模理论,研究不同轴向材料梯度分布的FGM楔形梁,通过数值仿真计算,分析讨论不同的转速、梯度分布规律以及变截面参数对系统动力学特性的影响. 结果表明,剪切效应对大高跨比的FGM楔形梁的变形影响较为明显,不容忽略;材料梯度分布规律和截面参数的选取均会对旋转FGM楔形梁的动力学响应和频率产生较大影响. 本文提出的考虑剪切效应的倾角刚柔耦合动力学模型是对以往非剪切模型的进一步完善,可应用于工程中的 Timoshenko梁结构的动力学问题求解.      

A shearable and thickness stretchable finite strain beam model for soft structures

Soft materials and structures have recently attracted lots of research interests as they provide paramount potential applications in diverse fields including soft robotics, wearable devices, stretchable electronics and biomedical engineering. In a previous work, an Euler–Bernoulli finite strain beam model with thickness stretching effect was proposed for soft thin structures subject to stiff constraint in the width direction. By extending that model to account for the transverse shear effect, a Timoshenko-type finite strain beam model within the plane-strain context is developed in the present work. With some kinematic hypotheses, the finite deformation of the beam is analyzed, constitutive equations are deduced from the theory of finite elasticity, and by employing the standard variational method, the equilibrium equations and associated boundary conditions are derived. In the limit of infinitesimal strain, the new model degenerates to the classical extensible and shearable elastica model. The corresponding incremental equilibrium equations and associated boundary conditions are also obtained. Based on the new beam model, analytical solutions are given for simple deformation modes, including uniaxial tension, simple shear, pure bending, and buckling under an axial load. Furthermore, numerical solution procedures and results are presented for cantilevered beams and simply supported beams with immovable ends. The results are also compared with the previously developed finite strain Euler–Bernoulli beam model to demonstrate the significance of transverse shear effect for soft beams with a small length-to-thickness ratio. The developed beam model will contribute to the design and analysis of soft robots and soft devices.     

Thermal effect on the deformation of a flexible beam with large kinematical driven overall motions

The research investigates the transient longitudinal and transverse deformation of a planar flexible beam with large overall motions in a temperature field. With the increase of temperature, longitudinal deformation is caused by the thermal expansion in axial direction. Due to the coupling of longitudinal and transverse deformation, the transverse deformation is induced, which is significant in cases where temperature increases rapidly in a very short period. Furthermore, the transverse temperature gradient, which is caused by the temperature variation in the transverse direction, may lead to transverse deformation. Considering the thermal strain, equations of motion of a flexible beam with arbitrary large overall motion are derived based on virtual work principle. The high order terms of the strain tensor are taken into account, such that the geometric nonlinear deformation terms are included in the dynamic equations. Simulation results of a rotating beam are shown to reveal the thermal effect and nonlinear effect on the dynamic performance of the beam.     

应变能中耦合项对旋转悬臂梁振动频率的影响

大范围运动悬臂梁的动力学建模问题对动力学特性分析及控制系统设计具有极其重要的作用. 当前研究多采用一次近似模型,其忽略了由轴向和横向变形所产生的应变能中的耦合项,然而这些项对动力学特性会产生影响. 通过讨论应变能的选取方式,计入了应变能中的耦合项;利用哈密尔顿原理建立结构的耦合振动模型;再借助瑞利—里兹法,以无大范围运动时的振型函数作为基本解组,得到了结构振动广义特征方程并求解. 通过数值算例对比分析,指出考虑应变能耦合项得到的频率与不考虑应变能耦合项得到的频率存在明显差别.      

An image encryption algorithm based on compound homogeneous hyper-chaotic system

A new dynamic model of a rotating flexible beam with a concentrated mass located in arbitrary position is derived based on the absolute nodal coordinate formulation, and its modal characteristics are investigated in this paper. To consider the concentrated mass at an arbitrary location of the beam, a Dirac’s delta function is used to express the mass per unit length of the beam. Based on the proposed dynamic model, the frequency analysis is performed. The nonlinear equation is transformed into the linear one via employing the linear perturbation analysis method. The stiffness matrix of static equilibrium of the system under the deformed condition is obtained, in which the effect of coupling between the longitudinal deformation and transversal deformation is included. This means even if only the chordwise bending equation is solved, the longitudinal vibration effect can be still considered. As we know, once the longitudinal deformation is large, it will significantly affect the chordwise bending vibration. So the proposed model in this paper is more accurate than the traditional dynamic models which are usually lack of the coupling terms between the longitudinal deformation and transversal deformation. In fact, the traditional dynamic models for the chordwise vibration analysis in the existing literature are usually linear due to neglecting the coupling terms, and consequently, they are only suitable for the modal characteristic analysis of a beam under small deformations. In order to get some general conclusions of the natural frequencies and mode shapes, the equation which governs the chordwise bending vibration of the rotating beam is transformed into a dimensionless form. The dynamic model presented in this paper is nonlinear and can be conveniently used to analyze the modal characteristics of a rotating flexible beam with large deformations. To demonstrate the power of the new dynamic model presented in this paper, the dynamic simulations involving the comparisons between the different frequencies obtained using the model proposed in this paper and the models in the existing literature and the investigating in frequency veering and mode shift phenomena are given. The simulation results show that the angular velocity of the flexible beam will give rise to the phenomena of the natural frequency loci veering and the associated mode shift which is verified in the previous studies. In addition, the phenomena of the natural frequency loci veering rather than crossing can be observed due to the changing of the magnitude of the concentrated mass or of the location of the concentrated mass which are found for the first time. Furthermore, there is an interesting phenomenon that the natural frequency loci will veer more than once due to different types of mode coupling between the bending and stretching vibrations of the rotating beam. At the same time, the mode shift phenomenon will occur correspondingly. Additionally, the characteristics of the vibration nodes are also investigated in this paper.     

弹性曲梁在机械和热载荷共同作用下的几何非线性模型及其数值解

基于Euler-Bernoulli梁的几何非线性理论,建立了弹性曲梁在任意分布机械载荷和热载荷共同作用下的几何非线性静平衡控制方程。该模型不仅计及了轴线伸长,同时也精确地考虑了梁的初始曲率对变形的影响以及轴向变形与弯曲变形之间的相互耦合效应。应用打靶法数值求解了半圆形曲梁在横向均匀升温作用下的非线性弯曲问题,数值比较了轴向伸长对曲梁变形的影响。     

带集中质量的旋转柔性曲梁动力学特性分析

本文对带集中质量的平面内旋转柔性曲梁动力学特性进行了研究.基于绝对节点坐标法推导出曲梁单元,其中该曲梁单元采用Green-Lagrangian应变,并根据曲梁变形前后的曲率变化和曲率的精确表达式计算了曲梁单元弹性力所作的虚功.通过虚功原理,利用$\delta$函数和中心刚体与悬臂曲梁之间的固支边界条件,建立了带集中质量的旋转柔性曲梁非线性动力学模型.基于该模型,本文仿真计算了悬臂曲梁的纯弯曲问题和带有刚柔耦合效应的旋转柔性曲梁动力学响应问题,以此分别讨论了所提出曲梁单元的收敛性和动力学模型的正确性.进一步应用D'Alembert原理,将旋转曲梁等效为带离心力的无旋转曲梁,通过线性摄动处理得到系统的特征方程,以此分别研究了旋转角速度、初始曲率和集中质量对曲梁动力学特性的影响.最后重点分析了旋转曲梁的频率转向和振型切换问题,并阐述了两者之间的相互关系.研究结果表明:随着旋转角速度的增大,曲梁的频率特性与直梁的频率特性相近,以及曲梁拉伸变形占主导的模态振型会提前.      

考虑尺度效应的夹层楔形多孔梁动力学分析

基于修正偶应力理论, 研究了具有大范围旋转中心刚体-功能梯度夹层Euler-Bernoulli楔形多孔柔性微梁系统的动力学特性.楔形梁是中间层为不完全功能梯度层, 两表层为均质材料的功能梯度夹层结构, 它可以减小传统夹层结构由于层与层之间材料属性的不同导致脱粘类型损伤的影响.采用假设模态法描述变形, 考虑具有捕捉动力刚化效应的非线性耦合项, 计及von Kármán几何非线性应变, 运用第二类Lagrange方程, 导出了适用于较大变形的高次刚柔耦合动力学方程.对在平面内做大范围运动的中心刚体-功能梯度夹层Euler-Bernoulli楔形多孔微梁的动力学特性进行了详细研究.研究表明: 功能梯度夹层楔形梁表层结构高度、旋转角速度、功能梯度幂指数、尺度参数、孔隙度以及各层结构的体积分数对系统的动力学特性都有很大的影响; 功能梯度夹层楔形梁综合了功能梯度直梁和楔形梁的特性, 其相对于功能梯度直梁的固有频率增大, 同时使得孔隙度对结构固有频率变化趋势的影响不再与功能梯度直梁相同; 由于柔性梁变形能中具有横向与轴向的耦合势能, 系统在稳态下的平衡位置发生了迁移现象; 系统随着尺度参数的变化发生了频率转向与振型转换.      

Geometric nonlinear dynamic analysis of curved beams using curved beam element

Instead of using the previous straight beam element to approximate the curved beam,in this paper,a curvilinear coordinate is employed to describe the deformations,and a new curved beam element is proposed to model the curved beam.Based on exact nonlinear strain-displacement relation,virtual work principle is used to derive dynamic equations for a rotating curved beam,with the effects of axial extensibility,shear deformation and rotary inertia taken into account.The constant matrices are solved numerically utilizing the Gauss quadrature integration method.Newmark and Newton-Raphson iteration methods are adopted to solve the differential equations of the rigid-flexible coupling system.The present results are compared with those obtained by commercial programs to validate the present finite method.In order to further illustrate the convergence and efficiency characteristics of the present modeling and computation formulation,comparison of the results of the present formulation with those of the ADAMS software are made.Furthermore,the present results obtained from linear formulation are compared with those from nonlinear formulation,and the special dynamic characteristics of the curved beam are concluded by comparison with those of the straight beam.     

Dynamic analysis of a rotating rigid-flexible coupled smart structure with large deformations

Based on Hamilton's principle,a new kind of fully coupled nonlinear dynamic model for a rotating rigid-flexible smart structure with a tip mass is proposed.The geometrically nonlinear effects of the axial,transverse displacement and rotation angle are considered by means of the first-order approximation coupling(FOAC)model theory, in which large deformations and the centrifugal stiffening effects are considered.Three kinds of systems are established respectively,which are a structure without piezoelectric layer,with piezoelectric layer in open circuit and closed circuit.Several simulations based on simplified models are presented to show the differences in characteristics between structures with and without the tip mass,between smart beams in closed and open circuit, and between the centrifugal effects in high speed rotating state or not.The last simulation calculates the dynamic response of the structure subjected to external electrical loading.     

基于初应力法的作大范围运动柔性梁动力学建模理论研究

研究了初应力法的作大范围运动柔性梁的建模理论.根据连续介质理论,考虑应变-位移中的非线性项,用一致质量有限元法对柔性梁进行离散,基于Jourdain速度变分原理导出定轴转动下大范围运动为自由的柔性梁刚-柔耦合动力学方程.从其刚柔耦合动力学方程出发,考虑在大范围运动已知情况下的结构动力学方程.通过引入准静态概念,把其结构动力学方程转化为准静态方程.对纵向和横向变形节点坐标进行坐标分离,解出与纵向变形相关的准静态方程,得到准静态时的纵向应力表达式,从而获得附加刚度项.并对此非惯性系下作大范围运动柔性梁的结构动力学方程进行数值仿真,对零次近似模型、一次近似模型、初应力法动力学模型的仿真结果进行分析,揭示三种模型的动力学性质的差异.     

考虑曲率纵向变形效应的大变形柔性梁刚柔耦合动力学建模与仿真

对在平面内做大范围转动的中心刚体柔性梁系统的动力学进行了研究,建立了考虑大变形效应的系统刚柔耦合动力学模型,并进行了动力学仿真.该动力学模型不但考虑了柔性梁横向弯曲变形和纵向变形(包含轴向拉伸变形和横向弯曲变形而引起的纵向缩短项),还考虑了纵向变形对曲率的影响,称为曲率纵向变形效应.在以往的研究中,柔性梁的横向弯曲变形能往往直接使用柔性梁横向弯曲变形来表达,并没有考虑纵向变形的影响.为了考虑柔性梁纵向变形对横向弯曲变形能的影响,在浮动坐标系下使用柔性梁参数方程形式的精确曲率公式来计算柔性梁的弯曲变形能.在此基础上建立了基于浮动坐标系的考虑曲率纵向变形效应的刚耦合动力学模型.论文给出了数值仿真算例,验证了本文所建的动力学模型既能适用于柔性梁的小变形问题,又能适用于大变形问题,且较现有高次刚柔耦合动力学模型更加适用于大变形问题的处理.论文还通过与能处理柔性梁大变形问题的绝对节点坐标法的比较,验证了模型的正确性.      

NONLINEAR DYNAMIC ANALYSIS OF FLEXIBLE MULTIBODY SYSTEM

The nonlinear dynamic equations of a multibody system composed of ?exible beams are derived by using the Lagrange multiplier method. The nonlinear Euler beam theory with inclusion of axial deformation e?ect is employed and its deformation ?eld is …     

非保守功能梯度弹性组合曲梁的精确模型及其数值解

基于直法线假设,采用可伸长梁的几何非线性理论,建立了功能梯度材料弹性组合曲梁受切线均布随从力作用下的静态大变形数学模型。该模型不仅计及了轴线伸长,同时也精确地考虑了梁的初始曲率对变形的影响以及轴向变形与弯曲变形之间的耦合效应。用打靶法数值求解了由金属和陶瓷两相材料所构成的一种FGM组合曲梁在沿轴线均布切向随动载荷作用下的非线性平面弯曲问题,给出了不同梯度指标下FGM弹性曲梁随载荷参数大范围变化的平衡路径,并与金属和陶瓷两种单相材料曲梁的相应特性进行了比较。     

旋转中心刚体-FGM梁刚柔热耦合动力学特性研究

对旋转中心刚体-功能梯度材料(functionally graded material,FGM)梁刚柔热耦合动力学特性进行研究.FGM梁为物理性能参数沿厚度方向呈幂律分布的欧拉伯努利梁.考虑柔性梁的横向弯曲变形和轴向拉伸变形, 并计入横向弯曲变形引起的纵向缩短,即非线性耦合变形量.考虑变截面空心梁在外部高温、内冷通道冷却情况下的热力耦合对系统动力学特性的影响,求解得到FGM梁沿厚度方向分布的温度场, 进而在本构关系中计入热应变.采用假设模态法对柔性梁变形场进行离散,运用第二类拉格朗日方程推导得到系统的刚柔热耦合动力学方程,并编制动力学仿真软件, 然后通过仿真算例对系统的动力学问题进行研究.结果表明:不同截面梁动力学响应差异较大, 因此需对实际系统合理建模;大范围运动已知时, 考虑热冲击载荷的FGM梁将有效抑制横向弯曲变形,而大范围运动恒定时随热冲击的叠加会出现高频振荡; 大范围运动未知时,外力矩和热冲击载荷相互作用产生热力耦合效应, 导致系统呈现高频振荡,同时与中心刚体大范围旋转运动产生刚柔热耦合效应.      

基于径向基点插值法的旋转Mindlin板高次刚柔耦合动力学模型

将无网格径向基点插值法(radial point interpolation method,RPIM)用于中心刚体?旋转柔性板的动力学分析.基于浮动坐标系方法和一阶剪切变形理论即Mindlin板理论,考虑剪切变形的影响,并计入板面内变形的非线性耦合变形项,采用径向基点插值法描述板的变形场,保留动能中有关非线性耦合变形项...     

柔性梁刚柔耦合碰撞动力学及变形传播研究

运用柔性多体系统刚柔耦合动力学理论,研究了作大范围回转运动柔性梁的碰撞动力学问题.考虑柔性梁的横向变形,以及横向变形引起的纵向缩短项即非线性耦合变形项.采用基于Hertz接触理论及非线性阻尼理论的非线性弹簧阻尼模型来求解碰撞过程中产生的碰撞力,运用第二类拉格朗日方程建立了系统的刚柔耦合碰撞动力学方程.编制仿真软件进行动力学仿真计算,得到了碰撞力和系统动力学响应,对比分析了不同动力学模型对系统动力学响应的影响.同时研究了碰撞导致的柔性梁横向变形传播的波动特性.     

Constrained relative motions in rotational mechanics

The dynamical effects of imposing constraints on the relative motions of component parts in a rotating mechanical system or structure are explored. It is noted that various simplifying assumptions in modeling the dynamics of elastic beams imply strain constraints, i.e., that the structure being modeled is rigid in certain directions. In a number of cases, such assumptions predict features in both the equilibrium and dynamic behavior which are qualitatively different from what is seen if the assumptions are relaxed. It is argued that many pitfalls may be avoided by adopting so-called geometrically exact models, and examples from the recent literature are cited to demonstrate the consequences of not doing this. These remarks are brought into focus by a detailed discussion of the nonlinear, nonlocal model of a shear-free, inextensible beam attached to a rotating rigid body. Here it is shown that linearization of the equations of motion about certain relative equilibrium configurations leads to a partial differential equation. Such spatially localized models are not obtained in general, however, and this leaves open questions regarding the way in which the geometry of a complex structure influences computational requirements and the possibility of exploiting parallelism in performing simulations. A general treatment of linearization about implicit solutions to equilibrium equations is presented and it is shown that this approach avoids unintended imposition of constraints on relative motions in the models. Finally, the example of a rotating kinematic chain shows how constraining the relative motions in a rotating mechanical system may destabilize uniformly rotating states.     

刚-柔耦合动力学系统的建模理论研究

刚－柔耦合动力学系统的传统的混合坐标方法是零次近似方法,在建模过程中,直接套用的结构动力学的小变形假设,忽略了变形位移的高次耦合变形量．本文对柔性梁建立较零次近似更精确的高次耦合动力学模型,从连续介质力学理论出发,在变形位移中,计及横向位移引起的轴向缩短,导出变形位移的二次耦合量．用一致质量有限元方法对梁进行离散,基于Jourdain速度变分原理导出大范围运动为自由的柔性梁的刚－柔耦合动力学方程．计算了柔性重力摆的角速度和摆端点的横向变形,揭示零次近似模型和耦合模型的刚－柔耦合动力学性质的根本差异．     

A robust hyper-elastic beam model under bi-axial normal-shear loadings

In this paper, a simple and robust constitutive model is proposed to simulate mechanical behaviors of hyper-elastic materials under bi-axial normal-shear loadings in the finite strain regime. The Mooney–Rivlin strain energy function is adopted to develop a two-dimensional (2D) normal-shear constitutive model within the framework of continuum mechanics. A motion field is first proposed for combined normal and shear deformations. The deformation gradient of the proposed field is calculated and then substituted into right Cauchy–Green deformation tensor. Constitutive equations are then derived for normal and shear deformations. They are two explicit coupled equations with high-level polynomial non-linearity. In order to examine capabilities of the developed hyper-elastic model, uniaxial tensile responses and non-linear stability behaviors of moderately thick straight and curved beams undergoing normal axial and transverse shear deformations are simulated and compared with experiments. Fused deposition modeling technique as a 3D printing technology is implemented to fabricate hyper-elastic beam structures from soft poly-lactic acid filaments. The printed specimens are tested under tensile/compressive in-plane and compressive out-of-plane forces. A finite element formulation along with the Newton–Raphson and Riks techniques is also developed to trace non-linear equilibrium path of beam structures in large defamation regimes. It is shown that the model is capable of predicting non-linear equilibrium characteristics of hyper-elastic straight and curved beams. It is found that the modeling of shear deformation and finite strain is essential toward an accurate prediction of the non-linear equilibrium responses of moderately thick hyper-elastic beams. Due to simplicity and accuracy, the model can serve in the future studies dealing with the analysis of hyper-elastic structures in which two normal and shear stress components are dominant.