Search NASA⌕ Search

SEARCH · Search NASA

Results for “large material deformation and rotation”

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.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

HOSS: an implementation of the combined finite-discrete element method

Nearly thirty years since its inception, the combined finite-discrete element method (FDEM) has made remarkable strides in becoming a mainstream analysis tool within the field of Computational Mechanics. FDEM was developed to effectively “bridge the gap” between two disparate Computational Mechanics approaches known as the finite and discrete element meth-ods. At Los Alamos National Laboratory (LANL) researchers developed the Hybrid Optimization Software Suite (HOSS) as a hybrid multi-physics platform, based on FDEM, for the simulation of solid material behavior complemented with the latest technological enhancements for full fluid–solid interaction. Furthermore, in HOSS, several newly developed FDEM algorithms have been implemented that yield more accurate material deformation formulations, inter-particle interaction solvers, and fracture and fragmentation solutions. Additionally, an explicit computational fluid dynamics solver and a novel fluid–solid interaction algorithms have been fully integrated (as opposed to coupled) into the HOSS’ solid mechanical solver, allowing for the study of an even wider range of problems. Advancements such as this are leading HOSS to become a tool of choice for multi-physics problems. Finally, HOSS has been successfully applied by a myriad of researchers for analysis in rock mechanics, oil and gas industries, engineering application (structural, mechanical and biomedical engineering), mining, blast loading, high velocity impact, as well as seismic and acoustic analysis. This paper intends to summarize the latest development and application efforts for HOSS.

58 GEOSCIENCES↗
NASA NTRS1998

Analysis of Curved Sandwich Panels Subjected to Combined Temperature Gradient and Mechanical Loads

The results of a detailed study of the nonlinear response of curved sandwich panels with composite face sheets and subjected to a temperature gradient through-the-thickness combined with mechanical loadings are presented. The analysis is based on a first-order shear-deformation Sanders-Budiansky type theory with the effects of large displacements, moderate rotations, transverse shear deformation and laminated anisotropic material behavior included. A mixed formulation is used with the fundamental unknowns consisting of the generalized displacements and the stress resultants of the panel. The nonlinear displacements, strain energy, principal strains, transverse shear stresses, transverse shear strain energy density, and their hierarchical sensitivity coefficients are evaluated. The hierarchical sensitivity coefficients measure the sensitivity of the nonlinear response to variations in the panel parameters, the effective properties of the face sheet layers and the core, and the micromechanical parameters. Numerical results are presented for cylindrical panels subjected to combined pressure loading, edge shortening or extension, edge shear and a temperature gradient through the thickness. The results show the effects of variations in the loading and the panel aspect ratio, on the nonlinear response and its sensitivity to changes in the various panel, effective layer and micromechanical parameters.

Noor, Ahmed K.↗
NASA NTRS1999

Analysis of Composite Panels Subjected to Thermo-Mechanical Loads

The results of a detailed study of the effect of cutout on the nonlinear response of curved unstiffened panels are presented. The panels are subjected to combined temperature gradient through-the-thickness combined with pressure loading and edge shortening or edge shear. The analysis is based on a first-order, shear deformation, Sanders-Budiansky-type shell theory with the effects of large displacements, moderate rotations, transverse shear deformation, and laminated anisotropic material behavior included. A mixed formulation is used with the fundamental unknowns consisting of the generalized displacements and the stress resultants of the panel. The nonlinear displacements, strain energy, principal strains, transverse shear stresses, transverse shear strain energy density, and their hierarchical sensitivity coefficients are evaluated. The hierarchical sensitivity coefficients measure the sensitivity of the nonlinear response to variations in the panel parameters, as well as in the material properties of the individual layers. Numerical results are presented for cylindrical panels and show the effects of variations in the loading and the size of the cutout on the global and local response quantities as well as their sensitivity to changes in the various panel, layer, and micromechanical parameters.

Noor, Ahmed K.↗
DOE OSTI2022

A finite micro-rotation material point method for micropolar solid and fluid dynamics with three-dimensional evolving contacts and free surfaces

This paper introduces an explicit material point method designed specifically for simulating the micropolar continuum dynamics in the finite deformation and finite microrotation regime. The material point method enables us to simulate large deformation problems while circumventing the potential mesh distortion without remeshing. To eliminate rotational motion damping and loss of angular momentum during the projection, we introduce the mapping for microinertia and angular momentum between particles and grids through the affine particle-in-cell approach. The microrotation and the curvature at each particle are updated through zero-order forward integration of the microgyration and its spatial gradient. We show that the microinertia and the angular momentum are conserved during the projections between particles and grids in our formulation. We verify the formulation and implementation by comparing with the analytical dispersion relation of micropolar waves under the small strain and small microrotation, as well as the analytical soliton solution for solids undergoing large deformation and large microrotation. Additionally we also demonstrate the capacity of the proposed computational framework to handle a wide spectrum of simulations that exhibit size effects in the geometrical nonlinear regime through three representative numerical examples, i.e., a cantilever beam torsion problem, a fragment-impact penetration problem, and a micropolar fluid discharging problem.

42 ENGINEERING↗
NASA NTRS2006

Simulations of Sea-Ice Dynamics Using the Material-Point Method

In recent years, the availability of large volumes of recorded ice motion derived from high-resolution SAR data has provided an amazingly detailed look at the deformation of the ice cover. The deformation is dominated by the appearance of linear kinematic features that have been associated with the presence of leads. These remarkable data put us in a position to begin detailed evaluation of current coupled mechanical and thermodynamic models of sea ice. This presentation will describe the material point method (MPM) for solving these model equations. MPM is a numerical method for continuum mechanics that combines the best aspects of Lagrangian and Eulerian discretizations. The material points provide a Lagrangian description of the ice that models convection naturally. Thus, properties such as ice thickness and compactness are computed in a Lagrangian frame and do not suffer from errors associated with Eulerian advection schemes, such as artificial diffusion, dispersion, or oscillations near discontinuities. This desirable property is illustrated by solving transport of ice in uniform, rotational and convergent velocity fields. Moreover, the ice geometry is represented by unconnected material points rather than a grid. This representation facilitates modeling the large deformations observed in the Arctic, as well as localized deformation along leads, and admits a sharp representation of the ice edge. MPM also easily allows the use of any ice constitutive model. The versatility of MPM is demonstrated by using two constitutive models for simulations of wind-driven ice. The first model is a standard viscous-plastic model with two thickness categories. The MPM solution to the viscous-plastic model agrees with previously published results using finite elements. The second model is a new elastic-decohesive model that explicitly represents leads. The model includes a mechanism to initiate leads, and to predict their orientation and width. The elastic-decohesion model can provide similar overall deformation as the viscous-plastic model; however, explicit regions of opening and shear are predicted. Furthermore, the efficiency of MPM with the elastic-decohesive model is competitive with the current best methods for sea ice dynamics. Simulations will also be presented for an area of the Beaufort Sea, where predictions can be validated against satellite observations of the Arctic.

Sulsky, D.↗
DOE OSTI2020

Deformation accommodating periodic computational domain for a uniform velocity gradient

Many multiscale methods for granular materials use periodic computational domains to consider particle scale interactions and then calculate the stress to drive the continuum scale calculations. For problems involving large material deformations, the computation domain often needs to be reinitialized because of the distortion, causing the loss of the history information of the system. This work introduces an algorithm to accommodate a large deformation of the material while maintaining the computational domain cuboid to avoid the domain reinitialization during the computation. The algorithm uses a rotating frame of reference, in which the velocity gradient can be represented by an upper triangular matrix. The deformation caused by the upper triangular matrix is treated by the image system implied in the periodicity to maintain the computational domain cuboid. The effect from the rotation of the reference frame is considered using the inertial forces. Finally, simulations of simple and pure shear motions are carried out to illustrate the algorithm.

42 ENGINEERING↗
NASA NTRSDate not supplied

TPSAS-NF1676L-32167-DND

Fiber Kinking Theory (FKT): Budiansky proposed that kink bands are a result of an initial (infinitely wide) band of misaligned fibers, nonlinear material behavior under shear loading, and accounting for large rotations of the fibers. CompDam is a NASA Langley CDM code that uses the deformation gradient, particular attention paid to kinematics of cracking.

Andrew Bergan↗
NASA NTRS1998

Mechanistic Models of Friction Stir Welding

Friction stir welding is a welding process developed at The Welding Institute (TWI) in England. The method uses very large strain plastic deformation of the material to join two pieces of metal together. The material is deformed using a tool which is forced between the two pieces which rotates causing a bond. Beyond this, very little is actually known although many people working in the field are willing to speculate on the detailed mechanisms involved. Some measurements made using sacrificial thermocouples at the weld joint indicate that the maximum temperature during the weld process is on the order of 370C - well below the melting temperature of the material. However, at this temperature, the material properties are highly temperature dependent, and the yield stress is approximately an order of magnitude less at this temperature than it is at room temperature. As expected, there are many interpretations of the physical mechanisms occurring during the weld process. Although there is very little published concerned with FSW, some of the anecdotal theories will be described. One describes the primary mechanism as frictional heating at the front of the tool caused by slip between the tool and the material. At elevated temperatures, the weld material becomes soft and deforms around the tool but not essentially altered by the tool rotation, similar to an extrusion. As the material meets again at the rear of the tool, the temperatures and pressures are sufficient to cause the material to bond. All other structures seen are secondary and unimportant. Another theory examined last summer at NASA's Marshall Space Flight Center (MSFC) was that there was no slip between the tool and the material resulting in a rotating mass of plastic weld material traveling at a variety of angular velocities - the greatest at the tool surface diminishing to zero at the outer edge of the plastic mass surrounding the tool. This conceptual model was followed by simplified calculations which showed that the balance of moments through the weld plug was not possible under steady state conditions and realistic temperature profiles. This led to some consideration of a quasi-steady oscillating process. Later when force measurements became available some models were modified and new ones were proposed.

Stewart, Michael B.↗
NASA NTRS1995

Load transfer in the stiffener-to-skin joints of a pressurized fuselage

Structural analyses are developed to determine the linear elastic and the geometrically nonlinear elastic response of an internally pressurized, orthogonally stiffened, composite material cylindrical shell. The configuration is a long circular cylindrical shell stiffened on the inside by a regular arrangement of identical stringers and identical rings. Periodicity permits the analysis of a unit cell model consisting of a portion of the shell wall centered over one stringer-ring joint. The stringer-ring-shell joint is modeled in an idealized manner; the stiffeners are mathematically permitted to pass through one another without contact, but do interact indirectly through their mutual contact with the shell at the joint. Discrete beams models of the stiffeners include a stringer with a symmetrical cross section and a ring with either a symmetrical or an asymmetrical open section. Mathematical formulations presented for the linear response include the effect of transverse shear deformations and the effect of warping of the ring's cross section due to torsion. These effects are important when the ring has an asymmetrical cross section because the loss of symmetry in the problem results in torsion and out-of-plane bending of the ring, and a concomitant rotation of the joint at the stiffener intersection about the circumferential axis. Data from a composite material crown panel typical of a large transport fuselage structure are used for two numerical examples. Although the inclusion of geometric nonlinearity reduces the 'pillowing' of the shell, it is found that bending is localized to a narrow region near the stiffener. Including warping deformation of the ring into the analysis changes the sense of the joint rotation. Transverse shear deformation models result in increased joint flexibility.

Johnson, Eric R.↗
NASA NTRS2018

Modeling Fiber Kinking at the Microscale and Mesoscale

A computational micromechanics (CMM) model is employed to interrogate the assumptions of a recently developed mesoscale continuum damage mechanics (CDM) model for fiber kinking. The CMM model considers an individually discretized three dimensional fiber and surrounding matrix accounting for nonlinearity in the fiber, matrix plasticity, fiber/matrix interface debonding, and geometric nonlinearity. Key parameters of the CMM model were measured through experiments. In particular, a novel experimental technique to characterize the in situ longitudinal compressive strength of carbon fibers through indentation of micropillars is presented. The CDM model is formulated on the basis of Budiansky's fiber kinking theory (FKT) with a constitutive deformation-decomposition approach to alleviate mesh size sensitivity. In contrast to conventional mesoscale CDM models that prescribe a constitutive response directly, the response of the proposed model is an outcome of material nonlinearity and large rotations of the fiber direction following FKT. Comparison of the predictions from the CMM and CDM models shows remarkable correlation in strength, post-peak residual stress, and fiber rotation, with less than 10% difference between the two models in most cases. Additional comparisons are made with several fiber kinking models proposed in the literature to highlight the efficacy of the two models. Finally, the CMM model is exercised in parametric studies to explore opportunities to improve the longitudinal compression strength of a ply through the use of nonconventional microstructures.

Herraez, Miguel↗
NASA NTRS1994

Global, quasi-steady dynamics of the distant solar wind 2: Deformation of the heliospheric current sheet

A three-dimensional (3-D) magnetohydrodynamic (MHD) numerical model is used to trace global deformations of the heliospheric current sheet (HCS) caused by large-scale dynamical interactions associated with corotating solar wind flows. Configurations incorporating the tilted-dipole geometry are investigated out past 30 AU for a variety of dipole tilt angles, alpha. Inclusion of the full, 3-D interplanetary dynamics allows north-south displacements and the east-west warping of the HCS by advective corotational effects to be accurately assessed for the first time. It is found that large-scale spatial correlations between velocity and density imposed at the coronal source (i.e., the geometric arrangement whereby regions of high-speed low-density material lying adjacent to areas of slow, dense flow interact obliquely under the influence of solar rotation) result in a distinctive pattern of deformation of the HCS. For an alpha = 30 deg tilted-dipole example, it is shown that typical zonal variations in radial velocity lead to significant folding of the HCS within about 5 AU of the Sun. By 10 AU, additional sharp bends appear near the latitudinal extremes of the HCS surface, where it is overtaken by shock fronts driven by 3-D corotating interaction regions (CIRs). Moreover, the model suggests that inside about 20 AU, major plasma structures are systematically organized, about the HCS, such that the greatest concentrations of material and magnetic field (the centroids of the 3-D CIR structures) are coincident with the folded crests of the HCS (near heliographic latitudes lambda = +/- alpha). Thus, in these circumstances many of the more interesting dynamical features inconveniently lie well away from the heliographic equator. At larger heliocentric distances, where neighboring CIRs begin to interact strongly, the warping of the HCS abates dramatically and the association between folds in the HCS and major field and density concentrations is weakened and ultimately breaks down. On the basis of the model calculations, quasi-steady deformations and related phenomena should figure prominently in the interplanetary medium, except during those periods when the magnetic dipole and solar spin axes are nearly aligned (alpha is approximately less than 10 deg). Even then, they may be of significance if there are any substantial local warps or kinks in the streamer belt.

Pizzo, V. J.↗
DOE OSTI2022

New large-strain FFT-based formulation and its application to model strain localization in nano-metallic laminates and other strongly anisotropic crystalline materials

This paper presents a new robust large-strain (LS) elasto-viscoplastic (EVP) formulation based on Fast Fourier Transforms (FFTs) for the prediction of the micro-mechanical response and microstructure evolution of polycrystalline and multiphase materials, with emphasis on the effect of strong crystallographic and/or morphologic anisotropy on localization of plastic deformation. In this work, the novel LS-EVPFFT formulation allows treatment of complex initial geometries and large deformations considering three grids of material points: a regular grid in the reference configuration, where FFTs can be performed; an irregular grid in the initial configuration, created by applying a stress-free displacement field to the reference regular grid; and an irregular grid in the current configuration, undergoing large strains and rotations as the material is loaded. Further numerical stability of the new formulation also required the use of a novel expression for the discrete modified Green’s operator, which reduces spurious field oscillations. After presenting and validating the new formulation by comparison with preexisting implementations and analytical solutions, LS-EVPFFT is applied to the prediction of slip and kink bands formation in polycrystalline columnar ice, and kink bands in single crystal zinc wires, showing good agreement with classic experiments. Finally, the model is used to study kink band formation during compression of Cu–Nb nano-metallic laminates (NMLs), in which accurate treatment of the complex geometry associated with the tortuosity of interfaces and large deformations become critical, showing consistency with corresponding micropillar experiments.

36 MATERIALS SCIENCE↗
NASA NTRS1979

Finite-strain large-deflection elastic-viscoplastic finite-element transient response analysis of structures

A method of analysis for thin structures that incorporates finite strain, elastic-plastic, strain hardening, time dependent material behavior implemented with respect to a fixed configuration and is consistently valid for finite strains and finite rotations is developed. The theory is formulated systematically in a body fixed system of convected coordinates with materially embedded vectors that deform in common with continuum. Tensors are considered as linear vector functions and use is made of the dyadic representation. The kinematics of a deformable continuum is treated in detail, carefully defining precisely all quantities necessary for the analysis. The finite strain theory developed gives much better predictions and agreement with experiment than does the traditional small strain theory, and at practically no additional cost. This represents a very significant advance in the capability for the reliable prediction of nonlinear transient structural responses, including the reliable prediction of strains large enough to produce ductile metal rupture.

Rodal, J. J. A.↗
NASA NTRS2015

Representing Matrix Cracks Through Decomposition of the Deformation Gradient Tensor in Continuum Damage Mechanics Methods

A method is presented to represent the large-deformation kinematics of intraply matrix cracks and delaminations in continuum damage mechanics (CDM) constitutive material models. The method involves the additive decomposition of the deformation gradient tensor into 'crack' and 'bulk material' components. The response of the intact bulk material is represented by a reduced deformation gradient tensor, and the opening of an embedded cohesive interface is represented by a normalized cohesive displacement-jump vector. The rotation of the embedded interface is tracked as the material deforms and as the crack opens. The distribution of the total local deformation between the bulk material and the cohesive interface components is determined by minimizing the difference between the cohesive stress and the bulk material stress projected onto the cohesive interface. The improvements to the accuracy of CDM models that incorporate the presented method over existing approaches are demonstrated for a single element subjected to simple shear deformation and for a finite element model of a unidirectional open-hole tension specimen. The material model is implemented as a VUMAT user subroutine for the Abaqus/Explicit finite element software. The presented deformation gradient decomposition method reduces the artificial load transfer across matrix cracks subjected to large shearing deformations, and avoids the spurious secondary failure modes that often occur in analyses based on conventional progressive damage models.

Leone, Frank A., Jr.↗
NASA NTRS1997

ZIP2DL: An Elastic-Plastic, Large-Rotation Finite-Element Stress Analysis and Crack-Growth Simulation Program

ZIP2DL is a two-dimensional, elastic-plastic finte element program for stress analysis and crack growth simulations, developed for the NASA Langley Research Center. It has many of the salient features of the ZIP2D program. For example, ZIP2DL contains five material models (linearly elastic, elastic-perfectly plastic, power-law hardening, linear hardening, and multi-linear hardening models), and it can simulate mixed-mode crack growth for prescribed crack growth paths under plane stress, plane strain and mixed state of stress conditions. Further, as an extension of ZIP2D, it also includes a number of new capabilities. The large-deformation kinematics in ZIP2DL will allow it to handle elastic problems with large strains and large rotations, and elastic-plastic problems with small strains and large rotations. Loading conditions in terms of surface traction, concentrated load, and nodal displacement can be applied with a default linear time dependence or they can be programmed according to a user-defined time dependence through a user subroutine. The restart capability of ZIP2DL will make it possible to stop the execution of the program at any time, analyze the results and/or modify execution options and resume and continue the execution of the program. This report includes three sectons: a theoretical manual section, a user manual section, and an example manual secton. In the theoretical secton, the mathematics behind the various aspects of the program are concisely outlined. In the user manual section, a line-by-line explanation of the input data is given. In the example manual secton, three types of examples are presented to demonstrate the accuracy and illustrate the use of this program.

Deng, Xiaomin↗
DOE OSTI2024

Three-dimensional continuum point cloud method for large deformation and its verification

This study presents a strong form based meshfree collocation method, which is named Continuum Point Cloud Method, to solve nonlinear field equations derived from classical mechanics for deformed bodies in three-dimensional Euclidean space. The method and its implementation are benchmarked against a nonlinear vector field using manufactured solutions. The analysis of mechanical fields firstly focuses on the study of St. Venant Kirchhoff and compressible neo-Hookean materials. Results for various initial boundary value problems are presented, including benchmark cases involving unidirectional tension and simple shear. Subsequently, the study concludes with an analysis of a displacement-controlled simulation of a compressible neo-Hookean material, specifically a bar that is pulled to 50% of its original length and rotated 90°. The pure tension case yields a 1.5% error in displacement between computed and expected values and a combined tension and torsion loading case provides further insight into material behavior under complex loading conditions. The resulting normal axial and transverse stress-strain curves are also presented. Lastly, the consistency and robustness of the proposed nonlinear numerical schemes are successfully demonstrated through various numerical experiments.

Compressible neo-Hookean materials↗
DOE OSTI2025

Meshfree simulation and prediction of recrystallized grain size in friction stir processed 316L stainless steel

Friction stir processing (FSP) is a promising solid-phase microstructural modification technique that can repair and enhance damaged stainless steel surfaces exposed to harsh environments. The quality of the repaired material is closely correlated to the recrystallized grain size in the stir zone (SZ), which is influenced by the thermomechanical conditions dictated by FSP process parameters. Thus, establishing a reliable relationship between these parameters and recrystallized grain size in the SZ is crucial for optimizing repair quality. However, existing experimental approaches often rely on indirect temperatures measured far from the SZ, along with rough strain rate estimations, which are imprecise and time-consuming. Meanwhile, existing mesh-based modeling methods usually face numerical challenges when dealing with the large material deformations inherent in FSP. Here, to address these issues, this study introduces a meshfree process model for FSP based on the smoothed particle hydrodynamics (SPH) method, aimed at predicting process conditions under different parameters. The model is validated using experimental data from 11 combinations of tool traverse and rotation speeds on 316 L stainless steel. Correlations between process parameters, material flow, temperature, strain, strain rate, and recrystallized grain size are revealed through SPH simulations and electron backscatter diffraction (EBSD) imaging. The results show that in situ SZ temperatures range from 1071 to 1322°C, which exceed the tool temperature by over 300°C. Furthermore, SZ temperature, strain rate, and grain size increase monotonically with higher tool temperature and faster traverse speed. A relationship is then established between the model-predicted Zener-Hollomon parameter and the recrystallized grain size based on EBSD data, expressed as ln(d) = -0.364 ln(Z) + 14.673. Finally, this relationship exhibits satisfactory accuracy with errors of less than 26.9% in predicting grain sizes at various SZ locations, which offers valuable insights for optimizing FSP repair processes for 316 L stainless steel.

316L stainless steel↗
DOE OSTI2020

Effect of grain structure and strain rate on dynamic recrystallization and deformation behavior: A phase field-crystal plasticity model

Herein, a phase field (PF)-crystal plasticity (CP) model has been developed to study the effects of initial grain structures, recrystallized grain orientation, and strain rates on dynamic recrystallization and deformation behavior under stress states typical of the Shear Assisted Processing and Extrusion (ShAPE) process. The model takes the geometrically necessary dislocation density as the nucleation criterion for recrystallization and minimizes the deformation energy as the driving force of recrystallized grain growth. Using this model, the effects of grain orientations and strain rates on dynamic recrystallization and deformation behavior were simulated. The results demonstrate that (1) a polycrystalline structure with the texture observed in ShAPE process has the lowest yield stress under the ShAPE stress conditions, (2) the recrystallized grains with the texture observed in ShAPE process largely soften the materials, and (3) higher shear strain rate or rotation speed results in larger magnitude of material softening.

36 MATERIALS SCIENCE↗