Closed-form Solution Research Articles

Configurational force method enables fracture assessment in soft materials

Configurational mechanics offers a framework for quantifying the tendency of defects to alter the material configuration. When applied to fracture mechanics, configurational forces can be used to quantify the propensity of cracks to propagate. An alternative, well-established approach involves analytical solutions for crack tip displacement fields. However, these solutions typically apply to a limited range of constitutive behaviors and oftentimes to the linear small strain regime. The ease of calculating configurational forces in a numerical Finite Element implementation, along with their applicability to soft fracture at large strains, motivates the study of their performance as a standalone fracture framework. In contrast to the majority of works that remain theoretical and numerical, our study includes a robust experimental approach to configurational forces at finite strains. We report tensile experiments on a soft elastomer with pre-cuts ante fracture initiation. In a first attempt to approach the J-integral via configurational forces, we explore the performance of the linear elastic fracture mechanics solutions on Pacman-shaped domains that reproduce the crack tip vicinity. Then, we implement the entire boundary value problem with three-dimensional simulations that replicate the empirical tensile deformation of the soft elastomer samples. Subsequently, the results are benchmarked against estimations of the J-integral obtained through a bespoke finite strain analytical crack tip solution. With the successful validation of the configurational force method at finite strains, we aim to establish a pipeline for the calculation of configurational forces in a standalone manner and circumventing the need for close-form analytical solutions.

Aeroelasticity analysis of anisogrid lattice sandwich cylindrical shells: Extended potential and piston theories

Based on the extended potential and piston aerodynamics theories, buckling and vibration analysis of sandwich cylindrical shells subjected to external airflow in both subsonic and supersonic regimes is presented. Sandwich cylindrical shells are considered to be made of same homogenous face sheets and an anisogrid composite lattice core with hexagonal cells. An equivalent single layer (ESL) theory in the framework of classical Love’s shell theory is used to formulate the global mechanical behavior of such shells. According to the extended potential and piston aerodynamic theories, aerodynamic pressure distributions are determined for steady, inviscid, irrotational compressible, and external airflow. The closed–form solutions of critical upstream pressure and natural frequency are developed for different boundary conditions using the Galerkin method. Finally, case studies are presented to compare between results of extended potential and piston theories and investigate the effects of upstream speed, geometric ratios, boundary conditions, and lattice specifications on various response quantities. Findings indicate that the differences in results of extended potential and piston theories are less than 3% in the supersonic regime with 2.3 ≤ M ∞ ≤ 2.7 for intermediate-length shells with a wide range of R . The piston theory is inadequate for analyzing aeroelasticity behavior of very long thin sandwich cylindrical shells.

Three-Dimensional Modeling and Inversion of Gravity Data Based on Topography: Urals Case Study

In this paper, the derivation of a concise closed form for the gravitational field of a polyhedron is presented. This formula forms the basis of the algorithm for calculating the gravitational field of an arbitrary shape body with high accuracy. Based on this algorithm, a method for gravity data inversion (creating density models of the Earth’s crust) has been developed. The algorithm can accept either regular or irregular polyhedron discretization for density model creation. The models are approximated with dense irregular grids, elements of which are polyhedrons. When performing gravity data inversion, we face three problems: topography with large amplitude, the sphericity of the planet, and a long computation time because of the large amount of data. In our previous works, we have already considered those problems separately but without explaining the details of the computation of the closed-form solution for a polyhedron. In this paper, we present for the first time a performance-effective numerical method for the inversion of gravity data based on topography. The method is based on closed-form expression for the gravity field of a spherical density model of the Earth’s crust with the upper topography layer, and provides great accuracy and speed of calculation. There are no restrictions on the model’s geometry or gravity data grid. As a case study, a spherical density model of the Earth’s crust of the Urals is created.

Periodicity and lifetime of orbits around elongated asteroids

This paper presents dynamics and control of orbits around an asteroid that is modelled as a uniformly rotating homogeneous rectangular parallelepiped (HRP). Aside from the non-uniform gravitational field of the asteroid, we include the solar radiation pressure (SRP) as a perturbing force. A simple model based on the polyhedral method is employed for calculating the gravity around the HRP. First, circular periodic orbits are identified both in the presence and absence of solar radiation pressure. Then, families of general non-circular periodic orbits in the presence of SRP are determined and subsequently analysed for stability. Next, for a general orbit in the equatorial plane, its stability against a catastrophic impact of an orbiting spacecraft onto the asteroid’s surface is checked in the presence of SRP. Orbits, where the solar radiation pressure is less/more dominant over the central body’s gravity, are identified with the aim of finding the best orbits for parking a spacecraft. Then, for long-period orbits, SRP is averaged and the time averaged evolution of the orbital elements is obtained. For a special set of initial condition, a closed form solution for the evolution of the orbital eccentricity is obtained. Using this solution, a lower limit for the life-time of a spacecraft that is left uncontrolled is derived in a closed form. Finally, active control strategies are developed for orbit-keeping ofspacecraft for any choice of orbital parameters. A linear-quadratic regulator (LQR) is designed to track a given Keplerian orbit with minimum fuel expenditure. Fuel expenditure on various orbits are calculated and compared to determine the orbits of cheapest operation. The work presented is applicable to small body (asteroids and comets) exploration missions which require online cancellation of the perturbing forces with minimal knowledge of the accurate shape, size, and composition of the small body prior to the mission.

Elastostatics of nonuniform miniaturized beams: Explicit solutions through a nonlocal transfer matrix formulation

A mathematically well-posed nonlocal model is formulated based on the variational approach and the transfer matrix method to investigate the size-dependent elastostatics of nonuniform miniaturized beams. The beams are composed of an arbitrary number of sub-beams with diverse material and geometrical properties, as well as small-scale size dependency. The model adopts a stress-driven nonlocal approach, a well-established framework in the Engineering Science community. The curvature of a sub-beam is defined through an integral convolution, considering the bending moments across all cross-sections of the sub-beam and a kernel function. The governing equations are solved and the deflections are derived in terms of some constants. The formulation uses local and interfacial transfer matrices, incorporating continuity conditions at cross-sections where sub-beams are joined, to define relations between constants in the solution of a generic sub-beam and those of the first sub-beam at the left end. The boundary conditions are then imposed to derive an explicit, closed-form solution for the deflection. The solution significantly simplifies the study of nonuniform beams with multiple sub-beams. The predictions of the model for two limiting cases, namely local nonuniform and nonlocal uniform beams, are in excellent agreement with the available literature data. The flexural behavior of nonuniform miniaturized beams, composed of two to five different sub-beams and subjected to different boundary conditions, is studied. The results are presented and discussed, emphasizing the effects of the material properties, nonlocalities, and lengths of the sub-beams on the deflection. It is demonstrated that the flexural response of nonlocal nonuniform beams is more complex than local counterparts. Unlike the local beams, dividing a nonlocal uniform beam into multiple sub-beams and then reconnecting them changes the overall stiffness of the beam. The study highlights the potential to design nonuniform miniaturized beams with specific configurations to control their flexural response effectively.

Approximate analytical solution of fractional-order generalized integro-differential equations via fractional derivative of shifted Vieta-Lucas polynomial

In this paper, we extend fractional-order derivative for the shifted Vieta-Lucas polynomial to generalized-fractional integro-differential equations involving non-local boundary conditions using Galerkin method as transformation technique and obtained N - \delta + 1 system of linear algebraic equations with \lambda i, i = 0, . . . , N unknowns, together with \delta non-local boundary conditions, we obtained (N + 1)- linear equations. The accuracy and effectiveness of the scheme was tested on some selected problems from the literature. Judging from the table of results and figures, we observed that the approximate solution corresponding to the problem that has exact solution in polynomial form gives a closed form solution while problem with non-polynomial exact solution gives better accuracy compared to the existing results.

Multiple solutions and conserved vectors of a shallow water wave equation arising in fluid mechanics; Lie group analysis

Investigation has shown that shallow water wave model equations constitute a class of hyperbolic partial differential equations describing the movement below a pressure surface in a fluid, and this makes them highly useful to be studied. Therefore, in this paper we analytically investigate a one-dimensional shallow water wave equation that is highly applicable in fluid dynamics. We leverage the robust nature of the Lie group technique to provide symmetries to the equation. Commutator and adjoint representations of symmetries were generated and presented in tabular form. Sequel to that, we perform reductions of the shallow water wave model under study via its point symmetries so that its exact solutions can be found. Thus, to achieve this, various combinations of the symmetries are made. Besides, we examine the travelling wave solutions of the model. The direct integration method is used to obtain elliptic integral functions as well as elliptic solutions in terms of the Weierstrass function. In addition, we utilize standard techniques such as (G′/G)-expansion, power series, and Kudryashov’s approaches to generate more closed-form solutions to the model. Other interesting solutions achieved include hyperbolic, trigonometric, rational, and exponential functions. In search of the physical meaning associated with the results, graphical demonstrations are included, and these are given in 2D, 3D, and density plots. Thereafter, the construction of conserved vectors for the aforementioned model is done with the use of the standard multiplier approach and Ibragimov’s theorem. These techniques give the opportunity to furnish various forms of conserved currents.

Hedging lookback-barrier option by Malliavin calculus in a mixed fractional Brownian motion environment

We present a Malliavin calculus approach to a mixed fractional Brownian motion option hedging model, that adequately describes, e.g., the hedging of a lookback-barrier option with the floating strike price. The Markovian setup and smooth stochastic differentials are necessary components in the payoff function for classical Δ-hedging of a contingent claim. This is in contrast to the Malliavin calculus approach, which may be used to any type of path-dependent options. Based on the fractional Clark–Ocone formula, an entirely probabilistic computation is developed to obtain the closed-form solution of the hedging portfolio for a lookback-barrier option with the floating strike price. Numerical experiments are performed to demonstrate the performance of our proposed hedging model, and we also conduct sensitivity analysis to investigate the correlation between model parameters and the hedging portfolio.

Amplifying Sine Unit: An Oscillatory Activation Function for Deep Neural Networks to Recover Nonlinear Oscillations Efficiently

Many industrial and real-life problems exhibit highly nonlinear periodic behaviors and the conventional methods may fall short of finding their analytical or closed-form solutions. Such problems demand cutting-edge computational tools with increased functionality and reduced cost. Recently, deep neural networks have gained massive research interest due to their ability to handle large data and universality to learn complex functions. In this work, we put forward a methodology based on deep neural networks with responsive layers structure to deal nonlinear oscillations in microelectromechanical systems. We incorporated some oscillatory and non-oscillatory activation functions such as growing cosine unit known as GCU, Sine, Mish and Tanh in our designed network to have a comprehensive analysis on their performance for highly nonlinear problems. Integrating oscillatory activation functions with deep neural networks definitely outperform in predicting the periodic patterns of underlying systems. To support oscillatory actuation for nonlinear systems, we have proposed a novel oscillatory activation function called Amplifying Sine Unit denoted as ASU which is more efficient than GCU for complex vibratory systems such as microelectromechanical systems. Experimental results show that the designed network with our proposed activation function ASU is more reliable and robust to handle the challenges posed by oscillations

Generalized variational framework with minimax optimization for parametric blind deconvolution

Blind deconvolution (BD), which aims to separate unknown convolved signals, is a fundamental problem in signal processing. Due to the ill-posedness and underdetermination of the convolution system, it is a challenging nonlinear inverse problem. This paper is devoted to the algorithmic studies of parametric BD, which is typically applied to recover images from ad hoc optical modalities. We propose a generalized variational framework for parametric BD with various priors and potential functions. By using the conjugate theory in convex analysis, the framework can be cast into a nonlinear saddle point problem. We employ the recent advances in minimax optimization to solve the parametric BD by the nonlinear primal-dual hybrid gradient method, with all subproblems admitting closed-form solutions. Numerical simulations on synthetic and real datasets demonstrate the compelling performance of the minimax optimization approach for solving parametric BD.

Optimal ordering based on shortage scenarios under presence of a low-quality substitute and supply uncertainty

Global trade has shifted from short-distance supply chains within countries to long-distance supply chains between countries. As a result, delivery lead time has become a key operational factor, and uncertainty causes stock-outs, resulting in massive losses. An important solution for shortages is product substitutes, but they can affect the quality of the end products. Moreover, even when demand is expected and planned for, unexpected last-minute demand makes operational conditions challenging. This paper considers a single-period problem of replenishing a high-quality product under both lead time and demand-related uncertainty. Stock-outs during and at the end of the replenishment period can be handled only by providing a low-quality substitute product that requires repeated replacements, therefore affecting pending demand. We provide detailed modeling of possible stock-out scenarios and show that, in presence of these two types of uncertainty, the optimal order quantity of a high-quality product can be determined using an upper and lower threshold and a base stock level equal to anticipated (planned) demand over the replenishment period. In particular, we derive a closed-form solution for the optimal order quantity under general probability distributions of lead time and demand and explicit solutions for uniform distributions. We also show the conditions under which optimal order quantities can be insensitive to the cost of the high-quality product.

Memory effect on the thermoelastic responses of a porous half-space with a novel fractional heat conduction model

Studies have shown that the existence of the anomalous diffusion of porous media leads to the abnormal heat conduction, which can be well described by fractional derivatives. In the present work, memory dependent derivative is introduced into the existing linear thermoelastic theory to analyze the transient thermoelastic responses of a porous half-space subjected to a thermal loading at the boundary. By using the method of Laplace transform and its numerical inversion, the closed form solutions of temperature, volume fraction, displacement, and stress are obtained. The effects of time delay, kernel function and a memory-dependent parameter on the physical fields are discussed. The results show that the memory dependent derivative plays an essential role in controlling the heat transfer process of porous materials.

Solving Nonlinear Fractional Differential Equations by Using Shehu Transform and Adomian Polynomials

The current article provides a detailed analysis of the solution of non-linear ordinary differential equations of fractional and non-fractional order in series forms using the Shehu transform (ST) and the Adomian decomposition method (ADM), also known as the Shehu transform Adomian decomposition method (STADM). Previously, these methods were used to solve differential equations of integer order as well as a very small number of ordinary and fractional differential equations (FDEs). The Caputo's operator is used in a number of well-known FDEs, including the logistic equation, the Van der Pole equation, and other non-fractional order differential equations like the nonlinear Bratu type equation. It is noted that all of the example issues had series solutions thanks to STADM. Plotting the graph for several series terms of the series solution demonstrates how the approximate solution tends to the closed form solution. In some example problems the impact of α is shown.

A closed-form structural stress-strain relation and implications on plastic design limits in codes and standards

The deformation of elastic-perfectly plastic beam or plate under combined through-thickness membrane and bending stresses (i.e., structural stress) is analyzed in this paper. Closed-form solutions of “plane-remains-plane” strain (i.e., structural strain) distribution as a function of applied stress and yield strength is derived. The relationship between structural stress and structural strain is also revealed for gaining insights on not only the plastic collapse conditions, but also through-thickness deformation behaviors leading to limit state. Based on the closed-form strain solutions, four deformation modes are identified: (1) elastic, (2) one-side yielding, (3) two-side yielding, and (4) plastic collapse. A membrane-bending stress interaction diagram is constructed to depict four zones that are associated with the four deformation modes in the normalized membrane and bending stress space. Closed-form solutions of residual stress, residual strain, and structural stress range-strain range relation after a loading-unloading cycle are also derived. Based on the parabolic plastic collapse stress limit curve, parabolic design stress limit curves with coherent safety margins for any combinations of membrane and bending stresses are constructed. The new insights gained in this study offer a more consistent definition of design stress limit than those used in the current codes and standards, e.g., ASME BP&V Code.

Efficient technique for evaluation of three-dimensional elastic–plastic fracture mechanics parameters based on equivalent distributed stress concept

This study proposes an efficient technique for the evaluation of the elastic–plastic fracture mechanics (EPFM) parameters of three-dimensional (3D) cracks based on the modified crack cohesive force (CCF) model. In the current method, the 3D crack’s EPFM parameters are analyzed by applying the pseudo crack face traction (CFT) to a two-dimensional (2D) center crack in an infinite plate. This pseudo traction is called the equivalent distributed stress (EDS). The EDS to be applied on the 2D center crack is determined in a manner that reproduces the relationship between crack length and stress intensity factor (SIF) of the intended 3D crack. Once EDS is determined, EPFM parameters for 3D cracks can be obtained instantly by using the center crack’s analytical closed-form solutions. This approach can compute the 3D crack’s EPFM parameters in dramatically less computation time than 3D finite element (FE) analysis. The effectiveness of the proposed method is validated by comparing the EPFM parameters of a 3D penny-shape crack obtained from this method against the reference analytical solutions and numerical FE analysis results. The parameters determined by the present method show excellent agreement with the reference solutions.

An Underwater Source Localization Method Using Bearing Measurements.

Angle-of-arrival (AOA) measurements are often used in underwater acoustical localization. Different from the traditional AOA model based on azimuth and elevation measurements, the AOA model studied in this paper uses bearing measurements. It is also often used in the Ultra-Short Baseline system (USBL). However, traditional acoustical localization needs additional range information. If the range information is unavailable, the closed-form solution is difficult to obtain only with bearing measurements. Thus, a localization closed-form solution using only bearing measurements is explored in this article. A pseudo-linear measurement model between the source position and the bearing measurements is derived, and considering the nonlinear relationship of the parameters, a weighted least-squares optimization equation based on multiple constraints is established. Different from the traditional two-step least-squares method, the semidefinite programming (SDP) method is designed to obtain the initial solution, and then a bias compensation method is proposed to further minimize localization errors based on the SDP result. Numerical simulations show that the performance of the proposed method can achieve Cramer-Rao lower bound (CRLB) accuracy. The field test also proves that the proposed method can locate the source position without range measurements and obtain the highest positioning accuracy.

Few-Shot Text Classification with an Efficient Prompt Tuning Method in Meta-Learning Framework

Meta-learning stands as a prevalent framework utilized in few-shot learning methods. Nonetheless, its efficacy hinges on substantial data availability during meta-training. Recent work adeptly tackled this hurdle by synergizing prompt tuning with the meta-learning paradigm, consequently attaining unparalleled performance on four benchmarks (FewRel, HuffPost, Reuters and Amazon). Nonetheless, the implementation efficacy of the previous method leaves room for enhancement, which is especially crucial when tuning larger language models. To this end, we introduce another expedited prompt tuning approach nested within the meta-learning framework. The novel approach normalizes the label information and sample information and uses the regression method to obtain the closed-form solution of each few-shot task, which significantly enhances inference speed, achieving a twofold improvement, while concurrently elevating average accuracy by [Formula: see text]% on the same benchmarks. Moreover, it demonstrates enhanced stability when faced with limited meta-training data, which is more applicable in many real scenarios where parallel data is rare. The source code is available to reproduce the results (http://github.com/Dr-Lv/EMPT).

MIMO radar partially correlated waveform design based on chirp rate diversity

Waveform design is a significant topic in multiple-input-multiple-output (MIMO) radar. In this paper, we address the problem of waveform design for colocated narrowband MIMO radar transmit beampattern synthesis. The existing methods have not comprehensively considered the transmit beampattern synthesis under multiple practical constraints, such as constant modulus constraint, Doppler tolerance, high range resolution and low auto-correlation sidelobe. We consider these practical constraints and propose a novel chirp rate diversity-based linear frequency modulated (CRD-LFM) waveform design for approximating the desired MIMO radar transmit beampattern. To tackle the problem, we improve the classical two-step approach. Firstly we design the covariance matrix using the existing closed-form solution. Then we model the problem of obtaining the waveforms from the covariance matrix as a nonlinear optimization problem, and we utilize the Sequential Quadratic Programming (SQP) algorithm to solve it. The simulations demonstrate that the waveforms obtained by this method can well approximate the desired transmit beampattern and have higher range resolution. Furthermore, we verify that the suppression of auto-correlation sidelobe can be achieved using amplitude weighting, which is easily realizable in practical engineering.

Surface-Based Probabilistic Fiber Tracking in Superficial White Matter.

The short association fibers or U-fibers travel in the superficial white matter (SWM) beneath the cortical layer. While the U-fibers play a crucial role in various brain disorders, there is a lack of effective tools to reconstruct their highly curved trajectory from diffusion MRI (dMRI). In this work, we propose a novel surface-based framework for the probabilistic tracking of fibers on the triangular mesh representation of the SWM. By deriving a closed-form solution to transform the spherical harmonics (SPHARM) coefficients of 3D fiber orientation distributions (FODs) to local coordinate systems on each triangle, we develop a novel approach to project the FODs onto the tangent space of the SWM. After that, we utilize parallel transport to realize the intrinsic propagation of streamlines on SWM following probabilistically sampled fiber directions. Our intrinsic and surface-based method eliminates the need to perform the necessary but challenging sharp turns in 3D compared with conventional volume-based tractography methods. Using data from the Human Connectome Project (HCP), we performed quantitative comparisons to demonstrate the proposed algorithm can more effectively reconstruct the U-fibers connecting the precentral and postcentral gyrus than previous methods. Quantitative validations were then performed on post-mortem MRIs to show the reconstructed U-fibers from our method more faithfully follow the SWM than volume-based tractography. Finally, we applied our algorithm to study the parietal U-fiber connectivity changes in autosomal dominant Alzheimer's disease (ADAD) patients and successfully detected significant associations between U-fiber connectivity and disease severity.

Analytical solutions for single and multiple scattering from rib-stiffened plates in water.

The interaction of an acoustic plane wave with a pair of plates connected by periodically spaced stiffeners in water is considered. The rib-stiffened structure is called a "flex-layer" because its low frequency response is dominated by bending stiffness. The quasi-static behavior is equivalent to a homogeneous layer of compressible fluid, which we identify as air for the purposes of comparison. In this way, an air layer is acoustically the same as a pair of thin elastic plates connected by a periodic spacing of ribs. At discrete higher frequencies, the flex-layer exhibits perfect acoustic transmission, the cause of which is identified as fluid-loaded plate waves propagating back and forth between the ribs. Both the low and finite frequency behavior of the flex-layer are fully explained by closed-form solutions for reflection and transmission. The analytical model is extended to two flex-layers in series, introducing new low and high frequency phenomena that are explained in terms of simple lumped parameter models.