SEARCH · Search NASA
Results for “Finite deformation theory”
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.
Acoustic Radiation Pressure
The theoretical foundation of acoustic radiation pressure in plane wave beams is reexamined. It is shown from finite deformation theory and the Boltzmann-Ehrenfest Adiabatic Principle that the Brillouin stress tensor (BST) is the radiation stress in Lagrangian coordinates (not Eulerian coordinates) and that the terms in the BST are not the momentum flux density and mean excess Eulerian stress but are simply contributions to the variation in the wave oscillation period resulting from changes in path length and true wave velocity, respectively, from virtual variations in the strain. It is shown that the radiation stress in Eulerian coordinates is the mean Cauchy stress (not the momentum flux density, as commonly assumed) and that Langevin's second relation does not yield an assessment of the mean Eulerian pressure, since the enthalpy used in the traditional derivations is a function of the thermodynamic tensions - not the Eulerian pressure. It is shown that the transformation between Lagrangian and Eulerian quantities cannot be obtained from the commonly-used expansion of one of the quantities in terms of the particle displacement, since the expansion provides only the difference between the value of the quantity at two different points in Cartesian space separated by the displacement. The proper transformation is obtained only by employing the transformation coefficients of finite deformation theory, which are defined in terms of the displacement gradients. Finite deformation theory leads to the result that for laterally unconfined, plane waves the Lagrangian and Eulerian radiation pressures are equal with the value (1/4)(2K) along the direction of wave propagation, where (K) is the mean kinetic energy density, and zero in directions normal to the propagation direction. This is contrary to the Langevin result that the Lagrangian radiation pressure in the propagation direction is equal to (2K) and the BST result that the Eulerian radiation pressure in that direction is the momentum flux density.
Finite deformation implementation of a mixed-mode single-integral type cohesive zone with reorienting surfaces of separation
To model material ductile failure and crack propagation, cohesive zone elements can be embedded along potential fracture paths in a finite element simulation. When damage criteria are met, elements in the mesh decohere, simulating the formation and propagation of a crack. In this paper, we present a novel computational algorithm based on finite deformation theory, essential to modeling crack initiation and growth in solids undergoing large deformations. This new algorithm was formulated within a Lagrangian frame of reference to extend previous cohesive zone algorithms to include modeling crack growth in finite deformation contexts. The local coordinate system, necessary for defining an embedded cohesive zone, is constructed based upon the current configuration and is updated within the nonlinear iteration process, thereby resulting in the convergence of the solution for a growing crack in a large deformation quasi-static setting. The model’s accuracy was demonstrated by comparing finite element model simulation results with the analytic case of a constant surface separation, as shown in the verification examples. The power and efficacy of the algorithm to capture large deformations during crack growth were then demonstrated with a double cantilever beam example case. It indicates that the model can be applied to a variety of physical circumstances for predicting crack initiation and growth with delamination and fracture.
Phase-Field Model Incorporating Large Inelastic Strain with Application to the Oxidation of High-Temperature Coating Systems
Structural metallic and ceramic materials at high temperatures often suffer from creep, fatigue and oxidation, which typically involves local large inelastic strain. The large strain is usually approached by Lagrangian-description-based finite deformation theory. Here we develop an efficient, Eulerian-description-based, Incremental Realization of Inelastic Deformation (IRID) algorithm to deal with large inelastic strain with small elastic strain. We apply this algorithm to simulate high temperature oxidation with severe volume expansion using a phase-field approach. The oxidation and deformation of bond coat can lead to spallation and failure of the thermal barrier coatings (TBC) and environmental barrier coatings (EBC), which is a critical concern in design of these coatings. We simulate the growth of thermally grown oxide (TGO) in these coating systems, with phase transformations and microstructure evolution incorporated. The simulation results are compared with relevant experimental reports in the literature.
Stress analysis of spacecraft parachutes using finite elements and large deformation theory
Spacecraft parachute stress analysis, using finite elements with nonlinear elastic properties to obtain shape and load distribution
Finite elasto-plastic deformation. I - Theory and numerical examples
It is demonstrated that the problem of elasto-plastic finite deformation is governed by a quasi-linear model irrespective of deformation magnitude. This feature follows from the adoption of a rate viewpoint toward finite deformation analysis in an Eulerian reference frame. Objectivity of the formulation is preserved by introduction of a frame-invariant stress rate. Equations for piece-wise linear incremental finite element analysis are developed by application of the Galerkin method to the instantaneously linear governing differential equations of the quasi-linear model. Numerical solution capability has been established for problems of plane strain and plane stress. The accuracy of the numerical analysis is demonstrated by consideration of a number of problems of homogeneous finite deformation admitting comparative analytic solution. It is shown that accurate, objective numerical solutions can be obtained for problems involving dimensional changes of an order of magnitude and rotations of a full radian.
On mixed and displacement finite element models of a refined shear deformation theory for laminated anisotropic plates
An improved plate theory that accounts for the transverse shear deformation is presented, and mixed and displacement finite element models of the theory are developed. The theory is based on an assumed displacement field in which the inplane displacements are expanded in terms of the thickness coordinate up to the cubic term and the transverse deflection is assumed to be independent of the thickness coordinate. The governing equations of motion for the theory are derived from the Hamilton's principle. The theory eliminates the need for shear correction factors because the transverse shear stresses are represented parabolically. A mixed finite element model that uses independent approximations of the displacements and moments, and a displacement model that uses only displacements as degrees of freedom are developed. A comparison of the numerical results for bending with the exact solutions of the new theory and the three-dimensional elasticity theory shows that the present theory (and hence the finite element models) is more accurate than other plate-theories of the same order.
Contact between a rigid cylinder and an orthotropic beam under initial stress
The contact behavior between a smooth rigid cylinder and a simply-supported orthotropic beam under uniaxial initial stresses is studied. The displacements are computed by superposing Mindlin plate solution with the solution obtained from Biot's theory of incremental deformation. Finite Fourier transforms are used in solving the equations. The effects of orthotropy and initial stresses on the contact stress distribution are investigated. An indentation law is established from the numerical results.
Smooth indentation of an initially stressed orthotropic beam
The contact behavior between a smooth rigid cylinder and a simply supported orthotropic beam under uniaxial initial stresses is studied. The displacements are computed by superposing Mindlin plate solution with the solution obtained from Biot's theory of incremental deformation. Finite Fourier transforms are used in solving the equations. A point matching technique is used to compute the contact stresses and the amount of indentation for a given contact length. The effects of orthotropy and initial stresses on the contact stress distribution are investigated. An indentation law is established from the numerical results.
On finite deformation elasto-plasticity
Lee (1969) has proposed a theory based on the decomposition of the total deformation gradient to an elastic and plastic part, and from it has concluded that the additive decomposition of the strain rates holds only approximately. Lubarda and Lee (1981) have declared that Lee's 'exact finite-deformation kinematics shows the almost universal assumption that the total velocity strain or rate of deformation is the sum of elastic and plastic rates to be in error'. Hence questions are raised regarding the validity of essentially all finite deformation elasto-plasticity theories. The present investigation is concerned with these questions. It is shown that the additive decomposition of the strain rates follows from all common finite elasto-plasticity concepts. Lee's theory is examined, and it is shown that this theory also leads to an additive strain rate decomposition, and therefore, his conclusion stems from misinterpretation. It is found that the elastic and the plastic strain rates considered by Lee do not correspond to the same configuration. They are, therefore, not compatible measures.
Postbuckling response and failure prediction of flat rectangular graphite-epoxy plates loaded in axial compression
The objectives of the study are to assess the capability of a first-order shear deformable degenerated shell finite element theory to predict the postbuckling response and failure modes of various graphite-epoxy panels loaded in axial compression. Two panels without holes and one with a hole are studied, and the resulting responses and failure modes correlated well with the experimental results. A progressive damage failure mechanism is applied in the nonlinear analysis, which proved successful in predictions of failure location, mode, and load.
Analysis of delamination related fracture processes in composites
This is a final report that summarizes the results achieved under this grant. The first major accomplishment is the development of the sublaminate modeling approach and shear deformation theory. The sublaminate approach allows the flexibility of considering one ply or groups of plies as a single laminated unit with effective properties. This approach is valid when the characteristic length of the response is small compared to the sublaminate thickness. The sublaminate approach was validated comparing its predictions with a finite element solution. A shear deformation theory represents an optimum compromise between accuracy and computational effort in delamination analysis of laminated composites. This conclusion was reached by applying several theories with increasing level of complexity to the prediction of interlaminar stresses and strain energy release rate in a double cracked-lap-shear configuration.
Exact solution of burst pressure for thick-walled pipes using the flow theory of plasticity
This paper develops an exact solution of burst pressure for defect-free, thick-walled pipes using the flow theory of plasticity in terms of the average shear stress yield criterion (or Zhu-Leis yield criterion). The pipe steel is assumed to obey the power-law strain hardening rule, and large plastic deformation is described by the finite strain theory. On this basis, internal pressure is obtained as a power series function of the effective strains on the inside and outside surfaces of the thick-walled pipe with use of the second Bernoulli numbers. At burst failure, the Zhu-Leis flow solution of burst pressure is determined as a power series solution, and the burst effective strains, the burst effective stresses, and the burst pressure are functions of the diameter ratio (D o / D i ), strain hardening exponent (n), and ultimate tensile strength (UTS). Similarly, the Tresca and von Mises flow solutions of burst pressure are also determined as a power series solution in terms of the Tresca and von Mises yield criteria. A general closed-from exact solution of burst pressure is then proposed for the three yield criteria, and the results showed that the proposed exact solution matches well with the power series solution for each yield criterion. Moreover, the von Mises flow solution is an upper bound prediction, the Tresca flow solution is a lower bound prediction, and the Zhu-Leis flow solution is an intermediate prediction that agrees well with the finite element analysis results of burst pressure for thick-walled pipes. Furthermore, two datasets of full-scale burst tests are then utilized to evaluate and validate the proposed flow solutions of burst pressure for both thin and thick-walled pipes.
Analysis of plane-plastic stress problems with axial symmetry in strain-hardening range
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.
A higher-order theory for geometrically nonlinear analysis of composite laminates
A third-order shear deformation theory of laminated composite plates and shells is developed, the Navier solutions are derived, and its finite element models are developed. The theory allows parabolic description of the transverse shear stresses, and therefore the shear correction factors of the usual shear deformation theory are not required in the present theory. The theory also accounts for the von Karman nonlinear strains. Closed-form solutions of the theory for rectangular cross-ply and angle-ply plates and cross-ply shells are developed. The finite element model is based on independent approximations of the displacements and bending moments (i.e., mixed finite element model), and therefore, only C sup o -approximation is required. The finite element model is used to analyze cross-ply and angle-ply laminated plates and shells for bending and natural vibration. Many of the numerical results presented here should serve as references for future investigations. Three major conclusions resulted from the research: First, for thick laminates, shear deformation theories predict deflections, stresses and vibration frequencies significantly different from those predicted by classical theories. Second, even for thin laminates, shear deformation effects are significant in dynamic and geometrically nonlinear analyses. Third, the present third-order theory is more accurate compared to the classical and firt-order theories in predicting static and dynamic response of laminated plates and shells made of high-modulus composite materials.
Nonlinear finite element models of laminated plates and shells
In this expository paper a first-order shear deformation laminate theory that accounts for the von Karman strains is described and its finite element model is developed. The element also accounts for Hill's elastic-perfectly plastic material model. Sample applications of the present element to linear and nonlinear bending of laminated plates and shells are presented. Numerical results showing the effects of boundary conditions, lamination scheme, and loading on the deflections and stresses are also presented. It is found that the effect of shear deformation and boundary conditions on deflections are significant.
Consistent linearization of micromorphic continuum theories
We outline the procedure of consistent linearization and apply it to the micromorphic, microstretch, and micropolar theories of continua. This yields tractable linear theories for nonlinearly elastic microstructured materials undergoing finite deformations. The results may be readily utilized in computational mechanics, stability and bifurcation analyses, and small-deformation problems in the context of these types of continua. Our results generalize those existing in current literature and facilitate their recovery upon incorporating appropriate kinematic and constitutive assumptions.
Compressive strength of flanges
The maximum compressive stress carried by a hinged flange is computed from a deformation theory of plasticity combined with the theory for finite deflections for this structure. The computed stresses agree well with those found experimentally. Empirical observation indicates that the results will also apply fairly well to the more commonly used flanges which are not hinged.