Search NASA⌕ Search

SEARCH · Search NASA

Results for “General Stress Response”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Modeling the Effective Elasticity of Anisotropic Porous Materials

The development and optimization of composite materials designed for thermal protection of NASA’s spacecraft require understanding their physical response to high-enthalpy environments. To predict their macro-scale properties and behavior, high-fidelity 3D simulations are performed at the microscale on realistic representations of these composites. The digital microstructures are generated either synthetically or through X-ray micro-computed tomography reconstructions. One of the main challenges in the prediction of the structural response of heatshield materials is the computation of the effective elasticity of the fibrous composite, as well as the understanding of the deformation and stresses generated at the microscale. These are driven by the fiber layout within the microstructure and the distribution of the infused matrix. In this effort, the micro-mechanical linear elastic behavior of fibrous ablators is modeled using a numerical method based on the Multi-Point Stress Approximation (MPSA) finite volume scheme, a generalization of the more commonly used Multi-Point Flux Approximation (MPFA) that was presented at the 10th Ablation Workshop. To predict the behavior of fibrous and woven architectures, algorithms that compute the local fiber orientation are used. The implementation of the MPSA was verified using analytical solutions, engineering test cases, and compared against legacy Finite Element Analysis (FEA) software. The stress analysis models were then applied to real geometries used by NASA in thermal protection systems such as fibrous preforms and woven materials and the results were compared to experimental data.

Elasticity↗
NASA NTRS1988

Data requirements to model creep in 9Cr-1Mo-V steel

Models for creep behavior are helpful in predicting response of components experiencing stress redistributions due to cyclic loads, and often the analyst would like information that correlates strain rate with history assuming simple hardening rules such as those based on time or strain. On the one hand, much progress has been made in the development of unified constitutive equations that include both hardening and softening through the introduction of state variables whose evolutions are history dependent. Although it is difficult to estimate specific data requirements for general application, there are several simple measurements that can be made in the course of creep testing and results reported in data bases. The issue is whether or not such data could be helpful in developing unified equations, and, if so, how should such data be reported. Data produced on a martensitic 9Cr-1Mo-V-Nb steel were examined with these issues in mind.

Swindeman, R. W.↗
NASA NTRS1996

Element Library for Three-Dimensional Stress Analysis by the Integrated Force Method

The Integrated Force Method, a recently developed method for analyzing structures, is extended in this paper to three-dimensional structural analysis. First, a general formulation is developed to generate the stress interpolation matrix in terms of complete polynomials of the required order. The formulation is based on definitions of the stress tensor components in term of stress functions. The stress functions are written as complete polynomials and substituted into expressions for stress components. Then elimination of the dependent coefficients leaves the stress components expressed as complete polynomials whose coefficients are defined as generalized independent forces. Such derived components of the stress tensor identically satisfy homogenous Navier equations of equilibrium. The resulting element matrices are invariant with respect to coordinate transformation and are free of spurious zero-energy modes. The formulation provides a rational way to calculate the exact number of independent forces necessary to arrive at an approximation of the required order for complete polynomials. The influence of reducing the number of independent forces on the accuracy of the response is also analyzed. The stress fields derived are used to develop a comprehensive finite element library for three-dimensional structural analysis by the Integrated Force Method. Both tetrahedral- and hexahedral-shaped elements capable of modeling arbitrary geometric configurations are developed. A number of examples with known analytical solutions are solved by using the developments presented herein. The results are in good agreement with the analytical solutions. The responses obtained with the Integrated Force Method are also compared with those generated by the standard displacement method. In most cases, the performance of the Integrated Force Method is better overall.

Kaljevic, Igor↗
NASA NTRS2003

Load Diffusion in Composite and Smart Structures

The research carried out here builds on our previous NASA supported research on the general topic of edge effects and load diffusion in composite structures. Further fundamental solid mechanics studies were carried out to provide a basis for assessing the complicated modeling necessary for the multi-functional large scale structures used by NASA. An understanding of the fundamental mechanisms of load diffusion in composite subcomponents is essential in developing primary composite structures. Some specific problems recently considered were those of end effects in smart materials and structures, study of the stress response of pressurized linear piezoelectric cylinders for both static and steady rotating configurations, an analysis of the effect of pre-stressing and pre-polarization on the decay of end effects in piezoelectric solids and investigation of constitutive models for hardening rubber-like materials. Our goal in the study of load diffusion is the development of readily applicable results for the decay lengths in terms of non-dimensional material and geometric parameters. Analytical models of load diffusion behavior are extremely valuable in building an intuitive base for developing refined modeling strategies and assessing results from finite element analyses.

Horgan, Cornelius O.↗
NASA NTRS1999

A Quantitative Study of Oxygen as a Metabolic Regulator

An acute reduction in oxygen (O2) delivery to a tissue is generally associated with a decrease in phosphocreatine, increases in ADP, NADH/NAD, and inorganic phosphate, increased rates of glycolysis and lactate production, and reduced rates of pyruvate and fatty acid oxidation. However, given the complexity of the human bioenergetic system and its components, it is difficult to determine quantitatively how cellular metabolic processes interact to maintain ATP homeostasis during stress (e.g., hypoxia, ischemia, and exercise). Of special interest is the determination of mechanisms relating tissue oxygenation to observed metabolic responses at the tissue, organ, and whole body levels and the quantification of how changes in tissue O2 availability affect the pathways of ATP synthesis and the metabolites that control these pathways. In this study, we extend a previously developed mathematical model of human bioenergetics to provide a physicochemical framework that permits quantitative understanding of O2 as a metabolic regulator. Specifically, the enhancement permits studying the effects of variations in tissue oxygenation and in parameters controlling the rate of cellular respiration on glycolysis, lactate production, and pyruvate oxidation. The whole body is described as a bioenergetic system consisting of metabolically distinct tissue/organ subsystems that exchange materials with the blood. In order to study the dynamic response of each subsystem to stimuli, we solve the ordinary differential equations describing the temporal evolution of metabolite levels, given the initial concentrations. The solver used in the present study is the packaged code LSODE, as implemented in the NASA Lewis kinetics and sensitivity analysis code, LSENS. A major advantage of LSENS is the efficient procedures supporting systematic sensitivity analysis, which provides the basic methods for studying parameter sensitivities (i.e., changes in model behavior due to parameter variation). Sensitivity analysis establishes relationships between model predictions and problem parameters (i.e., initial concentrations, rate coefficients, etc). It helps determine the effects of uncertainties or changes in these input parameters on the predictions, which ultimately are compared with experimental observations in order to validate the model. Sensitivity analysis can identify parameters that must be determined accurately because of their large effect on the model predictions and parameters that need not be known with great precision because they have little or no effect on the solution. This capability may prove to be important in optimizing the design of experiments, thereby reducing the use of animals. This approach can be applied to study the metabolic effects of reduced oxygen delivery to cardiac muscle due to local myocardial ischemia and the effects of acute hypoxia on brain metabolism. Other important applications of sensitivity analysis include identification of quantitatively relevant pathways and biochemical species within an overall mechanism, when examining the effects of a genetic anomaly or pathological state on energetic system components and whole system behavior.

Radhakrishnan, Krishnan↗
NASA NTRS1989

A stable numerical solution method in-plane loading of nonlinear viscoelastic laminated orthotropic materials

In response to the tremendous growth in the development of advanced materials, such as fiber-reinforced plastic (FRP) composite materials, a new numerical method is developed to analyze and predict the time-dependent properties of these materials. Basic concepts in viscoelasticity, laminated composites, and previous viscoelastic numerical methods are presented. A stable numerical method, called the nonlinear differential equation method (NDEM), is developed to calculate the in-plane stresses and strains over any time period for a general laminate constructed from nonlinear viscoelastic orthotropic plies. The method is implemented in an in-plane stress analysis computer program, called VCAP, to demonstrate its usefulness and to verify its accuracy. A number of actual experimental test results performed on Kevlar/epoxy composite laminates are compared to predictions calculated from the numerical method.

Gramoll, K. C.↗
NASA NTRS1997

Cycle-Powered Short Radius (1.9 m) Centrifuge: Effect of Exercise Versus Passive Acceleration on Heart Rate in Humans

In addition to extensive use of lower extremity physical exercise training as a countermeasure for the work capacity component of spaceflight deconditioning, some form of additional head-to-foot (+Gz) gravitational (orthostatic) stress may be required to further attenuate or prevent the signs and symptoms (nausea, vertigo, instability, fatigue) of the general reentry syndrome (GRS) that can reduce astronaut performance during landing. Orthostatic (head-to-foot) stress can be induced by standing, by lower body negative pressure, and by +Gz acceleration. One important question is whether acceleration training alone or with concurrent leg exercise would provide sufficient additive stimulation to attenuate the GRS. Use of a new human-powered centrifuge may be the answer. Thus, the purpose for this study was to compare heart rate (HR), i.e., a stress response during human-powered acceleration, in four men (35-62 yr) and two women (30-31 yr) during exercise acceleration versus passive acceleration (by an off-board operator) at 100% (maximal acceleration = A(max)), and at 25%, 50%, and 75% of A(max). Mean (+/-SE) A(max) was 43.7 +/- 1.3 rpm (+3.9 +/- 0.2Gz). Mean HR at exercise A(max) was 189 +/- 13 b/min (50-70 sec run time), and 142 +/- 22 b/min at passive A(max) (40-70 sec run time). Regression of mean HR on the various +Gz levels indicated explained variance (correlations squared) of r(exp 2) = 0.88 (exercise) and r(exp 2) = 0.96 (passive): exercise HR of 107 +/- 4 (25%) to 189 +/- 13 (100%) b/min were 43-50 b/min higher (p less than 0.05) than comparable passive HR of 64 +/- 2 to 142 +/- 22 b/min. Thus, exercise adds significant physiological stress during +Gz acceleration. Inflight use of this combined exercise and acceleration countermeasure may maintain work capacity as well as normalize acceleration and orthostatic tolerances which could attenuate or perhaps eliminate the GRS.

Greenleaf, J. E.↗
NASA NTRS1972

Study of effects of design details on structural response to acoustic excitation

The normal mode vibration characteristics of one-dimensional and two-dimensional panel arrays were investigated analytically and experimentally. The finite element displacement method was used to formulate the structural models. The structural models include a stiffness and consistent mass matrix for thin-walled open-section beams not previously reported, and a modification of a rectangular plate bending element to include a fundamental interior mode for the element as a generalized coordinate. Provision for adding lumped masses to represent accelerometers is also included. For the one-dimensional panel, the normal mode stress resultants are obtained by integrating the equilibrium equations directly. For the two-dimensional panel arrays, the solution for the stress resultants in the cover sheet was attempted by introducing stress functions dependent upon the displacement field of the plate element. This approach, although incomplete, is a natural extension of the one-dimensional analysis for plate elements. The exerimental program which was used to verify the analysis is described and experimental results are compared with the analysis.

Rudder, F. F., Jr.↗
NASA NTRS1983

Dynamic fracture toughnesses of reaction-bonded silicon nitride

The room-temperature dynamic fracture response of reaction-bonded silicon nitride is investigated using a hybrid experimental-numerical procedure. In this procedure, experimentally determined crack velocities are utilized to drive a dynamic finite-element code or dynamic finite-difference code in its generation mode in order to extract numerically the dynamic stress intensity factor of the fracturing specimen. Results show that the dynamic fracture toughness vs crack velocity relations of the two reaction-bonded silicon nitrides do not follow the general trend in those relations of brittle polymers and steel. A definite slow crack velocity during the initial phase of dynamic crack propagation is observed in reaction-bonded silicon nitride, which results in a nonunique dynamic fracture toughness vs crack velocity relation. In addition, it is found that a propagating crack will continue to propagate under a static stress intensity factor substantially lower than K(IC).

Kobayashi, A. S.↗
NASA NTRS1985

Effects of a high mean stress on the high cycle fatigue life of PWA 1480 and correlation of data by linear elastic fracture mechanics

High-cycle fatigue tests using 5-mm-diameter smooth specimens were performed on the single crystal alloy PWA 1480 (001 axis) at 70F (room temperature) in air and at 100F (538C) in vacuum (10 to the -6 power torr). Tests were conducted at zero mean stress as well as at high tensile mean stress. The results indicate that, although a tensile mean stress, in general, reduces life, the reduction in fatigue strength, for a given mean stress at a life of one million cycles, is much less than what is predicted by the usual linear Goodman plot. Further, the material appears to be significantly more resistant to mean stress effects at 1000F than at 70F. Metallographic examinations of failed specimens indicate that failures in all cases are initiated from micropores of sizes of the order of 30 to 40 microns. Since the macroscopic stress-strain response in all cases was observed to be linear elastic, linear elastic fracture mechanics (LEFM) analyses were carried out to determine the crack growth curves of the material assuming that crack initiation from a micropore (a sub o = 40 microns) occurs very early in life. The results indicate that the calculated crack growth rates at an R (defined as the ratio between minimum stress to maximum stress) value of zero are approximately the same at 70F as at 1000F. However, the calculated crack growth rates at other R ratios, both positive and negative, tend to be higher at 70F than at 1000F. Calculated threshold effects at large R values tend to be independent of temperature in the temperature regime studied. They are relatively constant with increasing R ratio up to a value of about 0.6, beyond which the calculated threshold stress intensity factor range decreases rapidly with increasing R ratios.

Majumdar, S.↗
NASA NTRS2015

Development of Long-Lifetime Pulsed Gas Valves for Pulsed Electric Thrusters

It is advantageous for gas-fed pulsed electric thrusters to employ pulsed valves so propellant is only flowing to the device during operation. The propellant utilization of the thruster will be maximized when all the gas injected into the thruster is acted upon by the fields produced by the electrical pulse. Gas that is injected too early will diffuse away from the thruster before the electrical pulse can act to accelerate the propellant. Gas that is injected too late will miss being accelerated by the already-completed electrical pulse. As a consequence, the valve must open quickly and close equally quickly, only remaining open for a short duration. In addition, the valve must have only a small amount of volume between the sealing body and the thruster so the front and back ends of the pulse are as coincident as possible with the valve cycling, with very little latent propellant remaining in the feed lines after the valve is closed. For a real mission of interest, a pulsed thruster can be expected to pulse at least 10(exp 10) - 10(exp 11) times, setting the range for the number of times a valve must open and close. The valves described in this paper have been fabricated and tested for operation in an inductive pulsed plasma thruster (IPPT) for in-space propulsion. In general, an IPPT is an electrodeless space propulsion device where a capacitor is charged to an initial voltage and then discharged, producing a high-current pulse through a coil. The field produced by this pulse ionizes propellant, inductively driving current in a plasma located near the face of the coil. Once the plasma is formed, it can be accelerated and expelled at a high exhaust velocity by the electromagnetic Lorentz body force arising from the interaction of the induced plasma current and the magnetic field produced by the current in the coil. The valve characteristics needed for the IPPT application require a fast-acting valve capable of a minimum of 10(exp 10) valve actuation cycles. Since even 10(exp 9) cycles is well above anything demonstrated, this lower value was selected as the design point for the present work. The valve seal must remain leak-tight throughout operation, and the body must maintain a low internal leakage at relatively high operating temperatures. The full set of design requirements used for this program are summarized in Table 1. In this work, we describe two pulsed gas valves that have been fabricated to have long lifetime and demonstrate the characteristics listed above. The first is a miniaturized, conventional electromagnet-based valve while the second is a piezoelectric-based valve design. The conventional valve, shown in Fig. 1, is opened by use of a solenoid electromagnetic actuator. When current is applied to the solenoid coil, magnetic forces pull the plunger away from the valve seat, allowing fluid to flow through the valve. Removal of electrical current permits the spring and fluid pressure to seat the plunger, halting the flow of fluid. The valve body is fabricated from 304L corrosion resistant steel (CRES) and while the parts that form the magnetic circuit are fabricated from 430 CRES. This material does not have optimum magnetic properties, but its corrosion resistance permits incorporation into a design without requiring an additional plating process. A viton O-ring compound (Parker V0884-75), selected for its mechanical strength at elevated temperatures, was used for the valve seat seal. The design was based solely on the use of analytical sizing calculations, as opposed to a more rigorous finite element analysis. While this valve is small and relatively lightweight, it does not represent a design that is optimized for mass and/or a given volume envelope. The piezoelectric valve is a "puller" valve design. Applying a voltage to the piezo crystal causes it to elongate and pull a pintle off the seat, opening the valve. The valve seal consists of the pintle with an external, spherically-formed tip fabricated from Inconel 625 and the seat fabricated from Inconel 718 and possessing a spherical surface that permits the pintle to align with the valve body centerline after each cycle, maintaining the circumferential metal-on-metal seal. The contact stress on the valve pintle is expected in the worst case to be 96.5 MPa (14 ksi), which is well below the work hardening value of 413.7 MPa (60 ksi) for Inconel 615. Results presented in this work will compare the performance of each valve, both to each other and to the valve design requirements given in Table 1. In general, both valves meet all the performance envelope requirements, but the piezoelectric valve consumes considerably less power and has a faster response time over the conventional solenoid-driven valves.

Burkhardt, Wendel M.↗
NASA NTRSDate not supplied

Progressive Failure Analysis of 3D Woven Composites via Multiscale Recursive Micromechanics

Multiscale progressive failure simulations have been performed for a 3D woven composite considering six length scales, spanning from the woven composite repeating unit cell to the matrix material containing voids. The Multiscale Recursive Micromechanics approach, which enables micromechanics models to call other micromechanics models (or themselves recursively) to consider finer and finer length scales, has been employed. The multiscale model uses both the generalized method of cells and Mori-Tanaka micromechanics theories and considers failure and damage at the lowest length scale. A simple subvolume elimination damage model, as well as the crack band model, are employed for the matrix and their predicted composite responses compared with and without a preexisting binder tow disbond. Results are also compared with experimental data for an AS4 carbon fiber/RTM6 epoxy matrix 3D orthogonal woven PMC, with good agreement in terms of global stiffness and global failure stress.

3D woven↗
NASA NTRS1985

Probabilistic Structural Analysis Theory Development

The objective of the Probabilistic Structural Analysis Methods (PSAM) project is to develop analysis techniques and computer programs for predicting the probabilistic response of critical structural components for current and future space propulsion systems. This technology will play a central role in establishing system performance and durability. The first year's technical activity is concentrating on probabilistic finite element formulation strategy and code development. Work is also in progress to survey critical materials and space shuttle mian engine components. The probabilistic finite element computer program NESSUS (Numerical Evaluation of Stochastic Structures Under Stress) is being developed. The final probabilistic code will have, in the general case, the capability of performing nonlinear dynamic of stochastic structures. It is the goal of the approximate methods effort to increase problem solving efficiency relative to finite element methods by using energy methods to generate trial solutions which satisfy the structural boundary conditions. These approximate methods will be less computer intensive relative to the finite element approach.

Burnside, O. H.↗
NASA NTRS1994

Evolution of protoplanetary disks with dynamo magnetic fields

The notion that planetary systems are formed within dusty disks is certainly not a new one; the modern planet formation paradigm is based on suggestions made by Laplace more than 200 years ago. More recently, the foundations of accretion disk theory where initially developed with this problem in mind, and in the last decade astronomical observations have indicated that many young stars have disks around them. Such observations support the generally accepted model of a viscous Keplerian accretion disk for the early stages of planetary system formation. However, one of the major uncertainties remaining in understanding the dynamical evolution of protoplanetary disks is the mechanism responsible for the transport of angular momentum and subsequent mass accretion through the disk. This is a fundamental piece of the planetary system genesis problem since such mechanisms will determine the environment in which planets are formed. Among the mechanisms suggested for this effect is the Maxwell stress associated with a magnetic field treading the disk. Due to the low internal temperatures through most of the disk, even the question of the existence of a magnetic field must be seriously studied before including magnetic effects in the disk dynamics. On the other hand, from meteoritic evidence it is believed that magnetic fields of significant magnitude existed in the earliest, PP-disk-like, stage of our own solar system's evolution. Hence, the hypothesis that PP disks are magnetized is not made solely on the basis of theory. Previous studies have addressed the problem of the existence of a magnetic field in a steady-state disk and have found that the low conductivity results in a fast diffusion of the magnetic field on timescales much shorter than the evolutionary timescale. Hence the only way for a magnetic field to exist in PP disks for a considerable portion of their lifetimes is for it to be continuously regenerated. In the present work, we present results on the self-consistent evolution of a turbulent PP disk including the effects of a dynamo-generated magnetic field.

Reyes-Ruiz, M.↗
NASA NTRS2011

The lagRST Model: A Turbulence Model for Non-Equilibrium Flows

This study presents a new class of turbulence model designed for wall bounded, high Reynolds number flows with separation. The model addresses deficiencies seen in the modeling of nonequilibrium turbulent flows. These flows generally have variable adverse pressure gradients which cause the turbulent quantities to react at a finite rate to changes in the mean flow quantities. This "lag" in the response of the turbulent quantities can t be modeled by most standard turbulence models, which are designed to model equilibrium turbulent boundary layers. The model presented uses a standard 2-equation model as the baseline for turbulent equilibrium calculations, but adds transport equations to account directly for non-equilibrium effects in the Reynolds Stress Tensor (RST) that are seen in large pressure gradients involving shock waves and separation. Comparisons are made to several standard turbulence modeling validation cases, including an incompressible boundary layer (both neutral and adverse pressure gradients), an incompressible mixing layer and a transonic bump flow. In addition, a hypersonic Shock Wave Turbulent Boundary Layer Interaction with separation is assessed along with a transonic capsule flow. Results show a substantial improvement over the baseline models for transonic separated flows. The results are mixed for the SWTBLI flows assessed. Separation predictions are not as good as the baseline models, but the over prediction of the peak heat flux downstream of the reattachment shock that plagues many models is reduced.

Lillard, Randolph P.↗
NASA NTRSDate not supplied

Computing the effective elasticity of anisotropic porous media from X-ray computed micro-tomography images

Development and optimization of composite materials designed for thermal protection of NASA’s spacecraft requires the understanding of their physical response to high-enthalpy environments. To predict their macro-scale properties and behavior, high-fidelity 3D simulations are performed at the micro-scale on realistic representations of these composites. The digital micro-structures are generated either synthetically or through X-ray micro-tomography reconstructions. One of the main challenges in the prediction of heatshield material structural response is the computation of the effective elasticity of the fibrous composite, as well as the understanding of the deformation and stresses generated at the micro-scale. These are driven by the fiber layout within the micro-structure and the distribution of the infused matrix. In this effort, the micro-mechanical linear elastic behavior of fibrous ablators is modeled through the use of a numerical method based on the Multi-Point Stress Approximation (MPSA) finite volume scheme. The MPSA, a generalization of the more commonly used Multi-Point Flux Approximation (MPFA), was discussed in a previous presentation at the 15th USCCN and in reference. To predict the behavior of fibrous and woven architectures, algorithms that compute the local fiber orientation are used. The implementation of the MPSA was verified using analytical solutions, engineering test cases and compared against legacy Finite Element Analysis (FEA) software. The stress analysis models were then applied to real geometries used by NASA in thermal protection systems such as fibrous preforms and woven materials and the results were compared to experimental data.

Federico Semeraro↗
NASA NTRS1971

Design study for a high reliability five-year spacecraft tape transport

Following the establishment of the overall transport concept, a study of all of the life limiting constraints associated with the transport were analyzed using modeling techniques. These design techniques included: (1) a response analysis from which the performance of the transport could be determined under operating conditions for a variety of conceptual variations both in a new and aged condition; (2) an analysis of a double cone guidance technique which yielded an optimum design for maximum guidance with minimum tape degradation; (3) an analysis of the tape pack design to eliminate spoking caused by negative tangential stress within the pack; (4) an evaluation of the stress levels experienced by the magnetic tape throughout the system; (5) a general review of the bearing and lubrication technology as applied to satellite recorders and hence the recommendation for using standard load carrying antifriction ball bearings; and (6) a kinetic analysis to determine the change in kinetic properties of the transport during operation.

Benn, G. S. L.↗
NASA NTRS2016

Plasticity Tool for Predicting Shear Nonlinearity of Unidirectional Laminates Under Multiaxial Loading

This study implements a plasticity tool to predict the nonlinear shear behavior of unidirectional composite laminates under multiaxial loadings, with an intent to further develop the tool for use in composite progressive damage analysis. The steps for developing the plasticity tool include establishing a general quadratic yield function, deriving the incremental elasto-plastic stress-strain relations using the yield function with associated flow rule, and integrating the elasto-plastic stress-strain relations with a modified Euler method and a substepping scheme. Micromechanics analyses are performed to obtain normal and shear stress-strain curves that are used in determining the plasticity parameters of the yield function. By analyzing a micromechanics model, a virtual testing approach is used to replace costly experimental tests for obtaining stress-strain responses of composites under various loadings. The predicted elastic moduli and Poisson's ratios are in good agreement with experimental data. The substepping scheme for integrating the elasto-plastic stress-strain relations is suitable for working with displacement-based finite element codes. An illustration problem is solved to show that the plasticity tool can predict the nonlinear shear behavior for a unidirectional laminate subjected to multiaxial loadings.

Wang, John T.↗