An explicit finite element method for dynamic analysis in three-medium coupling system and its application

In this paper, an explicit finite element method to analyze the dynamic responses of three-medium coupled systems with any terrain is developed on the basis of the numerical simulation of the continuous conditions on the bounda-ries among fluid saturated porous medium, elastic single-phase medium and ideal fluid medium. This method is a very effective one with the characteristic of high calculating speed and small memory needed because the formulae for this explicit finite element method have the characteristic of decoupling, and which does not need to solve sys-tem of linear equations. The method is applied to analyze the dynamic response of a reservoir with considering the dynamic interactions among water, dam, sediment and basement rock. The vertical displacement at the top point of the dam is calculated and some conclusions are given.

Coupled Vibration Analysis of Vehicle-Bridge System Based on Multi-Boby Dynamics

For establishing the refined numerical simulation model for coupled vibration between vehicle and bridge, the refined three-dimensional vehicle model is setup by multi-body system dynamics method, and finite element method of dynamic model is adopted to model the bridge. Taking Yujiang River Bridge on Nanning-Guangzhou railway line in China as study background, the?refined numerical simulation model of whole vehicle and whole bridge system for coupled vibration analysis is set up. The dynamic analysis model of the cable-stayed bridge is established by finite element method, and the natural vibration properties of the bridge are analyzed. The German ICE Electric Multiple Unit (EMU) train refined three-dimensional space vehicle model is set up by multi-system dynamics software SIMPACK, and the multiple non-linear properties are considered. The space vibration responses are calculated by co-simulation based on multi-body system dynamics and finite element method when the ICE EMU train passes the long span cable-stayed bridge at different speeds. In order to test if the bridge has the sufficient lateral or vertical rigidity and the operation stability is fine. The calculation results show: The operation safety can be guaranteed, and comfort?index is “excellent”. The bridge has sufficient rigidity, and vibration is in good condition.

Coupled method of finite and dynamic infinite elements for simulating wave propagation in elastic solids involving infinite domains

This paper deals with the coupled method of finite and dynamic infinite elements for simulating wave propagation in elastic and viscoelastic solids involving infinite domains.This method can be used to simultaneously simulate material complexities in the near field and the infinite extent of the far field.Based on the governing equations of wave motion in two-dimensional and three-dimensional elastic/viscoelastic solids,the mass and stiffness matrices of the dynamic infinite element have been derived.The proposed two-dimensional dynamic infinite element can be used to simulate both the P-wave and the SV-wave propagation within the element,while the proposed three-dimensional dynamic infinite element can be used to simultaneously simulate the Rayleigh wave,P-wave and S-wave propagation within the element.The related simulation results have demonstrated that the coupled method of finite and dynamic infinite elements can be accurately used to simulate,both physically and computationally,wave propagation in elastic/viscoelastic solids involving infinite domains.Thus,this method provides an advanced scientific tool for dealing with both scientific and engineering problems involving infinite domains.

RESEARCH ON SOLID-LIQUID COUPLING DYNAMICSOF PIPE CONVEYING FLUID

On the basis of Hamilton principle. the equation of sonlid-liquid coupling vibration of pipe conveying fluid is deduced. An asymmetrical sonlid-liquid coupling damp matrix and a symmetrical solid-liquid coupling Stiffness matrix are obtained. Using QR method , pipe’s nature frequencies are calculated. The curves of the first four orders of natural frequency-flow velocity of pipe waw given .The influence of flowing velocity ,pressure, solid-liquid coupling damp and solid-liquid coupling stiffness on natural frequency are discussed respectively.The dynamic respondence of the pipes for stepload with different flow velocity are calculated by Newmark method .It is found that,with the flow velocity increased, the nature frequency of the pipes reduced, increased,reduced again and so on.

Dynamic response analysis of a moored crane-ship with a flexible boom

The dynamic response of moored crane-ship is studied. Governing equations for the dynamic response of a crane-ship coupled with the pendulum motion of the payload are derived based on Lagrange’s equations. The boom is modeled based on finite element method, while the payload is modeled as a planar pendulum of point mass. The dynamic response was studied using numerical method. The calculation results show that the large-amplitude responses occur at wave periods near the natural period of the payload. Load swing angle is smaller for crane-ship with flexible boom, in comparison with rigid boom. The ship surge mo- tions have large vibrations for crane-ship with flexible boom, which were not observed for a rigid boom. The analysis identifies the significance of key parameters and reveals how the system design can be adjusted to avoid critical conditions.

Impacts of flexible obstructive working environment on dynamic performances of inspection robot for power transmission line

The rigid-flexible coupling dynamic modeling and simulation of an inspection robot were conducted to study the influences of the flexible obstructive working environment i.e. overhead transmission line on the robot＇s dynamic performance. First, considering the structure of the obstacles and symmetrical mechanism of the robot prototype, four basic subactions were abstracted to fulfill full-path kinematic tasks. Then, a multi-rigid-body dynamic model of the robot was built with Lagrange equation, whil^e a multi-flexible-body dynamic model of a span of lin~ was obtained by combining finite element method （FEM）, modal synthesis method and Lagrange equation. The two subsystem models were coupled under rolling along no-obstacle segment and overcoming obstacle poses, and these simulations of three subactions along different spans of line were performed in ADMAS. The simulation results, including the coupling vibration parameters and driving moment of joint motors, show the dynamic performances of the robot along ftexibile obstructive working path： in flexible obstructive working environment, the robot can fulfill the preset motion goals; it responses slower in more flexible path; the fluctuation of robot as well as driving moment of the corresponding joint in startup and brake region is greater than that in rigid environment; the fluctuation amplitude increases with increasing working environment flexibility.

Dynamic analysis on generalized linear elastic body subjected to large scale rigid rotations

The dynamic analysis of a generalized linear elastic body undergoing large rigid rotations is investigated. The generalized linear elastic body is described in kine- matics through translational and rotational deformations, and a modified constitutive relation for the rotational deformation is proposed between the couple stress and the curvature tensor. Thus, the balance equations of momentum and moment are used for the motion equations of the body. The floating frame of reference formulation is applied to the elastic body that conducts rotations about a fixed axis. The motion-deformation coupled model is developed in which three types of inertia forces along with their incre- ments are elucidated. The finite element governing equations for the dynamic analysis of the elastic body under large rotations are subsequently formulated with the aid of the constrained variational principle. A penalty parameter is introduced, and the rotational angles at element nodes are treated as independent variables to meet the requirement of C1 continuity. The elastic body is discretized through the isoparametric element with 8 nodes and 48 degrees-of-freedom. As an example with an application of the motion- deformation coupled model, the dynamic analysis on a rotating cantilever with two spatial layouts relative to the rotational axis is numerically implemented. Dynamic frequencies of the rotating cantilever are presented at prescribed constant spin velocities. The maximal rigid rotational velocity is extended for ensuring the applicability of the linear model. A complete set of dynamical response of the rotating cantilever in the case of spin-up maneuver is examined, it is shown that, under the ultimate rigid rotational velocities less than the maximal rigid rotational velocity, the stress strength may exceed the material strength tolerance even though the displacement and rotational angle responses are both convergent. The influence of the cantilever layouts on their responses and the multiple displacement trajectories observed in the floating frame is simultaneously investigated. The motion-deformation coupled model is surely expected to be applicable for a broad range of practical applications.

Numerical simulation of aircraft arresting process with an efficient full-scale rigid-flexible coupling dynamic model

The arresting process of carrier-based aircraft is widely recognized as a challenging task,characterized by the highest accident rate among all carrier-based aircraft operations.Dynamic simulation plays a crucial role in assessing the intricate responses of the arresting process,favoring the design of carrier-based aircraft.An efficient and accurate rigid-flexible coupling model for analyzing the dynamic response of the arresting process is proposed.By combining the dynamic characteristics of airframe,landing gear,arresting hook and arresting gear system,the rigid-flexible coupling dynamic model is established to reflect the relative motion of the coupling parts and arresting load.The dynamic model is verified through simulations of landing gear landing drops and by comparing the arresting simulation results with corresponding data in the US military standard.Additionally,simulations of the arresting process under different off-center distance and aircraft yaw angle are conducted to obtain the dynamic response of the aircraft during the arresting process.The result indicates that the rigid-flexible coupling dynamic model proposed is effective for analyzing the arresting dynamics response of carrier-based aircraft.The axial force of the arresting cable on both sides of the hook engagement point,pitch and yaw angle of aircraft are inconsistent under yaw and off-center arresting.The analysis method and obtained results provide valuable references for assessing the dynamic responses of carrier-based aircraft during arresting process and offer valuable in-sights in the design of carrier-based aircraft.

An Integrated Structural Strength Analysis Method for Spar Type Floating Wind Turbine

An integrated structural strength analysis method for a Spar type floating wind turbine is proposed in this paper,and technical issues related to turbine structure modeling and stress combination are also addressed.The NREL-5MW ＂Hywind＂ Spar type wind turbine is adopted as study object.Time-domain dynamic coupled simulations are performed by a fully-coupled aero-hydro-servo-elastic tool,FAST,on the purpose of obtaining the dynamic characteristics of the floating wind turbine,and determining parameters for design load cases of finite element calculation.Then design load cases are identified,and finite element analyses are performed for these design load cases.The structural stresses due to wave-induced loads and wind-induced loads are calculated,and then combined to assess the structural strength of the floating wind turbine.The feasibility of the proposed structural strength analysis method for floating wind turbines is then validated.

大跨度公铁两用斜拉-悬索协作体系桥断索动力响应

研究大跨度公铁两用斜拉-悬索协作体系桥在极端作用下发生断索时桥梁结构及桥上列车的动力响应,以G3铜陵长江公铁大桥为研究对象,建立全桥有限元模型;研究各个位置斜拉索、吊索断裂后桥梁结构的动力响应;研究在列车-桥梁耦合振动作用下,各断索组合工况中桥梁及列车的动力响应.结果表明:斜拉-悬索协作体系桥具有较高的结构刚度及安全冗余;在单根斜拉索、吊索断索后,可保持列车通行;在吊索区域单侧发生4根以上断索,可导致连续断索破坏;在斜拉索-吊索交替及斜拉索区域,单侧断索达到8根也不会发生连续断索破坏;断索动力放大系数取2.0合理,但仅对断索位置附近剩余吊索和斜拉索具有工程应用价值;当断索发生时,若列车经过,列车竖向加速度会发生突变,但仍在安全范围内;剩余斜拉索与吊索动应力增加主要来自于主梁作用,列车作用占比不超过13%.

龙羊峡重力拱坝模态识别与材料参数反演

基于龙羊峡重力拱坝实测振动数据,利用随机子空间法进行了大坝模态参数识别,建立了坝体-地基-库水系统的三维流固耦合有限元模型。根据识别出的一阶频率,分别考虑均质地基及分层地基,反演了坝体与基岩材料参数,并通过有限元软件ABAQUS进行了大坝的动力特征分析。结果表明:龙羊峡重力拱坝的一阶频率为2.778~3.196 Hz,反演得到坝体的弹性模量约为50 GPa~52 GPa;对于均值地基,其弹性模量的反演值为30 GPa;对于分层地基,2 490 m高程以下区域的弹性模量为34 GPa;实测振动数据反演得到的坝体弹性模量和地基弹性模量高于其设计值;龙羊峡重力拱坝的前3阶振型分别为正对称、反对称与正对称,其1阶频率和2阶频率相差较大,与重力坝自振频率分布特征相似。

悬臂同轴双转子系统不平衡耦合振动抑制试验研究

针对悬臂同轴双转子系统不平衡耦合振动突出问题,以某开式转子发动机输出端双转子系统为研究对象,对悬臂双转子系统振动耦合机理进行了理论推导、仿真分析和试验研究。首先,采用有限元法(FEM)建立了双转子系统等效动力学模型和耦合动力学方程,分析了悬臂同轴双转子系统耦合振动机理;然后,利用有限元软件对双转子系统进行了模态分析、临界转速分析以及不平衡响应分析,建立了自敏感系数概念以定量评价系统对不平衡振动的敏感度和耦合状态;最后,根据实际结构设计了悬臂同轴双转子实验台,进行了不同转速比下的不平衡振动抑制实验,并结合实验结果对仿真结果进行了验证。研究结果表明:模拟仿真提出的自敏感系数可作为双转子系统实际振动测点布置和不平衡振动抑制策略的理论依据;通过应用先平衡内转子后平衡外转子的不平衡振动分步抑振策略,可使双转子系统不平衡振动降幅达80%以上。该研究对双转子系统传感器测点布局和系统耦合振动抑制具有一定的指导意义。

基于LAGRANGE方程的深水钻井隔水管–水下井口系统动力分析

为了探究大质量、大刚度防喷器组(blowout preventers,BOPs)对深水钻井隔水管系统动态相应预测精度的影响,根据细长钻井隔水管与刚性防喷器组的结构特点,提出两者刚柔耦合的概念,采用能量法推导隔水管–防喷器组–水下井口系统的动能和势能,采用LAGRANGE方法建立耦合系统动力学理论模型,采用科学计算软件和Newmark-β直接积分法对动力学模型进行数值计算。以南海某深水钻井隔水管为例,开展基于耦合动力学模型的隔水管系统动态响应分析。结果表明,采用本文理论模型得到的隔水管不同位置的节点位移、单元弯矩、上部和下部挠性接头转角时程曲线、整体侧向位移包络线和弯矩包络线等与ABAQUS仿真结果均吻合良好,最大误差为8.8%。此方法可为深水钻井隔水管和水下井口系统动态分析提供参考。

大变幅加卸载下特厚煤层底板断层突水机理模拟研究

特厚煤层采场空间大、扰动范围广,强烈大变幅荷载易导致底板断层破裂加剧并诱发水害。数值模拟是揭示特厚煤层底板断层活化与突水机理的重要方法,准确反映大变幅加卸载下岩体破裂与裂隙水耦合特征是其合理性的关键。构建损伤变量与塑性应变、应力之间的相关关系,得到完整岩块的拉、加压卸载损伤演化方程;以平方拉剪应力与Benzeggagh-Kenane为初始、完全断裂准则,建立塑性位移与强度劣化关系,建立加卸载韧性断裂本构关系;基于实验数据建立贯通裂隙加卸载剪切本构关系。以基本方程与状态方程为基础,结合浸没边界方法,形成裂隙岩体水力学模拟理论。由此编制流体动力学-有限离散元CFD-FDEM耦合程序,模拟研究特厚煤层底板断层活化突水过程。结果表明:CFD-FDEM耦合程序可数值实现特厚煤层底板断层从(准)连续体到离散体转化,以及断层带裂隙水运移过程。底板断层采动破坏包络线呈W形,最深位于断层及其上盘(48.6m),最浅位于断层下盘(23m)。特厚煤层采场底板断层及其上盘受到较大超前集中应力,而后在采空区内大幅卸载,导致该位置出现显著二次破坏,并形成主要导水通道。研究成果为特厚煤层工作面底板断层水害防治提供理论支撑。

页岩水力压裂中多簇裂缝扩展的全耦合模拟

水平井和水力压裂是页岩气开发中的关键技术。对水力压裂中多簇裂缝同时扩展的物理过程进行了数值模拟。采用扩展有限元法(XFEM)模拟岩石中裂缝沿着任意路径扩展,采用有限体积法(FVM)模拟裂缝中流体的流动,并且考虑井筒中流体流动以及在各簇裂缝间的流量动态分配。通过牛顿迭代对全耦合物理过程进行数值求解,重点研究了初始长度不同的两条裂缝的扩展过程,证明较大的射孔摩阻能促进两条裂缝的同时扩展。并通过算例证明了本方法的精度和有效性。

液化场地桩-土地震相互作用振动台试验数值模拟

基于动力Biot理论将饱和砂土模拟为两相介质,正确考虑水和土颗粒之间变形关系,采用完全耦合u-p有限元公式模拟土体位移和孔压,选用多屈服面弹-塑性本构模型模拟黏土、砂土,砂层划分为20-8节点六面体单元,桩处理为梁-柱单元,建立三维有效应力振动台试验分析模型。通过振动台试验结果,验证该模型的正确性,显示该方法能够准确地模拟可液化场地桩-土-结构动力相互用的重要特征,这对于实际工程设计具有一定的参考意义。

轨枕空吊对轨枕动态性能的影响

为研究轨枕空吊对轨枕动态响应的影响,建立了车辆-轨道垂向耦合动力学模型.假设车辆模型为多刚体系统,用Euler梁模拟离散支撑的钢轨,轨枕视为弹性地基上的Euler梁,利用显式积分法求解车辆-轨道非线性动力学方程.结果表明:当轨枕部分悬空或完全悬空时,影响轨枕最大弯矩的位置;轨枕空吊数目越多,轨枕的弯矩越大;轨枕底座的完全悬空对邻近空吊区的轨枕动态行为影响较大.

搅拌摩擦焊接过程中材料流动形式

采用完全热力耦合模型对搅拌摩擦焊接过程进行模拟,并详细分析了搅拌摩擦焊接过程中的材料流动形式。结果表明,模型可以成功预测搅拌摩擦焊接过程材料流动和温度分布情况。通过对搅拌头周围材料流动的研究,分析了搅拌摩擦焊接过程中飞边现象形成的主要原因。研究了搅拌摩擦焊接构件不同厚度上材料的三维流动形式,通过与二维情况的比较证实,二维情况下的材料流动数值模拟结果对应于搅拌摩擦焊接构件靠近下表面部分的材料流动情况。从等效塑性应变的分布也能证实搅拌头轴肩对靠近上表面的材料行为具有明显影响,而对下表面附近材料行为影响较弱,从而说明二维情况对应三维情况靠近下表面的部分。

人-机-路环境下摩托车刚柔耦合系统动力学研究

目前摩托车系统动态特性研究一般采用多刚体动力学方法,没有考虑摩托车车架空间柔性体、人和路对摩托车动态特性的影响,致使其动态特性与实际存在较大的差别。把车架作为空间柔性体,考虑发动机和路面激励,建立了人-机-路环境下摩托车系统刚柔耦合动力学模型,对其加速、制动、转向等工况下系统的动态特性进行了仿真和评价,试验验证了理论分析的正确性;对摩托车通过国家B级路面和凸台两种工况下系统的平顺特性进行了仿真和道路试验,研究结论对摩托车整车开发及动态性能优化有一定指导意义。

浮体与系泊系统的耦合动力分析

浮体与系泊系统的耦合动力分析对各种新型海洋平台和海洋结构物的应用具有重要意义。在浮体和系泊系统的时域耦合分析中,研究并提出了一种预报浮体和系泊系统非线性运动响应和缆索张力的分析技术。考虑瞬时物面变化引起的各种非线性因素,在时域内建立浮体的非线性运动方程并求解,其中,入射波力和静恢复力按瞬时物面条件计算,辐射力和绕射力按浮体在瞬时平均湿表面为平衡位置做简谐振荡运动来计算;考虑系泊缆索受到的非线性流体动力及大位移,采用三维弹性梁理论将缆索模型化,把缆索看作是等拉伸、等弯曲刚度的细长结构,并采用非线性有限元法对缆索的运动方程进行求解;将浮体的分析方法和系泊缆索的分析方法结合,建立了浮体和系泊系统非线性运动响应和缆索张力的耦合计算方法;对缆索,使用高阶单元,只需要较少的单元就可以得到收敛解;对耦合动力分析方程,采用隐式、显式混合的数值格式,最终将方程化为显式格式进行求解,并且合理排列和控制求解次序,这样大大提高了计算效率。对系泊浮体的固有周期、各种工况下大幅非线性运动、缆索张力等进行了计算,并进行对比研究,可以看出,该文的方法既适用于缆索小变形条件下,也适用于缆索大变形条件下浮体与系泊系统的耦合动力分析。