Search results

This paper gives a bibliographical review of the finite element methods (FEMs) applied to the analysis of ceramics and glass materials. The bibliography at the end of the paper contains references to papers, conference proceedings and theses/dissertations on the subject that were published between 1977‐1998. The following topics are included: ceramics – material and mechanical properties in general, ceramic coatings and joining problems, ceramic composites, ferrites, piezoceramics, ceramic tools and machining, material processing simulations, fracture mechanics and damage, applications of ceramic/composites in engineering; glass – material and mechanical properties in general, glass fiber composites, material processing simulations, fracture mechanics and damage, and applications of glasses in engineering.

This study aims to comprehensively investigate the mixed-mode fracture behavior and mechanical properties of selective laser sintering (SLS) polyamide 12 (PA12) components, considering different build orientations and layer thicknesses. The primary objectives include the following. Conducting mixed-mode fracture and mechanical analyses on SLS PA12 parts. Investigating the influence of build orientation and layer thickness on the mechanical properties of SLS-printed components. Examining the fracture mechanisms of SLS-produced Arcan fracture and tensile specimens through experimental methods and finite element analyses.

The research used a combination of experimental techniques and numerical analyses. Tensile and Arcan fracture specimens were fabricated using the SLS process with varying build orientations (X, X–Y, Z) and layer thicknesses (0.1 mm, 0.2 mm). Mechanical properties, including tensile strength, modulus of elasticity and critical stress intensity factor, were quantified through experimental testing. Mixed-mode fracture tests were conducted using a specialized fixture, and finite element analyses using the J-integral method were performed to calculate fracture toughness. Scanning electron microscopy (SEM) was used for detailed morphological analysis of fractured surfaces.

The investigation revealed that the highest tensile properties were achieved in samples fabricated horizontally in the X orientation with a layer thickness of 0.1 mm. Additionally, parts manufactured with a layer thickness of 0.2 mm exhibited favorable mixed-mode fracture behavior. The results emphasize the significance of build orientation and layer thickness in influencing mechanical properties and fracture behavior. SEM analysis provided valuable insights into the failure mechanisms of SLS-produced PA12 components.

This study contributes to the field of additive manufacturing by providing a comprehensive analysis of the mixed-mode fracture behavior and mechanical properties of SLS-produced PA12 components. The investigation offers novel insights into the influence of build orientation and layer thickness on the performance of such components. The combination of experimental testing, numerical analyses and SEM morphological observations enhances the understanding of fracture behavior in additive manufacturing processes. The findings contribute to optimizing the design and manufacturing of high-quality PA12 components using SLS technology.

Finite element (FE) modeling of the human dentate mandible is the method of choice currently used for simulating structural fracture analyses in the mandibular region. A finite element model of a parasymphyseal fracture with an internal rigid fixation plate‐screw system has been developed to compare the effects of including frictionless/frictional contact boundary conditions at the fracture site. It is common practice to ignore contact boundary conditions in FE modeling of mandibular fractures due to the non‐linearities causing increased computational requirements. The stress distributions and displacements of the mandibular fracture region indicate a significant difference resulting from the introduction of realistic contact boundary conditions. These current findings suggest that even though the modeling of extreme situations, i.e. non‐contact modeling of unhealed fractures, may provide insight to non‐union problems, future mandibular fracture models should include frictional contact boundary conditions. This is in order to capture more realistic behavior of the system to be analyzed.

The purpose of this paper is to compare fracture‐fixation and bone‐healing promotion efficacies of an intramedullary (IM) nail‐type and an external osteosynthesis plate‐type femoral trochanteric‐fracture implants using the results of a combined multi‐body dynamics and finite element analyses. For both implants, fracture fixation was obtained using a dynamic hip blade which is anchored to the femur head on one end and is connected to the IM rod/plate on the other. The analysis was carried out for two pre‐fracture conditions of the femur: healthy and osteoporotic.

The musculoskeletal dynamics portion of the analysis was used to obtain realistic physiological loading conditions (i.e. muscle forces and joint reaction forces and moments) associated with four typical everyday activities of a patient, namely, walking, lunging, cycling and egress (i.e. exiting a passenger vehicle). The subsequent structural finite element analysis of the fractured femur/implant assembly was employed to quantify fracture‐fixation efficacy (as measured by the extents of lateral (found to be minor), flexural and torsional displacements of the two femur fragments) and the bone‐healing promotion efficacy (as quantified by the fraction of the fractured surface area which experienced desirable contact pressures).

The results obtained show that, in general, the IM‐rod type of implant out‐performs the osteosynthesis plate type of implant over a large range of scenarios involving relative importance of the bone‐healing promotion and fracture‐fixation efficacies, health condition of the femur and the activity level of the patient. More specifically, the more active the patient and the larger extent of osteoporosis in the femur, the more justifiable is the use of the IM‐rod type of implant.

The present approach enables assessment of the fracture‐fixation performance of orthopedic implants under physiologically realistic loading conditions.

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.

Describes a probabilistic methodology for fracture assessments of flawed structures constructed of ferritic steels using the research code WSTRESS. The probabilistic formulation for cleavage fracture implements a multiaxial form of the weakest link model which couples the macroscopic fracture behavior with a micromechanics model based on the statistics of microcracks. The Weibull stress, σ_w, emerges as a suitable near‐tip parameter to provide a connection between the microregime of failure and remote loading (J). WSTRESS builds on an iterative procedure to incorporate a 3‐D finite element description of the crack‐tip stress field and measured values of fracture toughness to calibrate the Weibull modulus, m, and the scale parameter, σ_u. Specific features of the code include statistical inference of Weibull parameters based on uncensored and censored models (with maximum likelihood method), construction of confidence intervals, several definitions for the near‐tip fracture process zone and other general facilities such as spatial integration of the stress field (to incorporate the random orientation of microcracks) and stochastic simulation of fracture data using the Monte Carlo method. The code also includes a convenient free‐form command language and a seamless interface with finite element results files stored in Patran binary or ASCII format.

The localized strain softening behaviour of concrete has been modelled by two approaches: (i) the stiffness degrading model based on the total stress‐strain constitutive relationship, and (ii) the tangent softening model based on the incremental stress‐strain relationship. The models are implemented using a new softening initiation criterion proposed for application in multi‐dimensional finite element analysis. Parametric analyses on plain concrete beams, tested experimentally by other researchers, have been carried out to investigate the required numerical efforts, the mesh objectivity, and the energy dissipation characteristics of the structures. The stiffness degrading model is very stable even when applied with relatively coarse finite element meshes. However, the computational demand of this model is relatively high. The combination of a total stress‐strain constitutive relationship to compute the element responses, and an incremental relationship to formulate the stiffness matrix, appears to be computationally efficient and stable, provided that adequately refined finite element mesh is used to model the structure.

A planar lattice network of beam elements is used to study the mechanisms of fracture in cement‐based materials. Beam properties are controlled by a nonlinear elastic fracture law which roughly accounts for three‐dimensionality of the material and fracture process. Special attention is given to modeling toughening mechanisms associated with aggregate‐matrix interface failure. The distributions of damage and fracture energy consumption are resolved at the material mesoscale and are shown to depend on strain gradient. An adaptive remeshing procedure is used to reduce computational cost and enable analyses of specimens of significantly differing scale, while keeping the lattice density constant. Larger process zones, higher specific fracture energies, and lower specific peak loads are obtained with increasing specimen size, in agreement with published test results. These computations provide information useful in developing refined macromodels for engineering analyses.

The unstable dynamic propagation of multistage hydrofracturing fractures leads to uneven development of the fracture network and research on the mechanism controlling this phenomenon indicates that the stress shadow effects around the fractures are the main mechanism causing this behaviour. Further studies and simulations of the stress shadow effects are necessary to understand the controlling mechanism and evaluate the fracturing effect.

In the process of stress-dependent unstable dynamic propagation of fractures, there are both continuous stress fields and discontinuous fractures; therefore, in order to study the stress-dependent unstable dynamic propagation of multistage fracture networks, a series of continuum-discontinuum numerical methods and models are reviewed, including the well-developed extended finite element method, displacement discontinuity method, boundary element method and finite element-discrete element method.

The superposition of the surrounding stress field during fracture propagation causes different degrees of stress shadow effects between fractures and the main controlling factors of stress shadow effects are fracture initiation sequence, perforation cluster spacing and well spacing. The perforation cluster spacing varies with the initiation sequence, resulting in different stress shadow effects between fractures; for example, the smaller the perforation cluster spacing and well spacing are, the stronger the stress shadow effects are and the more seriously the fracture propagation inhibition arises. Moreover, as the spacing of perforation clusters and well spacing increases, the stress shadow effects decrease and the fracture propagation follows an almost straight pattern. In addition, the computed results of the dynamic distribution of stress-dependent unstable dynamic propagation of fractures under different stress fields are summarised.

A state-of-art review of stress shadow effects and continuum-discontinuum methods for stress-dependent unstable dynamic propagation of multiple hydraulic fractures are well summarized and analysed. This paper can provide a reference for those engaged in the research of unstable dynamic propagation of multiple hydraulic structures and have a comprehensive grasp of the research in this field.

The interaction between hydraulic fracture (HF) and natural fracture (NF) in naturally fractured rocks is critical for hydraulic fracturing. This paper aims to focus on investigating the development of tensile and shear debonding zone on the NF caused by the stresses produced by HF, and the influence of NF’s debonding behavior on the interaction between HF and NF.

Theoretically, tensile and shear debonding modes of NF are considered, two dimensionless parameters are proposed to characterize the difficulty of tensile and shear failure of NF, respectively. Numerically, a finite element model combining the extended finite element method and cohesive zone method (CZM) is proposed to study NF’s debonding behavior and its influence on the interaction between HF and NF.

Both theoretical analysis and numerical simulation show the existence of two debonding modes. The numerical results also show that the HF can cross, offset or propagate along the NFs depending on the parameters’ value, resulting in different fracture network and stimulated reservoir volume. When they are large, the NF’s debonding area is small, HF tends to cross the NF and the fracture network is simple; when they are small, the NF’s debonding area is large, HF will propagate along the NF. In addition, HF is easier to propagate along with NF under tensile debonding mode while it is easier to pass through NF under shear debonding mode.

The theoretical and numerical considerations are taken into account in the influence of the debonding of NFs on the interaction between HFs and NFs and the influence on the formation of the fracture network.