Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix factorization”

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 307 records · Page 17

Overview of Potential Hazards in Electric Aircraft Charging Infrastructure

With increasing efforts in electrification of Advanced Air Mobility and electric aircraft, there is a growing need to install new infrastructure to support them. This report highlights potential hazards (non exhaustive) that the installation must be aware of and be prepared to mitigate with increased electrical equipment on site. The hazards can be due to natural causes as well as non natural causes. The report aims to provide starting guidelines for awareness of these hazards and to help plan for mitigation of the same.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Comprehensive Review of 2$β$ Decay Half-Lives

The double-beta (2β)-decay is the rarest nuclear physics process, and its experimental half-lives (T 1/2 ) exceed the age of the Universe from nine to fourteen orders of magnitude. Double-beta decay was observed, and its half-life was measured in 14 parent nuclei using direct, radiochemical, and geochemical methods. The decay observables are analyzed using the Evaluated Nuclear Structure Data File (ENSDF) procedures, and the recommended T 1/2 were deduced. Using the calculated values of phase factors, the effective nuclear matrix elements were extracted and compared with available data. Thousands of theoretical and experimental works have been dedicated to these topics in the last 85 years, and we present two data sets of recommended values to encapsulate the results.

2β-decay↗

On the invertibility of linear multivariable systems

Some results concerning invertibility of a class of linear, time-invariant systems are presented. It is shown that by an appropriate factorization of the transfer matrix of such a system, the problem of checking its invertibility can be reduced to that of checking the invertibility of a lower-order system. A sufficient condition for invertibility is also obtained.

Patel, R. V.↗

A multilevel approach for minimum weight structural design including local and system buckling constraints

A rational multilevel approach for minimum weight structural design of truss and wing structures including local and system buckling constraints is presented. Overall proportioning of the structure is achieved at the system level subject to strength, displacement and system buckling constraints, while the detailed component designs are carried out separately at the component level satisfying local buckling constraints. Total structural weight is taken to be the objective function at the system level while employing the change in the equivalent system stiffness of the component as the component level objective function. Finite element analysis is used to predict static response while system buckling behavior is handled by incorporating a geometric stiffness matrix capability. Buckling load factors and the corresponding mode shapes are obtained by solving the eigenvalue problem associated with the assembled elastic stiffness and geometric stiffness matrices for the structural system. At the component level various local buckling failure modes are guarded against using semi-empirical formulas. Mathematical programming techniques are employed at both the system and component level.

Schmit, L. A., Jr.↗

An exact zooming method for finite element analyses

An exact zooming technique which employs static condensation and exact structural reanalysis methods was developed. Successive application of static condensation reduces the system to one that is only associated with the degrees of freedom (dof) of the original model. Application of an exact static reanalysis technique permits the displacements at the dof of the original model that are contained in the zoomed portion of the structure to be obtained first. The response external to the zoom, as well as the response of additional dof within various levels of zooming, is computed. With the triangular factor of the stiffness matrix of the original system available, this approach involves only the solution of a system of equations of small order.

Hirai, I.↗

Recent developments in the design, testing and impact-damage tolerance of stiffened composite panels

The structural technology of laminated filamentary-composite stiffened-panel structures under combined in-plane and lateral loadings is discussed. Emphasis is on analyzing the behavior of the structures under load, determining appropriate structural proportions for weight efficient configurations, and effects of impact damage and geometric imperfections on structural performance. Experimental data on buckling of panels under in-plane compression validate the analysis and sizing methods, and illustrate structural performance and efficiency obtained from representative structures. It is shown that the strength of panels under in-plane compression can be degraded by low-velocity impact damage, and data are presented which indicate that the matrix is a significant factor influencing tolerance to impact damage.

Williams, J. G.↗

Task granularity studies on a many-processor CRAY X-MP

A hybrid granularity model is proposed for general concurrent solution. It is applied to the triangular factorization of a dense matrix ranging in size from 4 to 1024. Concurrency is achieved at two levels: (1) with small (micro) task granularity and (2) with large (blocked) task granularity. Relevance to a many-processor CRAY X-MP is demonstrated by simulation.

Calahan, D. A.↗

An implicit and stiffly stable finite element CFD algorithm for unsteady aerodynamics

A stable and accurate finite element CFD algorithm for hyperbolic/incompletely parabolic conservation law systems is described and verified. It combines a Taylor weak statement FEM, an optimal implicit Runge-Kutta time integration algorithm, and a matrix tensor product approximate factorization linear algebra procedure. The results of computational experiments show that the developed algorithm is robust.

Baker, A. J.↗

Microstructural non-uniformity as related to unidirectional compression performance

The compression performance of unidirectional composite laminates (UDCLs) is related not only to intrinsic or local strength, but also to a sample's ability to resist propagaton of initial damage. This paper introduces a simple shear lag model which is applied to a compression loaded UDCL with initial damage. Among the information obtained are axial overload and matrix layer shear strain factors in the neighborhood of the damage zone. These factors indicate that propagation may occur with relatively little initial damage and emphasize the need for a more complete understanding of compression behavior than currently exists.

Majors, Patrick S.↗

Expendable launch vehicle studies

Analytical support studies of expendable launch vehicles concentrate on the stability of the dynamics during launch especially during or near the region of maximum dynamic pressure. The in-plane dynamic equations of a generic launch vehicle with multiple flexible bending and fuel sloshing modes are developed and linearized. The information from LeRC about the grids, masses, and modes is incorporated into the model. The eigenvalues of the plant are analyzed for several modeling factors: utilizing diagonal mass matrix, uniform beam assumption, inclusion of aerodynamics, and the interaction between the aerodynamics and the flexible bending motion. Preliminary PID, LQR, and LQG control designs with sensor and actuator dynamics for this system and simulations are also conducted. The initial analysis for comparison of PD (proportional-derivative) and full state feedback LQR Linear quadratic regulator) shows that the split weighted LQR controller has better performance than that of the PD. In order to meet both the performance and robustness requirements, the H(sub infinity) robust controller for the expendable launch vehicle is developed. The simulation indicates that both the performance and robustness of the H(sub infinity) controller are better than that for the PID and LQG controllers. The modelling and analysis support studies team has continued development of methodology, using eigensensitivity analysis, to solve three classes of discrete eigenvalue equations. In the first class, the matrix elements are non-linear functions of the eigenvector. All non-linear periodic motion can be cast in this form. Here the eigenvector is comprised of the coefficients of complete basis functions spanning the response space and the eigenvalue is the frequency. The second class of eigenvalue problems studied is the quadratic eigenvalue problem. Solutions for linear viscously damped structures or viscoelastic structures can be reduced to this form. Particular attention is paid to Maxwell and Kelvin models. The third class of problems consists of linear eigenvalue problems in which the elements of the mass and stiffness matrices are stochastic. dynamic structural response for which the parameters are given by probabilistic distribution functions, rather than deterministic values, can be cast in this form. Solutions for several problems in each class will be presented.

Bainum, Peter M.↗

Scaling Laws Applied to a Modal Formulation of the Aeroservoelastic Equations

A method of scaling is described that easily converts the aeroelastic equations of motion of a full-sized aircraft into ones of a wind-tunnel model. To implement the method, a set of rules is provided for the conversion process involving matrix operations with scale factors. In addition, a technique for analytically incorporating a spring mounting system into the aeroelastic equations is also presented. As an example problem, a finite element model of a full-sized aircraft is introduced from the High Speed Research (HSR) program to exercise the scaling method. With a set of scale factor values, a brief outline is given of a procedure to generate the first-order aeroservoelastic analytical model representing the wind-tunnel model. To verify the scaling process as applied to the example problem, the root-locus patterns from the full-sized vehicle and the wind-tunnel model are compared to see if the root magnitudes scale with the frequency scale factor value. Selected time-history results are given from a numerical simulation of an active-controlled wind-tunnel model to demonstrate the utility of the scaling process.

Pototzky, Anthony S.↗

Direct Assembly of Modified Proteins on Carbon Nanotubes in an Aqueous Solution

Carbon nanotubes (CNTs) have superior mechanical and electrical properties that have opened up many potential applications. However, poor dispersibility and solubility, due to the substantial van der Waals attraction between tubes, have prevented the use of CNTs in practical applications, especially biotechnology applications. Effective dispersion of CNTs into small bundles or individual tubes in solvents is crucial to ensure homogeneous properties and enable practical applications. In addition to dispersion of CNTs into a solvent, the selection of appropriate solvent, which is compatible with a desired matrix, is an important factor to improve the mechanical, thermal, optical, and electrical properties of CNT-based fibers and composites. In particular, dispersion of CNTs into an aqueous system has been a challenge due to the hydrophobic nature of CNTs. Here we show an effective method for dispersion of both single wall CNTs (SWCNTs) and few wall CNTs (FWCNTs) in an aqueous buffer solution. We also show an assembly of cationized Pt-cored ferritins on the well dispersed CNTs in an aqueous buffer solution.

Kim, Jae-Woo↗

Improved Plate and Beam Models for Thermoviscoelastic Constitutive Modeling of Composites

The effective properties of composites are influenced by the time-dependent behavior of polymer matrices very sensitive to changes in temperature. Improved plate and beam models are required to efficiently design, and simulate composite structures when the long-term performance of large anisotropic composite structures is the matter of interest. In this work, mechanics of structure genome (MSG) is used to con-struct linear thermoviscoelastic plate and beam models that can homogenize three-dimensional heterogeneous materials made of constituents with time- and temperature-dependent behavior. The formulation derives the transient strain energy based on integral formulation for thermorheologically simple materials subject to finite temperature changes with the restriction that the strain is small. The reduced time parameter is introduced to relate the time-temperature dependency of the anisotropic material by means of master curves at reference conditions. The new formulation has been implemented in SwiftCompTM, a general-purpose multiscale constitutive modeling code based on MSG. Experimental data and three-dimensional direct numerical simulations of thin-ply high-strain composites (TP-HSC) using a commercial finite element analysis (FEA) package are conducted to verify the accuracy of SwiftCompTM results. The paper also analyzes the relationship between the shift factor of the polymer matrix and the temperature dependencies of the effective beam properties.

Finite element analysis↗

Investigating Low-Altitude Constellations of Ad-Hoc Lunar PNT System for Distributed Spacecraft Autonomy

In this study, we examine a low-altitude Lunar Position, Navigation, and Timing (LPNT) constellations and the localization performance of Centralized Extended Kalman Filter (CEKF) and Decentralized Extended Kalman Filter (DEKF) algorithms. The primary investigation involves a 100-node swarm operating at a 100 km altitude, in contrast to previous studies that examined a 21-node asset in a frozen-orbit at 5,500 km. The autonomous operation of large-scale swarm is based on two-way Inter-Satellite Link (ISL) measurements, which involve pseudoranges and relative velocities among swarm nodes. We perform a numerical assessment of the two filtering approaches, utilizing ‘fully sampled’ measurements from all available assets as well as ‘two ISL’ measurements where each spacecraft is restricted to only two antennas. This research includes an analysis of CEKF under 2-ISL constraints and evaluates the performance of DEKF in a 100-node swarm, which has not been explored in previous studies. In addition, we examine the impact of increasing the sampling frequency for DEKF, showing that the update cycle can be shortened from a 10-minute interval. A novel approach for ‘2-ISL limited’ DEKF will also be introduced, using a matching formulation that exhaustively enumerates all potential matches. This study provides valuable insights into large-scale distributed swarm operations, considering various filter configurations, sampling frequencies, matching strategies, and scalability of CEKF and DEKF for low-altitude LPNT applications. The Lunar PNT technology plays a key role in providing reliable and robust navigation services on the Moon's surface and the South pole, where the primary Lunar missions are planned. To support upcoming Lunar missions, including small satellites from NASA's Commercial Lunar Payload Services program, the Lunar PNT system must be adaptable to smaller platforms like CubeSats. Driven by the growing involvement of public and private exploration partnerships, the traditional low Earth orbit missions are shifting to beyond geosynchronous orbit [1]. These upcoming missions aim to foster a sustainable and innovative exploration program, in collaboration with commercial and international partners, to facilitate human expansion throughout the solar system and return new knowledge and opportunities to Earth [2]. As part of this trend, there are increasing efforts to utilize science missions in Lunar orbit to develop a non-dedicated and ad-hoc PNT network system. Two traditional approaches, the Deep Space Network (DSN) and the weak signal Global Positioning System (GPS), are established deep-space navigation technologies for missions beyond the geosynchronous orbit. Beginning in 1958, the DSN was developed to communicate with the Explorer 1 spacecraft based on the use of radiometric tracking in spacecraft navigation [3]. The DSN is capable of providing nearly unfettered coverage to spacecraft beyond low-Earth orbit (LEO), however, increased space mission volume has created concerns about future expectations of DSN usage for spacecraft navigation [4]. For cislunar mission applications, the position accuracy using DSN achieves 100 m (3σ) with at least three geometrically diverse ground stations when using radiometric tracking alone [5]. The DSN's dependence on Earth-based ground stations restricts its operational capabilities to periods of Earth visibility. This limitation, coupled with its poor localization performance, renders the DSN unsuitable for future lunar missions that demand continuous tracking and precise positioning. To satisfy the increasing requirements of DSN in Lunar applications, spacecrafts are also required to improve their onboard antenna power and efficiency of the transmission. However, there is an important aggregate cost trade between adding capabilities to every spacecraft and adding to a capacity on the ground that serves multiple spacecraft [6]. A weak GPS system can provide PNT service while the user spacecraft is bound to the Moon, leveraging a single, steerable high gain antenna with the relatively narrow beam which includes all the sources in its field of view [7]. However, the higher the altitude the receiver is above the GPS constellations, the poorer and the weaker are the relative geometry and the received signal powers, respectively, leading to a significant navigation accuracy reduction [8]. The transmitted power becomes weaker with increasing distance from the Earth as well as signals tracked from one of the side lobes of the GPS antenna pattern. As a results, the number of visible satellites and relative geometric condition of the GPS satellites at very high altitude drops dramatically and reduces the navigation solution accuracy. Therefore, the weak GPS system is also not an ideal way to provide PNT service to upcoming Lunar missions when considering its limited geometric condition and the recued navigation accuracy. Another navigation approach on the Moon is being developed, similar to the Global Navigation Satellite System (GNSS) on Earth, aiming to offer navigation service with continuous 24/7 coverage across the entire Lunar surface. For example, lunar communications relay and navigation systems (LCRNS) by NASA and Lunar navigation satellite systems (LNSS) by JAXA are designed to serve as dedicated Position, Navigation, and Timing (PNT) systems for the Moon. However, designing a dedicated LNSS and PNT service involves additional challenges, which are unique to the lunar environment, including limited payload capacity for the CubeSat platform, i.e., the size, weight, and power (SWaP) of the onboard clock, limited lunar ground monitoring stations, and limited financial investment as compared to the legacy Earth-GPS [9]. NASA’s focus on utilizing CubeSat platforms on the Moon leads to an alternative Lunar navigation platform that leverages the existing Lunar science and exploration assets. The small satellites used in Lunar missions can be used to create a low-cost, autonomous, ad-hoc, and on-demand mission-centric Lunar PNT swarm capable of providing PNT services to these low-cost lunar missions [10]. As upcoming Lunar missions will often operate at low-altitude about 30 km to 100 km for scientific observations and mapping purposes, the low-altitude orbital constellations could be employed to create an ad-hoc Lunar PNT system. However, several issues must be addressed, such as the instability of these orbits, which often require maintenance or are only suitable for short-duration missions, operating for fewer than 90 days. Additionally, at an altitude of 100 km, the satellites have a limited period during which they are above the horizon and capable of providing PNT service to users. The implementation of a non-dedicated, ad-hoc Lunar navigation constellation facilitates on-demand PNT services. A preliminary study of ad-hoc Lunar PNT system was conducted using 21 spacecraft in 5,5000 km altitude frozen orbits to test its feasibility and a basic performance of orbital asset localization among ad-hoc Lunar constellations in small satellites format [10]. These swarm assets are designed for autonomous localization with minimal Earth interaction, reducing dependency on bandwidth and ground resources. The design in [10] demonstrated the feasibility of a decentralized PNT approach, specifically employing a DEKF approach for state estimation, which helps minimize onboard operating costs. The DEKF method distributes computation across individual satellites, which lightens the computational load while maintaining accuracy in orbit ephemeris and clock offsets, similar to centralized systems [11]. In a follow-on study [12], each spacecraft was limited to 2 communications antennae, forcing the selection of measurements and scheduling spacecraft activities to perform the measurements. A matching algorithm is implemented to select the best measurements and schedule position estimation updates. The decentralized localization performance is also investigated with increasing levels of network degradation for swarm assets considering the impact of intermittent and permanent communication failure, to demonstrate the robustness and fidelity of the decentralized Lunar PNT service [13]. This study confirmed that the ad-hoc PNT constellations in frozen orbit are highly robust and resilient to communication failures. However, unlike frozen orbit swarm assets, the low-altitude satellites have a limited ground view at an altitude of 100 km, where the ad-hoc Lunar constellation consists of 98 low-altitude satellites, evenly distributed across seven circular polar orbital planes, alongside two satellites in a frozen orbit at an altitude of 5,500 km (Figure 1). Therefore, the number of satellites visible to ground users is significantly limited in low-altitude orbit constellations. As each visibility of a spacecraft remains intact for only a few ticks before it moves out of the field of view, the ground user encounters challenges in maintaining continuous navigation service, resulting in sparse availability and provision of Lunar PNT system. Consequently, service availability is primarily restricted to the Lunar South Pole region (Figure 2). Given these limitations and concerns, the localization performance of low-altitude swarm assets will be assessed in this study. We focus on the investigation of the localization performance of low-altitude swarm assets and ground users near the Lunar South Pole. The overall flow of the Lunar PNT simulation incorporates the DEKF approach of asset localization and the weighted least-squares approach in user localization (Figure 3). The autonomous Lunar PNT simulation is primarily implemented in MATLAB, where the DEKF based on the matching scheduler is implemented with Google’s OR-tools as a model builder and Gurobi optimization tool as a backend solver. The General Mission Analysis Tool (GMAT) is utilized to generate ephemeris data for swarm assets, and accounts for satellite orbital details, mass, and perturbations like solar radiation pressure and drag coefficients. Each ephemeris dataset is produced in the Moon International Celestial Reference Frame (ICRF) inertial coordinate system. For state estimation, the distributed swarm assets rely on two-way Inter-Satellite Link (ISL) measurements, which involve tracking pseudoranges and relative velocities between visible satellites and anchor nodes during each observation. Numerical evaluations of the decentralized localization process are conducted to demonstrate the feasibility of the low-altitude PNT system in providing reliable navigation services. The main approach involves using DEKF and CEKF to localize 100 satellites in low-altitude constellations, where the CEKF is implemented to serve as a baseline for comparing the performance of distributed algorithms. In both cases, we evaluate ‘fully sampled’ measurements from all available assets, and ‘two ISL’ measurements when spacecraft are constrained to have only two antennas. We test four estimation techniques: CEKF fully sampled, CEKF two ISL, DEKF fully sampled, and DEKF two ISL filters. As the DEKF update cycle is comprised of network setup, communication, and computations, a global broadcast network and 2-way ISL network setup will take from 4 to 6 minutes as maximum [12]. In this simulation, the DEKF update cycle is set to 10 minutes, including a 4-minute latency for obtaining and computing the actual measurement updates. We experiment an increased update cycle to demonstrate the feasibility and evaluate the impact on localization performance using various tuning values for measurement noise covariances (Figures 4 and 5). By comparing centralized and decentralized approaches using a matching algorithm, we analyze the influence of cross-correlation factors in the covariance matrix, assuming 100% reliability of all assets and measurements. The increased frequency and the adjustments of tuning parameters reveal distinct error patterns between the two scenarios. The localization accuracy of the swarm assets and ground users is assessed by taking the median error across 100 assets and one ground user (84.9°S, 137.5°E) over 7-day simulation period (Table 1). Since the user localization accuracy is significantly affected by the performance of the swarm assets, it is crucial to maintain high localization accuracy within the swarm. This study will continue to explore decentralized filtering for autonomous LPNT operations, with further investigation of an 'iterative' matching approach which enumerates every valid matching pair, planned for the following month.

Yeji Kim↗

Proton Scattering in NOvA Test Beam

The NuMI Off-axis $\nu_e$ Appearance (NOvA) experiment is a two-detector, long baseline, neutrino oscillation experiment, which aims to make a determination of the neutrino mass ordering, the octant of $\theta_{23}$, and measure possible charge-parity (CP) violation. Determining these parameters is a step towards answering some of the fundamental questions about neutrinos. Are neutrinos their own antiparticle? Could neutrinos be responsible for the matter-antimatter asymmetry of the universe? How do neutrinos get their mass? Answering these questions requires precise measurements of the parameters that govern how neutrinos behave, namely the mass squared splittings, mixing angles, and possible CP-violating phase factor in the PNMS matrix. Reaching high precision requires minimizing both statistical and systematic errors. As NOvA continues to accumulate data, the sizes of our statistical errors continue to shrink such that in a couple of years they will be comparable to our systematic errors. NOvA's current systematic errors are dominated by energy and calibration uncertainties. The Test Beam program was initiated to address these uncertainties by assessing the detector response in an environment where more parameters of the incoming particles are known. One of the main goals of the Test Beam program is to garner a better understanding of our systematic errors, particularly in energy measurement and detector response, so that the modeling of these parameters can be improved. I studied protons in the Test Beam Detector with momenta around \unit[1]{GeV/c}, the high end of the momentum range relevant to quasielastic neutrino interactions in NOvA. One question this sample can help to answer is if we correctly simulate the fraction of protons that have an inelastic scattering interaction versus ranging out. If protons inelastically scatter at different rates in our simulation and data, the correction applied to the hadronic system when reconstructing neutrino energy could be inaccurate, since events with inelastic scattering will have less visible energy in the detector. This is particularly important for quasielastic $\nu_\mu$ events, which have the lowest hadronic energy resolution, and therefore give us our best measurements of $\Delta m_{32}^2$ and $\theta_{23}$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Sparse Approximate Multifrontal Factorization with Butterfly Compression for High-Frequency Wave Equations

In this work, we present a fast and approximate multifrontal solver for large-scale sparse linear systems arising from finite-difference, finite-volume or finite-element discretization of high-frequency wave equations. The proposed solver leverages the butterfly algorithm and its hierarchical matrix extension for compressing and factorizing large frontal matrices via graph-distance guided entry evaluation or randomized matrix-vector multiplication-based schemes. Complexity analysis and numerical experiments demonstrate $\mathcal{O}(N\log^2 N)$ computation and $\mathcal{O}(N)$ memory complexity when applied to an $N\times N$ sparse system arising from 3D high-frequency Helmholtz and Maxwell problems.

97 MATHEMATICS AND COMPUTING↗

Component modes damping assignment methodology for articulated, multiflexible body structures

To simulate the dynamical motion of articulated, multiflexible body structures, one can use multibody simulation packages such as DISCOS. To this end, one must supply appropriate reduced-order models for all of the flexible components involved. The component modes projection and assembly model reduction (COMPARE) methodology is one way to construct these reduced-order component models, which when reassembled capture important system input-to-output mapping of the full-order model at multiple system configurations of interest. In conjunction, we must also supply component damping matrices which when reassembled generate a system damping matrix that has certain desirable properties. The problem of determining the damping factors of components' modes to achieve a given system damping matrix is addressed here. To this end, we must establish from first principles a matrix-algebraic relation between the system's modal damping matrix and the components' modal damping matrices. An unconstrained/constrained optimization problem can then be formulated to determine the component modes' damping factors that best satisfy that matrix-algebraic relation. The effectiveness of the developed methodology, called ModeDamp, has been successfully demonstrated on a high-order, finite element model of the Galileo spacecraft.

Lee, Allan Y.↗

Stability analysis of bridged cracks in brittle matrix composites

The bridging of matrix cracks by fibers is an important toughening mechanism in fiber reinforced brittle matrix composites. This paper presents the results of a nonlinear finite element analysis of the Mode-I propagation of a bridged matrix crack in a finite size specimen. The composite is modeled as an orthotropic continuum and the bridging due to the fibers is modeled as a distribution of tractions which resist crack opening. A critical stress intensity factor criterion is employed for matrix crack propagation while a critical crack opening condition is used for fiber failure. The structural response of the specimen (load-deflection curves) as well as the stress intensity factor of the propagating crack are calculated for various constituent properties and specimen configurations for both tensile and bending loading. By controlling the length of the bridged crack results are obtained which highlight the transition from stable to unstable behavior of the propagating crack.

Ballarini, Roberto↗