Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient stiffness”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A Computational Framework to design 3D stiffness gradient acoustic metamaterials for impedance matching

Acoustic waves play a crucial role in various applications, including medical imaging, non-destructive testing, and sonar systems. One of the significant challenges in these applications is impedance matching, which is essential for minimizing reflections and maximizing the transfer of acoustic energy between different media. Acoustic metamaterials offer a promising solution to this challenge. In addition to impedance control, gradient stiffness can enhance structural efficiency and enable spatial control of wave propagation, making it a valuable feature in acoustic metamaterial design. In this pa- per, we present our developed computational method to design 3D stiffness gradient acoustic metamaterials for impedance matching. The key steps in our approach include generating initial designs using a periodic covariance function to provide unit cells that are both periodic on the boundaries and randomly formed inside the unit cell. Furthermore, we integrated manufacturing constraints into the design process, ensuring that the structures are interconnected for fabrication. We propose two computational optimization algorithms: GenUnit, based on a non-dominated sorting genetic algorithm (NSGA-II), and MLMatch, which leverages differentiable machine learning. The two approaches are not separate contributions but complementary com- ponents of a unified framework. GenUnit requires no training data and directly interfaces with physics-based simulations, making it highly accurate but slower for large-scale exploration. In contrast, MLMatch is data-hungry during training but, once trained, enables near-instantaneous inference and broad design-space coverage. Together, they form a hybrid strategy: ML- Match rapidly explores the global design space, and GenUnit provides local refinement with high-fidelity accuracy. This balance between training cost, inference time, and precision is the motivation for including both methods in the same study. We applied this dual-algorithm framework to generate two metallic-based metamaterial designs that match the acoustic impedance of water while exhibiting a controlled gradient in stiffness (from stiff to soft). The stiffness gradient is particularly advantageous in applications where one side of the structure must interface with soft or sensitive surfaces, such as human tissue or delicate components. Here, this work paves the way for improved materials in various acoustic applications, particularly in ultrasound devices, by providing better impedance.

Metamaterial↗
NASA NTRS1997

Response of Composite Panels with Stiffness Gradients Due to Stiffener Terminations and Cutouts

The results of an analytical and experimental study of stiffened graphite-epoxy compression panels with terminated stiffeners are presented. The local stress gradients at the stiffener termination location are determined by finite element analysis. Three stiffener termination concepts are evaluated by analysis to determine the stiffener and skin laminate parameters that affect the panel response and failure. The effects of changing local skin laminate definitions, skin reinforcement details, and stiffener termination details on local stress gradients and load-path eccentricities are discussed. Analytical and test results are presented for panels with one terminated stiffener and for panels with one terminated stiffener and two unterminated stiffeners. The effects of a cutout in the skin of a panel with a terminated stiffener is also evaluated to determine the interaction between the stress gradients in the panel due to the cutout and those due to the terminated stiffener. The results of the study indicate that the critical failure modes of the panels initiate at the skin-stiffener interface near the end of the terminated stiffener.

Ambur, Damodar R.↗
DOE OSTI2021

Examination of stiff ion temperature gradient mode physics in simulations of DIII-D H-mode transport

This paper presents comparisons of predictions from two models of turbulent transport in tokamaks. The focus is on assessing how close a reduced but significantly less computationally demanding model (TGLF) does in matching the more detailed model (CGYRO) in predicting scalings with driving gradient ("stiffness"). The impacts of other physics mechanisms such as flow shear, fast ions, and electromagnetic fluctuations on these predictions is also studied. Overall, while there are clear differences in the scaling, these do not appear to significantly impact the fidelity of the reduced model in prediction the experimental conditions studied relative to the more advanced model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Mechanically graded granular scaffolds for osteochondral tissue engineering

Engineered scaffolds designed to approximate the mechanical microenvironment of the osteochondral unit often address this complexity using discrete, two-phase architectures that introduce mechanical discontinuities and interfacial stress concentrations rather than a contiguous stiffness transition. To address this challenge, we created a photoannealed polyethylene glycol (PEG) granular scaffold with a spatially controlled stiffness gradient within a cell-permissive, macroporous architecture. Stiffness was dictated by photoannealing microgels using a photomask. We tuned void volume and available surface area by varying microgel diameter and tested how mesenchymal stromal cells (MSCs) interpret local mechanical environments. MSCs exhibited position-dependent differences in morphology, cytoskeletal structure, matrix deposition, and lineage-specific gene expression within the gradient scaffolds. Softer regions supported rounded cell morphology and deposition of a glycosaminoglycan-rich matrix, whereas stiffer regions promoted cell elongation, increased cytoskeletal tension, and expression of mineral-associated markers. Gradients formed from smaller microgels magnified these spatial responses by increasing cellular confinement and adhesion site availability. Disruption of actomyosin contractility eliminated these regional differences, demonstrating that MSCs rely on tension-dependent mechanotransduction to interpret the gradient. These findings reveal that coupling microgel architecture with continuous stiffness transitions provides a tractable platform to study multiscale mechanobiologic regulation and spatially guide osteochondral tissue formation.

Biological and medical sciences↗
DOE OSTI2022

Examination of stiff ion temperature gradient mode physics in simulations of DIII-D H-mode transport

A systematic evaluation of gyrokinetic and gyrofluid model predictions of ion temperature gradient (ITG) stability and transport using parameters from DIII-D high confinement mode (H-mode) plasmas has been performed. The nonlinear CGYRO code is used to make the gyrokinetic predictions, and the quasilinear TGLF model for the corresponding gyrofluid predictions. The assessments are made at three radii (normalized toroidal flux ρtor = 0.4, 0.55, and 0.7) in three different plasma scenarios with varying levels of neutral beam heating and torque. For each of the nine cases (3 radii × 3 scenarios) considered, ITG turbulence is found to be the dominant long-wavelength instability and transport mechanism. The inclusions of both transverse magnetic fluctuations and dynamic fast beam ions are stabilizing for all cases considered, with strongest effects seen at ρor = 0.4 where the fast ion population and normalized plasma pressure β = 2μ0nT/B2 are highest. The further inclusion of parallel magnetic fluctuations does not have a meaningful impact on the ITG turbulence in these scenarios, but does destabilize (in combination with fast ions) new high-frequency instabilities at ρtor = 0.4 in the high power scenarios. In each case the linear and nonlinear ITG critical gradients are predicted to be lower than the measured ITG scale lengths and their associated uncertainties. Inclusion of equilibrium flow shear in the transport predictions generally leads to an upshift in effective critical gradient rather than a qualitative change in the predicted stiffness, with stronger responses typically seen in the gyrokinetic predictions than in the gyrofluid results. However, in most cases these upshifted gradients still remain below the measured values and their uncertainties. Although the predicted critical gradients are below the measured gradients, both models predicted flux-matching gradients consistent with measured values in six of the nine cases considered, with no clear systematic over- or underprediction. Thus, while the experimental ion temperature profiles do not appear to be closely pinned to the ITG critical gradient, both gyrokinetic and gyrofluid models are able to accurately match the measured gradients reasonably well in most cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Investigating profile stiffness and critical gradients in shaped TCV discharges using local gyrokinetic simulations of turbulent transport

The experimental observation made on the TCV tokamak of a significant confinement improvement in plasmas with negative triangularity ( δ < 0) compared to those with standard positive triangularity has been interpreted in terms of different degrees of profile stiffness (Sauter et al. Phys. Plasmas 21 055906, 2014) and/or different critical gradients. Employing the Eulerian gyrokinetic code GENE (Jenko et al. Phys. Plasmas 7 1904, 2000), profile stiffness and critical gradients are studied under TCV relevant conditions. For the considered experimental discharges, trapped electron modes (TEMs) and electron temperature gradient (ETG) modes are the dominant microinstabilities, with the latter providing a significant contribution to the non-linear electron heat fluxes near the plasma edge. Two series of simulations with different levels of realism are performed, addressing the question of profile stiffness at various radial locations. Retaining finite collisionality, impurities and electromagnetic effects, as well as the physical electron-to-ion mass ratio are all necessary in order to approach the experimental flux measurements. However, flux-tube simulations are unable to fully reproduce the TCV results, pointing towards the need to carry out radially nonlocal (global) simulations, i.e. retaining finite machine size effects, in a future study. Some conclusions about the effect of triangularity can nevertheless be drawn based on the flux-tube results. In particular, the importance of considering the sensitivity to both temperature and density gradient is shown. The flux tube results show an increase of the critical gradients towards the edge, further enhanced when δ < 0, and they also appear to indicate a reduction of profile stiffness towards plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Role of fast ion pressure in the isotope effect in JET L-mode plasmas

This paper presents results of JET ITER-like wall L-mode experiments in hydrogen and deuterium (D) plasmas, dedicated to the study of the isotope dependence of ion heat transport by determination of the ion critical gradient and stiffness by varying the ion cyclotron resonance heating power deposition. When no strong role of fast ions in the plasma core is expected, the main difference between the two isotope plasmas is determined by the plasma edge and the core behavior is consistent with a gyro-Bohm scaling. When the heating power (and the fast ion pressure) is increased, in addition to the difference in the edge region, also the plasma core shows substantial changes. The stabilization of ion heat transport by fast ions, clearly visible in D plasmas, appears to be weaker in H plasmas, resulting in a higher ion heat flux in H with apparent anti-gyro-Bohm mass scaling. The difference is found to be caused by the different fast ion pressure between H and D plasmas, related to the heating power settings and to the different fast ion slowing down time, and is completely accounted for in non-linear gyrokinetic simulations. The application of the TGLF quasi-linear model to this set of data is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Impact of ionization and transport on pedestal density structure in DIII-D and Alcator C-Mod

Abstract This paper investigates the role of ionization on the pedestal structure using both measurements and modeling for H-mode plasmas on DIII-D and Alcator C-Mod to enhance our ability to predict pedestal behavior in future pilot plants. The impact of the neutral penetration depth on the pedestal density is investigated using dimensionally matching hydrogen and deuterium DIII-D H-mode discharges at low and high electron density. The DIII-D Lyman- α diagnostic measurements show that hydrogen neutrals penetrate deeper inside the plasma on both the high field and low field side, while the pedestal electron density structure is similar for both isotopes. However, as the opaqueness increases we observe that the pedestal density gradient becomes stiff, similar to prior observations on DIII-D and C-Mod (Mordijck 2020 Nuclear Fusion 60 082006). In addition, these results also confirm prior measured and modeled poloidal asymmetries in neutral densities, indicating that to make transport predictions, 2D neutral modeling is necessary. The first direct validation of SOLPS-ITER for the measured brightness, emissivity and neutral densities for three different confinement regimes on C-Mod is introduced. The SOLPS-ITER model shows good agreement, within the constrains of the model for all regimes. In addition, a comparison of SOLPS-ITER modeling for DIII-D and C-Mod shows that as opaqueness increases, the role of divertor fueling and thus poloidal asymmetries in the neutral density profiles decreases. Based on these experimental and modeling results we estimate the size of a potential particle pinch using typical values for the diffusion coefficient for both DIII-D and C-Mod H-mode discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Numerical analysis of bipropellant combustion in liquid thrust chambers by an Eulerian-Eulerian approach

The liquid thrust chambers performance (LTCP) code is used for parametric studies of flow and combustion in liquid rocket engines. Multiphase flow equations are solved in an Eulerian-Eulerian framework, and multistep finite rate chemistry is incorporated. The discretization scheme is fully implicit and is based on the total variation diminishing (TVD) scheme, which is accurate, robust, very efficient and capable of handling steep gradients and stiff chemistry. Effects of injection velocity and chamber size have been considered, and the effect of group combustion on the evaporation rate has been studied for a dense spray.

Dang, A. L.↗
NASA NTRS1971

A study of the flexible body dynamics of the RAE-A satellite

Spacecraft flight data was compared with simulated data computed by the RAE and IMP dynamics simulator computer programs. There was not sufficient data available to provide definitive data on antenna boom position to the Radio Astronomy experimenter. There was, however, enough data to indicate the trend of antenna boom stiffness and temperature gradient. It was possible to obtain a good fit of simulated data to measured spacecraft attitude data when the antenna booms were extended to 450 feet. The value of the RAE-A antenna boom stiffness appears to be significantly less than the preflight ground measured value. It was confirmed that an improved simulation can be obtained by using values of the root angles of the antenna booms and by using an angle of the damper plane of motion, relative to the plane of the antennas, which differed from preflight ground measurements.

Source record↗
DOE OSTI2019

MR elastography of liver at 3 Tesla: comparison of gradient-recalled echo (GRE) and spin-echo (SE) echo-planar imaging (EPI) sequences and agreement across stiffness measurements

Purpose: To compare 2D gradient-recalled echo (GRE) and 2D spin-echo (SE) echo-planar imaging (EPI) MR elastography (MRE) for measurement of hepatic stiffness in adult patients with known or suspected liver disease at 3 Tesla. Materials and methods: Three hundred and eighty-seven consecutive patients underwent MRE of the liver at 3 Tesla with 2D-GRE and 2D-SE-EPI sequences. ‘Mean liver stiffness (LS)’ calculated by averaging 3 ROIs in the right lobe, ‘Maximum LS’ calculated by an ROI in the right lobe; and ‘Freehand LS’ calculated by an ROI in the entire liver were measured by two independent readers. Inter-observer and inter-class variability in stiffness measurements were assessed. Stiffness values were correlated with degree of liver fibrosis (METAVIR scores) in 97 patients who underwent biopsy. The diagnostic performance was compared by a receiver-operating characteristic analysis. Results: The technical failure rate was 2.8% for 2D-SE-EPI (11/387) and 4.1% for 2D-GRE (16/387, 9 had R{sub 2}* > 80 s{sup −1} indicating iron overload). There is high reproducibility for both GRE and SE-EPI variants (ICC = 0.84–0.94 for both GRE and SE-EPI MRE). The highest sensitivity, specificity, and accuracy of differentiating mild fibrosis (F0–F2) from advanced fibrosis (F3–F4) are 0.84 (GRE Freehand measurement), 0.92 (GRE Maximum stiffness measurement), and 0.88 (GRE Freehand measurement), respectively. Conclusions: High intra-class correlation and intra-reader correlation are seen on measured hepatic stiffness for both 2D-GRE and 2D-SE-EPI MRE. 2D-SE-EPI has lower failure rate. Diagnostic performance of both sequences is equivalent, with highest sensitivity for 2D-GRE Freehand stiffness measurement, and highest specificity 2D-GRE Maximum stiffness measurement.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2022

A Polymer Canvas with the Stiffness of the Bone Matrix to Study and Control Mesenchymal Stem Cell Response

Reproducing in vitro the complex multiscale physical features of human tissues creates novel biomedical opportunities and fundamental understanding of cell–environment interfaces and interactions. While stiffness has been recognized as a key driver of cell behavior, systematic studies on the role of stiffness have been limited to values in the KPa–MPa range, significantly below the stiffness of bone. Here, in this work, a platform enabling the tuning of the stiffness of a biocompatible polymeric interface up to values characteristic of human bone is reported, which are in the GPa range, by using extremely thin polymer films on glass and cross-linking the films using ultraviolet (UV) light irradiation. It is shown that a higher stiffness is related to better adhesion, proliferation, and osteogenic differentiation, and that it is possible to switch on/off cell attachment and growth by solely tuning the stiffness of the interface, without any surface chemistry or topography modification. Since the stiffness is tuned directly by UV irradiation, this platform is ideal for rapid and simple fabrication of stiffness patterns and gradients, thus representing an innovative tool for combinatorial studies of the synergistic effect of tissue environmental cues on cell behavior, and creates new opportunities for next-generation biosensors, single-cell patterning, and lab-on-a-chip devices.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Machine Learning Predicts the Timing and Shear Stress Evolution of Lab Earthquakes Using Active Seismic Monitoring of Fault Zone Processes

Abstract Machine learning (ML) techniques have become increasingly important in seismology and earthquake science. Lab‐based studies have used acoustic emission data to predict time‐to‐failure and stress state, and in a few cases, the same approach has been used for field data. However, the underlying physical mechanisms that allow lab earthquake prediction and seismic forecasting remain poorly resolved. Here, we address this knowledge gap by coupling active‐source seismic data, which probe asperity‐scale processes, with ML methods. We show that elastic waves passing through the lab fault zone contain information that can predict the full spectrum of labquakes from slow slip instabilities to highly aperiodic events. The ML methods utilize systematic changes in P‐wave amplitude and velocity to accurately predict the timing and shear stress during labquakes. The ML predictions improve in accuracy closer to fault failure, demonstrating that the predictive power of the ultrasonic signals improves as the fault approaches failure. Our results demonstrate that the relationship between the ultrasonic parameters and fault slip rate, and in turn, the systematically evolving real area of contact and asperity stiffness allow the gradient boosting algorithm to “learn” about the state of the fault and its proximity to failure. Broadly, our results demonstrate the utility of physics‐informed ML in forecasting the imminence of fault slip at the laboratory scale, which may have important implications for earthquake mechanics in nature.

58 GEOSCIENCES↗
DOE OSTI2022

Fibrous anisotropy and mineral gradients within the radula stylus of chiton: Controlled stiffness and damage tolerance in a flexible biological composite

Over hundreds of millions of years, organisms have evolved architected structures via precise control over hierarchically assembled components, including the integration of dissimilar materials. One such example is found in the radula system of chitons, intertidal mollusks that feed on algae growing on the rock. Their radula consists of multiple rows of ultrahard teeth, each integrated with a foldable belt-like substrate via a stiff, yet flexible stylus, which is essential for efficient rasping during the feeding process. Here, we investigate the nano and micro-scale components and architectures as well as regional mechanical properties of the stylus, and their subsequent role during the rasping of Cryptochiton stelleri. Three important factors were determined to contribute to the regio-specific stiffness of the stylus: the presence of mineral components, highly oriented chitinous fibers, and a chemically cross-linked protein matrix. All these factors are varied throughout the stylus. There is a high mineral content on the trailing edge close to the tooth and a cross-linked matrix on the leading edge, both with orientational specific oriented chitin fibers that provide force transduction to the tooth. Conversely, there is a significant lack of mineral or cross-linked matrix in the proximal end as well as a low degree of fiber orientation, resulting in a flexible region that can accommodate torsion and flexure during rasping. Understanding the graded composite structure of the stylus and applying this unique design to various engineering fields such as soft robotics, biotechnology, and the medical industry, can inspire the production of high-performance materials.

Lee, Jung-Eun↗
DOE OSTI2019

Magnetic resonance elastography SE-EPI vs GRE sequences at 3T in a pediatric population with liver disease

Purpose: The goal of our study is to compare hepatic stiffness measures using gradient-recalled echo (GRE) versus spin-echo echo planar imaging (SE-EPI)-based MR Elastography (MRE) at 3T used to measure hepatic stiffness in a patients with suspected liver diseases. Materials and methods: This retrospective study included 52 patients with liver disease who underwent a 3T MRE exam including both an investigational SE-EPI-based technique and a product GRE-based technique. Regions of interest (ROI) were placed on the elastograms to measure elastography-derived liver stiffness as well as the area included within the ROIs. The mean liver stiffness values and area of ROIs were compared. Results: The mean liver stiffness was 3.72 kilopascal (kPa) ± 1.29 using GRE MRE and 3.78 kPa ± 1.13 using SE-EPI MRE. Measurement of liver stiffness showed excellent agreement between the two pulse sequences with a mean bias of − 0.1 kPa (range − 1.8 to 1.7 kPa) between sequences. The mean measurable ROI area was higher with SE-EPI (313.8 cm{sup 2} ± 213.8) than with the GRE technique (208.6 cm{sup 2} ± 114.8), and the difference was statistically significant (P < 0.05). Conclusions: Our data shows excellent agreement of measured liver stiffness between GRE and SE-EPI-based sequences at 3T. Our results show the advantage of a SE-EPI MRE sequence in terms of image quality, ROI size and acquisition time with equivalent liver stiffness measurements as compared to GRE-MRE sequence.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2025

Analysis of the HAE activity in the TJ-II stellarator using a Landau closure model

The aim of this study is to analyze the stability of helical Alfvén eigenmodes (HAEs) in TJ-II discharges and the stabilizing effect of the energetic particles generated by the neutral beam injector (NBI) on pressure gradient-driven modes (PGDMs). HAE and PGDM stability is studied using the linear version of the gyro-fluid code FAR3d and the continuous structure by the STELLGAP code. First, Alfvén eigenmode (AE) and PGDM activity observed in the experiments is reproduced by the simulations, identifying unstable m/n = 4/7 − 2/3 and 7/12 − 5/8 HAEs triggered around ρ = 0.66 showing a frequency of 209 and 204 kHz, respectively, as well as 5/3 PGDM. Next, a parametric study is performed with respect to the thermal ion density and iota profile in the middle-outer plasma region to verify the robustness of the simulation results with respect to the uncertainty of experimental profiles. The analysis confirms that experimental uncertainty does not cause large deviations in the simulation results, showing the destabilization of the same HAEs for all the configurations tested. The simulations also indicate the decay of the 5/3 PGDM growth rate as the energetic particle (EP) population in the plasma increases, consistent with the experiment. Stability analysis of the n = 3, 7, 11, n = 5, 9, 13, n = 6, 10, 14, and n = 8, 12 helical families is performed with respect to the NBI operational regime for different EP energies, β as well as deposition profiles. The most unstable configuration is the radially localized on-axis NBI operation (stiff EP density profile gradients nearby the magnetic axis). Using the simulation model that reproduces the observed Alfvén activity, we extend the study to analyze NBI performance within a theoretical framework. It shows that increasing NBI voltage (which raises EP energy) leads to a degradation in NBI performance for a given power (related to EP β and their density). To achieve better NBI operation, higher voltage must be balanced with lower injection power, ensuring stable AEs while keeping the same EP β.

Plasma confinement↗
NASA NTRS1990

EZDESIT: A Computer Program for Structural Element Sizing and Vehicle Weight Prediction

A method is presented for formulating stiffness terms and thermal coefficients of fiber reinforced composite, stiffened panels for input to finite element analysis (FEA). The method is robust enough to handle panels with general cross sectional shapes, including those which are unsymmetric and/or unbalanced. New thermal coefficients are introduced to quantify panel response from through-the-thickness temperature gradients. Equations are defined for stiffness, thermal expansion, and thermal bending that consider the full complement of membrane, bending, membrane-bending coupling, and transverse shear behavior. A technique of implementing this capability with a single plane of shell finite elements using the MSC/NASTRAN™ FEA program is revealed. Finally, an example of a composite, hat-stiffened panel is included to demonstrate errors that occur when an unsymmetric panel is symmetrically formulated as traditionally done. These erroneous values and the correct ones produced from the presented method are listed for the panel's stiffness matrices, thermal expansion and bending vectors, and the forces and moments that are thermally induced from in-plane and through-the-thickness temperature gradients.

Jeffrey A Cerro↗
NASA NTRS1992

Composite Stiffened Panels: Part I - Stiffness, Thermal Expansion, and Thermal Bending Formulation for Finite Element Analysis

A method is presented for formulating stiffness terms and thermal coefficients of fiber-reinforced composite, stiffened panels for input to finite element analysis (FEA). The method is robust enough to handle panels with general cross sectional shapes, including those which are unsymmetric and/or unbalanced. New thermal coefficients are introduced to quantify panel response from through-the-thickness temperature gradients. Equations are defined for stiffness, thermal expansion, and thermal bending that consider the full complement of membrane, bending, membrane-bending coupling, and transverse shear behavior. A technique of implementing this capability with a single plane of shell finite elements using the MSC/NASTRAN(TradeMark) FEA program is revealed. Finally, an example of a composite, hat-stiffened panel is included to demonstrate errors that occur when an unsymmetric panel is symmetrically formulated as traditionally done. These erroneous values and the correct ones produced from the presented method are listed for the panel's stiffness matrices, thermal expansion and bending vectors, and the forces and moments that are thermally induced from in-plane and through-the-thickness temperature gradients.

Craig S Collier↗