A new node-pairs contact algorithm is proposed to deal with a composite material or bi-material interface crack face contact and friction problem （e.g., resistant coating and thermal barrier coatings） subjected to c...A new node-pairs contact algorithm is proposed to deal with a composite material or bi-material interface crack face contact and friction problem （e.g., resistant coating and thermal barrier coatings） subjected to complicated load conditions. To decrease the calculation scale and calculation errors, the local Lagrange multipliers are solved only on a pair of contact nodes using the Jacobi iteration method, and the constraint modification of the tangential multipliers are required. After the calculation of the present node-pairs Lagrange multiplier, it is turned to next contact node-pairs until all node-pairs have finished. Compared with an ordinary contact algorithm, the new local node-pairs contact algorithm is allowed a more precise element on the contact face without the stiffness matrix singularity. The stress intensity factors （SIFs） and the contact region of an infinite plate central crack are calculated and show good agreement with those in the literature. The contact zone near the crack tip as well as its influence on singularity of stress fields are studied. Furthermore, the frictional contacts are also considered and found to have a significant influence on the SIFs. The normalized mode-II stress intensity factors KII for the friction coefficient decrease by 16% when f changes from 1 to 0.展开更多
地震破裂过程的数值模拟通常是假定破裂发生在先存断层上.但是,实际地震过程中,除了沿着先存断层产生自发破裂外,很有可能会产生新的破裂.为此,本文将传统有限元与扩展有限元方法(Extended Finite Element Method,XFEM)相结合,首先利用...地震破裂过程的数值模拟通常是假定破裂发生在先存断层上.但是,实际地震过程中,除了沿着先存断层产生自发破裂外,很有可能会产生新的破裂.为此,本文将传统有限元与扩展有限元方法(Extended Finite Element Method,XFEM)相结合,首先利用传统有限元方法模拟先存断层上的自发破裂过程,随后采用扩展有限元方法模拟不同介质(内摩擦角分别为0°、22.5°和45°)中,由于先存断层破裂造成的断层尖端应力集中导致的新生破裂过程.当内摩擦角为0°时,2条新生破裂与先存断层走向之间的夹角分别为0°和90°;当内摩擦角为22.5°时,该夹角变为11.25°和78.75°;当内摩擦角为45°时,该夹角变为22.5°和67.5°.在这三种情形下,新生的2条破裂,相互垂直,但破裂特征(破裂长度及破裂速度)不同,与先存断层之间的夹角越小的新生破裂,其扩展距离越远,扩展速度也更快(甚至为超剪切破裂).此外,在模拟结果中,岩石强度越低,断层扩展的长度越长,扩展速度也越快.本研究对于地震破裂全过程的数值模拟、深刻认识地震震源过程、断层生长以及地震灾害评估等有着重要的科学意义.展开更多
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.展开更多
Based on the existing experiment results, the fretting fatigue contact geometry of a riveted two aluminum specimen was studied using the finite element method. The contact stress fields of the inner and outer contact ...Based on the existing experiment results, the fretting fatigue contact geometry of a riveted two aluminum specimen was studied using the finite element method. The contact stress fields of the inner and outer contact edges on the two specimen’s up and down surface under different contact friction coefficient and the fatigue loads were analyzed, the influences of the contact friction coefficient and remote stress on crack initiation were discussed. The results were in well agreement with the observations of the existing experiments, that is, the crack initiated places of the first aluminum specimen change from the area of 900 to 450, and the crack initiated places of the second aluminum specimen change from the area of 900 to 1350 with the increase of the friction coefficient and the remote stress.展开更多
The basic equations of the representation of plasticity at an interface crack by an inclined strip yield superdislocation model are derived. With a special combination of the material properties and small-scale yieldi...The basic equations of the representation of plasticity at an interface crack by an inclined strip yield superdislocation model are derived. With a special combination of the material properties and small-scale yielding case, the problem is reduced to an algebraic equation in an unknown, the ratio of the plastic zone size. Some discussions on the limitations of this model are presented.展开更多
The seismic stability of a cracked dam was examined in this study. Geometric nonlinearity and large deformations, as well as the contact condition at the crack site, were taken into consideration. The location of pene...The seismic stability of a cracked dam was examined in this study. Geometric nonlinearity and large deformations, as well as the contact condition at the crack site, were taken into consideration. The location of penetrated cracks was first identified using the concrete plastic-damage model based on the nonlinear finite element method （FEM）. Then, the hard contact algorithm was used to simulate the crack interaction in the normal direction, and the Coloumb friction model was used to simulate the crack interaction in the tangential direction. After verification of numerical models through a case study, the seismic stability of the Koyna Dam with two types of penetrated cracks is discussed in detail with different seismic peak accelerations, and the collapse processes of the cracked dam are also presented. The results show that the stability of the dam with two types of penetrated cracks can be ensured in an earthquake with a magnitude of the original Koyna earthquake, and the cracked dam has a large earthquake-resistant margin. The failure processes of the cracked dam in strong earthquakes can be divided into two stages： the sliding stage and the overturning stage. The sliding stage ends near the peak acceleration, and the top block slides a long distance along the crack before the collapse occurs. The maximum sliding displacement of the top block will decrease with an increasing friction coefficient at the crack site.展开更多
For a compression-shear mixed mode interface crack, it is difficult to solve the stress and strain fields considering the material viscosity, the crack-tip singularity, the frictional effect, and the mixed loading lev...For a compression-shear mixed mode interface crack, it is difficult to solve the stress and strain fields considering the material viscosity, the crack-tip singularity, the frictional effect, and the mixed loading level. In this paper, a mechanical model of the dynamic propagation interface crack for the compression-shear mixed mode is proposed using an elastic-viscoplastic constitutive model. The governing equations of propagation crack interface at the crack-tip are given. The numerical analysis is performed for the interface crack of the compression-shear mixed mode by introducing a displacement function and some boundary conditions. The distributed regularities of stress field of the interface crack-tip are discussed with several special parameters. The final results show that the viscosity effect and the frictional contact effect on the crack surface and the mixed-load parameter are important factors in studying the mixed mode interface crack- tip fields. These fields are controlled by the viscosity coefficient, the Mach number, and the singularity exponent.展开更多
基金supported by the National Basic Research Program of China(Grant No.2012CB026200)the National Natural Science Foundation of China(Grant No.50878048)
文摘A new node-pairs contact algorithm is proposed to deal with a composite material or bi-material interface crack face contact and friction problem （e.g., resistant coating and thermal barrier coatings） subjected to complicated load conditions. To decrease the calculation scale and calculation errors, the local Lagrange multipliers are solved only on a pair of contact nodes using the Jacobi iteration method, and the constraint modification of the tangential multipliers are required. After the calculation of the present node-pairs Lagrange multiplier, it is turned to next contact node-pairs until all node-pairs have finished. Compared with an ordinary contact algorithm, the new local node-pairs contact algorithm is allowed a more precise element on the contact face without the stiffness matrix singularity. The stress intensity factors （SIFs） and the contact region of an infinite plate central crack are calculated and show good agreement with those in the literature. The contact zone near the crack tip as well as its influence on singularity of stress fields are studied. Furthermore, the frictional contacts are also considered and found to have a significant influence on the SIFs. The normalized mode-II stress intensity factors KII for the friction coefficient decrease by 16% when f changes from 1 to 0.
文摘地震破裂过程的数值模拟通常是假定破裂发生在先存断层上.但是,实际地震过程中,除了沿着先存断层产生自发破裂外,很有可能会产生新的破裂.为此,本文将传统有限元与扩展有限元方法(Extended Finite Element Method,XFEM)相结合,首先利用传统有限元方法模拟先存断层上的自发破裂过程,随后采用扩展有限元方法模拟不同介质(内摩擦角分别为0°、22.5°和45°)中,由于先存断层破裂造成的断层尖端应力集中导致的新生破裂过程.当内摩擦角为0°时,2条新生破裂与先存断层走向之间的夹角分别为0°和90°;当内摩擦角为22.5°时,该夹角变为11.25°和78.75°;当内摩擦角为45°时,该夹角变为22.5°和67.5°.在这三种情形下,新生的2条破裂,相互垂直,但破裂特征(破裂长度及破裂速度)不同,与先存断层之间的夹角越小的新生破裂,其扩展距离越远,扩展速度也更快(甚至为超剪切破裂).此外,在模拟结果中,岩石强度越低,断层扩展的长度越长,扩展速度也越快.本研究对于地震破裂全过程的数值模拟、深刻认识地震震源过程、断层生长以及地震灾害评估等有着重要的科学意义.
基金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.
文摘Based on the existing experiment results, the fretting fatigue contact geometry of a riveted two aluminum specimen was studied using the finite element method. The contact stress fields of the inner and outer contact edges on the two specimen’s up and down surface under different contact friction coefficient and the fatigue loads were analyzed, the influences of the contact friction coefficient and remote stress on crack initiation were discussed. The results were in well agreement with the observations of the existing experiments, that is, the crack initiated places of the first aluminum specimen change from the area of 900 to 450, and the crack initiated places of the second aluminum specimen change from the area of 900 to 1350 with the increase of the friction coefficient and the remote stress.
基金Visiting Research Fellow in Institute of Industrial Science.University of Tokyo,supported by Japan Society for the Promotion of Science.
文摘The basic equations of the representation of plasticity at an interface crack by an inclined strip yield superdislocation model are derived. With a special combination of the material properties and small-scale yielding case, the problem is reduced to an algebraic equation in an unknown, the ratio of the plastic zone size. Some discussions on the limitations of this model are presented.
基金supported by the National Basic Research Program of China (973 Program, Grant No. 2007CB714104)the National Natural Science Foundation of China (Grant No. 50779011)the Innovative Project for Graduate Students of Jiangsu Province (Grant No. CX10B_202Z)
文摘The seismic stability of a cracked dam was examined in this study. Geometric nonlinearity and large deformations, as well as the contact condition at the crack site, were taken into consideration. The location of penetrated cracks was first identified using the concrete plastic-damage model based on the nonlinear finite element method （FEM）. Then, the hard contact algorithm was used to simulate the crack interaction in the normal direction, and the Coloumb friction model was used to simulate the crack interaction in the tangential direction. After verification of numerical models through a case study, the seismic stability of the Koyna Dam with two types of penetrated cracks is discussed in detail with different seismic peak accelerations, and the collapse processes of the cracked dam are also presented. The results show that the stability of the dam with two types of penetrated cracks can be ensured in an earthquake with a magnitude of the original Koyna earthquake, and the cracked dam has a large earthquake-resistant margin. The failure processes of the cracked dam in strong earthquakes can be divided into two stages： the sliding stage and the overturning stage. The sliding stage ends near the peak acceleration, and the top block slides a long distance along the crack before the collapse occurs. The maximum sliding displacement of the top block will decrease with an increasing friction coefficient at the crack site.
基金Project supported by the National Natural Science Foundation of China(No.11302054)the Fundamental Research Funds for the Central Universities(No.HEUCF130216)
文摘For a compression-shear mixed mode interface crack, it is difficult to solve the stress and strain fields considering the material viscosity, the crack-tip singularity, the frictional effect, and the mixed loading level. In this paper, a mechanical model of the dynamic propagation interface crack for the compression-shear mixed mode is proposed using an elastic-viscoplastic constitutive model. The governing equations of propagation crack interface at the crack-tip are given. The numerical analysis is performed for the interface crack of the compression-shear mixed mode by introducing a displacement function and some boundary conditions. The distributed regularities of stress field of the interface crack-tip are discussed with several special parameters. The final results show that the viscosity effect and the frictional contact effect on the crack surface and the mixed-load parameter are important factors in studying the mixed mode interface crack- tip fields. These fields are controlled by the viscosity coefficient, the Mach number, and the singularity exponent.
