A Finite Element Study of Crack Behavior for Carbon Na-notube Reinforced Bone Cement

Polymethylmethacrylate (PMMA) bone cement is a polymeric material that is widely used as a structural orthopedic material. However, it is not an ideal material for bone grafting due to its fragility. Carbon nanotubes (CNTs) have been introduced in order to reinforce PMMA resulting in a composite material which exhibits improved tensile properties, increased fatigue resistance and fracture toughness. This improvement is potentially due to bridging and arresting cracks as well as absorption of energy. In this study, a two-dimensional finite element model is presented for the fracture analysis of PMMA-CNT composite material. Instead of the classical single fiber model, the present work considers an ensemble of CNTs interacting with a pre-existing crack. Casca is used to produce a two dimensional mesh and the fracture analysis is performed using Franc 2D. The model is subjected to uni-axial loading in the transverse plane and the interaction between the crack and CNTs is evaluated by determining the stress intensity factor in the vicinity of the crack tips. The effects of geometric parameters of the CNTs and the material structural heterogeneity on crack propagation trajectory are investigated. Furthermore, the effects of CNT diameter, wall thickness and elastic mismatch between the matrix and the nanotubes on crack growth are studied. The results illustrate that the CNTs repel cracks during loading as they act as barriers to crack growth. As a result, the incorporation of CNTs into PMMA reduces crack growth but more importantly increases the fracture resistance of bone cement.

COMBINED DELAUNAY TRIANGULATION AND ADAPTIVE FINITE ELEMENT METHOD FOR CRACK GROWTH ANALYSIS

The paper presents the utilization of the adaptive Delaunay triangulation in the finite element modeling of two dimensional crack propagation problems, including detailed description of the proposed procedure which consists of the Delaunay triangulation algorithm and an adaptive remeshing technique. The adaptive remeshing technique generates small elements around crack tips and large elements in the other regions. The resulting stress intensity factors and simulated crack propagation behavior are used to evaluate the effectiveness of the procedure. Three sample problems of a center cracked plate, a single edge cracked plate and a compact tension specimen, are simulated and their results assessed.

ANALYSIS OF ELECTRIC BOUNDARY CONDITION EFFECTS ON CRACK PROPAGATION IN PIEZOELECTRIC CERAMICS

There are three types of cracks: impermeable crack, permeable crack and conducting crack, with different electric boundary conditions on faces of cracks in piezoelectric ceramics, which poses difficulties in the analysis of piezoelectric fracture problems. In this paper, in contrast to our previous FEM formulation, the numerical analysis is based on the use of exact electric boundary conditions at the crack faces, thus the common assumption of electric impermeability in the FEM analysis is avoided. The crack behavior and elasto-electric fields near a crack tip in a PZT-5 piezoelectric ceramic under mechanical, electrical and coupled mechanical- electrical loads with different electric boundary conditions on crack faces are investigated. It is found that the dielectric medium between the crack faces will reduce the singularity of stress and electric displacement. Furthermore, when the permittivity of the dielectric medium in the crack gap is of the same order as that of the piezoelectric ceramic, the crack becomes a conducting crack, the applied electric field has no effect on the crack propagation.

Crack propagation simulation in brittle elastic materials by a phase field method

To overcome the difficulties of re-meshing and tracking the crack-tip in other computational methods for crack propagation simulations,the phase field method based on the minimum energy principle is introduced by defining a continuous phase field variable(x)∈[0,1]to characterize discontinuous cracks in brittle materials.This method can well describe the crack initiation and propagation without assuming the shape,size and orientation of the initial crack in advance.In this paper,a phase field method based on Miehe's approach[Miehe et al.,Comp.Meth.App.Mech.Eng.(2010)]is applied to simulate different crack propagation problems in twodimensional(2D),isotropic and linear elastic materials.The numerical implementation of the phase field method is realized within the framework of the finite element method(FEM).The validity,accuracy and efficiency of the present method are verified by comparing the numerical results with other reference results in literature.Several numerical examples are presented to show the effects of the loading type(tension and shear),boundary conditions,and initial crack location and orientation on the crack propagation path and force-displacement curve.Furthermore,for a single edge-cracked bi-material specimen,the influences of the loading type and the crack location on the crack propagation trajectory and force-displacement curve are also investigated and discussed.It is demonstrated that the phase field method is an efficient tool for the numerical simulation of the crack propagation problems in brittle elastic materials,and the corresponding results may have an important relevance for predicting and preventing possible crack propagations in engineering applications.

A MULTI-COUPLED FINITE ELEMENT ANALYSIS OF RESISTANCE SPOT WELDING PROCESS

A two-dimensional axisymmetric finite element model is developed to analyze the transient thermal and mechanical behaviors of the Resistance Spot Welding （RSW） process using commercial software ANSYS. Firstly a direct-coupled electrical-thermal Finite Element Analysis （FEA） is performed to analyze the transient thermal characteristics of the RSW process. Then based on the thermal results a sequential coupled thermo-elastic-plastic analysis is conducted to determine the mechanical features of the RSW process. The thermal history of the whole process and the temperature distribution of the weldment are obtained through the analysis. The mechanical features, including the distributions of the contact pressure at both the faying surface and the electrode-workpiece interface, the stress and strain distributions in the weldment and their changes during the RSW process, the deformation of the weldment and the electrode displacement are also calculated.

Finite element simulation of stress intensity factors in elastic-plastic crack growth

A finite element program developed elastic-plastic crack propagation simulation using Fortran language. At each propagation step, the adaptive mesh is automatically refined based on a posteriori h-type refinement using norm stress error estimator. A rosette of quarter-point elements is then constructed around the crack tip to facilitate the prediction of crack growth based on the maximum normal stress criterion and to calculate stress intensity factors under plane stress and plane strain conditions. Crack was modelled to propagate through the inter-element in the mesh. Some examples are presented to show the results of the implementation.

Fatigue Crack Propagation Analysis of Orthotropic Steel Bridge with Crack Tip Elastoplastic Consideration

Due to the complex structure and dense weld of the orthotropic steel bridge deck(OSBD),fatigue cracks are prone to occur in the typical welding details.Welding residual stress(WRS)will cause a plastic zone at the crack tip.In this paper,an elastoplastic constitutive model based on the Chaboche kinematic hardening model was introduced,and the extended finite element method(XFEM)was used to study the influence of material elastoplasticity and crack tip plastic zone on the law of fatigue crack propagation.By judging the stress state of the residual stress field at the crack tip and selecting different crack propagation rate models to investigate the crack propagation law when plastic deformation was considered,the propagation path and propagation rate of fatigue crack of the OSBD were obtained.The results show that,whether the residual stress field is considered or not,the plastic deformation at the crack tip will not cause the obvious closure of the fatigue crack at the U-rib toe during the crack propagation process,but will significantly affect the crack propagation path.When material plasticity is considered,the propagation angle of fatigue crack at the U-rib toe basically remains unchanged along the short-axis direction of the initial crack,but is going up along the long-axis direction,and the crack tip plastic zone inhibits the propagation of the crack tip on one side.Compared with linear elastic materials,the crack propagation law considering material plasticity is more consistent with that in actual bridge engineering.In terms of the propagation rate,if the residual stress field is not considered,the fatigue crack propagation rate at U-rib toe with plasticity considered is slightly higher than that without plasticity considered,because plastic deformation will affect the amplitude of energy release rate.When considering the WRS field,the fatigue crack propagation rate at U-rib toe is increased due to the combined actions of plastic deformation and stress ratio R.

Stable Fatigue Crack Propagation of 16MnR Steel

A research on the stable fatigue crack propagation of 16MnR steel is investigated systematically in this paper.First,control experiments of 16MnR with compact tension specimen is conducted to study the effect of R-ratios,specimen thickness and notch sizes.The experiments show that the fatigue crack growth(FCG)rate in stable propagation was insensitive to these factors.Then,the stress intensity factor(SIF)is computed and compared by displacement interpolation method,J integral and interaction integral method respectively.The simulation shows that optimization on the mesh density and the angle of singular element improved the computational efficiency and accuracy of SIF and the interaction integral method has an obvious advantage on stability.Finally,the FCG rate is modeled by the Jiang fatigue damage criterion and the extended finite element method(XFEM)respectively.The simulation results of FCG rate are in line with experiments data and indicate that XFEM method is more accurate than Jiang fatigue damage method.

Using Extended Finite Element Method for Computation of the Stress Intensity Factor, Crack Growth Simulation and Predicting Fatigue Crack Growth in a Slant-Cracked Plate of 6061-T651 Aluminum

The 6061-T651 aluminium alloy is one of the most common aluminium alloys for marine components and general structures. The stress intensity factor (SIF) is an important parameter for estimating the life of the cracked structure. In this paper, the stress intensity factors of a slant-cracked plate, which is made of 6061-T651 aluminum, have been calculated using extended finite element method (XFEM) and finite element method (FEM) in ABAQUS software and the results were compared with theoretical values. Numerical values obtained from these two methods were close to the theoretical values. In simulations of crack growth at different crack angles, the crack propagation angle values were closer to the theoretical values in XFEM method. Also, the accuracy and validity of fatigue crack growth curve were much closer to the theoretical graph in XFEM than the FEM. Therefore, in this paper the capabilities of XFEM were realized in analyzing issues such as cracks.

Crack propagation of Mg-Li alloys by impact load

Under impact load,using energy estimation and updated finite element method,the crack propagation of Mg-Li alloys was studied,and the relations between high speed impact,time interval and crack propagation velocity were deduced. The effects of impact pattern characteristic parameters such as impact range and time interval on impact fatigue crack growth and propagation rates were discussed. The results show that with the development of crack length,the crack propagation rate along the central sample alloy swiftly increases to the maximum. Then,the crack of other parts has different rate variations. The calculation and simulation results fit well with the experimental ones,which verifies the validity of energy estimation and updated finite element analysis. Impact loading usually induces materials and accelerates the fatigue crack growth.

Crack Propagation Path in Two-Directionally Graded Composites Subjected to Mixed-Mode Ⅰ+Ⅱ Loading

Crack propagation path in two-directionally graded composites was investigated by the finite element method.A graded extended finite element method(XFEM)was employed to calculate displacement and stress fields in cracked graded structures.And a post-processing subroutine of interaction energy integral was implemented to extract the mixed-mode stress intensity factors(SIFs).The maximum hoop stress(MHS)criterion was adopted to predict crack growth direction based on the assumption of local homogenization of asymptotic crack-tip fields in graded materials.Effects of material nonhomogeneous parameters on crack propagation paths were also discussed in detail.It is shown that the present method can provide relatively accurate predictions of crack paths in two-directionally graded composites.Crack propagates in the decreasing direction of effective Young′s modulus.The shape and steepness of property gradient perpendicular to the crack surface have great influences on crack paths.Through redesigning material property reasonably,crack growth in graded material can be changed to improve mechanical behaviours of cracked structures.

Early Crack Propagation Behavior of Laser Metal Deposited Ti-6Al-4V Alloy Under High Cycles Fatigue Loading

The crack initiation and early propagation are of great significance to the overall fatigue life of material.In order to investigate the anisotropic fracture behavior of laser metal deposited Ti-6Al-4V alloy(LMD Ti64)during the early stage,the fourpoint bending fatigue test was carried out on specimens of three different directions,as well as the forged specimens.The results indicate the anisotropic crack initiation and early propagation of LMD Ti64.The direction perpendicular to the deposition direction exhibits a better fatigue resistance than the other two.The crack initiation position and propagation path are dominated by the microstructure in the vicinity of U-notch.LMD Ti64 has a typical small crack effect,and the early crack propagation velocities in three directions are similar.Affected by the slip system of LMD Ti64,secondary cracks frequently occur,which are often found to have an angle of 60°to the main crack.The electron backscatter diffraction analysis indicates that LMD Ti64 has preferred orientations,i.e.,strong 0001//Z texture and 001//Z texture.Their crystallographic orientation will change as the direction of columnarβgrains turns over,resulting in the fatigue anisotropy of LMD Ti64 in crack initiation and early crack propagation process.

FINITE ELEMENT MODELLING OF COMPLEX 3D STATIC AND DYNAMIC CRACK PROPAGATION BY EMBEDDING COHESIVE ELEMENTS IN ABAQUS

This study proposes an algorithm of embedding cohesive elements in Abaqus and develops the computer code to model 3D complex cragk propagation in quasi-brittle materials in a relatively easy and efficient manner. The cohesive elements with softening traction-separation relations and damage initiation and evolution laws are embedded between solid elements in regions of interest in the initial mesh to model potential cracks. The initial mesh can consist of tetrahedrons, wedges, bricks or a mixture of these elements. Neither remeshing nor objective crack propagation criteria are needed. Four examples of concrete specimens, including a wedgesplitting test, a notched beam under torsion, a pull-out test of an anchored cylinder and a notched beam under impact, were modelled and analysed. The simulated crack propagation processes and load-displacement curves agreed well with test results or other numerical simulations for all the examples using initial meshes with reasonable densities. Making use of Abaqus＇s rich pre/post- processing functionalities and powerful standard/explicit solvers, the developed method offers a practical tool for engineering analysts to model complex 3D fracture problems.

Mesoscale fracture behavior of Longmaxi outcrop shale with different bedding angles:Experimental and numerical investigations

The mechanical properties and fracturing mechanism of shale containing beddings are critically important in shale gas exploitation and wellbore stability.To investigate the effects of shale bedding on crack behavior and fracturing mechanism,scanning electron microscope(SEM)with a loading system was employed to carry out three-point bending tests on Longmaxi outcrop shale.The crack initiation and propagation of Longmaxi shale were observed and recorded by taking photos during loading.The cracking paths were extracted to calculate the crack length through a MATLAB program.The peak load,fracture toughness and fracture energy all increase with the bedding angle from 0°to 90°.The crack length and energy were also found to increase with the bedding angle in the range of 0°-600 and then drop slightly.The fracturing mechanism of shale includes the main crack affected by the bedding angle and disturbed by randomly distributed particles.The main cracking path was accompanied by several microcrack branches which could form an interconnected crack system.When the main crack encounters larger sedimentary particles,it will deflect around the particles and then restore to the initial direction.A numerical technique using extended finite element method(XFEM)coupled with anisotropic cohesive damage criteria was developed,which is able to capture the dependence of crack propagations on bedding angle and sedimentary particles.This study sheds light on understanding and predicting mesoscale fracture behavior of shale with different bedding angles.

Simulation Research on Stress Intensity Factors of Different Crack Aspect Ratios on Hollow Axles

Because of the wicked service environment of the high speed train, it is possible that the hollow axle of the train may encounter the foreign object damage and form a sharp notch. Under the fatigue loading a crack can initiate from the notch and propagate to failure. It is noted that the stress intensity factor is the control parameter of the crack propagating, for the purpose of getting the more exact propagation characteristics, the stress intensity factor is studied mainly. The service loads of hollow axles are defined, and the stress distribution of hollow axles is obtained according to the load spectrum. The semi-ellipse crack configuration is defined with three parameters： the aspect ratio, the relative depth and the relative location along the crack front. Quarter point 20-node isoparametric degenerate singular elements are used for the region near the crack tip. The finite element model of crack extension of hollow axle is created, and the crack front is dispersed which can realize orthogonal extension. Based on this the stress intensity factors of crack front were calculated, and the distribution rules of the stress intensity factors of different initial crack shapes are obtained. The conclusions are compared with that of the analytic method and they agree with each other very well, and the calculating results show that there is a close relationship between the stress intensity factor and the initial crack shape. For a round crack the stress intensity factor at the surface point increases faster than the one at the center point with the crock propagation. However, for a narrow crack, the results are in contrast with that of a round one. So, all the cracks with different shapes propagate toward to a similar shape, and they grow at this shape to end. The study may contribute to the crack propagate characteristics research.

Application of scaled boundary finite element method in static and dynamic fracture problems

The scaled boundary finite element method （SBFEM） is a recently developed numerical method combining advantages of both finite element methods （FEM） and boundary element methods （BEM） and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors （SIFs） directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.

A novel virtual node method for polygonal elements

A novel polygonal finite element method （PFEM） based on partition of unity is proposed, termed the virtual node method （VNM）. To test the performance of the present method, numerical examples are given for solid mechanics problems. With a polynomial form, the VNM achieves better results than those of traditional PFEMs, including the Wachspress method and the mean value method in standard patch tests. Compared with the standard triangular FEM, the VNM can achieve better accuracy. With the ability to construct shape functions on polygonal elements, the VNM provides greater flexibility in mesh generation. Therefore, several fracture problems are studied to demonstrate the potential implementation. With the advantage of the VNM, the convenient refinement and remeshing strategy are applied.

Hydration Process and Crack Tendency of Concrete Based on Resistivity and Restrained Shrinkage Crack

Hydration process, crack potential and setting time of concrete grade C30, C40 and C50 were monitored by using a non-contact electrical resistivity apparatus, a novel plastic ring mould and penetration resistance methods, respectively. The results show the highest resistivity of C30 at the early stage until a point when C50 accelerated and overtook the others. It has been experimentally confirmed that the crossing point of C30 and C50 corresponds to the final setting time of C50. From resistivity derivative curve, four different stages were observed upon which the hydration process is classified; these are dissolution, induction, acceleration and deceleration periods. Consequently, restrained shrinkage crack and setting time results demonstrated that C50 set and cracked the earliest. The cracking time of all the samples occurred within a reasonable experimental period thus the novel plastic ring is a convenient method for predicting concrete＇s crack potential. The highest inflection time（t_i） obtained from resistivity curve and the final setting time（t_f） were used with crack time（t_c） in coming up with mathematical models for the prediction of concrete＇s cracking age for the range of concrete grade considered. Finally, an ANSYS numerical simulation supports the experimental findings in terms of the earliest crack age of C50 and the crack location.

Numerical simulation of plate rigid restraint cracking tests based on cohesive element model

Cohesive element is developed from the Dugdal-Barenblatt model in the field of fracture mechanics. The mechanical constitutive relation of cohesive element can be artificially assumed depending on the specific applications. It has been successfully applied in the study of crystal plasticity/brittle fracture process and decohesion between delaminations. In this paper, tensile experiments of large steel plate with different length of pre-existing cracks are conducted. Based on commercial software ABAQUS, cohesive element is adopted to simulate the tensile tests, and appropriate parameter values are obtained by fitting displacement-force curves. Using these parameters, a numerical method is presented by applying cohesive element to thermo-elastic-plastic finite element method （TEP-FEM） to simulate plate rigid restraint cracking （PRRC） tests. By changing constitutive relation of cohesive element, dimensions of the model and welding conditions, the influence of welding restraint intensity and welding conditions on the crack propagation are discussed, respectively. Three types of welding cold cracking are simulated. Significant influence of welding cold cracking on resistant stress in welding line is captured by this numerical method.

Towards fully automatic modelling of the fracture process in quasi-brittle and ductile materials:a unified crack growth criterion

Fully automatic finite element(FE) modelling of the fracture process in quasi-brittle materials such as concrete and rocks and ductile materials such as metals and alloys,is of great significance in assessing structural integrity and presents tre-mendous challenges to the engineering community. One challenge lies in the adoption of an objective and effective crack propagation criterion. This paper proposes a crack propagation criterion based on the principle of energy conservation and the cohesive zone model(CZM) . The virtual crack extension technique is used to calculate the differential terms in the criterion. A fully-automatic discrete crack modelling methodology,integrating the developed criterion,the CZM to model the crack,a simple remeshing procedure to accommodate crack propagation,the J2 flow theory implemented within the incremental plasticity framework to model the ductile materials,and a local arc-length solver to the nonlinear equation system,is developed and im-plemented in an in-house program. Three examples,i.e.,a plain concrete beam with a single shear crack,a reinforced concrete(RC) beam with multiple cracks and a compact-tension steel specimen,are simulated. Good agreement between numerical predictions and experimental data is found,which demonstrates the applicability of the criterion to both quasi-brittle and ductile materials.