Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Strain calculation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,135 records · Page 63

Elevated Temperature Deformation of Fe-39.8Al and Fe-15.6Mn-39.4Al

The elevated temperature compressive properties of binary Fe-39.8 at % Al and Fe-15.6Mn-39.4Al have been measured between 1000 and 1300 K at strain rates between 10(exp 7) and 10(exp 3)/ s. Although the Mn addition to iron aluminide did not change the basic deformation characteristics, the Mn-modified alloy was slightly weaker. In the regime where deformation of FeAl occurs by a high stress exponent mechanism (n = 6), strength increases as the grain size decreases at least for diameters between approx. 200 and approx. 10 microns. Due to the limitation in the grain size-flow stress-temperature-strain rate database, the influence of further reductions of the grain size on strength is uncertain. Based on the appearance of subgrains in deformed iron aluminide, the comparison of grain diameters to expected subgrain sizes, and the grain size exponent and stress exponent calculated from deformation experiments, it is believed that grain size strengthening is the result of an artificial limitation on subgrain size as proposed by Sherby, Klundt and Miller.

Whittenberger, J. Daniel↗

A Study of the Free Oscillations of the Earth

Published observations on the toroidal oscillations of the earth are critically reviewed. A supplementary analysis of the record obtained by the Lamont strain seismometer is presented. Eleven toroidal modes are identified, and it is concluded that the periods are known to within 1 per cent. A perturbation scheme involving the ratio of the angular velocity of the earth to the resonant frequency is used in calculating the effects due to the rotation of the earth on the resonant frequency. The free oscillations are viewed as a superposition of traveling waves. In a nonrotating system two traveling waves combine to produce a stationary standing wave. In a rotating system, the rotation distinguishes between waves that travel in the direction of rotation and those that travel in the opposite direction. Rotation removes a degeneracy and results in a splitting of a spectral peak of order l into 2 times l plus 1 peaks. The fractional displacement in frequency for the lowest-order toroidal oscillations is 1/206 and of the same order as the Q of the peak, so that splitting will probably not be observed in the toroidal oscillations. Viewed locally, rotation causes a particle to precess about a direction parallel to the axis of rotation. This precession will cause a variation of amplitude with time if the motion is recorded by an instrument with an anisotropic response function. Care is therefore needed in studying the time decay of a given spectral peak. Rotation also couples the normal coordinates so that a motion that is initially purely horizontal will develop a vertical component. It is expected that vertical seismometers should record particle motion with the toroidal frequencies. The perturbations of the toroidal oscillations due to core-mantle interaction are treated in detail. An exact expression is obtained for the rate of energy dissipated by a finitely conducting plate oscillating across a magnetic field. The energy dissipated at the core-mantle boundary due to viscous and hydromagnetic coupling is shown to be insignificant as compared with the energy dissipated within the mantle. The toroidal magnetic field leaking into the lower mantle combines with the dipole field, resulting in a stress on the mantle, tending to stiffen the lower boundary. The stress is of sufficient magnitude to produce a displacement toward higher frequency in the lower-order toroidal oscillations. Observations on the (sub 0) T (sub 2) oscillations lead to an estimate of the toroidal magnetic field in the lower mantle. A calculation of elastic energy in the low-order oscillations suggests a value of 10 (sup 18) ergs per cycles per hour for the energy density at low frequencies in the Chilean earthquake. Each mode of oscillation has a characteristic radial distribution of elastic energy associated with it. This distribution determines which parts of the earth contribute most heavily in determining a particular resonant frequency. The distribution of energy for the lower 17 modes for a homogeneous and a Gutenberg model earth is calculated. The resonant frequencies for models of the earth based on the Gutenberg and Lehmann distribution of elastic properties are presented. It is shown that the Gutenberg model earth fits the observations more closely than the Lehmann model and that a slight alteration of the Gutenberg model gives a significantly better fit to the observations. The alteration involves a lower shear-wave velocity in the lower mantle while the Gutenberg velocity distribution is maintained in the upper mantle. Various studies of the earth's oscillations coupled with surface-wave investigations substantiate Gutenberg's hypothesis of a layer of low velocity in the upper mantle. The physical conditions required for the formation of a region of low velocity are examined in detail. The results confirm Birch's earlier statement that a temperature gradient in excess of 6 degrees to 7 degrees per kilometer is needed to produce a decrease in velocity. The low-velocity layer does not require that the temperature approach or exceed the melting temperature. If tile upper mantle is homogeneous, the region of lower velocity should commence at the base of the crust and extend to 150 kilometers under the oceans and about 100 kilometers under continental regions. The distribution of thermal conductivity and radioactivity consistent with the low-velocity layer is also considered.

MacDonald, Gordon J. F.↗

Predicting Unreinforced Fabric Mechanical Behavior with Recurrent Neural Networks

Unreinforced woven fabrics are widely employed in various high-performance applications, including parachute deployment systems, airbags, and ballistic armor. The analysis of such materials is inherently complex due to the multiscale structure of these materials, and the dependence of macroscale behavior on changes that occur at lower scales. Previously, NASA’s Multiscale Analysis Tool (NASMAT) showed its capability in predicting unreinforced fabric behavior at the macroscale by capturing finite rotations that occur at the mesoscale. Though effective, the tool can face high computational cost for large, complex problems, motivating the need for the development of a surrogate model that can capture the same behavior. A recurrent neural network (RNN) was developed and trained on virtual NASMAT data to mimic the physics-based solutions while improving the computational runtime. The architecture of the RNN to best simulate the fabric behavior was carefully crafted based on heuristic knowledge of predicting physics-based temporal data, manual hyperparameter case studies, and Hyperband optimization.. The resultant model was able to predict a variety of stress-strain curves for fabrics with different mesoscale geometries, and was further validated by comparing to experimental data for the K706 style Kevlar plain-weave fabric, demonstrating the ability of the model to effectively capture the geometric changes in the fabric without explicitly calculating them, as is done in NASMAT. Furthermore, the tool showed its ability to improve on the runtime by a factor of 10 for fabric solutions compared to the multiscale tool, which would further enable the simulation of complex loading scenarios on unreinforced fabrics.

Fabric↗

Effect of chemical substitution and external strain on phase stability and ferroelectricity in two dimensional M 2 CT 2 MXenes

Two dimensional ferroelectric materials are gaining increasing attention for use in ultrathin electronic devices owing to the presence of a spontaneous polarization down to one or two monolayers. However, such materials are difficult to identify, especially those with out-of-plane electric polarizations. Previous work predicted that a metastable ferroelectric phase exists in the 2D MXene Sc 2 CO 2 , while further studies have predicted that this phase exists in other MXene chemistries. However, questions remain about the origin of ferroelectricity, the stability of this phase relative to other competing phases, and the effect of external stimuli in these materials. In this work, we use density functional theory calculations to investigate 12 M 2 CT 2 MXenes (M = transition metal, T = surface terminating group) and determine which have the ferroelectric phase as their ground state. We compute these materials’ polarizations, densities of states, phonon band structures, Bader charges, and Born effective charges in the ferroelectric phase to elucidate the reasons for its stabilization. We demonstrate that this ferroelectric phase can be preferentially stabilized in non-ferroelectric MXenes through full chemical substitution of Sc or O, alloying of the Sc sites, or application of epitaxial strain. Finally, we show that these materials have excellent piezoelectric properties as well. This work provides a detailed understanding of ferroelectric MXenes and show how the number of 2D ferroelectric materials can be increased through chemical substitution or application of external stimuli.

36 MATERIALS SCIENCE↗

Towards cubic symmetry for Ir 4 + : Structure and magnetism of the antifluorite K 2 IrBr 6

Here, crystal structure, electronic state of Ir 4 + , and magnetic properties of the antifluorite compound K 2 IrBr 6 are studied using high-resolution synchrotron x-ray diffraction, resonant inelastic x-ray scattering (RIXS), thermodynamic and transport measurements, and ab initio calculations. The crystal symmetry is reduced from cubic at room temperature to tetragonal below 170 K and eventually to monoclinic below 122 K. These changes are tracked by the evolution of the noncubic crystal-field splitting Δ measured by RIXS. Nonmonotonic changes in Δ are ascribed to the competing effects of the tilt, rotation, and deformation of the IrBr 6 octahedra as well as tetragonal strain on the electronic levels of Ir 4 + . The Néel temperature of T N = 11.9 K exceeds that of the isostructural K 2 IrCl 6 , and the magnitude of frustration on the fcc spin lattice decreases. We argue that the replacement of Cl by Br weakens electronic correlations and enhances magnetic couplings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Nucleation and growth of rolling contact failure of 440C bearing steel

A 'two-body' elasto-plastic finite element model of 2-dimensional rolling and rolling-plus-sliding was developed to treat the effect of surface irregularities. The model consists of a smooth cylinder in contact with a semi-infinite half-space that is either smooth or fitted with one of 0.4 microns deep or 7 microns deep groove, or a 0.4 microns high ridge-like asperity. The model incorporates elastic-linear-kinematic hardening-plastic (ELKP) and non-linear-kinematic hardening-plastic (NLKP) material constitutive relations appropriate for hardened bearing steel and the 440C grade. The calculated contact pressure distribution is Hertzian for smooth body contact, and it displays intense, stationary, pressure spikes superposed on the Hertzian pressure for contact with the grooved and ridged surface. The results obtained for the 0.4 microns deep groove compare well with those reported by Elsharkawy and Hamrock for an EHD lubricated contact. The effect of translating the counterface on the half space as opposed to indenting the half space with the counter face with no translation is studied. The stress and strain values near the surface are found to be similar for the two cases, whereas they are significantly different in the subsurface. It is seen that when tiny shoulders are introduced at the edge of the groove in the finite element model, the incremental plasticity and residual stresses are significantly higher in the vicinity of the right shoulder (rolling direction is from left to right) than at the left shoulder. This may explain the experimental observation that the spall nucleation occurs at the exit end of the artificially planted indents. Pure rolling calculations are compared with rolling + sliding calculations. For a coefficient of friction, mu = 0.1, the effect of friction is found to be small. Efforts were made to identify the material constitutive relations which best describe the deformation characteristics of the bearing steels in the initial few cycles. Elastic-linear-kinematic hardening-plastic (ELKP) material constitutive relations produce less net plastic deformation in the initial stages for a given stress, than seen in experiments. A new set of constitutive relations: non-linear-kinematic hardening-plastic (NLKP) was used. This material model produces more plasticity than the ELKP model and shows promise for treating the net distortions in the early stages. Techniques for performing experimental measurements that can be compared with the finite element calculations were devised. The measurements are being performed on 9mm-diameter, 440C steel cylindrical rolling elements in contact with 12.5 mm-diameter, 52100 steel balls in a 3-ball-rod fatigue test machine operating at 3600 RPM. Artificial, 7 microns deep, indents were inserted on the running track of the cylindrical rolling elements and profilometer measurements of these indents made, before and after the rolling. These preliminary measurements show that the indents are substantially deformed plastically in the process of rolling. The deformations of the groove calculated with the finite element model are comparable to those measured experimentally.

Gupta, V.↗

In-Situ Blade Strain Measurements of a Crossflow Turbine Operating in a Tidal Flow

This data was collected between October 25 and December 12 of 2022 at the University of New Hampshire (UNH) and Atlantic Marine Energy Center (AMEC) turbine deployment platform (TDP). The goal was to collect blade strain data from a crossflow turbine operating in a tidal flow. A table in ('Deployment Schedule.PNG') outlines the various dates when each instrument was operational, and more details can be found via literature listed in 'Related Publications'.txt. This dataset includes zipped folders for each instrument containing data in .csv files for the relevant duration specific to each instrument, along with separate README file for each measurement. Some instrument files are quite large and can pose a challenge for a visual spreadsheet editor to open. A processing software like MATLAB or Python is recommended. All data contained in this submission is unfiltered/unprocessed data unless otherwise noted in the README file. Blade strain was measured using 8 foil-based strain gauges along the span of a single turbine blade. Water currents were measured using Acoustic Doppler Current Profilers (ADCP's) and Acoustic Doppler Velocimeters (ADV's) both upstream and downstream of the turbine for inflow, wake and turbulence measurements. Electrical power output was measured using the Voltsys rectifier. Shaft speed was calculated based on the Voltsys measurements of the permanent magnet three phase generator AC generation frequency, coupled directly to the cross flow turbine under test (i.e., no gear box). Platform motions were captured using a Yost IMU (inertial measurement unit). Turbine thrust loading was measured using a reaction arm about the turbine deployment platform spanning beam, where two bi-directional load cells were connected to the system via a pinned connection. The TDP is a floating structure moored on the Portsmouth facing side of Memorial Bridge pier #2, which spans the Piscataqua River between Portsmouth, NH and Kittery, ME. The Piscataqua River connects the Great Bay Estuary to the Gulf of Maine and the river currents are dominated by tidal flow with water velocities exceeding 2.5 m/s during spring ebb tides at this site which were previously characterized by Chancey 2019. The turbine under test was a modified New Energy Corporation (Calgary, CA) model EVG-025 4-blade H-Darrius type vertical axis cross flow turbine that rotates in the clockwise direction with a rotor diameter of 3.2m and blade length of 1.7m. The hydro-foil profile was a NACA 0021 with a 10 inch chord length and a blade preset pitch angle of +4deg with a positive angle corresponding with the toe in direction. The standard EVG-025 has a rotor diameter of 3.4m and its rated power output is 25kW at 3 m/s. The rotor diameter was reduced to accommodate the size of the existing TDP moon-pool. A single blade of this turbine was further modified to accommodate 8 full-bridge strain gauges (Bharath et al 2023, Bichanich et al 2024). For power performance and other relevant details on the turbine and its characteristics, see O'Byrne 2022.

16 TIDAL AND WAVE POWER↗

Hybrid interatomic potential for Sn

To design materials for extreme applications, it is important to understand and predict phase transitions and their influence on material properties under high pressures and temperatures. Atomistic modeling can be a useful tool to assess these behaviors. However, this can be difficult due to the lack of fidelity of the interatomic potentials in reproducing this high pressure and temperature extreme behavior. Here, in this work, a hybrid EAM-R—which is the combination of embedded atom method (EAM) and rapid artificial neural network potential—for Tin (Sn) is described which is capable of accurately modeling the complex sequence of phase transitions between different metallic polymorphs as a function of pressure. This hybrid approach ensures that a basic empirical potential like EAM is used as a lower energy bound. By using the final activation function, the neural network contribution to energy must be positive, assuring stability over the whole configuration space. This implementation has the capacity to reproduce density functional theory results at 6 orders of magnitude slower than a pair potential for molecular dynamics simulation, including elastic and plastic characteristics and relative energies of each phase. Using calculations of the Gibbs free energy, it is demonstrated that the potential precisely predicts the experimentally observed phase changes at temperatures and pressures across the whole phase diagram. At 10.2 GPa, the present potential predicts a first-order phase transition between body-centered tetragonal (BCT) β-Sn and another polymorph of BCT-Sn. This structure transforms into body-centered cubic near the experimentally reported value at 33 GPa. Thus, the Sn potential developed in this paper can be used to study complex deformation mechanisms under extreme conditions of high pressure and strain rates unlike existing potentials. Moreover, the framework developed in this paper can be extended for different material systems with complex phase diagrams.

36 MATERIALS SCIENCE↗

Design of biphenylene-derived tunable dirac materials

The exploration of carbon allotropes has unveiled a series of two-dimensional (2D) materials with unique electronic and mechanical properties, yet the need for stable structures with tailored electronic properties persists. Here, in this study, we introduce a new class of 2D carbon allotropes derived from the biphenylene network (BPN), incorporating acetylenic linkages to tune their structural and electronic characteristics. Through density functional theory calculations, we identified ten novel BPN-derived structures that exhibit both energetic and dynamic stability, confirmed by cohesive energy and phonon spectrum analyses. Among them, BPN-02 and BPN04 are metallic, featuring critically-tilted type-III Dirac cones under ~ 5 % biaxial strain, while BPN-22 is a semiconductor with a band gap of 0.95 eV and exhibits highly anisotropic carrier mobility. Additionally, these structures demonstrate significant anisotropy in their elastic properties, further distinguishing them from other 2D carbon materials like graphene. Our findings suggest that these novel BPN-based structures have strong potential for next-generation electronic and optoelectronic applications, providing new avenues for the design and synthesis of advanced carbon materials.

2D carbon allotropes↗

Initial demonstration of automated fuel performance modeling with 1977 EBR-II metallic fuel pins using BISON code with FIPD and IMIS databases

Using the BISON fuel performance code, simulations were conducted using an automated process to read initial and operating conditions from the Fuels Irradiation and Physics Database (FIPD) and Integral Fast Reactor materials information system (IMIS) database, which contains metallic fuel data from the Experimental Breeder Reactor-II (EBR-II). This work demonstrates use of an integrated framework to access the vast majority of EBR-II experimental fuel pin data to support rapid development of fuel performance models for next-generation metallic fuel systems. With this capability, validation for fuel qualification can be performed rapidly. Between IMIS and FIPD, there is enough information to conduct 1977 unique EBR-II metallic fuel pin histories from 24 different experiments, at varying levels of detail between the two databases. Each of these histories includes a high-resolution power history, flux history, coolant channel flow rates, and coolant channel temperatures. Fission gas release (FGR), cumulative damage fraction (CDF), fuel axial swelling, cladding profilometry, and burnup were all simulated in BISON. The results were compared to post-irradiation examination (PIE) results for the initial demonstration of automated BISON modeling. BISON simulations conducted with IMIS and FIPD were in rough agreement with PIE measurements and calculations. Cladding profilometry, FGR, and fuel axial swelling were found to be in rough agreement with PIE measurements, depending on the physics used within the BISON input files. Here, the mechanical contact solver chosen was found to significantly impact axial fuel swelling and cladding strain predictions. CDF values were assessed to see whether pin failure may have been predicted (CDF ≥ 1). This work suggests that continued development of an automated tool for BISON should focus on inclusion of the Fast Flux Test Facility (FFTF) experimental data for a larger database for metallic fuel, improved physical models to better capture fuel performance, such as fuel-cladding interactions, and a more detailed comparison with available PIE data to further the BISON model development.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Operando TEM study of a working copper catalyst during ethylene oxidation

Active catalysts are typically metastable, and their surface state depends on the gas-phase chemical potential and reaction kinetics. To gain relevant insights into structure-performance relationships, it is essential to investigate catalysts under their operational conditions. Here, we use operando TEM combining real-time observations with online mass spectrometry (MS) to study a Cu catalyst during ethylene oxidation. We identify three distinct regimes characterized by varying structures and states that show different selectivities with temperature, and elucidate the reaction pathways with the aid of theoretical calculations. Our findings reveal that quasi-static Cu 2 O at low temperatures is selective towards ethylene oxide (EO) and acetaldehyde (AcH) via an oxometallacycle (OMC) pathway. In the dynamic Cu 0 /Cu 2 O oscillation regime at medium temperatures, partially reduced and strained oxides decrease the activation energies associated with partial oxidation. At high temperatures, the catalyst is predominantly Cu 0 , partially covered by a monolayer Cu 2 O. While Cu 0 is extremely efficient in dehydrogenation and eventual combustion, the monolayer oxide favors direct EO formation. These results challenge conclusions drawn from ultra-high vacuum studies that suggested metallic copper would be a selective epoxidation catalyst and highlight the need for operando study under realistic conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

X-ray diffraction at the National Ignition Facility

We report details of an experimental platform implemented at the National Ignition Facility to obtain in situ powder diffraction data from solids dynamically compressed to extreme pressures. Thin samples are sandwiched between tamper layers and ramp compressed using a gradual increase in the drive-laser irradiance. Pressure history in the sample is determined using high-precision velocimetry measurements. Up to two independently timed pulses of x rays are produced at or near the time of peak pressure by laser illumination of thin metal foils. The quasi-monochromatic x-ray pulses have a mean wavelength selectable between 0.6 Å and 1.9 Å depending on the foil material. The diffracted signal is recorded on image plates with a typical 2 θ x-ray scattering angle uncertainty of about 0.2° and resolution of about 1°. Analytic expressions are reported for systematic corrections to 2 θ due to finite pinhole size and sample offset. A new variant of a nonlinear background subtraction algorithm is described, which has been used to observe diffraction lines at signal-to-background ratios as low as a few percent. Variations in system response over the detector area are compensated in order to obtain accurate line intensities; this system response calculation includes a new analytic approximation for image-plate sensitivity as a function of photon energy and incident angle. This experimental platform has been used up to 2 TPa (20 Mbar) to determine the crystal structure, measure the density, and evaluate the strain-induced texturing of a variety of compressed samples spanning periods 2–7 on the periodic table.

47 OTHER INSTRUMENTATION↗

Single pair of Weyl nodes in the spin-canted structure of EuCd 2 As 2

We report time reversal symmetry breaking Weyl semimetals are unique among Weyl materials in allowing the minimal number of Weyl points, thus offering the clearest signatures of the associated physics. Here we present neutron diffraction, density-functional theory, and transport measurement results which indicate that EuCd 2 As 2 , under ambient field, strain, and pressure, is such a material with a single pair of Weyl points. Our work reveals a magnetic structure (magnetic space group C2'/m') with Eu moments pointing along the [210] direction in plane and canted ~ 30° out of plane. Density functional theory calculations using this structure show that the observed canting drastically alters the relevant electronic bands, relative to the in-plane order, leading to a single set of well-defined Weyl points. Furthermore, we find the canting angle can tune the distance of the Weyl points above the Fermi level, with the smallest distance at low canting angles. Finally, transport measurements of the anomalous Hall effect and longitudinal magnetoresistance exhibit properties indicative of a chiral anomaly, thus supporting the neutron scattering and DFT results suggesting EuCd 2 As 2 is close to the ideal situation of the Weyl hydrogen atom.

36 MATERIALS SCIENCE↗

Lamellar fluctuations melt ferroelectricity.

We consider a standard Ginzburg-Landau model of a ferroelectric whose electrical polarization is coupled to gradients of elastic strain. At the harmonic level, such flexoelectric interaction is known to hybridize acoustic and optic phonon modes and lead to phases with modulated lattice structures that precede the state with spontaneously broken inversion symmetry. Here, we use the self-consistent phonon approximation to calculate the effects of thermal and quantum polarization fluctuations on the bare hybridized modes to show that such long-range modulated order is unstable at all temperatures. We discuss the implications for the nearly ferroelectric SrTiO3 and KTaO3, and we propose that these systems are melted versions of an underlying modulated state that is dominated by nonzero momentum thermal fluctuations except at the very lowest temperatures.

Guzman-Verri, G. G.↗

Phase Demodulation by Frequency Chirping in Coherence Microwave Photonic Interferometry

This paper presents a signal processing method to demodulate the optical interference phase of cascaded individual optical fiber intrinsic Fabry-Perot interferometric (IFPI) sensors in a coherent microwave-photonic interferometry (CMPI) distributed sensing system. The new method utilizes the chirp effect of electro-optic modulator (EOM) to create a quasi-quadrature optical interference phase shift between two adjacent pulses which correspond to two adjacent reflection points in the time domain. The phase shift can be controlled by adjusting the bias voltage that is applied to the EOM. The interference phase is calculated by elliptically fitting the phase shift. The interference phase change is proportional to the optical path difference (OPD) change of the interferometer, and the sign can be used to differentiate the increase or decrease of the OPD. The method is demonstrated for distributed strain sensing, showing good linearity, high resolution and large dynamic range.

42 ENGINEERING↗

Variation of refractive index in strained In(x)Ga(1-x)As-GaAs heterostructures

In(x)Ga(1-x)As-GaAs heterostructures and strained-layer superlattices can be used as optical waveguides. For such applications it is important to know explicitly the refractive index variation with mismatch strain and with alloying in the ternary layer. Starting from the Kramers-Kronig integral dispersion relations, a model has been developed from which the refractive index change in the ternary layer of In(x)Ga(1-x)As-GaAs heterojunctions can be calculated. The results are presented and discussed. The expected changes in a superlattice have been qualitatively predicted.

Das, U.↗

Manufacturing stresses and strains in filament wound cylinders

Tests were performed to verify a previously developed model for simulating the manufacturing process of filament wound cylinders. The axial and hoop strains were measured during cure inside a filament wound Fiberite T300/976 graphite-epoxy cylinder. The measured strains were compared to those computed by the model. Good agreements were found between the data and the model, indicating that the model is a useful representation of the process. For the conditions of the test, the manufacturing stresses inside the cylinder were also calculated using the model.

Calius, E. P.↗

Probabilistic nonlinear finite element analysis of composite structures

A probabilistic finite element analysis procedure for laminated composite shells is developed. A total Lagrangian finite element formulation, employing a degenerated three-dimensional laminated composite shell element with the full Green-Lagrange strains and first-order shear deformable kinematics, is used. The first-order second-moment technique for probabilistic finite element analysis of random fields is employed, and results are presented in the form of mean and variance of the structural response. Reliability calculations are made by using the first-order reliability method combined with sensitivity derivatives from the finite element analysis. Both ply-level and micromechanics-level random variables are incorporated, the latter by means of the Aboudi micromechanics model. Two sample problems are solved to verify the accuracy of the procedures developed and to quantify the variability of certain material type/structure combinations. In general, the procedure is quite effective in determining the response statistics and reliability for linear and geometric nonlinear behavior of laminated composite shells.

Engelstad, S. P.↗