基于ABAQUS平台考虑T应力的Ⅰ型裂纹扩展模拟开发

工程结构在制造工艺过程中或使用期间会产生裂纹,对结构断裂路径的预测和研究是防治工程安全问题发生的重要手段。在考虑裂纹尖端应力场常数项T应力的基础上对传统的最大周向应力准则(Maximum tangential stress criterion,MTS)和最小应变能密度因子准则(Minimum strain energy density criterion,SED)进行修正,采用Python语言对ABAQUS的前、后处理和有限元计算模块进行二次开发,通过计算最优解的粒子群算法(Particle swarm optimization,PSO)将修正后的准则编入裂纹自动扩展程序脚本中。利用上述二次开发程序对初始纯Ⅰ型裂纹的扩展路径进行模拟,结果表明:采用ABAQUS脚本程序模拟结果与相关文献实验结果吻合,表明了程序的有效性,进而实现考虑T应力的多种断裂准则对裂纹扩展路径的预测;当T应力值处于一定范围内时,修正的MTS准则无法预测裂纹发生的偏转现象,扩展路径呈直线,此时可采用修正的SED准则进行预测。

基于辛几何方法测定Ⅰ型应力强度因子

为得到高精度的Ⅰ型应力强度因子,提出了一种基于辛几何方法测定裂纹尖端Ⅰ型应力强度因子的方法。通过在裂纹尖端构建扇形区域,分离辛几何方法下裂纹尖端的全状态向量,结合裂纹尖端的位移场,计算出裂纹的应力强度因子。利用Ⅰ型裂纹尖端位移场的理论解,研究了扇形区域的特征点个数、圆心角大小以及半径长度对应力强度因子精确度的影响。研究结果表明:特征点个数对该方法求解Ⅰ型应力强度因子的精度的影响较大,特征点个数大于37时,相对误差稳定趋向于0;扇形圆心角对应力强度因子的精度影响较小,但是相同精度要求下,不同的扇形圆心角需要的特征点个数不同;扇形半径过小时,由于断裂过程中裂纹尖端塑性区域的影响,应力强度因子的解不稳定且相对误差较大。对比相对误差可以发现该方法计算得到的Ⅰ型应力强度因子的精度较高,且该方法只需得到裂纹尖端的位移场就可以精确得到应力强度因子,解决了以往求解体系在复杂应力等情况下求解步骤复杂,导致无法精确得到Ⅰ型应力强度因子的难题。

浅析OBE理念下的《汽车构造Ⅰ》一流课程建设

为体现以学生为中心、产出导向和持续改进教育理念以及课程的高阶性、创新性以及挑战度,文章以《汽车构造Ⅰ》课程为例,按照线上线下混合式一流课程的建设标准,遵循OBE理念,加强课程总体设计,融入课程思政元素,整合优化教学内容,设计“AEEA”的三阶线上线下混合式教学模式,开展多种教学方法的应用,建立过程化、多元化、混合式的课程考核评价体系,以期为高校专任教师开展一流课程建设与实施应用提供思路。

Boundary Element Analysis forModeⅢCrack Problems of Thin-Walled Structures from Micro-to Nano-Scales

This paper develops a new numerical framework for modeⅢcrack problems of thin-walled structures by integrating multiple advanced techniques in the boundary element literature.The details of special crack-tip elements for displacement and stress are derived.An exponential transformation technique is introduced to accurately calculate the nearly singular integral,which is the key task of the boundary element simulation of thin-walled structures.Three numerical experiments with different types of cracks are provided to verify the performance of the present numerical framework.Numerical results demonstrate that the present scheme is valid for modeⅢcrack problems of thin-walled structures with the thickness-to-length ratio in the microscale,even nanoscale,regime.

2D-SiC_(f)/SiC复合材料层间Ⅰ型断裂试验及表征

二维编织碳化硅纤维增强碳化硅复合材料(2D-SiC_(f)/SiC)在航空领域中得到广泛使用,然而该材料层间强度低,使其易萌生层间裂纹,引起分层破坏。为此,本工作采用楔形双悬臂梁法(W-DCB)和悬臂梁法(DCB)开展了层间Ⅰ型断裂试验,获得了2D-SiC_(f)/SiC的层间裂纹驱动的加载数据,得到了其裂纹端口张开力及张开位移变形曲线。在试验加载过程,通过光学显微镜监测了视觉裂纹扩展过程,探究了2D-SiC_(f)/SiC的层间I型裂纹扩展规律。结合理论分析和裂纹视觉特征解释了加载曲线拐点及其他特征点的断裂力学含义。利用扫描电子显微镜分析了2D-SiC_(f)/SiC的层间断面特征,揭示了断面分层裂纹扩展机制。结果表明:W-DCB方法测量的2D-SiC_(f)/SiC层间Ⅰ型初始能量释放率与DCB方法等效;2D-SiC_(f)/SiC层间Ⅰ型断裂过程中,裂纹端口变形曲线的多峰性不符合经典线弹性断裂力学预测的加载峰后特征,反映了2D-SiC_(f)/SiC层间约束关系的复杂性;层间断面为结构性非完全损伤,发生了局部纤维桥连现象。

压载荷下Ⅰ-Ⅱ复合型中心裂纹板弹塑性屈曲载荷研究

为了研究含裂纹结构承载能力,采用有限元方法建立压载荷下Ⅰ-Ⅱ复合型中心裂纹板模型,计算裂纹板弹塑性屈曲载荷,讨论相对裂纹长度(a/w)、裂纹倾角(θ)、板宽(w)、板厚(t)、边界条件对弹塑性屈曲载荷及其变化规律的影响。对于上下夹持中心裂纹板,随着裂纹长度和裂纹倾角的增大,板宽越大,中心裂纹板屈曲载荷对裂纹越敏感,中心裂纹板承载能力受到裂纹长度和裂纹倾角的影响越明显。弹塑性屈曲载荷随着板厚的增大而增大,并且随着裂纹倾角的增大,板厚越大,屈曲载荷的下降程度越显著。上下简支中心裂纹板与上下夹持中心裂纹板的屈曲载荷变化规律基本一致,且上下夹持中心裂纹板的屈曲载荷显著大于上下简支中心裂纹板的屈曲载荷。给出不同边界条件下弹塑性屈曲载荷与相对裂纹长度、裂纹倾角的回归关系式。通过中心裂纹板弹塑性屈曲试验发现,采用压缩力学性能计算弹塑性屈曲载荷更接近试验结果,试验屈曲载荷的变化规律与有限元模拟计算结果一致。研究结果对含裂纹结构的承载能力评价具有参考价值。

冲击荷载下充填节理岩体Ⅰ型裂纹动态扩展特性研究

为研究冲击荷载作用下充填节理岩体Ⅰ型裂纹动态扩展特性,利用分离式霍普金森杆(split Hopkinson pressure bar,简称SHPB)试验系统及裂纹扩展测试系统,对石膏、水泥及云石胶3种充填节理的侧开三角形单边裂纹(singlecleavage triangle,简称SCT)岩体试件进行动态冲击压缩试验,对比分析裂纹动态扩展过程及冲击破坏模式,并通过试验-数值法研究了Ⅰ型裂纹动态扩展过程中的应力强度因子及能量释放率的变化规律。研究表明:充填节理岩体主要有3种破坏模式,分别为预制裂纹扩展贯穿整个试件、裂纹贯穿后发生充填节理破坏、充填节理先破坏并阻滞裂纹贯穿;节理充填体的破坏与充填体强度、胶结力及应变率等相关;裂纹扩展速度自起裂处振荡增长,并在贯穿充填节理前某个位置达到最大,而应力强度因子及能量释放率在节理附近达到最大;充填节理刚度退化会抑制裂纹的扩展,并造成能量释放率的急剧增大;岩性组合差异加剧岩体试件材料属性的失配,导致应力强度因子及能量释放率随裂纹穿过节理后出现不同程度降低。

Detection of internal crack growth in polyethylene pipe using guided wave ultrasonic testing

Despite the success of guided wave ultrasonic inspection for internal defect detection in steel pipes,its application on polyethylene(PE)pipe remains relatively unexplored.The growth of internal cracks in PE pipe severely affects its pressure-holding capacity,hence the early detection of internal cracks is crucial for effective pipeline maintenance strategies.This study extends the scope of guided wave-based ultrasonic testing to detect the growth of internal cracks in a natural gas distribution PE pipe.Laboratory experiments and a finite element model were planned to study the wave-crack interaction at different stages of axially oriented internal crack growth with a piezoceramic transducer-based setup arranged in a pitch-catch configuration.Mode dispersion analysis supplemented with preliminary experiments was performed to isolate the optimal inspection frequency,leading to the selection of the T(0,1)mode at 50-kHz for the investigation.A transmission index based on the energy of the T(0,1)mode was developed to trace the extent of simulated crack growth.The findings revealed an inverse linear correlation between the transmission index and the crack depth for crack growth beyond 20%crack depth.

(Ⅰ+Ⅱ)MIXED MODE HYDROGEN INDUCED CRACKING

Hydrogen induced cracking(HIC)of 0.3% C,1% Cr,1% Mn,1% Si high strength steel has been studied under simple mode Ⅰ,mode Ⅱ and(Ⅰ+Ⅱ)mixed mode loading conditions.Af- ter being hydrogen-charged in IN H_2SO_4 solution,the material behaved hydrogen embrittlement in all the cases studied.The threshold K_(ⅡH)/K_(ⅡX) of HIC under mode Ⅱ load- ing was 0.27,which was nearly the same as that K(ⅠH)/K_(ⅠX)=0.29 under mode Ⅰ loading. While the thresholds of-HIC under(Ⅰ+Ⅱ)mixed mode loading were 0.36,0.41 and 0.37 cor- responding to the K_Ⅱ/K_Ⅰ ratio of 0.27,0.4 and O.81.The results show that simple mode Ⅰ or mode Ⅱ loading is more susceptible to hydrogen embrittlement than(Ⅰ+Ⅱ)mixed mode. For explaining the experimental results,the effects of triaxial stress as well as plastic deformation ahead of crack tip has been discussed.

THE ANALYSES OF THE THREE-DIMENSIONAL STRESS STRUCTURE NEAR THE CRACK TIP OF MODE I CT SPECIMENS IN ELASTICPLASTIC STATE(Ⅰ)—THE ANALYSES OF CONSTRAINT PARAMETERS AND FRACTURE PARAMETERS

In the present paper, three dimensional analyses of some general constraint parameters and fracture parameters near the crack tip of Mode I CT specimens in two different thicknesses are carried out by employing ADINA program. The results reveal that the constraints along the thickness direction are obviously separated into two parts: the keeping similar high constraint field (Z1) and rapid reducing constraints one(Z2). The two fields are experimentally confiremed to correspond to the smooth region and the shear lip on the fracture face respectively. So the three dimensional stress structure of Mode I specimens can be derived through discussing the two fields respectively. The distribution of the Crack Tip Opening Displacement (CTOD) along the thickness direction and the three dimensional distribution of the void growth ratio (Vg) near the crack tip are also obtained. The two fracture parameters are in similar trends along the thickness direction, and both of them can reflect the effect of thickness and that of the loading level to a certain degree.

A PENNY-SHAPED CRACK IN A TRANSVERSELY ISOTROPIC PIEZOELECTRIC SOLID:MODES Ⅱ AND Ⅲ PROBLEMS

An exact analysis of the modes Ⅱ and Ⅲ problems of a penny- shaped crack in a transversely isotropic piezoelectric medium is performed in this paper.The potential theory method is employed based on the general solution of three-dimensional piezoelasticity and the four harmonics involved are represented by one complex potential.Previous results in potential theory are then utilized to obtain the exact solution that is expressed in terms of elementary functions.Comparison is made between the current results with those published and good agreement is obtained.

MECHANISM ANALYSIS OF THICKNESS EFFECT ON MIXED MODE Ⅰ/Ⅱ FRACTURE OF LC4-CS ALUMINUM ALLOY

Mixed mode Ⅰ/Ⅱ fracture erperiments of LC4-CS aluminum alloy were conductedby using tension--shear specimens with thicknesses of 2, 4, 8 and 14mm. Fracturemechanisms of thickness effect on mixed mode Ⅰ/Ⅱ fracture were first examined fromfracture surface morphology to correlate with the macroscopic fracture behavior andstress state. It is found that specimen thickness has a strong influence on mixed modefracture. As thickness varies from thin to thick the macroscopic fracture surfacesappear the characteristics of plane stress state (2mm, 4mm--thick specimen), three--dimensional stress state (8mm--thick specimens), and plane strain state (14mm--thickspecimens), respectively. The specimens of all kinds of thicknesses are typical of ten-sile type failure under mode Ⅰ loading condition and shear type failure under mode Ⅱloading condition. Two distinct features coexist on the fracture surfaces under mixedmode loading conditions, and the corresponding proportion varies with loading mix-ity. Void--growth processes are the failure mechanism in both predominately tensile-and shears--type fractures. The size and depth of dimples on the fracture surface varygreatly with thickness. Therefore, it is extraordinary necessary to take into accountthe thickness effect when a mixed mode fracture criterion is being established.

In-plane shear(ModeⅡ) crack sub-critical propagation of rock at high temperature

In-plane shear crack sub-critical propagation of rock at high temperature was studied by finite element method and shear-box(i.e.compression-shear) test with newly designed electrically conductive adhesive method.Numerical and experimental results show that the normalized shear(Mode Ⅱ) stress intensity factors,K ⅡT/KT0 is decreased as the temperature increases because high temperature can improve stress distribution at crack tip and reduce the Mode Ⅱ stress intensity factor.Microscopic features of fractured surface are of little pits and secondary micro-cracks in the vicinity(1.5-4.0 mm) of the crack tip.The chevron-shape secondary cracks gradually merge in the length of about 4-5 mm and disappear along the direction of crack propagation.Stable shear crack propagation time is increased with the increasing temperature while the stable shear crack propagation rate is decreased with the increasing temperature,since high temperature can increase the shear(Mode Ⅱ) fracture toughness and prevent the crack growth.It is necessary to ensure the ligament of specimen long enough to measure the maximum unstable crack propagation rate of rock.

Effect of holes on in-plane shear(Mode Ⅱ) crack sub-critical propagation of rock

Shear-box(i.e.compression-shear) test and newly designed electrically conductive adhesive method were used to measure shear crack sub-critical propagation time and rate of sandstone specimen.Different cubic specimens with and without holes were tested to study the effect of holes on the shear crack sub-critical propagation.Numerical and experimental results show that three independent variables of hole,the interval distance S,the distance between the center of hole and the crack tip L,and hole radius R,have different contribution to the ratio of stress intensity factor of the specimen with holes to that of the specimen without hole,KⅡ/KⅡ0.Increasing S and decreasing L and R will result in the decrease of KⅡ/KⅡ0 and help crack arrest.The weight relation of the independent variables for KⅡ/KⅡ0 is S>L>R.The specimen DH3 with the largest value of S and the smallest values of L and R has the longest sub-critical crack propagation time and the smallest sub-critical crack propagation rate.Adding two suitable holes symmetrically to the original crack plane in rock specimen is considered to be a potential method for crack arrest of rock.

EXACT SOLUTIONS OF NEAR CRACK LINE FIELDS FOR MODE I CRACK UNDER PLANE STRESS CONDITION IN AN ELASTIC－PERFECTLY PLASTIC SOLID

The near crack line analysis method has been used in the present paper,The classical small scale yielding conditions have been completely abandoned in the analyses and one inappropriate matching condition used to be used at the elasticplastic boundary has been corrected.The reasonable solution of the plastic stresses near the crack line region has been established.By matching the plastic stresses with the exact elastic stresses at the elastic-plastic boundary,the plastic stresses the length of the plastic zone and the unit normal vector of the elastic-plastic boundary near the crock line region have been obtained for a mode I crack under uniaxial tension,as well as a mode I crack under biaxial tension,which shows that for both conditions the plastic stress componentsσy, and σsy.he length of the plastic zone and the unit normal vector of the elastic-plastic boundary are quite the same while the plastic stress σs is different.

Interface crack problems for mode Ⅱ of double dissimilar orthotropic composite materials

The fracture problems near the similar orthotropic composite materials are interface crack tip for mode Ⅱ of double disstudied. The mechanical models of interface crack for mode Ⅱ are given. By translating the governing equations into the generalized hi-harmonic equations, the stress functions containing two stress singularity exponents are derived with the help of a complex function method. Based on the boundary conditions, a system of non-homogeneous linear equations is found. Two real stress singularity exponents are determined be solving this system under appropriate conditions about bimaterial engineering parameters. According to the uniqueness theorem of limit, both the formulae of stress intensity factors and theoretical solutions of stress field near the interface crack tip are derived. When the two orthotropic materials are the same, the stress singularity exponents, stress intensity factors and stresses for mode II crack of the orthotropic single material are obtained.

A HALF PLANE CRACK UNDER THREE-DIMENSIONAL COMBINED MODE IMPACT LOADING

The dynamic stress intensity factor history for a half plane crack in an otherwise unbounded elastic body, with the crack faces subjected to a traction distribution consisting of two pairs of suddenly-applied shear line loads is consid- ered. The analytic expression for the combined mode stress intensity factors as a function of time is obtained. The method of solution is based on the application of integral transforms and the Wiener-Hopf technique. Some features of the solutions are discussed and graphical numerical results are presented.

A THEORETICAL MODEL OF MATRIX CRACKING IN COMPOSITE LAMINATES UNDER THERMOMECHANICAL LOADING

A theoretical approach is presented for analyzing the ply crackingin general symmetric lami- nates subjected to any combination ofin-plane mechanical loading and uniform temperature changes. Theequivalent constraint model proposed by the authors in a previouswork is used to account for the cracking in- teraction betweenlaminae in the laminates. By using a superposition schemce and thestress field solutions the energy release rate for a ply cracking isexplicitly as a function of stiffness reduction parameters of thelaminates. The ratio of mode Ⅰ to mode Ⅱ is introduced formconstruction of the fracture criterion. The effects of the laminateparameters and the crack spacing on the energy release rate and themode mixity are illustrated. Finally, the model is used to predictthe thermomechanical load for the first-ply-cracking.

A DISCUSSION ABOUT THE EFFECTIVENESS ON USING A VARIATE-TRANSFORMATION TO FIND OUT SOLUTIONS OF MODE Ⅲ CRACK PROBLEMS IN POWER HARDENING MEDIA

By the use of the transformations of physical plane to strain plane and physical planeto stress plane,an analytic expression of the asymptotic solution near a modeⅢcrack tipin a power hardening medium can be obtained.In this paper the effectiveness of thetransformation is discussed.Analytical results show that the transformation is effectiveexcept for a special limit case of power hardening media—the ideal plastic materials.

ROUGHNESS-INDUCED SHEAR RESISTANCE OF MODE Ⅱ CRACK GROWTH

It is obtained in this paper that the fatigue threshold value of mode H was 1.9 times of that of mode Ⅰ in dual-phase steel(DPS),and the normal stress intensity factor range oJ mode Ⅱ branch crack tip was 2.2 times of that of mode Ⅰ.Above results illustrate that the resistance of mode Ⅱ crack growth was higher than that of mode Ⅰ,the former resulting from roughness-induced shear resistance,the latter,crack closure. The mode Ⅱ component can play two important roles in near-threshold fatigue crack growth:(1)increasing crack tip plasticity which accelerates the crack growth and(2)intro- ducing crack surface contact and rubbing to reduce the crack propagation rate.By means of crack closure,the quantity of shear resistance was easily solved in this paper.The friction shear stress strength factor range of mode Ⅱ,K_,is still much higher than the closure stress strength factor range of mode Ⅰ,K_(Ⅰ,cl).This illustrated that the roughness enlarged the second role and played a role of shielding crack tip from mode Ⅱ crack.