Search results

The purpose of this paper is to present special nine‐node quadrilateral elements to discretize the un‐cracked boundary and the inclined surface crack in a transversely isotropic cuboid under a uniform vertical traction along its top and bottom surfaces by a three‐dimensional (3D) boundary element method (BEM) formulation. The mixed‐mode stress intensity factors (SIFs), K_I, K_II and K_III, are calculated.

A 3D dual‐BEM or single‐domain BEM is employed to solve the fracture problems in a linear anisotropic elastic cuboid. The transversely isotropic plane has an arbitrary orientation, and the crack surface is along an inclined plane. The mixed 3D SIFs are evaluated by using the asymptotical relation between the SIFs and the relative crack opening displacements.

Numerical results show clearly the influence of the material and crack orientations on the mixed‐mode SIFs. For comparison, the mode‐I SIF when a horizontal rectangular crack is embedded entirely within the cuboid is calculated also. It is observed that the SIF values along the crack front are larger when the crack is closer to the surface of the cuboid than those when the crack is far away from the surface.

The FORTRAN program developed is limited to regular surface cracks which can be discretized by the quadrilateral shape function; it is not very efficient and suitable for irregular crack shapes.

The evaluation of the 3D mixed‐mode SIFs in the transversely isotropic material may have direct practical applications. The SIFs have been used in engineering design to obtain the safety factor of the elastic structures.

This is the first time that the special nine‐node quadrilateral shape function has been applied to the boundary containing the crack mouth. The numerical method developed can be applied to the SIF calculation in a finite transversely isotropic cuboid within an inclined surface crack. The computational approach and the results of SIFs are of great value for the modeling and design of anisotropic elastic structures.

Wavelets are being increasingly used for damage diagnostics. The purpose of this paper is to present an algorithm that uses the wavelet transform for detecting mixed-mode, also known as combined mode, cracks in large truss structures.

The mixed-mode crack is modeled by superposing two damage modes, and this model is combined with a finite element model of the truss. The natural modes of the truss are processed through the wavelet transform and then used to determine the damage location. The influence of multiple parameters such as truss geometry, crack geometry, number of truss members, orientation of truss members, etc. is investigated as part of the study.

The proposed damage detection algorithm is found to be successful in detecting single mode as well as mixed-mode cracks even in the presence of significant end effects, and even when a relatively coarse sampling of natural modes is used. Results from multiple simulations that involve three commonly used truss structures are presented. A correlation between damage severity and the magnitude of wavelet coefficients is observed.

The proposed algorithm is found to be successful in accurately detecting damage, but direct determination of damage severity is found to be challenging.

Since the most of structures and structural components suffers from cyclic loadings, the study on the fatigue failure due to the crack growth has a great importance. The purpose of this paper is to present a three-dimensional fatigue crack growth simulation of embedded cracks using s-version finite element method (SFEM). Using the numerical results, the validity of the fitness-for-service (FFS) code evaluation method is verified.

In this paper, three-dimensional fatigue crack propagation analysis of embedded cracks is performed using the SFEM. SFEM is a numerical analysis method in which the shape of the structure is represented by a global mesh, and cracks are modeled by local meshes independently. The independent global and local meshes are superimposed to obtain the displacement solution of the problem simultaneously.

The fatigue crack growth of arbitrary shape of cracks is slow compared to that of the simplified circular crack and the crack approximated based on the FFS code of the Japan Society of Mechanical Engineers (JSME). The results tell us that the FFS code of JSME can provide a conservative evaluation of the fatigue crack growth and the residual life time.

This paper presents a three-dimensional fatigue crack growth simulation of embedded cracks using SFEM. Using this method, it is possible to apply mixed mode loads to complex shaped cracks that are closer to realistic conditions.

This study provides an analysis of patch repair for cracked aircraft structures. Delamination is a type of damage that affects the patch's behavior. The purpose of this study is to assess the influence of delamination on repair performance.

An analytical and numerical study using the finite element method was conducted for a cracked plate repaired with a patch containing a pre-existing delamination defect. The method for defining the contact pair surfaces and modeling the delamination interaction within the patch interface is specified using the virtual crack closure technique (VCCT) approach.

The efficiency of the repair is measured in terms of the J-integral. The effects of delamination initiation, mechanical loading, crack length and patch stacking sequences are presented. It is noted that in mode I, delamination propagation is only significant at node A. The numerical results are in good agreement with those of the analytical solution found in the literature. It is observed that the patch's behavior is strongly dependent on loading, crack size and stacking sequences in terms of reducing the structure's lifespan, especially in the presence of delamination.

The numerical modeling presented by the VCCT approach is highly valuable for studying delamination evolution. The influence of loading, crack size and stacking sequences on repair performance is discussed in this work.

This paper aims to provide SIF calculation method for engineering application.

In this paper, the stress intensity factors (SIFs) calculation method is applied to the anisotropic Ni-based single crystal film cooling holes (FCHs) structure.

Based on contour integral, the anisotropic SIFs analysis finite element method (FEM) in Ni-based single crystal is proposed. The applicability and mesh independence of the method is assessed by comparing the calculated SIFs using mode of plate with an edge crack. Anisotropic SIFs can be calculated with excellent accuracy using the finite element contour integral approach. Then, the effect of crystal orientation and FCHs interference on the anisotropic SIFs is clarified. The SIFs of FCH edge crack in the [011] orientated Ni-based single crystal increases faster than the other two orientations. And the SIF of horizontal interference FCHs edge crack is also larger than that of the inclined interference one.

The SIFs of the FCH edge crack in the turbine air-cooled blade are innovatively computed using the sub-model method. Both the Mode I and II SIFs of FCHs edge crack in blade increase with crack growing.

Behavior of mode I crack tip in fiber metal laminate (FML) differs from that in homogeneous or plain specimen made of metal used in the laminate due to the load transfer effect in the laminate caused by property mismatch between dissimilar material layers. The purpose of this paper is to present a finite element investigation on the characteristics of crack tip in monotonically loaded and residually stressed FML.

Crack tip characteristics are assessed by: the sizes of various zones that form at the tip; and crack tip energy release rates. The same are found by modeling two types of Glare laminates under monotonic tension with different crack orientations in SSY regime – Type I and Type II. Residual stresses are externally introduced in the models. Delaminations are modeled by cohesive elements. Crack tip zone sizes are measured from finite element solutions. Values of J integrals are computed over cyclic paths near the crack tips. Identically cracked and loaded plain aluminum alloy specimens are also modeled for comparison.

The sizes of crack tip zones in Glare laminates are found to be different than those in plain specimens. Process zone is observed to form at crack tip in Type I laminate whereas it does not develop in Type II laminate, the reverse being true in plain specimens. Values of J integrals near crack tips are also found to deviate from those in plain specimens, higher in Type I laminate due to crack tip stress amplification and lower in Type II laminate due to stress reduction. Crack orientation decides the amplification or shielding effect in the laminate.

There is scope for validating the numerical results reported in the paper by theoretical models.

The method to quantify crack tip shielding and amplification is presented that shall be useful in checking the structural integrity/safety of the laminate during actual service conditions.

Shielding and amplification effects are explicitly described and illustrated in the paper. Suitability of using J integrals over paths crossing non-homogeneous and property mismatched material layers is tested. Use of cohesive zone method that is readily applicable in finite element procedures and is relatively simple, fast and reasonably accurate is also demonstrated.

Based on the non-local piezoelectricity theory, this paper is concerned with two collinear permeable Mode-I cracks in piezoelectric materials subjected to the harmonic stress wave. The paper aims to discuss this issue.

According to the Fourier transformation, the problem is formulated into two pairs of dual integral equations, in which the unknown variables are the displacement jumps across the crack surfaces.

Finally, the dynamic non-local stress and the dynamic non-local electric displacement fields near the crack tips are obtained. Numerical results are provided to illustrate the effects of the distance between the two collinear cracks, the lattice parameter and the circular frequency of the incident waves on the entire dynamic fields near the crack tips, which play an important role in designing new structures in engineering.

Different from the classical solutions, the present solution exhibits no stress and electric displacement singularities at the crack tips in piezoelectric materials. It is found that the maximum stress and maximum electric displacement can be used as a fracture criterion.

Adhesively bonded joints are gaining importance in the structural joining processes competing against welding and bolting processes. However, long-term behaviour of adhesively bonded joints is still an open question. Due to the increasing interest in adhesively bonded joints, mainly in the transports industry, there is a need to deep the knowledge about the fatigue behaviour of adhesive joints with metallic substrates allowing the development of reliable joints to resist cyclic loadings. The paper aims to discuss these issues.

An experimental research aiming at characterizing the fatigue behaviour of adhesively bonded aluminium substrates is presented in this paper, covering both fatigue crack propagation and global S-N behaviours. Double cantilever beam (DCB), end notch flexure (ENF) and double lap joints (DLJ) specimens built using the AA6061T651 substrate and epoxy adhesive were used to evaluate the pure modes I and II fatigue crack propagation rates and the S-N fatigue behaviours.

DCB and ENF specimens allowed the formulation of pure modes I and II fatigue crack propagation laws including the propagation thresholds. DLJs showed higher static shear strength than recommended by the manufacturer for aluminium substrates, but fatigue resistance of the DLJs was lower than suggested by the manufacturer. The fatigue damage process in the DLJs was dominated by a fatigue crack initiation process.

A consistent fatigue research on adhesively bonded aluminium substrates is presented covering in the same study aspects of fatigue crack propagation and fatigue crack initiation. Data reduction schemes involving both numerical and analytical procedures were followed. Proposed work constitutes a rigorous basis for future fatigue prediction models developments.

The purpose of this paper is to study theoretically the ability of the prestressed foam core composite sandwich Split Cantilever Beam (SCB) for generating mixed-mode II/III crack loading conditions (the mode II fracture was provided by prestressing the beam using imposed transverse displacements).

The concepts of linear-elastic fracture mechanics were used. The fracture behavior was studied in terms of the strain energy release rate. For this purpose, a three-dimensional finite element model of the prestressed sandwich SCB was developed. The virtual crack closure technique was applied in order to analyze the strain energy release rate mode components distribution along the crack front.

It was found that the distribution is non-symmetric. The analysis revealed that a wide mixed-mode II/III ratios range can be generated by varying the magnitude of the imposed transverse displacement. The influence of the sandwich core material on the mixed-mode II/III fracture behavior was investigated. For this purpose, three sandwich beam configurations with different rigid cellular foam core were simulated. It was found that the strain energy release rate decreases when the foam core density increases.

For the first time, a mixed-mode II/III fracture study of foam core composite sandwich beam is performed.

It has been usual to prefer an enrichment pattern independent of the mesh when applying singular functions in the Generalized/eXtended finite element method (G/XFEM). This choice, when modeling crack tip singularities through extrinsic enrichment, has been understood as the only way to surpass the typical poor convergence rate obtained with the finite element method (FEM), on uniform or quasi-uniform meshes conforming to the crack. Then, the purpose of this study is to revisit the topological enrichment strategy in the light of a higher-order continuity obtained with a smooth partition of unity (PoU). Aiming to verify the smoothness' impacts on the blending phenomenon, a series of numerical experiments is conceived to compare the two GFEM versions: the conventional one, based on piecewise continuous PoU's, and another which considers PoU's with high-regularity.

The stress approximations right at the crack tip vicinity are qualified by focusing on crack severity parameters. For this purpose, the material forces method originated from the configurational mechanics is employed. Some attempts to improve solution using different polynomial enrichment schemes, besides the singular one, are discussed aiming to verify the transition/blending effects. A classical two-dimensional problem of the linear elastic fracture mechanics (LEFM) is solved, considering the pure mode I and the mixed-mode loadings.

The results reveal that, in the presence of smooth PoU's, the topological enrichment can still be considered as a suitable strategy for extrinsic enrichment. First, because such an enrichment pattern still can treat the crack independently of the mesh and deliver some advantage in terms of convergence rates, under certain conditions, when compared to the conventional FEM. Second, because the topological pattern demands fewer degrees of freedom and impacts conditioning less than the geometrical strategy.

Several outputs are presented, considering estimations for the J–integral and the angle of probable crack advance, this last computed from two different strategies to monitoring blending/transition effects, besides some comments about conditioning. Both h- and p-behaviors are displayed to allow a discussion from different points of view concerning the topological enrichment in smooth GFEM.