A nonlinear planar beam formulation with stretch and shear deformations under end forces and moments

A new nonlinear planar beam formulation with stretch and shear deformations is developed in this work to study equilibria of a beam under arbitrary end forces and moments. The slope angle and stretch strain of the centroid line, and shear strain of cross-sections, are chosen as dependent variables in this formulation, and end forces and moments can be either prescribed or resultant forces and moments due to constraints. Static equations of equilibria are derived from the principle of virtual work, which consist of one second-order ordinary differential equation and two algebraic equations. These equations are discretized using the finite difference method, and equilibria of the beam can be accurately calculated. For practical, geometrically nonlinear beam problems, stretch and shear strains are usually small, and a good approximate solution of the equations can be derived from the solution of the corresponding Euler–Bernoulli beam problem. The bending deformation of the beam is the only important one in a slender beam, and stretch and shear strains can be derived from it, which give a theoretical validation of the accuracy and applicability of the nonlinear Euler–Bernoulli beam formulation. Relations between end forces and moments and relative displacements of two ends of the beam can be easily calculated. This formulation is powerful in the study of buckling of beams with various boundary conditions under compression, and can be used to calculate post-buckling equilibria of beams. Higher-order buckling modes of a long slender beam that have complex configurations are also studied using this formulation.     

NUMERICAL STUDY OF THE FLOW PAST A CYLINDER EXCITED TRANSVERSELY TO THE INCIDENT STREAM. PART 1: LOCK-IN ZONE,HYDRODYNAMIC FORCES AND WAKE GEOMETRY

The numerical study of the flow past a circular cylinder forced to oscillate transversely to the incident stream is presented herein, at a fixed Reynolds number equal to 106. The finite element technique was favoured for the solution of the Navier–Stokes equations, in the formulation where the stream function and the vorticity are the field variables. The cylinder oscillation frequency ranged between 0·80 and 1·20 of the natural vortex-shedding frequency, and the oscillation amplitude extended up to 50% of the cylinder diameter. Since the resolution of the characteristics of synchronized wakes is the focus of the study, the first task is the determination of the boundary of the lock-in region. The computation revealed that, when the cylinder oscillation frequency exceeds the frequency of the natural shedding of vortices, the flow is not absolutely periodic at subsequent cycles but a quasiperiodic flow pattern occurs, which creates difficulty in the determination of the lock-in boundary. The time histories of the drag and lift forces for various oscillation parameters are presented, while the vorticity contours were favoured for the numerical flow visualization. The hydrodynamic forces, the phase angle between the lift force and the cylinder displacement, and the parameters of the wake geometry when steady state was reached, are presented in cumulative diagrams. These diagrams indicate the effect of the oscillation parameters on the hydrodynamic forces and on the wake geometry.     

Prediction of unwinding behaviors and problems of cables from inner-winding spool dispensers

Wire-guided control technologies are widely used to increase the targeting accuracy of advanced military weapons through the use of unwinding dispensers to guarantee that unwinding occurs without any problems, such as tangling or cutting. In this study, the transient behaviors of cables unwinding from inner-winding cylindrical spool dispensers are investigated. The cable is withdrawn from the spool dispenser at a constant velocity through a fixed point located along the axis of the spool dispenser. And when the cable is flown out of the dispenser, because several dynamic forces such as inertial forces, Coriolis forces, centrifugal forces, tensile forces, and fluid-resistance forces act on the cables, the cables exhibit highly nonlinear and complex unwinding behaviors which are called as unwinding balloons. For predicting these complex unwinding motions, the governing equations of motion for cables unwinding from a cylindrical spool dispenser of the inner-winding type are derived using the extended Hamilton’s principles for an open system in which mass can be transported at each boundary. Modified finite difference methods are used to discretize the spatial variables of the derived nonlinear partial differential equations. The time responses of the unwinding cables are calculated using Newmark time integration methods. Finally, in order to present numerical examples, an inner-winding spool dispenser that can be unwound up to a length of 50 km is designed by using simple geometrical relationships. The behaviors during unwinding from the designed inner spool dispenser are presented and modifications of the spool dispenser for avoiding unwinding problems are proposed by using the suggested numerical methods.     

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.     

Three-Dimensional Large Deformation Analysis of the Multibody Pantograph/Catenary Systems

To accurately model the nonlinear behavior of the pantograph/catenary systems, it is necessary to take into consideration the effect of the large deformation of the catenary and its interaction with the nonlinear pantograph system dynamics. The large deformation of the catenary is modeled in this investigation using the three-dimensional finite element absolute nodal coordinate formulation. To model the interaction between the pantograph and the catenary, a sliding joint that allows for the motion of the pan-head on the catenary cable is formulated. To this end, a non-generalized arc-length parameter is introduced in order to be able to accurately predict the location of the point of contact between the pan-head and the catenary. The resulting system of differential and algebraic equations formulated in terms of reference coordinates, finite element absolute nodal coordinates, and non-generalized arc-length and contact surface parameters are solved using computational multibody system algorithms. A detailed three-dimensional multibody railroad vehicle model is developed to demonstrate the use of the formulation presented in this paper. In this model, the interaction between the wheel and the rail is considered. For future research, a method is proposed to deal with the problem of the loss of contact between the pan-head and the catenary cable.     

Calculation of laminar recirculating flows using a local non-staggered grid refinement system

A grid-embedding technique for the solution of two-dimensional incompressible flows governed by the Navier-Stokes equations is presented. A finite volume method with collocated primitive variables is employed to ensure conservation at the interfaces of embedding grids as well as global conservation. The discretized equations are solved simultaneously for the whole domain, providing a strong coupling between regions of different refinement. The formulation presented herein is applicable to uniform or non-uniform Cartesian meshes. The method was applied to the solution of two scalar transport equations, to cavity flows driven by body and shear forces and to a sudden plane contraction flow. The numerical predictions are compared with the exact solutions when available and with experimental data. The results show that neither the convergence rate nor the stability of the method is affected by the presence of embedded grids. Embedded grids provide a better distribution of grid nodes over the computational domain and consequently the solution accuracy was improved. The grid-embedding technique proved also that significant savings in computing time could be achieved.     

Some Results of the Finite Deflection Analysis of Clamped Skew Sandwich Plates

The suitability of Galerkin's method for the solution of the problem of the finite deflection analysis of clamped skew sandwich plates is studied. The five coupled nonlinear governing differential equations for sandwich plates are transformed into nonlinear algebraic equations by using Galerkin's method of error minimization. These equations are then solved using an iterative algorithm suggested by Brown. Comparisons of the results of the present analysis with available solutions show good agreement. Numerical results are presented for skew sandwich plates for a wide range of values of the core modulus for different skew angles and aspect ratios. Simplicity in formulation and computation is the advantage of the method as compared with other methods of nonlinear analysis. Computing time and memory requirements in a digital computer are relatively very small, which makes the method attractive.     

Nonlinear dynamic response and active vibration control of the viscoelastic cable with small sag

IntroductionCablesareveryefficientstructuralmembersandhencehavebeenwidelyusedinmanylong_spanstructures,includingcable_supportbridges,guyedtowersandcable_supportroofs.Sincecablesarelight,veryflexibleandlightlydamped ,structuresutilizingcables,i.e .,cable_structuresystems,usuallyhavevariousdynamicproblems.Theirmodelsarethereforeverimportantinpredictingandcontrollingtheirresponses.Inthelastdecade,thenonlineardynamicvibrationandstabilitybehaviorofcablesandcable_structureshavedrawntheattentionofman…     

Control structure interaction in the nonlinear analysis of flexible mechanical systems

The effect of the control structure interaction on the feedforward control law as well as the dynamics of flexible mechanical systems is examined in this investigation. An inverse dynamics procedure is developed for the analysis of the dynamic motion of interconnected rigid and flexible bodies. This method is used to examine the effect of the elastic deformation on the driving forces in flexible mechanical systems. The driving forces are expressed in terms of the specified motion trajectories and the deformations of the elastic members. The system equations of motion are formulated using Lagrange's equation. A finite element discretization of the flexible bodies is used to define the deformation degrees of freedom. The algebraic constraint equations that describe the motion trajectories and joint constraints between adjacent bodies are adjoined to the system differential equations of motion using the vector of Lagrange multipliers. A unique displacement field is then identified by imposing an appropriate set of reference conditions. The effect of the nonlinear centrifugal and Coriolis forces that depend on the body displacements and velocities are taken into consideration. A direct numerical integration method coupled with a Newton-Raphson algorithm is used to solve the resulting nonlinear differential and algebraic equations of motion. The formulation obtained for the flexible mechanical system is compared with the rigid body dynamic formulation. The effect of the sampling time, number of vibration modes, the viscous damping, and the selection of the constrained modes are examined. The results presented in this numerical study demonstrate that the use of the driving forees obtained using the rigid body analysis can lead to a significant error when these forces are used as the feedforward control law for the flexible mechanical system. The analysis presented in this investigation differs significantly from previously published work in many ways. It includes the effect of the structural flexibility on the centrifugal and Coriolis forces, it accounts for all inertia nonlinearities resulting from the coupling between the rigid body and elastic displacements, it uses a precise definition of the equipollent systems of forces in flexible body dynamics, it demonstrates the use of general purpose multibody computer codes in the feedforward control of flexible mechanical systems, and it demonstrates numerically the effect of the selected set of constrained modes on the feedforward control law.     

Nonlinear axisymmetric deformations of an elastic tube under external pressure

The problem of the finite axisymmetric deformation of a thick-walled circular cylindrical elastic tube subject to pressure on its external lateral boundaries and zero displacement on its ends is formulated for an incompressible isotropic neo-Hookean material. The formulation is fully nonlinear and can accommodate large strains and large displacements. The governing system of nonlinear partial differential equations is derived and then solved numerically using the C++ based object-oriented finite element library Libmesh. The weighted residual-Galerkin method and the Newton-Krylov nonlinear solver are adopted for solving the governing equations. Since the nonlinear problem is highly sensitive to small changes in the numerical scheme, convergence was obtained only when the analytical Jacobian matrix was used. A Lagrangian mesh is used to discretize the governing partial differential equations. Results are presented for different parameters, such as wall thickness and aspect ratio, and comparison is made with the corresponding linear elasticity formulation of the problem, the results of which agree with those of the nonlinear formulation only for small external pressure. Not surprisingly, the nonlinear results depart significantly from the linear ones for larger values of the pressure and when the strains in the tube wall become large. Typical nonlinear characteristics exhibited are the “corner bulging” of short tubes, and multiple modes of deformation for longer tubes.     

On the dynamics of rotating circular cable structures

The combined effect of gravity and a centrifugal field on the dynamics of a heavy cable with negligible bending stiffness is treated. The inextensible length of the cable is preassigned in such a way that the cross-sectional shape of the cable structure is circular if the contiguration rotates about its central axis with a constant speed and no external forces act. Under the additional influence of gravity, complicated nonlinear weight-excited vibrations occur.To understand the variety of vibrational phenomena. The dynamic system is studied by some experiments first. In order to explain the experimental results theoretically, the governing nonlinear boundary value problem is derived next. Subsequently, an appropriate variational formulation for application of a Rayleigh Ritz procedure is suggested. The resulting nonlinear ordinary differential equations approximately deseribe a part of the observed vibrational behaviour.Both the experiments and the calculations demonstrate that, for high speeds, nonlinear weight-excited vibrations about the circular reference configuration occur. On the other hand, for low velocitites, periodie and even chaos-like snap-through phenomena appear.     

Stress resultants plasticity with general closest point projection

In this paper, general closest point projection algorithm is derived for the elastoplastic behavior of a cross-section of a beam finite element. For given section deformations, the section forces (stress resultants) and the section tangent stiffness matrix are obtained as the response for the cross-section. Backward Euler time integration rule is used for the solution of the nonlinear evolution equations. The solution yields the general closest projection algorithm for stress resultants plasticity model. Algorithmic consistent tangent stiffness matrix for the section is derived. Numerical verification of the algorithms in a mixed formulation beam finite element proves the accuracy and robustness of the approach in simulating nonlinear behavior.     

斜拉桥拉索-阻尼器系统非线性响应分析

考虑索的抗弯刚度、垂度及几何非线性的影响,得出了索一阻尼器系统的空间非线性振动偏微分方程,用中心差分法将微分方程在空间内离散,导出了系统的非线性振动常微分方程组。结合Newmark法及虚拟力法提出了一种用于求解非线性振动瞬态响应的杂交分析算法。并以典型的斜拉桥拉索为研究对象,给出了数值算例,并与Runge—Kutta直接积分法进行了比较,说明了杂交算法的准确性及有效性。     

Large Amplitude Three-Dimensional Free Vibrations of Inclined Sagged Elastic Cables

The nonlinear characteristics in the large amplitude three-dimensionalfree vibrations of inclined sagged elastic cables are investigated. Amodel formulation which is not limited to cables having smallsag-to-span ratios and takes into account the axial deformation effectis considered. Based on a multi-degree-of-freedom cable model, a finitedifference discretization is employed within a numerical solution of thegoverning equations of three-dimensional coupled motion. Variousnumerical examples of arbitrarily inclined sagged cables with initialout-of-plane or in-plane motions are carried out for the case of aspecified end tension. The major findings consist of highlighting theextent of two-and three-dimensional nonlinear couplings, the occurrenceof nonlinear dynamic tensions, and the meaningfulness of modaltransition phenomena ensuing from the activation of various internalresonance conditions. The influence of cable inclination on thenonlinear dynamic behavior is also evaluated. Comprehensive discussionand comparison of large amplitude free vibrations of horizontal andinclined sagged cables are presented.     

Absolute nodal coordinate formulation of large-deformation piezoelectric laminated plates

The Absolute Nodal Coordinate Formulation (ANCF) has been initiated in 1996 by Shabana (Computational Continuum Mechanics, 3rd edn., Cambridge: Cambridge University Press, 2008). It introduces large displacements of planar and spatial finite elements relative to the global reference frame without using any local frame. A sub-family of beam, plate and cable finite elements with large deformations are proposed and employed the 3D theory of continuum mechanics. In the ANCF, the nodal coordinates consist of absolute position coordinates and gradients that can be used to define a unique rotation and deformation fields within the element. In contrast to other large deformation formulations, the equations of motion contain constant mass matrices as well as zero centrifugal and Coriolis inertia forces. The only nonlinear term is a vector of elastic forces. This investigation concerns a way to generate new finite element in the ANCF for laminated composite plates. This formulation utilizes the assumption that the bonds between the laminae are thin and shear is non-deformable. Consequently, the Equivalent Single Layer, ESL model, is implemented. In the ESL models, the laminate is assumed to deform as a single layer, assuming a smooth variation of the displacement field across the thickness. In this paper, the coupled electromechanical effect of Piezoelectric Laminated Plate is imposed within the ANCF thin plate element, in such a way as to achieve the continuity of the gradients at the nodal points, and obtain a formulation that automatically satisfies the principle of work and energy. Convergence and accuracy of the finite-element ANCF Piezoelectric Laminated Plate is demonstrated in geometrically nonlinear static and dynamic test problems, as well as in linear analysis of natural frequencies. The computer implementation and several numerical examples are presented in order to demonstrate the use of the formulation developed in this paper. A comparison with the commercial finite element package COMSOL MULTIPHYSICS () is carried out with an excellent agreement.     

A theory of finite elastic consolidation

Presented in this paper is a general theory describing the consolidation of a porous elastic soil. The formulation allows for the occurrence of finite geometry changes and finite elastic strains during the consolidation process. The governing equations have been cast in a rate form and the laws which determine deformation and pore fluid flow, i.e. Hooke's law and Darcy's law, are presented in a frame indifferent manner. A numerical technique is described that provides an approximate solution to the governing equations. The theory and the solution technique are illustrated by several examples of practical interest.     

Hydrodynamic and Geometric Stiffening Effects on the Out-of-Plane Waves of Submerged Cables

This study focuses on the relative importance of two sources of nonlinearities affecting submerged cable response. The first of these is the added fluid damping offered by the surrounding medium while the second is the geometric stiffening offered by the cable through finite extensions of its centerline. The contribution of each nonlinear effect, taken separately and in tandem, is evaluated herein through the study of structural waves that form in the (out-of-plane) direction normal to the cable equilibrium plane.Numerical solutions are pursued herein using a finite difference algorithm which is brought to bear on two nonlinear cable/fluid models including: (1)~a nonlinear submerged cable model in which hydrodynamic drag is the sole nonlinear mechanism (referred to herein as the 'nonlinear drag model'); and (2)~a nonlinear submerged cable model in which hydrodynamic drag and geometric stiffening are both active nonlinear mechanisms (the 'nonlinear elastic-drag model'). Numerical solutions for propagating cable waves are developed for the case of a long suspension subjected to a concentrated harmonic excitation source. Conclusions are subsequently drawn regarding the spatial decay of the resulting out-of-plane waves and the dynamic cable tension induced by these waves. The effect of these two nonlinear mechanisms is further explored through the analysis of two additional, linear models: (3)~a simple linear taut string model without drag (the 'simple model'); and (4)~a linear taut string model with linear drag (the 'linear drag model'). The results of all models are critically compared and the range of validity of the linear/cable fluid models are assessed.     

On a pure complementary energy principle and a force-based finite element formulation for non-linear elastic cables

A complementary-dual force-based finite element formulation is proposed for the geometrically exact quasi-static analysis of one-dimensional hyperelastic perfectly flexible cables lying in the two-dimensional space. This formulation employs as approximate functions the exact statically admissible force fields, i.e., those that satisfy the equilibrium differential equations in strong form, as well as the equilibrium boundary conditions. The formulation relies on a principle of total complementary energy only expressed in terms of force fields, being therefore called a pure principle. Under the assumption of stress-unilateral behavior, this principle can be regarded as being dual to the principle of minimum total potential energy, corresponding therefore to a maximum principle. Some numerical applications, including cables suspended from two and three points at the same level or at different levels, with both Hookean and Neo-Hookean material behaviors, are presented. As it will be shown, in contrast to the standard two-node displacement-based formulation derived from the principle of minimum total potential energy, the proposed dual force-based formulation is capable of providing the exact solution of a given problem only using a single finite element per cable. Both the proposed principle of pure complementary energy and its corresponding force-based finite element formulation can be easily extended to the case of cables lying in the three-dimensional space.     

Finite Deflection of Skew Sandwich Plates on Elastic Foundations by the Galerkin Method

ABSTRACT

Application of the Galerkin method to various fluid and structural mechanics problems that are governed by a single linear or nonlinear differential equation is well known [1-5]. Recently, the method has been extended to finite element formulations [6-10], In this paper the suitability of the Galerkin method for solution of large deflection problems of plates is studied. The method is first applied to investigate large deflection behavior of clamped isotropic plates on elastic foundations. After validity of the method is established, it is then extended to analyze problems of large deflection of clamped skew sandwich plates, both with and without elastic foundations. The plates are considered to be subjected to uniformly distributed loads. The governing differential equations for the sandwich plate in terms of displacements in Cartesian coordinates are first established and then transformed into skew coordinates. The nonlinear differential equations of the plates are then transformed into nonlinear algebraic equations, using the Galerkin method. These equations are solved using a Newton-Raphson iterative procedure. The parameters considered herein for large deflection behavior of skew sandwich plates are the aspect ratio of the plate, Poisson's ratio, skew angle, shearing stiffnesses of the core, and foundation moduli. Numerical results are presented for skew sandwich plates for various skew angles and aspect ratios. Simplicity and quick convergence are the advantages of the method, in comparison with other much more laborious numerical methods that require extensive computer facilities.