Mechanical analysis of deepwater drilling riser system based on multibody system dynamics

A multibody system including a drilling riser system,tensioners and a floating platform is key equipment for offshore oil and gas drilling.Most of the previous studies only focus on the drilling riser system rather than the multibody system.Mechanical characteristics of the deepwater drilling riser system cannot be analyzed accurately in a simplified model.Therefore,a three-dimensional multibody analysis program is developed.The static and dynamic characteristics of the deepwater drilling riser system under different platform motions are analyzed based on the developed program.The results show that the static displacement of the riser system with tensioners is smaller than that without tensioners,which means the tensioners can suppress the deformation of the riser system.Under surge and sway motions of the platform,the dynamic displacement of the riser system with tensioners is also smaller than that without tensioners due to the tensioner suppression effect.Besides,the heave motion induces a uniform axial vibration of the riser system,while roll and pitch motions excite the riser system to vibrate laterally.Compared with the stress amplitude due to surge and sway motions,the stress amplitude of the riser system due to heave,roll and pitch motions is relatively small but cannot be neglected.

Computational dynamics of soft machines

Soft machine refers to a kind of mechanical system made of soft materials to complete sophisticated missions, such as handling a fragile object and crawling along a narrow tunnel corner, under low cost control and actuation. Hence, soft machines have raised great challenges to computational dynamics. In this review article, recent studies of the authors on the dynamic modeling, numerical simulation, and experimental validation of soft machines are summarized in the framework of multibody system dynamics. The dynamic modeling approaches are presented first for the geometric nonlinearities of coupled overall motions and large deformations of a soft component, the physical nonlinearities of a soft component made of hyperelastic or elastoplastic materials, and the frictional contacts/impacts of soft components, respectively. Then the computation approach is outlined for the dynamic simulation of soft machines governed by a set of differential-algebraic equations of very high dimensions, with an emphasis on the efficient computations of the nonlinear elastic force vector of finite elements. The validations of the proposed approaches are given via three case studies, including the locomotion of a soft quadrupedal robot, the spinning deployment of a solar sail of a spacecraft, and the deployment of a mesh reflector of a satellite antenna, as well as the corresponding experimental studies. Finally, some remarks are made for future studies.

Numerical and experimental studies on impact dynamics of a planar flexible multibody system

In this paper a computational methodology on impact dynamics of the flexible multibody system is presented. First, the floating frame of reference approach and nodal coordinates on the basis of finite element formulation are used to describe the kinematics of planar deformable bodies. According to the kinematic description of contact conditions, the contact constraint equations of planar flexible bodies are derived. Based on the varying topology technique the impact dynamic equations for a planar multibody system are established. Then the initial conditions of the equations in each contact stage are determined according to the discontinuity theory in continuum mechanics. The experiments between the aluminum rods are performed to check the correctness of the proposed method. Through the comparison between the numerical and experimental results the proposed method is validated. Experimental results also show that the impulse momentum method cannot accurately predict the complex impact dynamic phenomena and the continuous model may lead to a serious error when used to simulate the impact problems with significant wave propagation effects.

AN AUTOMATIC CONSTRAINT VIOLATION STABI- LIZATION METHOD FOR DIFFERENTIAL/ ALGEBRAIC EQUATIONS OF MOTION IN MULTIBODY SYSTEM DYNAMICS

A new automatic constraint violation stabilization method for numerical integration of Euler_Lagrange equations of motion in dynamics of multibody systems is presented. The parameters α,β used in the traditional constraint violation stabilization method are determined according to the integration time step size and Taylor expansion method automatically. The direct integration method, the traditional constraint violation stabilization method and the new method presented in this paper are compared finally.

An approach for the coupled simulation of machining processes using multibody system and smoothed particle hydrodynamics algorithms

The prediction accuracy of a simulation method is limited by its theoretical background. This fact can lead to disadvantages regarding the simulation quality when investigating systems of high complexity, e.g. containing components showing a fairly different behavior. To overcome this limitation, co-simulation approaches are used more and more, combining the advantages of different simulation disciplines. That is why we propose a new strategy for the dynamic simulation of cutting processes. The method couples Lagrangian particle methods, such as the smoothed particle hydrodynamics （SPH） method, and multibody system （MBS） tools using co-simulations. We demonstrate the capability of the new approach by providing simulation results of an orthogonal cutting process and comparing them with experimental data. @ 2013 The Chinese Society of Theoretical and Applied Mechanics. [doi：10.1063/2.1301305]

Fast simulation of flexible multibody dynamics with electric-hydraulic drive system

Hardware in the loop simulation （HILS） has been investigated in the field of the multibody dynamics （MBD）, which combined the MBD simulation with the actual mechanical system. The fast simulation is necessary for the HILS system in order to require the real time simulation. This paper presents a fast simulation technique using the domain decomposition method with the iteration in the flexible multibody system in which flexible linkage system and electro-hydraulic drive system are coupled with each other. C 2013 The Chinese Society of Theoretical and Applied Mechanics.[doi：10.1063/2.1301301]

A METHOD FOR SOLVING THE DYNAMICS OF MULTIBODY SYSTEMS WITH RHEONOMIC AND NONHOLONOMIC CONSTRAINTS

The solution of the dynamic problem of multibody systems subject to rheonomic and nonholonomic constraints is achieved by applying singular value decomposition of the constraint matrix and projections of the dynamic equations of the systems along the feasible and unfeasible directions of the constraints. Formula to solve the constraint reaction forces and a method to avoid the violation of the constraints are also given.The solution does not rely on coordinates used to describe the systems and is computational efficitive example is finally presnted.

Dynamic Investigations on the Wear Behavior of a 3D Revolute Joint Considering Time-Varying Contact Stiffness: Simulation and Experiment

The existence of the relative radial and axial movements of a revolute joint’s journal and bearing is widely known.The three-dimensional(3D)revolute joint model considers relative radial and axial clearances;therefore,the freedoms of motion and contact scenarios are more realistic than those of the two-dimensional model.This paper proposes a wear model that integrates the modeling of a 3D revolute clearance joint and the contact force and wear depth calculations.Time-varying contact stiffness is first considered in the contact force model.Also,a cycle-update wear depth calculation strategy is presented.A digital image correlation(DIC)non-contact measurement and a cylindricity test are conducted.The measurement results are compared with the numerical simulation,and the proposed model’s correctness and the wear depth calculation strategy are verified.The results show that the wear amount distribution on the bearing’s inner surface is uneven in the axial and radial directions due to the journal’s stochastic oscillations.The maximum wear depth locates where at the bearing’s edges the motion direction of the follower shifts.These find-ings help to seek the revolute joints’wear-prone parts and enhance their durability and reliability through improved design.

A Geometrically Exact Triangular Shell Element Based on Reproducing Kernel DMS-Splines

Tomodel amultibody systemcomposed of shell components,a geometrically exact Kirchho-Love triangular shell element is proposed.The middle surface of the shell element is described by using the DMS-splines,which can exactly represent arbitrary topology piecewise polynomial triangular surfaces.The proposed shell element employs only nodal displacement and can automatically maintain C1 continuity properties at the element boundaries.A reproducing DMS-spline kernel skill is also introduced to improve computation stability and accuracy.The proposed triangular shell element based on reproducing kernel DMS-splines can achieve an almost optimal convergent rate.Finally,the proposed shell element is validated via three static problems of shells and the dynamic simulation of aexible multibody system undergoing both overall motions and large deformations.

A GENERAL PROCEDURE TO CAPTURE THE "DYNAMIC STIFFNESS

A general procedure to capture the 'dynanmic Stiffness' is presented in this paper. The governing equations of motion are formulated for an arbitrary flexible body in large overall motion based on Kane's equations . The linearization is performed peroperly by means of geometrically nonlinear straindisplacement relations and the nonlinear expression of angular velocity so that the 'dynamical stiffness' terms can be captured naturally in a general tcase. The concept and formulations of partial velocity and angular velocity arrays of Huston's method are extended to the flexible body and form the basis of the analysis. The validity and generality of the procedure presented in the paper are verified by numerical results of its application in both the beam and plate models.

Simulation of Restraint System Performance upon the Occupant's Response During Impact

Aim To study the influence of restraint system performance upon the occupant's response during impact, and provide a scientific base for occupant restraint system design. Methods \ In the light of basic theory of multibody system dynamics and impact dynamics on the basis of classical theory of impact, R W method is adopted to construct the vehicle occupant system model consisting of fourteen rigid bodies, thirty seven DOFs and slip joints for the simulation. A software named SVC3D(3 dimensional simulation of vehicle crash) is developed in the FORTRAN language. Results\ The results of simulation have a good coincidence with those of tests and the restraint system with low elongation webbing and equipped with pretensioner provides better restraint effect for the occupant. Conclusion\ The model of vehicle occupant multibody system and SVC3D are suitable for use. Occupant should be belted with low elongation webbing to a certain degree and occupant restraint system should be equipped with pretensioner.

A Study on the Dynamic Analysis of A Tracked Vehicle for Ocean Mining on the Deep Seabed

A study is presented on the dynamic analysis of a tracked vehicle for mining on the deep seabed of very soft soil. Equations for the interaction between the track and extremely soft seabed are employed to develop a track/soil interaction module called TVAS. The vehicle is modeled as a multibody dynamic system by the use of a multibody dynamic analysis program. The module developed is cooperated with the multibody dynamic analysis program with a user-defined subroutine. The dynamic behavior and the conceptual design of the mining vehicle on the deep seabed are investigated.

Multibody dynamic analysis using a rotation-free shell element with corotational frame

Rotation-free shell formulation is a simple and effective method to model a shell with large deformation. Moreover, it can be compatible with the existing theories of finite element method. However, a rotation-free shell is seldom employed in multibody systems. Using a derivative of rigid body motion, an efficient nonlinear shell model is proposed based on the rotation-free shell element and corotational frame. The bending and membrane strains of the shell have been simplified by isolating deformational displacements from the detailed description of rigid body motion. The consistent stiffness matrix can be obtained easily in this form of shell model. To model the multibody system consisting of the presented shells, joint kinematic constraints including translational and rotational constraints are deduced in the context of geometric nonlinear rotation-free element. A simple node-to-surface contact discretization and penalty method are adopted for contacts between shells. A series of analyses for multibody system dynamics are presented to validate the proposed formulation. Furthermore,the deployment of a large scaled solar array is presented to verify the comprehensive performance of the nonlinear shell model.

Balance Mechanism Design of Single Cylinder Engine Based on Continuous Mass Distribution of Connecting Rod

For a single cylinder engine, the total unbalanced inertial forces occur in the engine block, which results in engine’s vibration and deteriorated noise. In order to eliminate the unbalanced forces, counterweight and primary balance shaft should be attached to the cylinder block so that engine durability and ride comfortability may be further improved. Traditionally one third of connecting rod assembly’s mass is treated as reciprocating mass, and two thirds as rotating mass when designing balance mechanism. In this paper, a new method based on the multibody dynamics simulation is introduced to separate the reciprocating mass and rotating mass of connecting rod assembly. The model consists of crankshaft, connecting rod, piston and the simulation is performed subsequently. According to the simulation results of the main bearing loads, the reciprocating mass and rotating mass are separated. Finally a new balance mechanism is designed and simulation results show that it completely balances inertial forces to improve the engine’s noise vibration and harshness performance.

Adaptive control of track tension estimation using radial basis function neural network

Track tension is a major factor influencing the reliability of a track.In order to reduce the risk of track peel-off,it is necessary to keep track tension constant.However,it is difficult to measure the dynamic tension during off-road operation.Based on the analysis of the relation and external forces depending on free body diagrams of the idler,idler arm,road wheel and road arm,a theoretical estimation model of track tension is built.Comparing estimation results with multibody dynamics simulation results,the rationality of track tension monitor is validated.By the aid of this monitor,a track tension control system is designed,which includes a self-tuning proportional-integral-derivative(PID)controller based on radial basis function neural network,an electro-hydraulic servo system and an idler arm.The tightness of track can be adjusted by turning the idler arm.Simulation results of the vehicle starting process indicate that the controller can reach different expected tensions quickly and accurately.Compared with a traditional PID controller,the proposed controller has a stronger anti-disturbance ability by amending control parameters online.

A recursive formulation based on corotational frame for flexible planar beams with large displacement

A forward recursive formulation based on corotational frame is proposed for flexible planar beams with large displacement.The traditional recursive formulation has been successfully used for flexible mutibody dynamics to improve the computational efficiency based on floating frame,in which the assumption of small strain and deflection is adopted.The proposed recursive formulation could be used for large displacement problems based on the corotational frame.It means that the recursive scheme is used not only for adjacent bodies but also for adjacent beam elements.The nodal relative rotation coordinates of the planar beam are used to obtain equations with minimal generalized coordinates in present formulation.The proposed formulation is different from absolute nodal coordinate formulation and the geometrically exact beam formulation in which the absolute coordinates are used.The recursive scheme and minimal set of dynamic equations lead to a high computational efficiency in numerical integration.Numerical examples are carried out to demonstrate the accuracy and validity of this formulation.For all of the examples,the results of the present formulation are in good agreement with results obtained using commercial software and the published results.Moreover,it is shown that the present formulation is more efficient than the formulation in ANSYS based on GEBF.

Modeling joints with clearance and friction in multibody dynamic simulation of automotive differentials

Defects in kinematic joints can sometimes highly influence the simulation response of the whole multibody system within which these joints are included. For instance, the clearance, the friction, the lubrication and the flexibility affect the transient behaviour, reduce the component life and produce noise and vibration for classical joints such as prismatics, cylindrics or universal joints.In this work, a new 3D cylindrical joint model which accounts for the clearance, the misalignment and the friction is presented. This formulation has been used to represent the link between the planet gears and the planet carrier in an automotive differential model. C 2013 The Chinese Society of Theoretical and Applied Mechanics. [doi：10.1063/2.1301303]

A discrete model of a rope with bending stiffness or viscous damping

A discrete model of a rope is developed and used to simulate the plane motion of the rope fixed at one end.Actually,two systems are presented,whose members are rigid but non-ideal joints involve elasticity or dissipation.The dissipation is reflected simply by viscous damping model, whereas the bending stiffness conception is based on the classical curvature-bending moment relationship for beams and simple geometrical formulas.Equations of motion are derived and their complexity is discussed from the computational point of view.Since modified extended backward differentiation formulas（MEBDF）of Cash are implemented to solve the resulting initial value problems,the technique scheme is outlined.Numerical experiments are performed and influences of the elasticity and damping on behaviour of the model are analyzed.Basic energy principles are used to verify the obtained results.

A new subregion mesh method for the investigation of the elastic-plastic impact in flexible multibody systems

Impact processes between flexible bodies often lead to local stress concentration and wave propagation of high frequency. Therefore, the modeling of flexible multibody systems involving impact should consider the local plastic deformation and the strict requirements of the spatial discretization. Owing to the nonlinearity of the stiffness matrix, the reduction of the element number is extremely important. For the contact-impact problem, since different regions have different requirements regarding the element size, a new subregion mesh method is proposed to reduce the number of the unnecessary elements. A dynamic model for flexible multibody systems with elastic-plastic contact impact is established based on a floating frame of reference formulation and complete Lagrange incremental nonlinear finite-element method to investigate the effect of the elastic-plastic deformation as well as spatial discretization. Experiments on the impact between two bodies are carried out to validate the correctness of the elastic-plastic model. The proposed formulation is applied to a slider-crank system with elastic-plastic impact.

Dynamic Response of Functional Virtual Prototyping of Motorcycle under Combined Excitations

This paper developed 3D product models of motorcycle and engine by UGNX as well as virtual prototyping by ADAMS program with road roughness generated by MATLAB. Under the straight-line running condition, the dynamic responses of motorcycle multibody system to both road and engine excitations were compared with those to only road excitation in terms of vertical acceleration response, amplitude frequency response and power spectral density. The comparisons of simulation data showed that the response due to flat road excitation was around 20 Hz, while that to the combined excitations was in a wide frequency band, of which the major components focused on 10 Hz, 15 Hz, 35 Hz ,70 Hz, 100 Hz and even higher frequencies, reflecting the characteristics of engine excitation based on its unbalanced inertia force and torque. It is concluded that the high fidelity virtual prototyping can simulate the dynamics of motorcycle product well in investigating the vibration and ride comfort performance.