Two dimensional wave front shape induced in a homogeneously strained elastic body by a point perturbing body force.
Two dimensional wavefront shape induced in finitely strained elastic body by impulsive point body force
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Two dimensional wavefront shape induced in finitely strained elastic body by impulsive point body force
A simple method is developed for solving plane-plastic-stress problems with axial symmetry in the strain-hardening range which is based on the deformation theory of plasticity employing the finite-strain concept. The equations defining the problems are first reduced to two simultaneous nonlinear differential equations involving two dependent variables: (a) the octahedral shear strain, and (b) a parameter indicating the ratio of principal stresses. By multiplying the load and dividing the radius by an arbitrary constant, it is possible to solve these problems without iteration for any value of the modified load. The constant is determined by the boundary condition. This method is applied to a circular membrane under pressure, a rotating disk without and with a central hole, and an infinite plate with a circular hole. Two materials, inconel x and 16-25-6, the octahedral shear stress-strain relations of which do not follow the power law, are used. Distributions of octahedral shear strain, as well as of principal stresses and strains, are obtained. These results are compared with the results of the same problems in the elastic range.
An elastic-plastic (incremental and small strain) finite element analysis was used with a crack growth criterion to study crack initiation, stable crack growth, and instability under monotonic loading to failure of metallic materials. The crack growth criterion was a critical crack-tip-opening displacement (CTOD) at a specified distance from the crack tip, or equivalently, a critical crack-tip-opening angle (CTOA). Whenever the CTOD (or CTOA) equaled or exceeded a critical value, the crack was assumed to grow. Single values of critical CTOD were used in the analysis to model crack initiation, stable crack growth, and instability for 7075-T651 and 2024-T351 aluminum alloy compact specimens. Calculated and experimentally measured CTOD values at initiation agreed well for both aluminum alloys. These critical CTOD values were also used to predict failure loads on center-crack tension specimens and a specially-designed three-hole-crack tension specimen made of the two aluminum alloys and of 304 stainless steel. All specimens were 12.7 mm thick. Predicted failure loads for 7075-T651 aluminum alloy and 304 stainless steel specimens were generally within + or - 15 percent of experimental failure loads, whereas the predicted failure loads for 2024-T351 aluminum alloy specimens were generally within + or - 5 percent of the experimental loads.
Higher order correction terms for the stiffness and inertia matrices associated with a triangular plane stress-strain finite dynamic element are developed in detail. Numerical results are presented which indicate that the adoption of these matrices along with a suitable quadratic matrix eigenproblem solver effects a significant economy in the free vibration solution of structure when compared with the analysis based on the usual finite element procedure. Finally, a FORTRAN IV computer program listing of the various relevant matrices is given.
In order to explain the results of the testing conducted in Part I and relate the subscale test specimen performance to case performance during prelaunch loading, a comprehensive analysis program was conducted. Finite element models of both full-scale case and subscale test coupons were made. Validated axisymmetric models were used for the case analyses. Modeling of coupons was accomplished using a generalized plane strain finite element model in most cases, although some three-dimensional finite element modeling was used to investigate the importance of three-dimensional effects on the coupon performance. The large number of plies and complexity of the test coupons used limited the finite element grids to a through-the-thickness discretization of one element per ply. Delaminations observed in the helical ply dropoff region were modeled using frictionless sliders. Comparisons of predicted strains and displacements were made with the test data from both cases and coupons to validate the finite element models. Based on observations of the test results and finite element predictions, a failure criterion for both coupons and cases was postulated and then validated by comparing performance predictions to test results. Finally, the validated failure criterion was used to predict the performance of the first two flight cases under prelaunch loads.
The objective is to study the assumed natural strain (ANS) formulation of finite elements from a variational standpoint. The study is based on two hybrid extensions of the Reissner-type functional that uses strains and displacements as independent fields. One of the forms is a genuine variational principle that contains an independent boundary traction field, whereas the other one represents a restricted variational principle. Two procedures for element level elimination of the strain field are discussed, and one of them is shown to be equivalent to the inclusion of incompatible displacement modes. Also, the 4-node C(exp 0) plate bending quadrilateral element is used to illustrate applications of this theory.
The paper is concerned with the development of constitutive relations for large nonisothermal elastic-viscoplastic deformations for metals. The kinematics of elastic-plastic deformation, valid for finite strains and rotations, is presented. The resulting elastic-plastic uncoupled equations for the deformation rate combined with use of the incremental elasticity law permits a precise and purely deductive development of elastic-viscoplastic theory. It is shown that a phenomenological thermodynamic theory in which the elastic deformation and the temperature are state variables, including few internal variables, can be utilized to construct elastic-viscoplastic constitutive equations, which are appropriate for metals. The limiting case of inviscid plasticity is examined.
An improved hybrid particle-finite element method has been developed for hypervelocity impact simulation. The method combines the general contact-impact capabilities of particle codes with the true Lagrangian kinematics of large strain finite element formulations. Unlike some alternative schemes which couple Lagrangian finite element models with smooth particle hydrodynamics, the present formulation makes no use of slidelines or penalty forces. The method has been implemented in a parallel, three dimensional computer code. Simulations of three dimensional orbital debris impact problems using this parallel hybrid particle-finite element code, show good agreement with experiment and good speedup in parallel computation. The simulations included single and multi-plate shields as well as aluminum and composite shielding materials. at an impact velocity of eleven kilometers per second.
The objective of this paper is to develop a constitutive model for finitely deformed viscoelastic-viscoplastic materials and a micromechanics approach to homogenizing composites consisting of such materials. The development of the constitutive model involves establishing a thermodynamic framework based on finite strain theory, developing a viscoelasticity and a viscoplasticity model based on the thermodynamic framework, developing a radial return algorithm based on a classic framework, and deriving a closed-from incremental constitutive relation in the spatial configuration. The development of the micromechanics approach involves pulling-back the above constitutive relation to the material configuration, formulating a variational statement with the resulting constitutive relation, discretizing variational statement in a finite-dimensional space, and solving the discretized variational statement using an Euler–Newton predictor–corrector method. The constitutive model is calibrated via monotonic uniaxial compression tests on a polymer, and the calibrated model is validated by comparing its predictions with the cyclic test data. It is shown capable of characterizing viscoelasticity, viscoplasticity, and complex loading paths. The micromechanics approach’s capabilities are demonstrated through homogenizing a unidirectional fiber-reinforced composite, subjected to uniaxial, biaxial, and shear loading, at different strain rates. It is demonstrated to be capable of handling rate dependence and complex loading paths. The present framework can be further improved by implementing more sophisticated viscoelasticity and viscoplasticity models in future work.
The accuracy of the finite difference method in the solution of linear elasticity problems that involve either a stress discontinuity or a stress singularity is considered. Solutions to three elasticity problems are discussed in detail: a semi-infinite plane subjected to a uniform load over a portion of its boundary; a bimetallic plate under uniform tensile stress; and a long, midplane symmetric, fiber reinforced laminate subjected to uniform axial strain. Finite difference solutions to the three problems are compared with finite element solutions to corresponding problems. For the first problem a comparison with the exact solution is also made. The finite difference formulations for the three problems are based on second order finite difference formulas that provide for variable spacings in two perpendicular directions. Forward and backward difference formulas are used near boundaries where their use eliminates the need for fictitious grid points.
A consistent step-wise approach is presented to investigate the damage mechanism in composite bonded skin/stringer constructions under uniaxial and biaxial (in-plane/out-of-plane) loading conditions. The approach uses experiments to detect the failure mechanism, computational stress analysis to determine the location of first matrix cracking and computational fracture mechanics to investigate the potential for delamination growth. In a first step, tests were performed on specimens, which consisted of a tapered composite flange, representing a stringer or frame, bonded onto a composite skin. Tests were performed under monotonic loading conditions in tension, three-point bending, and combined tension/bending to evaluate the debonding mechanisms between the skin and the bonded stringer. For combined tension/bending testing, a unique servohydraulic load frame was used that was capable of applying both in-plane tension and out-of-plane bending loads simultaneously. Specimen edges were examined on the microscope to document the damage occurrence and to identify typical damage patterns. For all three load cases, observed failure initiated in the flange, near the flange tip, causing the flange to almost fully debond from skin. In a second step, a two dimensional plane-strain finite element model was developed to analyze the different test cases using a geometrically nonlinear solution. For all three loading conditions, computed principal stresses exceeded the transverse strength of the material in those areas of the flange where the matrix cracks had developed during the tests. In a third step, delaminations of various lengths were simulated in two locations where delaminations were observed during the tests. The analyses showed that at the loads corresponding to matrix ply crack initiation computed strain energy release rates exceeded the values obtained from a mixed mode failure criterion in one location, Hence. Unstable delamination propagation is likely to occur as observed in the experiments.
A consistent step-wise approach is presented to investigate the damage mechanism in composite bonded skin/stringer constructions under uniaxial and biaxial (in-plane/out-of-plane) loading conditions. The approach uses experiments to detect the failure mechanism, computational stress analysis to determine the location of first matrix cracking and computational fracture mechanics to investigate the potential for delamination growth. In a first step, tests were performed on specimens, which consisted of a tapered composite flange, representing a stringer or frame, bonded onto a composite skin. Tests were performed under monotonic loading conditions in tension, three-point bending, and combined tension/bending to evaluate the debonding mechanisms between the skin and the bonded stringer. For combined tension/bending testing, a unique servohydraulic load frame was used that was capable of applying both in-plane tension and out-of-plane bending loads simultaneously. Specimen edges were examined on the microscope to document the damage occurrence and to identify typical damage patterns. For all three load cases, observed failure initiated in the flange, near the flange tip, causing the flange to almost fully debond from the skin. In a second step, a two-dimensional plane-strain finite element model was developed to analyze the different test cases using a geometrically nonlinear solution. For all three loading conditions, computed principal stresses exceeded the transverse strength of the material in those areas of the flange where the matrix cracks had developed during the tests. In a third step, delaminations of various lengths were simulated in two locations where delaminations were observed during the tests. The analyses showed that at the loads corresponding to matrix ply crack initiation computed strain energy release rates exceeded the values obtained from a mixed mode failure criterion in one location. Hence, unstable delamination propagation is likely to occur as observed in the experiments.
The assumed natural strain formulation of finite elements is interpreted from a variational standpoint. The approach is based on hybrid extensions of the Reissner-type functional which uses the strains and displacements as independent fields. Consideration is restricted to linear elasticity. The four-node C(0) plate-bending quadrilateral is used as a specific example to illustrate the application of the present interpretation. A key finding is that any change in the strain-displacement interpolation from the variationally consistent interpolation must be associated in some way to the addition of incompatible displacement modes.
The dynamic analysis of three dimensional elastic beams which experience large rotational and large deformational motions are examined. The beam motion is modeled using an inertial reference for the translational displacements and a body-fixed reference for the rotational quantities. Finite strain rod theories are then defined in conjunction with the beam kinematic description which accounts for the effects of stretching, bending, torsion, and transverse shear deformations. A convected coordinate representation of the Cauchy stress tensor and a conjugate strain definition is introduced to model the beam deformation. To treat the beam dynamics, a two-stage modification of the central difference algorithm is presented to integrate the translational coordinates and the angular velocity vector. The angular orientation is then obtained from the application of an implicit integration algorithm to the Euler parameter/angular velocity kinematical relation. The combined developments of the objective internal force computation with the dynamic solution procedures result in the computational preservation of total energy for undamped systems. The present methodology is also extended to model the dynamics of deployment/retrieval of the flexible members. A moving spatial grid corresponding to the configuration of a deployed rigid beam is employed as a reference for the dynamic variables. A transient integration scheme which accurately accounts for the deforming spatial grid is derived from a space-time finite element discretization of a Hamiltonian variational statement. The computational results of this general deforming finite element beam formulation are compared to reported results for a planar inverse-spaghetti problem.
Finite expansion of spherical cavity in locking medium using finite strain theory
Elastic solids fourth order anharmonic equation of state from finite strain theory, resolving ambiguities by Mie-Gruneisen equation
Theoretical equations of state in geophysics, considering systematics approach to laboratory data, seismic velocity profiles, finite strain and atomistic approach
Wavefront singularities for the displacement functions, associated with the radiation of linear elastic waves from a point source embedded in a finitely strained two-dimensional elastic solid, are examined in detail. It is found that generally the singularities are of order d to the -1/2 power, where d measures distance away from the front. However, in certain exceptional cases singularities of order d to the -n power, where n = 1/4, 2/3, 3/4, may be encountered.