Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear structural analysis”

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.

1,203 records · Page 67
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Aerodynamic simulation on massively parallel systems

This paper briefly addresses the computational requirements for the analysis of complete configurations of aircraft and spacecraft currently under design to be used for advanced transportation in commercial applications as well as in space flight. The discussion clearly shows that massively parallel systems are the only alternative which is both cost effective and on the other hand can provide the necessary TeraFlops, needed to satisfy the narrow design margins of modern vehicles. It is assumed that the solution of the governing physical equations, i.e., the Navier-Stokes equations which may be complemented by chemistry and turbulence models, is done on multiblock grids. This technique is situated between the fully structured approach of classical boundary fitted grids and the fully unstructured tetrahedra grids. A fully structured grid best represents the flow physics, while the unstructured grid gives best geometrical flexibility. The multiblock grid employed is structured within a block, but completely unstructured on the block level. While a completely unstructured grid is not straightforward to parallelize, the above mentioned multiblock grid is inherently parallel, in particular for multiple instruction multiple datastream (MIMD) machines. In this paper guidelines are provided for setting up or modifying an existing sequential code so that a direct parallelization on a massively parallel system is possible. Results are presented for three parallel systems, namely the Intel hypercube, the Ncube hypercube, and the FPS 500 system. Some preliminary results for an 8K CM2 machine will also be mentioned. The code run is the two dimensional grid generation module of Grid, which is a general two dimensional and three dimensional grid generation code for complex geometries. A system of nonlinear Poisson equations is solved. This code is also a good testcase for complex fluid dynamics codes, since the same datastructures are used. All systems provided good speedups, but message passing MIMD systems seem to be best suited for large miltiblock applications.

Haeuser, Jochem↗
NASA NTRS2012

The Transition from Complex Crater to Peak-Ring Basin on the Moon: New Observations from the Lunar Orbiter Laser Altimeter (LOLA) Instrument

Impact craters on planetary bodies transition with increasing size from simple, to complex, to peak-ring basins and finally to multi-ring basins. Important to understanding the relationship between complex craters with central peaks and multi-ring basins is the analysis of protobasins (exhibiting a rim crest and interior ring plus a central peak) and peak-ring basins (exhibiting a rim crest and an interior ring). New data have permitted improved portrayal and classification of these transitional features on the Moon. We used new 128 pixel/degree gridded topographic data from the Lunar Orbiter Laser Altimeter (LOLA) instrument onboard the Lunar Reconnaissance Orbiter, combined with image mosaics, to conduct a survey of craters >50 km in diameter on the Moon and to update the existing catalogs of lunar peak-ring basins and protobasins. Our updated catalog includes 17 peak-ring basins (rim-crest diameters range from 207 km to 582 km, geometric mean = 343 km) and 3 protobasins (137-170 km, geometric mean = 157 km). Several basins inferred to be multi-ring basins in prior studies (Apollo, Moscoviense, Grimaldi, Freundlich-Sharonov, Coulomb-Sarton, and Korolev) are now classified as peak-ring basins due to their similarities with lunar peak-ring basin morphologies and absence of definitive topographic ring structures greater than two in number. We also include in our catalog 23 craters exhibiting small ring-like clusters of peaks (50-205 km, geometric mean = 81 km); one (Humboldt) exhibits a rim-crest diameter and an interior morphology that may be uniquely transitional to the process of forming peak rings. Comparisons of the predictions of models for the formation of peak-ring basins with the characteristics of the new basin catalog for the Moon suggest that formation and modification of an interior melt cavity and nonlinear scaling of impact melt volume with crater diameter provide important controls on the development of peak rings. In particular, a power-law model of growth of an interior melt cavity with increasing crater diameter is consistent with power-law fits to the peak-ring basin data for the Moon and Mercury. We suggest that the relationship between the depth of melting and depth of the transient cavity offers a plausible control on the onset diameter and subsequent development of peak-ring basins and also multi-ring basins, which is consistent with both planetary gravitational acceleration and mean impact velocity being important in determining the onset of basin morphological forms on the terrestrial planets.

Baker, David M. H.↗
NASA NTRS2006

Verification of Numerical Solutions for the Deployment of the Highly Nonlinear MARSIS Antenna Boom Lenticular Joints

The MARSIS antenna booms are constructed using lenticular hinges between straight boom segments in a novel design which allows the booms to be extremely lightweight while retaining a high stiffness and well defined structural properties once they are deployed. Lenticular hinges are elegant in form but are complicated to model as they deploy dynamically and require highly specialized nonlinear techniques founded on carefully measured mechanical properties. Results from component level testing were incorporated into a highly specialized ADAMS model which employed an automated damping algorithm to account for the discontinuous boom lengths formed during the deployment. Additional models with more limited capabilities were also developed in both DADS and ABAQUS to verify the ADAMS model computations and to help better define the numerical behavior of the models at the component and system levels. A careful comparison is made between the ADAMS and DADS models in a series of progressive steps in order to verify their numerical results. Different trade studies considered in the model development are outlined to demonstrate a suitable level of model fidelity. Some model sensitivities to various parameters are explored using subscale and full system models. Finally, some full system DADS models are exercised to illustrate the limitations of traditional modeling techniques for variable geometry systems which were overcome in the ADAMS model.

Mars↗
NASA NTRS2014

Buckling Testing and Analysis of Honeycomb Sandwich Panel Arc Segments of a Full-Scale Fairing Barrel: 8-ply Out-of-Autoclave Facesheets - Part 3

Four honeycomb sandwich panels, representing 1/16th arc segments of a 10 m diameter barrel section of the heavy lift launch vehicle, were manufactured under the NASA Composites for Exploration program and the NASA Constellation Ares V program. Two configurations were chosen for the panels: 6-ply facesheets with 1.125 in. honeycomb core and 8-ply facesheets with 1.000 in. honeycomb core. Additionally, two separate carbon fiber/epoxy material systems were chosen for the facesheets: inautoclave IM7/977-3 and out-of-autoclave T40-800B/5320-1. Smaller 3- by 5-ft panels were cut from the 1/16th barrel sections. These panels were tested under compressive loading at the NASA Langley Research Center. Furthermore, linear eigenvalue and geometrically nonlinear finite element analyses were performed to predict the compressive response of the 3- by 5-ft panels. This manuscript summarizes the experimental and analytical modeling efforts pertaining to the panel composed of 8-ply, T40-800B/5320-1 facesheets (referred to as Panel C). To improve the robustness of the geometrically nonlinear finite element model, measured surface imperfections were included in the geometry of the model. Both the linear and nonlinear, two-dimensional (2-D) and three-dimensional (3-D), models yield good qualitative and quantitative predictions. Additionally, it was predicted correctly that the panel would fail in buckling prior to failing in strength.

heavy lift launch vehicles↗
NASA NTRS2021

Application of a Developmental Composite Material Model to Predict the Crush Response of Two Energy Absorbers

In 2012, a consortium was formed with the goal of creating a new composite material model capable of predicting the wide range of properties, accumulated damage behavior, and the many different types of failure in composites under impact loading. The material model was developed for execution in the commercially available nonlinear, explicit transient dynamic finite element code, LS-DYNA®. This material model incorporates three sub-models: deformation, damage, and failure. In addition, the model accounts for strain rate and temperature effects and relies heavily on the input of tabulated material response data. The model is designated *MAT_COMPOSITE_TABULATED_PLASTICITY_DAMAGE, or *MAT_213. Initially, *MAT_213 was developed for use with solid elements only; however, a thin shell element formulation for *MAT_213 has been adapted recently. The objective of this project was to find a suitable modeling example to investigate the capabilities and performance of *MAT_213. In 2012, two composite energy absorbers were designed and evaluated at NASA Langley Research Center through multi-level testing and simulation. The first was a conical-shaped energy absorber, designated the conusoid, which consisted of four layers of hybrid carbon-Kevlar® plain-weave fabric oriented at [+45°/-45°/-45°/+45°] with respect to the vertical direction. The second was a sinusoidal-shaped energy absorber, designated the sinusoid, which consisted of hybrid carbon-Kevlar® plain-weave fabric face sheets, two layers for each face sheet oriented at ±45° with respect to the vertical direction, and a closed-cell ELFOAM® P200 polyisocyanurate foam core. Finite element models were developed of the energy absorbers and simulations were performed using LS-DYNA®. In this paper, the development of a *MAT_213 model of a hybrid carbon-Kevlar® plain-weave fabric is presented. Next, comparisons with material characterization tests are presented. Then, test-analysis results are documented for each energy absorber as comparisons of time-history responses, as well as predicted and experimental structural deformations and progressive damage under impact loading using the *MAT_213 material model. Since a prior *MAT_58, or *MAT_LAMINATED_COMPOSITE_FABRIC material model was used in previous simulations of the energy absorbers, comparisons are made between *MAT_58 and *MAT_213 model predictions with test data. Finally, the paper includes a comprehensive list of “lessons learned,” in which a series of parametric studies are documented that were performed to investigate specific issues related to the material model. These “lessons learned” are included in hopes that they may help future *MAT_213 users.

crashworthiness, LS-DYNA transient dynamic simulat↗
NASA NTRS2014

Buckling Testing and Analysis of Honeycomb Sandwich Panel Arc Segments of a Full-Scale Fairing Barrel: Comparison of In- and Out-of-Autoclave Facesheet Configurations

Four honeycomb sandwich panels, representing 1/16th arc segments of a 10-m diameter barrel section of the Heavy Lift Launch Vehicle, were manufactured and tested under the NASA Composites for Exploration and the NASA Constellation Ares V programs. Two configurations were chosen for the panels: 6-ply facesheets with 1.125 in. honeycomb core and 8-ply facesheets with 1.0 in. honeycomb core. Additionally, two separate carbon fiber/epoxy material systems were chosen for the facesheets: in-autoclave IM7/977-3 and out-of-autoclave T40-800b/5320-1. Smaller 3 ft. by 5 ft. panels were cut from the 1/16th barrel sections and tested under compressive loading. Furthermore, linear eigenvalue and geometrically nonlinear finite element analyses were performed to predict the compressive response of each 3 ft. by 5 ft. panel. To improve the robustness of the geometrically nonlinear finite element model, measured surface imperfections were included in the geometry of the model. Both the linear and nonlinear models yielded good qualitative and quantitative predictions. Additionally, it was correctly predicted that the panel would fail in buckling prior to failing in strength. Furthermore, several imperfection studies were performed to investigate the influence of geometric imperfections, fiber angle misalignments, and three-dimensional effects on the compressive response of the panel.

epoxy matrix composites↗
NASA NTRS2016

Development of Multi-Physics Dynamics Models for High-Frequency Large-Amplitude Structural Response Simulation

An analytic approach is demonstrated to reveal potential pyroshock-driven dynamic effects causing power losses in the Thermo-Electric (TE) module bars of the Mars Science Laboratory (MSL) Multi-Mission Radioisotope Thermoelectric Generator (MMRTG). This study utilizes high-fidelity finite element analysis with SIERRA/PRESTO codes to estimate wave propagation effects due to large-amplitude suddenly-applied pyro shock loads in the MMRTG. A high fidelity model of the TE module bar was created with approximately 30 million degrees-of-freedom (DOF). First, a quasi-static preload was applied on top of the TE module bar, then transient tri-axial acceleration inputs were simultaneously applied on the preloaded module. The applied input acceleration signals were measured during MMRTG shock qualification tests performed at the Jet Propulsion Laboratory. An explicit finite element solver in the SIERRA/PRESTO computational environment, along with a 3000 processor parallel super -computing framework at NASA AMES, was used for the simulation. The simulation results were investigated both qualitatively and quantitatively. The predicted shock wave propagation results provide detailed structural responses throughout the TE module bar, and key insights into the dynamic response (i.e., loads, displacements, accelerations) of critical internal spring/piston compression systems, TE materials, and internal component interfaces in the MMRTG TE module bar. They also provide confidence on the viability of this high-fidelity modeling scheme to accurately predict shock wave propagation patterns within complex structures. This analytic approach is envisioned for modeling shock sensitive hardware susceptible to intense shock environments positioned near shock separation devices in modern space vehicles and systems.

Shock Waves↗
NASA NTRS1999

A fluid--structure interaction finite element analysis of pulsatile blood flow through a compliant stenotic artery

A new model is used to analyze the fully coupled problem of pulsatile blood flow through a compliant, axisymmetric stenotic artery using the finite element method. The model uses large displacement and large strain theory for the solid, and the full Navier-Stokes equations for the fluid. The effect of increasing area reduction on fluid dynamic and structural stresses is presented. Results show that pressure drop, peak wall shear stress, and maximum principal stress in the lesion all increase dramatically as the area reduction in the stenosis is increased from 51 to 89 percent. Further reductions in stenosis cross-sectional area, however, produce relatively little additional change in these parameters due to a concomitant reduction in flow rate caused by the losses in the constriction. Inner wall hoop stretch amplitude just distal to the stenosis also increases with increasing stenosis severity, as downstream pressures are reduced to a physiological minimum. The contraction of the artery distal to the stenosis generates a significant compressive stress on the downstream shoulder of the lesion. Dynamic narrowing of the stenosis is also seen, further augmenting area constriction at times of peak flow. Pressure drop results are found to compare well to an experimentally based theoretical curve, despite the assumption of laminar flow.

Non-NASA Center↗
NASA NTRS1993

Internal strain analysis of ceramics using scanning laser acoustic microscopy

Quantitative studies of material behavior characteristics are essential for predicting the functionality of a material under its operating conditions. A nonintrusive methodology for measuring the in situ strain of small dimeter (to 11 microns) ceramic fibers under uniaxial tensile loading and the local internal strains of ceramics and ceramic composites under flexural loading is introduced. The strain measurements and experimentally observed mechanical behavior are analyzed in terms of the microstructural development and fracture behavior of each test specimen evaluated. Measurement and analysis of Nicalon silicon carbide (SiC) fiber (15 microns diameter) indicate that the mean elastic modulus of the individual fiber is 185.3 GPa. Deviations observed in the experimentally determined elastic modulus values between specimens were attributed to microstructural variations which occur during processing. Corresponding variations in the fracture surface morphology were also observed. The observed local mechanical behavior of a lithium alumino-silicate (LAS) glass ceramic, a LAS/SiC monofilament composite, and a calcium alumino-silicate (CAS)/SiC fully reinforced composite exhibits nonlinearities and apparent hysteresis due to the subcritical mechanical loading. Local hysteresis in the LAS matrices coincided with the occurrence of multiple fracture initiation sites, localized microcracking, and secondary cracking. The observed microcracking phenomenon was attributed to stress relaxation of residual stresses developed during processing, and local interaction of the crack front with the microstructure. The relaxation strain and stress predicted on apparent mechanical hysteresis effects were defined and correlated with the magnitude of the measured fracture stress for each specimen studied. This quantitative correlation indicated a repeatable measure of the stress at which matrix microcracking occurred for stress relief of each material system. Stress relaxation occurred prior to the onset of steady state cracking conditions. The relaxation stress occurred at 18.5 percent of the fracture stress in LAS and 11.0 percent of the yield stress in CAS/SiC. The relaxation stress ratio was dependent upon the dominant fracture mode of the LAS/SiC specimens. Relaxation stress ratios greater than 0.30 were observed for specimens which fractured due to shear at the fiber matrix interface; specimens which fracture due to tensile cracking had relaxation stress ratios less than 0.30. The stress relaxation ratio appeared to be a specific characteristic of the glass ceramic material. The measured stress relaxation for LAS indicated a measure of the inherent residual stresses in the material due to processing and suggested localized toughening mechanisms for brittle material structures.

Kent, Renee M.↗
NASA NTRS1994

Strong gravitational lensing statistics as a test of cosmogonic scenarios

Gravitational lensing statistics can provide a direct and powerful test of cosmic structure formation theories. Since lensing tests, directly, the magnitude of the nonlinear mass density fluctuations on lines of sight to distant objects, no issues of 'bias' (of mass fluctuations with respect to galaxy density fluctuations) exist here, although lensing observations provide their own ambiguities of interpretation. We develop numerical techniques for generating model density distributions with the very large spatial dynamic range required by lensing considerations and for identifying regions of the simulations capable of multiple image lensing in a conservative and computationally efficient way that should be accurate for splittings significantly larger than 3 seconds. Applying these techniques to existing standard Cold dark matter (CDM) (Omega = 1) and Primeval Baryon Isocurvature (PBI) (Omega = 0.2) simulations (normalized to the Cosmic Background Explorer Satellite (COBE) amplitude), we find that the CDM model predicts large splitting (greater than 8 seconds) lensing events roughly an order-of-magnitude more frequently than the PBI model. Under the reasonable but idealized assumption that lensing structrues can be modeled as singular isothermal spheres (SIS), the predictions can be directly compared to observations of lensing events in quasar samples. Several large splitting (Delta Theta is greater than 8 seconds) cases are predicted in the standard CDM model (the exact number being dependent on the treatment of amplification bias), whereas none is observed. In a formal sense, the comparison excludes the CDM model at high confidence (essentially for the same reason that CDM predicts excessive small-scale cosmic velocity dispersions.) A very rough assessment of low-density but flat CDM model (Omega = 0.3, Lambda/3H(sup 2 sub 0) = 0.7) indicates a far lower and probably acceptable level of lensing. The PBI model is consistent with, but not strongly tested by, the available lensing data, and other open models would presumably do as well as PBI. These preliminary conclusions and the assumptions on which they are based can be tested and the analysis can be applied to other cosmogonic models by straightforward extension of the work presented here.

Cen, Renyue↗
NASA NTRS2006

Nonlinear Spring Finite Elements for Predicting Mode I-Dominated Delamination Growth in Laminated Structure with Through-Thickness reinforcement

One particular concern of polymer matrix composite laminates is the relatively low resistance to delamination cracking, in particular when the dominant type of failure is mode I opening. One method proposed for alleviating this problem involves the insertion pultruded carbon pins through the laminate thickness. The pins, known as z-pins, are inserted into the prepreg laminate using an ultrasonic hammer prior to the curing process, resulting in a field of pins embedded normal to the laminate plane as illustrated in Figure. 1. Pin diameters range between 0.28-mm to 0.5-mm and standard areal densities range from 0.5% to 4%. The z-pins are provided by the manufacturer, Aztex(Registered TradeMark) , in a low-density foam preform, which acts to stabilize orientation of the pins during the insertion process [1-3]. Typical pin materials include boron and carbon fibers embedded in a polymer matrix. A number of methods have been developed for predicting delamination growth in laminates reinforced with z-pins. During a study on the effect of z-pin reinforcement on mode I delamination resistance, finite element analyses of z-pin reinforced double cantilever beam (DCB) specimens were performed by Cartie and Partridge [4]. The z-pin bridging stresses were modeled by applying equivalent forces at the pin locations. Single z-pin pull-out tests were performed to characterize the traction law of the pins under mode I loading conditions. Analytical solutions for delamination growth in z-pin reinforced DCB specimens were independently derived by Robinson and Das [5] and Ratcliffe and O'Brien [6]. In the former case, pin bridging stresses were modeled using a distributed load and in the latter example the bridging stresses were discretely modeled by way of grounded springs. Additionally, Robinson and Das developed a data reduction strategy for calculating mode I fracture toughness, G(sub Ic), from a z-pin reinforced DCB specimen test [5]. In both cases a traction law similar to that adopted by Cartie and Partridge was used to represent z-pin failure under mode I loading conditions. In the current work spring elements available in most commercial finite element codes were used to model z-pins. The traction law used in previous analyses [4-6] was employed to represent z-pin damage. This method is intended for and is limited to simulating z-pins in composite laminate structure containing mode I-dominated delamination cracking. The current technique differs from previous analyses in that spring finite elements (available in commercial codes) are employed for simulating zpins, reducing the complexity of the analysis construction process. Furthermore, the analysis method can be applied to general structure that experiences mode I-dominated delamination cracking, in contrast to existing analytical solutions that are only applicable to coupon DCB specimens.

Ratcliffe, James G.↗
NASA NTRSDate not supplied

Nonlinear Buckling Analysis of Tow-Steered Composite Cylinders with Cutouts

The buckling and postbuckling behavior of two composite tow-steered shells with cutouts of different sizes is assessed using nonlinear finite element (FE) analysis and compared to experimental measurements. The cylindrical shells are manufactured using an automated fiber placement system, where the shells’ fiber orientation angles vary continuously around the shell circumference from ±10 degrees on the crown and keel to ±45 degrees on the sides. One shell features thickness variations due to tow overlaps that result from application of all 24 tows during each pass of the fiber placement system. The second shell uses the system’s tow drop/add capability to achieve a more uniform wall thickness without overlaps. Unreinforced cutouts of two different sizes—the first smaller cutout representing a passenger door on a commercial aircraft and the second larger cutout a cargo door—were machined into each of the two cylinders resulting in a total of four test cases. These cylinders were tested in axial compression and buckled elastically in previous work and are now analyzed using nonlinear FE models to compare bifurcation buckling loads as well as the load-displacement response in the prebuckling and postbuckling regimes. For all four shells analyzed, the prebuckling stiffness, buckling load, and deformation mode sequence throughout the loading-unloading cycle is accurately reproduced by the models. In particular, the shells first buckle locally around the cutouts in a stable (super-critical) manner with only a slight decrease in axial stiffness, which occurs due to the favorable load redistribution facilitated by tow steering. The shells then buckle globally in an unstable (sub-critical) manner with diamond-shaped buckles forming to the left and right of the cutouts. The buckling load of all shells with cutouts is at least 82% of the buckling load of the pristine shells without cutouts. Overall, the ability to sustain local buckling phenomena, and the relatively small reductions in global buckling load compared to pristine shells without cutouts, demonstrates the great potential of using tow steering to mitigate the adverse effects of cutouts in axially-compressed shell structures.

composites↗
NASA NTRSDate not supplied

Nonlinear Buckling Analysis of Tow-Steered Composite Cylinders with Cutouts

The buckling and postbuckling behavior of two composite tow-steered shells with cutouts of different sizes is assessed using nonlinear finite element (FE) analysis and compared to experimental measurements. The cylindrical shells are manufactured using an automated fiber placement system, where the shells’ fiber orientation angles vary continuously around the shell circumference from 10 degrees on the crown and keel to 45 degrees on the sides. One shell features thickness variations due to tow overlaps that result from application of all 24 tows during each pass of the fiber placement system. The second shell uses the system’s tow drop/add capability to achieve a more uniform wall thickness without overlaps. Unreinforced cutouts of two different sizes—the first smaller cutout representing a passenger door on a commercial aircraft and the second larger cutout a cargo door—were machined into each of the two cylinders resulting in a total of four test cases. These cylinders were tested in axial compression and buckled elastically in previous work and are now analyzed using nonlinear FE models to compare bifurcation buckling loads as well as the load-displacement response in the prebuckling and postbuckling regimes. For all four shells analyzed, the prebuckling stiffness, buckling load, and deformation mode sequence throughout the loading-unloading cycle is accurately reproduced by the models. In particular, the shells first buckle locally around the cutouts in a stable (super-critical) manner with only a slight decrease in axial stiffness, which occurs due to the favorable load redistribution facilitated by tow steering. The shells then buckle globally in an unstable (sub-critical) manner with diamond-shaped buckles forming to the left and right of the cutouts. The buckling load of all shells with cutouts is at least 82% of the buckling load of the pristine shells without cutouts. Overall, the ability to sustain local buckling phenomena, and the relatively small reductions in global buckling load compared to pristine shells without cutouts, demonstrates the great potential of using tow steering to mitigate the adverse effects of cutouts in axially-compressed shell structures.

Buckling↗
NASA NTRS2021

Micromechanics-based Modeling of Laminated SiC/SiC Ceramic Matrix Composites

The behavior and response of ceramic matrix composites (CMCs), in particular silicon carbide fiber reinforced silicon carbide matrix (SiC/SiC), is affected by many factors such as variation of fiber volume fraction, residual stresses resulting from processing of the composites at high temperature, random microstructures, and the presence of matrix flaws (e.g., voids, pores, cracks etc.) as well as general material nonlinearity and heterogeneity that occurs randomly in a composite. Residual stresses arising from the phase change of constituents are evaluated in this paper and it is shown that they do influence composite strength and need to be properly accounted for. Additionally, the microstructures (location of fiber centers, coating thickness etc.) of advanced CMCs are usually disordered (or random) and fiber diameter and strength typically have a distribution. They rarely resemble the ordered fiber packing (square, rectangular, or hexagonal) that is generally assumed in micromechanics-based models with periodic boundary conditions for computational expediency. These issues raise the question of how should one model such systems effectively? Can an ordered hexagonal packed repeating unit cell (RUC) accurately represent the random microstructure behavior? How many fibers need to be included to enable accurate representation? Clearly, the number of fibers within an RUC must be limited to insure a balance between accuracy and efficiency. NASA’s in-house micromechanics-based code MAC/GMC provides a framework to analyze such RUCs for the overall composite behavior and the FEAMAC computer code provides linkage of MAC/GMC to the commercial FEA code, ABAQUS. The appropriate level of discretization of the RUC as well as the analysis method employed, i.e., Generalized Method of Cells (GMC) or High Fidelity Generalized Method of Cells (HFGMC), is investigated in this paper in the context of a unidirectional as well as a cross-ply laminated CMC. Results including effective composite properties, proportional limit stress (an important design parameter) and fatigue are shown utilizing both GMC as well as HFGMC. Finally, a few multiscale analyses are performed on smooth bar test coupons as well as test coupons with features such as open-hole and double notches using FEAMAC. Best practices and guidance are provided to take these phenomena into account and keep a proper balance between fidelity (accuracy) and efficiency. Following these guidelines can account for important physics of the problem and provide significant advantages when performing large multiscale composite structural analyses.

Ceramic Matrix Composites↗
NASA NTRS2019

In Situ Observation of Growth Dynamics in DECLIC Directional Solidification Insert Onboard ISS: DSI-R Flight Campaign

The study of solidification microstructure formation is of utmost importance for materials design and processing, as solid-liquid interface patterns largely govern mechanical and physical properties. Pattern selection occurs under dynamic conditions of growth in which the initial morphological instability evolves nonlinearly and undergoes a reorganization process. The dynamic and nonlinear nature of this instability renders in situ observation of the interface an invaluable tool to gain knowledge on the time-evolution of the interface pattern. Transparent organic analogs, which solidify like metallic alloys, allow direct visualization of interface dynamics. Extensive ground-based studies of both metallic and organic bulk samples have established the presence of significant convection during solidification processes that alters the formation of interfacial microstructures. A reduced-gravity environment is therefore mandatory for fluid flow elimination in bulk samples. In the framework of the CNES project MISOL3D (MIcrostrutures de SOLidification 3D) and the NASA projects DSIP (Dynamical Selection of 3D Interface Patterns), SPADES (SPAtiotemporal Evolution of three-dimensional DEndritic array Structures) and CAMUS (ComputAtional Studies of MicrostrUcture Formation During Alloy Solidification in Microgravity), we participated in the development of the Directional Solidification Insert (DSI) of the DEvice for the study of Critical Liquids and Crystallization (DECLIC). The DECLIC-DSI is dedicated to in situ and real-time characterization of solid-liquid interface patterns during directional solidification of transparent alloys in diffusive transport regime. Between April 2010 and March 2011, the first ISS campaign (DSI) explored the entire range of microstructures resulting in unprecedented observations. A second campaign (DSI-R), performed between October 2017 and December 2018, in which the insert contained an alloy of higher solute concentration, allowed to complete the benchmark database. The increase of solute concentration resulted in well-developed dendritic patterns at lower velocities (lower interface curvature and larger tip radius). The microstructure resulting from dendritic growth is dominant in metallurgy so that it is fundamental to understand the mechanisms of its formation. The main aims of this experimental campaign are to understand: the mechanisms of the cell to dendrite transition, the fundamental mechanisms of sidebranching formation, the dependence of dendrite tip shapes on growth conditions, the interaction of primary array and secondary sidebranches, and the influence of subgrain boundaries on the spatiotemporal organization of the array structure. Preparation, analysis and interpretation of the experiments performed onboard ISS are considerably enhanced by experiments performed on ground using thin-samples (Pr. Trivedi’s group) and phase-field simulations of microstructure formation in a diffuse growth regime (Pr. Karma’s group). In this summary, we will present an initial assessment of the results obtained during the DSI-R campaign. Nomenclature Vp: pulling velocity G: thermal gradient L: solidified length rtip: tip curvature radius Acronyms/Abbreviations CADMOS: Centre d'Aide au Développement des activités en Micro-pesanteur et des Opérations Spatiales CNES: Centre National d’Études Spatiales DECLIC: Device for the study of Critical LIquids and Crystallization DSI: Directional Solidification Insert DSI-R: Directional Solidification Insert - Refurbish ISS: International Space Station NASA: National Aeronautics and Space Administration PF: phase-field SCN: Succinonitrile 3D: three-dimensional

Mota, F.L.↗