长条形基坑带顶部支撑排桩支护体系空间变形研究

提出了一种建立具有固定双面约束多点摩擦的多体系统动力学方程的方法. 用笛卡尔坐标阵 描述系统的位形,根据局部方法的递推关系建立系统的约束方程,应用第一类Lagrange方程 建立该系统的动力学方程,使得具有摩擦的约束面的法向力与Lagrange乘子一一对应,便于 摩擦力的分析与计算,并用矩阵形式给出了摩擦力的广义力的一般表达式. 应用增广法将微 分-代数方程组转化为常微分方程组,并用分块矩阵的形式给出,以便于方程的编程与计算. 给出了一种改进的试算法,可提高计算效率. 最后给出了一个算例,应用试算法和RK法对算 例进行了数值仿真.     

双面约束多点摩擦多体系统的建模和数值方法

提出了一种建立具有固定双面约束多点摩擦的多体系统动力学方程的方法. 用笛卡尔坐标阵描述系统的位形,根据局部方法的递推关系建立系统的约束方程,应用第一类Lagrange方程建立该系统的动力学方程,使得具有摩擦的约束面的法向力与Lagrange乘子一一对应,便于摩擦力的分析与计算,并用矩阵形式给出了摩擦力的广义力的一般表达式. 应用增广法将微分-代数方程组转化为常微分方程组,并用分块矩阵的形式给出,以便于方程的编程与计算.给出了一种改进的试算法,可提高计算效率. 最后给出了一个算例,应用试算法和RK法对算例进行了数值仿真.      

二维系统传递矩阵法

以薄板横向振动的振动特性问题为例提出了二维系统传递矩阵法. 从质点和无质量梁的传递 矩阵出发，建立了按列排列的薄板子结构的传递矩阵. 用二维系统传递矩阵法获得了整个板 的总体传递方程，从而可得到板在任意一种边界条件下的特征方程. 数值求解了两种情况下 薄板横向振动的固有振动特性. 计算结果表明，用该方法可用于研究类似薄板的二维系统的 动力学问题，且无需建立系统的总体动力学方程.     

挠性多体系统动力学

本文以递推形式给出了挠性多体系统的动力学方程。此方程具有如下特点:(1)对整个系统而言,方程数目多,但递推的每一步只需要解低维的方程组,从而避免了高维矩阵的求逆,减少了计算量。(2)对不同类型的问题,方程形式具有很高的相似性,便于程序编制,各个需要求逆的矩阵都是对称正定的,这样减少了计算量和计算机内存单元。(3)这种方程形式处理非树形和带有约束条件的系统很方便。     

漂浮基空间机器人系统基于增广变量思想的改进变结构滑模控制

讨论载体位置与姿态均不受控制的漂浮基空间机器人系统的控制问题.首先导出了空间机器人欠驱动形式的系统动力学方程.之后借助于增广变量法,证明可以将上述系统动力学方程及系统增广广义Jacobi矩阵分别表示为一组适当选择的组合惯性参数的线性函数.以此为基础,根据具有较强鲁棒性的变结构滑模控制理论,设计了一种空间机器人惯性空间期望轨迹跟踪的改进变结构滑模控制方案.与传统变结构滑模控制相比,所提控制方案通过一次离线预估控制律中相关矩阵元素的上下限,从而避免了实时控制过程中重复计算系统动力学方程中科氏力、离心力项的麻烦,因此有效减少了计算量,更适用于机载计算机运算能力有限的空间机器人控制系统实时应用.仿真运算,证实了方法的有效性.     

多体系统传递矩阵法研究进展

作为一种多体系统动力学新方法, 多体系统传递矩阵法由于其无需系统总体动力学方程和快速计算的特点, 已被广泛用于各种多管火箭、自行火炮、舰炮等复杂大型机械系统动力学分析与设计. 本文介绍了该方法的研究进展, 包括: 线性多体系统传递矩阵法、多体系统离散时间传递矩阵法、二维系统传递矩阵法、受控多体系统传递矩阵法、多体系统传递矩阵法和通常动力学方法的混合方法等, 给出了该方法解决自行火炮、多管火箭武器多体系统动力学的重大工程应用实例.      

多体系统动力学方程违约修正的数值计算方法

多体系统动力学方程为微分代数方程,一般将其转化成常微分方程组进行数值计算,在数值积分的过程中约束方程的违约会逐渐增大.本文对具有完整、定常约束的多体系统,在修改的带乘子Lagrange正则形式的方程的基础上,根据Baumgarte提出的违约修正的方法,给出了一种多体系统微分代数方程违约修正法和系统的动力学方程的矩阵表达式.通过对曲柄-滑块机构的数值仿真,计算结果表明本文给出的方法在计算精度和计算效率上好于Baumgarte提出的两种违约修正的方法.     

多体系统动力学Kane方法的改进

基于Kane方法, 针对约束多体系统, 建立了一种新型的自动组集系统动力学方程的方法. 首先提出偏速度矩阵和偏角速度矩阵的概念, 将各体对系统广义惯性力的贡献用简洁、统一的数学形式表达. 然后引入各个运动学变量的递推关系以提高建模效率.最后对新型的Kane方程进行扩展, 用于处理多体系统中的运动约束. 该算法适用于任意多体系统, 建立的动力学模型不含待定乘子, 维数与系统广义速率相同, 利于控制系统设计.对带有闭环约束的空间多机械臂系统的数值仿真验证了方法的正确性.      

弹性动力学高阶核无关快速多极边界元法

基于核无关的快速多极方法, 发展了一种弹性动力学问题的快速、高精度边界元分析方法. 采用基于二次曲面单元的Nystr?m 离散, 将边界积分方程转化为求和形式, 可以方便地进行加速计算;由于采用二次元, 边界元分析精度很高. 将一种新型快速多极方法用于Nystr?m 边界元法的加速计算, 该方法的数值实现简便、不依赖于积分方程基本解的表达式, 因此通用性很好;该方法还具有最优的计算量和存储量、精度高且可以控制. 结合Nystr?m 边界元系数矩阵和快速多极方法转换矩阵的特点, 提出一种大幅度降低边界元内存消耗的策略. 数值结果表明, 该方法无论在分析精度, 还是计算速度和内存消耗上, 都大大优于同类方法, 是一种快速、通用的工程弹性动力学问题大规模数值分析方法.      

非线性动力系统线性模型数值计算的Taylor变换法

将非线性动力系统化为连续变化的线性系统,并导出任意自治或非自治非线性动力系统的瞬时线性化方程,该线性方程的连续变化描述了系统的全部复杂动力行为.进一步采用Taylor变换法求解系统的线性化方程,得到一种非线性动力系统数值计算的新方法,避免了指数矩阵展开的乘积运算.计算实例表明该方法在不增加计算机时的前提下,精度高于传统的Houbolt,Wilson-θ及Newmark-β等方法.计算了Duffing方程和van　Pol方程的混沌及周期特性.     

A HIGH ORDER KERNEL INDEPENDENT FAST MULTIPOLE BOUNDARY ELEMENT METHOD FOR ELASTODYNAMICS

基于核无关的快速多极方法, 发展了一种弹性动力学问题的快速、高精度边界元分析方法. 采用基于二次曲面单元的Nyström 离散, 将边界积分方程转化为求和形式, 可以方便地进行加速计算;由于采用二次元, 边界元分析精度很高. 将一种新型快速多极方法用于Nyström 边界元法的加速计算, 该方法的数值实现简便、不依赖于积分方程基本解的表达式, 因此通用性很好;该方法还具有最优的计算量和存储量、精度高且可以控制. 结合Nyström 边界元系数矩阵和快速多极方法转换矩阵的特点, 提出一种大幅度降低边界元内存消耗的策略. 数值结果表明, 该方法无论在分析精度, 还是计算速度和内存消耗上, 都大大优于同类方法, 是一种快速、通用的工程弹性动力学问题大规模数值分析方法.     

带约束平面多刚体系统动力学方程的隐式算法

在多体系统动力学正则方程的基础上建立了平面多体系统正则方程的隐式数值算法。利用平面运动的特性，对正则方程进行了简化，导出了该方程的Ｊａｃｏｂｉ矩阵的一般表达式，给出了Ｒｕｎｇｅ－Ｋｕｔａ多体系统动力学方程隐式数值计算方法。算例表明，该方法是一种计算速度和精度均理想的数值方法。     

质点和多刚体系统的广义动静法

将经典动静法的形式作适当处理后可用于推导一般完整或非完整力学系统的动力学方程,而避免使用动力学函数。这种方程可以采用矩阵形式写出以便面向计算机,文中以树状多刚体系统及操作机器人为例子说明其应用。     

含摩擦滑移铰平面多刚体动力学的HLCP方法

给出了一种基于水平线性互补问题(HLCP)的双边约束滑移铰含摩擦平面多刚体系统动力学的数值计算方法.首先从系统的约束方程出发建立了每个双边约束中两个约束面的法向约束力的互补关系;并利用摩擦余量的概念给出了库仑摩擦定律的互补表达式.然后针对事件驱动算法将该类非光滑多体系统动力学方程中双边约束中法向约束力切换和系统"stick-slip"运动状态切换的判断统一成HLCP的求解,通过定义一组变维数矩阵给出了HLCP形式的一般表达式,便于编程计算.最后通过一个两自由度多体系统的算例验证了该方法的有效性.     

小阻尼转子系统动力学的本征解

本文把小阻尼转子系统的动力学方程化为一种新形式，用摄动法展开该方程，导致一个反对称矩阵的广义本征值问题，因此获得了高精度的本征解。文中给出了一个简例。     

多体系统Lagrange方程数值算法的研究进展

Lagrange方法是建立多体系统动力学方程的普遍方法之一,其方程的形式为常微分方程组或微分 - 代数方程组,数值计算与数值分析是研究多体系统动力学特性的重要方法.本文简要介绍了多体系统动力学方程的第一、二类Lagrange方程和修正的Lagrange方程的基本形式及这些方程的正则形式,着重介绍了正则方程在数值计算中的特点,就多体系统Lagrange方程的隐式算法、辛算法和多体系统动力学特性的数值分析方法(包括数值仿真、Poincar'e映射和Lyapunov指数的计算方法)的研究现状进行了综述.     

一种直接求解的动态载荷识别方法

基于直接求解系统方程的位移响应,提出了一种时域内动态载荷的识别方法.该方法从中心差分格式的有限元离散形式出发,进行分块矩阵变换,将载荷识别问题转化为一个虚构系统方程的响应求解问题.通过数值算例,用两种不同形式的载荷检验了该方法的识别能力,分析了有无零相位滤波器对识别结果的影响.结果表明:所提出的载荷识别方法在位移响应数据含有噪声的情况下,能有效稳定地实现载荷识别;且计算简单快捷,抗噪能力强,不需进行模态分解,不存在模态截取等问题.     

考虑弹簧阻尼作动器解析雅可比矩阵的多刚体动力学分析

弹簧-阻尼-作动器(spring-damper-actuator,SDA)是多体系统中常见的力元,在工程领域中有着广泛的应用.采用绝对坐标方法建立的多体系统动力学控制方程通常是复杂的非线性微分-代数方程组.为了保证数值解的精度和稳定性,通常需要采用隐式算法求解动力学方程,而雅可比矩阵的计算在隐式数值求解过程中至关重要.对于含有SDA的多体系统,SDA造成的附加雅可比矩阵是与广义坐标和广义速度相关的高度非线性函数.目前的很多研究工作专注于广义力向量的计算,然而对附加雅克比矩阵的计算则少有关注.针对含SDA的多刚体系统进行动力学分析,首先基于Newmark算法研究其在动力学方程求解中的雅可比矩阵的构成形式;然后推导SDA的广义力向量对应的附加雅可比矩阵,其中包括广义力向量对广义坐标和对广义速度的偏导数矩阵.最后通过两个数值算例研究附加雅可比矩阵对动力学分析收敛性的影响;数值分析表明:当SDA的刚度、阻尼和作动力数值较大时,SDA导致的附加雅可比矩阵对数值解的收敛性有重要影响;当考虑SDA对应的附加雅可比矩阵时,动力学分析可以以较少的迭代步实现收敛,从而减少分析时间.     

轴向移动悬臂梁高效动力学建模及频率响应分析

轴向移动梁动力学问题具有广泛的工程应用背景,如:机械手、机床主轴、武器身管等.计算轴向移动梁动力学响应是评估结构动力学性能以及最终指导结构设计的一个重要手段.采用Rayleigh-Ritz法、拉格朗日方程推导了轴向移动悬臂梁时变动力学方程.选取幂级数函数构造试函数对轴向移动系统动力问题进行求解.幂级数函数良好的积分与微分性能,使得推导容易以矩阵的形式快速进行,便于符号运算软件直接生成MATLAB程序.由于MATLAB基本数据单位为矩阵,符号软件生成的程序只需经过简单修改便可进行动力学计算.大大缩短了轴向移动梁从建模到动力学分析的时间,过程十分高效.通过四组算例,将本文方法计算得到的动力学响应与文献数据进行对比,对该方法准确性进行了验证,并给出了幂级数函数拟合阶数的选取原则.以此为基础,研究了轴向移动梁的频率响应特性.分为考虑重力与忽略重力两种情况,讨论了轴向振动幅度对其频率响应特性的影响.     

AN EFFICIENT DYNAMIC MODELING METHOD OF AN AXIALLY MOVING CANTILEVER BEAM AND FREQUENCY RESPONSE ANALYSIS

轴向移动梁动力学问题具有广泛的工程应用背景,如：机械手、机床主轴、武器身管等.计算轴向移动梁动力学响应是评估结构动力学性能以及最终指导结构设计的一个重要手段.采用Rayleigh-Ritz法、拉格朗日方程推导了轴向移动悬臂梁时变动力学方程.选取幂级数函数构造试函数对轴向移动系统动力问题进行求解.幂级数函数良好的积分与微分性能,使得推导容易以矩阵的形式快速进行,便于符号运算软件直接生成MATLAB程序.由于MATLAB基本数据单位为矩阵,符号软件生成的程序只需经过简单修改便可进行动力学计算.大大缩短了轴向移动梁从建模到动力学分析的时间,过程十分高效.通过四组算例,将本文方法计算得到的动力学响应与文献数据进行对比,对该方法准确性进行了验证,并给出了幂级数函数拟合阶数的选取原则.以此为基础,研究了轴向移动梁的频率响应特性.分为考虑重力与忽略重力两种情况,讨论了轴向振动幅度对其频率响应特性的影响.