Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plastic strains”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Distributed fiber optic strain sensing of bending deformation of a well mockup in the laboratory

Well integrity is critical to the safety and success of subsurface energy exploration and management, as leakage of fluids from subsurface reservoirs is often induced by compromised wells. This study investigates bending deformation of a mockup of an oil/gas well that is subjected to loads expected in deviated wells under reservoir compaction and assesses the feasibility of utilizing distributed fiber optic strain sensing to monitor the deformation. Here, a well mockup, which consists of outer and inner steel pipes with a cemented annulus, is tested under three-point bending loading and its strain and curvature development is monitored by Brillouin optical time domain reflectometry/analysis (BOTDR/A) as well as optical frequency domain reflectometry (OFDR). The primary objective of this research is to assess the strain sensing performance of newly fabricated fiber optic cables and to identify key cable characteristics which could improve the quality of distributed strain measurements with BOTDR/A. Results show that the tight-buffered cable is best suited for strain sensing as its maximum errors in the strain measurement were -36% and -24% against conventional sensors at the maximum elastic and plastic bending loads, respectively, whereas those of the non-tight-buffered cable were -45% and -71%, respectively. Similar trends were obtained in the bending curvature measurement. The detailed design of the tight-buffered cable is presented to elucidate key characteristics of such a cable, which will facilitate accurate distributed strain sensing in oil and gas wells.

42 ENGINEERING↗
DOE OSTI2022

Two–Photon Polymerized Shape Memory Microfibers: A New Mechanical Characterization Method in Liquid

Two-photon polymerization (TPP) is widely used to create 3D micro- and nanoscale scaffolds for biological and mechanobiological studies, which often require the mechanical characterization of the TPP fabricated structures. To satisfy physiological requirements, most of the mechanical characterizations need to be conducted in liquid. However, previous characterizations of TPP fabricated structures are all conducted in air due to the limitation of conventional micro- and nanoscale mechanical testing methods. In this study, a new experimental method is reported for testing the mechanical properties of TPP-printed microfibers in liquid. The experiments show that the mechanical behaviors of the microfibers tested in liquid are significantly different from those tested in air. By controlling the TPP writing parameters, the mechanical properties of the microfibers can be tailored over a wide range to meet a variety of mechanobiology applications. In addition, it is found that, in water, the plasticly deformed microfibers can return to their predeformed shape after tensile strain is released. The shape recovery time is dependent on the size of microfibers. The experimental method represents a significant advancement in mechanical testing of TPP fabricated structures and may help release the full potential of TPP fabricated 3D tissue scaffolds for mechanobiological studies.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Data-driven surrogates for high dimensional models using Gaussian process regression on the Grassmann manifold

This paper introduces a surrogate modeling scheme based on Grassmannian manifold learning to be used for cost-efficient predictions of high-dimensional stochastic systems. The method exploits subspace-structured features of each solution by projecting it onto a Grassmann manifold. This point-wise linear dimensionality reduction harnesses the structural information to assess the similarity between solutions at different points in the input parameter space. The method utilizes a solution clustering approach in order to identify regions of the parameter space over which solutions are sufficiently similarly such that they can be interpolated on the Grassmannian. In this clustering, the reduced-order solutions are partitioned into disjoint clusters on the Grassmann manifold using the eigen-structure of properly defined Grassmannian kernels and, the Karcher mean of each cluster is estimated. Then, the points in each cluster are projected onto the tangent space with origin at the corresponding Karcher mean using the exponential mapping. For each cluster, a Gaussian process regression model is trained that maps the input parameters of the system to the reduced solution points of the corresponding cluster projected onto the tangent space. Using this Gaussian process model, the full-field solution can be efficiently predicted at any new point in the parameter space. In certain cases, the solution clusters will span disjoint regions of the parameter space. In such cases, for each of the solution clusters we utilize a second, density-based spatial clustering to group their corresponding input parameter points in the Euclidean space. The proposed method is applied to two numerical examples. Here, the first is a nonlinear stochastic ordinary differential equation with uncertain initial conditions where the surrogate is used to predict the time history solution. The second involves modeling of plastic deformation in a model amorphous solid using the Shear Transformation Zone theory of plasticity, where the proposed surrogate is used to predict the full strain field of a material specimen under large shear strains.

42 ENGINEERING↗
DOE OSTI2021

Irradiation hardening and ductility loss of Eurofer97 steel variants after neutron irradiation to ITER-TBM relevant conditions

Ten Eurofer97 steel variants, produced by non-standard fabrication-processing routes and modified alloying chemistries, were studied by neutron irradiations in the high flux isotope reactor. The irradiations were performed to ITER-TBM relevant conditions of ~255–350 °C, 2.94–3.24 dpa. We quantified the irradiation-induced degradation of the steels using mechanical property tests. All the steels suffered from irradiation hardening, where a significant increase in Vickers microhardness and yield stress (σ YS ) occurred, accompanied with severe loss of tensile elongation. The extent of hardening was material dependent. For T irr = 300±30 °C, most steels showed σ YS increase in the range of ~30% to as high as ~66%, except for a low temperature tempered steel with σYS increase below 15%. Despite large losses in elongation, most failures were ductile. Significant post-necking ductility was retained with reduction in area (RA) between 65–75%, but <50% for low temperature tempered steels. The ultimate tensile stress to yield stress (σ UTS /σ YS ) ratios decreased significantly after irradiation, highlighting irradiation-induced strain hardening capacity reduction. No major effect of irradiation on the plastic instability stress (σ PIS ) and true fracture stress of the steels was observed. By comparing the tensile stresses in true stress units and with literature, the results suggest that RAFM steel designing should target materials with a large separation between σ PIS and σ YS , to ensure the materials can maintain large work hardening and uniform deformation capability after irradiation. The tensile data of the steels additionally revealed a compelling evidence of an inverse trend between the change in RA and increase in σ YS of the neutron irradiated Eurofer97 type steels.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Residual strain-stress in manganese steel railway crossing determined by synchrotron and laboratory X-rays

The residual stress within the nose of a worn manganese steel railway crossing was investigated. The macroscopic residual stress at the running surface was determined in areas with and without visual faults by laboratory X-rays. To get information about the penetration depth of the residual stress, synchrotron X-ray diffraction was used to map the distribution of local (within grain) residual strains below the running surface. The laboratory X-ray measurements revealed significant compressive stress on the rail running surface. The synchrotron measurements documented that local strains exist through the entire depth investigated. Relations between residual stress, plastic deformation, and crack formation are analysed, leading to a discussion on potentials/limitations of the two experimental methods for investigations of residual stress in large engineering components.

differentia aperture synchrotron diffraction↗
DOE OSTI2023

L1 0 Ordering in MnAl and FeNi Influenced by Magnetic Field and Strain

Due to various materials supply chain challenges, magnets free of constrained elements are attracting increasing interest. Magnetic materials such as rare-earth free FeNi and MnAl have been receiving considerable attention due to the high magnetocrystalline anisotropy and other associated magnetic properties derived from their unique chemically ordered tetragonal crystal structure, denoted as the L1 0 structure. However, synthesis of L1 0 FeNi has had limited success due to the extremely low atomic mobilities of Fe and Ni. In this work, isostructural MnAl was first studied as proxy to understand the L1 0 ordering process. Here, evidence of L1 0 ordering in FeNi derived from TEM studies is presented, where ordering was facilitated by the application of strain and magnetic field provided during thermal treatment of a severely plastically deformed FeNi alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ordering effects on deformation substructures and strain hardening behavior of a CrCoNi based medium entropy alloy

Here, a CrCoNi based medium entropy alloy with small additions of Ti, Al and Nb (denoted as (CrCoNi) 93 Al 4 Ti 2 Nb) in the as-quenched condition, exhibits tensile properties comparable to those of the equiatomic CrCoNi alloy at room temperature. Dark field transmission electron microscopy (TEM), atomic resolution high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) together with atom probe tomography (APT) show that spatially-localized long range ordering (LRO) L1 2 domains exist in this alloy. The evolution of deformation substructure with plastic deformation in this alloy was characterized using electron backscatter diffraction (EBSD), electron channeling contrast imaging (ECCI) and STEM based techniques including the recently developed weak beam dark field STEM imaging. Plastic deformation occurs by the slip of a/2<110>dislocations, which are narrowly dissociated into Shockley partial dislocations on {111} slip planes. Their dissociation distances in the (CrCoNi) 93 Al 4 Ti 2 Nb alloy are much smaller than the widths of the corresponding partials in the equiatomic CrCoNi alloy due to one or more of the minor alloying elements (Al, Ti, Nb). Dislocation slip in this alloy has a pronounced planar character. The leading dislocations in slip bands glide as pairs due to the existence of LRO domains. Multipoles were formed through the slip of dislocations with opposite signs on adjacent {111} slip planes. Those multipoles serve as building blocks for the formation of subgrain structures consisting of fine slip bands. The distances between slip bands were continuously refined during plastic deformation and dynamic refinement of slip bands plays a crucial role in strain hardening. The effects of LRO domains on planar dislocation slip, the deactivation of deformation twinning and strain hardening of this alloy are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Investigation of phase-transformation path in TiZrHf(VNbTa)x refractory high-entropy alloys and its effect on mechanical property

Refractory high-entropy alloys (RHEAs) show a great potential as structural materials due to their remarkable properties, such as high strength and thermal stability. Despite these advantages, the RHEAs still suffer from the limited ductility at room temperature. To overcome the strength-ductility trade-off, the new RHEA design strategy, transformation-induced plasticity (TRIP) HEAs, have been developed by applying metastability-engineering approach. With the newly developed RHEAs design strategy, we designed metastable RHEAs, TiZrHfVxNbxTax (x = 1.0, 0.5 and 0.1) with decrease of the β-stabilizing elements (V, Nb, and Ta) to promote the transformation-induced ductility and work-hardening capability. The proper addition of the β-stabilizing elements (V, Nb, and Ta) into the ternary equi-molar TiZrHf alloy (TiZrHfV 0.1 Nb 0.1 Ta 0.1 ) leads to stress/strain-induced phase transformation from the body-centered cubic (BCC) to hexagonal-close-packed (HCP) α' and/or orthorhombic α" phases. Here, detailed studies on the phase transformation sequence during the deformation were carried out using X-ray diffraction (XRD), electron back-scattered diffraction (EBSD), and transmission-electron microscope (TEM), which revealed the stepwise phase transformation that is responsible for the strong strain hardening behavior. Furthermore, we further suggested the new transformation-induced-plasticity (TRIP) HEA design strategy, applying the Bo-Md approach.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Evaluating the grain-scale deformation behavior of a single-phase FCC high entropy alloy using synchrotron high energy diffraction microscopy

We report although the deformation behavior of high-entropy alloys (HEAs) has been extensively studied at the macroscale, many important properties have yet to be explored for these alloys at the microscale, thus hampering accurate prediction of damage and failure processes. Synchrotron high-energy diffraction microscopy (HEDM) and fast-Fourier transform-based crystal plasticity modeling was conducted to investigate the three-dimensional (3D) grain-resolved micromechanical response for approximately 1,900 constituent grains within a single-phase FCC HEA up to 1% applied strain. The evolution of grain-resolved elastic strains, lattice reorientations, and maximum resolved shear stresses (mRSS) were evaluated to quantify elastic, yield, and fully plastic behavior. Overall, the initial critical resolved shear stress (CRSS), determined via in situ HEDM and companion modeling, was found to be > 20% higher than estimated using the classical polycrystalline Taylor factor (M = 3.06). However, a descriptive parameter based on the average grain-resolved Taylor factor ($\overline{M}$) was found to show excellent agreement with plastic yielding behavior observed within HEDM datasets. Noticeable deviations in HEDM lattice reorientations compared to both EVP-FFT simulations and classical predictions for FCC polycrystals were discovered, highlighting the complexity in correlating local lattice reorientations, Taylor, and Schmid factors with plastic response for this material at the grain-scale. Therefore, it is anticipated that the overall trends and parameter identification of 3D grain-resolved properties in this study can serve as an important foundation for continued mesoscale investigation on both well-established and newly developed Cantor-like HEAs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A new coupled multiphase flow–finite strain deformation–fault slip framework for induced seismicity

Production of hydrocarbons and water from subsurface reservoirs are known to cause permanent deformation of the reservoir and seismicity along faults both of which are detrimental to sustainable development of natural resources. Most of the prior studies on understanding fluid flow-induced plasticity and seismicity have focused on one or the other phenomenon due to the numerical difficulty associated with simultaneous modeling of the two failure phenomena because stress and deformation evolve non-linearly in both plasticity and seismicity. However, in reservoirs undergoing long-term production, plastic failure can alter the stress paths of points on a fault such that the onset, location, and magnitude of actual seismic events can no longer be predicted by a poroelastic simulation due to inaccurate stress and deformation history. We present a computational framework for coupled multiphase flow, finite strain poroplastic deformation, and dynamic fault slip and use it to understand the impact of plastic deformation on the onset, location, and magnitude of induced fault slip events. We evaluate the impact of plasticity on reservoir pressure, deformation, induced stress, and fault slip by comparing infinitesimal strain elastic and finite strain poro-elastoplastic models. For real-world applications, we consider different scenarios where the reservoir is either mechanically weaker or stronger than the caprock. We analyze the stress evolution as a function of the change in reservoir pressure to understand the role of contrast in reservoir and caprock elastic moduli on geomechanical stability. The results show that the poroplastic reservoir exhibits larger vertical deformation and delayed slip than the poroelastic reservoir after the same amount of oil production. For the same amount of pressure drop, a reservoir with a smaller modulus than the caprock displays a larger vertical displacement and an earlier onset of both plastic failure and fault slip. For a reservoir with a larger modulus than the caprock, vertical displacement is larger on the reservoir top boundary and smaller on the ground surface, and a higher pressure drop is needed to induce plastic failure and fault slip.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Bayesian calibration of a physics-based crystal plasticity and damage model

In this work, we present a model parameter calibration procedure for a physics-based crystal plasticity model. The calibration process utilizes a powerful statistics-based Bayesian calibration method. Calibration of the crystal plasticity parameters makes use of experimentally-measured data, i.e. compressive stress–strain response, from $\langle$100$\rangle$ and $\langle$123$\rangle$ single crystal copper dynamically loaded via Kolsky bar tests. The calibration of damage parameters is achieved using experimentally-measured free-surface velocity history data from plate impact test on $\langle$100$\rangle$ and $\langle$110$\rangle$ single crystal copper, which generates shock compression followed by dynamic tensile failure. A validation assessment is then carried out by comparing the calibrated model predictions and experimental measurements of the dynamic tensile damage generated in an impacted bicrystal copper plate. Lastly, a model-informed rationale for the experimentally-observed dependence of the spatial distribution of ductile damage (porosity) on crystallography is provided.

42 ENGINEERING↗
DOE OSTI2023

Progression of creep deformation from grain boundaries to grain interior in Al-Cu-Mn-Zr alloys

Creep mechanisms are studied in θ'-Al 2 Cu-strengthened Al-Cu-Mn-Zr alloys at 300 and 350°C for (i) ACMZ, a base alloy without further alloying elements and (ii) RR350, a commercial alloy with additions of Ni and Co forming distinct grain-boundary precipitates. At high stresses, creep is dominated by dislocations bypassing θ' precipitates within grains via the Orowan mechanism, as evidenced by (i) very high stress exponent (n~20-25) and (ii) α-Al and θ' lattice strains (measured via in-situ neutron diffraction) evolving during creep in a manner consistent with load transfer from the plastically-deforming α-Al matrix to elastically-deforming θ' precipitates. At intermediate stresses, both alloys exhibit a n~3 regime, where α-Al and θ' lattice strains scale near-linearly with applied stress while remaining largely unaffected by strain accumulation, indicating that Orowan looping or dislocation pile-up around θ' is now inactive within the grains. Rather, dislocation motion occurs solely in θ'-precipitate-free zones (θ'-PFZ) where high dislocation densities are observed via TEM after creep deformation. Plastic flow at θ'-PFZ and/or localized pipe diffusion are expected to enable grain-boundary sliding (GBS), which is proposed as the rate-limiting mechanism in the n~3 regime. Ni/Co-rich precipitates at RR350 grain-boundaries, with negligible θ'-PFZ around them, share load (as determined via neutron diffraction) with the α-Al matrix more effectively than θ-Al 2 Cu precipitates at ACMZ grain-boundaries, with wide surrounding θ'-PFZ. So, high creep resistance in the n~3 GBS regime of RR350 is enabled by coarsening-resistant grain-boundary precipitates, forming without concomitant development of weak θ'-PFZ, which effectively share load with the grains.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Elucidating the temperature dependence of TRIP in Q&P steels using synchrotron X-Ray diffraction, constituent phase properties, and strain-based kinetics models

Understanding the deformation-induced martensitic transformation (DIMT) is critical for interpreting the structure-property relationships that govern the performance of transformation-induced plasticity (TRIP) assisted steels. However, modern TRIP-assisted steels often exhibit DIMT kinetics that are not easily captured by existing empirical models based on bulk tensile strain. We address this challenge by combined bulk uniaxial tensile tests and in-situ high energy synchrotron X-ray diffraction, which resolved the phase volume fractions, stress-strain response, and microstructure evolution of each constituent phase. A modification of the Olson-Cohen model is implemented, which describes the martensitic transformation kinetics as a function of the estimated partitioned strain in austenite, rather than the bulk tensile strain. This DIMT kinetic model is used as a framework to clarify the root cause of an insufficiently understood toughness trough reported for TRIP-assisted steels during deformation at elevated temperatures. Here, the importance of the temperature-dependent toughness is discussed, based on the opportunity to modify deformation processes to tailor the DIMT kinetics and mechanical properties during forming and in service.

36 MATERIALS SCIENCE↗
DOE OSTI2022

High strain rate deformation of aged TRIP Ti-10V-2Fe-3Al (wt.%) examined by in-situ synchrotron X-ray diffraction

Transformation Induced Plasticity (TRIP) is a promising avenue for tailoring the work hardening response of metastable $β$ titanium (Ti) alloys. Here we show that aged TRIP Ti-10V-2Fe-3Al (wt.%) maintains higher elongations and flow stresses as strain rate increases, if phase stability and microstructural characteristics are tuned. Low temperature aging influences the matrix $β$ phase stability by $ω$ phase precipitation, which affords a promising way to impact the TRIP effect and obtain desirable mechanical properties, ranging from high damping capacity to good strength/ductility combinations. Although TRIP is active during quasi-static and dynamic testing up to 2000 s –1 , increasing aging time and/or strain rate reduces the overall propensity for the TRIP effect and extent of transformation, which occurs rapidly just at the onset of yielding. Importantly, TRIP with $ω$ phase precipitation provides interesting alloying, microstructure, and property design strategies for engineering applications like lightweight protective structures, where high strains and the need for energy absorption are encountered.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Atomistic Investigation of Plastic Deformation and Dislocation Motion in Uranium Mononitride

Uranium mononitride (UN) is a promising advanced nuclear fuel due to its high thermal conductivity and high fissile density. However, many aspects of its mechanical behavior, particularly at reactor-relevant conditions, remain unclear. In this study, molecular dynamics (MD) simulations were employed to investigate the deformation behavior and dislocation motion in UN. We found that the Kocevski potential predicts the principal slip system as $\frac{1}{2}$ $\langle110\rangle${110}, aligning with experimental data. On the other hand, the Tseplyaev potential predicts slip to primarily occur on $\frac{1}{2}$ $\langle110\rangle${111}. MD simulations of stress–strain behavior were used to estimate the nanoindentation hardness, revealing that the Kocevski potential accurately predicts hardness even though it fails to model dynamic plasticity. Complete dislocation mobility functions have been fitted for the edge and screw dislocations in both the thermally activated and phonon-drag regimes. The 300 K linear mobility of the edge dislocation using the Tseplyaev potential was found to be 817 Pa -1 ·s -1 , whereas that of the screw dislocation using the Kocevski potential was found to be 4546 Pa -1 ·s -1 . At intermediate stresses, we observed that the subsonic steady-state motion of the edge dislocation in UN is intermittently interrupted by velocity jumps, reaching the average sound velocity. Finally, the threshold Schmid stress is calculated as 179–197 MPa, which gives an upper-limit estimate of the uniaxial yield stress of polycrystalline UN of 548–603 MPa. These findings, including the fitted dislocation mobility function, provide essential input for future plasticity and dislocation dynamics models of nuclear fuels.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mobility of soil-biodegradable nanoplastics in unsaturated porous media affected by protein-corona

Soil-biodegradable plastic has been increasingly used as mulches in agriculture, which provides not only agronomical benefits but also in situ disposal and biodegradation options. However, soil-biodegradable plastic mulches inevitably fragment into micro- and nanoplastics during biodegradation, which can reside in soils or migrate into deep soils, where they may not degrade readily due to reduced microbial activity. To date, little is known about the transport of soil-biodegradable micro- and nanoplastics in soils. Here, in this study, we studied the transport of soil-biodegradable nanoplastics (∼200 nm) made of polybutylene adipate co-terephthalate (PBAT) in unsaturated sand (proxy for soil). Specifically, we studied the mobility of pristine and weathered PBAT nanoplastics in the absence and presence of proteins (positively charged lysozyme and negatively charged bovine serum albumin, pH = 7.7). We found that (1) both pristine and the weathered PBAT nanoplastics were mobile; (2) positively charged lysozyme formed protein-coronas around PBAT nanoplastics and inhibited the transport; and (3) decreased water saturation promoted the retention of PBAT nanoplastics via physical straining. These results suggest that soil-biodegradable nanoplastics fragmented from soil-biodegradable plastic mulches are mobile and may readily migrate into deep soil layers, but positively charged proteins and unsaturated flow would prevent such transport via formation of protein-corona and physical straining.

Yu, Yingxue↗
DOE OSTI2020

Progressive collapse analysis of stiffened panel structure with consideration of actual welding distortion and residual stresses

To date, inclusion of welding distortion and residual stresses in determining the strength of marine structures is based on empirical formulations due to welding complexity. In this study, a more feasible numerical model was developed from the perspective of inherent strain theory. The accuracy of this simplified numerical model was confirmed by transient thermo-elastic–plastic finite element analysis using a novel dynamic mesh refining and iterative substructure method which is 50 times faster than conventional finite element method. Welding distortion by the two numerical approaches was validated with experimental measurement. The collapse analyses indicated that the inherent strain model predicted a compressive strength and load-displacement curve with high accuracy, as validated by the transient thermo-mechanical model. The global stiffness and the ultimate strength of the welded structure would be overestimated by about 22% if no welding effects were considered. The proposed computation procedure for determining the ultimate strength of welded structures is accurate, simple, and practical and enables improvement of ultimate strength through optimization of welding joint design to reduce the welding residual stress and distortion.

Huang, Hui↗
DOE OSTI2020

Notch-tensile behavior of Al0.1CrFeCoNi high entropy alloy

Notch-tensile behavior of an Al0.1CrFeCoNi high entropy alloy (Al0.1-HEA) was studied using V-notch geometry and with local strain mapping using digital image correlation (DIC). A notch-strength ratio of 1.51 indicated notch strengthening. Further analysis of stress-strain response supplemented with microstructural analysis revealed that, while the presence of notch results in strengthening due to work hardening by twinning induced plasticity, the notch also contributes strongly to geometrical softening in the non-uniform ductility regime, and accounts for the onset of failure. The strain localization behavior of Al0.1-HEA due to the presence of V-notch was compared with three conventional alloys: Inconel 625 nickel-based superalloy, 304 stainless steel and Ti–6Al–4V titanium alloy. The study revealed that the nature of notch widening with increasing strain was dependent on material characteristics. The extent of notch widening impacted the local strain field and stress distribution, thereby influencing the propensity for crack initiation and growth. The experimental results were verified by finite element analysis.

Sinha, Subhasis↗