Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shear deformation”

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 1,207 records · Page 67

Stress-controlled medium-amplitude oscillatory shear (MAOStress) of PVA–Borax

We report the first-ever complete measurement of MAOStress material functions, which reveal that stress can be more fundamental than strain or strain rate for understanding linearity limits as a function of Deborah number. The material used is a canonical viscoelastic liquid with a single dominant relaxation time: polyvinyl alcohol (PVA) polymer solution cross-linked with tetrahydroborate (Borax) solution. We outline experimental limit lines and their dependence on geometry and test conditions. These MAOStress measurements enable us to observe the frequency dependence of the weakly nonlinear deviation as a function of stress amplitude. The observed features of MAOStress material functions are distinctly simpler than MAOStrain, where the frequency dependence is much more dramatic. The strain-stiffening transient network model was used to derive a model-informed normalization of the nonlinear material functions that accounts for their scaling with linear material properties. Moreover, we compare the frequency dependence of the critical stress, strain, and strain-rate for the linearity limit, which are rigorously computed from the MAOStress and MAOStrain material functions. While critical strain and strain-rate change by orders of magnitude throughout the Deborah number range, critical stress changes by a factor of about 2, showing that stress is a more fundamental measure of nonlinearity strength. This work extends the experimental accessibility of the weakly nonlinear regime to stress-controlled instruments and deformations, which reveal material physics beyond linear viscoelasticity but at conditions that are accessible to theory and detailed simulation.

Mechanics↗

Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions

Many materials of interest to Sandia transition from fluid to solid or have regions of both phases coexisting simultaneously. Currently there are, unfortunately, no material models that can accurately predict this material response. This is relevant to applications that "birth stress" related to geoscience, nuclear safety, manufacturing, energy production and bioscience. Accurately capturing solidification and residual stress enables fully predictive simulations of the evolving front shape or final product. Accurately resolving flow of proppants or blood could reduce environmental impact or lead to better treatments for heart attacks, thrombosis, or aneurism. We will address a science question in this proposal: When does residual stress develop during the critical transition from liquid to solid and how does it affect material deformation? Our hypothesis is that these early phases of stress development are critical to predictive simulation of material performance, net shape, and aging. In this project, we use advanced constitutive models with yield stress to represent both fluid and solid behavior simultaneously. The report provides an abbreviated description of the results from our LDRD "Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions," since we have written four papers that document the work in detail and which we reference. We give highlights of the work and describe the gravitationally driven flow visualization experiment on a model yield stress fluid, Carbopol, at various concentrations and flow rates. We were able to collapse the data on a single master curve by showing it was self-similar. We also describe the Carbopol rheology and the constitutive equations of interest including the Bingham-Carreau-Yasuda model, the Saramito model, and the HB-Saramito model including parameter estimation for the shear and oscillatory rheology. We present several computational models including the 3D moving mesh simulations of both the Saramito models and Bingham-Carreau-Yasuda (BCY) model. We also show results from the BCY model using a 3D level set method and two different ways of handling reduced order Hele-Shaw modeling for generalized Newtonian fluids. We present some first ever two-dimensional results for the modified Jeffries Kamani-Donley-Rogers constitutive equation developed during this project. We include some recent results with a successful Saramito-level set coupling that allows us to tackle problems with complex geometries like mold filling in a thin gap with an obstacle, without the need for remeshing or remapping. We report on some experiments for curing systems where fluorescent particles are used to track material flow. These experiments were carried out in an oven on Sylgard 184 as a model polymerizing system. We conclude the report with a summary of accomplishments and some thoughts on follow-on work.

36 MATERIALS SCIENCE↗

Microstructural evolution and intermetallic formation in Zn-3Mg (wt%) powder mixture processed by high-pressure torsion

Severe plastic deformation (SPD) techniques have been used extensively over the past 40 years for producing strong metals and alloys. High-pressure torsion (HPT) is one of the most promising SPD techniques for achieving high strength through nanoscale grain refinement and phase transformation. In this research, a mixture of pure zinc (Zn) and magnesium (Mg) powders, Zn-3Mg (wt%), was HPT-processed under a pressure of 6 GPa for 1, 5, 10, 20, and 30 turns at room temperature to achieve a high strength biodegradable material. In order to understand the effects of pre-consolidation on the resulting microstructure and hardness, HPT processing was performed on loose powders placed in the die and also on a pre-compacted powder mixture and the resulting HPT disks were characterized by X-ray diffraction, scanning electron microscopy, atom probe tomography, and Vickers microhardness. In both cases, the HPT disk microstructures contained nanoscale grains, and stable and metastable strain-induced intermetallics, but an unusual softening appeared at large shear strains. Grain size, grain morphology, and the formation of different intermetallics were analyzed to explain the unusual hardness distribution, and it was found that an inverse Hall-Petch relationship between hardness and grain size exists. It is suggested that thermally-activated phenomena such as grain boundary sliding contributed to the strain-induced softening of this nano-structured biomaterial due to its low melting point. Finally, the current results are compared with those for HPT-processed cast alloys and hybrids of the same composition.

36 MATERIALS SCIENCE↗

Comparison of the Mechanical Properties of a Conjugated Polymer Deposited Using Spin Coating, Interfacial Spreading, Solution Shearing, and Spray Coating

The mechanical properties of π-conjugated (semiconducting) polymers are a key determinant of the stability and manufacturability of devices envisioned for applications in energy and healthcare. These properties—including modulus, extensibility, toughness, and strength—are influenced by the morphology of the solid film, which depends on the method of processing. To date, the majority of work done on the mechanical properties of semiconducting polymers has been performed on films deposited by spin coating, a process not amenable to the manufacturing of large-area films. Here, we compare the mechanical properties of thin films of regioregular poly(3-heptylthiophene) (P3HpT) produced by three scalable deposition processes—interfacial spreading, solution shearing, and spray coating—and spin coating (as a reference). Our results lead to four principal conclusions. (1) Spray-coated films have poor mechanical robustness due to defects and inhomogeneous thickness. (2) Sheared films show the highest modulus, strength, and toughness, likely resulting from a decrease in free volume. (3) Interfacially spread films show a lower modulus but greater fracture strain than spin-coated films. (4) The trends observed in the tensile behavior of films cast using different deposition processes held true for both P3HpT and poly(3-butylthiophene) (P3BT), an analogue with a higher glass transition temperature. Grazing incidence X-ray diffraction and ultraviolet–visible spectroscopy reveal many notable differences in the solid structures of P3HpT films generated by all four processes. Herein, while these morphological differences provide possible explanations for differences in the electronic properties (hole mobility), we find that the mechanical properties of the film are dominated by the free volume and surface topography. In field-effect transistors, spread films had mobilities more than 1 magnitude greater than any other films, likely due to a relatively high proportion of edge-on texturing and long coherence length in the crystalline domains. Overall, spread films offer the best combination of deformability and charge-transport properties.

36 MATERIALS SCIENCE↗

In situ measurement of three-dimensional intergranular stress localizations and grain yielding under elastoplastic axial-torsional loading

The three-dimensional grain-averaged response of solid bar samples under non-proportional (NP) elastoplastic axial-torsional loading was investigated using in situ high energy diffraction microscopy (HEDM) and companion crystal plasticity finite element (CPFE) modeling. Important stress metrics including applied shear (σ θZ ) and axial (σ ZZ ) stress tensor components, stress and stress deviator tensor invariants (I 1 , J 2 , and J 3 ), von Mises equivalent stress (σ$^{grain}_{VM}$), maximum resolved shear stress (mRSS), stress triaxiality (η), and lode angle parameter ($\barθ$) values were tracked for ~300 grains under two different loading conditions: (1) Torsion-dominated loading (low NP) and (2) Tension-torsion loading (high NP) in equiatomic NiCoCr, a representative multicomponent face-centered cubic (FCC) superalloy. Overall, significant stress localizations existed within both samples as evidenced by the radial dependence of grain-resolved σ θZ , σ$^{grain}_{VM}$, and J 2 ; by comparison, I 1 , J 3 , η, and $\barθ$ metrics did not show discernible trends within the volume. These stress localizations reveal a complex interplay between axial and shear stress components (e.g., stress coupling) resulting in grain yielding near the sample surface largely driven by shear stress, whereas internal grain yielding was largely accommodated by axial stress. Grain-resolved stress localization trends were described well by the CPFE model, although some discrepancies in magnitude occurred, particularly for volumetric stress metrics (I 1 and η) due to initial type II residual stress distributions. The superposition of initial residual stress states onto CPFE grain-resolved data significantly improved model accuracy for η. This suggests that residual stresses more strongly influence the simulation of volumetric rather than deviatoric (yield) stress metrics.

36 MATERIALS SCIENCE↗

Effect of Microscopic Damage Events on Static and Ballistic Impact Strength of Triaxial Braid Composites

In previous work, the ballistic impact resistance of triaxial braided carbon/epoxy composites made with large flat tows (12k and 24k) was examined by impacting 2 X2 X0.125" composite panels with gelatin projectiles. Several high strength, intermediate modulus carbon fibers were used in combination with both untoughened and toughened matrix materials. A wide range of penetration thresholds were measured for the various fiber/matrix combinations. However, there was no clear relationship between the penetration threshold and the properties of the constituents. During some of these experiments high speed cameras were used to view the failure process, and full-field strain measurements were made to determine the strain at the onset of failure. However, these experiments provided only limited insight into the microscopic failure processes responsible for the wide range of impact resistance observed. In order to investigate potential microscopic failure processes in more detail, quasi-static tests were performed in tension, compression, and shear. Full-field strain measurement techniques were used to identify local regions of high strain resulting from microscopic failures. Microscopic failure events near the specimen surface, such as splitting of fiber bundles in surface plies, were easily identified. Subsurface damage, such as fiber fracture or fiber bundle splitting, could be identified by its effect on in-plane surface strains. Subsurface delamination could be detected as an out-of-plane deflection at the surface. Using this data, failure criteria could be established at the fiber tow level for use in analysis. An analytical formulation was developed to allow the microscopic failure criteria to be used in place of macroscopic properties as input to simulations performed using the commercial explicit finite element code, LS-DYNA. The test methods developed to investigate microscopic failure will be presented along with methods for determining local failure criteria that can be used in analysis. Results of simulations performed using LS-DYNA will be presented to illustrate the capabilities and limitations for simulating failure during quasi-static deformation and during ballistic impact of large unit cell size triaxial braid composites.

Littell, Justin D.↗

Variational Asymptotic Homogenization of Finitely Deformed Viscoelastic-Viscoplastic Composites

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.

Finite element analysis↗

Sculpting crystals one Burgers vector at a time: Toward colloidal lattice robot swarms

Plastic deformation of crystalline materials with isotropic particle attractions proceeds by the creation and migration of dislocations under the influence of external forces. If dislocations are produced and migrated under the action of local forces, then material shape change can occur without the application of surface forces. We investigate how particles with variable diameters can be embedded in colloidal monolayers to produce dislocations on demand. Additionally, we find in simulation that when embedded clusters of variable diameter particles are taken through multiple cycles of swelling and shrinking, large cumulative plastic slip is produced by the creation and biased motion of dislocation pairs in the solid for embedded clusters of particular geometries. In this way, dislocations emitted by these clusters (biased “dislocation emitters”) can be used to reshape colloidal matter. Our results are also applicable to larger-scale swarms of robotic particles that organize into dense ordered two-dimensional (2D) arrangements. We conclude with a discussion of how dislocations fulfill for colloids the role sought by “metamodules” in lattice robotics research and show how successive applications of shear as a unit operation can produce shape change through slicing and swirling.

36 MATERIALS SCIENCE↗

The dynamics of oceanic fronts. Part 1: The Gulf Stream

The establishment and maintenance of the mean hydrographic properties of large scale density fronts in the upper ocean is considered. The dynamics is studied by posing an initial value problem starting with a near surface discharge of buoyant water with a prescribed density deficit into an ambient stationary fluid of uniform density. The full time dependent diffusion and Navier-Stokes equations for a constant Coriolis parameter are used in this study. Scaling analysis reveals three independent length scales of the problem, namely a radius of deformation or inertial length scale, Lo, a buoyance length scale, ho, and a diffusive length scale, hv. Two basic dimensionless parameters are then formed from these length scales, the thermal (or more precisely, the densimetric) Rossby number, Ro = Lo/ho and the Ekman number, E = hv/ho. The governing equations are then suitably scaled and the resulting normalized equations are shown to depend on E alone for problems of oceanic interest. Under this scaling, the solutions are similar for all Ro. It is also shown that 1/Ro is a measure of the frontal slope. The governing equations are solved numerically and the scaling analysis is confirmed. The solution indicates that an equilibrium state is established. The front can then be rendered stationary by a barotropic current from a larger scale along-front pressure gradient. In that quasisteady state, and for small values of E, the main thermocline and the inclined isopycnics forming the front have evolved, together with the along-front jet. Conservation of potential vorticity is also obtained in the light water pool. The surface jet exhibits anticyclonic shear in the light water pool and cyclonic shear across the front.

Kao, T. W.↗

An Axial-Torsional, Thermomechanical Fatigue Testing Technique

A technique for conducting strain-controlled, thermomechanical, axial-torsional fatigue tests on thin-walled tubular specimens was developed. Three waveforms of loading, namely, the axial strain waveform, the engineering shear strain waveform, and the temperature waveform were required in these tests. The phasing relationships between the mechanical strain waveforms and the temperature and axial strain waveforms were used to define a set of four axial-torsional, thermomechanical fatigue (AT-TMF) tests. Real-time test control (3 channels) and data acquisition (a minimum of 7 channels) were performed with a software program written in C language and executed on a personal computer. The AT-TMF testing technique was used to investigate the axial-torsional thermomechanical fatigue behavior of a cobalt-base superalloy, Haynes 188. The maximum and minimum temperatures selected for the AT-TMF tests were 760 and 316 C, respectively. Details of the testing system, calibration of the dynamic temperature profile of the thin-walled tubular specimen, thermal strain compensation technique, and test control and data acquisition schemes, are reported. The isothermal, axial, torsional, and in- and out-of-phase axial-torsional fatigue behaviors of Haynes 188 at 316 and 760 C were characterized in previous investigations. The cyclic deformation and fatigue behaviors of Haynes 188 in AT-TMF tests are compared to the previously reported isothermal axial-torsional behavior of this superalloy at the maximum and minimum temperatures.

Kalluri, Sreeramesh↗

Bulk nanocrystalline Al–Mg–Y alloys with amorphous grain boundary complexions display high strength and compressive plasticity

Although nanocrystalline alloys regularly exhibit high strengths, their use in structural applications often face challenges due to sample size limitations, unstable microstructures, and the limited ability to plastically deform. The incorporation of amorphous grain boundary complexions has been proposed to address these issues, by simultaneously stabilizing nanocrystalline grain structures for scale-up processing and improving alloy toughness. In the present study, the mechanical behavior of bulk nanocrystalline Al–Mg–Y is examined with macroscale compression testing, probing a length scale that is relevant to real-world structural applications. Bulk samples were fabricated via a simple powder metallurgy approach, with different hot-pressing temperatures and durations employed for consolidation in order to investigate microstructural and property evolution. All of the specimens contained primary face-centered cubic Al and secondary Al 4 C 3 and Al 3 Y phases, with the Al 3 Y particles exhibiting two populations of small equiaxed and larger elongated particles. Appreciable plasticity was measured along with high ultimate stresses over 800 MPa due to the presence of amorphous grain boundary complexions. Microstructural characterization of fracture surfaces revealed that the area fraction of dimpled regions increased with longer hot-pressing time. Most importantly, the elongated Al 3 Y particles formed regular cellular patterns with increasing hot-pressing time, delaying shear localization and significantly enhancing plasticity. The hierarchy present in the microstructure of the Al–Mg–Y alloy, from amorphous grain boundary complexions to secondary phases, gives rise to excellent bulk mechanical properties, which are attractive for structural applications.

Bulk nanocrystalline alloy↗

Effects of anisotropic shale creep on the stress and permeability evolution of a geological nuclear waste repository

To ensure long-term safety and performance, geological nuclear waste repositories require low-permeability barriers such as bentonite buffers and/or shale host rock. Shale is not only known for its low permeability but also for its trend to undergo time-dependent deformation (i.e., creep), which could heal damage, but the effects of shale creep on the long-term performance of nuclear waste repositories have not been clearly understood. In particular, the anisotropic nature of shale (i.e., bedding) could have a significant effect on its creep behavior, and consequently, on the long-term performance of nuclear waste repositories. In this research, numerical simulations were carried out with the objective of showing the effects of anisotropic shale creep on the stress and permeability evolution of a generic geological nuclear waste repository in shale. The TOUGH-FLAC simulator, a thermo-hydromechanically (THM) coupled numerical code, was used for the simulations. To achieve the objective, comparisons were performed between the results of anisotropic shale creep simulations and those of different simulation cases, namely, no creep (i.e., elastic), isotropic creep, and long-term creep shale cases. Results of the comparisons show that the elastic and isotropic creep shale cases respectively led to the overestimation and underestimation of stress and permeability in the repository, whereas the long-term creep shale case, which accumulated greater creep in later periods than in earlier periods, helped to keep large shear and tensile stresses from developing while maintaining compressive spherical stress, resulting in consistently low permeability levels. These results indicate that performance assessments with elastic and isotropic creep formation models will provide the upper and lower bound estimates of stress and permeability, while more reasonable estimates will be given by an anisotropic creep formation model, and that shale with long-term creep characteristics will be beneficial in many aspects of the safety and performance of nuclear waste repositories.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Synchrotron CT imaging of lattice structures with engineered defects

Understanding mechanical failure, crack propagation, and compressive behavior at the micrometer scale is essential for tailoring material properties for structural performance in cellular materials. Typically, modeling of traditional polymer foam materials is clouded by the lack of control in material morphology and its inherent stochastic structure. Additive manufacturing with sub-micrometer resolution provides a direct path for experimenters to specifically tailor structures needed by modelers to explicitly probe mechanical performance. Using laboratory-based 3D X-ray computed tomographic imaging (CT), the examination of deformation and damage provides a critical path to understand how these soft materials behave. Additionally, synchrotron CT yields realistic information at higher strain rates to directly validate the robustness of our finite element modeling. For this study, nanolithographic printing was employed to generate a series of engineered lattices with increasing levels of defects through the random removal of ligaments. These structures were mechanically tested and imaged with laboratory-based microCT. Additionally, synchrotron experiments were conducted in which the structures were imaged in 3D at 14 Hz during compression at a 0.4 s –1 strain rate. These 3D images show the changes in the structure as the ligaments bend, buckle and fracture in real time. This technique provides a robust framework for developing our methodologies and future exploration of engineered structures.

36 MATERIALS SCIENCE↗

The Effect Stacking Fault Segregation and Phase Transformations have on Creep Strength in Ni-based Superalloys

In this study, two commercially similar polycrystalline Ni-based disk superalloys (LSHR and ME3) were creep tested at 760C and 552MPa to approximately 0.3% plastic strain. LSHR consistently displayed superior creep properties at this stress/temperature regime even though the microstructural characteristics between the two alloys were comparable. High resolution structural and chemical analysis, however, revealed significant differences between the two alloys among active gamma prime shearing modes involving superlattice intrinsic and extrinsic stacking faults. In ME3, Co and Cr segregation and Ni and Al depletion were observed along the intrinsic faults - revealing a gamma prime -to- gamma phase transformation. Conversely in LSHR, an alloy with a higher W content, Co and W segregation was observed along the intrinsic faults. This observation combined with scanning transmission electron microscopy (STEM) simulations confirm a gamma prime-to-D019 phase transformation along the intrinsic faults in LSHR. Using experimental observations and density functional theory calculations, a novel local phase transformation strengthening mechanism is proposed that could be further utilized to improve the high temperature creep capabilities of Ni-base disk alloys.

Superalloy↗

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work is limited to very simple models. In this project, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. The proposed work is based on previously funded NASA work that allowed us to develop a two-dimensional numerical method for boiling flows and to demonstrate the ability of the method to simulate film boiling. While numerical simulations of multi-fluid flows have been advanced in a major way during the last five years, or so, similar capability for flows with phase change are still in their infancy. Although the feasibility of the proposed approach has been demonstrated, it has yet to be extended and applied to fully three-dimensional simulations. Here, a fully three-dimensional, parallel, grid adaptive code will be developed. The numerical method will be used to study nucleate boiling in microgravity, with particular emphasis on two aspects of the problem: 1) Examination of the growth of bubbles at a wall nucleation site and the instabilities of rapidly growing bubbles. Particular emphasis will be put on accurately capturing the thin wall layer left behind as a bubble expands along a wall, on computing instabilities on bubble surfaces as bubbles grow, and on quantifying the effects of both these phenomena on heat transfer; and 2) Examination of the effect of shear flow on bubble growth and heat transfer.

Tryggvason, Gretar↗

Incorporation of Plasticity and Damage Into an Orthotropic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased usage in the aerospace and automotive industries. While there are several composite material models currently available within commercial transient dynamic finite element codes, several features have been identified as being lacking in the currently available material models that could substantially enhance the predictive capability of the impact simulations. A specific desired feature pertains to the incorporation of both plasticity and damage within the material model. Another desired feature relates to using experimentally based tabulated stress-strain input to define the evolution of plasticity and damage as opposed to specifying discrete input properties (such as modulus and strength) and employing analytical functions to track the response of the material. To begin to address these needs, a combined plasticity and damage model suitable for use with both solid and shell elements is being developed for implementation within the commercial code LS-DYNA. The plasticity model is based on extending the Tsai-Wu composite failure model into a strain-hardening based orthotropic plasticity model with a non-associative flow rule. The evolution of the yield surface is determined based on tabulated stress-strain curves in the various normal and shear directions and is tracked using the effective plastic strain. The effective plastic strain is computed by using the non-associative flow rule in combination with appropriate numerical methods. To compute the evolution of damage, a strain equivalent semi-coupled formulation is used, in which a load in one direction results in a stiffness reduction in multiple coordinate directions. A specific laminated composite is examined to demonstrate the process of characterizing and analyzing the response of a composite using the developed model.

Impact↗

Incorporation of Plasticity and Damage Into an Orthotropic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased usage in the aerospace and automotive industries. While there are several composite material models currently available within commercial transient dynamic finite element codes, several features have been identified as being lacking in the currently available material models that could substantially enhance the predictive capability of the impact simulations. A specific desired feature pertains to the incorporation of both plasticity and damage within the material model. Another desired feature relates to using experimentally based tabulated stress-strain input to define the evolution of plasticity and damage as opposed to specifying discrete input properties (such as modulus and strength) and employing analytical functions to track the response of the material. To begin to address these needs, a combined plasticity and damage model suitable for use with both solid and shell elements is being developed for implementation within the commercial code LS-DYNA. The plasticity model is based on extending the Tsai-Wu composite failure model into a strain-hardening based orthotropic plasticity model with a non-associative flow rule. The evolution of the yield surface is determined based on tabulated stress-strain curves in the various normal and shear directions and is tracked using the effective plastic strain. The effective plastic strain is computed by using the non-associative flow rule in combination with appropriate numerical methods. To compute the evolution of damage, a strain equivalent semi-coupled formulation is used, in which a load in one direction results in a stiffness reduction in multiple coordinate directions. A specific laminated composite is examined to demonstrate the process of characterizing and analyzing the response of a composite using the developed model.

Impact↗

Manifestations of Dynamic Strain Aging in Soft-Oriented NiAl Single Crystals

The tensile and compressive properties of six NiAl-base single-crystal alloys have been investigated at temperatures between 77 and 1200 K. The normalized critical resolved shear stresses (CRSS/E) and work-hardening rates (Theta/E) for these alloys generally decreased with increasing temperature. However, anomalous peaks or plateaus for these properties were observed in conventional purity (CPNiAl), Si-doped (NiAl-Si), C-doped low Si (UF-NiAl1), and Mo-doped (NiAl-Mo) alloys at intermediate temperatures (600 to 1000 K). This anomalous behavior was not observed in high-purity, low interstitial material (HP-NiAl). Low or negative strain-rate sensitivities (SRS) also were observed in all six alloys in this intermediate temperature range. Coincident with the occurrence of negative strain-rate sensitivities was the observation of serrated stress-strain curves in the CPNiAl and NiAl-Si alloys. These phenomena have been attributed to dynamic strain aging (DSA). Chemical analysis of the alloys used in this study suggests that the main specie responsible for strain aging in NiAl is C but indicate that residual Si impurities can enhance the strain aging effects. The corresponding dislocation microstructures at low temperatures (300 to 600 K) were composed of well-defined cells. At intermediate temperatures (600 to 900 K), either poorly defined cells or coarse bands of localized slip, reminiscent of the vein structures observed in low-cycle fatigue specimens deformed in the DSA regime, were observed in conventional purity, Si-doped, and in Mo-doped alloys. In contrast, a well-defined cell structure persisted in the low interstitial, high-purity alloy. At elevated temperatures (greater than or equal to 1000 K), more uniformly distributed dislocations and sub-boundaries were observed in all alloys. These observations are consistent with the occurrence of DSA in NiAl single-crystal alloys at intermediate temperatures.

Weaver, M. L.↗