Approximation Of Stress Field Research Articles

Photoelastic study of a double edge notched plate for determination of the Williams series expansion

In this work, digital photoelasticity method is applied for assessment of the crack tip linear fracture mechanics parameters for a plate with double edge notches and different other crack configurations. The overarching objective of the study is to obtain the coefficients of the Williams series expansion for the stress and displacement fields in the vicinity of the crack tip by the digital photoelasticity technique for the double edge notched plate. The digital image processing tool for experimental data obtained from the photoelasticity experiments is developed and utilized. The digital image processing tool is based on the Ramesh approach but allows us to scan the image in any direction and to analyse the image after any number of logical operations. In the digital image processing isochromatic fringe analysis, the optical data contained in the transmission photoelastic isochromatics were converted into text file and then the points of isochromatic fringes with minimum light intensity were used for evaluating fracture mechanics parameters. The multi-parameter stress field approximation is used. The mixed mode fracture parameters, especially stress intensity factors (SIF) are estimated for specimen configurations like double edge notches and inclined center crack using the proposed algorithm based on the classical over-deterministic method. The effects of higher-order terms in the Williams expansion were analysed for different cracked specimens. It is shown that the higher order terms are needed for accurate characterization of the stress field in the vicinity of the crack tip. The experimental SIF values estimated using the proposed method are compared with analytical / finite element analysis (FEA) results, and are found to be in good agreement.

Photoelastic evaluation of stress fields and coefficients of multi-parameter asymptotic expansion of the crack-tip stress field

In this work, digital photoelasticity method is appliedfor estimation of the crack tip fracture mechanics parameters for a plate with double edge notches and different crack configurations.The digital image processing tool for experimental data obtained from the photoelasticity experiments is developedand utilized. The digital image processing tool is based on the Ramesh approach but allows us to scan the image in any direction and to analyse the image after any number of logical operations. In the digital image processing fringe analysis, the optical data contained in the transmission photoelasticisochromatics were converted into text file and then the points of isochromatic fringes with minimum light intensity were used for evaluating fracture mechanics parameters. The multi-parameter stress field approximation is used. The mixed mode fracture parameters, especially stress intensity factors (SIF) are estimated for specimen configurations like double edge notches and inclined center crack using the proposed algorithm based on the classical over-deterministic method. The effects of higher-order terms in the Williams expansion were analysed for different cracked specimens. It is shown that the higher order terms are needed for accurate characterization of the stress field in the vicinity of the crack tip. The experimental SIF values estimated using the proposed method are compared with analytical/finite element analysis (FEA) results, and are found to be in good agreement.

A Stress Approach Model of Functionally Graded Shells

This paper presents a model called SAM-FG (stress approach model of functionally graded shells) for linear elastic, thin, and moderately thick shells made of functionally graded materials. The model is an extension of the SAM-H model, originally created for homogeneous shells. Assuming that the material is orthotropic and that one of its orthotropic directions is the thickness direction, the extension consists in considering that the 3D compliance tensor may depend on the thickness coordinate. The model starts with a tunable polynomial approximation of the 3D stress field that contains the same generalized forces as SAM-H. This stress approximation verifies the 3D equilibrium equations and the stress boundary conditions at the faces of the shell. As in SAM-H, 5 generalized displacements appear in SAM-FG. By applying the Hellinger–Reissner functional and Reissner’s variational method, the generalized forces, strains, and equations in SAM-FG turn out to be the same as in SAM-H, except for the generalized constitutive equations. To prove the accuracy of the model, SAM-FG is first applied to a simply supported, functionally graded plate and its results are compared to other models. To validate the model for shell-like structures, SAM-FG results are compared to those obtained with solid finite element calculations for three case studies of structures subjected to an internal pressure. The first one deals with a hollow sphere made of an isotropic functionally graded material. The second case considers a hollow cylinder made of an orthotropic functionally graded material. In the last case, a catenoid with an isotropic functionally graded material is studied. In all cases, the mean displacements are correctly predicted, even if the main purpose of the SAM-FG model is not to calculate these fields accurately. The stress field approximations are very accurate, and since the implementation of the shell model in a finite element code would imply 5 degrees of freedom per node, SAM-FG is a good alternative to solid finite element calculations for the structural analysis of functionally graded shells with a reasonable computational cost.

Convolutional Neural Network Algorithm for Stress Field Approximation for a Uniaxially Loaded Plate with Holes

ABSTRACT Knowing the high stress concentration and its location around a hole is of practical importance in designing engineering structures, effective parameters can be selected in order to achieve minimum stress concentration around the hole(s). The visualisation of the stress distribution is done experimentally or numerically using finite element analysis (FEA). This study involved the approximation of the von Mises stress field distribution around regular holes in finite metallic plates, based on deep learning 2D encoder-decoder convolutional neural network (CNN). The CNN based stress field prediction is achieved by nonlinear 2D image training data obtained from ABAQUS, the geometric parameter of flat plate and the hole(s) are the input to CNN and the von Mises stress filed distribution is the output, assuming the plate considered is finite, isotropic, plane stress state, linearly elastic and uniaxial loading condition with total number of trainable parameters for the neural network is 3,437,143. The prediction accuracy achieved is 94.18% with a loss of 0.03% between the ground truth and the predicted images Considerable care was taken to minimize the complexity of the CNN architecture and to make it interpretable and to achieve high accuracy with available low data. The proposed method can predict acceptable accurate solution to a problem with the geometries not included in the training dataset. This approach can be directly applied as it provides immediate feedback for real-time design iterations at the early stage of design. Keywords: Convolutional neural network, finite element analysis, isotropic plate, machine learning, regular hole, stress concentration, uniaxial loading

Propagation conditions of an eccentric crack in a semi-circular disk loaded in I+II mixed mode

Multi-parameter fracture mechanics concept is applied to investigate crack propagation conditions during a three-point bending (3PB) test performed on an asymmetrical semi-circular disk under I+II mixedmode loading. The asymmetry is represented by an eccentricity of the angled crack as well as by different distances of both supports during the 3PB test from the middle of the specimen. The Williams series expansion is used for the crack-tip stress field approximation and fracture response of the crack is assessed by means of the combination of the numerical and analytical solution. Comparison between the individual configurations is carried out and discussion on the importance of the higher-order terms considered in the Williams expansion for various cases of asymmetry is presented.

Utilization of Williams’ Power Series for Estimation of Crack Behavior under Mixed-Mode Loading

Multi-parameter fracture mechanics concept has been applied to investigate crack behavior under mixed-mode loading, particularly in a semi-circular bending disc. The so-called Williams’ series expansion is used for the crack-tip stress field approximation. It has been shown that application of the generalized fracture mechanics concept can be crucial for materials with specific fracture behavior, such as elastic-plastic or quasi-brittle one, when fracture occurs not only in the very vicinity of the crack tip, but also in a more distant surrounding. Then, considering the higher-order terms of the Williams’ expansion in fracture criteria can be helpful. The attention is devoted to the analysis of the influence of various distances between the supports during the three-point bending test on the results of the further crack propagation direction.

Crack propagation in mixed-mode specimens described via multi-parameter fracture mechanics

The main objective of the paper is to employ a multi-parameter fracture mechanics concept to describe crack propagation through a specimen loaded in mixed-mode. This concept in particular was used because it has been shown that application of the generalized fracture mechanics concept can play a key role for materials with specific fracture behaviour, i.e. when fracture processes occur not only in the very vicinity of the crack tip, but also at larger distances from it. Two mixed-mode (I+II) geometries for the investigation of crack behaviour are presented here. The Williams series expansion is used for crack-tip stress field approximation. Then, considering the higher-order terms of the Williams expansion with regard to maximum tangential stress criterion can provide better estimates of the crack deflection angle. The coefficients of the Williams expansion were determined by means of the over-deterministic method for the purposes of this work. This analysis was performed for each cracked configuration, which is very time-consuming and makes the analysis very extensive. The crack propagation angle obtained by means of the generalized fracture criterion is discussed in detail. It was found that single-parameter fracture mechanics is sufficient when applied close to the crack tip and when mode I of loading prevails, while multi-parameter fracture mechanics can be recommended at larger distances from the crack tip and for configurations where mode II becomes dominant.

Multi-parameter fracture mechanics

A mixed-mode geometry has been chosen to investigate a crack propagation using the multi-parameter fracture mechanics concept. The so-called Williams’ series expansion is used for the crack-tip stress field approximation. It has been shown that application of the generalized fracture mechanics concept can be crucial for materials with specific fracture behaviour, such as elastic-plastic or quasi-brittle one, when fracture occurs not only in the very vicinity of the crack tip, but also in a more distant surrounding. Then, considering the higher-order terms of the Williams’ expansion in fracture criteria (describing the crack stability and/or crack propagation direction) can bring more precise results. The coefficients of the Williams’ expansion must be calculated numerically (for instance by means of the over-deterministic method in this work) for each cracked configuration, which is very time-consuming, and the analysis is very extensive even for a few basic cracked specimen configurations. On the other hand, a suitable choice of the geometrical configuration of the cracked disc enables performing experiments only on the specimens that could prove the theory about the importance of using the higher-order terms.

Williams expansion-based approximation of the stress field in an Al 2024 body with a crack from optical measurements

A study on the approximation of the stress field in the vicinity of crack tip in a compact tension specimen made from Al 2024-T351 is presented. Crack tip stress tensor components are expressed using the linear elastic fracture mechanics (LEFM) theory in this work, more precisely via its multi-parameter formulation, i.e. by Williams power series (WPS). Determination of coefficients of terms of this series is performed using a least squares-based regression technique known as over deterministic method (ODM) for which displacements data obtained experimentally via optical measurements are taken as inputs. The stress fields reconstructed based on the displacement data obtained experimentally by means of optical measurements are verified by means of the stress field approximations derived for the normalized CT specimen via hybrid elements.

Analysis of accuracy of Williams series approximation of stress field in cracked body – influence of area of interest around crack-tip on multi-parameter regression performance

A study on the accuracy of an approximation of the stress field in a cracked body is presented. Crack-tip stress tensor is expressed using the linear elastic fracture mechanics (LEFM) theory in this work, more precisely via its multi-parameter formulation, i.e. by Williams power series (WPS). Determination of coefficients of terms of this series is performed using a least squares-based regression technique known as over-deterministic method (ODM) for which results from finite element (FE) method computation are usually taken as inputs. Main attention is paid to a detailed analysis of a suitable selection of FE nodes whose results serve as the inputs to the employed method. Two different ways of FE nodal selection are compared – nodes selected from the crack tip vicinity lying at a ring of a certain radius versus nodes selected more or less uniformly from a specified part of the test specimen body. Comparison of these approaches is made with the help of procedures developed by the authors which enable both the determination of the coefficients of terms of the analytical WPS approximation of the stress field based on the FE results and the backward reconstruction of the field (again using WPS) from those determined terms’ coefficients/functions. The wedge-splitting test (WST) specimen with a crack is taken as example for the study.

Analytical Approximation of Crack-Tip Stress Field: Study on Efficient Determination of Coefficients of Higher Order Terms of Williams Power Expansion

The paper presents a study on accuracy of the multi-parameter approximation of the stress field in a cracked body using Williams power series calibrated according to results of FEM computation. Main attention is paid to a detailed analysis of suitable selection of FE nodes, whose results serve as inputs to the over-deterministic method (ODM) employed for determination of the coefficients of terms of the Williams expansion describing the crack-tip field. Two different ways of FE nodal selection are compared – nodes selected from the crack tip vicinity versus nodes selected from a specific part of the test specimen body in a specific way. Comparison is made with the usage of the authors’ own developed procedures enabling both the determination of the coefficients of the analytical approximation of stress field based on the FE results and the backward reconstruction of the field from those determined terms’ coefficients/functions. The wedge-splitting test (WST) specimen with a crack is taken as example for the study.

Multi-parameter approximation of the stress field in a cracked body in the more distant surroundings of the crack tip

The accuracy of the multi-parameter approximation of the stress fields in a cracked body is studied in the paper. This approximation, which uses the Williams power series expansion (WE), is intended to be used to estimate the extent of the nonlinear zone, which plays a significant role within tensile failure and the fatigue assessment of non-brittle materials. The characteristics of this zone could be potentially incorporated into methods of determining the true values of fracture parameters and the fatigue behaviour descriptors of materials exhibiting nonlinear failure. Considering the fact that in the case of elastic–plastic and especially quasi-brittle materials the size of this zone is substantial in comparison to specimen dimensions, it is necessary to consider a large region around the crack tip for this task. An automatic routine implemented as a Java application was created to determine the values of coefficients of the higher order terms of the WE that describe crack-tip fields. These values are calculated using the Over-Deterministic Method (ODM), which is applied to the results of the finite element (FE) analysis of an arbitrary mode I test geometry. Furthermore, another Java application developed by the authors provides an analytical reconstruction of the crack-tip stress field by means of the truncated WE, and enables detailed analysis of the crack-tip stress field approximation. The developed procedures simplify the analysis of the description of mechanical fields at a greater distance from the crack tip considerably. The presented research is focused on the optimisation of the selection of FE nodal results entering the ODM procedure used to determine the values of coefficients of the higher order terms of the WE. The aim is to improve the accuracy of the approximation.

Estimation of the crack propagation direction in a mixed-mode geometry via multi-parameter fracture criteria

The presented work introduces a numerical parametric study on the crack propagation direction under mixed-mode conditions (mode I + II). It is conducted for the geometry of an eccentric asymmetric fourpoint bending of a single edge notched beam specimen; various levels of mode-mixity are ensured by modifications in the crack length and crack eccentricity. The direction of crack propagation is estimated semianalytically using both the maximum tangential stress criterion and the strain energy density criterion (implemented as a procedure within the used finite element computational code) as well as numerically (from verification reasons). Multi-parameter fracture mechanics is employed in the presented work for precise analytical evaluation of the stress field in the cracked specimen. This theory is based on description of the stress and deformation fields in the cracked body by means of their approximation using several initial terms of the Williams power series. Recent studies show that utilization of only first term of the series, which corresponds to the stress intensity factor (SIF), the single controlling parameter for the crack initiation and propagation assessment in brittle materials, is insufficient in many crack problems. It appears also in this study that the higher-order terms of the asymptotic crack-tip field are of great relevance for the conducted analysis, similarly to a number of other fracture phenomena (near-crack-tip stress field approximation, non-linear zone extent estimation, etc.).

Multi-parameter approximation of stress field in a cracked specimen using purpose-built Java applications

A quality of multi-parameter approximation of the stress and displacement fields around a crack tip in a non-brittle material test specimen is studied in the paper. The stress field approximation using Williams power series is intended to be utilized for estimation of the nonlinear zone extent which potentially plays a role within methods for determination of true values of fracture parameters of materials exhibiting nonlinear failure. Considering the fact that in the case of elastic-plastic and especially quasi-brittle materials the size of this zone is substantial in comparison to the specimen dimensions, it is necessary to take a large region around the crack tip into account for this task. An automatic utility created to determine the values of coefficients of the higher order terms of Williams power series by usage of over-deterministic method applied on results of finite element analysis of arbitrary mode I test geometry is one of the tested procedures. The second one provides a backward reconstruction of the crack-tip stress field analytically by means of truncated Williams expansion. The developed procedures are based on the support of Java programming language and considerably simplify analyses of the mechanical fields’ description in a farther distance from the crack-tip. The presented research is focused on optimization of selection of FE nodal results for improvement of accuracy of the approximation.

単結晶Ｎｉ基超合金中のセル状変質域形成に対するマイクロメカニクス的アプローチ

A cellular microstructure may be nucleated in single crystal Ni-based superalloys when they have experienced a local damage associated with a plastic deformation during their in-service and post-service periods of gas turbine. In this work, an analytical method to investigate the morphological evolution in the cellular microstructure is proposed. The method is essentially based on the Eshelby's micromechanics theory, and it has been extended so as to be applied for material system containing high volume fraction inclusions, by employing the average stress field approximation by Mori and Tanaka. This method enables us to estimate a stable shape of precipitate in the material system, which must be influenced by many factors : e.g. precipitate volume fraction ; Young's modulus ratio and lattice misfit between matrix and precipitate ; external stress field in multi-axial state ; and heterogeneity of plastic strain between matrix and precipitate. A series of numerical calculations were summarized on the stable shape maps. The application of the method to predict the γ' rafting in superalloys during creep showed that the proposed method supplied reasonable explanations to all the experimental observations, and that the heterogeneity of plastic strain between matrix and precipitates may play a significant role in the shape stability of the precipitate. Furthermore, the method could successfully estimate the morphology in the cellular microstructure formed around the spherical indentation in CMSX-4 single crystal Ni-based superalloy. The present method will be a useful tool for failure analysis of superalloy components and for an estimation of in-service stress state.

A higher-order hybrid-mixed harmonic shell-of-revolution element

An accurate higher-order hybrid-mixed C 0 harmonic shell-of-revolution element is proposed. Unlike existing hybrid-mixed shell-of-revolution elements, the present element introduces additional nodeless degrees only for displacement field interpolation in order to enhance the solution convergence rate. In the course of developing the present element, we address some fundamental issues such as the effect of the nodeless degrees and the role of the stress field approximation consistent with the displacement field. Through several numerical examples, the use of the additional nodeless degrees is shown to yield very accurate results only with a few elements.

Fixed grid finite elements in elasticity problems

A fixed grid representation of the finite element domain is used to solve elasticity problems. An integration between CAD and FEM systems is made by using this fixed grid representation. Some considerations about FG approximation and stiffness matrix generation are presented. The element stiffness matrix for an element in the mesh is obtained as a factor of a standard element stiffness matrix. A least squares local approximation of stress field is used to calculate the stress at boundary elements. A problem with an analytical solution is used to measure the displacement and stress error. Several meshes and several geometrical configurations of the problem are used to test the reliability of the fixed grid finite element method.

Asymptotic analysis of stability for prismatic solids under axial loads

This work addresses the stability of axially loaded prismatic beams with any simply connected crosssection. The solids obey a general class of rate-independent constitutive laws, and can sustain finite strains in either compression or tension. The proposed method is based on multiple scale asymptotic analysis, and starts with the full Lagrangian formulation for the three-dimensional stability problem, where the boundary conditions are chosen to avoid the formation of boundary layers. The calculations proceed by taking the limit of the beam's slenderness parameter, ε (ε 2 ≡ area/length 2), going to zero, thus resulting in asymptotic expressions for the critical loads and modes. The analysis presents a consistent and unified treatment for both compressive (buckling) and tensile (necking) instabilities, and is carried out explicitly up to o( ε 4) in each case. The present method circumvents the standard structural mechanics approach for the stability problem of beams which requires the choice of displacement and stress field approximations in order to construct a nonlinear beam theory. Moreover, this work provides a consistent way to calculate the effect of the beam's slenderness on the critical load and mode to any order of accuracy required. In contrast, engineering theories give accurately the lowest order terms ( O( ε 2)—Euler load—in compression or O(1)—maximum load—in tension) but give only approximately the next higher order terms, with the exception of simple section geometries where exact stability results are available. The proposed method is used to calculate the critical loads and eigenmodes for bars of several different cross-sections (circular, square, cruciform and L-shaped). Elastic beams are considered in compression and elastoplastic beams are considered in tension. The O( ε 2) and O( ε 4) asymptotic results are compared to the exact finite element calculations for the corresponding three-dimensional prismatic solids. The O( ε 4) results give significant improvement over the O( ε 2) results, even for extremely stubby beams, and in particular for the case of cross-sections with commensurate dimensions.

Acoustoelasticity: Ultrasonic stress field reconstruction

Based on the theory of acoustoelasticity, a new ultrasonic stress reconstruction method—the generalized acoustic ratio (GAR) technique—is developed for locally plane structures and orthotropic materials. For given transit times of the three wave modes and the shear wave polarization angle, the local plane stress tensor is uniquely determined. The GAR technique yields accurate stress estimates with relatively small temperature sensitivity. Based on calibration constants from three uniaxial specimens, the entire stress field in a compact tension specimen is reconstructed. The results are in very good agreement with stress predictions from an elastoplastic finite element analysis. To further improve the measurements, a numerical technique, the stress field approximation (SFA) technique, is developed. The SFA technique uses a smooth local bicubic spline approximation and aims at improving the overall stress field estimate by enforcing the equilibrium equations, the stress boundary conditions and symmetry conditions. Numerical results show that both the average error and its spread are indeed reduced.

Formulation of the variational boundary element method with improved accuracy

A formulation which uses Lagrange multipliers to “smooth out” discontinuities in the stress field approximation at the joints of the boundary elements, and which is in fact a hybrid formulation of the finite-element method [1, 2], is considered. The improvement to the a posteriori estimate approximation is analysed.