SEARCH · Engineering Papers

Results for “Finite element contact”

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.

SourceResearch connectionFrom yearThrough year

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

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Finite element analysis of steady and transiently moving/rolling nonlinear viscoelastic structure. Part 1: Theory

In a three part series of papers, a generalized finite element analysis scheme is developed to handle the steady and transient response of moving/rolling nonlinear viscoelastic structure. This paper considers the development of the moving/rolling element strategy, including the effects of large deformation kinematics and viscoelasticity modelled by fractional integro-differential operators. To improve the solution strategy, a special hierarchical constraint procedure is developed for the case of steady rolling/translating as well as a transient scheme involving the use of a Grunwaldian representation of the fractional operator. In the second and third parts of the paper, 3-D extensions are developed along with transient contact strategies enabling the handling of impacts with obstructions. Overall, the various developments are benchmarked via comprehensive 2- and 3-D simulations. These are correlated with experimental data to define modelling capabilities.

Padovan, Joe↗

NASA NTRS2009

Thermal-Structural Analysis of PICA Tiles for Solar Tower Test

Thermal protection materials used in spacecraft heatshields are subjected to severe thermal and mechanical loading environments during re-entry into earth atmosphere. In order to investigate the reliability of PICA tiles in the presence of high thermal gradients as well as mechanical loads, the authors designed and conducted solar-tower tests. This paper presents the design and analysis work for this tests series. Coupled non-linear thermal-mechanical finite element analyses was conducted to estimate in-depth temperature distribution and stress contours for various cases. The first set of analyses performed on isolated PICA tile showed that stresses generated during the tests were below the PICA allowable limit and should not lead to any catastrophic failure during the test. The tests results were consistent with analytical predictions. The temperature distribution and magnitude of the measured strains were also consistent with predicted values. The second test series is designed to test the arrayed PICA tiles with various gap-filler materials. A nonlinear contact method is used to model the complex geometry with various tiles. The analyses for these coupons predict the stress contours in PICA and inside gap fillers. Suitable mechanical loads for this architecture will be predicted, which can be applied during the test to exceed the allowable limits and demonstrate failure modes. Thermocouple and strain-gauge data obtained from the solar tower tests will be used for subsequent analyses and validation of FEM models.

Agrawal, Parul↗

DOE OSTI2026

Non-conformal interface-cohesive modeling with the shifted boundary method

The accurate simulation of boundary- and interface-dominated problems on complex geometries remains challenging when boundary- or interface-fitted meshes are difficult to generate, particularly for curved boundaries, polycrystalline microstructures, and dense interface networks. The Shifted Boundary Method (SBM) alleviates this meshing burden by shifting the enforcement of boundary conditions from the true boundary to a nearby surrogate boundary and recovering the effect of the true boundary through geometric correction terms, thereby enabling standard finite element spaces on non-boundary-fitted meshes. In this report, we develop a general shiftedboundary and shifted-interface framework within the open-source MOOSE framework. We first present a general SBM implementation for complex geometries on non-boundary-fitted meshes. We then adopt the Shifted Interface Method (SIM) for internal interfaces and develop a unified shifted-interface treatment in which the interface law is enforced on a surrogate interface and the effect of the true interface is recovered through shifted jumps, fluxes, and tractions. This perspective brings scalar thermal-contact and vector-valued cohesive-zone mechanics into a single framework, the latter realized as the Shifted Cohesive Zone Method (SCZM) and coupled with history-dependent constitutive models from NEML2. We further extend the MOOSE mesh infrastructure to support cohesive-zone calculations on distributed meshes. The framework is verified and demonstrated through three progressive studies: Poisson’s equation on a smoothed starshaped domain, a manufactured thermal-contact problem on a non-interface-fitted mesh, and a two-dimensional polycrystalline representative volume element combining crystal plasticity with cohesive grain-boundary interfaces. Across these studies, the shifted formulations reproduce boundary- and interface-fitted reference solutions with high fidelity, indicating that the proposed framework provides an accurate and efficient route to boundary- and interface-dominated simulations on arbitrary geometries without requiring fitted meshes.

Yang, Cheng-Hau↗

NASA NTRS1995

Resolution of the buoyancy in the 8-foot high temperature tunnel combustor

Currently, the 8-Foot High Temperature Tunnel (8-Ft. HTT) combustor produces a good profile at only one point (2000 psia and 3650 R with oxygen enrichment). Air is enriched with oxygen (liquid) so that the combustor product gas will contain the volumetric amount of oxygen normally found in air. The oxygen enriched air has a large fraction that is not reacted and flows through the outer periphery of the fuel injector. This ring of cold air in addition to the relatively cold walls of the combustor set up buoyancy forces that produce a segregation of relatively cool gases at the bottom of the combustor exit. The basic problem is to produce a test gas that has uniform properties at all combustor conditions. The combustor temperature may be as high as 3700 R or as low as 2000 R. Combustor pressures can be as high as 3500 psia (no oxygen enrichment) and as low as 600 psia. The segregation is most severe with oxygen enriched air, since its temperature is lower and its density is high. The combustor is lined with nickel 201 and can be operated at about 1600 R maximum. A global mixing process is desired that produces an acceptable profile of temperature, species, and velocity at the exit of the combustor. The ultimate goal is a temperature profile with about 100 R variance and about 2 percent variance in oxygen. The exit total temperature must not be lowered significantly by the mixing apparatus or mechanisms employed. If immersed bodies are used, they must also be kept very hot. All combustor wall modifications must be able to survive the heat and structural conditions of the varied operating conditions. Our approach to resolving this issue is being conducted in three stages: (1) Consider mixing exclusively, (2) Resolve the heat transfer concerns resulting from the chosen mixing strategy, and (3) Solve the material and structural problems resulting from stages (1) and (2). Since the 8-Ft. HTT is unavailable for experimentation, the study is conducted exclusively with computational fluid dynamic (CPD) codes (Fluent/Uns and Rampant 3.1) using unstructured grid through body fitted coordinates. Both CFD codes are general purpose Navier-Stokes solution packages that can solve integral conservation equations for conservation of mass, momentum and energy. The governing equations are discretized using a control-volume finite-element method on unstructured triangular 2-D grids. In the interest of time, a 3-D tetrahedral grid was used to check the 2-D results on one mixing strategy and the 2-D results were confirmed. Preliminary results indicate that excellent mixing can be achieved with a body placed in the center of the flowing hot fluid with a minor modification to the combustor wall similar to a model positioned in the test section of a wind tunnel. The concept here, is to create longitudinal vortices strong enough to bring both fluids into intimate contact with each other near the wall where the cold fluid resides. However, there may be a trade-off in the choice of mixing strategy, heat transfer and structural requirements. Currently, we are examining the more promising geometries for heat transfer concerns and developing strategies for the material of construction of the center body.

Loney, Norman W.↗

NASA NTRS2005

Structural Benchmark Tests of Composite Combustion Chamber Support Completed

A series of mechanical load tests was completed on several novel design concepts for extremely lightweight combustion chamber support structures at the NASA Glenn Research Center (http://www.nasa.gov/glenn/). The tests included compliance evaluation, preliminary proof loadings, high-strain cyclic testing, and finally residual strength testing of each design (see the photograph on the left). Loads were applied with single rollers (see the photograph on the right) or pressure plates (not shown) located midspan on each side to minimize the influence of contact stresses on corner deformation measurements. Where rollers alone were used, a more severe structural loading was produced than the corresponding equal-force pressure loading: the maximum transverse shear force existed over the entire length of each side, and the corner bending moments were greater than for a distributed (pressure) loading. Failure modes initiating at the corner only provided a qualitative indication of the performance limitations since the stress state was not identical to internal pressure. Configurations were tested at both room and elevated temperatures. Experimental results were used to evaluate analytical prediction tools and finite-element methodologies for future work, and they were essential to provide insight into the deformation at the corners. The tests also were used to assess fabrication and bonding details for the complicated structures. They will be used to further optimize the design of the support structures for weight performance and the efficacy of corner reinforcement.

Krause, David L.↗

NASA NTRS2018

High Strain Rate Response of Adhesively Bonded Fiber-Reinforced Composite Joints – A Computational Study to Guide Experimental Design

Adhesively bonded carbon fiber-reinforced epoxy composite laminates are widely used in aerospace applications. During a high energy impact event, these laminates are often subjected to high strain rate loading. However, the influence of high strain rate loading on the response of these composite joints is not well understood. Computational finite element (FE) modeling and simulations are conducted to guide the design of high strain rate experiments. Two different experimental designs based on split Hopkinson bar were numerically modeled to simulate Mode I and Mode II types loading in the composite. In addition, the computational approach adopted in this study helps in understanding the high strain rate response of adhesively bonded composite joints subjected to nominally Mode I and Mode II loading. The modeling approach consists of a ply-level 3D FE model, a progressive damage constitutive model for the composite material behavior and a cohesive tie-break contact element for interlaminar delamination.

Ravindran, Suraj↗

NASA NTRS2019

Fast Steering Mirror Disturbance Effects on Overall System Optical Performance for the Large Ultraviolet/optical/infrared Surveyor (LUVOIR) Concept Using a Non-Contact Vibration Isolation and Precision Pointing System

As the optical performance requirements of space telescopes get more stringent, the need to analyze all possible error sources early in the mission design becomes critical. One large telescope with tight performance requirements is the Large Ultraviolet / Optical / Infrared Surveyor (LUVOIR) concept. The LUVOIR concept includes a 15-meter-diameter segmented-aperture telescope with a suite of serviceable instruments operating over a range of wavelengths between 100nm to 2.5um. Using an isolation architecture that involves no mechanical contact between the telescope and the host spacecraft structure allows for tighter performance metrics than current space-based telescopes being flown. Because of this separation, the spacecraft disturbances can be greatly reduced and disturbances on the telescope payload contribute more to the optical performance error. A portion of the optical performance error comes from the disturbances generated from the motion of the Fast Steering Mirror (FSM) on the payload. Characterizing the effects of this disturbance gives insight into the specifications on the FSM needed to achieve the tight optical performance requirements of the overall system. Through analysis of the LUVOIR finite element model and linear optical model given a range of input disturbances at the FSM, the optical performance of the telescope and recommendations for FSM specifications can be determined. The LUVOIR observatory control strategy consists of a multi-loop control architecture including the spacecraft Attitude Control System (ACS), Vibration Isolation and Precision Pointing System (VIPPS), and FSM. This paper focuses on the control loop containing the FSM disturbances and their effects on the telescope optical performance.

Sacks, Lia W.↗

NASA NTRS2008

Informatics and High Throughput Screening of Thermophysical Properties

The combination of computer-aided experiments with computational modeling enables a new class of powerful tools for materials research. A non-contact method for measuring density, thermal expansion, and creep of undercooled and high-temperature materials has been developed, using electrostatic levitation and optical diagnostics, including digital video. These experiments were designed to take advantage of the large volume of data (many gigabytes/experiment, terabytes/campaign) to gain additional information about the samples. For example, using sub-pixel interpolation to measure about 1000 vectors per image of the sample's surface allows the density of an axisymmetric sample to be determined to an accuracy of about 200 ppm (0.02%). A similar analysis applied to the surface shape of a rapidly rotating sample is combined with finite element modeling to determine the stress-dependence of creep in the sample in a single test. Details of the methods for both the computer-aided experiments and computational models will be discussed.

Hyers, Robert W.↗

NASA NTRS2005

Effect of Impact Location on the Response of Shuttle Wing Leading Edge Panel 9

The objective of this paper is to compare the results of several simulations performed to determine the worst-case location for a foam impact on the Space Shuttle wing leading edge. The simulations were performed using the commercial non-linear transient dynamic finite element code, LS-DYNA. These simulations represent the first in a series of parametric studies performed to support the selection of the worst-case impact scenario. Panel 9 was selected for this study to enable comparisons with previous simulations performed during the Columbia Accident Investigation. The projectile for this study is a 5.5-in cube of typical external tank foam weighing 0.23 lb. Seven locations spanning the panel surface were impacted with the foam cube. For each of these cases, the foam was traveling at 1000 ft/s directly aft, along the orbiter X-axis. Results compared from the parametric studies included strains, contact forces, and material energies for various simulations. The results show that the worst case impact location was on the top surface, near the apex.

Lyle, Karen H.↗

NASA NTRS2015

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

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

Burkhardt, Wendel M.↗

NASA NTRS2009

Advances in Non-contact Measurement of Creep Properties

Our team has developed a novel approach to measuring creep at extremely high temperatures using electrostatic levitation (ESL). This method has been demonstrated on niobium up to 2300 C, while ESL has melted tungsten (3400 C). High-precision machined spheres of the sample are levitated in the NASA MSFC ESL, a national user facility, and heated with a laser. The laser is aligned off-center so that the absorbed photons transfer their momentum to the sample, causing it to rotate at up to 250,000+ RPM. The rapid rotation loads the sample through centripetal acceleration, causing it to deform. The deformation of the sample is captured on high-speed video, which is analyzed by machine-vision software from the University of Massachusetts. The deformations are compared to finite element models to determine the constitutive constants in the creep relation. Furthermore, the noncontact method exploits stress gradients within the sample to determine the stress exponent in a single test. This method was validated in collaboration with the University of Tennessee for niobium at 1985 C, with agreement within the uncertainty of the conventional measurements. A similar method is being employed on Ultra-High-Temperature ZrB2- SiC composites, which may see application in rocket nozzles and sharp leading edges for hypersonic vehicles.

Hyers, Robert W.↗

DOE OSTI2022

Enhancing high-fidelity nonlinear solver with reduced order model

Abstract We propose the use of reduced order modeling (ROM) to reduce the computational cost and improve the convergence rate of nonlinear solvers of full order models (FOM) for solving partial differential equations. In this study, a novel ROM-assisted approach is developed to improve the computational efficiency of FOM nonlinear solvers by using ROM’s prediction as an initial guess. We hypothesize that the nonlinear solver will take fewer steps to the converged solutions with an initial guess that is closer to the real solutions. To evaluate our approach, four physical problems with varying degrees of nonlinearity in flow and mechanics have been tested: Richards’ equation of water flow in heterogeneous porous media, a contact problem in a hyperelastic material, two-phase flow in layered porous media, and fracture propagation in a homogeneous material. Overall, our approach maintains the FOM’s accuracy while speeding up nonlinear solver by 18–73% (through suitable ROM-assisted FOMs). More importantly, the proximity of ROM’s prediction to the solution space leads to the improved convergence of FOMs that would have otherwise diverged with default initial guesses. We demonstrate that the ROM’s accuracy can impact the computational efficiency with more accurate ROM solutions, resulting in a better cost reduction. We also illustrate that this approach could be used in many FOM discretizations (e.g., finite volume, finite element, or a combination of those). Since our ROMs are data-driven and non-intrusive, the proposed procedure can easily lend itself to any nonlinear physics-based problem.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2006

NASA Tech Briefs, May 2006

T-Shaped Emitter Metal Structures for HBTs Rigorous Estimation of SNR of a PSK Communication Link Advanced Ka-Band Transceiver With Monopulse Tracking EMI Filters for Low-Temperature Applications Lightweight Electronic Camera for Research on Clouds Pilot Weather Advisor System Waveguide Power-Amplifier Module for 80 to 150 GHz Better Back Contacts for Solar Cells on Flexible Substrates Topics covered include:Tunable, Highly Stable Lasers for Coherent Lidar; Improved Photon-Emission-Microscope System; Program Synthesizes UML Sequence Diagrams; Aspect-Oriented Subprogram Synthesizes UML Sequence Diagrams; Updated Computational Model of Cosmic Rays Near Earth; Software for Alignment of Segments of a Telescope Mirror; Simulation of Dropping of Cargo With Parachutes; DAVE-ML Utility Program; Robust Control for the Mercury Laser Altimeter; Thermally Stable Piezoelectric and Pyroelectric Polymers; Combustion Synthesis of Ca3(PO4)2 Net-Shape Surgical Implants; Stochastic Representation of Chaos Using Terminal Attractors; Two High-Temperature Foil Journal Bearings; Using Plates To Represent Fillets in Finite-Element Modeling; Repairing Chipped Silicide Coatings on Refractory Metal Substrates; Simplified Fabrication of Helical Copper Antennas; Graded-Index "Whispering-Gallery" Optical Microresonators; Optical Profilometers Using Adaptive Signal Processing; Manufacture of Sparse-Spectrum Optical Microresonators; Exact Tuning of High-Q Optical Microresonators by Use of UV; Automation for "Direct-to" Clearances in Air-Traffic Control; Improved Traps for Removing Gases From Coolant Liquids; and Lunar Constellation of Frozen Elliptical Inclined Orbits.

Source record↗

DOE OSTI2022

Large-scale parametric modeling of spent nuclear fuel dynamics in the 30 cm package drop scenario

Packages used to transport spent nuclear fuel (SNF) are required by the U.S. Code of Federal Regulations 10 CFR 71.71 to demonstrate satisfactory performance during a drop scenario. While the CFR is meant to ensure safe package function, it does not evaluate survival of the SNF within. The U.S. Department of Energy Spent Fuel and Waste Science and Technology program is working on closing the knowledge gap related to the response of SNF to external mechanical loads, including the hypothetical 30 cm package drop scenario in the CFR. In support of this effort, LS-DYNA finite element simulations were developed by Pacific Northwest National Laboratory (PNNL) to model generic drop scenarios at both the package and fuel assembly level. The models were validated against one-third scale package and full scale fuel assembly drop test data and were exercised to predict fuel cladding strains in a narrow range of model configurations. This work describes a large-scale parametric study conducted by PNNL using the previously developed and validated PWR finite element model, with the addition of a new generic BWR assembly model. The motivation for the parametric study was to characterize the broad range of SNF responses in the 30 cm package drop scenario. This was accomplished by varying the drop orientation, fuel assembly type (17x17 PWR and 10x10 BWR), burnup, cladding temperature, spacer grid buckling load, package mass, impact limiter stiffness, and mechanical gap conditions within the basket. A MATLAB framework was developed to automate LS-DYNA model generation and execution on PNNL institutional computing resources. In total, over 2000 simulations were performed. For each simulation, the SNF response was quantified in terms of permanent grid deformation, fuel rod contact pressure, and strains within the fuel rods, guide tubes, and water rods. The results provide valuable insight into the range of responses that could be reasonably expected from SNF in the hypothetical drop scenario, as well as the sensitivity to each input parameter. The results of this parametric study are a key component of the testing and modeling strategy the Spent Fuel and Waste Science and Technology program is using to close the external loads knowledge gap.

Kadooka, Kevin↗

DOE OSTI2022

Enhanced critical axial tensile strain limit of CORC ® wires: FEM and analytical modeling

Abstract Conductor on Round Core (CORC ® ) cables and wires are composed of spiraled high-temperature superconducting (HTS) rare-earth barium copper oxide (REBCO) tapes, wound in multiple layers, and can carry very high currents in background magnetic fields of more than 20 T. They combine isotropic flexibility and high resilience to electromagnetic and thermal loads. The brittle nature of HTS tapes limits the maximum allowable axial tensile strain in superconducting cables. An intrinsic tensile strain above about 0.45% will introduce cracks in the REBCO layer of straight HTS tapes resulting in irreversible damage. The helical fashion at which the REBCO tapes are wound around the central core allows tapes to experience only a fraction of the total axial tensile strain applied to the CORC ® wire. As a result, the critical strain limit of CORC ® wires can be increased by a factor of more than 10 that of REBCO tapes. Finite element (FE) and analytical models are developed to predict the performance of CORC ® wires under axial tensile strain. A parametric analysis is carried out by varying the winding angle, the Poisson’s ratio of the CORC ® wire core, the core diameter, and the tape width. The results show that a small variation in winding angle can have a significant impact on the cable’s axial tensile strain tolerance. While the radial contraction of the helically wound tapes in a CORC ® wire under axial tensile strain depends on its winding angle, it is mostly driven by the Poisson’s ratio of the central core, affecting the tape strain state and thus its performance. Contact pressure from multiple layers within the CORC ® wire also affects the CORC ® wire performance. The FE model can be used to optimize the cable design for specific application conditions, resulting in an irreversible strain limit of CORC ® cables and wires as high as 7%.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1988

Three-dimensional analysis of a postbuckled embedded delamination

Delamination growth caused by local buckling of a delaminated group of plies was investigated. Delamination growth was assumed to be governed by the strain energy release rates, G(1), G(2) and G(3). The strain energy release rates were calculated using a geometrically nonlinear three-dimensional finite element analysis. The program is described and several checks of the analysis are discussed. Based on a limited parametric study, the following conclusions were reached: (1) the problem is definitely mixed mode (in some cases G(1) is larger than G(2), for other cases the opposite is true); (2) in general, there is a large gradient in the strain energy release rates along the delamination front; (3) the locations of maximum G(1) and G(2) depend on the delamination shape and the applied strain; (4) the mode 3 component was negligible for all cases considered; and (5) the analysis predicted that parts of the delamination would overlap. The results presented did not impose contact constraints to prevent overlapping. Further work is needed to determine the effects of allowing the overlapping.

Whitcomb, John D.↗

NASA NTRS1989

Three-dimensional analysis of a postbuckled embedded delamination

Whitcomb, John D.↗

NASA NTRS1992

Higher order finite element analysis of thick composite laminates

A higher order, sub-parametric, laminated, 3D solid finite element was used for the analysis of very thick laminated composite plates. The geometry of this element is defined by four nodes in the X-Y plane which define a prism of material through the thickness of the laminate. There are twenty-four degrees of freedom at each node; translations at the upper and lower surfaces of the laminate in each of the three coordinate directions, and the derivatives of these translations with respect to each coordinate. This choice of degrees of freedom leads to displacement and strain compatibility at the corners. Stacking sequence effects are accounted for by explicitly integrating the strain energy density through the thickness of the element. The laminated solid element was combined with a gap-contact element to analyze thick laminated composite lugs loaded through flexible pins. The resulting model accounts for pin bending effects that produce non-uniform bearing stresses through the thickness of the lug. A thick composite lug experimental test program was performed, and provided data that was used to validate the analytical model. Two lug geometries and three stacking sequences were tested.

Goering, J.↗