大变形复合材料薄板多体系统动力学建模

本文对大变形复合材料薄板的多体系统动力学建模方法进行研究。基于Kirchhoff假设,法线与中面保持垂直,从格林应变的表达式出发,建立了面内应变和曲率与绝对位置坐标和斜率的关系,在此基础上推导了广义弹性力阵和弹性力阵对广义坐标的导数阵,用绝对节点坐标方法建立了大变形复合材料薄板多体系统的动力学方程,用广义法和和牛顿迭代法求解微分-代数混合方程。对外载荷作用下的复合材料薄板进行数值仿真,通过与ANSYS的仿真结果进行对比,验证了本文建模方法的准确性和快速收敛性。最后,将建模方法应用于复合材料太阳帆板展开机构的数值仿真,分析了不同铺层情况下驱动力和约束力的振动特性。     

多体系统动力学微分-代数方程广义-α 投影法

高效、稳定的微分-代数方程数值求解方法是多体系统动力学领域的关键问题之一。该文针对多体系统动力学指标3微分-代数方程,对目前多体系统动力学中引入的隐式时域逐步积分方法进行了深入研究,提出了适用于一般质量矩阵的广义-α -S法,并结合约束投影方法,构造了广义-α -S投影法。该方法既能较好地保持系统总能量,又能较高程度地同时满足位移约束、速度级约束和加速度级约束,并且在步长较大时可稳定求解,计算效率较高。     

复合柔性结构航天器动力学建模研究

柔性航天器动力学建模的传统方法是采用混合坐标法,针对中心刚体带大型柔性附件类的航天器,这种方法在理论建模和工程应用方面都获得了极大的成功。在中心刚体加柔性附件类航天器柔性动力学研究成果基础上,通过计及柔性体与柔性体连接点间的复合位移变形,利用混合坐标法建立了复合柔性结构航天器动力学模型,其软件系统DASFA 2.0已初步用于工程分析设计。     

考虑附加质量的旋转柔性梁的随机动力学分析

摘要：研究了带有附加质量的中心刚体-柔性悬臂梁系统在参数具有随机性时作大范围运动的动力响应问题。基于假设模态法和Lagrange方程建立了带有附加质量的中心刚体-柔性悬臂梁系统的一次近似耦合随机动力学方程,利用混沌多项式结合高效回归法将其转化为完全隐式纯微分方程,求解方程得到柔性悬臂梁变形位移响应的数字特征。最后,通过数值仿真对物理参数和几何参数具有随机性的系统进行动力特性研究。仿真结果表明：利用随机参数的动力学模型能客观地反映出系统的动力学行为;部分随机参数的分散性对柔性体动力响应的影响不可忽视。     

梁动力学问题重心有理插值配点法

提出数值求解梁动力学问题的高精度重心有理插值配点法。采用重心有理插值张量积形式近似梁在任意时刻及位置挠度,运用配点法获得梁动力学问题控制方程与初边值条件的离散代数方程组。利用微分矩阵与矩阵张量积运算记号,将离散后代数方程组写成简洁矩阵形式。通过置换法施加边界条件及初始条件求解代数方程组,获得梁动力学问题在计算节点处位移值。数值算例表明,重心有理插值配点法具有算式简单、计算节点适应性好、程序实施方便、计算精度高等优点。     

热载荷作用下大变形柔性梁刚柔耦合动力学分析

从非线性应变-位移关系式出发,用虚功原理建立了热载荷作用的柔性梁的热传导方程和旋转刚体-梁系统的刚-柔耦合动力学方程.由于考虑了刚度阵的高次变形项,适用于大变形问题.对温度、弹性变形和刚体运动变量联合求解.研究了热流引起的温度梯度对弹性变形和刚体转动的影响,以及在大变形情况下的几何非线性效应.     

刚柔耦合多体系统动力学模型降阶

提出了基于模态综合法的刚柔耦合多体系统动力学模型降阶方法。该方法用自然坐标法和绝对节点坐标法分别描述刚柔耦合多体系统中的刚体构件和柔性体构件,同时用Craig-Bampton方法对柔性体模型进行减缩。对于刚体构件与柔性体构件之间只存在线性约束的情况,建立了消除线性约束的刚柔耦合多体系统动力学方程。最后,为了验证的该方法的有效性,对刚柔耦合双摆进行了研究。仿真结果表明:适当选择模态就可以在满足计算精度的同时减少计算时间,提高计算效率。     

轮式悬架移动2连杆柔性机械手动力学研究与仿真

对轮式悬架2连杆柔性移动机械手(平面)进行了系统的动力学研究.该轮式移动机械手由带有弹性-阻尼悬架系统的移动载体和柔性机械手所组成,并假定移动载体以恒速通过不规则路面.采用经典瑞利-里兹(Rayleigh-Ritz)法和浮动坐标法描述机械手弹性变形与参考运动间的动力学耦合问题,综合利用拉格朗日原理和牛顿-欧拉方程并在笛卡尔坐标系下 ,以矩阵、矢量简洁的形式构建了该移动柔性机械手系统的完整动力学模型.最后采用数值的方法给出了该动力学模型正解的仿真结果.通过与刚体模型、刚柔混合模型仿真结果的比较,证实了该柔体系统存在动力学耦合现象.     

关于微分-代数系统的Runge-Kutta迭代方法的研究

为了求解非自制指标-1的微分-代数系统,我们研究基于Runge-Kutta方法的动力学迭代过程,得到相关的非线性微分-代数方程的收敛理论,这类迭代过程具有一般性和灵活性,且沿着时间域网格点可以选取不同的插值函数.     

无网格法在中心刚体-旋转柔性梁系统动力学分析中的应用EI北大核心CSCD

将无网格点插值法、径向基点插值法、光滑节点插值法用于中心刚体-旋转柔性梁的动力学分析。基于浮动坐标系方法,考虑梁的纵向拉伸变形和横向弯曲变形,并计入横向弯曲变形引起的纵向缩短,即非线性耦合项,运用第二类Lagrange方程推导得到作大范围运动的中心刚体-旋转柔性梁系统的动力学方程。将无网格法的仿真结果与有限元法和假设模态法进行比较分析,表明其作为一种柔性体离散方法在中心刚体-旋转柔性梁的刚柔耦合多体系统动力学的研究中具有可推广性。     

Direct differentiation method for sensitivity analysis based on transfer matrix method for multibody systems

On one hand, the new version of transfer matrix method for multibody systems (NV‐MSTMM), has been proposed by formulating transfer equations of elements in acceleration level instead of position level as in the original discrete time transfer matrix method of multibody systems to study multibody system dynamics. This new formulation avoids local linearization and allows using any integration algorithms. On the other hand, sensitivity analysis is an important way to improve the optimization efficiency of multibody system dynamics. In this paper, a totally novel direct differentiation method based on NV‐MSTMM for sensitivity analysis of multibody systems is developed. Based on direct differentiation method, sensitivity analysis matrix for each kind of element is established. By assembling transfer matrices and sensitivity analysis matrices based on differentiation law of multiplication, the sensitivity analysis equation of overall transfer equation is deduced. The computing procedure of the proposed method is also presented. All these improvements as well as three numerical examples show that the direct differentiation method based on NV‐MSTMM is suitable for optimizing the dynamic sensitivity in multi–rigid‐body systems.     

广义多体系统动力学的速度变换矩阵综合法

文中提出了广义多体系统和速度变换矩阵的概念,提出了一种新的加速度变换关系,以带不定乘子的拉格朗日方程为基础推导得到了求解复杂系统动力学问题的一种新方法,即广义多体系统的速度变换矩阵综合法。利用该方法,可根据无耦合广义体的动力学参数和系统的速度变换矩阵直接获得广义多体系统的动力学方程,其中不含拉格朗日不定乘子和约束反力,且方程中逆矩阵求解的维数等于系统的自由度数,因而有利于提高计算效率。该方法主要面向计算机实现程式化的算法,系统的动力方程可以由计算机自动完成运算,从而避免了繁琐的解析推导工作。     

A novel nonsmooth dynamics method for multibody systems with friction and impact based on the symplectic discrete format

As multibody systems often involve unilateral constraints, nonsmooth phenomena, such as impacts and friction, are common in engineering. Therefore, a valid nonsmooth dynamics method is highly important for multibody systems. An accuracy representation of multibody systems is an important performance indicator of numerical algorithms, and the energy balance can be used efficiently evaluate the performance of nonsmooth dynamics methods. In this article, differential algebraic equations (DAEs) of a multibody system are constructed using the D'Alembert's principle, and a novel nonsmooth dynamics method based on symplectic discrete format is proposed. The symplectic discrete format can maintain the energy conservation of a conservative system; this property is expected to extend to nonconservative systems with nonsmooth phenomena in this article. To evaluate the properties of the proposed method, several numerical examples are considered, and the results of the proposed method are compared with those of Moreau's midpoint rule. The results demonstrate that the solutions obtained using the proposed method, which is based on the symplectic discrete format, can realize a higher solution accuracy and lower numerical energy dissipation, even under a large time step.     

An efficient high-precision recursive dynamic algorithm for closed-loop multibody systems

As most closed-loop multibody systems do not have independent generalized coordinates, their dynamic equations are differential/algebraic equations (DAEs). In order to accurately solve DAEs, a usual method is using generalized α-class numerical methods to convert DAEs into difference equations by differential discretization and solve them by the Newton iteration method. However, the complexity of this method is O(n²) or more in each iteration, since it requires calculating the complex Jacobian matrix. Therefore, how to improve computational efficiency is an urgent problem. In this paper, we modify this method to make it more efficient. The first change is in the phase of building dynamic equations. We use the spatial vector note and the recursive method to establish dynamic equations (DAEs) of closed-loop multibody systems, which makes the Jacobian matrix have a special sparse structure. The second change is in the phase of solving difference equations. On the basis of the topology information of the system, we simplify this Jacobian matrix by proper matrix processing and solve the difference equations recursively. After these changes, the algorithm complexity can reach O(n) in each iteration. The algorithm proposed in this paper is not only accurate, which can control well the position/velocity constraint errors, but also efficient. It is suitable for chain systems, tree systems, and closed-loop systems.     

微分对接条件对次谐共振影响的研究

根据模态综合法并借助二阶非完整约束系统的Routh型方程,建立了考虑全部代数和微分对接条件下,非线性转子—支承系统的运动微分方程。然后采用一种新的等效线性化技术,求解系统的次谐共振,大大简化了分析过程并提高了计算效率。通过对对接条件之作用的定量研究,结果表明:微分对接条件对次谐共振影响较大,而由不独立的微分对接条件转化得来的代数对接条件,对次谐共振影响很小;考虑微分对接条件,能明显提高方法的收敛速度。因此,应当考虑微分对接条件,尽管这会增加推导工作量和综合后矩阵的阶数,但并不影响计算效率。     

Interval uncertain method for multibody mechanical systems using Chebyshev inclusion functions

This study proposes a new uncertain analysis method for multibody dynamics of mechanical systems based on Chebyshev inclusion functions The interval model accounts for the uncertainties in multibody mechanical systems comprising uncertain‐but‐bounded parameters, which only requires lower and upper bounds of uncertain parameters, without having to know probability distributions. A Chebyshev inclusion function based on the truncated Chebyshev series, rather than the Taylor inclusion function, is proposed to achieve sharper and tighter bounds for meaningful solutions of interval functions, to effectively handle the overestimation caused by the wrapping effect, intrinsic to interval computations. The Mehler integral is used to evaluate the coefficients of Chebyshev polynomials in the numerical implementation. The multibody dynamics of mechanical systems are governed by index‐3 differential algebraic equations (DAEs), including a combination of differential equations and algebraic equations, responsible for the dynamics of the system subject to certain constraints. The proposed interval method with Chebyshev inclusion functions is applied to solve the DAEs in association with appropriate numerical solvers. This study employs HHT‐I3 as the numerical solver to transform the DAEs into a series of nonlinear algebraic equations at each integration time step, which are solved further by using the Newton–Raphson iterative method at the current time step. Two typical multibody dynamic systems with interval parameters, the slider crank and double pendulum mechanisms, are employed to demonstrate the effectiveness of the proposed methodology. The results show that the proposed methodology can supply sufficient numerical accuracy with a reasonable computational cost and is able to effectively handle the wrapping effect, as cosine functions are incorporated to sharpen the range of non‐monotonic interval functions. Copyright © 2013 John Wiley & Sons, Ltd.     

A Lagrangian approach to the non-causal inverse dynamics of flexible multibody systems: The three-dimensional case

This paper addresses the problem of end-point trajectory tracking in flexible multibody systems through the use of inverse dynamics. A global Lagrangian approach is employed in formulating the system equations of motion, and an iterative procedure is proposed to achieve end-point trajectory tracking in three-dimensional, flexible multibody systems. Each iteration involves firstly, a recursive inverse kinematics procedure wherein elastic displacements are determined in terms of the rigid body co-ordinates and Lagrange multipliers, secondly, an explicit computation of the inverse dynamic joint actuation, and thirdly, a non-recursive forward dynamic analysis wherein generalized co-ordinates and Lagrange multipliers are determined in terms of the joint actuation and desired end-point co-ordinates. In contrast with the recursive methods previously proposed, this new method is the most general since it is suitable for both open-chain and closed-chain configurations of three-dimensional multibody systems. The algorithm yields stable, non-casual actuating joint torques and associated Lagrange multipliers that account for the constraint forces between flexible multibody components.     

基于二阶矩阵微分方程的机械振动系统线性二次型调节器设计

本文讨论机械振动系统线性二次型状态调节器（LQR）问题，直接针对系统二阶运动微分方程，性能指标为一个依赖于二阶导数的泛函。由欧拉-拉格朗日方程得出一个系统矩阵增广的二阶线性微分方程，指出该方程稳定的特征对就是最优控制振动系统闭环特征对，并给出求解最优控制状态反馈矩阵的方法，另外，由本文方法还可得出基于速度和加速度反馈的最优控制反馈矩阵。这里不涉及求解代数矩阵Riccati方程。     

受控线性多体系统传递矩阵法

受控线性多体系统的稳态运动是多体系统动力学的重要研究内容之一。本文以简单的受控多体系统为例，建立了受控线性多体系统传递矩阵法，能方便快捷地求解受控线性多体系统的稳态运动。建立了控制力作用下集中质量和弹簧阻尼铰的扩展传递矩阵和扩展传递方程。将牛顿法和受控线性多体系统传递矩阵法的计算结果进行了比较。实例表明，受控线性多体系统传递矩阵法不仅能用于受控线性多体系统的动力学分析，而且完整地保持了传递矩阵法的所有优点。     

色噪声与确定性谐波联合激励下Bouc-Wen动力系统响应的统计线性化方法

提出了一种用于求解色噪声和确定性谐波联合作用下单自由度Bouc?Wen系统响应的统计线性化方法.基于系统响应可分解为确定性谐波和零均值随机分量之和的假定,将原滞回运动方程等效地化为两组耦合的且分别以确定性和随机动力响应为未知量的非线性微分方程.利用谐波平衡法求解确定性运动方程,利用统计线性化方法求解色噪声激励下的随机运...