Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex loading”

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,189 records · Page 66

Quantifying the Solar Energy Resource for Puerto Rico

After Hurricane Maria, multiple U.S. Department of Energy laboratories studied the state of the electric grid in Puerto Rico and analyzed grid resilience and grid integration of renewable energy. As part of the work done at the National Renewable Energy Laboratory, researchers created new solar resource data, conducted a technical potential and supply curve analysis, and studied the interannual variability of the solar resource. A new methodology was developed to downscale solar resource data from the National Solar Radiation Data Base (NSRDB) from a 4-km x 4-km spatial and 30-minute temporal resolution to a 2 km x 2 km and 5-minute resolution. This methodology primarily used simple physical principles to develop high-resolution cloud properties which were then used to compute solar radiation. The high-resolution datasets were validated against ground measurements and the error metrics were found to be similar to the original lower resolution dataset. Using 20 years of downscaled data from the NSRDB multi-year capacity factors for photovoltaics (PV) were developed for both single-axis tracking and fixed latitude-tilt configurations. Use of the multi-year data provides the ability to understand variability in capacity factors due to variability in weather over a long period of time. For Puerto Rico the coastal regions were found to have significant higher capacity factors than inland. Using land-use and terrain information a technical potential analysis was conducted for Puerto Rico. This analysis restricted single PV plant development to a maximum of 100 MW nameplate capacity. The nameplate capacity for each municipality were then determined. Based on our assumptions, 56 of the 78 total municipalities of Puerto Rico contain some level of solar capacity. Most of the interior municipalities did not have any capacity because of the geographic exclusions used in this study. The lowest capacity for a PV plant observed in a municipality was 10 MW. The maximum capacity within a county was 2,000 MW. Further a supply curve analysis was conducted by taking the results of the technical potential and quantifying system and transmission costs. The levelized cost of energy (LCOE) was calculated for each theoretical PV plant site, and the levelized cost of transmission was added to the LCOE to produce a total cost estimate for each site. The results of the supply curve analysis allow for a relative comparison of the cost for integrating new PV capacity into the grid. This analysis indicates that cheaper total LCOE sites tend to be larger in capacity. The total capacity in this study was found to be far beyond the maximum peak load for the island. However, this study does not consider the economic and market potential for development. The cumulative capacity presented in this study assumes that the best locations are developed first and ignores the complex decision paths for new power plant development. Therefore, this analysis can only be treated as illustrative. Finally this study investigates the impact of inter-annual variability of resource using a variety of metrices including probability of exceedance and variation in capacity factor and LCOE. This study demonstrates that the capacity factor or LCOE could vary by over 10% year to year. This clearly indicates the risks involved in using any particular year of data and clearly points to the use of multi-year data to reduce some of the risks related to variability in weather.

14 SOLAR ENERGY↗

Voltage-spike analysis for a free-running parallel inverter

Unwanted and sometimes damaging high-amplitude voltage spikes occur during each half cycle in many transistor saturable-core inverters at the moment when the core saturates and the transistors switch. The analysis shows that spikes are an intrinsic characteristic of certain types of inverters even with negligible leakage inductance and purely resistive load. The small but unavoidable after-saturation inductance of the saturable-core transformer plays an essential role in creating these undesired thigh-voltage spikes. State-plane analysis provides insight into the complex interaction between core and transistors, and shows the circuit parameters upon which the magnitude of these spikes depends.

Lee, F. C. Y.↗

Aerospace technology development of three types of solid state remote power controllers for 120VDC with current ratings of five, and thirty amperes, one type having current limiting

The first generation of remote power controllers (RPC) developed included: a 5-ampere design (Type 1), capable of limiting maximum overload current to 15 amperes for .1 sec; and 5-ampere noncurrent (Type 2) and 30-ampere noncurrent (Type 3) limiting designs, both with selectable instant trip levels for high-current overload. Each design provides overcurrent protection through an inverse I squared T trip-out function with an automatic reset option and demonstrates step-applied fault capability with a 4000-ampere surge, fast-risetime (low-inductance) power source. They also meet MIL - STD - 461A specification for electromagnetic interference. The second generation RPCs traded off specification compliance for reduction in cost and complexity for the Type 1 and 2 designs and give comparable or improved performance in most areas. The noncurrent limiting RPC proved to be a more economical and feasible method of overload protection for certain load types.

Baker, D. E.↗

Session on High Speed Civil Transport Design Capability Using MDO and High Performance Computing

Since the inception of CAS in 1992, NASA Langley has been conducting research into applying multidisciplinary optimization (MDO) and high performance computing toward reducing aircraft design cycle time. The focus of this research has been the development of a series of computational frameworks and associated applications that increased in capability, complexity, and performance over time. The culmination of this effort is an automated high-fidelity analysis capability for a high speed civil transport (HSCT) vehicle installed on a network of heterogeneous computers with a computational framework built using Common Object Request Broker Architecture (CORBA) and Java. The main focus of the research in the early years was the development of the Framework for Interdisciplinary Design Optimization (FIDO) and associated HSCT applications. While the FIDO effort was eventually halted, work continued on HSCT applications of ever increasing complexity. The current application, HSCT4.0, employs high fidelity CFD and FEM analysis codes. For each analysis cycle, the vehicle geometry and computational grids are updated using new values for design variables. Processes for aeroelastic trim, loads convergence, displacement transfer, stress and buckling, and performance have been developed. In all, a total of 70 processes are integrated in the analysis framework. Many of the key processes include automatic differentiation capabilities to provide sensitivity information that can be used in optimization. A software engineering process was developed to manage this large project. Defining the interactions among 70 processes turned out to be an enormous, but essential, task. A formal requirements document was prepared that defined data flow among processes and subprocesses. A design document was then developed that translated the requirements into actual software design. A validation program was defined and implemented to ensure that codes integrated into the framework produced the same results as their standalone counterparts. Finally, a Commercial Off the Shelf (COTS) configuration management system was used to organize the software development. A computational environment, CJOPT, based on the Common Object Request Broker Architecture, CORBA, and the Java programming language has been developed as a framework for multidisciplinary analysis and Optimization. The environment exploits the parallelisms inherent in the application and distributes the constituent disciplines on machines best suited to their needs. In CJOpt, a discipline code is "wrapped" as an object. An interface to the object identifies the functionality (services) provided by the discipline, defined in Interface Definition Language (IDL) and implemented using Java. The results of using the HSCT4.0 capability are described. A summary of lessons learned is also presented. The use of some of the processes, codes, and techniques by industry are highlighted. The application of the methodology developed in this research to other aircraft are described. Finally, we show how the experience gained is being applied to entirely new vehicles, such as the Reusable Space Transportation System. Additional information is contained in the original.

Rehder, Joe↗

Electronic Performance Support for Operational Systems: A Case Study of the Link Monitor and Control Operator Assistant

For complex operational systems, help needs to come from the inside out. It is often not realistic to call a help desk for problems that need immediate attention, especially for tasks that put a heavy cognitive load on the system operator. This session addresses the issues associated with providing electronic performance support for operational systems, including situations where the system is already fielded and can only change through evolution rather than revolution. We present a case study based on our experiences in developing the Link Monitor and Control Operator Assistant for NASA's Deep Space Network (DSN). The goals of the Operator Assistant are to improve the operability of the system and increase the efficiency of mission operations.

Hill, Randall W., Jr.↗

SSME Investment in Turbomachinery Inducer Impeller Design Tools and Methodology

Within the rocket engine industry, SSME turbomachines are the de facto standards of success with regard to meeting aggressive performance requirements under challenging operational environments. Over the Shuttle era, SSME has invested heavily in our national inducer impeller design infrastructure. While both low and high pressure turbopump failures/anomaly resolution efforts spurred some of these investments, the SSME program was a major benefactor of key areas of turbomachinery inducer-impeller research outside of flight manifest pressures. Over the past several decades, key turbopump internal environments have been interrogated via highly instrumented hot-fire and cold-flow testing. Likewise, SSME has sponsored the advancement of time accurate and cavitating inducer impeller computation fluid dynamics (CFD) tools. These investments together have led to a better understanding of the complex internal flow fields within aggressive high performing inducers and impellers. New design tools and methodologies have evolved which intend to provide confident blade designs which strike an appropriate balance between performance and self induced load management.

Zoladz, Thomas↗

Spray-On Foam Insulations for Launch Vehicle Cryogenic Tanks

Spray-on foam insulation (SOFI) has been developed for use on the cryogenic tanks of space launch vehicles beginning in the 1960s with the Apollo program. The use of SOFI was further developed for the Space Shuttle program. The External Tank (ET) of the Space Shuttle, consisting of a forward liquid oxygen tank in line with an aft liquid hydrogen tank, requires thermal insulation over its outer surface to prevent ice formation and avoid in-flight damage to the ceramic tile thermal protection system on the adjacent Orbiter. The insulation also provides system control and stability with throughout the lengthy process of cooldown, loading, and replenishing the tank. There are two main types of SOFI used on the ET: acreage (with the rind) and closeout (machined surface). The thermal performance of the seemingly simple SOFI system is a complex of many variables starting with the large temperature difference of from 200 to 260 K through the typical 25-mm thickness. Environmental factors include air temperature and humidity, wind speed, solar exposure, and aging or weathering history. Additional factors include manufacturing details, launch processing operations, and number of cryogenic thermal cycles. The study of the cryogenic thermal performance of SOFI under large temperature differentials is the subject of this article. The amount of moisture taken into the foam during the cold soak phase, termed Cryogenic Moisture Uptake, must also be considered. The heat leakage rates through these foams were measured under representative conditions using laboratory standard liquid nitrogen boiloff apparatus. Test articles included baseline, aged, and weathered specimens. Testing was performed over the entire pressure range from high vacuum to ambient pressure. Values for apparent thermal conductivity and heat flux were calculated and compared with prior data. As the prior data of record was obtained for small temperature differentials on non-weathered foams, analysis of the different methods is provided. Recent advancements and applications of SOFI systems on future launch vehicles and spacecraft are also addressed.

Fesmire, J. E.↗

Material Characterization and Geometric Segmentation of a Composite Structure Using Microfocus X-Ray Computed Tomography Image-Based Finite Element Modeling

This study utilizes microfocus x-ray computed tomography (CT) slice sets to model and characterize the damage locations and sizes in thermal protection system materials that underwent impact testing. ScanIP/FE software is used to visualize and process the slice sets, followed by mesh generation on the segmented volumetric rendering. Then, the local stress fields around several of the damaged regions are calculated for realistic mission profiles that subject the sample to extreme temperature and other severe environmental conditions. The resulting stress fields are used to quantify damage severity and make an assessment as to whether damage that did not penetrate to the base material can still result in catastrophic failure of the structure. It is expected that this study will demonstrate that finite element modeling based on an accurate three-dimensional rendered model from a series of CT slices is an essential tool to quantify the internal macroscopic defects and damage of a complex system made out of thermal protection material. Results obtained showing details of segmented images; three-dimensional volume-rendered models, finite element meshes generated, and the resulting thermomechanical stress state due to impact loading for the material are presented and discussed. Further, this study is conducted to exhibit certain high-caliber capabilities that the nondestructive evaluation (NDE) group at NASA Glenn Research Center can offer to assist in assessing the structural durability of such highly specialized materials so improvements in their performance and capacities to handle harsh operating conditions can be made.

Abdul-Aziz, Ali↗

NASA Phase2 Unlimited Rights Final Report Research and Technology for Aerospace Propulsion Systems (RTAPS) Task Order: ERA Advanced Core Compressor Technology Program

As part of the NASA Integrated Systems Research Program (ISRP), Environmentally Responsible Aviation (ERA) Project, NASA GRC and GE Aviation have partnered to investigate the technology barriers associated with improved fuel economy of large gas turbine engines. Crucial to improving fuel economy is the increase in pressure ratio (PR) of the core high pressure compressor (HPC). A key enabler to increasing HPC pressure ratios are new technology front stage airfoil designs that have higher efficiency, higher pressure ratios and improvements in front stage part-power operability. The purpose of this work is to design and test two high pressure ratio, high efficiency front block designs. The baseline design, hereafter referred to as Build 1, is similar to what would be found in a high efficiency 30:1 PR class HPC. A redesigned front block with improved blading will also be tested and hereafter be referred to as Build 2. The goal of Build 2 is to increase stage aero loading (increase stage pressure ratio) beyond Build 1 while maintaining equivalent efficiency and stall margin. The detailed inter-stage and traverse data from the two builds will be used to improve understanding of the complex flow physics at high-speed and part-speed conditions that impact performance and operability. By addressing, understanding and solving the challenges associated with aerodynamic losses and stage matching better design strategies can be incorporated in future core compressors with higher efficiency, thereby enabling improvement in mission fuel burn while still meeting operability requirements.

Axial Compressors↗

NSUF BRR Cask Shipment Summary Report

On September 16, 2021 the Idaho National Laboratory shipped a Battelle Research Reactor (BRR) cask loaded with five Nuclear Science User Facilities (NSUF) experiments and one Accident Tolerant Fuels (ATF) experiment from the Advanced Test Reactor (ATR) to the Hot Fuel Examination Facility (HFEF) at the INL Materials and Fuels Complex (MFC) for post-irradiation examination.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

TRUST Sensors in Environments: Thermocouples (SE-TC), Release FY25

The Delivery Environments Testbeds to Reduce Uncertainty in Simulations and Tests (TRUST) project is a broad project intended to analyze simplified problems experimentally and with modeling and simulation. The purpose of analyzing these simplified problems is to extend solution methods to more complex problems, as well as understand deficiencies and gaps in knowledge of methods currently used in more complex analyses. The TRUST project encompasses several smaller testbeds intended to isolate individual phenomena. The testbed under consideration in this report is the Sensors in Environments: Thermocouples testbed. In previous years, the purpose of this testbed was to quantify uncertainty of thermocouple sensors. To accomplish this, an aluminum plate was placed in a thermal chamber and subject to various types of thermal loading. Thermocouples were placed in various locations on the aluminum plate in various configurations (e.g., embedded in the plate, placed under Kapton tape), and an effort was made to quantify uncertainty in these measurements. Finite element simulations were performed to investigate how sensitive these measurements were to parameters such as the boundary conditions on the plate and material properties. However, a fundamental source of uncertainty in this analysis was the convective heat transfer from the plate. Convective heat transfer is a complex physical phenomenon comprised of a number of interacting sub-processes and is difficult to predict accurately a priori. As such, the main purpose of this testbed in FY25 was to better understand, both experimentally and numerically, the convective heat transfer from the plate. This is a highly applicable problem to several more complex problems, as convective heat transfer occurs in nearly all problems where a body is moving through air. Numerically, this required a two-step approach. First, the air flow in the thermal chamber was in vestigated using computational fluid dynamics. The commercial solver Fluent was used to perform these simulations. From these simulations, a heat transfer coefficient over the surface of the plate was calculated. This heat transfer was then used as boundary conditions for finite element heat transfer simulations within the plate, which were performed using Abaqus. Significant effort was devoted to automating the handoff between these two solvers. Experimentally, previous thermocouple results in the plate were used to validate the time-dependent thermal profiles produced from Abaqus. Further experimental efforts were performed both to help validate the Fluent simulations and to inform its boundary conditions. For example, hot-wire anemometers were used to measure the velocity in the chamber, which would be particularly useful in understanding the chamber inlet velocity. Thermocouple measurements were also taken in the chamber, instead of only on the plate, to serve as validation evidence for the Fluent simulations. Numerical results showed that the Fluent to Abaqus workflow matched previous plate thermocouple measurements well. This type of handoff is useful for more complex experiments, or those that are not able to be examined in as great of detail as this testbed, as it was performed without any experimental input. Experimental results, however, were more mixed. The anemometers proved unreliable, with inconsistent measurements across all anemometers, even at locations that were nearly identical. On the other hand, the thermocouples provided a relatively rich view of the temperature field in the chamber.

42 ENGINEERING↗

Monolithic Flexure Pre-Stressed Ultrasonic Horns

High-power ultrasonic actuators are generally assembled with a horn, backing, stress bolt, piezoelectric rings, and electrodes. The manufacturing process is complex, expensive, difficult, and time-consuming. The internal stress bolt needs to be insulated and presents a potential internal discharge point, which can decrease actuator life. Also, the introduction of a center hole for the bolt causes many failures, reducing the throughput of the manufactured actuators. A new design has been developed for producing ultrasonic horn actuators. This design consists of using flexures rather than stress bolts, allowing one to apply pre-load to the piezoelectric material. It also allows one to manufacture them from a single material/plate, rapid prototype them, or make an array in a plate or 3D structure. The actuator is easily assembled, and application of pre-stress greater than 25 MPa was demonstrated. The horn consists of external flexures that eliminate the need for the conventional stress bolt internal to the piezoelectric, and reduces the related complexity. The stress bolts are required in existing horns to provide prestress on piezoelectric stacks when driven at high power levels. In addition, the manufacturing process benefits from the amenability to produce horn structures with internal cavities. The removal of the pre-stress bolt removes a potential internal electric discharge point in the actuator. In addition, it significantly reduces the chances of mechanical failure in the piezoelectric stacks that result from the hole surface in conventional piezoelectric actuators. The novel features of this disclosure are: 1. A design that can be manufactured from a single piece of metal using EDM, precision machining, or rapid prototyping. 2. Increased electromechanical coupling of the horn actuator. 3. Higher energy density. 4. A monolithic structure of a horn that consists of an external flexure or flexures that can be used to pre-stress a solid piezoelectric structure rather than a bolt, which requires a through hole in the piezoelectric material. 5. A flexure system with low stiffness that accommodates mechanical creep with minor reduction in pre-stress.

Sherrit, Stewart↗

Controlling the size of non-axisymmetric magnetic footprints using resonant magnetic perturbations

The structure of the non-axisymmetric heat load distribution at the divertor plates is determined not only by the toroidal but also from the poloidal spectrum of non-axisymmetric field perturbations. Whether they are intrinsic, like error fields, or they are applied through 3D coils, the non-axisymmetric fields produce complex 3D edge magnetic topologies (footprints) that alter the properties of the heat and particle flux distributions on the divertor target plates. In this manuscript, a study of the impact of applied 3D field poloidal spectrum on the footprint size and structure is done for the DIII-D tokamak using the resistive MHD code M3D-C1 coupled with the field line tracing code TRIP3D. To resolve the impact of the poloidal spectrum of the magnetic perturbation, the relative phase of the two rows of in-vessel 3D coils used to produce both a n = 2 and a n = 3 perturbation is varied, where n is the toroidal harmonic of the magnetic perturbation. This shows that the largest footprint is predicted when the relative phase of the two rows is close to zero, which is also where the resonant coupling with the plasma is maximized. These results suggest that it will be challenging to decouple the footprint size from the requisite resonant coupling for RMP–ELM control. The correlation between the measured heat load and particle flux distributions at the outer divertor plates in DIII-D and the magnetic measurements is in good agreement with the predicted dependence of the magnetic footprint size on the amplitude of the resonant component of the plasma response.

3D fields↗

Evolution of global and local deformation in additively manufactured octet truss lattice structures

Additively manufactured lattice truss structures, often referred to as architected cellular materials, present significant advantages over conventional structures due to their unique characteristics such as high strength-to-weight ratios and surface area-to-volume ratios. These geometrically complex structures, however, come with concomitant challenges for qualification and inspection. In this study, compression testing interrupted with micro-computed tomography inspection was conducted to monitor the evolution of global and local deformation throughout the loading process of 304L stainless steel octet truss lattice structures. Both two- and three-dimensional image analysis techniques were leveraged to characterize geometric heterogeneities resulting from the laser powder bed fusion manufacturing process as well as track the structure throughout deformation. Variations from model-predicted behavior resulting from these heterogeneities are considered relative to the predicted and actual responses of the structures during compression to better understand, model, and predict the octet truss lattice structure compression response.

36 MATERIALS SCIENCE↗

Plate motion in sheared granular fault system

Plate motion near the fault gouge layer, and the elastic interplay between the gouge layer and the plate under stick-slip conditions, is key to understanding the dynamics of sheared granular fault systems. Here, a two-dimensional implementation of the combined finite-discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly simulate a sheared granular gouge fault system. In this work we focus on investigating the influence of normal load, driving shear velocity and plate stiffness on the velocities and displacements in the direction parallel to the shear direction (x-direction) measured at locations on the upper and lower plates just adjacent to the gouge. The simulations show that during slip phases the magnitudes of the measured velocities on the upper and lower plates are proportional to the normal load and may be inversely proportional to the square root of the plate's shear modulus. Whereas, the driving shear velocity does not show distinct influence on the measured velocities. Additionally, large slip velocities are generally associated with large macroscopic friction coefficient drops. For the models subjected to smaller normal loads, larger shear velocities and with stiffer shear plates, the same magnitude of slip velocity could cause a larger drop of macroscopic friction coefficient. During stick phases, the velocities of the upper and lower plates are respectively slightly greater and slightly smaller than half of the driving shear velocity and are both in the same direction of shear. The shear strain rate of the gouge is calculated from this velocity difference between the upper and lower plate during stick phases and thus the gouge effective shear modulus can be calculated. The results show that the gouge effective shear modulus increases proportionally with normal load, while the influence of shear velocity and plate stiffness on gouge effective shear modulus is minor. The simulations address the dynamics of a laboratory-scale fault gouge system and may aid in revealing the complexities of earthquake frictional dynamics.

58 GEOSCIENCES↗

PFSA-Ionomer Adsorption to C and Pt/C Particles in Fuel-Cell Inks

Catalyst inks used to make fuel-cell electrodes consist of Pt/C catalyst particles and a perfluorosulfonic acid (PFSA) ionomer dispersed in water/alcohol solvent mixtures. PFSA ionomer in the ink adsorbs to the surface of the catalyst particles, dictating the dispersion colloid properties. Following adsorption, the subsequent distribution of excess nonadsorbed ionomer in the ink then governs the final structure of the electrode. Here, we characterize the adsorption of the PFSA ionomer onto Pt/C catalyst particles. PFSA adsorption is largely irreversible. Adsorbed sulfonic-acid moieties impart a negative charge on the catalyst surface, causing electrostatic repulsion between the free ionomer in solution and the ionomer-covered Pt/C particle surface. The amount of adsorption is limited by the resulting electrostatic charge that grows as more ionomer adsorbs, and the catalyst surface becomes more negatively charged. Attenuating electrostatic repulsion by increasing the ink ionic strength promotes ionomer adsorption. Electrostatically limited adsorption is observed, irrespective of the solvent water/n-propanol ratio or the catalyst particle porosity and Pt loading. Experimentally measured ionomer adsorption isotherms are well predicted by a Smoluchowski-based kinetic adsorption model, in which the electrostatic energy barrier for adsorption is predicted from DLVO theory. These findings help to unravel the complex phenomena within these colloidal dispersions, allowing for subsequent tailoring of inks to optimize fuel-cell electrode structure and performance.

Adsorption↗

Material Model Parameters Optimization in Liquid Mercury Target Dynamics Simulation With Machine Learning Surrogates

A pulsed spallation target is subjected to very short (∼0.7μs) but intense loads (23.3 kJ) from repeated proton pulses, which knock away neutrons from the mercury atoms’ nuclei for a wide range application in physics, engineering, medicine, petroleum exploration, biology, chemistry, etc. The effect of this pulsed loading on the stainless-steel target module which contains the flowing mercury target material is difficult to predict not only due to its short but intense explosive-like physical reaction, but also the nonlinear material behavior of the liquid mercury in the structure. Injecting small helium bubbles in the mercury has been an efficient method of mitigating the pressure wave at high power level stage. However, prediction of the resultant loading on the target is more difficult when helium gas is intentionally injected into the mercury. A 2-phase material model that incorporates the Rayleigh-Plesset (R-P) model is expected to address this complex multi-physics dynamics problem by including the bubble dynamics in the liquid mercury. A parameter sensitivity study was firstly employed to understand their impact on the simulation strains. The investigated parameters included E, μ, γ, σ, n, VFgas, and gas cumulative volume curve control parameters a and b. Verification and validation results from sparse polynomial expansions (SPE) method and directional Gaussian smoothing (DGS) optimization show that the surrogate model had training error of ∼7% and validation error of ∼15%, indicating that machine learning methods and surrogate models can help optimize the uncertain parameters in the complex 2-phase material model. This approach is expected to fill the knowledge gap between unknown liquid-gas mixture material model and measured vessel strain responses.

Lin, Lianshan↗

Accurate interlaminar stress recovery from finite element analysis

The accuracy and robustness of a two-dimensional smoothing methodology is examined for the problem of recovering accurate interlaminar shear stress distributions in laminated composite and sandwich plates. The smoothing methodology is based on a variational formulation which combines discrete least-squares and penalty-constraint functionals in a single variational form. The smoothing analysis utilizes optimal strains computed at discrete locations in a finite element analysis. These discrete strain data are smoothed with a smoothing element discretization, producing superior accuracy strains and their first gradients. The approach enables the resulting smooth strain field to be practically C1-continuous throughout the domain of smoothing, exhibiting superconvergent properties of the smoothed quantity. The continuous strain gradients are also obtained directly from the solution. The recovered strain gradients are subsequently employed in the integration o equilibrium equations to obtain accurate interlaminar shear stresses. The problem is a simply-supported rectangular plate under a doubly sinusoidal load. The problem has an exact analytic solution which serves as a measure of goodness of the recovered interlaminar shear stresses. The method has the versatility of being applicable to the analysis of rather general and complex structures built of distinct components and materials, such as found in aircraft design. For these types of structures, the smoothing is achieved with 'patches', each patch covering the domain in which the smoothed quantity is physically continuous.

Tessler, Alexander↗