Airy Stress Function Research Articles

Investigating teeth root stresses of various spur gears using stress function and FEM

AbstractThe gear is the most efficient and reliable portion of the transmission system. In machine tools, gears are utilized to handle heavy loads. These gears are used indefinitely if the prerequisites are met. Gear tooth bending stress for numerous spur gears has been examined in this research. The influence of changing the pressure angles on the spur gear tooth root stress has been investigated. The purpose of this study is to look for a new way of finding the stresses generated in the roots of the gear teeth as well as the effect of changing pressure angles on these stresses. Solidworks 2021 uses parametric formulation to model gears, analyzed for deformation and maximum bending stress in ANSYS. The airy stress function is used to conduct the nominal root tooth stress. The enhancement percentages of modified cases are 5.92%, 1.70%, 10.70%, 17.22%, 10.86%, 24.32%, 20.91%, and 26.15% for the cases 2, 3, 4, 5, 6, 7, 8, and 9, respectively according to the stress function findings.

The impact of periodic axial loads on nonlinear dynamic instability behavior of Inconel 625 pipes

In various engineering fields like aerospace and aircraft structures or marine and offshore platforms, constitutive material of critical components should be made of specific materials that can work properly in the required workspace. Such materials must have excellent properties such as high mechanical strength as well as great resistance to corrosion, oxidation, and creep. Inconel 625 is a superalloy that is chosen as constitutive material of critical components due to its great abilities. On the other hand, since investigating Inconel 625 pipe has not been done yet, different mechanical characteristics of using structures made of Inconel 625 should be assessed. Additionally, doing so would be necessary to gather information for current industrial affairs and also future investigations. Therefore, the nonlinear dynamic instability response of axially loaded Inconel 625 pipes is investigated in the current article. The pipe structure is modeled via the Donnell shell theory and nonlinear von Kármán theory. The motion equations of pipes are established by applying the Hamiltonian approach. Then, in order to alter the nonlinear derived partial differential equations into the Mathieu-Hill equation, both Navier's solution and Airy stress function are implemented. Additionally, the amplitudes of steady-state oscillation of the Inconel 625 pipe are determined by employing Bolotin's method. Eventually, the impacts of various effective parameters on the nonlinear dynamic behaviors of Inconel 625 pipes are evaluated. The results indicate static and dynamic load factors possess a remarkable effect on instability exciting areas and steady-state vibration amplitudes of pipe. Moreover, the dynamic instability response of the pipe is dependent on the radius-to-thickness and length-to-radius ratios, and also how the ratios are affected depends on the wave number.

The emergence of Airy stress function in two-dimensional disordered packings of particles

The packing of hard-core particles in contact with their neighbors is a statically determinate problem. The probability functional method allows deriving the full system of equations for the stress tensor components analytically. For an isotropic and homogeneous two-dimensional packing, we derive the classical Euler–Cauchy and Navier equations; the latter allows expressing the stress tensor of the system in terms of the Airy stress function.

An innovative and simple approach for predicting tunnelling-induced ground movements semi-analytically

ABSTRACT A semi-analytical elastic solution for predicting tunnelling-induced ground movements in clayey media is presented in this study. This study modelled the non-uniform ground movements (or non-uniform ground loss) around the tunnel cavity using the Gaussian distribution function. A new approach is proposed to model the tunnel peripheral displacements other than the simplistic trigonometric functions. It is shown that by modelling the ground loss in an elliptical fashion, (i) better prediction of ground surface settlements is obtained compared to the circular ground loss in some cases, while (ii) in other cases, circular ground loss gives better predictions. This indicates that both these approaches should be used to predict the ground surface settlements. Further, the ground displacements observed in East-West Metro, Kolkata, are discussed as a case study. The semi-analytically obtained displacements give good conformity with the observed field measurements in East-West Metro. This solution is much simpler to follow and can be easily employed in a code compared to the complex variable method or airy stress function method.

Airy stress function for proposed thermoelastic triangular elements

Airy stress function for proposed thermoelastic triangular elements

A developed hybrid experimental–analytical method for thermal stress analysis of a deep U-notched plate

A general framework of hybridizing experimentally-recorded load-induced thermal information by means of thermoelastic stress analysis (TSA) of a loaded structure with the Michell solution of Airy stress function in isotropic linear elasticity is proposed for in-plane stresses determination. The capability of the proposed hybrid method in separating the TSA signals into individual stresses is demonstrated by stress-analyzing a deep U-notched aluminum plate without neither knowing the entire geometry and distant loading/boundary conditions nor using supplementary experimental information. Even though no experimental data were processed at, and adjacent to, the edges of the U-notched plate, the individual stresses are determined throughout the surface of the plate. Prior processing actual experimental TSA data, the accuracy and the stability of the proposed hybrid method through adding artificial noise scatters in the processed simulated data with the number of retained terms in the generalized Airy stress expressions are firstly investigated. The effects of the superimposed noise scatters as well as the number of employed data in the evaluated stresses are also assessed. Finally, the reliability of the determined hybrid-experimental stresses is supported through a comparison with finite-element predictions and strain gauge measurements. The ease of implementation, the accuracy and the precision of the evaluated stress, and the versatility to the type of employed data make the proposed hybrid method more attractive and superior than other hybrid methods.

Numerical solution of problem of equilibrium of parallelepiped in stresses

Usually, planar and spatial problems of the theory of elasticity in stresses are solved using the Airy or Maxwell stress function, respectively. Solving a plane problem usually comes down to solving a biharmonic equation. In the spatial case, it can be solved by the Castigliano variational method. This work is devoted to the mathematical and numerical solution of the spatial problem of elasticity theory directly concerning stresses. In this case, the boundary value problem consists of three equilibrium equations and three Beltrami-Michell equations with the corresponding boundary conditions. At the same time, equilibrium equations are considered as missing boundary conditions. Symmetric finite-difference equations are constructed, and the problems of equilibrium of a parallelepiped under the action of domed and uniformly distributed loads are solved numerically. The reliability of the results obtained is substantiated by comparing the numerical results with the wellknown Filonenko-Borodich solution and with the results of solving the boundary value problem of the equilibrium of the parallelepiped with respect to displacements.

Interactive Exploration of Tension-Compression Mixed Shells

Achieving a pure-compression stress state is considered central to the form-finding of shell structures. However, the pure-compression assumption restricts the geometry of the structure's plan in that any free boundary edges cannot bulge outward. Allowing both tension and compression is essential so that overhanging leaves can stretch out toward the sky. When performing tension-compression mixed form-finding, a problem with boundary condition (BC) compatibility arises. Since the form-finding equation is hyperbolic, boundary information propagates along the asymptotic lines of the stress function. If conflicting BC data is prescribed at either end of an asymptotic line, the problem becomes ill-posed. This requires a user of a form-finding method to know the solution in advance. By contrast, pure-tension or pure-compression problems are elliptic and always give solutions under any BCs sufficient to restrain rigid motion. To solve the form-finding problem for tension-compression mixed shells, we focus on the Airy's stress function, which describes the stress field in a shell. Rather than taking the stress function as given, we instead treat both the stress function and the shell as unknowns. This doubles the solution variables, turning the problem to one that has an infinity of different solutions. By enforcing equilibrium in the shell interior and prescribing the correct matching pairs of BCs to both the stress function and the shell, a stress function and shell can be simultaneously found such that equilibrium is satisfied everywhere in the shell interior and thus automatically has compatible BCs by construction. The problem of a potentially over-constrained form-finding is thus avoided by expanding the solution space and creating an under-determined problem. By varying inputs and repeatedly searching for stress function-shell pairs that fall within the solution space, a user is allowed to interactively explore the possible forms of tension-compression mixed shells under the given plan of the shell.

Advancing contact of a 2D elastic curved beam indented by a rigid pin with clearance

A two-dimensional analytical solution is presented for stresses and displacements in an elastic curved beam forming an incomplete ring in frictionless and unbonded contact with a rigid pin loaded by a point force and in the presence of clearance. The circular beam is modeled as an incomplete elastic thick ring, constrained at both ends and in a plane stress state. The stress and displacement fields within the beam are derived from a biharmonic Airy stress function, according to the Michell solution in polar coordinates. The mixed boundary value problem is reduced to a set of dual series equations and then to a non-homogeneous linear system of infinite equations, which is then solved by truncation. The non-linear relations between the applied load and the contact angle or the pressure distribution are obtained by using an inverse method. The analytical results are compared with finite element predictions for a pin-lug connection and a reasonable agreement is observed for several typical geometries. The peaks of contact pressure and von Mises equivalent stress and their location within the curved beam are evidenced.

Nonlinear dynamic stability analysis of axial impact loaded structures via the nonlocal strain gradient theory

In engineering applications, there has been an overwhelming tendency towards portability, miniaturization, and integration in recent years. To link the intrinsic size dependency feature of the small-scale structure with its structural stability, the nonlocal strain gradient theory, which captures the size effect in a more general size-dependent continuum-based model, is introduced to explore the nonlinear dynamic stability behaviour of nanoplates. Four types of axial impact loading configurations, namely, sinusoidal, exponential, rectangular, and damping, are considered. Some practical factors, such as Winkler-Pasternak elastic foundation and damping, are taken into account in the analysis. The equations of motion for the size-dependent initially imperfect plate are derived in the framework of the first-order shear deformation plate theory in conjunction with the Von Kármán nonlinear terms. Then the Airy stress function corresponding to simply supported nanoplate is introduced; then, by applying the Galerkin method, the obtained differential equations are addressed by the fourth-order Runge-Kutta algorithm. Subsequently, the specific value of the critical dynamic buckling load is determined by the Volmir criterion. Organic solar cells (OSCs), a type of emerging solar-to-electrical energy conversion nanodevice, are used as an illustrative example within the existing framework. The effects of size dependency in conjunction with the pulse load configuration, the initial imperfection, the elastic foundation, as well as the damping ratio on the nonlinear dynamic buckling behaviour of the OSC are thoroughly investigated.

On designing plane-faced funicular gridshells

Concepts from isotropic geometry, Timoshenko’s shell theory, Airy stress functions and Maxwell’s reciprocal diagrams are combined in the design of plane-faced funicular gridshells. The notions of self-Airy and mixed-Airy gridshells are introduced, with an emphasis on self-tied gridshells. This paper extends the work of J.C. Maxwell for 2D pin-jointed trusses to 2.5D gridshells with the addition of a new reciprocal figure called the slope diagram. The form, force and slope diagrams are combined by a mixed area calculation to produce another new figure called the Maxwell-Mondrian diagram. A powerful new design process leveraging the relationship between the gridshell geometry and the Airy stress function is presented.

Calculating local geomembrane strains from gravel particle indentations with thin plate theory

A new general method is presented to calculate local strains in geomembranes from the deformed shape imposed by overlying coarse gravel particles under vertical pressure. Past methods assume that the geomembrane attains its deformed shape by only deforming vertically and hence neglect the effect of lateral displacements on strain. The new method treats the geomembrane as a thin plate with mid surface components of displacements in three directions (x, y and z). Lateral components of displacements (those in the x and y directions) are related to vertical displacements (z direction) by large-strain-displacement relationships and compatibility of strains. Normal and shear strains in the lateral directions are calculated using Airy's stress function and a linear elastic constitutive law. Bending and torsional strains calculated from curvature and are added to the mid surface strains to find strains on the top and bottom surfaces. The method was validated against data sets with known three-dimensional displacements and strain generated by finite element analysis. The application of the new method to calculate local strains in a geomembrane from the deformed shape obtained from a protection-layer-assessment physical test is illustrated.

Couple stresses and discrete potentials in the vertex model of cellular monolayers.

The vertex model is widely used to simulate the mechanical properties of confluent epithelia and other multicellular tissues. This inherently discrete framework allows a Cauchy stress to be attributed to each cell, and its symmetric component has been widely reported, at least for planar monolayers. Here, we consider the stress attributed to the neighbourhood of each tricellular junction, evaluating in particular its leading-order antisymmetric component and the associated couple stresses, which characterise the degree to which individual cells experience (and resist) in-plane bending deformations. We develop discrete potential theory for localised monolayers having disordered internal structure and use this to derive the analogues of Airy and Mindlin stress functions. These scalar potentials typically have broad-banded spectra, highlighting the contributions of small-scale defects and boundary layers to global stress patterns. An affine approximation attributes couple stresses to pressure differences between cells sharing a trijunction, but simulations indicate an additional role for non-affine deformations.

Numerical Solution of the Two-Dimensional Elasticity Problem in Strains

Usually, the boundary value problems of the theory of elasticity are formulated with respect to displacements, and are reduced to the well-known Lame equations. Strains and stresses can be calculated from displacements as a solution to Lame's equation. Also known are the Beltrami Mitchell equations, which make it possible to formulate the boundary value problem of the theory of elasticity with respect to stresses. Currently, the boundary value problems of the theory of elasticity in stresses are studied in more detail in the two-dimensional case, and usually solved numerically with the introduction of the Airy stress function. But, the direct solution of boundary value problems of elasticity theory with respect to stresses requires further researches. This work, similarly to the boundary value problem in stresses, is devoted to the formulation and numerical solution of boundary value problems of the theory of elasticity with respect to deformations. The proposed boundary value problem consists of six Beltrami-Mitchell-type equations depending on strains and three equations of the equilibrium equation expressed with respect to deformations. As boundary conditions, in addition to the usual conditions for surface forces, three additional conditions are also introduced based on the equilibrium equations. The boundary value problem is considered in detail for a rectangular area. The discrete analogue of the boundary value problem is composed by the finite difference method. The convergence of difference schemes and an iterative method for their solution are studied. Software has been developed in the C++ environment for solving boundary value problems in the theory of elasticity and deformation. A number of boundary value problems on the deformation of a rectangular plate are solved numerically under various boundary conditions. The reliability of the obtained results is substantiated by comparing the numerical results, with the exact solution, as well as with the known solutions of the plate tension problems with parabolic and uniformly distributed edge loads.

Form-finding of shells containing both tension and compression using the Airy stress function

Pure-compression shells have been the central topic in the form-finding of shells. This study investigated tension-compression mixed-type shells by utilizing a NURBS-based isogeometric form-finding approach that analyses Airy stress functions to expand the possible plan geometry. A complete set of smooth version graphic statics tools is provided to support the analyses. The method is validated using examples with known solutions, and a further example demonstrates the possible forms of shells that the proposed method permits. Additionally, a guideline to configure a proper set of boundary conditions is presented through the lens of asymptotic lines of the stress functions.

Wavelet solution for hygrothermomechanical bending of initially defected plate undergoing large deformation on nonlinear elastic foundation

An efficient wavelet approach has been developed for geometrically nonlinear bending analysis of largely deformed plate with arbitrary initial curvature on nonlinear Winkler–Pasternak foundation loaded with sinusoidally hygrothermal stresses. Translational and rotational constraints have been analyzed and equivalently transformed into the Cauchy-type boundary conditions of Airy stress function, which are efficiently approached by interpolating Coiflet-type wavelet. A newly hygrothermomechanical bending model of plate subjected to the adiabatically impermeable edges or with boundaries of sinusoidal temperature and humidity has been firstly formulated. The highly coupled and nonlinear governing equations have been decomposed into linear differential terms by homotopy transformation and successfully solved by the Coiflet-type wavelet Galerkin method. The wavelet strategy has been validated superior in dealing with extremely nonlinear bending and high-precision reconstitution of initial deflection in hygrothermal environment, while the obtained solutions are in excellent agreement with published results. Effect of initial deflection, hygrothermal load and elastic foundation has been investigated, which implies the hygrothermal effect on plate in linear bending is prominent but gradually negligible in large deformation, while the nonlinear effects of elastic foundation become significant by the increase of deflection but not obvious at smally deformed state.

Elastic-Plastic Stress Analysis of Steel Fiber Reinforced Composite Plates Under Axial Load

Composite materials, obtained by combining two or more materials; It is defined as a new type of material with high strength, high rigidity and lightness. Composite plates are structural elements that are used in machines and structures under different loads, consist of at least two types of materials and can be produced in various constructions. In this study, elastic-plastic stress analysis of polymer matrix continuous fiber reinforced composite plate under axial load was solved with Airy Stress Function proposed as a 5th order non-uniform polynomial to solve the elasticity problem. Polyethylene matrix composite reinforced with steel fibers was taken as the plate material and the material was accepted as ideal elastic-plastic. Tsai-Hill Yield Criterion was used for the plastic solution. According to the results of the analysis, as the fiber angle increased in the composite plate, the plastic stress limit decreased, the increase in the fiber angle decreased the plastic stress limit, and the decrease in the plastic stress limit caused the residual stresses to increase.

Analytic formulation of elastic field around edge dislocation adjacent to slanted free surface.

Explicit and tractable formulation of the internal stress field around edge dislocations is indispensable for considering the mechanics of fine crystalline solids, because the motion of edge dislocations in a slanted direction with respect to the free surface often plays a vital role in the plastic deformation of the solids under loading. In this study, we formulated an analytical solution for the stress distribution that occurs around edge dislocations embedded in a semi-infinite elastic medium. This formulation is based on the image force method and the Airy stress function method; it describes the variation in the stress distribution with changes in the slanted angle between the traction-free flat surface of the medium and the Burgers vector of the edge dislocation. Furthermore, our analytical solution shows that the attractive force acting on the edge dislocation due to the presence of the free surface is always perpendicular to the surface, regardless of the relative angle of the Burgers vector with the surface.

Nonlinear and chaotic vibrations of FG double curved sandwich shallow shells resting on visco-elastic nonlinear Hetenyi foundation under combined resonances

In this study, the nonlinear and chaotic instability of functionally graded (FG) double curved shallow sandwich shells resting on a viscoelastic Hetenyi foundation, under simultaneous effect of in-plane and transverse excitations is studied. Employing the third-order Reddy theory and von-Karman relations and the Hamilton principle, the partial differential equations of motion under movable SS boundary conditions are derived. Introducing the trigonometric Airy stress function and applying the Galerkin’s method, the equations are reduced to a set of nonlinear ODEs with time. Then, under forcing resonance conditions and using the perturbation method, the modulation equations at the stationary conditions are derived and solved numerically. Then stability of non-trivial solutions for the resonance amplitude corresponding to the presence of limit cycle oscillation is investigated. Then by extracting the characteristic resonance amplitude curves, the effect of various parameters, including: frequency detuning parameter, in-plane excitation amplitude, linear, shear and nonlinear stiffness and damping parameters of the foundation on nonlinear response are analyzed. Finally, by extracting two degrees of freedom system time response, bifurcation and chaotic characteristic curves of the problem, the conditions for occurrence the periodic, double periodic, multi-periodic and chaotic behaviors under the simultaneous effect of parametric and internal resonances are studied comprehensively.

Vibration and buckling analyses of functionally graded plates based on refined plate theory using airy stress function

In this paper, an approach based on the refined plate theory and Airy stress function has been proposed to investigate the vibration and buckling behaviors of functionally graded (FG) rectangular plates. The significant feature of the proposed approach is considering only two unknowns in the displacement field in which the contribution due to shear and bending to the total transverse displacement are clarified. Using the extended Hamilton’s principle and defining an Airy stress function corresponding to the compatibility equation, the equations of motion, which do not explicitly include the in-plane displacements, are derived. The accuracy and effectiveness of the current model is shown by comparing the natural frequencies and buckling loads of various FG rectangular plates calculated by the proposed approach with even three-dimensional (3-D) and quasi-3D solutions. Besides, the exact dynamic response of a square FG plate due to a harmonic central force is investigated using modal analysis. This approach is capable of handling quasi-3D models , by selecting proper functions in the displacement field to consider the thickness stretching effect. By doing this, behavior of FG square plates with various coefficient functions are compared with each other. Therefore, through implementing the Airy stress function, the number of variables is reduced, in turn, simplifying the dynamic model. Consequently, it can be used for wide range study of static and dynamic behaviors of FG plates.