A Nonlinear Spectral Finite Element Model for Analysis of Wave Propagation in Solid with Internal Friction and Dissipation

A geometrically non-linear Spectral Finite Flement Model (SFEM) including hysteresis, internal friction and viscous dissipation in the material is developed and is used to study non-linear dissipative wave propagation in elementary rod under high amplitude pulse loading. The solution to non-linear dispersive dissipative equation constitutes one of the most difficult problems in contemporary mathematical physics. Although intensive research towards analytical developments are on, a general purpose cumputational discretization technique for complex applications, such as finite element, but with all the features of travelling wave (TW) solutions is not available. The present effort is aimed towards development of such computational framework. Fast Fourier Transform (FFT) is used for transformation between temporal and frequency domain. SFEM for the associated linear system is used as initial state for vector iteration. General purpose procedure involving matrix computation and frequency domain convolution operators are used and implemented in a finite element code. Convergnence of the spectral residual force vector ensures the solution accuracy. Important conclusions are drawn from the numerical simulations. Future course of developments are highlighted.

[1]  D. Flanagan,et al.  PRONTO 3D: A three-dimensional transient solid dynamics program , 1989 .

[2]  R. Taylor,et al.  Lagrange constraints for transient finite element surface contact , 1991 .

[3]  J. H. Biffle,et al.  WULFF: a set of computer programs for the large displacement dynamic response of three-dimensional solids , 1976 .

[4]  R. D. Mindlin,et al.  Influence of rotary inertia and shear on flexural motions of isotropic, elastic plates , 1951 .

[5]  Peter Wriggers,et al.  A note on tangent stiffness for fully nonlinear contact problems , 1985 .

[6]  V. K. Luk,et al.  An analytical model on penetration of eroding long rods into metallic targets , 1991 .

[7]  J. O. Hallquist,et al.  User's manual for DYNA2D: an explicit two-dimensional hydrodynamic finite-element code with interactive rezoning , 1982 .

[8]  J. C. Simo,et al.  Consistent tangent operators for rate-independent elastoplasticity☆ , 1985 .

[9]  K. Metzinger,et al.  Bobbin stresses generated by wire winding , 1991 .

[10]  J. C. Simo,et al.  Algorithms for static and dynamic multiplicative plasticity that preserve the classical return mapping schemes of the infinitesimal theory , 1992 .

[11]  Martin W. Heinstein,et al.  An analysis of smoothed particle hydrodynamics , 1994 .

[12]  A. Maradudin,et al.  An Introduction To Applied Anisotropic Elasticity , 1961 .

[13]  J. M. Crissman,et al.  Creep failure and fracture of polyethylene in uniaxial extension , 1979 .

[14]  J. Dienes On the analysis of rotation and stress rate in deforming bodies , 1979 .

[15]  H. Huang Transient Interaction of Plane Acoustic Waves with a Spherical Elastic Shell , 1969 .

[16]  David J. Benson,et al.  Sliding interfaces with contact-impact in large-scale Lagrangian computations , 1985 .

[17]  L. A. Schoof,et al.  EXODUS II: A finite element data model , 1994 .

[18]  R. Landel,et al.  The Strain‐Energy Function of a Hyperelastic Material in Terms of the Extension Ratios , 1967 .

[19]  Samuel W. Key,et al.  HONDO II: a finite element computer program for the large deformation dynamic response of axisymmetric solids. [HONDO] , 1978 .

[20]  Melanie L. Lenard Convergence conditions for restarted conjugate gradient methods with inaccurate line searches , 1976, Math. Program..

[21]  T. Belytschko,et al.  A uniform strain hexahedron and quadrilateral with orthogonal hourglass control , 1981 .

[22]  W. Benz,et al.  Applications of Smooth Particle Hydrodynamics (SPH) to astrophysical problems , 1987 .

[23]  Ted Belytschko,et al.  Eigenvalues and stable time steps for the bilinear Mindlin plate element , 1985 .

[24]  David J. Benson,et al.  A single surface contact algorithm for the post-buckling analysis of shell structures , 1990 .

[25]  M. M. Carroll,et al.  Foundations of Solid Mechanics , 1985 .

[26]  Thomas L. Geers,et al.  Excitation of a fluid‐filled, submerged spherical shell by a transient acoustic wave , 1993 .

[27]  L. P. Bindeman,et al.  Assumed strain stabilization of the eight node hexahedral element , 1993 .

[28]  B. J. Thorne A damage model for rock fragmentation and comparison of calculations with blasting experiments in granite , 1990 .

[29]  Noboru Kikuchi,et al.  Penalty/finite-element approximations of a class of unilateral problems in linear elasticity , 1981 .

[30]  M. J. Forrestal,et al.  Penetration into soil targets , 1992 .

[31]  R. D. Krieg,et al.  Accuracies of Numerical Solution Methods for the Elastic-Perfectly Plastic Model , 1977 .

[32]  C. M. Stone,et al.  SANCHO: a finite element computer program for the quasistatic, large deformation, inelastic response of two-dimensional solids , 1985 .

[33]  D. Tzou,et al.  A spherical cavity-expansion penetration model for concrete targets , 1997 .

[34]  Anthony Ralston,et al.  Mathematical Methods for Digital Computers , 1960 .

[35]  Dimitri P. Bertsekas,et al.  Constrained Optimization and Lagrange Multiplier Methods , 1982 .

[36]  Bjarne Stroustrup,et al.  C++ Programming Language , 1986, IEEE Softw..

[37]  P. Alart,et al.  A generalized Newton method for contact problems with friction , 1988 .

[38]  R. F. Kulak,et al.  Accurate Numerical Solutions for Elastic-Plastic Models , 1979 .

[39]  M. Horne,et al.  Theory of Plasticity , 1955, Nature.

[40]  Edward L. Wilson,et al.  Numerical methods in finite element analysis , 1976 .

[41]  D. B. Longcope,et al.  Penetration into ductile metal targets with rigid spherical-nose rods , 1995 .

[42]  Jürg Dual,et al.  High-resolution analysis of the complex wave spectrum in a cylindrical shell containing a viscoelastic medium. Part II. Experimental results versus theory , 1997 .

[43]  Mazen R. Tabbara,et al.  Spherical cavity-expansion forcing function in PRONTO 3D for application to penetration problems , 1997 .

[44]  J. Monaghan Why Particle Methods Work , 1982 .

[45]  John E. Dennis,et al.  Numerical methods for unconstrained optimization and nonlinear equations , 1983, Prentice Hall series in computational mathematics.

[46]  A. Malkin,et al.  Some conditions for rupture of polymer liquids in extension , 1997 .

[47]  H. D. Hibbitt,et al.  Finite element analysis of sheet forming processes , 1990 .

[48]  Satya N. Atluri,et al.  Non-linear analysis of wave propagation using transform methods , 1993 .

[49]  P. Wriggers,et al.  FINITE ELEMENT FORMULATION OF LARGE DEFORMATION IMPACT-CONTACT PROBLEMS WITH FRICTION , 1990 .

[50]  J. Monaghan,et al.  Shock simulation by the particle method SPH , 1983 .

[51]  M. J. Forrestal,et al.  Penetration of 6061-T651 Aluminum Targets With Rigid Long Rods , 1988 .

[52]  Thomas J. R. Hughes,et al.  Improved numerical dissipation for time integration algorithms in structural dynamics , 1977 .

[53]  L. E. Malvern Introduction to the mechanics of a continuous medium , 1969 .

[54]  T. L. Warren,et al.  Effects of strain hardening and strain-rate sensitivity on the penetration of aluminum targets with spherical-nosed rods , 1998 .

[55]  P M Campbell,et al.  Some New Algorithms for Boundary Value Problems in Smooth Particle Hydrodynamics , 1989 .

[56]  T. A. Duffey,et al.  Non-normal impact of earth penetrators , 1996 .

[57]  Duan Li,et al.  On Restart Procedures for the Conjugate Gradient Method , 2004, Numerical Algorithms.

[58]  S. J. Plimpton,et al.  Molecular Dynamics Simulations of Short-Range Force Systems on 1024-Node Hypercubes , 1990, Proceedings of the Fifth Distributed Memory Computing Conference, 1990..

[59]  Gregory D. Sjaardema Overview of the Sandia National Laboratories Engineering Analysis Code Access System , 1999 .

[60]  J. Lysmer,et al.  Finite Dynamic Model for Infinite Media , 1969 .

[61]  Jürg Dual,et al.  High-resolution analysis of the complex wave spectrum in a cylindrical shell containing a viscoelastic medium. Part I. Theory and numerical results , 1997 .

[62]  M. E. Kipp,et al.  WONDY V: a one-dimensional finite-difference wave-propagation code , 1982 .

[63]  D. L. Hicks,et al.  Stability analysis of WONDY (a hydrocode based on the artificial viscosity method of von Neumann and Richtmyer) for a special case of Maxwell’s law , 1978 .

[64]  Srinivasan Gopalakrishnan,et al.  A spectral finite element for analysis of wave propagation in uniform composite tubes , 2003 .

[65]  A. P. Gilkey,et al.  EXODUS: A finite element file format for pre- and postprocessing , 1988 .

[66]  J. W. Swegle,et al.  Conservative smoothing versus artificial viscosity , 1994 .

[67]  D. Roy Mahapatra,et al.  A spectral finite element model for analysis of axial–flexural–shear coupled wave propagation in laminated composite beams , 2003 .

[68]  C. Y. Cheung,et al.  Experimental study of pure bitumens in tension, compression, and shear , 1997 .

[69]  L. M. Taylor,et al.  Some computational aspects of large deformation, rate-dependent plasticity problems , 1983 .

[70]  Jerry I. Lin,et al.  Explicit algorithms for the nonlinear dynamics of shells , 1984 .

[71]  R. Hill The mathematical theory of plasticity , 1950 .

[72]  G. Strang,et al.  The solution of nonlinear finite element equations , 1979 .

[73]  Z. Mroz,et al.  Associated and non-associated sliding rules in contact friction problems. , 1978 .

[74]  J. Daniel The Conjugate Gradient Method for Linear and Nonlinear Operator Equations , 1967 .

[75]  T. Hughes,et al.  Finite rotation effects in numerical integration of rate constitutive equations arising in large‐deformation analysis , 1980 .

[76]  Richard W. Connelly,et al.  The Use of Extensional Rheometry to Establish Operating Parameters for Stretching Processes , 1982 .

[77]  J. W. Swegle,et al.  Stabilizing S.P.H. with conservative smoothing , 1994 .

[78]  V. Rich Personal communication , 1989, Nature.

[79]  Mark L. Wilkins,et al.  Impact of cylinders on a rigid boundary , 1973 .

[80]  Amy P. Gilkey,et al.  ALGEBRA — A Program That Algebraically Manipulates the Output of a Finite Element Analysis (EXODUS Version) , 2000 .

[81]  N. Mott,et al.  The theory of indentation and hardness tests , 1945 .

[82]  J. C. Simo,et al.  Algorithmic symmetrization of coulomb frictional problems using augmented lagrangians , 1993 .

[83]  E. Lee,et al.  JWL equation of state coefficients for high explosives , 1973 .

[84]  D. P. Flanagan,et al.  GENESIS finite element mesh file format , 1986 .

[85]  D. B. Jr. Longcope,et al.  Oblique penetration modeling and correlation with field tests into a soil target , 1996 .

[86]  R. Ogden Recent Advances in the Phenomenological Theory of Rubber Elasticity , 1986 .

[87]  Donald C. Bogue,et al.  The extensional and failure properties of polymer melts , 1975 .

[88]  Ted Belytschko,et al.  A three-dimensional impact-penetration algorithm with erosion , 1987 .

[89]  Peter Wriggers,et al.  A note on finite‐element implementation of pressure boundary loading , 1991 .

[90]  H. Parisch A consistent tangent stiffness matrix for three‐dimensional non‐linear contact analysis , 1989 .

[91]  L. Taylor,et al.  A finite element model for the analysis of tailored pulse stimulation of boreholes , 1983 .

[92]  C. M. Reeves,et al.  Function minimization by conjugate gradients , 1964, Comput. J..

[93]  Sjaardema NUMBERS: A collection of utilities for pre- and postprocessing two- and three-dimensional EXODUS finite element models , 1989 .

[94]  W. Wunderlich,et al.  Mechanics of Structures. Variational and Computational Methods , 1993 .

[95]  N. Bloembergen,et al.  Theory and Experimental Results , 1948 .

[96]  Amy Gilkey,et al.  BLOT: A mesh and curve plot program for the output of a finite element analysis , 1989 .

[97]  James L. White,et al.  Investigation of failure during elongational flow of polymer melts , 1977 .

[98]  Lawrence D. Cloutman,et al.  An evaluation of smoothed particle hydrodynamics , 1991 .

[99]  D. Grady,et al.  Strain-rate dependent fracture initiation , 1980 .

[100]  J. Monaghan,et al.  Kernel estimates as a basis for general particle methods in hydrodynamics , 1982 .

[101]  J. Daniel,et al.  A conjugate gradient approach to nonlinear elliptic boundary value problems in irregular regions , 1974 .

[102]  B. M. Dobratz,et al.  LLNL Explosives Handbook, Properties of Chemical Explosives and Explosive Simulants , 1985 .

[103]  Anil Chaudhary,et al.  A solution method for static and dynamic analysis of three-dimensional contact problems with friction , 1986 .

[104]  Francis H. Harlow,et al.  PIC and its progeny , 1987 .

[105]  R. Beyer,et al.  Soviet Physics—JETP , 1960 .

[106]  J. M. Kennedy,et al.  Hourglass control in linear and nonlinear problems , 1983 .

[107]  B. Balachandran,et al.  SPECTRAL ANALYSES OF NON-LINEAR INTERACTIONS , 1996 .

[108]  Sia Nemat-Nasser,et al.  Phenomenological modelling of rate-dependent plasticity for high strain rate problems , 1989 .

[109]  D. Owen,et al.  A combined finite‐discrete element method in transient dynamics of fracturing solids , 1995 .

[110]  M. J. Forrestal,et al.  Penetration of Strain-Hardening Targets With Rigid Spherical-Nose Rods , 1991 .

[111]  Richard Chait,et al.  Factors influencing the strength differential of high strength steels , 1972 .

[112]  Nathan M. Newmark,et al.  A Method of Computation for Structural Dynamics , 1959 .

[113]  J. H. Biffle JAC3D -- A three-dimensional finite element computer program for the nonlinear quasi-static response of solids with the conjugate gradient method; Yucca Mountain Site Characterization Project , 1993 .

[114]  W. D. Wightman Philosophical Transactions of the Royal Society , 1961, Nature.

[115]  B. S. Altman,et al.  An empirical equation for penetration depth of ogive-nose projectiles into concrete targets , 1994 .

[116]  Peter G. Bergmann,et al.  Non-Linear Field Theories , 1949 .

[117]  Joseph J Monaghan,et al.  An introduction to SPH , 1987 .

[118]  R. D. Richtmyer,et al.  Difference methods for initial-value problems , 1959 .

[119]  L. Lucy A numerical approach to the testing of the fission hypothesis. , 1977 .

[120]  Joel S. Kuszmaul A Technique For Predicting Fragmentation And Fragment Sizes Resulting From Rock Blasting , 1987 .

[121]  J. J. Rushchitsky,et al.  Interaction of Waves in Solid Mixtures , 1999 .

[122]  R. D. Richtmyer,et al.  A Method for the Numerical Calculation of Hydrodynamic Shocks , 1950 .

[123]  James F. Doyle,et al.  Wave Propagation in Structures , 1989 .

[124]  T. Belytschko,et al.  Large displacement, transient analysis of space frames , 1977 .

[125]  James L. White,et al.  Experimental study of elongational flow and failure of polymer melts , 1978 .

[126]  A. Douglas,et al.  Finite elements for structural analysis , 1983 .

[127]  M. Kruskal,et al.  New similarity reductions of the Boussinesq equation , 1989 .

[128]  J. Oden,et al.  Contact Problems in Elasticity: A Study of Variational Inequalities and Finite Element Methods , 1987 .

[129]  K. F. Riley,et al.  Mathematical methods for the physical sciences , 1975 .

[130]  J. T. Holden On the finite deflections of thin beams , 1972 .

[131]  Mark O. Neal,et al.  Contact‐impact by the pinball algorithm with penalty and Lagrangian methods , 1991 .

[132]  Tunc Geveci,et al.  Advanced Calculus , 2014, Nature.

[133]  A. Giannakopoulos The return mapping method for the integration of friction constitutive relations , 1989 .

[134]  C. M. Stone,et al.  SANTOS - a two-dimensional finite element program for the quasistatic, large deformation, inelastic response of solids , 1997 .

[135]  B. Budiansky,et al.  Elastic moduli of a cracked solid , 1976 .

[136]  A. H. Marchertas,et al.  Nonlinear Finite-Element Method for Plates and Its Application to Dynamic Response of Reactor Fuel Subassemblies , 1974 .

[137]  R. Glowinski,et al.  Augmented Lagrangian and Operator-Splitting Methods in Nonlinear Mechanics , 1987 .

[138]  J. C. Simo,et al.  A continuum-based finite element formulation for the implicit solution of multibody, large deformation-frictional contact problems , 1993 .

[139]  M. J. Forrestal,et al.  Penetration of concrete targets with ogive-nose steel rods , 1998 .

[140]  J. Monaghan Particle methods for hydrodynamics , 1985 .

[141]  Sjaardema GJOIN: A program for merging two or more GENESIS databases , 1992 .

[142]  Peter A. Clarkson,et al.  Nonclassical symmetry reductions for the Kadomtsev-Petviashvili equation , 1991 .

[143]  A. P. Gilkey,et al.  GEN3D: A GENESIS database 2D to 3D transformation program , 1989 .

[144]  W. C. Mills-Curran,et al.  SUPES: a Software Utilities Package for the Engineering Sciences , 1986 .

[145]  Jerrold E. Marsden,et al.  The Hamiltonian structure of nonlinear elasticity: The material and convective representations of solids, rods, and plates , 1988 .

[146]  P. Perzyna Fundamental Problems in Viscoplasticity , 1966 .

[147]  Martin W. Heinstein,et al.  A general-purpose contact detection algorithm for nonlinear structural analysis codes , 1993 .

[148]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[149]  D. Flanagan,et al.  An accurate numerical algorithm for stress integration with finite rotations , 1987 .

[150]  H. S. Morgan,et al.  A phenomenological constitutive model for low density polyurethane foams , 1987 .

[151]  Pierre E. Dupont,et al.  A WAVE APPROACH TO ESTIMATING FREQUENCY-DEPENDENT DAMPING UNDER TRANSIENT LOADING , 2000 .

[152]  R. G. Whirley,et al.  An assessment of numerical algorithms for plane stress and shell elastoplasticity on supercomputers , 1989 .

[153]  J. J. Anza,et al.  A mixed finite element model for the elastic contact problem , 1989 .

[154]  Larry D. Libersky,et al.  Smooth particle hydrodynamics with strength of materials , 1991 .

[155]  C. Kelley Iterative Methods for Linear and Nonlinear Equations , 1987 .

[156]  Tod A. Laursen,et al.  The convected description in large deformation frictional contact problems , 1994 .

[157]  L. M. Taylor,et al.  Microcrack-induced damage accumulation in brittle rock under dynamic loading , 1986 .

[159]  Clifford Ambrose Truesdell,et al.  A first course in rational continuum mechanics , 1976 .