界面超高速自相似动态分层分析：反平面情况

对材料界面超高速自相似动态分层的反平面问题进行了解析分析。分层模拟为界面裂纹由零长度自相似扩展，扩展速度为蹭音速或超音速。首先考虑运动集中载荷作用下界面动态分层的情况，利用界面裂纹自相似扩展的运动位错模型将问题归结为奇异积分方程，并求得解析解，分析了裂纹尖端的应力奇性，获得了动应力强度因子。最后，利用叠加原理给出了ｘ＾ｎ型载荷作用下界面动态分层的解。     

大直径SHPB实验中的高温加载技术及其应用

为研究材料的高温动态力学行为，提出一套由自主设计的温控系统和100 mm SHPB装置组成的高温SHPB实验系统，采用ANSYS软件对界面热传导及其对实验结果的影响进行了计算分析，论证了该实验技术的可靠性，并对混凝土的高温动态力学性能进行了研究。结果表明:在大直径合金钢材质SHPB装置上对混凝土等热惰性材料进行高温冲击实验，冷接触时间临界值为1.00 s，本文中提出的高温加载技术可将冷接触时间控制在0.50 s以内，实验技术可靠；同一加载速率下，随着温度从常温升到1 000 ℃，高温混凝土的动态应力应变曲线呈现出塑性变化趋势，动态抗压强度先提高后降低，动态峰值应变则不断增大。     

界面裂纹尖端的动态奇异特性

本文研究了界面裂纹尖端的动态应力场的奇异特性.引入尖端无摩擦接触的界面裂纹模型并采用具有运动边界的控制积分方程.证明了在动态界面裂纹尖端仅存在平方根奇异的应力场.数值结果表明接触区中的正应力确保持为压应力.为表现界面裂纹的动态特性,给出了应力强度因子和裂纹面接触区尺寸的数值结果.     

复杂界面(界面层)条件下的弹性波传播问题研究综述

界面（界面层）广泛存在于各种材料和结构中,并具有千差万别的形态．各种界面条件下的波动问题有着重要的理论与实际意义．综述了完好粘接界面、弱连接界面（或界面层）以及接触界面等复杂界面模型下弹性波传播问题的研究现状,主要集中于界面模型的建立、波传播问题的研究方法及主要结论．并提出了值得进一步研究的问题．     

温度循环载荷对柔性接头界面损伤的影响

固体火箭发动机在生产、运输和储存的过程中会受到环境湿热和老化载荷的作用,导致柔性接头界面的力学性能产生退化.为研究柔性接头界面力学性能退化对界面损伤的影响规律,基于双线性内聚力模型建立了一种描述界面力学性能退化的数学模型.以某柔性接头为研究对象,为便于研究,以温度循环载荷代替环境老化载荷开展了界面损伤分析,并采用ABAQUS 6.14进行了仿真.结果 表明,与后法兰粘接的界面损伤程度最大.当温度循环周期达到672次时,界面13和界面14的损伤程度分别由初始状态的2.9％和4.3％增加到5.2％和8.2％,增幅分别高达81.2％和91.7％.同时,温度循环载荷会加速柔性接头的界面损伤.该界面力学性能退化模型可为柔性接头老化研究提供参考.     

应用半权函数法计算界面裂纹的应力强度因子

应用半权函数法求解双材料界面裂纹的应力强度因子，得到以半权函数对参考位移与应力加权积分的形式表示的应力强度因子。针对特征值为复数λ的双材料界面裂纹裂尖应力和位移场，设置与之对应特征值为-λ的位移函数，即半权函数。半权函数的应力函数满足平衡方程，应力应变关系，界面的连续条件以及在裂纹面上面力为0；半权函数与裂纹体的几何尺寸无关，对边界条件没有要求。由功的互等定理得到应力强度因子KⅠ和KⅡ的积分形式表达式。本文计算了多种情况下界面裂纹应力强度因子的算例，与文献结果符合得很好。由于裂尖应力的振荡奇异性已经在积分中避免，只需考虑绕裂尖远场的任意路径上位移和应力，即使采用该路径上较粗糙的参考解也可以得到较精确的结果。     

固体炸药爆轰与惰性介质相互作用的一种扩散界面模型

提出一种保持热力学一致性的扩散界面模型，用来数值模拟固体炸药爆轰与惰性介质的相互作用问题。基于混合网格内各组分物质间可以达到力学平衡状态而不能达到热学平衡状态的假设，由混合网格能量守恒以及压力相等条件，推导出每种组分物质的体积分数演化方程。由此获得的扩散界面模型包括组分物质的质量守恒方程、混合物质的动量及总能量守恒方程，同时包括组分物质的体积分数演化方程和混合物质的压力演化方程。该扩散界面模型的主要特点是考虑了化学反应以及热学非平衡的影响。提出的扩散界面模型在物质界面附近不会出现物理量的非物理振荡现象、适用于任意表达形式的物质状态方程以及任意数目的惰性介质。     

界面裂纹萌生与扩展的分子动力学模拟

运用分子动力学模拟方法研究了裂纹在界面端处萌生与沿界面扩展的临界条件. 模拟考虑了一双相材料的3种模型，即构成90°/90°和 90°/180°夹角的两个界面端和一个界面裂纹. 模拟采用了包含原子区域与连续区域的并发型多尺度模型，即在界面端尖端和裂纹尖端附近 采用分子动力学(MD)方法，MD区域之外则按照线弹性有限元方法分析. 结果表明，在断裂启动时刻，3个模型沿界面的最大应力均达到界面理想强度；而且，其界 面能恰好足以克服界面材料的本征内聚能. 因此，界面端裂纹萌生与沿界面扩展的断裂条件可以通过界面理想强度和内聚能联系起来. 并基于模拟计算结果提出了界面断裂启动的统一准则.     

钛/钢复合板爆炸焊接实验

以3mm 厚的TA2钛板和26mm 厚的正火态Q345R为材料,通过爆炸焊接实验,对钛/钢复合板 爆炸焊接的动态参数进行了研究。结合复合板结合界面特征、复合板结合强度(剪切强度)以及界面波的金相 组织,讨论了钛/钢爆炸焊接时获得高强度结合和规则的界面结合波状形态的条件。对于3mm 厚TA2与 26mm厚正火态Q345R,该条件是动态碰撞角17,动态碰撞速度vp760m/s。根据界面波及基板轧制 金相组织形态,分析了形成界面波的机理,认为射流阻碍复板连续碰撞基板是形成界面波的一个主要原因。     

各向异性双材料界面裂纹扩展裂尖场

应用界面断裂力学理论和Stroh方法,研究了广义平面变形下动态裂纹沿着各向异性双材料界面扩展时的裂尖奇异应力及动态应力强度因子.双材料界面的动态裂尖区域特性主要由两个实矩阵W和D确定,且裂尖奇异应力和动态应力强度因子可以由包含这两个矩阵的柯西奇异积分方程确定,同时给出了动态应力强度因子和能量释放率的显示表达式.算例得出当裂纹以小速度扩展时,裂尖振荡因子ε与静态时几乎相同,当界面裂纹扩展速度接近瑞利波速时,ε趋于无穷大;同时得出应力强度因子及能量释放率随裂纹扩展速度的变化关系.     

基于分子动力学的结合材料界面破坏准则

结合材料的破坏通常都是从界面或其附近发生的,但界面破坏的机理及其评价准则尚未十分清楚.采用分子动力学模拟方法,可以对结合材料的界面破坏过程进行模拟,从而获得结合材料的界面应力和界面破坏之间的关系.界面破坏可以分为奇异应力场作用下的破坏,和界面应力集中引起的破坏两种.虽然在分子动力学模拟中采用了高度简化的界面模型,但对界面破坏过程的模拟,仍可以帮助人们获得结合材料界面破坏过程的规律性认识.分别模拟远场作用下界面上存在初始裂纹和界面附近存在初始裂纹两种情况下的界面破坏,根据分子动力学模拟结果,提出了一个结合材料界面破坏的准则.     

双材料界面断裂力学模型与实验方法

纤维增强聚合物(FRP)质轻、高强, 可提高结构的刚度、强度、抗震性能和耐久性, 近年来在结构加固及工程改造 中得到广泛应用. FRP与传统复合材料之间形成双材料黏结界面, 界面断裂特性是决定双材 料结构性能的关键因素. 对双材料界面裂纹尖端应力场理论、界面裂纹模型、黏结界面 I型、II型及混合型断裂试验及理论研究现状进行综合评述和分析. 界面模型主要有经典梁/ 板理论和刚性节点模型、考虑剪切变形的双亚层理论和半刚性节点模型、基于双亚层理论的 柔性节点模型、考虑剪切变形的多层亚层理论和多亚层柔性节点模型、弹性地基梁模型以及 黏聚模型. 还介绍了双材料界面断裂力学在FRP-混凝土研究中的应用.     

An accurate interface reconstruction method using piecewise circular arcs

A novel piecewise circular interface construction (PCIC) method for accurate reconstruction of interface in a two‐phase flow problem is proposed. This is under the framework of a fixed grid, volume of fluid approach applied on a two‐dimensional semistaggered structured grid. Fluid interface in each mixed cell is represented using a geometric template of piecewise circular arc. Data corresponding to arc center coordinates and radius are first predicted using curve fitting methods and corrected with the help of volume fraction constraints. Further corrections are carried out to achieve function (c0) continuity at cell boundaries. The proposed method does not require additional calculations for the determination of curvature (for calculation of surface tension force), since it is obtained as part of reconstruction process itself. For dynamic interface construction, simple analytical expressions are derived to construct edge matched flux polygons. Area of intersection of flux polygons with area covered by primary fluid is determined to effect geometric advection across a PCIC interface. Accuracy of this method is demonstrated by the reconstruction of standard static and dynamically evolving interface problems. Accuracy levels superior to most interface reconstruction methods using PLIC and schemes using higher order curves are established. Finally, the capability to handle a complex two‐phase flow problem simulation viz the four‐vortex flow field, where interface undergoes breakage and coalescence, is also demonstrated.     

Anti-plane waves near an interface between two piezoelectric half-spaces

We show the existence of certain new waves that can propagate near an interface between two half-spaces of different piezoelectric ceramics, where the interface is modeled by a membrane with the surface/interface elasticity [1]. The current configuration can be reduced into a number of well-known results as special cases, such as Love wave, Bleustein and Gulyaev wave. Together with our previous work for the imperfect interface [2], a full range of consideration of the interface affecting the anti-plane waves is now completed.     

Fracture behavior of a bonded magneto-electro-elastic rectangular plate with an interface crack

In this paper the anti-plane problem for an interface crack between two dissimilar magneto-electro-elastic plates subjected to anti-plane mechanical and in-plane magneto-electrical loads is investigated. The interface crack is assumed to be either magneto-electrically impermeable or permeable, and the position of the interface crack is arbitrary. The finite Fourier transform method is employed to reduce the mixed boundary-value problem to triple trigonometric series equations. The dislocation density functions and proper replacement of the variables are introduced to reduce these series equations to a standard Cauchy singular integral equation of the first kind. The resulting integral equation together with the corresponding single-valued condition is approximated as a system of linear algebra equations which can be easily solved. Field intensity factors and energy release rates are determined numerically and discussed in detail. Numerical results show the effects of crack configuration and loading combination parameters on the fracture behaviors of crack tips according to energy release rate criterion. The study of this problem is expected to have applications to the investigation of dynamic fracture properties of magneto-electro-elastic materials with cracks.     

Cohesive fracture simulation and failure modes of FRP–concrete bonded interfaces

A work-of-fracture method using three-point bend beam (3PBB) specimen, commonly employed to determine the fracture energy of concrete, is adapted to evaluate the mode-I cohesive fracture of fiber reinforced plastic (FRP) composite–concrete adhesively bonded interfaces. In this study, a bilinear damage cohesive zone model (CZM) is used to simulate cohesive fracture of FRP–concrete bonded interfaces. The interface cohesive process damage model is proposed to simulate the adhesive–concrete interface debonding; while a tensile plastic damage model is used to account for the cohesive cracking of concrete near the bond line. The influences of the important interface parameters, such as the interface cohesive strength, concrete tensile strength, critical interface energy, and concrete fracture energy, on the interface failure modes and load-carrying capacity are discussed in detail through a numerical finite element parametric study. The results of numerical simulations indicate that there is a transition of the failure modes controlling the interface fracture process. Three failure modes in the mode-I fracture of FRP–concrete interface bond are identified: (1) complete adhesive–concrete interface debonding (a weak bond), (2) complete concrete cohesive cracking near the bond line (a strong bond), and (3) a combined failure of interface debonding and concrete cohesive cracking. With the change of interface parameters, the transition of failure modes from interface debonding to concrete cohesive cracking is captured, and such a transition cannot be revealed by using a conventional fracture mechanics-based approach, in which only an energy criterion for fracture is employed. The proposed cohesive damage models for the interface and concrete combined with the numerical finite element simulation can be used to analyze the interface fracture process, predict the load-carrying capacity and ductility, and optimize the interface design, and they can further shed new light on the interface failure modes and transition mechanism which emulate the practical application.     

Crack growth of an interface crack between two dissimilar magneto-electro-elastic materials under anti-plane mechanical and in-plane electric magnetic impact

This paper examines the dynamic response of an interface crack between two dissimilar magneto-electro-elastic materials subjected to the mechanical and electric magnetic impacts. The magneto-electric impermeable boundary conditions are adopted. Laplace and Fourier transforms and dislocation density functions are employed to reduce the mixed boundary value problem to Cauchy singular integral equations in Laplace transform domain, which are solved numerically. Lots of numerical results are given graphically in time domain. The effects of electric impact loading and magnetic impact loading on dynamic energy density factors are discussed. Crack growth and propagation is predicted. The study of this problem is expected to have applications to the investigation of dynamic fracture properties of magneto-electro-elastic materials with cracks.     

Application of bimaterial interface corner failure mechanics to silicon/glass anodic bonds

Motivated by the existence of a universal singular stress field at bimaterial interface corners, a fair amount of work has been performed to support the use of the corresponding critical stress intensities to correlate fracture initiation. The approach is in the spirit of interface fracture mechanics but applicable to a different class of problems, specifically, when a crack does not previously exist (or cannot be detected, at least economically), and when subsequent crack propagation does not necessarily occur along the interface. Here we further progress toward the development, understanding, and application of the approach, both experimentally and theoretically, for a series of silicon/glass anodically bonded structures. To this end we designed and fabricated two series of silicon/glass anodically bonded bimaterial specimens with different interface corner geometries that commonly arise from different silicon etching technologies. Offset three-point flexure tests were performed that resulted in brittle fracture that initiated at the interface corner. From a rigorous stress analysis at the interface corner, we determined the order of the stress singularities and the angular variation of the stress fields. We computed the corresponding stress intensities via full-field finite element analyses of the silicon/glass specimens loaded in offset three-point flexure. Measured fracture data show that although the failure stress varies significantly with bond size, the corresponding critical stress intensity of the dominant mode is constant, thus providing support for its use as a fracture initiation criterion. In the light of both the stress analysis and the measured fracture data, we discuss the effect of mode mixity (loosely shearing versus opening) and show that it has little influence on the results for the specimens and loading considered in this study. Via an idealized model of a small crack, either interfacial or extending into one of the adherends, we study the effects of geometrical perturbations at the interface corner on the stress state, and discuss implications for fracture analysis and interpretation of fracture data. We also explore the prediction of the crack initiation angle and achieve reasonable success with a simple criterion based on the maximum circumferential stress near the uncracked interface corner.