一种O(n)算法复杂度的递推绝对节点坐标法研究

相比于传统的浮动坐标法,绝对节点坐标法(absolute nodal coordinate formulation,ANCF)在处理柔性体非线性大变形问题上具有显著优势,但是对于ANCF的求解目前主要依据拉格朗日方程等分析力学原理建立微分代数方程(differential algebraic equation,DAE)进行,其算法复杂度为O(n2)或O(n3)(n为系统自由度),且求解过程存在位置或速度的违约问题.据此,研究了一种O(n)算法复杂度的递推绝对节点坐标法(recursive absolute nodal coordinate formulation,RANCF).该方法采用ANCF描述大变形柔性体,借鉴铰接体递推算法(articulatedbody algorithm,ABA)思路建立多柔体系统逐单元的运动学和动力学递推关系,得到微分形式的系统动力学方程(ordinary differential equation,ODE).在ODE方程中,系统广义质量阵为三对角块矩阵,通过恰当的矩阵处理,可以得到逐单元求解该方程的递推算法.在此基础上,给出了RANCF算法的详细流程,并对流程中每个步骤进行了细致的算法效率分析,证明了RANCF是算法复杂度为O(n)的高效算法.RANCF方法保留了ANCF对大转动、大变形多柔体系统精确计算的优点,同时极大地提升了算法效率,特别在处理高自由度复杂多柔体系统中具有显著优势.并且该方法采用ODE求解,无DAE的违约问题,因此具有更高的算法精度.最后,在算例部分,通过MSC.ADAMS仿真软件、能量守恒测试、算法复杂度曲线对RANCF的正确性、计算精度和计算效率进行了验证.     

基于绝对节点坐标法的平面梁有限变形下变形重构

现有的对有限变形条件下柔性结构变形重构的研究往往单纯基于曲率与应变间的几何关系, 同时忽略了被测体的纵向变形及其与弯曲变形的耦合效应. 为得到一种更加精确且能借助现有的力学工具进行应用方向扩展的变形重构方法, 以平面梁为对象, 借鉴变形重构逆有限元法的思想, 将平面梁的变形重构问题视作一类最优化问题. 首先, 通过引入绝对节点坐标法(absolute nodal coordinate formulation, ANCF)对柔性结构大变形下非线性的平面梁应变?位移关系进行精确描述, 构造了一种逆梯度缩减ANCF平面索梁单元. 然后, 对此逆ANCF单元进行改进, 在简化节点自由度的同时通过引入罚函数确保单元节点处的曲率连续性, 既保证了本问题的适定性, 也提升了最终解的精确性. 最后, 基于该单元利用Newton法构造了平面梁有限变形下变形重构问题的两种求解算法, 即逐单元算法和多单元整体算法, 以实现不同需求下的稳定求解. 数值仿真结果表明, 本方法在大变形条件下的变形重构误差小于1%, 而且在测点较少的情况下依然保持较高的精度, 同时验证了本方法的收敛性与计算效率.      

柔性多体系统动力学绝对节点坐标方法研究进展

阐述了多体系统动力学理论的研究背景,指出了多种传统的柔性多体系统动力学研究方法的不足.系统地从4个方面回顾了柔性多体系统动力学绝对节点坐标方法诞生十几年以来的研究进展,即:单元研究进展、系统动力学方程求解数值算法研究进展、非线性材料多体系统动力学研究进展以及相关的应用研究进展.最后提出了值得进一步研究的问题.      

基于绝对节点坐标法的弹性线方法研究

相比于浮动坐标系法,绝对节点坐标法(absolute nodal coordinate formulation, ANCF)在处理柔性体非线性大变形问题上具有显著优势, ANCF将单元节点坐标定义在全局坐标系下,采用斜率矢量代替节点转角坐标,具有常数质量阵,不存在科氏离心力等优点,然而弹性力阵为非线性项,其求解将比较耗时且占用资源.据此,在弹性力求解方法中,引入弹性线方法 (elastic line method, ELM),该方法将格林–拉格朗日应变张量定义在中心线上,采用曲率公式来定义弯曲应变,转角公式来定义扭转应变.同时采用有限元法对三维柔性梁位移场进行离散,求解梁单元常数质量阵、广义刚度阵、广义力阵,进而得到单元的动力学方程,通过转换矩阵得到三维梁的动力学方程.接着从理论上指出连续介质力学方法 (continuum mechanics method, CMM)和弹性线方法在求解弹性力上的不同点,并编制动力学仿真软件.最后分别采用连续介质力学方法和弹性线方法对柔性单摆以及履带式车辆的动力学问题进行仿真分析,结果表明:弹性线方法能在保证精度的前提下有效提高计算效率.     

弹性力学轴对称问题的有限元线法

给出了解弹性力学空间轴对称问题的有限元线法的基本理论。该法包括了２－４条结线的等参数单元，沿结线方向的两点边值问题采用插值矩阵法解之。算例表明，本法具有良好的收敛性和较高的计算精度。     

线性黏弹性球面波的特征线分析

基于ZWT黏弹性本构方程建立了体现高应变率效应的黏弹性球面波的控制方程组,包含5个偏微分方程,解5个未知量v、r、、r和。采用特征线法,问题转化为解3族特征线上的5个常微分方程,物理上图像清晰,数学上易于求解。特征线数值分析显示,黏弹性球面波的衰减和弥散效应超过线弹性球面波。球面扩散引起的环向拉应力是导致介质拉伸破坏的主因。进一步还针对强间断黏弹性球面波得出其衰减特性的解析表达式,表明这种更强的衰减特性是几何扩散效应和本构黏性效应两者共同作用的后果。     

基于节点密度的自适应EFG法拓扑优化方法研究

针对无网格法拓扑优化中设计变量庞大、计算效率低的问题,提出了基于节点密度的自适应EFG拓扑优化方法.该方法以节点相对密度为设计变量,建立了以结构柔度最小为优化目标的自适应EFG拓扑优化模型,给出了拓扑优化中自适应节点加密准则、节点加密方法和新增节点设计变量更新规则.算例结果表明,采用论文提出的自适应EFG拓扑优化方法进行结构拓扑优化时,可以减少结构分析和优化的设计变量数目,提高优化效率.     

Kirchhoff弹性杆分析动力学的准坐标表达

作为DNA的力学模型,依据Kirchhoff动力学比拟思想建立的弹性细杆的分析力学方法已从静力学深入到动力学。由于静力学平衡微分方程与刚体动力学相当,因此,弹性细杆动力学的分析力学方程必是以弧坐标和时间为双自变量的偏微分方程。以横截面的形心速度以及弯扭度和角速度沿主轴的分量为准速度,定义了准坐标,导出了准坐标的微分和变分运算的交换关系。从Hamilton原理出发,利用准坐标的微分和变分运算的交换关系,导出了Kirchhoff弹性杆动力学准坐标下的Boltzmann-Hamel方程,并由此导出Lanrange方程。指出了Boltzmann-Hamel方程显式即为弹性杆动力学的Kirchhoff方程。定义关于弧坐标和时间的正则变量和Hamilton函数,导出Boltzmann-Hamel方程的正则形式。本文结果是以弹性杆静力学和刚性杆动力学为其特例。作为例子,建立了垂挂的在重力作用下作平面运动的弹性细杆的动力学微分方程以说明本文方法的应用。     

弹性力学的复变量数值流形方法

数值流形方法通过引入数学和物理双重网格,将插值域和积分域分别定义在两个不同的覆盖上来完成系统能量泛函积分运算. 当采用高阶函数构造位移函数时,广义节点自由度将大大增加. 在求解系统的平衡方程中,运算量是与自由度的三次方成正比的,因此数值流形方法的计算量是较大的. 为此,在复变量理论的基础上,采用一维基函数建立二维问题的逼近试函数,然后将其应用于弹性力学的数值流形方法,提出了复变量数值流形方法,推导了弹性力学的复变量数值流形方法的公式. 与传统的数值流形方法相比,复变量数值流形方法具有计算量小、精度高的优点.      

INTEGRATION OF NON-UNIFORM RATIONAL B-SPLINES GEOMETRY AND RATIONAL ABSOLUTE NODAL COORDINATES FORMULATION FINITE ELEMENT ANALYSIS

This investigation is intended to develop a computer procedure for the integration of NURBS geometry and the rational absolute nodal coordinate formulation （RANCF） finite element analysis. A linear transformation is given that can be used to convert the NURBS curve to RANCF cable element mesh retaining the same geometry and the same degree of continuity, including the discussion of continuity control and mesh refinement. The green strain tensor is used to establish the nonlinear dynamic equations with numerical examples to demonstrate the use of the procedure in the dynamic analysis of flexible bodies.     

完备化ANCF薄板单元及在钢板弹簧动力学建模中的应用

绝对节点坐标方法已在多体系统动力学研究中广泛应用, 但常用来描述板壳类结构的薄板单元, 由于梯度不完备而无法直接用于带有初始弯曲参考构型的柔性体变形描述. 为避免全参数板单元建立车辆钢板弹簧模型时存在的严重截面闭锁问题, 拟采用薄板单元用于板簧建模. 为此, 探索了将现有绝对节点坐标薄板单元纳入一般连续介质力学弹性力表达的方法, 采用中面上单位法向量作为单元厚度方向的梯度向量, 从而得到了完备化的薄板单元及其描述初始弯曲构型时消除初应变的方法. 在此基础上通过定义簧片的未变形构型, 在钢板弹簧中引入可控的预应力, 实现对钢板弹簧装配过程的准确模拟. 通过数值算例验证了本方法的正确性. 建立了车辆钢板弹簧模型, 通过建立在簧片上的局部坐标系实现接触点的跨单元搜索, 并采用惩罚函数法和平滑化的库伦摩擦模型施加簧片间的接触力. 引入参考节点的概念建立了整合车身与吊耳及其机构运动关系的刚柔耦合模型.}}      

基于节点独立变量的连续体结构动态拓扑优化

提出基于节点独立变量的连续体结构动态拓扑优化方法.以动态结构响应量最小或最大为目标,体积比为约束,建立了动态结构拓扑优化模型.基于数字图像处理过滤技术得到清晰、边界光滑和体现网格无关性的优化结果.通过二维结构数值算例对理论方法进行验证.结果表明,该方法在连续体结构动态拓扑优化设计中具有可行性和有效性.     

含分数阻尼特性元件的多体系统动力学研究

在绝对节点坐标体系下研究了具有分数导数阻尼特性元件的多体系统动力学建模、求解问题. 采用基于绝对节点坐标的无闭锁效应剪变梁单元离散柔性构件,建立了含常数质量矩阵的系统动力学方程, 并采用数值耗散可控的广义a方法求解. 通过数值算例计算,对比研究了算法参数与阻尼项的分数指数对系统动力学响应的影响规律.该方法可以进一步扩展到众多工程实际问题研究中.      

基于Green函数的动态载荷区间识别方法研究

基于位移响应时间序列,提出了对结构动态载荷时间历程进行识别的区间方法.载荷在时域内可用一系列脉冲响应函数来表示,系统的响应是载荷与单位脉冲响应函数(Green函数)的卷积分.采用区间向量定量化不确定但有界的位移响应时间序列,并将Green函数离散,构造区间线性方程组,并利用区间摄动法进行求解,从而得到待识别载荷时间历程的区间估计.对不同测点下识别得到的载荷区间进行集员识别,改进识别结果精度.最后通过对十杆桁架结构和悬臂梁结构的动态载荷识别,验证了本文方法的有效性和适用性.     

基于遗传算法的弹性地基加肋板肋梁无网格优化分析

基于遗传算法及一阶剪切理论, 提出一种弹性地基上加肋板肋条位置优化的无网格方法. 首先, 通过一系列点来离散平板及肋条, 并用弹簧模拟弹性地基, 从而得到加肋板的无网格模型; 其次, 基于一阶剪切理论及移动最小二乘近似原理导出位移场, 求出弹性地基加肋板总势能; 再次, 根据哈密顿原理导出结构的弯曲控制方程, 并通过完全转换法处理边界条件; 最后, 引入遗传算法和改进遗传算法, 以肋条的位置为设计变量、弹性地基板的中心点挠度最小值为目标函数, 对肋条位置进行优化达到地基板控制点挠度最小的目的. 以不同参数、载荷布置形式的弹性地基加肋板为例, 与ABAQUS有限元结果及文献解进行比较. 研究表明, 采用所提出的无网格模型可有效求解弹性地基上加肋板弯曲问题, 结果易收敛, 同时基于遗传算法与改进混合遗传算法所提出的无网格优化方法均可有效优化弹性地基加肋板肋条位置, 后者计算效率相对较高, 只进行了三次迭代便可获得稳定的最优解, 此外在优化过程中肋条位置改变时只需要重新计算位移转换矩阵, 又避免了网格重构.      

弹性薄板分析的条形传递函数方法

提出一种用于矩形弹性薄板变形分析的条形传递函数方法．一个矩形区域首先沿某一个方向被剖分成若干个条形子域,分割这些子域的直线称为结线,在结线上定义位移函数,它是结线坐标的一维函数,结线的两个端点称为结点．为适应复杂边界条件,在边界结线上定义若干结点,该结线的位移函数用结点位移参数插值表示．每个条形子域的变形用结线位移函数和适当的插值函数（形函数）表示．结线位移函数和结点位移参数满足的平衡微分方程及代数方程由变分原理给出     

弹性力学的实时神经计算原理与数值仿真

针对现代结构分析的特点，提出了基于神经网络结构的弹性力学分析原理；给出求解该问题的网络———改进的Ｈｏｐｆｉｅｌｄ和ＴＨ网络．提出用ＢＰ网络来实现单元刚度矩阵的实时计算．最后，对两个简单结构的弹性力学神经计算进行了数值仿真