Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Functional coupling”

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 919 records · Page 51

APIS: Honeybee Foraging Task Assignment for Use in Uncertain and Unreliable Environments

Multiagent Cyber-Physical-Human (CPH) systems in realistic environments operate under uncertain conditions. Communication among agents, aimed at reducing the uncertainty, is itself subject to uncertainty. We propose to manage uncertainties in autonomous, long-duration operations of multiagent systems via a modified Honeybee Foraging (HBF) behavioral scheme. The resulting system, Autonomous Persistent Intelligent Swarm (APIS),incorporates two new behaviors to ameliorate informational uncertainty. When “scouting”, agents are tasked based on informational quality and reliability rather than solely on priorities. When “dancing”, agents are tasked to rendezvous with other dancing agents to exchange information at close range, where successful communication is guaranteed. When coupled with uncertainty-aware modeling across agents, these behaviors improve situational awareness and resilience of the system, enabling it to function under more uncertain conditions arising during long-duration missions.

Autonomous systems↗

APIS: Honeybee Foraging Task Assignment for Use in Uncertain and Unreliable Environments

Multiagent Cyber-Physical-Human (CPH) systems in realistic environments operate under uncertain conditions. Communication among agents, aimed at reducing the uncertainty, is itself subject to uncertainty. We propose to manage uncertainties in autonomous, long-duration operations of multiagent systems via a modified Honeybee Foraging (HBF) behavioral scheme. The resulting system, Autonomous Persistent Intelligent Swarm (APIS),incorporates two new behaviors to ameliorate informational uncertainty. When “scouting”, agents are tasked based on informational quality and reliability rather than solely on priorities. When “dancing”, agents are tasked to rendezvous with other dancing agents to exchange information at close range, where successful communication is guaranteed. When coupled with uncertainty-aware modeling across agents, these behaviors improve situational awareness and resilience of the system, enabling it to function under more uncertain conditions arising during long-duration missions.

Autonomous systems↗

Assessing Potential Geophysical and Environmental Impacts from Frequent Rocket Launch Missions at Kennedy Space Center

Kennedy Space Center (KSC) in Florida has been utilized for space missions over many years with a gradually increasing number of rocket launches. There have been multiple studies where air-coupled acoustic waves and infrasound originating from launched rockets were used for operational purposes, such as locating booster trajectories as a function of changing atmospheric conditions. However, no study has utilized the acoustoelastic waves as a signal source for subsurface seismic investigations. We conducted a dispersion analysis using the seismic energy recorded from the Artemis I rocket launch at KSC in November 2022, and from these results we generated depth-sensitivity kernels at different wave frequencies. The kernels were compared with sedimentary core data to verify boundaries of carbonate layers above the Floridan Aquifer System. Accumulative information of bedrock-sediment boundary across the sedimentary platform could especially provide geo-structural evidence that manifests the configuration of coastal features. Continuous dispersion analysis of seismic recordings from consecutive rocket launches also has potential to identify non-stationary environmental effects, such as reorientation of sedimentary structures and fluctuation of the saltwater/groundwater lens from gravitational tides, which may affect erosional susceptibility of coastal features.

Han Byul Woo↗

Assessing Potential Geophysical and Environmental Impacts from Frequent Rocket Launch Missions at Kennedy Space Center

Kennedy Space Center (KSC) in Florida has been utilized for space missions over many years with a gradually increasing number of rocket launches. There have been multiple studies where air-coupled acoustic waves and infrasound originating from launched rockets were used for operational purposes, such as locating booster trajectories as a function of changing atmospheric conditions. However, no study has utilized the acoustoelastic waves as a signal source for subsurface seismic investigations. We conducted a dispersion analysis using the seismic energy recorded from the Artemis I rocket launch at KSC in November 2022, and from these results we generated depth-sensitivity kernels at different wave frequencies. The kernels were compared with sedimentary core data to verify boundaries of carbonate layers above the Floridan Aquifer System. Accumulative information of bedrock-sediment boundary across the sedimentary platform could especially provide geo-structural evidence that manifests the configuration of coastal features. Continuous dispersion analysis of seismic recordings from consecutive rocket launches also has potential to identify non-stationary environmental effects, such as reorientation of sedimentary structures and fluctuation of the saltwater/groundwater lens from gravitational tides, which may affect erosional susceptibility of coastal features.

Rocket Launch Seismicity↗

An historical survey of computational methods in optimal control.

Review of some of the salient theoretical developments in the specific area of optimal control algorithms. The first algorithms for optimal control were aimed at unconstrained problems and were derived by using first- and second-variation methods of the calculus of variations. These methods have subsequently been recognized as gradient, Newton-Raphson, or Gauss-Newton methods in function space. A much more recent addition to the arsenal of unconstrained optimal control algorithms are several variations of conjugate-gradient methods. At first, constrained optimal control problems could only be solved by exterior penalty function methods. Later algorithms specifically designed for constrained problems have appeared. Among these are methods for solving the unconstrained linear quadratic regulator problem, as well as certain constrained minimum-time and minimum-energy problems. Differential-dynamic programming was developed from dynamic programming considerations. The conditional-gradient method, the gradient-projection method, and a couple of feasible directions methods were obtained as extensions or adaptations of related algorithms for finite-dimensional problems. Finally, the so-called epsilon-methods combine the Ritz method with penalty function techniques.

Polak, E.↗

Coupled Marangoni-Benard/Rayleigh-Benard Instability with Temperature Dependent Viscosity

The onset of convection induced by coupled surface tension gradient and buoyancy forces is investigated with temperature dependent viscosity. Both surface tension and viscosity are assumed to vary linearly with temperature. The limiting case, Ma=0, is the buoyancy driven convection problem typically referred to as the Rayleigh-Benard stability problem. The other limiting case, Ra=0, is the surface tension (gradient) driven flow problem referred to as the Marangoni-Benard problem. The equations and boundary conditions obtained from the linear analysis are solved numerically as a generalized eigenvalue problem. Neutral stability curves for different viscosity slopes have been generated for the Marangoni-Benard and Rayleigh-Benard problems. It is shown that the curves can be collapsed to a single curve by appropriately scaling the results for each of the limiting cases. The critical Marangoni number is determined as a function of the slope of the viscosity temperature variation, epsilon, for different values of the Rayleigh number. When the viscosity decreases linearly with temperature, the coupled buoyancy-surface tension problem, including the limiting cases of Ra=0 and Ma=0, is found to be more stable than the constant viscosity case.

Skarda, J. Raymond Lee↗

Incorporating Plant Phenology Dynamics in a Biophysical Canopy Model

The Multi-Layer Canopy Model (MLCan) is a vegetation model created to capture plant responses to environmental change. Themodel vertically resolves carbon uptake, water vapor and energy exchange at each canopy level by coupling photosynthesis, stomatal conductance and leaf energy balance. The model is forced by incoming shortwave and longwave radiation, as well as near-surface meteorological conditions. The original formulation of MLCan utilized canopy structural traits derived from observations. This project aims to incorporate a plant phenology scheme within MLCan allowing these structural traits to vary dynamically. In the plant phenology scheme implemented here, plant growth is dependent on environmental conditions such as air temperature and soil moisture. The scheme includes functionality that models plant germination, growth, and senescence. These growth stages dictate the variation in six different vegetative carbon pools: storage, leaves, stem, coarse roots, fine roots, and reproductive. The magnitudes of these carbon pools determine land surface parameters such as leaf area index, canopy height, rooting depth and root water uptake capacity. Coupling this phenology scheme with MLCan allows for a more flexible representation of the structure and function of vegetation as it responds to changing environmental conditions.

environmental changes↗

A spontaneously arising mutation in connexin32 with repeated passage of FRTL-5 cells coincides with increased growth rate and reduced thyroxine release

In this study we examine changes in the cellular properties of FRTL-5 cells as a function of passage number, with particular emphasis on gap junction expression, karyotype, morphology, growth rate and thyroxine (T(4)) release. Early passage FRTL-5 follicular cells transfer dye through gap junctions from injected cell(s) to third-order neighboring cells and beyond within their respective follicles and have immuno-detectable connexin32 (Cx32) type gap junctional plaques in their lateral contacting plasma membranes. By contrast, FRTL-5 cells established as monolayers, or as follicles from cultures passed more than 15 times, did not transfer microinjected Lucifer Yellow dye to contiguous neighboring cells and did not express any immuno-detectable rat thyroid specific connexins (Cx43, Cx32 or Cx26). Western blots confirmed that total, membrane and cytosolic Cx32 protein was present only in early pass follicular cultures. To better understand the passage-dependent loss of Cx32 expression, RT-PCR primers were made to the most unique sequences of the rat Cx32 molecule, the cytoplasmic and carboxyl-terminal regions. These primers were used to screen FRTL-5 RNA from cultures of various passage numbers. The results revealed that later passage cultures had a single base deletion in the middle of the Cx32 cytoplasmic loop region at nucleotide position 378. This base deletion was in the middle position of the codon for amino acid 116, which is normally a CAC (histidine) but read with the frame shift was a CCC (proline). The four amino acids that followed this deletion were also altered with the fourth one becoming UAA, the ochre translation stop codon. This premature stopping of translation resulted in a truncation of 60% of the protein, which included the remaining cytoplasmic loop, third and fourth transmembrane regions and the carboxyl-terminus. The later passage cultures did not produce a carboxyl-terminal RT-PCR product, indicating that the mRNA was also truncated. These regions of the Cx32 molecule contain the sequences and epitopes to which probes and antibodies are directed, and as such alterations of these regions with repeated passage explains reports by others that FRTL-5 cells do not express Cx32, and implies that cultures used for these assessments were passed more than 15 times. To determine if genetic or epigenetic abnormalities existed in FRTL-5 cells we performed chromosome spreads from various passage cultures. FRTL-5 cells have been reported to be diploid and more recently non-diploid; however, we found them to be fully tetraploid. This tetraploidy appears to be unstable in that later passes are tetraploid plus two or three extra chromosomes. There were no obvious translocations, breaks or large-scale interstitial deletions of any chromosomes in the FRTL-5 cultures tested. As FRTL-5 cells were repeatedly passed their morphology changed. Monolayer areas spread from beneath the follicles, and the follicles became flattened in appearance. These physical changes were coincident with dramatically increased growth rates. Early cultures (passed 3-12 times) divided on average every 49+/-1 h, whereas later passes (passes 20-25) divided every 28+/-3 h. To correlate these changes with a measure of thyroid function we assayed T(4) output. Early passage follicular cultures incubated for 6 h with sodium iodide, released on average 5.27+/- 0.33 ng/ml of T(4)/100 follicles. Later passes, or early passes treated with heptanol to down-regulate Cx32, released an average of 3.84+/-0.50 ng/ml of T(4)/100 follicles. There was a 27% difference in T(4) release between early follicular cultures, that were coupled by Cx32, and late or down-regulated early follicular cultures, that were uncoupled (P<0.0001). Collectively, the physical changes documented in this study were coincident with the loss of functional Cx32. This implies a relationship between the loss of intercellular communication and changes in morphogenic appearance, growth rate and reduced thyroid function and supports the previously postulated, tumor-suppressor role for Cx32. FRTL-5 cultures from low passage numbers are an excellent model of primary thyroid cells. However, many reports in the literature ascribe features to FRTL-5 cells that are mutually inconsistent. These differences may be resolved in the future by addressing the passage number and the conditional differences of the cultures being studied.

NASA Discipline Radiation Health↗

Spatial orientation of the vestibular system

1. A simplified three-dimensional state space model of visual vestibular interaction was formulated. Matrix and dynamical system operators representing coupling from the semicircular canals and the visual system to the velocity storage integrator were incorporated into the model. 2. It was postulated that the system matrix for a tilted position was a composition of two linear transformations of the system matrix for the upright position. One transformation modifies the eigenvalues of the system matrix while another rotates the pitch and roll eigenvectors with the head, while maintaining the yaw axis eigenvector approximately spatially invariant. Using this representation, the response characteristics of the pitch, roll, and yaw eye velocity were obtained in terms of the eigenvalues and associated eigenvectors. 3. Using OKAN data obtained from monkeys and comparing to the model predictions, the eigenvalues and eigenvectors of the system matrix were identified as a function of tilt to the side or of tilt to the prone positions, using a modification of the Marquardt algorithm. The yaw eigenvector for right-side-down tilt and for downward pitch cross-coupling was approximately 30 degrees from the spatial vertical. For the prone position, the eigenvector was computed to be approximately 20 degrees relative to the spatial vertical. For both side-down and prone positions, oblique OKN induced along eigenvector directions generated OKAN which decayed to zero along a straight line with approximately a single time constant. This was verified by a spectral analysis of the residual sequence about the straight line fit to the decaying data. The residual sequence was associated with a narrow autocorrelation function and a wide power spectrum. 4. Parameters found using the Marquardt algorithm were incorporated into the model. Diagonal matrices in a head coordinate frame were introduced to represent the direct pathway and the coupling of the visual system to the integrator. Model simulations predicted the behavior of yaw and pitch OKN and OKAN when the animal was upright, as well as the cross-coupling in the tilted position. The trajectories in velocity space were also accurately simulated. 5. There were similarities between the monkey eigenvectors and human perception of the spatial vertical. For side-down tilts and downward eye velocity cross-coupling, there was only an Aubert (A) effect. For upward eye velocity cross-coupling there were both Muller (E) and Aubert (A) effects. The mean of the eigenvectors for upward and downward eye velocities overlay human 1 x g perceptual data.(ABSTRACT TRUNCATED AT 400 WORDS).

NASA Discipline Neuroscience↗

DGEN Aeropropulsion Research Turbofan (DART) Core/Combustor-Noise Infinite-Tube-Probe (ITP) Transfer Function

Direct measurements of turbofan engine core unsteady pressure is complicated by the extreme thermal environment within the engine, preventing currently available transducers from operating. It is necessary to use a remote measurement configuration such as an infinite-tube-probe (ITP) configuration in which the transducer is teed into a line pneumatically coupled with the location of interest and the other end of the ITP contains an "infinite" waveguide designed to prevent any reflections. This configuration has become a standard method for experiments concerning core noise. In order to relate the spectral pressure measurements obtained remotely with the ITP, a transfer (frequency response) function needs to be calculated based on a known pressure field. This is accomplished using a normal impedance tube at Glenn Research Center, where a flush-mounted microphone is used as a reference to the ITP signal in the calculation of the ITP's transfer function.

Boyle, Devin K.↗

Analytical study of the liquid phase transient behavior of a high temperature heat pipe

The transient operation of the liquid phase of a high temperature heat pipe is studied. The study was conducted in support of advanced heat pipe applications that require reliable transport of high temperature drops and significant distances under a broad spectrum of operating conditions. The heat pipe configuration studied consists of a sealed cylindrical enclosure containing a capillary wick structure and sodium working fluid. The wick is an annular flow channel configuration formed between the enclosure interior wall and a concentric cylindrical tube of fine pore screen. The study approach is analytical through the solution of the governing equations. The energy equation is solved over the pipe wall and liquid region using the finite difference Peaceman-Rachford alternating direction implicit numerical method. The continuity and momentum equations are solved over the liquid region by the integral method. The energy equation and liquid dynamics equation are tightly coupled due to the phase change process at the liquid-vapor interface. A kinetic theory model is used to define the phase change process in terms of the temperature jump between the liquid-vapor surface and the bulk vapor. Extensive auxiliary relations, including sodium properties as functions of temperature, are used to close the analytical system. The solution procedure is implemented in a FORTRAN algorithm with some optimization features to take advantage of the IBM System/370 Model 3090 vectorization facility. The code was intended for coupling to a vapor phase algorithm so that the entire heat pipe problem could be solved. As a test of code capabilities, the vapor phase was approximated in a simple manner.

Roche, Gregory Lawrence↗

Coupled elastic response of open section laminated composite beams subject to generalized beam loading

A generalized beam theory for arbitrary open section laminated composite beams subjected to generalized beam loading is developed using a strength of materials approach. At present, the theory is limited to statistically determinate beams. Solution of the general system of equations yields strain and curvature distributions as functions of the axial coordinate x and the contour coordinate s. Average stresses on a given cross section or ply stresses in a local coordinate system are determined using these distributions. Data obtained for thin-walled open section composite beams reveal that elastic coupling occurs among curvatures, shearing strains, and axial strain, depending on the laminate stiffness parameters. Laminate ply layup does not affect the location of the shear center in the thin-wall approximation, but may have a significant impact for thicker walled open sections. The average axial stress distribution is found to be distinctly different from the isotropic distribution, but the average shear stress distribution is the same as that in the isotropic case.

Zvarick, Albert G.↗

Point Spread Function of ASTRO-H Soft X-Ray Telescope (SXT)

ASTRO-H (Hitomi) satellite equips two Soft X-ray Telescopes (SXTs), one of which (SXT-S) is coupled to Soft-X-ray Spectrometer (SXS) while the other (SXT-I) is coupled to Soft X-ray Imager (SXI). Although SXTs are lightweight of approximately 42 kgmodule1 and have large on-axis effective area (EA) of approximately 450 cm(exp 2) at 4.5 keV module(sub 1) by themselves, their angular resolutions are moderate approximately 1.2 arcmin in half power diameter. The amount of contamination into the SXS FOV (3.05 times 3.05 arcmin(exp 2) from nearby sources was measured in the ground-based calibration at the beamline in Institute of Space and Astronautical Science. The contamination at 4.5 keV were measured with sources distant from the SXS center by one width of the FOV in perpendicular and diagonal directions, that is, 3 and 4.5 arcmin-off, respectively. The average EA of the contamination in the four directions with the 3 and 4.5 arcmin-off were measured to be 2 and 0.6% of the on-axis EA of 412 cm (exp) for the SXS FOV, respectively. The contamination from a source distant by two FOV widths in a diagonal direction, that is, 8.6 arcmin-off was measured to be 0.1% of the on-axis at 4.5 keV. The contamination amounts were also measured at 1.5 keV and 8.0 keV which indicated that the ratio of the contamination EA to that of on-axis hardly depended on the source energy. The off-axis SXT-I images from 4.5 to 27 arcmin were acquired at intervals of -4.5 arcmin for the SXI FOV of 38 times 38 arcmin(exp 2). The image shrinked as the off-axis angle increased. Above 13.5 arcmin of off-angle, a stray appeared around the image center in the off-axis direction. As for the on-axis image, a ring-shaped stray appeared at the edge of SXI of approximately 18 arcmin distant from the image center.

Soft X-ray spectrometer↗

Modeling of Inner Magnetosphere Coupling Processes

The Ring Current (RC) is the biggest energy player in the inner magnetosphere. It is the source of free energy for Electromagnetic Ion Cyclotron (EMIC) wave excitation provided by a temperature anisotropy of RC ions, which develops naturally during inward E B convection from the plasmasheet. The cold plasmasphere, which is under the strong influence of the magnetospheric electric field, strongly mediates the RC-EMIC wave-particle-coupling process and ultimately becomes part of the particle and energy interplay. On the other hand, there is a strong influence of the RC on the inner magnetospheric electric and magnetic field configurations and these configurations, in turn, are important to RC dynamics. Therefore, one of the biggest needs for inner magnetospheric research is the continued progression toward a coupled, interconnected system with the inclusion of nonlinear feedback mechanisms between the plasma populations, the electric and magnetic fields, and plasma waves. As we clearly demonstrated in our studies, EMIC waves strongly interact with electrons and ions of energies ranging from approx.1 eV to approx.10 MeV, and that these waves strongly affect the dynamics of resonant RC ions, thermal electrons and ions, and the outer RB relativistic electrons. As we found, the rate of ion and electron scattering/heating in the Earth's magnetosphere is not only controlled by the wave intensity-spatial-temporal distribution but also strongly depends on the spectral distribution of the wave power. The latter is also a function of the plasmaspheric heavy ion content, and the plasma density and temperature distributions along the magnetic field lines. The above discussion places RC-EMIC wave coupling dynamics in context with inner magnetospheric coupling processes and, ultimately, relates RC studies with plasmaspheric and Superthermal Electrons formation processes as well as with outer RB physics.

Khazanov, George V.↗

Analytic computation of energy derivatives - Relationships among partial derivatives of a variationally determined function

Formulas are presented relating Taylor series expansion coefficients of three functions of several variables, the energy of the trial wave function (W), the energy computed using the optimized variational wave function (E), and the response function (lambda), under certain conditions. Partial derivatives of lambda are obtained through solution of a recursive system of linear equations, and solution through order n yields derivatives of E through order 2n + 1, extending Puley's application of Wigner's 2n + 1 rule to partial derivatives in couple perturbation theory. An examination of numerical accuracy shows that the usual two-term second derivative formula is less stable than an alternative four-term formula, and that previous claims that energy derivatives are stationary properties of the wave function are fallacious. The results have application to quantum theoretical methods for the computation of derivative properties such as infrared frequencies and intensities.

King, H. F.↗

A Modified Wall Matching Treatment to Account for Local Solid to Fluid Thermal Coupling

The wall-matching methodology of Wilcox is modified to include a solid-wall, thermal-conduction model. This coupled fluid-thermal-structure model is derived assuming that the wall thermal-structure behavior is locally one-dimensional and that structural deformations, due to thermally induced stresses, are not significant. The one-dimensional coupled fluid-thermal-structure model is derived such that the wall temperature is removed as an independent boundary condition variable. The one-dimensional coupled fluid-thermal-structure model is also derived for the general case of an arbitrary mixture of thermally prefect gases and a wall of arbitrary thickness and conductivity by using a compressible, streamwise-pressure-gradient-corrected, wall-matching function and Fourier's law of heat conduction. The resulting model was implemented in the VULCAN CFD code as a new boundary condition type. VULCAN was then used to simulate a two-dimensional Mach 6 wind tunnel facility nozzle flow to demonstrate/validate the one-dimensional coupled fluid-thermal-structure model. The nozzle internal-wall surface temperature and heat transfer distributions computed using the one-dimensional coupled fluid-thermal-structure model are compared to wall temperature and heat transfer distributions from an iterative multi-dimensional analysis obtained by coupling the VULCAN CFD code and the MSC/NASTRAN-thermal code. The one-dimensional coupled fluid-thermal-structure model analysis is shown to be very robust and in excellent agreement with the multi-dimensional iteratively coupled analysis. It is also shown that the one-dimensional analysis can be used as an initial guess for the multi-dimensional iteratively coupled analysis.

White, Jeffery A.↗

A Modified Wall Matching Treatment to Account for Local Solid to Fluid Thermal Coupling

The wall-matching methodology of Wilcox is modified to include a solid-wall, thermal-conduction model. This coupled fluid-thermal-structure model is derived assuming that the wall thermal-structure behavior is locally one-dimensional and that structural deformations, due to thermally induced stresses, are not significant. The one-dimensional coupled fluid-thermal-structure model is derived such that the wall temperature is removed as an independent boundary condition variable. The one-dimensional coupled fluid-thermal-structure model is also derived for the general case of an arbitrary mixture of thermally prefect gases and a wall of arbitrary thickness and conductivity by using a compressible, streamwise-pressure-gradient-corrected, wall-matching function and Fourier's law of heat conduction. The resulting model was implemented in the VULCAN CFD code as a new boundary condition type. VULCAN was then used to simulate a two-dimensional Mach 6 wind tunnel facility nozzle flow to demonstrate/validate the one-dimensional coupled fluid-thermal-structure model. The nozzle internal-wall surface temperature and heat transfer distributions computed using the one-dimensional coupled fluid-thermal-structure model are compared to wall temperature and heat transfer distributions from an iterative multi-dimensional analysis obtained by coupling the VULCAN CFD code and the MSC/NASTRAN-thermal code. The one-dimensional coupled fluid-thermal-structure model analysis is shown to be very robust and in excellent agreement with the multi-dimensional iteratively coupled analysis. It is also shown that the one-dimensional analysis can be used as an initial guess for the multi-dimensional iteratively coupled analysis.

White, Jeffery A.↗

Incorporation of Failure 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 use in the aerospace and automotive communities. The aerospace community has identified several key capabilities which are currently lacking in the available material models in commercial transient dynamic finite element codes. To attempt to improve the predictive capability of composite impact simulations, a next generation material model is being developed for incorporation within the commercial transient dynamic finite element code LS-DYNA. The material model, which incorporates plasticity, damage and failure, utilizes experimentally based tabulated input to define the evolution of plasticity and damage and the initiation of failure as opposed to specifying discrete input parameters such as modulus and strength. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. For the damage model, a strain equivalent formulation is used to allow for the uncoupling of the deformation and damage analyses. In the damage model, a semi-coupled approach is employed where the overall damage in a particular coordinate direction is assumed to be a multiplicative combination of the damage in that direction resulting from the applied loads in various coordinate directions. For the failure model, a tabulated approach is utilized in which a stress or strain based invariant is defined as a function of the location of the current stress state in stress space to define the initiation of failure. Failure surfaces can be defined with any arbitrary shape, unlike traditional failure models where the mathematical functions used to define the failure surface impose a specific shape on the failure surface. In the current paper, the complete development of the failure model is described and the generation of a tabulated failure surface for a representative composite material is discussed.

Impact↗