Search results

A review is made of methods for calculating parameters characterizing crack tip behaviour in non‐linear materials. Convenient methods of calculating J‐integral type quantities are reviewed, classified broadly into two groups, as domain integrals and virtual crack extension techniques. In addition to considerations of how such quantities may be calculated by finite elements, assessment methods of conducting the actual incremental analyses are described.

This study aims to contribute to the numerical modelling of drop impact on a flip-chip component assembled on printed circuit boards using solder micro-bumps. This contribution is based on the introduction of non-linear fracture mechanics in the numerical approach.

The integration of non-linear fracture mechanics into the numerical approach requires the proposal and validation of several simplifying assumptions. Initially, a dynamic 3D model was simplified to a dynamic 2D model. Subsequently, the dynamic 2D model is replaced with an equivalent static 2D model. The equivalent static 2D model was used to perform calculations considering the non-linear fracture mechanics. A crack was modelled in the critical bump. The J-integral was used as a comparative parameter to study the effects of crack length, crack position and chip thickness on the fracture toughness of the solder bump.

The different simplifying assumptions were validated by comparing the results obtained by the various models. Numerical results showed a high risk of failure at the critical solder bump in a zone close to the intermetallic layer. The obtained results were in agreement with the post-test observations using the “Dye and Pry” methods.

The originality of this study lies in the introduction of non-linear fracture mechanics to model the mechanical response of solder bumps during drop impact. This study led to some interesting conclusions, highlighting the advantage of introducing non-linear fracture mechanics into the numerical simulations of microelectronic components during a drop impact.

The purpose of this paper is to develop basic principles of deterministic structural integrity assessment of a component with a crack- or notch-like defect by including safety factors against fracture and plastic collapse in criteria equations of linear and nonlinear fracture mechanics.

The safety factors against fracture are calculated by demanding that the applied critical stress should not be less than the yield stress of the material for a component with a crack or a notch of the acceptable size. Structural integrity assessment of the engineering components damaged by crack- or notch-like defects is discussed from view point of the failure assessment diagram (FAD). The methodology of the FAD has been employed for the structural integrity analysis and assessment of acceptable sizes of throw-thickness notch in a plate under tension and surface longitudinal notch-like defects in a pressure vessel.

Basic equations have been presented to calculate the safety factor against fracture for critical values of the stress intensity factor, crack tip opening displacement (CTOD), the J-integral and the FAD as well as to estimate an acceptable (safe) region for an engineering component with a crack- or notch-like defect of the acceptable size. It was shown that safety factors against fracture depend on both the safety factor against plastic collapse and employed fracture mechanics criterion. The effect of crack/notch tip constraint is incorporated into criteria equations for the calculation of safety factors against fracture.

The deterministic method of fracture mechanics is recommended for structural integrity assessment of a component with a crack- or notch-like defect by including safety factors against fracture and plastic collapse in criteria equations of linear and nonlinear fracture mechanics.

The purpose of this paper is to further develop the truss‐like discrete element method (DEM) in order to make it suitable to deal with damage and fracture problems.

Finite and boundary elements are the best developed methods in the field of numerical fracture and damage mechanics. However, these methods are based on a continuum approach, and thus, the modelling of crack nucleation and propagation could be sometimes a cumbersome task. Besides, discrete methods possess the natural ability to introduce discontinuities in a very direct and intuitive way by simply breaking the link between their discrete components. Within this context, the present work extends the capabilities of a truss‐like DEM via the introduction of three novel features: a tri‐linear elasto‐plastic constitutive law; a methodology for crack discretization and the computation of stress intensity factors; and a methodology for the computation of the stress field components from the unixial discrete‐element results.

Obtained results show the suitability and the performance of the proposed methodologies to solve static and dynamic crack problems (including crack propagation) in brittle and elasto‐plastic materials. Computed results are in good agreement with experimental and numerical results reported in the bibliography.

This paper demonstrates the versatility of the truss‐like DEM to deal with damage mechanics problems. The approach used in this work can be extended to the implementation of time‐dependent damage mechanisms. Besides, the capabilities of the discrete approach could be exploited by coupling the truss‐like DEM to finite and boundary element methods. Coupling strategies would allow using the DEM to model the regions of the problem where crack nucleation and propagation occurs, while finite or boundary elements are used to model the undamaged regions.

The scope of the truss‐like DEM has been extended. New procedures have been introduced to deal with elastoplastic‐crack problems and to improve the post processing of the stress results.

Contact problems are among the most difficult ones in mechanics. Due to its practical importance, the problem has been receiving extensive research work over the years. The finite element method has been widely used to solve contact problems with various grades of complexity. Great progress has been made on both theoretical studies and engineering applications. This paper reviews some of the main developments in contact theories and finite element solution techniques for static contact problems. Classical and variational formulations of the problem are first given and then finite element solution techniques are reviewed. Available constraint methods, friction laws and contact searching algorithms are also briefly described. At the end of the paper, a bibliography is included, listing about seven hundred papers which are related to static contact problems and have been published in various journals and conference proceedings from 1976.

Mode I static fracture behavior of polymer composites is studied using the tapered double cantilever beam test method. A non‐linear three‐dimensional finite element model is developed to analyze the test data. The fracture toughness is evaluated using a J‐integral approach. A non‐uniform distribution of the J‐integral value along the crack front is obtained with maximum at the mid‐plane of the specimen. It is shown that taking into account the damage induced non‐linear behavior improves the fracture toughness. This is explained with increased strain energy dissipation as a result of the non‐linear behavior.

This paper gives a bibliographical review of the finite element and boundary element parallel processing techniques from the theoretical and application points of view. Topics include: theory – domain decomposition/partitioning, load balancing, parallel solvers/algorithms, parallel mesh generation, adaptive methods, and visualization/graphics; applications – structural mechanics problems, dynamic problems, material/geometrical non‐linear problems, contact problems, fracture mechanics, field problems, coupled problems, sensitivity and optimization, and other problems; hardware and software environments – hardware environments, programming techniques, and software development and presentations. The bibliography at the end of this paper contains 850 references to papers, conference proceedings and theses/dissertations dealing with presented subjects that were published between 1996 and 2002.

The purpose of this paper is to systematically investigate the fluid lag phenomena and its influence in the hydraulic fracturing process, including all stages of fluid-lag evolution, the transition between different stages and their coupling with dynamic fracture propagation under common conditions.

A plane 2D model is developed to simulate the complex evolution of fluid lag during the propagation of a hydraulic fracture driven by an impressible Newtonian fluid. Based on the finite element method, a fully implicit solution scheme is proposed to solve the strongly coupled rock deformation, fluid flow and fracture propagation. Using the proposed model, comprehensive parametric studies are performed to examine the evolution of fluid lag in various geological and operational conditions.

The numerical simulations predict that the lag ratio is around 5% or even lower at the beginning stage of hydraulic fracture under practical geological conditions. With the fracture propagation, the lag ratio keeps decreasing and can be ignored in the late stage of hydraulic fracturing for typical parameter combinations. On the numerical aspect, whether the fluid lag can be ignored depends not only on the lag ratio but also on the minimum mesh size used for fluid flow. In addition, an overall mixed-mode fracture propagation factor is proposed to describe the relationship between diverse parameters and fracture curvature.

In this study, relatively simple physical models such as linear elasticity for solid, Newtonian model for fluid and linear elasticity fracture mechanics for fracture are used. The current model does not account for such effects like leak off, poroelasticity and softening of rock formations, which may also visibly affect the fluid lag depending on specific reservoir conditions.

This study helps to understand the effect of fluid lag during hydraulic fracturing processes and provides numerical experience in dealing with the fluid lag with finite element simulation.

The purpose of this paper is to perform a theoretical analysis of non-linear delamination fracture in cantilever beam opened notch (CBON) configuration. It is assumed that the non-linear mechanical behavior of the CBON can be described by using a stress-strain curve with power-law hardening.

The fracture analysis is carried-out by applying the integration contour independent J-integral. For this purpose, a model based on the technical beam theory is used. Equation is derived for determination of the CBON specimen curvature in elastic-plastic stage of deformation. The equation is solved by using the MatLab program system. Solutions of the J-integral are obtained at linear-elastic as well as elastic-plastic behavior of the CBON. The influence of the power-law exponent on the non-linear fracture is evaluated.

The analysis reveals that the J-integral value increases when the exponent of the power-law increases. The solution obtained here is very useful for parametric analyses of the non-linear fracture behavior, since the simple formulas derived capture the essentials of the fracture response.

Beside for parametric investigations, the solution obtained here can also be applied for calculating the critical J-integral value at non-linear behavior using experimentally determined critical fracture load at the onset of crack growth from the initial crack tip position in the CBON configuration.

An analysis is performed of the non-linear fracture in the CBON configuration by applying the J-integral approach, assuming that the mechanical response can be modeled using a stress-strain curve with power-law hardening.

This paper presents a mathematical model for the analysis of cohesive fracture propagation through a non‐homogeneous porous medium. Governing equations are stated within the frame of Biot's theory, accounting for the flow through the solid skeleton, along the fracture and across its sides toward the surrounding medium. The numerical solution is obtained in a 2D context, exploiting the capabilities of an efficient mesh generator, and requires continuous updating of the domain as the fractures enucleate and propagate. It results that fracture paths and their velocity of propagation, usually assumed as known, are supplied directly by the model without introducing any simplifying assumption.