A comparative analysis of methods for solving equations in the nonlinear elasticity theory

An approach to solving a variational equation used in geometrically and physically nonlinear problems of deformable body mechanics is considered. This approach is based on the continuation of a solution with respect to the loading parameter. Large systems of nonlinear ordinary differential equations arise in such problems. Usually, these systems are solved by the Euler methods. It is proposed to use the Runge-Kutta and multistep methods and to estimate the total computational cost. A dependence of numerical errors on the number of integration steps is obtained. An optimal method for solving nonlinear problems is chosen on the basis of this dependence.     

带孔隙损伤的弹性固体的广义变分原理

应用弹性微结构理论,建立了具广义力场带孔隙损伤线弹性固体的基本模型．应用变积方法,同时分别建立了带孔隙损伤弹性固体四类和两类变量的广义变分原理,这些变分原理对应着带孔隙损伤弹性固体微分方程和初值边值条件．应用弹性微结构理论,建立了带孔隙损伤的弹性Timoshenko 梁的基本方程,得到带孔隙损伤的弹性Timoshenko 梁两类变量的广义变分原理．这些广义变分原理为近似求解带孔隙损伤的弹性问题提供了有效途径．     

Generalized variational principles of the viscoelastic body with voids and their applications

From the Boltzmann‘ s constitutive law of viscoelastic materials and the linear theory of elastic materials with voids, a constitutive model of generalized force fields for viscoelastic solids with voids was given. By using the variational integral method, the convolution-type functional was given and the corresponding generalized variational principles and potential energy principle of viscoelastic solids with voids were presented. It can be shown that the variational principles correspond to the differential equations and theinitial and boundary conditions of viscoelastic body with voids. As an application, a generalized variational principle of viscoelastic Timoshenko beams with damage was obtained which corresponds to the differential equations of generalized motion and the initial and boundary conditions of beams. The variational principles provide a way for solving problems of viscoelastic solids with voids.     

Analysis of right closed thin-walled prismatic shells supported by stringers with geometric and physical nonlinearity taken into account

We consider thin-walled right-angle closed prismatic shells with rigid contour of the transverse cross-section. Such shells underlie the schemes used in the analysis of various thin-walled spatial structures. The use of nonlinear physical and geometric relations in the computations permits numerically obtaining the strength margin of the corresponding structures. In the present paper, we propose methods for obtaining a boundary value problem and analyzing such shells with nonlinear factors taken into account; the problem is presented as a system of linear differential equations with variable coefficients. We show that, within the approach proposed, this boundary value problem has a fixed structure independent of the special form of nonlinearity. The entire variety of problems of static analysis of right-angle prismatic shells with nonlinear factors taken into account can be reduced to solving this boundary value problem. Methods for taking a specific nonlinearity into account are treated as various methods for obtaining expressions for the variable coefficients in the matrices of the boundary value problem. We present methods for solving this boundary value problem numerically; these methods are independent of the specific form of the nonlinearity.     

Computational dynamics of a 3D elastic string pendulum attached to a rigid body and an inertially fixed reel mechanism

A high fidelity model is developed for an elastic string pendulum, one end of which is attached to a rigid body while the other end is attached to an inertially fixed reel mechanism which allows the unstretched length of the string to be dynamically varied. The string is assumed to have distributed mass and elasticity that permits axial deformations. The rigid body is attached to the string at an arbitrary point, and the resulting string pendulum system exhibits nontrivial coupling between the elastic wave propagation in the string and the rigid body dynamics. Variational methods are used to develop coupled ordinary and partial differential equations of motion. Computational methods, referred to as Lie group variational integrators, are then developed, based on a finite element approximation and the use of variational methods in a discrete-time setting to obtain discrete-time equations of motion. This approach preserves the geometry of the configurations, and leads to accurate and efficient algorithms that have guaranteed accuracy properties that make them suitable for many dynamic simulations, especially over long simulation times. Numerical results are presented for typical examples involving a constant length string, string deployment, and string retrieval. These demonstrate the complicated dynamics that arise in a string pendulum from the interaction of the rigid body motion, elastic wave dynamics in the string, and the disturbances introduced by the reeling mechanism. Such interactions are dynamically important in many engineering problems, but tend be obscured in lower fidelity models.     

Construction of optimal laws of variation in the angular momentum vector of a dynamically symmetric rigid body

We consider the problem of construction of optimal laws of variation in the angular momentum vector of a dynamically symmetric rigid body so as to ensure the transition of the rigid body from an arbitrary initial angular position to the required final angular position. For the functionals to be minimized, we use combined performance functionals, one of which characterizes the expenditure of time and of the squared modulus of the angular momentum vector in a given proportion, while the other characterizes the expenditure of time and momentum of the modulus of the angular momentum vector necessary to change the rigid body orientation. The control (the vector of the rigid body angular momentum) is assumed to be bounded in the modulus. The problem is solved by using Pontryagin’s maximum principle and the quaternion differential equation [1, 2] relating the vector of the dynamically symmetric rigid body angular momentum to the quaternion of orientation of the coordinate system rotating with respect to the rigid body about its dynamical symmetry axis at an angular velocity proportional to the angular momentum vector projection on the axis. The use of such a model of rotational motion leads to the problem of optimal control with the moving right end of the trajectory and significantly simplifies the analytic study of the problem of construction of optimal laws of variation in the angular momentum vector, because this model explicitly exploits the body angular momentum quaternion (control) instead of the rigid body absolute angular velocity quaternion. We construct general analytic solutions of the differential equations for the boundary-value problems which form systems of nine nonlinear differential equations. It is shown that the process of solving the differential boundary-value problems is reduced to solving two scalar algebraic transcendental equations.     

Experimental determination of the hydrodynamic coefficients for a cylindrical cavity with radial baffles

Some results are presented of an experimental determination of the hydrodynamic coefficients of the equations of the disturbed motion of a rigid body having a cavity in the form of a right circular cylinder with a flat bottom, which is partially filled with a liquid, for the case of radial baffles.     

The pendulum-slosh problem: Simulation using a time-dependent conformal mapping

Suspending a rectangular vessel which is partially filled with fluid from a single rigid pivoting pole produces an interesting theoretical model with which to investigate the dynamic coupling between fluid motion and vessel rotation. The exact equations for this coupled system are derived with the fluid motion governed by the Euler equations relative to the moving frame of the vessel, and the vessel motion governed by a modified forced pendulum equation. The nonlinear equations of motion for the fluid are solved numerically via a time-dependent conformal mapping, which maps the physical domain to a rectangle in the computational domain with a time dependent conformal modulus. The numerical scheme expresses the implicit free-surface boundary conditions as two explicit partial differential equations which are then solved via a pseudo-spectral method in space. The coupled system is integrated in time with a fourth-order Runge–Kutta method. The starting point for the simulations is the linear neutral stability contour discovered by Turner et al. (2015, Journal of Fluid & Structures 52, 166–180). Near the contour the nonlinear results confirm the instability boundary, and far from the neutral curve (parameterized by longer pole lengths) nonlinearity is found to significantly alter the vessel response. Results are also presented for an initial condition given by a superposition of two sloshing modes with approximately the same frequency from the linear characteristic equation. In this case the fluid initial conditions generate large nonlinear vessel motions, which may have implications for systems designed to oscillate in a confined space or on the slosh-induced-rolling of a ship.     

Studies of nonlinear deformation and stability of oval and elliptic cylindrical shells under axial compression

The stability of noncircular shells, in contrast to that of circular ones, has not been studied sufficiently well yet. The publications about circular shells are counted by thousands, but there are only several dozens of papers dealing with noncircular shells. This can be explained on the one hand by the fact that such shells are less used in practice and on the other hand by the difficulties encountered when solving problems involving a nonconstant curvature radius, which results in the appearance of variable coefficients in the stability equations. The well-known solutions of stability problems were obtained by analytic methods and, as a rule, in the linear approximation without taking into account the moments and nonlinearity of the shell precritical state, i.e., in the classical approximation. Here we use the finite element method in displacements to solve the problem of geometrically nonlinear deformation and stability of cylindrical shells with noncircular contour of the transverse cross-section. We use quadrilateral finite elements of shells of natural curvature. In the approximations to the element displacements, we explicitly distinguish the displacements of elements as rigid bodies. We use the Lagrange variational principle to obtain a nonlinear system of algebraic equations for determining the unknown nodal finite elements. We solve the system by a step method with respect to the load using the Newton-Kantorovich linearization at each step. The linear systems are solved by the Kraut method. The critical loads are determined with the use of the Silvester stability criterion when solving the nonlinear problem. We develop an algorithm for solving the problem numerically on personal computers. We also study the nonlinear deformation and stability of shells with oval and elliptic transverse cross-section in a wide range of variations in the ovalization and ellipticity parameters. We find the critical loads and the shell buckling modes. We also examine how the critical loads are affected by the strain nonlinearity and the ovalization and ellipticity of shells.     

On electromagnetic forming processes in finitely strained solids: Theory and examples

The process of electromagnetic forming (EMF) is a high velocity manufacturing technique that uses electromagnetic (Lorentz) body forces to shape sheet metal parts. EMF holds several advantages over conventional forming techniques: speed, repeatability, one-sided tooling, and most importantly considerable ductility increase in several metals. Current modeling techniques for EMF processes are not based on coupled variational principles to simultaneously account for electromagnetic and mechanical effects. Typically, separate solutions to the electromagnetic (Maxwell) and motion (Newton) equations are combined in staggered or lock-step methods, sequentially solving the mechanical and electromagnetic problems.The present work addresses these issues by introducing a fully coupled Lagrangian (reference configuration) least-action variational principle, involving magnetic flux and electric potentials and the displacement field as independent variables. The corresponding Euler-Lagrange equations are Maxwell's and Newton's equations in the reference configuration, which are shown to coincide with their current configuration counterparts obtained independently by a direct approach. The general theory is subsequently simplified for EMF processes by considering the eddy current approximation.Next, an application is presented for axisymmetric EMF problems. It is shown that the proposed variational principle forms the basis of a variational integration numerical scheme that provides an efficient staggered solution algorithm. As an illustration a number of such processes are simulated, inspired by recent experiments of freely expanding uncoated and polyurea-coated aluminum tubes.     

Solution of the optimal turn problem for a spherically symmetric rigid body with arbitrary boundary conditions in the class of generalized conical motions

We consider the problem of time- and energy consumption-optimal turn of a rigid body with spherical mass distribution under arbitrary boundary conditions on the angular position and angular velocity of the rigid body. The optimal turn problem is modified in the class of generalized conical motions, which allows one to obtain closed-form solutions for equations of motion with arbitrary constants. Thus, solving the optimal control boundary value problem is reduced to solving a system of nonlinear algebraic equations for the constants. Numerical examples are considered to illustrate the proximity between the solutions of the traditional and modified problems of optimal turn of a rigid body.     

Dynamics of a rigid cylinder near a plane boundary in the radiation field of an acoustic wave

An approach is described for investigation of the interaction between a rigid body and a viscous fluid boundary under acoustic wave propagation. The influence of the liquid on the rigid body is determined as a mean force, which is a constant in the time component of the hydrodynamic force. This enables the use of a previously developed technique for calculation of pressure in a compressible viscous liquid. The technique takes into account the second-order terms with respect to the wave field parameters and is based on investigation of a system of initially nonlinear hydromechanics equations that can be simplified with respect to the wave motion parameters of the liquid. It has proven possible to retain the second-order terms for determination of stresses in the liquid without having to solve the system of nonlinear equations. The stresses can be expressed in terms of parameters found in the solution of the linearized equations of the compressible viscous liquid. In this way, the solution of linearized equations is expressed in terms of a scalar and vector potentials. The problem statement is derived for a rigid cylinder located near a rigid flat wall under the effects of a wave propagating perpendicular to the wall. The solution for this particular example is obtained.     

Impact of a spherical rigid body on the surface of a cavity in a compressible liquid: An axisymmetric problem

The shock interaction of a spherical rigid body with a spherical cavity is studied. This nonstationary mixed boundary-value problem with an unknown boundary is reduced to an infinite system of linear Volterra equations of the second kind and the differential equation of motion of the body. The hydrodynamic and kinematic characteristics of the process are obtained __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 1, pp. 11–19, January 2008.     

Instability of sloshing motion in a vessel undergoing pivoted oscillations

Suspending a rectangular vessel partially filled with an inviscid fluid from a single rigid pivoting rod produces an interesting physical model for investigating the dynamic coupling between the fluid and vessel motion. The fluid motion is governed by the Euler equations relative to the moving frame of the vessel, and the vessel motion is given by a modified forced pendulum equation. The fully nonlinear, two-dimensional, equations of motion are derived and linearised for small-amplitude vessel and free-surface motions, and the natural frequencies of the system analysed. It is found that the linear problem exhibits an unstable solution if the rod length is shorter than a critical length which depends on the length of the vessel, the fluid height and the ratio of the fluid and vessel masses. In addition, we identify the existence of 1:1 resonances in the system where the symmetric sloshing modes oscillate with the same frequency as the coupled fluid/vessel motion. The implications of instability and resonance on the nonlinear problem are also briefly discussed.     

Oscillations of a vessel containing a viscous fluid

The oscillations of a rigid body having a cavity partially filled with an ideal fluid have been studied in numerous reports, for example, [1–6]. Certain analogous problems in the case of a viscous fluid for particular shapes of the cavity were considered in [6, 7]. The general equations of motion of a rigid body having a cavity partially filled with a viscous liquid were derived in [8]. These equations were obtained for a cavity of arbitrary form under the following assumptions: 1) the body and the liquid perform small oscillations (linear approximation applicable); 2) the Reynolds number is large (viscosity is small). In the case of an ideal liquid the equations of [8] become the previously known equations of [2–6]. In the present paper, on the basis of the equations of [8], we study the free and the forced oscillations of a body with a cavity (vessel) which is partially filled with a viscous liquid. For simplicity we consider translational oscillations of a body with a liquid, since even in this case the characteristic mechanical properties of the system resulting from the viscosity of the liquid and the presence of a free surface manifest themselves.The solutions are obtained for a cavity of arbitrary shape. We then consider some specific cavity shapes.     

大型偏心结构吊装欠约束问题研究

吊装施工过程中被吊模块的水平度是作业要求的重要指标,通常需要增加配重调平。传统有限元方法需要补充约束以消除单元刚体位移,且需要重复计算平衡方程来求解调平载荷,效率不高。将模块的运动分解为随动坐标系的整体运动以及相对该坐标系的弹性变形,可将欠约束问题化为多体系统的静平衡问题。基于虚功率原理推导了吊装平顺时刻的节点力平衡方程以及相应的切线刚度矩阵,并将配重表示为基础配重与载荷系数相乘的形式。通过对节点力平衡方程求导,得到一组以载荷系数为自变量的微分方程,通过求解微分方程并结合水平度判据,可快速搜寻满足水平度要求的载荷系数。数值算例表明,该方法在解决偏心模块吊装欠约束问题方面具有明显的优势,在确定配重载荷方面具有较快的速度和合理的精度。     

Off-center impact of an elastic column by a rigid mass

Based on the Timoshenko beam model the equations of motion are obtained for large deflection of off-center impact of a column by a rigid mass via Hamilton's principle. These are a set of coupled nonlinear partial differential equations. The Newmark time integration scheme and differential quadrature method are employed to convert the equations into a set of nonlinear algebraic equations for displacement components. The equations are solved numerically and the effects of weight and velocity of the rigid mass and also off-center distance on deformation of the column are studied.     

A Variational Approach to Dynamics of Flexible Multibody Systems

Abstract

This paper presents a variational formulation of constrained dynamics of flexible multibody systems, using a vector-variational calculus approach. Body reference frames are used to define global position and orientation of individual bodies in the system, located and oriented by position of its origin and Euler parameters, respectively. Small strain linear elastic deformation of individual components, relative to their body reference frames, is defined by linear combinations of deformation modes that are induced by constraint reaction forces and normal modes of vibration. A library of kinematic couplings between flexible and/or rigid bodies is defined and analyzed. Variational equations of motion for multibody systems are obtained and reduced to mixed differential-algebraic equations of motion. A space structure that must deform during deployment is analyzed, to illustrate use of the methods developed     

Dynamics of Slider-Crank Mechanisms with Flexible Supports and Non-Ideal Forcing

Transient and steady state dynamic response of a class of slider-crank mechanisms is investigated. Specifically, the class of mechanisms examined involves rigid members but compliant supporting bearings. Moreover, the mechanisms are subjected to non-ideal forcing. Namely, both the driving and the resisting loads are expressed as a function of the angular coordinate describing the crank rotation. First, an appropriate set of equations of motion is derived by applying Lagrange's equations. These equations are strongly nonlinear due to the large rigid body rotation of the crank and the connecting rod, as well as due to the nonlinearities associated with the bearing action and the form of the driving and the resisting loads. Consequently, the dynamics of the resulting dynamical system is examined by solving the equations of motion numerically. More specifically, transient response is captured by direct integration, while determination of complete branches of steady state response is achieved by applying appropriate numerical methodologies. Initially, mechanisms whose crankshaft is supported by bearings with rolling elements and linear stiffness characteristics are examined. Then, numerical results are presented for rolling element bearings with nonlinear stiffness characteristics. Finally, the study is focused on mechanisms supported by hydrodynamic bearings. In all cases, the attention is focused on investigating the influence of the system parameters on its dynamics. Moreover, models with constant crank angular velocity are first analysed, since they provide valuable insight into some aspects of the system dynamics. Eventually, the emphasis is shifted to the general case of non-ideal forcing, originating from the dependence of the driving and the resisting moments on the crankshaft motion.