Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ELASTIC SHELL”

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.

498 records · Page 28

Impact Simulation of Composite Panels for Aerospace Applications

The aerospace research community has been working towards the development of a comprehensive framework for the design of composite structures at the component, subassembly, assembly, and system levels. Such a framework requires accurate, efficient and easy-to-use mathematical and computational models that are in sync with physical reality. The research work in this paper is designed to help an aerospace structural engineer move towards this goal by simulating dynamic impact to composite panels. First, the existing orthotropic visco-elastic visco-plastic material model in LS-DYNA commercial finite element program, MAT_COMPOSITE_TABULATED_PLASTICITY_DAMAGE (MAT_213) is enhanced with the development and use of rate-dependent damage data. Using state-of-the-art experimental procedures, impact validation test data are created. Second, the versatility of the constitutive model is examined by its ability to simulate a low velocity (27 ft/s, 8 m/s) impact test, and two high velocity impact tests (385–417 ft/s, 117–127 m/s) that straddle the ballistic limit of the projectile and the composite panels. Third, these impact events are modeled using three different finite element types – under-integrated and fully-integrated solid elements, and fully-integrated thin shell elements. Insight into accuracy versus numerical efficiency is gained by comparing the performance of solid versus thin shell elements. The numerical results show that accurate predictions of these impact events can be obtained with reasonable effort and care.

Simulations↗

A general low frequency acoustic radiation capability for NASTRAN

A new capability called NASHUA is described for calculating the radiated acoustic sound pressure field exterior to a harmonically-excited arbitrary submerged 3-D elastic structure. The surface fluid pressures and velocities are first calculated by coupling a NASTRAN finite element model of the structure with a discretized form of the Helmholtz surface integral equation for the exterior fluid. After the fluid impedance is calculated, most of the required matrix operations are performed using the general matrix manipulation package (DMAP) available in NASTRAN. Far field radiated pressures are then calculated from the surface solution using the Helmholtz exterior integral equation. Other output quantities include the maximum sound pressure levels in each of the three coordinate planes, the rms and average surface pressures and normal velocities, the total radiated power and the radiation efficiency. The overall approach is illustrated and validated using known analytic solutions for submerged spherical shells subjected to both uniform and nonuniform applied loads.

Everstine, G. C.↗

Hydro-isostatic deflection and tectonic tilting in the central Andes: Initial results of a GPS survey of Lake Minchin shorelines

Sufficiently large lake loads provide a means of probing rheological stratification of the crust and upper mantle. Lake Minchin was the largest of the late Pleistocene pluvial lakes in the central Andes. Prominent shorelines, which formed during temporary still-stands in the climatically driven lake level history, preserve records of lateral variations in subsequent net vertical motions. At its maximum extent the lake was 140 m deep and spanned 400 km N-S and 200 km E-R. The load of surficial water contained in Lake Minchin was sufficient to depress the crust and underlying mantle by 20-40 m, depending on the subjacent rheology. Any other differential vertical motions will also be recorded as departures from horizontality of the shorelines. We recently conducted a survey of shoreline elevations of Lake Minchin with the express intent of monitoring the hydro-isostatic deflection and tectonic tilting. Using real-time differential Global Positioning System (GPS), we measured topographic profiles across suites of shorelines at 15 widely separated locations throughout the basin. Horizontal and vertical accuracies attained are roughly 30 and 70 cm, respectively. Geomorphic evidence suggests that the highest shoreline was occupied only briefly (probably less than 200 years) and radiocarbon dates on gastropod shells found in association with the shore deposits constrain the age to roughly 17 kyr. The basin-side pattern of elevations of the highest shoreline is composed of two distinct signals: (27 +/- 1) m of hydro-isostatic deflection due to the lake load, and a planar tilt with east and north components of (6.8 +/- 0.4) 10(exp -5) and 9-5.3 +/- 0.3) 10(exp -5). This rate of tilting is too high to be plausibly attributed to steady tectonism, and presumably reflects some unresolved combination of tectonism plus the effects of oceanic and lacustrine loads on a laterally heterogeneous substrate. The history of lake level fluctuations is still inadequately known to allow detailed inferences of crust and mantle rheology. However, it is already clear that the effective elastic plate thickness is closer to 40 km than the 60-70 km crustal thickness in the central Andes and the effective viscosity is less than 5 10(exp 20) Pa s.

Bills, Bruce G.↗

Formulation and Performance Evaluation of Epoxy Sealant Systems for Double-Shell Tank Bottom Refurbishment

The performance of epoxy sealants used in the refurbishment of double-shell tank (DST) systems requires balancing processability, thermomechanical stability, and adhesion to cementitious substrates. This study incorporates Heloxy 8 as a reactive diluent into Westlake 862 epoxy to tailor workability and cured-state properties. Rheological time-sweep analysis demonstrates that increasing the diluent content significantly reduces complex viscosity and extends workability, thereby improving pumpability and flow for large-area applications. However, the targeted 2-hour processing window is not fully achieved. Differential scanning calorimetry (DSC) confirms that all formulations cure at room temperature to glass transition temperatures ( T g ) at least 20 °C above the maximum DST operating temperature (27 °C), thereby ensuring service in the glassy regime. Dynamic mechanical analysis (DMA) reveals formulation-dependent reductions in tan delta and increases in storage modulus, indicating increasingly elastic and mechanically stable networks with diluent incorporation. Pull-off adhesion testing shows that modified formulations (70–90% Westlake epoxy) exhibit significantly higher adhesion strengths than the unmodified system. Grout cohesive failure indicates that interfacial bonding exceeds substrate strength. Collectively, these results demonstrate that controlled reactive diluent incorporation enables optimization of processing behavior, interfacial adhesion, and thermomechanical performance, supporting the suitability of the modified epoxy systems as durable sealant layers for cementitious barrier applications in hazardous waste containment infrastructure.

Differential scanning calorimetry↗

Coupled BE/FE/BE approach for scattering from fluid-filled structures

NASHUA is a coupled finite element/boundary element capability built around NASTRAN for calculating the low frequency far-field acoustic pressure field radiated or scattered by an arbitrary, submerged, three-dimensional, elastic structure subjected to either internal time-harmonic mechanical loads or external time-harmonic incident loadings. Described here are the formulation and use of NASHUA for solving such structural acoustics problems when the structure is fluid-filled. NASTRAN is used to generate the structural finite element model and to perform most of the required matrix operations. Both fluid domains are modeled using the boundary element capability in NASHUA, whose matrix formulation (and the associated NASTRAN DMAP) for evacuated structures can be used with suitable interpretation of the matrix definitions. After computing surface pressures and normal velocities, far-field pressures are evaluated using an asymptotic form of the Helmholtz exterior integral equation. The proposed numerical approach is validated by comparing the acoustic field scattered from a submerged fluid-filled spherical thin shell to that obtained with a series solution, which is also derived here.

Everstine, Gordon C.↗

Hierarchic models for laminated plates

Structural plates and shells are three-dimensional bodies, one dimension of which happens to be much smaller than the other two. Thus, the quality of a plate or shell model must be judged on the basis of how well its exact solution approximates the corresponding three-dimensional problem. Of course, the exact solution depends not only on the choice of the model but also on the topology, material properties, loading and constraints. The desired degree of approximation depends on the analyst's goals in performing the analysis. For these reasons models have to be chosen adaptively. Hierarchic sequences of models make adaptive selection of the model which is best suited for the purposes of a particular analysis possible. The principles governing the formulation of hierarchic models for laminated plates are presented. The essential features of the hierarchic models described models are: (1) the exact solutions corresponding to the hierarchic sequence of models converge to the exact solution of the corresponding problem of elasticity for a fixed laminate thickness; and (2) the exact solution of each model converges to the same limit as the exact solution of the corresponding problem of elasticity with respect to the laminate thickness approaching zero. The formulation is based on one parameter (beta) which characterizes the hierarchic sequence of models, and a set of constants whose influence was assessed by a numerical sensitivity study. The recommended selection of these constants results in the number of fields increasing by three for each increment in the power of beta. Numerical examples analyzed with the proposed sequence of models are included and good correlation with the reference solutions was found. Results were obtained for laminated strips (plates in cylindrical bending) and for square and rectangular plates with uniform loading and with homogeneous boundary conditions. Cross-ply and angle-ply laminates were evaluated and the results compared with those of MSC/PROBE. Hierarchic models make the computation of any engineering data possible to an arbitrary level of precision within the framework of the theory of elasticity.

Actis, Ricardo Luis↗

Analysis of Element Size and Element Type on SS304L Material Model Performance

Finite element solutions are only valid if the information inputted into the models is a valid assumption. This is why it is important to determine if the model behaves as expected especially in the material response. Material models that consider deformation in the elastic and plastic regions are desirable. As a material deforms plastically, hardening will occur until the ultimate strength is reached due to the reduction of the cross sectional area in a uniaxial tensile test environment. A phenomenological approximation of this behavior which is known as power law hardening (plh) can be expressed as σ = σ y + K $ϵ^{n}_{p}$ where K and n are empirically determined constants, and σ y is the proportional limit of the material. The first objective of this study is to verify the behavior for a calibrated plh material model for selective deviatoric hexahedral elements when applied to a different material stock. The Sierra SM input for this calibrated model can be found in the appendix. The second objective of this study is to determine the differences in failure characteristics for considerably large shell elements versus that of fine hexahedral elements.

36 MATERIALS SCIENCE↗

Multimaterial 3D Printing in Activating Bath Enables In Situ Polymerization of Thermosets with Intricate Geometries and Diverse Elastic Behaviors

Polydicyclopentadiene, p(DCPD), is a high‐performance thermoset valued for its exceptional toughness, strength, and stiffness. When copolymerized with 1,5‐cyclooctadiene (COD), its mechanical properties can be tuned from glassy to rubbery at room temperature. While frontal polymerization enables a rapid and energy‐efficient route to 3D print DCPD‐based materials, challenges such as ink shelf life and gravitational distortion, especially in direct ink writing of soft COD‐rich formulations, must be considered. Here, a complementary chemical strategy is presented, embedded 3D printing, that enables localized in situ polymerization of printed DCPD/COD inks within a reactive support matrix. The matrix provides both physical support and a reservoir of chemical activator, which diffuses into the ink, activates a latent bis(N‐heterocyclic carbene) Ru precatalyst, and initiates ring‐opening metathesis polymerization. Curing begins at the ink–matrix interface and propagates inward via diffusion, stabilizing the interface and preventing capillary‐driven deformation regardless of the matrix yield stress. This approach eliminates the need for cold storage, external curing, or photoinitiation, significantly expanding the processing window. Using this method, diverse thermosetting and elastomeric architectures are fabricated with features as small as 5 µm and aspect ratios of 100, including interlinked chains, shallow spherical shells exhibiting snap‐through buckling, and hair‐like fin arrays inaccessible through traditional techniques.

chemical activation↗

Acoustic Communication in Deep Ice at Ocean Worlds

One of NASA priorities is the in-situ exploration of ocean worlds in the solar system where potentially there might be life under the ice shell. This requires reaching the ocean below extremely cold through significant deep ice. Jupiter’s moon, Europa, is such a challenging body, where it is estimated to have a 40 km thick ice shell. An approach for reaching the ocean has been conceived using a melting probe Cryobot concept that has been studied for a potential future mission. A lander is assumed to be the platform from which the Cryobot would be deployed. The ice penetrating vehicle concept consists of a cylindrical, narrow-body probe that encases a radioisotope heat/power source that would be used to do the penetration by melting through the icy crust. The baseline design of the probe includes a suite of science instruments to analyze the ice during descent and the liquid ocean underneath. For communication, a set of fiber optic wire as well as wireless RF in the very cold porous top layer is assumed, and then acoustic modules would be used for the communication in warmer denser ice over distance of 25 km between the modules. In addition to the acoustic communication modules, a sonar is part of the concept, for obstacle avoidance. The focus of this paper is on the use of elastic waves in the 1kHz range.

Lih, Shyh-Shiuh↗

Lowest-energy broad 𝛼-cluster resonances in 19 F

There is strong astrophysical interest in the structure of 19 F states near the α-decay threshold, as their properties are important for understanding the development of α clustering in the 20 Ne region. The emergence of clustered states, and more generally of states that couple strongly to reaction channels near their decay thresholds, is also a topic of current interest in theoretical nuclear physics. In this work, we determine the parameters of broad low-spin states in 19 F near the α-decay threshold and present a theoretical study of these states. The analysis is restricted to ℓ = 0 and 1 resonances in the α + 15 N system close to the α threshold in 19 F. Excitation functions for 15 N(α,α) elastic scattering were measured using the Thick Target Inverse Kinematics (TTIK) method, and these new data, together with older high–energy-resolution measurements, were analyzed within the R-matrix framework. The nuclear structure of 19 F was calculated using configuration-interaction methods with a recently developed effective interaction Hamiltonian. As a result, we identify a sequence of α-cluster resonances in 19 F and map the distribution of clustering strength, which is relevant for astrophysical reaction modeling and for the theoretical understanding of many-body dynamics and the emergence of clustering in loosely bound or unstable nuclei. Furthermore, the work advances theoretical insight into the origins of clustering and highlights open questions for future theoretical and experimental studies.

19F↗

Nuclear properties with semilocal momentum-space regularized chiral interactions beyond N 2 LO

We present a comprehensive investigation of few-nucleon systems as well as light and medium-mass nuclei up to A = 48 using the current Low Energy Nuclear Physics International Collaboration two-nucleon interactions in combination with the third-order (N 2 LO) three-nucleon forces. To address the systematic overbinding of nuclei starting from A ≈ 10 found in our earlier study utilizing the N 2 LO two- and three-nucleon forces, we take into account higher-order corrections to the two-nucleon potentials up through fifth order in chiral effective field theory. Here, the resulting Hamiltonian can be completely determined using the A = 3 binding energies and selected nucleon-deuteron cross sections as input. It is then shown to predict other nucleon-deuteron scattering observables and spectra of light p-shell nuclei, for which a detailed correlated truncation error analysis is performed, in agreement with experimental data. Moreover, the predicted ground state energies of nuclei in the oxygen isotopic chain from 14 O to 26 O as well as 40 Ca and 48 Ca show a remarkably good agreement with experimental values, given that the Hamiltonian is fixed completely from the A ≤ 3 data, once the fourth-order (N 3 LO) corrections to the two-nucleon interactions are taken into account. On the other hand, the charge radii are found to be underpredicted by ≈10% for the oxygen isotopes and by almost 20% for 40 Ca and 48 Ca.

20 ≤ A ≤ 38↗

Residual Strength Analyses of Monolithic Structures

Finite-element fracture simulation methodology predicts the residual strength of damaged aircraft structures. The methodology uses the critical crack-tip-opening-angle (CTOA) fracture criterion to characterize the fracture behavior of the material. The CTOA fracture criterion assumes that stable crack growth occurs when the crack-tip angle reaches a constant critical value. The use of the CTOA criterion requires an elastic- plastic, finite-element analysis. The critical CTOA value is determined by simulating fracture behavior in laboratory specimens, such as a compact specimen, to obtain the angle that best fits the observed test behavior. The critical CTOA value appears to be independent of loading, crack length, and in-plane dimensions. However, it is a function of material thickness and local crack-front constraint. Modeling the local constraint requires either a three-dimensional analysis or a two-dimensional analysis with an approximation to account for the constraint effects. In recent times as the aircraft industry is leaning towards monolithic structures with the intention of reducing part count and manufacturing cost, there has been a consistent effort at NASA Langley to extend critical CTOA based numerical methodology in the analysis of integrally-stiffened panels.In this regard, a series of fracture tests were conducted on both flat and curved aluminum alloy integrally-stiffened panels. These flat panels were subjected to uniaxial tension and during the test, applied load-crack extension, out-of-plane displacements and local deformations around the crack tip region were measured. Compact and middle-crack tension specimens were tested to determine the critical angle (wc) using three-dimensional code (ZIP3D) and the plane-strain core height (hJ using two-dimensional code (STAGS). These values were then used in the STAGS analysis to predict the fracture behavior of the integrally-stiffened panels. The analyses modeled stable tearing, buckling, and crack branching at the integral stiffener using different values of critical CTOA for different material thicknesses and orientation. Comparisons were made between measured and predicted load-crack extension, out-of-plane displacements and local deformations around the crack tip region. Simultaneously, three-dimensional capabilities to model crack branching and to monitor stable crack growth of multiple cracks in a large thick integrally-stiffened flat panels were implemented in three-dimensional finite element code (ZIP3D) and tested by analyzing the integrally-stiffened panels tested at Alcoa. The residual strength of the panels predicted from STAGS and ZP3D code compared very well with experimental data. In recent times, STAGS software has been updated with new features and now one can have combinations of solid and shell elements in the residual strength analysis of integrally-stiffened panels.

Forth, Scott↗