电子封装元件的高阶层离散均匀化分析

将均匀化理论应用于具有非完全(单层内)周期性微结构的倒装焊底充胶电子封装元件，建立了高阶逐层离散层板模型，用解析法分析热载荷下结构的温度应力. 计算结果与有限元解的比较表明，该分析模型和方法是有效的，而且比较简便. 算例分析结果显示，胶层厚度、焊点密度、胶与焊点材料的模量比和体积比，对于焊点温度应力有明显影响.     

高密度封装倒装焊的新简化模型及疲劳寿命分析

高密度封装结构中存在大量焊球，在进行有限元建模时，考虑到焊球数量多且尺寸小的特点，需简化焊球层模型．本文针对高密度封装结构中倒装焊及底填胶，提出了一种新的简化模型，采用非重点部位简化和材料均匀化相结合的方法，将内部焊球层和底填胶简化为均匀层，只保留外圈几层焊球．通过建立代表性体积单元，计算得到均匀化层的材料参数．讨论了外层焊球圈数对焊球层危险点应力的影响，发现保留外圈两层焊球就能得到非常精确的结果．利用新简化模型计算了高密度封装结构的疲劳寿命，所获得的结果与未均匀化模型结果误差在0.3%以内．     

基于有限断裂法和比例边界有限元法的裂纹扩展模拟

基于有限断裂法和比例边界有限元法提出了一种裂缝开裂过程模拟的数值模型。采用基于有限断裂法的混合断裂准则作为起裂及扩展的判断标准,当最大环向应力和能量释放率同时达到其临界值时,裂缝扩展。结合多边形比例边界有限元法,可以半解析地求解裂尖区域附近的应力场和位移场,在裂尖附近无需富集即可获得高精度的解。计算能量释放率时,只需将裂尖多边形内的裂尖位置局部调整,无需改变整体网格的分布,网格重剖分的工作量降至最少。裂缝扩展步长通过混合断裂准则确定,避免了人为假设的随意性,并可以实现裂缝变步长扩展的模拟,更符合实际情况。通过对四点剪切梁的复合型裂缝扩展过程的模拟,对本文模型进行了验证,并应用于重力坝模型的裂缝扩展模拟,计算结果表明,本文提出的模型简单易行且精度较高。     

线弹性广义界面层模型的研究

为了用界面层方法分析一般结构中的裂纹扩展问题,提出了广义界面层的概念.基于能量释放率G相等的原则,对界面层的本构关系进行了线弹性等效.在此基础上,重点推导了含Ⅰ型裂纹双悬梁(DCB)的载荷-位移函数,并借助Abaqus软件对双悬梁的分层问题进行数值模拟.理论推导和数值模拟的结果证明了广义界面层的线弹性等效是可行的.     

焊点/底充胶夹层等效弹性常数的简化分析

将高阶逐层离散层板模型和均匀化理论相结合,应用于具有非完全(单层内)周期性微结构的多层结构,用解析法分析了电子封装结构中由焊点和底充胶构成的非均质夹层的等效弹性常数．计算结果与已有结果进行了比较,验证了文中分析方法的有效性和简便性．     

橡胶—钢粘接界面断裂韧性实验研究

采用双悬臂夹层梁试样,对橡胶—钢粘接界面的I型边缘裂纹的扩展情况及断裂韧性进行实验研究.通过声发射技术的监测将裂纹扩展划分为无损伤、开裂、失稳扩展三个阶段,把线弹性断裂力学中的柔度和能量释放率引入到橡胶—钢粘接界面断裂韧性的计算中,由实验测出裂纹的张开位移计算出夹层梁的柔度曲线,进而得到了能量释放率曲线,并对断裂韧性的变化规律进行了分析.     

双拉试样裂纹扩展过程中的应力强度因子计算

论文将使用一种界面单元来解决二维裂纹的静态扩展问题.这种界面单元基于虚拟裂纹闭合法,利用商业有限元软件ABAQUS的用户自定义单元UEL功能,发展为界面断裂单元,计算应变能释放率(GⅠ和GⅡ).在裂纹尖端的两个节点间设置一个特殊刚度的弹簧,并引入哑节点计算裂纹尖端后面的张开位移和裂纹尖端前面的虚拟裂纹扩展量.采用这种单元计算应变能释放率时不需要使用奇异单元或折叠单元,不会出现收敛问题,也不需要复杂的后续处理.因此,采用这种断裂单元分析二维裂纹扩展问题是方便的、高效率的,而且也能得到可靠的精度.     

基于非线性弹簧界面的复合材料分层损伤预测

基于功能等效的原理,用自定义本构的非线性弹簧模拟复合材料层间界面,弹簧的刚度随着界面层的牵引力-位移曲线变化,可以表征界面材料性能的线性和非线性退化过程。对复合材料单一型分层和混合型分层损伤演化过程进行了模拟分析,研究了单元尺寸、界面强度等参数对模拟结果的影响机理。结果表明:非线性弹簧界面单元能够准确地模拟分层损伤的起始和扩展过程,使用非线性弹簧模拟界面可以减小计算模型规模,与内聚力单元相比计算效率提高约10倍;裂纹每扩展一个单元长度所需要的能量与单元尺寸成正比,单元尺寸越大,裂纹扩展所需要的能量越大;界面强度越低,初始裂纹尖端张开过程越平缓,模型的收敛性越好,可以通过降低弹簧单元界面强度来减小模型计算规模。     

金属材料三维断裂韧度厚度效应的有限元模拟

随着金属材料大壁厚结构件在工程中的广泛应用,对其断裂韧度的厚度效应研究具有重要的科学意义和工程价值。本研究基于有限元和实验相结合的方法,对金属材料断裂韧度的厚度效应进行预测。首先,通过一组薄壁厚金属材料标准三点弯曲试验得到试样失效时的临界载荷值,并利用内聚力模型与基于虚拟裂纹闭合技术的裂纹扩展模拟方法得到裂纹扩展时的单元临界能量释放率。随后,以此临界能量释放率作为裂纹扩展的启裂准则门槛值,通过有限元计算得到不同试样厚度下裂纹启裂时的裂尖断裂参数随着厚度的变化规律。最后,为了验证有限元模拟结果的准确性,本研究进行了另外两组不同厚度下三点弯曲试样的断裂韧度试验,并将试验结果与有限元结果进行了对比,验证了有限元所模拟的断裂韧度厚度效应的准确性。本研究旨在,通过薄壁厚三点弯曲试样的实验结果结合有限元模拟工作,即可实现金属材料断裂韧度的整个厚度效应曲线,为任意厚度下金属材料断裂韧度预测提供一种可靠的研究方法,有益于缩减试验成本,为大壁厚工程结构件的失效预测提供依据。     

石墨烯裂纹扩展行为研究

基于分子动力学方法对含预制裂纹石墨烯进行扶手椅向拉伸断裂模拟。使用连续介质理论结合分子动力学计算石墨烯能量释放率,确定石墨烯能量释放率GIC为10.25 J/m2;应力强度因子KIC为 3.33MPam^1/2。进一步对影响石墨烯裂纹扩展速率的因素-初始裂纹长度与加载速率进行讨论。结果表明：裂纹初始长度与加载率会在一定程度上影响石墨烯中裂纹扩展速率。裂纹扩展速率会随着初始裂纹长度的增加而降低;但随着初始裂纹长度的增加,裂纹扩展速率对其敏感度降低。裂纹扩展速率会随着加载率的升高而增大。 初始裂纹长度与加载率对裂纹扩展速率的影响有一定的关联性,加载率的升高会降低裂纹扩展速率对初始裂纹长度变化的敏感度。在此基础上确定了石墨烯中裂纹扩展极限速率为8350 m/s。关联性,加载率的升高会降低裂纹扩展速率对初始裂纹长度变化的敏感度。在此基础上确定了石墨烯中裂纹扩展极限速率为8350 m/s。     

Quasicontinuum simulation of single crystal nano-plate with a mixed-mode crack

The quasicontinuum (QC) method is employed to simulate a nickel single crystal nano-plate with a mixed-mode crack. Atomic stresses near the crack tip are fitted according to the elastoplastic fracture mechanics equations. It is found that the atomic stress fields neighboring the crack tip are also singular and controlled by the atomic stress intensity factors. And then the critical energy release rates for brittle and ductile fracture are computed and compared in order to predict crack propagation or dislocation emission. Four possible slip directions at the crack tip are pointed out. Finally, the slip direction around the crack tip is determined by the shear stress and it is well consistent with the atomic pictures from the QC simulation.     

冲击载荷下扩展裂纹尖端动态能量释放率分布的焦散线分析

借助高速摄影捕捉裂纹瞬态扩展过程，利用动态焦散线研究了含有裂纹的三点弯曲梁在冲击载荷作用下扩展裂纹尖端的动态能量释放率分布规律；综合分析了裂纹扩展时间、长度、速度，以及扩展裂纹尖端动态应力强度因子与它的变化关系，表明了动态能量释放率在裂纹扩展过程中的驱动作用。     

Dynamic crack tip fields and dynamic crack propagation characteristics of anisotropic material

Based on mechanics of anisotropic material, the dynamic crack propagation problem of I/II mixed mode crack in an infinite anisotropic body is investigated. Expressions of dynamic stress intensity factors for modes I and II crack are obtained. Components of dynamic stress and dynamic displacements around the crack tip are derived. The strain energy density theory is used to predict the dynamic crack extension angle. The critical strain energy density is determined by the strength parameters of anisotropic materials. The obtained dynamic crack tip fields are unified and applicable to the analysis of the crack tip fields of anisotropic material, orthotropic material and isotropic material under dynamic or static load. The obtained results show Crack propagation characteristics are represented by the mechanical properties of anisotropic material, i.e., crack propagation velocity M and fiber direction α. In particular, the fiber direction α and the crack propagation velocity M give greater influence on the variations of the stress fields and displacement fields. Fracture angle is found to depend not only on the crack propagation but also on the anisotropic character of the material.     

On the fracture toughness of ferroelastic materials

The toughness enhancement due to domain switching near a steadily growing crack in a ferroelastic material is analyzed. The constitutive response of the material is taken to be characteristic of a polycrystalline sample assembled from randomly oriented tetragonal single crystal grains. The constitutive law accounts for the strain saturation, asymmetry in tension versus compression, Bauschinger effects, reverse switching, and strain reorientation that can occur in these materials due to the non-proportional loading that arises near a propagating crack. Crack growth is assumed to proceed at a critical level of the crack tip energy release rate. Detailed finite element calculations are carried out to determine the stress and strain fields near the growing tip, and the ratio of the far field applied energy release rate to the crack tip energy release rate. The results of the finite element calculations are then compared to analytical models that assume the linear isotropic K-field solution holds for either the near tip stress or strain field. Ultimately, the model is able to account for the experimentally observed toughness enhancement in ferroelastic ceramics.     

Effect of orthotropy on crack propagation

The tendency of moving cracks to spread along the preferred directions of material anisotropy is treated. Depending on the velocity of crack propagation, the change of material properties in orthogonal planes is shown to affect the bifurcation characteristics. The problem is reduced to a system of dual integral equations that can be solved in a standard fashion. Of particular interest is the dynamic stress field near the tip of a moving crack in an orthotropic material. Although the 1√r stress singularity is preserved with r being the radial distance measured from the crack tip, the angular variations of the stresses are dependent on crack speed and material anisotropy. The possibility of crack bifurcation is examined by application of the strain energy density criterion for several composite systems. Crack branching is found to be enhanced by material anisotropy, a phenomenon that is not uncommon in composite materials.     

Investigation on thermomechanical fracture in the framework of configurational forces

A configurational force approach is developed for providing a fresh look onto classical aspects of thermomechanical fracture. The theoretical framework is based on the finite deformation and makes no restrictions on the material response. The integral form of configurational force balance at the crack tip is constructed, and the concentrated configurational body force is decomposed into the inertial and internal parts. The energy release rate is evaluated through the generalized second law of thermodynamics applicable to configurational force system. The theoretical investigation shows that the negative of the projection of the internal configurational force concentrated at the crack tip along the direction of crack propagation plays the role of the energy release rate and acts directly in response to crack propagation. This finding enables us to deal with the thermomechanical fracture problems in material space.     

矩形压电体中反平面裂纹的电弹性场

基于线性压电理论,本文获得了含有中心反平面裂纹的矩形压电体中的奇异应力和电场。利用Fourier积分变换和Fourier正弦级数将电绝缘型裂纹问题化为对偶积分方程,并进一步归结为易于求解的第二类Fred-holm积分方程。获得了裂纹尖端应力、应变、电位移和电场的解析解,求得了裂纹尖端场的强度因子及能量释放率。分析了压电矩形体的几何尺寸对它们的影响。结果表明,对于电绝缘型裂纹,裂纹尖端附近的各个场变量都具有-1/2阶的奇异性,能量释放率与电荷载的方向及大小有关,并且有可能为负值。     

基于DIC的爆炸加载下脆性材料裂纹扩展规律的试验研究

通过数字图像相关法（DIC），应用PMMA对爆炸加载条件下脆性材料的裂纹扩展规律进行了试验研究。基于对称性试验模型，实现了裂纹尖端位置和应变场信息的同步记录。以此为基础，通过对比分析获知，主应变场应变值最大点不能作为裂纹尖端的判断依据。并以动态裂纹扩展速度为参量，应用断裂动力学和最小二乘牛顿迭代法，计算出了考虑惯性效应的Ⅰ-Ⅱ混合型裂纹的应力强度因子:K_Ⅰ和K_Ⅱ值会随着裂纹扩展方向改变而发生突变；K_Ⅰ最大值为2.63 MPa·m1/2，最小值为0.89 MPa·m1/2；其整体变化趋势表明，爆炸加载条件下脆性材料裂纹扩展随能量积聚和释放呈循环阶梯式递减发展。     

An energy-consistent fracture model for ferroelectrics

The fracture behavior of ferroelectrics has been intensively studied in recent decades, though currently a widely accepted fracture mechanism is still lacking. In this work, enlightened by previous experimental observations that crack propagation in ferroelectrics is always accompa-nied by domain switching, we propose a micromechanical model in which both crack propagation and domain switch-ing are controlled by energy-based criteria. Both electric energy and mechanical energy can induce domain switching, while only mechanical energy can drive crack propagation. Furthermore, constrained domain switching is considered in this model, leading to the gradient domain switching zone near the crack tip. Analysis results show that stress-induced ferroelastic switching always has a toughening effect as the mechanical energy release rate serves as the driving force for both fracture and domain switching. In compari-son, the electric-field-induced switching may have either a toughening or detoughening effect. The proposed model can qualitatively agree with the existing experimental results.     

无网格与有限元的耦合在动态断裂研究中的应用

提出将无网格Galerkin法与有限元耦合的方法用于分析动态裂纹扩展问题,只在裂尖附近区域沿裂纹扩展方向布置无网格结点,而在其他区域采用一般的有限元,区域交界处的结点采用MLS方法插值,然后将求得的结点值再分配到有限单元的相关结点上,保证了无网格区域和有限元区域的交界处位移的连续。避免了网格的再生成,同时也克服了单纯使用无网格Galerkin法所带来的边界条件难处理及计算效率较低的缺点。数值算例显示这种方法是有效的。