Ultrafast Measurements of Mode-Specific Deformation Potentials of Bi 2 Te 3 and Bi 2 Se 3
Not Available
SEARCH · Search NASA
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.
Not Available
Level density ρ(E,Q) is derived within the micro-macroscopic approximation (MMA) for a system of strongly interacting Fermi particles with the energy E and additional integrals of motion Q , in line with several topics of the universal and fruitful activity of A. S. Davydov. Within the extended Thomas Fermi and semiclassical periodic orbit theory beyond the Fermi-gas saddle-point method, we obtain ρ ∝ I ν (S)/S ν , where I ν (S) is the modified Bessel function of the entropy S . For small shell-structure contribution, one finds ν = κ/2 + 1, where κ is the number of additional integrals of motion. This integer number is a dimension of Q, Q = { N, Z , …} for the case of two-component atomic nuclei, where N and Z are the numbers of neutrons and protons, respectively. For much larger shell structure contributions, one obtains ν = κ /2 + 2. The MMA level density ρ reaches the well-known Fermi gas asymptote for large excitation energies and the finite micro-canonical combinatoric limit for low excitation energies. Further, the additional integrals of motion can also be the projection of the angular momentum of a nuclear system for nuclear rotations of deformed nuclei, number of excitons for collective dynamics, and so on. Fitting the MMA total level density ρ( E , Q) for a set of the integrals of motion Q = { N, Z }, to experimental data on a long nuclear isotope chain for low excitation energies, one obtains the results for the inverse level-density parameter K , which differs significantly from those of neutron resonances due to shell, isotopic asymmetry, and pairing effects.
SWIFT (Silicon-water in Familiar Template) foams are silicone foams that have a continuous porosity without having any obvious spherical or ellipsoidal pore shapes. Like other foamed rubber, SWIFT foams can experience large reversible deformations. Because of their microstructure, SWIFT foams have a stiffering behavior in compression that is very abrupt compared to most other foamed rubbers. Work by Benedikt reported by Miller showed that fitting compressive test data with the CHIPFoam model in its current form is impractical and inaccurate. This report describes an extension made to CHIPFoam to allow the modeling of SWIFT foams. While the extension allows the modeling of SWIFT foams, it also extends the range of porosity of foams that can be modelled with CHIPFoam to include very high porosity (>99%) foamed rubber. The previous model was unstable for porosities greater than 70%. The changes made to CHIPFoam include the introduction of a critical porosity at which the compressible Danielsson strain energy function that couples volumetric and deviatoric mechanical behavior begins to increase with compression. This new feature also required the introduction of a third term in the multiplicative decomposition of the relative volume, J . The organization of this report is described here. The first section, SWIFT Foam Structure and Mechanical Behavior, describes the structure of SWIFT foams as contrasted to other foamed rubbers. It also contrasts the mechanical behavior of SWIFT foams to other more traditional foamed rubbers. The second section, CHIPFoam Modifications, describes the specific modifications done to be able to modify the delayed rapid compressive stiffening seen in SWIFT foams. The third section, Results, demonstrates the use of the modified model to fit compressive loading response data from a SWIFT foam sample. The fourth section, Conclusions and Further Work, puts the work in perspective and presents likely further work.
This paper presents the first measurements of the azimuthal anisotropy coefficients $v_{n}$, which quantify the $n^{\mathrm{th}}$-order Fourier modulation of charged-particle azimuthal distributions, for $n=2$-4 in $\sqrt{s_{\mathrm{NN}}}=5.36$ TeV $\mathrm{^{16}O}+\mathrm{^{16}O}$ and $\mathrm{^{20}Ne}+\mathrm{^{20}Ne}$ collisions recorded with the ATLAS detector at the Large Hadron Collider in 2025. The $v_{n}$ coefficients are measured as a function of transverse momentum ($p_{\mathrm{T}}$), collision centrality, and event multiplicity. They are extracted using two complementary methods: two-particle correlations with a template-fit subtraction of short-range non-flow contributions, and four-particle subevent cumulants, which intrinsically suppress non-flow effects and provide sensitivity to flow fluctuations. The results show a clear hierarchy $v_{2} > v_{3} > v_{4}$ and a non-monotonic dependence on $p_{\mathrm{T}}$, reaching a maximum around 2 GeV, consistent with trends observed in heavy-ion collisions. Detailed comparisons between the two collision systems reveal an enhanced $v_{2}$ in central $\mathrm{^{20}Ne}+\mathrm{^{20}Ne}$ collisions, consistent with theory expectations based on the predicted prolate deformation of neon nuclei, in contrast to the slightly tetrahedral structure predicted for oxygen. The four-particle cumulant results highlight strong event-by-event fluctuations and provide the greatest sensitivity to nuclear shape effects. These measurements can place new constraints on the initial geometry and the hydrodynamic response in light-ion collisions, offering valuable input for models of nuclear structure.
We introduce a bosonic ambitwistor string theory in AdS space. Even though the theory is anomalous at the quantum level, one can nevertheless use it in the classical limit to derive a novel formula for correlation functions of boundary CFT operators in arbitrary space-time dimensions. The resulting construction can be treated as a natural extension of the CHY formalism for the flat-space S-matrix, as it similarly expresses tree-level amplitudes in AdS as integrals over the moduli space of Riemann spheres with punctures. These integrals localize on an operator-valued version of scattering equations, which we derive directly from the ambitwistor string action on a coset manifold. As a testing ground for this formalism we focus on the simplest case of ambitwistor string coupled to two cur- rent algebras, which gives bi-adjoint scalar correlators in AdS. In order to evaluate them directly, we make use of a series of contour deformations on the moduli space of punctured Riemann spheres and check that the result agrees with tree level Witten diagram computations to all multiplicity. We also initiate the study of eigenfunctions of scattering equations in AdS, which interpolate between conformal partial waves in different OPE channels, and point out a connection to an elliptic deformation of the Calogero-Sutherland model.
Biofouling disrupts surface functionality and integrity of engineered substrates. A variety of natural materials such as plant leaves and insect wings have evolved sophisticated physical mechanisms capable of preventing biofouling. Over the past decade, several reports have pinpointed nanoscale surface topography as an important regulator of the surface adhesion and growth of bacteria. Although artificial nanoengineered features have been used to create bactericidal materials that kill adhered bacteria, functional surfaces capable of synergistically providing anti-biofouling and bactericidal properties remain to be developed. Furthermore,fundamental questions pertaining to the need for intrinsic hydrophobicity to achieve bactericidal performance or the crucial role played by structure length scale (nano vs. micro), remain to be answered. Here, we demonstrate highly scalable, cost effective, and efficient nanoengineered multifunctional surfaces that possess both anti-biofouling and bactericidal properties on industrially relevant copper (Cu) and aluminum (A1) substrates. We characterize biofouling and bactericidal performance using a combination of scanning electron microscopy (SEM), atomic force microscopy (AFM), live-dead bacterial staining and imaging, as well as solution phase measurements of bacterial viability. SEM results showed that nanostructures created on both Cuand Al were capable of physical deformation of adhered E. coli . Bacterial viability measurements on both Cu and Al indicated a complex interaction between the anti-biofouling and bactericidal nature of these materials and their surface topography, chemistry, and structure. We found that nano-length structures, as compared to micro-length, provide improved bactericidal properties,and that increased hydrophobicity greatly decreased the of adhered bacteria while also modestly increasing the surfaces killing capacity. This study provides additional insights into design guidelines for materials that are not only bactericidal, but also anti-biofouling, using a simple and economic method.
Both separate and combined wind tunnel and vibration shaker tests were conducted on two structural panels representative of the Shuttle orbiter in the NASA LaRC 8-foot transonic pressure tunnel (TPT) to determine the effects of combined loads on the thermal protection system (TPS). The primary objective of this test was to provide a combined full-scale load environment and realistic time history of the dynamic pressures, Mach numbers (through transonic), and dynamic structural responses of these panels. The panels were selected from orbiter locations where interactive load sources such as aerodynamic shock waves, turbulent boundary layers, strut-induced vorticity, and substrate deformation combined to provide high bonding loads between the TPS and the orbiter structure. The test panels for this program were highly instrumented with static and dynamic pressure gages, accelerometers, deflectometers, strain gages, Schileren and high-speed photography, and special instrumentation necessary to determine TPS/structure interface loads and tile motions. Two test specimens of each orbiter panel were utilized. Both were high-fidelity representations of the selected orbiter location. The first panel (calibration panel) for each test was utilized to calibrate the tunnel and shaker system as well as provide a data-gathering source for applied static and dynamic loads. The second panel (test panel) for each test was subjected to flight time exposures only so that TPS characteristics as a function of flight mission exposure could be determined. These tests were the last in a series that was suce6ssfully completed and necessary to certify the TPS prior to the first launching of the Shuttle.
When plasma-beta greater than 1 then the gas pressure dominates over the magnetic pressure. This ratio as a function along the coronal magnetic field lines varies from beta greater than 1 in the photosphere at the base of the field lines, to beta much less than 1 in the mid-corona, to beta greater than 1 in the upper corona. Almost all magnetic field extrapolations do not or cannot take into account the full range of beta. They essentially assume beta much less than 1, since the full boundary conditions do not exist in the beta greater than 1 regions. We use a basic parametric representation of the magnetic field lines such that the field lines can be manipulated to match linear features in the EUV and SXR coronal images in a least squares sense. This research employs free-form deformation mathematics to generate the associated coronal magnetic field. In our research program, the complex magnetic field topology uses Parametric Transformation Analysis (PTA) which is a new and innovative method to describe the coronal fields that we are developing. In this technique the field lines can be viewed as being embedded in a plastic medium, the frozen-in-field-line concept. As the medium is deformed the field lines are similarly deformed. However the advantage of the PTA method is that the field line movement represents a transformation of one magnetic field solution into another magnetic field solution. When fully implemented, this method will allow the resulting magnetic field solution to fully match the magnetic field lines with EUV/SXR coronal loops by minimizing the differences in direction and dispersion of a collection of PTA magnetic field lines and observed field lines. The derived magnetic field will then allow beta greater than 1 regions to be included, the electric currents to be calculated, and the Lorentz force to be determined. The advantage of this technique is that the solution is: (1) independent of the upper and side boundary conditions, (2) allows non-vanishing magnetic forces, and (3) provides a global magnetic field solution, which contains high- and low-beta regimes and maximizes the similarity between the field lines structure and all the coronal images of the region. The coronal image analysis is crucial to the investigation and for the first time these images can be exploited to derive the coronal magnetic field in a well-posed mathematical formulation. This program is an outgrowth of an investigation in which an extrapolated potential field was required to be "inflated" in order to have the field lines match the Yohkoh/SXT images. The field lines were radially stretched resulting in a better match to the coronal loops of an active region. The PTA method of radial and non-radial deformations of field lines to provide a match to the EUV/SXR images will be presented.
Conducting biological research in space requires consideration be given to isolating appropriate control parameters. For in vitro cell cultures, numerous environmental factors can adversely affect data interpretation. A biological response attributed to microgravity can, in theory, be explicitly correlated to a specific lack of weight or gravity-driven motion occurring to, within or around a cell. Weight can be broken down to include the formation of hydrostatic gradients, structural load (stress) or physical deformation (strain). Gravitationally induced motion within or near individual cells in a fluid includes sedimentation (or buoyancy) of the cell and associated shear forces, displacement of cytoskeleton or organelles, and factors associated with intra- or extracellular mass transport. Finally, and of particular importance for cell culture experiments, the collective effects of gravity must be considered for the overall system consisting of the cells, their environment and the device in which they are contained. This does not, however, rule out other confounding variables such as launch acceleration, on orbit vibration, transient acceleration impulses or radiation, which can be isolated using onboard centrifuges or vibration isolation techniques. A framework is offered for characterizing specific cause-and-effect relationships for gravity-dependent responses as a function of the above parameters.
Ionizable polymers form dynamic networks with domains controlled by two distinct energy scales, ionic interactions and van der Waals forces; both evolve under elongational flows during their processing into viable materials. A molecular level insight of their nonlinear response, paramount to controlling their structure, is attained by fully atomistic molecular dynamics simulations of a model ionizable polymer, polystyrene sulfonate. As a function of increasing elongational flow rate, the systems display an initial elastic response, followed by an ionic fraction-dependent strain hardening, stress overshoot, and eventually strain-thinning. As the sulfonation fraction increases, the chain elongation becomes more heterogeneous. Finally, flow-driven ionic assembly dynamics that continuously break and reform control the response of the system.
The development of more efficient catalytic processes and new approaches to control catalytic activity and selectivity are central to the realization of more selective, atom economic, and energy efficient routes to value added chemicals and polymers. The reactivity and selectivity of a transition metal catalyst is intimately related to the ligand-sphere geometry and, in many cases, the ideal ligand geometry for one step of a catalytic cycle is poorly matched to the ideal ligand geometry for another, resulting in sub-optimal efficiency. Macroscopic mechanical forces are both large, potentially much larger than interatomic forces, and are directional and localized to an extent that differentiates them from other forms of energy input such as heat or light. As such, mechanical force represents a heretofore untapped approach to modulate catalyst geometry, with the potential to reversibly modulate catalyst geometry on the timescale of catalytic turnover or monomer enchainment. This project has addressed the fundamental challenges in material-to-molecule strain coupling associated with the development of a new class of mechanically responsive catalysts (mechanocatalysts) in which active organotransition metal catalysts are strategically embedded in a flexible polymer network such that application of external mechanical force (stretching or deformation) leads to modulation of catalyst geometry, and hence reactivity and selectivity. Our efforts during the tenure of this grant were directed toward the elucidation of force-reactivity relationships of elementary transformations that occur within the first coordination sphere of a transition metal complex employing stiff stilbene photoswitches tethered to a flexible bidentate phosphine ligand derived from MeOBiphep as molecular force probes which provide a range of compressive and extension forces to the coupled transition metal complex depending on the geometry of the stiff stilbene and length of the tethering chains. During the tenure of this grant, we have quantified the rate of C(sp 2 )-C(sp 2 ) reductive elimination from platinum(II) diaryl complexes containing bis(phosphine) force probe ligands as a function of mechanical force; compressive forces decreased the rate of reductive elimination whereas extension forces increased the rate relative to the strain-free MeOBiphep complex with a 3.4-fold change in rate over a ~290 pN range of restoring forces. In a similar manner, we have quantified the rate of oxidative addition of bromobenzene to low-ligated palladium(0) complexes containing force probe ligands as a function of mechanical force; compressive forces increase the rate of oxidative addition, whereas tensile forces decrease the rate with a ~6 fold change in rate across ~340 pN of force applied to the complexes. In both cases, experimental and computational analyses argue strongly against any significant force-induced perturbation of ground state geometry within the first coordination sphere of the reactant complexes. Rather, the force/rate behavior observed for these transformations across these ranges of forces is attributed to the coupling of force to the nuclear motion comprising the reaction coordinates for reductive elimination and oxidative addition. These results together inform the development of catalysts whose activity can be tuned by an external force that is adjusted within a catalytic cycle and suggest opportunities to experimentally map geometry changes associated with reactions in transition metal complexes and potential strategies for force-modulated catalysis.
The lift force on a deformable liquid sphere moving in steady, plane Poiseuille-Stokes flow and subjected to an external body force is calculated. The results are obtained by seeking a solution to Stokes' equations for the motion of the liquids inside and outside the slightly perturbed sphere surface, as expansions valid for small values of the ratio of the Weber number to the Reynolds number. When the ratio of the drop and external fluid viscosities is small, the lift exerted on a neutrally buoyant drop is found to be approximately one-tenth of the magnitude of the force reported by Wohl and Rubinow acting on the same drop in unbounded Poiseuille flow in a tube. The resultant trajectory of the drop is calculated and displayed as a function of the external body force.
The previously developed dynamic theory of growth and detachment of vapour bubbles under conditions of nucleate pool boiling is applied to study motion and deformation of a bubble evolving at a single nucleation site. The bubble growth is presumed to be thermally controlled, and two components of heat transfer to the bubble are accounted of: the one from the bulk of surrounding liquid and the one due to heat conduction across a liquid microlayer formed underneath the bubble. Bubble evolution is governed by the buoyancy and an effective surface tension force, both the forces making the bubble centre of mass move away from the wall and, thus, assisting its detachment. Buoyancy-controlled and surface-tension-controlled regimes are considered separately in a meticulous way. The duration of the whole process of bubble evolution till detachment, the rate of growth, and the bubble departure size are found as functions of time and physical and operating parameters. Some repeatedly observed phenomena, such as an influence of gravity on the growth rate, are explained. Inferences of the model agree qualitatively with available experimental evidence, and conclusions pertaining to the dependence on gravity of the bubble radius at detachment and the whole time of the bubble development when being attached to the wall are confirmed quantitatively.
Unreinforced woven fabrics have been implemented in a variety of high performance applications, including body armor, deployable structures, and as the reinforcement material in composites. Multiscale modeling techniques have significantly improved the capabilities of simulation-based tools to capture fabric mechanics efficiently and accurately, but often lack in their prediction of failure and require pairing with finite element analysis (FEA) software, limiting their application to the design of ‘fit-for-purpose’ materials. NASA’s Multiscale Analysis Tool (NASMAT) is a standalone multiscale program that has been traditionally used in the analysis of reinforced composites materials. More recently, it has been amended to simulate unreinforced fabric behavior by allowing the geometric state of the tows to change with applied loading due to the lack of a reinforcement material, such as the matrix seen in composites. Previous work has shown the ability of NASMAT to capture nonlinear macroscale behavior by predicting geometric changes in the state of each subcell as a function of the applied loading and allowing each subcell in the analysis to rotate according to these predicted changes, as well as predict nonlinear behavior due to the fiber breakage failure mode. In this work, the capability of predicting the onset and propagation of failure in plain woven fabrics in NASMAT is presented for the yarn pullout failure mode, which occurs when a fabric is loaded at an off-axis angle relative to the warp of weft tow direction. Yarn pullout behavior is initiated by determining the applied load in which the shear resistance of the contact area between yarn families is overcome. When failure is initiated, yarn pullout is determined to have occurred when the applied displacement, calculated from global strain, exceeds the deformed position of a given contact points between yarn families, determined from pin-joint kinematics. Contact points where pullout has occurred contribute to a global damage parameter used to modify the homogenized stiffness of the fabric, resulting in nonlinear behavior observed at the macroscale. The off-axis loading behavior and yarn pullout failure theory have been developed and implemented into NASMAT such that users can simulate off-axis tensile behavior of fabrics in a single, standalone multiscale tool. Simulations are compared to uniaxial tensile tests at various off-axis angles to demonstrate the capability of the tool in its prediction of both the onset of failure at each off-axis angles and the stress-strain behavior as failure progresses.
Dislocation slip, deformation twinning, phase transformations, and fast fracture are energy dissipation mechanisms that accommodate mechanical loading in materials. The energetically unfavorable formation of dislocation attributes to the ease of cracking and the low damage tolerance observed in superhard ceramics, notably boron carbide. Here, this work demonstrates that room temperature dislocation slip can be enabled in boron carbide by altering its chemistry through Al doping. The activation of dislocation slip is mechanistically explained by quantum mechanics simulations and electron microscopy, which indicate that strain energy is released through basal icosahedral slip facilitated by icosahedral rotation and chain bond breaking and reconfiguring. The new insight gained through this work suggests that atomic doping could be an effective strategy to tune deformation mechanisms in boron carbide, which provides a significant potential for limiting amorphization and catastrophic failure, and opens a new strategy to enhance damage tolerance in brittle ceramics.
The Mechanical Threshold Stress (MTS) model provides excellent predictive capabilities for the material constitutive response for a wide range of temperatures and strain rates. However, the MTS model fails to capture the rapidly increasing yield stress at high strain rate behavior as the deformation controlling mechanism transitions from thermal activation to drag mechanisms, only capturing the linear behavior. Further, the model typically over predicts the flow stress behavior at yield and post yield due to its use of a constant work hardening rate parameter derived from the stress–strain response at constant saturation stress. An alternative approach to fitting portions of the MTS model is investigated and mathematical models are developed to address these issues. The results show that with appropriate experimental data, the mechanical threshold stress and work hardening rate parameters within the MTS model can quite easily and accurately be modified to extend applicability to high strain rate behavior and more accurately model the initial flow stress behavior at early work hardening rates without modification of the functions core to the MTS model itself.
Natural gas furnaces are the most common space heating equipment in the U.S. residential and commercial building markets. However, current residential natural gas condensing furnaces generate substantial acidic condensate as well as significant emissions of sulfur oxides (SOx), nitrogen oxides (NOx), carbon monoxide (CO), hydrocarbons (HC), and methane (CH 4 ) contributing to environmental degradation of air, water, and soil. This report describes a novel solution to reduce the environmental impact of natural gas condensing furnaces based on the technology of a monolithic acidic gas reduction (AGR) catalyst for SOx trapping, NOx redox to nitrogen, and oxidation of formic acid, CO, HC, and CH 4 . The AGR technology offers the following benefits: (1) a neutral furnace condensate with a pH of ~7, allowing its safe release into the sewer system thus eliminating a second drainage system; (2) trapping and removing nearly all SOx emissions; (3) NOx emissions nearly at nearly 1-2 ng/J, more than 95% lower than new emissions standards in California; (4) the use of a low-cost heat exchanger as a condensing heat exchanger (HX) since the condensate is not acidic, avoiding the need for expensive stainless steel alloys; and (5) unburnt fuel energy recovery to boost efficiency.The AGR component and AGR-enabled furnace performance were broadly tested to determine their effects on long-term reliability and durability, as well as SOx storage and regeneration activity. The AGR regeneration does not impair the performance in achieving neutral condensate and ultra-low NOx emissions, and the AGR catalyst subjected to regeneration activities continued to function well and achieved slightly better annual fuel utilization efficiency (AFUE). The 400-hour reliability and durability test of the retrofitted condensing furnace with the AGR component shows that the furnace unit still achieves a neutral furnace condensate with a pH of ~7 and enables 0~3 ng/J of NOx emissions. However, the 400-hour operation slightly degraded the AFUE because of soot particle accumulation caused by frequent incomplete combustion owing to inappropriate condensate drainage during testing. Thus, proper condensate drainage is critical for AGR-enabled furnaces. Furthermore, neutron computed tomography was employed to survey the aged AGR component and demonstrate high-resolution 2D and 3D representations for the nondestructive diagnosis of the AGR component. The tomography showed that the AGR component did not deform or suffer broken AGR channels. A new AGR catalyst with low precious metal loading was preliminarily explored to identify a pathway of optimizing AGR material loading and maximizing acidic gas reduction at low cost. The new AGR component can reduce precious metal loading by 38% and still achieve neutral condensate and ultralow NOx emissions. The furnace with the AGR component of low Pt/Rh loading enables a maximum AFUE of 97%, which is meaningfully higher than the original furnace. Long-duration testing for the furnace enabled with the low precious metal loading AGR component will be vital in future research. Although the current work demonstrates a proof of concept for the AGR-enabled furnace, the AGR assembly needs to be optimized and integrated into the design of new OEM furnace products. Furthermore, the AGR technology can be applied not only for residential gas furnaces, but also for commercial rooftop units, gas heat pumps, gas-fired water heaters, combustion boilers, and other systems.
Finite-element fracture simulation methodology predicts the residual strength of damaged aircraft structures. The methodology uses the critical crack-tip-opening-angle (CTOA) fracture criterion to characterize the fracture behavior of the material. The CTOA fracture criterion assumes that stable crack growth occurs when the crack-tip angle reaches a constant critical value. The use of the CTOA criterion requires an elastic- plastic, finite-element analysis. The critical CTOA value is determined by simulating fracture behavior in laboratory specimens, such as a compact specimen, to obtain the angle that best fits the observed test behavior. The critical CTOA value appears to be independent of loading, crack length, and in-plane dimensions. However, it is a function of material thickness and local crack-front constraint. Modeling the local constraint requires either a three-dimensional analysis or a two-dimensional analysis with an approximation to account for the constraint effects. In recent times as the aircraft industry is leaning towards monolithic structures with the intention of reducing part count and manufacturing cost, there has been a consistent effort at NASA Langley to extend critical CTOA based numerical methodology in the analysis of integrally-stiffened panels.In this regard, a series of fracture tests were conducted on both flat and curved aluminum alloy integrally-stiffened panels. These flat panels were subjected to uniaxial tension and during the test, applied load-crack extension, out-of-plane displacements and local deformations around the crack tip region were measured. Compact and middle-crack tension specimens were tested to determine the critical angle (wc) using three-dimensional code (ZIP3D) and the plane-strain core height (hJ using two-dimensional code (STAGS). These values were then used in the STAGS analysis to predict the fracture behavior of the integrally-stiffened panels. The analyses modeled stable tearing, buckling, and crack branching at the integral stiffener using different values of critical CTOA for different material thicknesses and orientation. Comparisons were made between measured and predicted load-crack extension, out-of-plane displacements and local deformations around the crack tip region. Simultaneously, three-dimensional capabilities to model crack branching and to monitor stable crack growth of multiple cracks in a large thick integrally-stiffened flat panels were implemented in three-dimensional finite element code (ZIP3D) and tested by analyzing the integrally-stiffened panels tested at Alcoa. The residual strength of the panels predicted from STAGS and ZP3D code compared very well with experimental data. In recent times, STAGS software has been updated with new features and now one can have combinations of solid and shell elements in the residual strength analysis of integrally-stiffened panels.