Taylor stochastic finite element method (SFEM) is applied to analyze the uncertainty of plane multiple cracks stress intensity factors (SIFs) considering the uncertainties of material properties, crack length, and loa...Taylor stochastic finite element method (SFEM) is applied to analyze the uncertainty of plane multiple cracks stress intensity factors (SIFs) considering the uncertainties of material properties, crack length, and load. The stochastic finite element model of plane multiple cracks are presented. In this model, crack tips are meshed with six-node triangular quarter-point elements; and other area is meshed with six-node triangular elements. The partial derivatives of displacement and stiffness matrix with resp...展开更多
Effects of alloy elements on the microstructure and crack resistance of Fe-C-Cr weld surfacing layer were investigated. The results show that microstructures of the layer mainly consist of carbides and austenite matri...Effects of alloy elements on the microstructure and crack resistance of Fe-C-Cr weld surfacing layer were investigated. The results show that microstructures of the layer mainly consist of carbides and austenite matrix. Increasing C and Cr contents impair the crack resistance of the layer due to increased amount of brittle carbides. The addition of Ni, Nb or Mo improves the crack resistance of Fe-C-Cr weld surfacing layer by increasing the amount of austenite and forming fine NbC or M 7C 3 carbides in the layer. But, the excessive Nb (>2.50wt%) or Mo (>1.88wt%) impairs the crack resistance of the layer, which has relation with increased carbides or carbide coarsening and austenite matrix solid solution strengthening. The proper combination of C, Cr, Ni, Nb and Mo can further improve not only the crack resistance of Fe-C-Cr weld surfacing layer but also the erosion resistance as a result of fine NbC and M 7C 3 carbides distributing uniformly in austenite matrix. The optimal layer compositions are 3.05wt%C, 20.58wt%Cr, 1.75wt%Ni, 2.00wt%Nb and 1.88wt%Mo.展开更多
Rare earth elements（REE） are a set of seventeen chemical elements in the periodic table, which have been found in a variety classes of minerals in the world. The global demand for REE and their compounds has increas...Rare earth elements（REE） are a set of seventeen chemical elements in the periodic table, which have been found in a variety classes of minerals in the world. The global demand for REE and their compounds has increased in high-tech applications over the past decades due to their chemical, catalytic, electrical, magnetic, and optical properties. After mining, comminution and ore beneficiation, hydrometallurgical and pyrometallurgical processes are carried out to recover REE compounds. The concentrate cracking is a process wherein the structure of a REE mineral is modified in order to dissolve REE in a weak acid solution. This would be accomplished through various techniques such as alkaline cracking and acid baking. The cracked concentrate is then leached to dissolve REE in a solution and subsequently recover those via neutralization, precipitation or solvent extraction methods. This paper reviewed the cracking and leaching methods applied to REE concentrates to produce an intermediate product for further REE separation.展开更多
Self-deformation cracking is the cracking caused by thermal deformation, autogenous volume deformation or shrinkage deformation. In this paper, an extended finite element calculation method was deduced for concrete cr...Self-deformation cracking is the cracking caused by thermal deformation, autogenous volume deformation or shrinkage deformation. In this paper, an extended finite element calculation method was deduced for concrete crack propagation under a constant hydration and hardening condition during the construction period, and a corresponding programming code was developed. The experimental investigation shows that initial crack propagation caused by self-deformation loads can be analyzed by this program. This improved algorithm was a preliminary application of the XFEM to the problem of the concrete self-deformation cracking during the hydration and hardening period. However, room for improvement exists for this algorithm in terms of matching calculation programs with mass concrete temperature fields containing cooling pipes and the influence of creep or damage on crack propagation.展开更多
The effect of rare earth elements on the thermal cracking resistance of high speed steel （HSS） rolls was investigated. Laser rapid heating was used for thermal fatigue experiments. Thermal cracks and microstructure ...The effect of rare earth elements on the thermal cracking resistance of high speed steel （HSS） rolls was investigated. Laser rapid heating was used for thermal fatigue experiments. Thermal cracks and microstructure were observed using metalloscopy and scanning electron microscopy. The results showed that thermal cracks initiated from the interface between the matrix and eutectic carbides （including M6C and M7C3 type carbides）,and propagated along the interface between the two phases. MC type carbides enriched with vanadium could prevent the propagation of thermal cracks. The presence of rare earth elements decreased the quantity of big eutectic carbides,and proportionally increased spherical and rod-shaped MC type carbide content. HSS0 （0.00% RE） had approximately three times the thermal cracking density of HSS3 （0.12 wt.% RE）. Rare earth elements were shown to significantly improve the microstructure and thermal cracking resistance of HSS rolls.展开更多
A three-dimensional thermo-mechanical coupled finite element model is built up to simulate the phenomena of dynamical contact and frictional heating of crack faces when the plate containing the crack is excited by hig...A three-dimensional thermo-mechanical coupled finite element model is built up to simulate the phenomena of dynamical contact and frictional heating of crack faces when the plate containing the crack is excited by high-intensity ultrasonic pulses. In the finite element model, the high-power ultrasonic transducer is modeled by using a piezoelectric thermal-analogy method, and the dynamical interaction between both crack faces is modeled using a contact-impact theory. In the simulations, the frictional heating taking place at the crack faces is quantitatively calculated by using finite element thermal-structural coupling analysis, especially, the influences of acoustic chaos to plate vibration and crack heating are calculated and analysed in detail. Meanwhile, the related ultrasonic infrared images are also obtained experimentally, and the theoretical simulation results are in agreement with that of the experiments. The results show that, by using the theoretical method, a good simulation of dynamic interaction and friction heating process of the crack faces under non-chaotic or chaotic sound excitation can be obtained.展开更多
The formation and propagation of the popular off-corner subsurface cracks in bloom continuous casting were investigated through thermo-mechanical analysis using three coupled thermo-mechanical models. A two-dimensiona...The formation and propagation of the popular off-corner subsurface cracks in bloom continuous casting were investigated through thermo-mechanical analysis using three coupled thermo-mechanical models. A two-dimensional thermo-elasto-visco-plastic finite element model was developed to predict the mould gap evolution, temperature profiles and deformation behavior of the solidified shell in the mould region. Then, a three-dimensional model was adopted to calculate the shell growth, tempera- ture history and the development of stresses and strains of the shell in the following secondary cooling zones. Finally, another three-dimensional model was used to analyze the stress distributions in the straightening region, The results showed that the off-corner cracks in the shell originated from the mould owing to the tensile strain developed in the crack sensitive regions of the solidification front, and they could be driven deeper by the possible severe surface temperature rebound and the extensive tensile stress in the secondary cooling zone, especially upon the straightening operation of the bloom casting. It is revealed that more homogenous shell temperature and thickness can be obtained through optimization of mould corner radius, casting speed and secondary cooling scheme, which help to decrease stress and strain concentration and therefore prevent the initiation of the cracks.展开更多
The Finite Element Method of Lines （FEMOL） is a semi-analytic approach and takes a position between FEM and analytic methods. First, FEMOL in Fracture Mechanics is presented in detail. Then, the method is applied to...The Finite Element Method of Lines （FEMOL） is a semi-analytic approach and takes a position between FEM and analytic methods. First, FEMOL in Fracture Mechanics is presented in detail. Then, the method is applied to a set of examples such as edge-crack plate, the central-crack plate, the plate with cracks emanating from a hole under tensile or under combination loads of tensile and bending. Their dimensionless stress distribution, the stress intensify factor （SIF） and crack opening displacement （COD） are obtained, and comparison with known solutions by other methods are reported. It is found that a good accuracy is achieved by FEMOL. The method is successfully modified to remarkably increase the accuracy and reduce convergence difficulties. So it is a very useful and new tool in studying fracture mechanics problems.展开更多
The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good ac...The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good accuracy can be obtained with relatively coarse girds. In particular, application to the tension specimen shows very good agreement with the evaluation of stress intensity factors, which is better than the results of other methods. This implies a considerable potential for using this method in the 3D analysis of finite geometry solids and suggests a possible extension of this technique to nonlinear material behavior.展开更多
In this paper, a method is developed for determining the effective stiffness of the cracked component. The stiffness matrix of the cracked component is integrated into the global stiffness matrix of the finite element...In this paper, a method is developed for determining the effective stiffness of the cracked component. The stiffness matrix of the cracked component is integrated into the global stiffness matrix of the finite element model of the global platform for the FE calculation of the structure in any environmental conditions. The stiffness matrix equation of the cracked component is derived by use of the finite variation principle and fracture mechanics. The equivalent parameters defining the element that simulates the cracked component are mathematically presented, and can be easily used for the FE calculation of large scale cracked structures together with any finite element program. The theories developed are validated by both lab tests and numerical calculations, and applied to the evaluation of crack effect on the strength of a fixed platform and a self-elevating drilling rig.展开更多
Based on the crack tip field expansion of the Reissner plate, a special high order bending crack tip element is developed, and the element stiffness matrix is given in the explicit form, which is especially convenien...Based on the crack tip field expansion of the Reissner plate, a special high order bending crack tip element is developed, and the element stiffness matrix is given in the explicit form, which is especially convenient for engineering analyses. A numerical example is presented and compared with previous results to demonstrate the efficiency and accuracy of the special element.展开更多
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 ...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.展开更多
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 ...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.展开更多
A simple and effective boundary element method for stress intensity factor calculation for crack problems in a plane elastic plate is presented. The boundary element method consists of the constant displacement discon...A simple and effective boundary element method for stress intensity factor calculation for crack problems in a plane elastic plate is presented. The boundary element method consists of the constant displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity elements proposed by YAN Xiangqiao. In the boundary element implementation the left or the right crack-tip displacement discontinuity element was placed locally at the corresponding left or right each crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. Test examples (i. e. , a center crack in an infinite plate under tension, a circular hole and a crack in an infinite plate under tension) are included to illustrate that the numerical approach is very simple and accurate for stress intensity factor calculation of plane elasticity crack problems. In addition, specifically, the stress intensity factors of branching cracks emanating from a square hole in a rectangular plate under biaxial loads were analysed. These numerical results indicate the present numerical approach is very effective for calculating stress intensity factors of complex cracks in a 2-D finite body, and are used to reveal the effect of the biaxial loads and the cracked body geometry on stress intensity factors.展开更多
It has been postulated that, with tensile loading conditions, micro-cracks onthin hard film act as stress concentrators enhancing plastic deformation of the substrate materialin their vicinity. Under favorable conditi...It has been postulated that, with tensile loading conditions, micro-cracks onthin hard film act as stress concentrators enhancing plastic deformation of the substrate materialin their vicinity. Under favorable conditions the localized plastic flow near the cracks may turninto macroscopic plastic strain thus affects the plasticity behaviors of the substrate. Thisphenomenon is analyzed quantitatively with finite element method with special attention focused onthe analysis and discussion of the effects of plastic work hardening rate, film thickness and crackdepth on maximum plastic strain, critical loading stress and the size of the local plasticdeformation zone. Results show that micro-cracks on thin hard film have unnegligible effects on theplasticity behaviors of the substrate material under tensile loading.展开更多
The Saint-Venant torsion problems of a cylinder with curvilinear cracks were considered and reduced to solving the boundary integral equations only on cracks. Using the interpolation models for both singular crack tip...The Saint-Venant torsion problems of a cylinder with curvilinear cracks were considered and reduced to solving the boundary integral equations only on cracks. Using the interpolation models for both singular crack tip elements and other crack linear elements, the boundary element formulas of the torsion rigidity and stress intensity factors were given. Some typical torsion problems of a cylinder involving a straight, kinked or curvilinear crack were calculated. The obtained results for the case of straight crack agree well with those given by using the Gauss-Chebyshev integration formulas, which demonstrates the validity and applicability of the present boundary element method.展开更多
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 co...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.展开更多
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 applica...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.展开更多
Spherical indentation of ceramic coatings with metallic interlayer was performed by means of axisymmetric finite element analysis （FEA）. Two typical ceramic coatings with relatively high and low elastic modulus depo...Spherical indentation of ceramic coatings with metallic interlayer was performed by means of axisymmetric finite element analysis （FEA）. Two typical ceramic coatings with relatively high and low elastic modulus deposited on aluminum alloy and carbon steel were considered. The fracture mechanics of the ceramic coatings mechanisms due to occurrence of surface ring cracks extending traverse the coating thickness under spherical indentation are investigated within the framework of linear fracture mechanics. The J-integral associated to such cracks was computed. The evolution of J-integral vs the crack length and the indentation depth was studied. The effects of the interlayer, the coating and the substrate on the J-integral evolution were discussed. The results show that a suitable metallic interlayer can improve the fracture resistance of the coating systems under the same indentation conditions through reducing the J-integral.展开更多
The extended finite element method is used to analyze a plate with two parallel edge cracks impacted by a cylindrical projectile. The influence of the impact speed, crack length, plate thickness and notch tip radius o...The extended finite element method is used to analyze a plate with two parallel edge cracks impacted by a cylindrical projectile. The influence of the impact speed, crack length, plate thickness and notch tip radius on the crack initiation and propagation is studied. Dynamics equations are solved by an implicit time integration scheme which is unconditionally stable. Very good agreement is achieved between numerical predictions and experimental results. The critical velocity of the crack initiation under different conditions is examined. The influence of the crack length is greater than that of the impact speed, plate thickness and notch tip radius.展开更多
基金National Natural Science Foundation of China (10577015)Chinese Aeronautics Foundation (2006ZD53050, 03B53008)
文摘Taylor stochastic finite element method (SFEM) is applied to analyze the uncertainty of plane multiple cracks stress intensity factors (SIFs) considering the uncertainties of material properties, crack length, and load. The stochastic finite element model of plane multiple cracks are presented. In this model, crack tips are meshed with six-node triangular quarter-point elements; and other area is meshed with six-node triangular elements. The partial derivatives of displacement and stiffness matrix with resp...
文摘Effects of alloy elements on the microstructure and crack resistance of Fe-C-Cr weld surfacing layer were investigated. The results show that microstructures of the layer mainly consist of carbides and austenite matrix. Increasing C and Cr contents impair the crack resistance of the layer due to increased amount of brittle carbides. The addition of Ni, Nb or Mo improves the crack resistance of Fe-C-Cr weld surfacing layer by increasing the amount of austenite and forming fine NbC or M 7C 3 carbides in the layer. But, the excessive Nb (>2.50wt%) or Mo (>1.88wt%) impairs the crack resistance of the layer, which has relation with increased carbides or carbide coarsening and austenite matrix solid solution strengthening. The proper combination of C, Cr, Ni, Nb and Mo can further improve not only the crack resistance of Fe-C-Cr weld surfacing layer but also the erosion resistance as a result of fine NbC and M 7C 3 carbides distributing uniformly in austenite matrix. The optimal layer compositions are 3.05wt%C, 20.58wt%Cr, 1.75wt%Ni, 2.00wt%Nb and 1.88wt%Mo.
基金supported by the Research Initiation Grant of Queen’s University,Canada(#378140)Ontario Centers of Excellence VIP fund(#25394)
文摘Rare earth elements（REE） are a set of seventeen chemical elements in the periodic table, which have been found in a variety classes of minerals in the world. The global demand for REE and their compounds has increased in high-tech applications over the past decades due to their chemical, catalytic, electrical, magnetic, and optical properties. After mining, comminution and ore beneficiation, hydrometallurgical and pyrometallurgical processes are carried out to recover REE compounds. The concentrate cracking is a process wherein the structure of a REE mineral is modified in order to dissolve REE in a weak acid solution. This would be accomplished through various techniques such as alkaline cracking and acid baking. The cracked concentrate is then leached to dissolve REE in a solution and subsequently recover those via neutralization, precipitation or solvent extraction methods. This paper reviewed the cracking and leaching methods applied to REE concentrates to produce an intermediate product for further REE separation.
基金Funded by the National Natural Science Foundation of China(Nos.51409264,51509020,51209219)Zhejiang Provincial Natural Science Foundation(No.LQ15E090003)
文摘Self-deformation cracking is the cracking caused by thermal deformation, autogenous volume deformation or shrinkage deformation. In this paper, an extended finite element calculation method was deduced for concrete crack propagation under a constant hydration and hardening condition during the construction period, and a corresponding programming code was developed. The experimental investigation shows that initial crack propagation caused by self-deformation loads can be analyzed by this program. This improved algorithm was a preliminary application of the XFEM to the problem of the concrete self-deformation cracking during the hydration and hardening period. However, room for improvement exists for this algorithm in terms of matching calculation programs with mass concrete temperature fields containing cooling pipes and the influence of creep or damage on crack propagation.
基金Project supported by the National Natural Science Foundation of China (51071136)
文摘The effect of rare earth elements on the thermal cracking resistance of high speed steel （HSS） rolls was investigated. Laser rapid heating was used for thermal fatigue experiments. Thermal cracks and microstructure were observed using metalloscopy and scanning electron microscopy. The results showed that thermal cracks initiated from the interface between the matrix and eutectic carbides （including M6C and M7C3 type carbides）,and propagated along the interface between the two phases. MC type carbides enriched with vanadium could prevent the propagation of thermal cracks. The presence of rare earth elements decreased the quantity of big eutectic carbides,and proportionally increased spherical and rod-shaped MC type carbide content. HSS0 （0.00% RE） had approximately three times the thermal cracking density of HSS3 （0.12 wt.% RE）. Rare earth elements were shown to significantly improve the microstructure and thermal cracking resistance of HSS rolls.
基金Project supported by the National Natural Science Foundation of China (Grant No. 10574073)
文摘A three-dimensional thermo-mechanical coupled finite element model is built up to simulate the phenomena of dynamical contact and frictional heating of crack faces when the plate containing the crack is excited by high-intensity ultrasonic pulses. In the finite element model, the high-power ultrasonic transducer is modeled by using a piezoelectric thermal-analogy method, and the dynamical interaction between both crack faces is modeled using a contact-impact theory. In the simulations, the frictional heating taking place at the crack faces is quantitatively calculated by using finite element thermal-structural coupling analysis, especially, the influences of acoustic chaos to plate vibration and crack heating are calculated and analysed in detail. Meanwhile, the related ultrasonic infrared images are also obtained experimentally, and the theoretical simulation results are in agreement with that of the experiments. The results show that, by using the theoretical method, a good simulation of dynamic interaction and friction heating process of the crack faces under non-chaotic or chaotic sound excitation can be obtained.
文摘The formation and propagation of the popular off-corner subsurface cracks in bloom continuous casting were investigated through thermo-mechanical analysis using three coupled thermo-mechanical models. A two-dimensional thermo-elasto-visco-plastic finite element model was developed to predict the mould gap evolution, temperature profiles and deformation behavior of the solidified shell in the mould region. Then, a three-dimensional model was adopted to calculate the shell growth, tempera- ture history and the development of stresses and strains of the shell in the following secondary cooling zones. Finally, another three-dimensional model was used to analyze the stress distributions in the straightening region, The results showed that the off-corner cracks in the shell originated from the mould owing to the tensile strain developed in the crack sensitive regions of the solidification front, and they could be driven deeper by the possible severe surface temperature rebound and the extensive tensile stress in the secondary cooling zone, especially upon the straightening operation of the bloom casting. It is revealed that more homogenous shell temperature and thickness can be obtained through optimization of mould corner radius, casting speed and secondary cooling scheme, which help to decrease stress and strain concentration and therefore prevent the initiation of the cracks.
文摘The Finite Element Method of Lines （FEMOL） is a semi-analytic approach and takes a position between FEM and analytic methods. First, FEMOL in Fracture Mechanics is presented in detail. Then, the method is applied to a set of examples such as edge-crack plate, the central-crack plate, the plate with cracks emanating from a hole under tensile or under combination loads of tensile and bending. Their dimensionless stress distribution, the stress intensify factor （SIF） and crack opening displacement （COD） are obtained, and comparison with known solutions by other methods are reported. It is found that a good accuracy is achieved by FEMOL. The method is successfully modified to remarkably increase the accuracy and reduce convergence difficulties. So it is a very useful and new tool in studying fracture mechanics problems.
文摘The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good accuracy can be obtained with relatively coarse girds. In particular, application to the tension specimen shows very good agreement with the evaluation of stress intensity factors, which is better than the results of other methods. This implies a considerable potential for using this method in the 3D analysis of finite geometry solids and suggests a possible extension of this technique to nonlinear material behavior.
文摘In this paper, a method is developed for determining the effective stiffness of the cracked component. The stiffness matrix of the cracked component is integrated into the global stiffness matrix of the finite element model of the global platform for the FE calculation of the structure in any environmental conditions. The stiffness matrix equation of the cracked component is derived by use of the finite variation principle and fracture mechanics. The equivalent parameters defining the element that simulates the cracked component are mathematically presented, and can be easily used for the FE calculation of large scale cracked structures together with any finite element program. The theories developed are validated by both lab tests and numerical calculations, and applied to the evaluation of crack effect on the strength of a fixed platform and a self-elevating drilling rig.
文摘Based on the crack tip field expansion of the Reissner plate, a special high order bending crack tip element is developed, and the element stiffness matrix is given in the explicit form, which is especially convenient for engineering analyses. A numerical example is presented and compared with previous results to demonstrate the efficiency and accuracy of the special element.
基金supported by EPSRC UK(No.EP/F00656X/1)Xiangting Su's one-year visit to the University of Liverpoosupported by the China Scholarship Council and the National Natural Science Foundation of China(No.50579081).
文摘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.
文摘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.
基金Project supported by the National Natural Science Foundation of China (No. 10272037)
文摘A simple and effective boundary element method for stress intensity factor calculation for crack problems in a plane elastic plate is presented. The boundary element method consists of the constant displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity elements proposed by YAN Xiangqiao. In the boundary element implementation the left or the right crack-tip displacement discontinuity element was placed locally at the corresponding left or right each crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. Test examples (i. e. , a center crack in an infinite plate under tension, a circular hole and a crack in an infinite plate under tension) are included to illustrate that the numerical approach is very simple and accurate for stress intensity factor calculation of plane elasticity crack problems. In addition, specifically, the stress intensity factors of branching cracks emanating from a square hole in a rectangular plate under biaxial loads were analysed. These numerical results indicate the present numerical approach is very effective for calculating stress intensity factors of complex cracks in a 2-D finite body, and are used to reveal the effect of the biaxial loads and the cracked body geometry on stress intensity factors.
基金This project is supported by National Natural Science Foundation of China(No.59705009).
文摘It has been postulated that, with tensile loading conditions, micro-cracks onthin hard film act as stress concentrators enhancing plastic deformation of the substrate materialin their vicinity. Under favorable conditions the localized plastic flow near the cracks may turninto macroscopic plastic strain thus affects the plasticity behaviors of the substrate. Thisphenomenon is analyzed quantitatively with finite element method with special attention focused onthe analysis and discussion of the effects of plastic work hardening rate, film thickness and crackdepth on maximum plastic strain, critical loading stress and the size of the local plasticdeformation zone. Results show that micro-cracks on thin hard film have unnegligible effects on theplasticity behaviors of the substrate material under tensile loading.
文摘The Saint-Venant torsion problems of a cylinder with curvilinear cracks were considered and reduced to solving the boundary integral equations only on cracks. Using the interpolation models for both singular crack tip elements and other crack linear elements, the boundary element formulas of the torsion rigidity and stress intensity factors were given. Some typical torsion problems of a cylinder involving a straight, kinked or curvilinear crack were calculated. The obtained results for the case of straight crack agree well with those given by using the Gauss-Chebyshev integration formulas, which demonstrates the validity and applicability of the present boundary element method.
文摘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.
文摘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.
文摘Spherical indentation of ceramic coatings with metallic interlayer was performed by means of axisymmetric finite element analysis （FEA）. Two typical ceramic coatings with relatively high and low elastic modulus deposited on aluminum alloy and carbon steel were considered. The fracture mechanics of the ceramic coatings mechanisms due to occurrence of surface ring cracks extending traverse the coating thickness under spherical indentation are investigated within the framework of linear fracture mechanics. The J-integral associated to such cracks was computed. The evolution of J-integral vs the crack length and the indentation depth was studied. The effects of the interlayer, the coating and the substrate on the J-integral evolution were discussed. The results show that a suitable metallic interlayer can improve the fracture resistance of the coating systems under the same indentation conditions through reducing the J-integral.
基金supported by the National Natural Science Foundation of China(Nos.11272096 and 11472086)the Research Fund for the Doctoral Program of Higher Education of China(No.20112304110015)
文摘The extended finite element method is used to analyze a plate with two parallel edge cracks impacted by a cylindrical projectile. The influence of the impact speed, crack length, plate thickness and notch tip radius on the crack initiation and propagation is studied. Dynamics equations are solved by an implicit time integration scheme which is unconditionally stable. Very good agreement is achieved between numerical predictions and experimental results. The critical velocity of the crack initiation under different conditions is examined. The influence of the crack length is greater than that of the impact speed, plate thickness and notch tip radius.
