Search results

The purpose of this study is to develop an efficient numerical procedure for simulating the effect of printing orientation, as one of the primary sources of anisotropy in 3D-printed components, on their fracture properties.

The extended finite element method and the cohesive zone model (XFEM-CZM) are used to develop subroutines for fracture simulation. The ability of two prevalent models, i.e. the continuous-varying fracture properties (CVF) model and the weak plane model (WPM), and a combination of both models (WPM-CVF) are evaluated to capture fracture behavior of the additively manufactured samples. These models are based on the non-local and local forms of the anisotropic maximum tangential stress criterion. The numerical models are assessed by comparing their results with experimental outcomes of 16 different configurations of polycarbonate samples printed using the material extrusion technique.

The results demonstrate that the CVF exaggerates the level of anisotropy, and the WPM cannot detect the mild anisotropy of 3D-printed parts, while the WPM-CVF produces the best results. Additionally, the non-local scheme outperforms the local approach in terms of finite element analysis performance, such as mesh dependency, robustness, etc.

This paper provides a method for modeling anisotropic fracture in 3D-printed objects. A new damage model based on a combination of two prevalent models is offered. Moreover, the developed subroutines for implementing the non-local anisotropic fracture criterion enable a reliable crack propagation simulation in media with varying degrees of complication, such as anisotropy.

Due to its inexpensive production costs, low stress concentration and maintenance-friendliness, the adhesive bonded pipe joint is frequently utilized for pipe connection. However, further theoretical analysis is needed to understand the debonding failure mechanism of such bonded pipe joints under axial tension.

In this study, based on the bi-linear cohesive zone model, the integrated closed-form solutions were derived by considering the axial stiffness ratio and failure stage to determine the relative interfacial slip, interfacial shear stress and relationship of tension–displacement in the bonded pipe joint.

Additionally, solutions for the critical bonded length and the ultimate load capacity were put forth. Besides, the numerical study was conducted to verify the theoretical solutions regarding the load–displacement relationship. The interfacial shear stress distribution at different failure stages was presented to understand the interfacial shear stress transmission and debonding process. The effect of bonded length on the ultimate load and ductility of pipe joints was also discussed.

The findings in this study can give a reference for the design of bonded pipe joints in their actual engineering applications.

Pipelines are subject to pits, holes and cracks after staying in service for a while, especially in harsh environments. To repair the pipelines, composite materials are used, due to composite materials' low cost, high-corrosion resistance and easy handling. This study aims to investigate the reliability of the blister test for evaluating the bonding strength of multiwall carbon nanotube (MWCNT) on woven carbon-reinforced epoxy.

Flexural, hardness and Izod impact tests were used to evaluate MWCNT effect on the epoxy by adding different amounts, 0.2, 0.4, 0.6, 0.8 and 1 wt. %, of MWCNT, to be compared with pure epoxy.

The results showed that 0.8 wt.% gives the highest strength. The experimental results of 0.8 wt.% MWCNT reinforced carbon composite was compared with the finite element model under blister test, and the results showed high similarities.

Evaluation of the reliability and the advantages of MWCNT considering the high aspect ratio and high tensile strength, which is more than 15 times compared to steel, MWCNT enhances the strength, stiffness and toughness of epoxy used as a matrix in repairing pipelines, which leads to an increase in the resistance of composite materials against oil internal pressure before delamination.

This study implements the fourth-order phase field method (PFM) for modeling fracture in brittle materials. The weak form of the fourth-order PFM requires C¹ basis functions for the crack evolution scalar field in a finite element framework. To address this, non-Sibsonian type shape functions that are nonpolynomial types based on distance measures, are used in the context of natural neighbor shape functions. The capability and efficiency of this method are studied for modeling cracks.

The weak form of the fourth-order PFM is derived from two governing equations for finite element modeling. C⁰ non-Sibsonian shape functions are derived using distance measures on a generalized quad element. Then these shape functions are degree elevated with Bernstein-Bezier (BB) patch to get higher-order continuity (C¹) in the shape function. The quad element is divided into several background triangular elements to apply the Gauss-quadrature rule for numerical integration. Both fourth-order and second-order PFMs are implemented in a finite element framework. The efficiency of the interpolation function is studied in terms of convergence and accuracy for capturing crack topology in the fourth-order PFM.

It is observed that fourth-order PFM has higher accuracy and convergence than second-order PFM using non-Sibsonian type interpolants. The former predicts higher failure loads and failure displacements compared to the second-order model due to the addition of higher-order terms in the energy equation. The fracture pattern is realistic when only the tensile part of the strain energy is taken for fracture evolution. The fracture pattern is also observed in the compressive region when both tensile and compressive energy for crack evolution are taken into account, which is unrealistic. Length scale has a certain specific effect on the failure load of the specimen.

Fourth-order PFM is implemented using C¹ non-Sibsonian type of shape functions. The derivation and implementation are carried out for both the second-order and fourth-order PFM. The length scale effect on both models is shown. The better accuracy and convergence rate of the fourth-order PFM over second-order PFM are studied using the current approach. The critical difference between the isotropic phase field and the hybrid phase field approach is also presented to showcase the importance of strain energy decomposition in PFM.

In this paper, an attempt is made to obtain buckling loads, ultimate bearing capacity and other required structural characteristics of grid structure panels. The numerical method for post-buckling behavior analysis of panels involving multiple invisible damages is also presented.

In this paper, two bidirectional stiffened composite panels are manufactured and tested. Multiple discrete invisible damages are introduced in different positions of the stringers, and the experimental and simulation investigation of buckling and post-buckling were carried out on the damaged stiffened panels.

The simulation load–displacement curves are compared with the experimental results, and it is found that the simulation model can well predict the occurrence of buckling and failure loads. The strain curve shows that the rate of strain change at the damaged site is greater than that at the undamaged site, which reflects that the debond is more likely occurred at the damaged site. The simulation verifies that the panel is usually crushed due to matrix compression and fiber–matrix shear.

In this paper, post-buckling tests and numerical simulations of bidirectional stiffened composite panels with impact damage were carried out. Two panels with four longitudinal stringers and two transverse stringers were manufactured and tested. The buckling and post-buckling characteristics of the grid structure are obtained, and the failure mechanism of the structure is explained. This is helpful for the design of wall panel structure.

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.