橡胶弹簧非线性刚度的有限元解

本文探讨了用双重非线性有限元手段，求解橡胶弹簧的非线性刚度。橡胶的本构关系采用了广义虎克定律和基于应变能的超弹性材料的本构关系联合使用的方案，用单元内的静水压力ｐ，来修正广义虎克定律的弹性常数Ｅ和Ｇ。还提出了带刚性域的超单元和波前凝聚法两种新的办法，获得橡胶弹簧的刚度，使内存用量和ＣＰＵ时间大幅度下降。因而几何、物理双重非线性橡胶弹簧的计算，能在微机上实现。文内选用英国ＤＵＮＬＯＰ橡胶公司橡胶元件     

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

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

含随机参数的非线性黏弹性有限元计算

进行了粗粒土与结构接触面单调和循环加载试验，基于宏细观测量结果, 扩展了 损伤概念以 描述该类接触面在受载过程中的物态演化, 及由于物态演化导致的力学特性从初始状态到最终 稳定状态的连续变化过程. 揭示了接触面损伤的细观物理基础主要是接触面内土的颗粒破碎 和剪切压密这两种物态演化；指出接触面的剪胀体应变可以划分为可逆性和不可逆性剪胀体 应变两部分，其中不可逆性剪胀体应变可作为接触面损伤发展的宏观量度，因此其归一化 形式可作为一种损伤因子的定义；提出了建立粗粒土与结构接触面一种损伤本构关系的基本思路.     

船舶碰撞数值仿真的附加质量模型

根据船舶碰撞的运动滞后和局部损伤特性 ,采用附加质量处理流体 结构耦合作用 ,用详细的有限元模型表达撞击船首和被撞船侧的直接涉撞区结构 ,而将非碰撞区的船体结构视为刚体 ,形成一个附加质量模型用于碰撞仿真计算。附加质量模型不仅建模工作量小 ,而且CPU时间也大大减少。算例表明 ,附加质量模型与流固耦合模型的分析结果具有良好的一致性 ,可以满足一般的工程精度要求。     

超磁致伸缩换能器的磁滞非线性动力学仿真

在超磁致伸缩材料输出端位移假设的基础上,建立了以输出顶杆为研究对象、考虑材料内部磁滞非线性及预压弹簧特性的超磁致伸缩换能器非线性动力学模型; 并将Simulink仿真系统应用到超磁致伸缩换能器系统动力学仿真中. 仿真结果表明: 换能器模型系统为稳定周期运动系统,并具有滞回非线性特性;弹簧非线性项中的平方项是影响换能器模型系统的主要因素. 经验证,数值计算结果与文献给出的试验测试结果吻合较好.     

多孔硅橡胶材料有限元计算的切线本构矩阵

为进行多孔硅橡胶材料非线性有限元计算，针对用伸长比表示的应变能函数，推导了多孔硅橡胶材料非线性有限元计算中切线本构关系的具体形式。并采用便于编程的矩阵记法，给出了非线性有限元计算中克希荷夫应力和格林应变之间的切线关系式。     

非线性碰摩力对碰摩转子分叉与混沌行为的影响

研究了具有非线性碰摩力的转子局部碰摩的分叉与混沌运动,利用计算机仿真对某发动机转子的碰摩故障进行了数值模拟,讨论了转子系统参数的变化对转子混池运动状态的影响,并与碰摩实验结果进行了比较,发现了具有非线性碰摩力的转子局部碰摩转子系统的各种多周期运动和混沌运动及其演变过程。     

橡胶楔体与刚性缺口接触的大变形分析

利用一种新的橡胶材料应变能函数，对橡胶楔体与刚性缺口接触大变形问题进行了分析。得到了接触尖点附近变形的奇异性特征，给出了奇异性指数与材料常数、橡胶楔体角度、刚性缺口角度之间的关系式。同时编制了大变形有限元程序，计算得到了与理论解一致的结论。     

不可压缩流动的移动粒子半隐式法数值模拟

介绍了一种新的移动粒子半隐式方法。该方法采用Lagrange描述,用一系列离散粒子代替流体,通过半隐式求解和时间推进法来进行计算,弥补了传统网格法对复杂形状、大变形、高速撞击等情况下网格划分和重构过于繁琐的缺点,在模拟不可压缩流动问题中得到了较好的效果。     

有限元应力场可视化的一种新方法

提出一种画二维和主等参单元彩色云图的新方法。首先，根据单元覆盖象素的数量在正方形母元上均匀布置足够数量的点。第二，用单元的结点坐标、结点应力以及等参扦值函数，插值计算这些点的总体坐标及应力。第三，根据预先设立的颜色与应力之间的函数关系，求出每个点对应的象素的颜色。最后，逐个把单元覆盖的象素设置为求出的颜色。与消陷技术以及画家算法结合起来，这种方法可以清晰画出三维应力场的云图。实例表明，本文方法简单、准确，易于实现。     

Elasto-plastic analysis of Reissner-Mindlin plates by the use of hybrid stress FEM with layered model

Based on a _mR -type variational formulation featuring a cross-depth layered model in conjunction with a mechanical sub-element for simulating the material constitution, the cross-depth plasticity development of the Reissner-Mindlin plate is investigated by following the loading process. A 4-node quadrilateral hybrid-stressc ⁰-continuous plate-bending element HPT-9 is formulated. Numerical examinations demonstrate its remarkable characteristic behavior in being free from spurious kinematic mode, capable of alleviating locking difficulties as the thin plate limit is approached and providing numerical results with remarkable accuracy and computational efficiency over Spilker's counterpart LH₄. An elasto-plastic analysis of Reissner-Mindlin plates has justified the validity and effectiveness of the present scheme in depicting the cross-depth plasticity development following the loading process.The Project is Supported by National Natural Science Foundation of China.     

Dynamics of repelling soliton collisions in coupled Schrödinger equations

The system of Coupled Nonlinear Schrödinger's Equations (CNLSE) is solved numerically by means of a conservative difference scheme. Values of the cross-modulation parameter, α₂, are chosen to induce repelling collisions. The resulting waveforms are fitted to soliton profiles and studied via internal parameters including phase velocity to gain insight into the dynamics of the repelling collisions. A decaying oscillation of the post-collision profiles is observed and investigated in terms of varying α₂ and relative initial velocity of colliding profiles.     

Tension of a cylindrical membrane partially stretched over a rigid cylinder

We study a contact problem with friction for a hyperelastic long thin-walled tube. One end of the tube is placed over an immovable, rough, rigid cylinder and an axial force is applied to another end. We assume the deformation of the tube is finite and axisymmetric. The tube is modeled by a semi-infinity cylindrical membrane. The axial force tends to a constant value at large distances from the inclusion. The membrane is made of an incompressible, homogeneous, isotropic elastic material. A contact between the membrane and the rigid cylinder is with a dry friction. The membrane will not slide off the cylinder only by friction and at a sufficient contact area. The friction is described by Coulomb’s law. We study a minimum length of the membrane which is in contact with the rigid cylinder and is needed to hold the membrane on the rigid cylinder. We obtain an explicit solution for the Bartenev–Khazanovich (Varga) strain–energy function and numerical results for the Mooney–Rivlin and Fung models.     

有限元网格自组织神经网络的自动生成

利用多个自组织神经网络的同时训练,可以一次性自动地生成计算区域内与计算区域边界上所有网格节点,尔后按FGT法的基本思想将节点连接成所需的三角形或四边形单元。由于在自组织神经网络中采用了改善的目标函数,节点的分布可以实现自适应调节以反映网格疏密分布的要求。文末几个算例验收证明本算法具有自适应性,适用于凸域、凹域、多连通域等多种情况。     

弹性接触问题的常刚度迭代方法及其加速

通过引入接触间隙荷载，在荷载项中反映接触状态变化的影响，避免了重新形成接触刚度矩阵，构造了弹性接触问题的常刚度迭代格式。然后借助于不动点迭代的Steffensen加速技术给出了一种求解弹性接触问题的常刚度迭代加速算法。文中最后给出两个算例表明了本文方法的有效性。     

弹性与弹塑性问题的有限元与边界元耦合解法

本文提出了一种新的弹性与弹塑性问题的对称耦合解法。根据分区广义变分原理,直接导出问题的求解方程式。通过典型算例,验证了该方法的有效性,本文建议的方法与超单元形式的耦合法相比,在理论上比较直接,在计算上更为经济。