膜-基复合材料界面剪应力的三维半解析计算

 采用影响系数法对膜-基复合材料的界面剪应力三维半解析进行 了分析研究.利用三维有限元方法对薄膜的影响系数进行了计算. 将 基体作为半无限大体，利用其平面边界作用单位力时的位移场解析 解，得到基体的影响系数. 结果表明，对膜-基复合材料界面剪应力 进行三维半解析计算，克服了完全用三维有限元对其进行计算的限 制，为该类问题的分析提供了新的途径.     

纤维增强韧性基体界面力学行为

分析了纤维增强韧性基体的界面力学行为及其失效机理,按剪滞理论和应变理化规律研究微复合材料的弹塑性变形和应力状态,讨论了幂硬化和线性硬化基体的弹塑性变形和界面应力分布,并给出纤维应力和位移的表达式。按最大剪应力强度理论建立了纤维／基体界面失效准则,推导出弹塑性界面失效的平均剪应力随纤维埋入长度的变化关系。     

TiN薄膜的应力状态对摩擦学性能的影响

用Ｘ射线衍射仪测定了在５２１００钢基体上离子束增强沉积ＴｉＮ膜、等离子体化学气相沉积ＴｉＮ膜和离子镀ＴｉＮ膜的应力状态，分析了不同工艺方法制取的ＴｉＮ薄膜的应力形成的影响因素，比较了３种薄膜在不同载荷和摩擦速度条件下的摩擦学性能，分析了膜－基界面两侧应力状态对膜－基结合力、薄膜的耐磨性能和磨损机理的影响．结果表明：３种ＴｉＮ／５２１００钢试样在薄膜内的应力均为压应力，但在界面附近基体一侧的应力状态是随着工艺方法的不同而不同，３种膜的硬度和膜－基结合力都依次下降，而其内应力与膜－基应力的差值则是依次增大，分别为２６９．０ＭＰａ，６６０．５ＭＰａ和１０６３．３ＭＰａ，因而前者显示出最高的膜－基结合力和最佳的摩擦学性能；而后２种膜则显示出依次渐差的膜－基结合力和摩擦学性能     

薄膜涂层材料界面纯剪破坏标准试验法的开发

界面强度是薄膜涂层材料的最重要的性能指标之一，目前尚缺乏有效的测量方法。测量界面强度的主要困难在于寻求一种便于试验的试件形状和加载方式，使得界面上能够产生不同的应力状态，即在不同的剥离应力和剪应力比的状态下发生破坏。本文采用有限元数值分析（ABAQUS），研究了几种具有简单几何形状的试件的界面应力，并基于大量的数值试验，设计了剪应力为破坏支配因素的试件形状和加载方式，并且给出了便于应用的最大界面剪应力的经验估算公式。该经验公式可以适用于各种材料组合和薄膜涂层的厚度。研究结果表明，通过对薄膜涂层材料试件的基体引入缺口以产生应力集中，进行普通的四点弯曲试验，可以进行剪应力占支配地位的界面破坏试验。利用本文提出的试件形状和相应的最大界面剪应力经验公式，可以通过破坏试验简单地得到界面的剪切强度。     

微纳米尺度膜/基结构界面性能的实验研究

为了研究铜膜/有机玻璃结构的界面性能,首先对沉积在有机玻璃基底上300nm厚的铜膜进行了单轴压缩实验,部分区域的薄膜因屈曲而脱离基底。选择在膜/基粘接良好区域、膜/基脱粘区域分别进行等位移纳米压痕实验。利用膜/基粘接良好区域处硬度/弹性模量与压痕位移的关系来确定膜/基结构的临界脱粘位移。基于宏观力学中表征界面性能的能量法,利用两个区域等位移的塑性功差值来确定界面能量释放率。研究结果表明:当压痕位移约450nm时,膜/基结构开始出现界面脱粘,实验测得铜膜/有机玻璃结构的界面能量释放率值在6.81~10.32J/m2之间。     

模型复合材料弹塑性界面应力分析

由纤维增强弹塑性基体所产生的界面具有弹塑性力学行为。考虑到一般材料的塑性变形都遵循幂硬化规律,对模型复合材料的界面进行弹性和应变硬化状态下的变形规律及其应力分析。以纤维拔出试验为研究模型,将界面分成弹性区和塑性区。利用界面应力剪滞理论,分别建立弹性区和塑性区的界面力学基本方程。选择适当的位移函数满足基本方程及埋入纤维的边界条件,再按位移函数求出弹性区和塑性区的界面剪应力。推导出平均界面剪应力与纤维     

纤维增强陶瓷基复合材料疲劳迟滞回线模型研究

纤维增强陶瓷基复合材料初始加载到疲劳峰值应力时, 基体出现裂纹, 纤维/基体界面发生脱粘. 在疲劳载荷作用下, 纤维相对基体在界面脱粘区往复滑移使得陶瓷基复合材料出现疲劳迟滞现象. 建立了纤维陶瓷基复合材料疲劳迟滞回线细观力学模型, 采用断裂力学方法确定了初始加载纤维/基体界面脱粘长度、卸载界面反向滑移长度与重新加载新界面滑移长度, 分析了4种不同界面滑移情况的疲劳迟滞回线. 假设正交铺设与编织陶瓷基复合材料疲劳迟滞回线主要受0°铺层、轴向纱线内纤维/基体界面滑移的影响, 预测了单向、正交铺设与编织陶瓷基复合材料在不同峰值应力与不同循环的疲劳迟滞回线, 与试验结果吻合.      

INVESTIGATION ON FATIGUE HYSTERESIS LOOPS MODELS OF FIBRE-REINFORCED CERAMIC-MATRIX COMPOSITES

纤维增强陶瓷基复合材料初始加载到疲劳峰值应力时, 基体出现裂纹, 纤维/基体界面发生脱粘. 在疲劳载荷作用下, 纤维相对基体在界面脱粘区往复滑移使得陶瓷基复合材料出现疲劳迟滞现象. 建立了纤维陶瓷基复合材料疲劳迟滞回线细观力学模型, 采用断裂力学方法确定了初始加载纤维/基体界面脱粘长度、卸载界面反向滑移长度与重新加载新界面滑移长度, 分析了4种不同界面滑移情况的疲劳迟滞回线. 假设正交铺设与编织陶瓷基复合材料疲劳迟滞回线主要受0°铺层、轴向纱线内纤维/基体界面滑移的影响, 预测了单向、正交铺设与编织陶瓷基复合材料在不同峰值应力与不同循环的疲劳迟滞回线, 与试验结果吻合.     

含曲线型膜基界面的高分子基金属薄膜延展性能

柔性电子中联接电子元器件的互联金属导线多以附着在高分子基底上的薄膜形式存在。由于此类膜基体系在服役过程中需要承受相对较大的变形,如何改进高分子基金属薄膜的延展性能成为制约柔性电子技术发展的关键问题之一。以往的研究通过对高分子基底进行酸碱腐蚀、喷砂等表面糙化处理,虽然可以有效提高膜基结合性能,但却很少考虑基底表面糙化处理对提高膜基体系延展性能的影响。本文首先实验研究了在含糙化表面的聚酰亚胺基底上附着Cu膜的延展性能,结果表明,提高基底表面粗糙度能够显著降低Cu膜在拉伸条件下的裂纹密度。由于膜基体系表面裂纹的扩展与薄膜表面拉伸正应力分布相关,后者将直接影响薄膜的延展性,采用有限元方法模拟计算了基底表面糙化处理后,金属薄膜在拉伸状态下的应力分布。在计算模型中,膜基界面被处理成正弦曲线形式的理想化界面,并考虑了金属薄膜的外表面为平直状和曲线状两种情况。结果显示,曲线型界面可显著改变后一种情况下金属薄膜在拉伸状态下的表面正应力分布,从而达到抑制金属表面裂纹的扩展以及降低裂纹密度的作用。最后,采用内聚力模型模拟膜基界面,研究了在拉伸条件下曲线型界面的损伤分布情况。结果表明,相对于平直界面,曲线型界面不易发生如界面损伤和界面裂纹扩展的破坏,而且振幅波长比越大的曲线型界面越不容易发生破坏。     

纤维增强金属基复合材料整体蠕变性能的细观力学分析

高温下金属基复合材料的蠕变主要由基体蠕变和界面扩散蠕变两部分构成,以往的研究中常常只考虑其中一种蠕变机理,从而导致得到的规律具有较大的局限性．本文提出了一种可预测金属基复合材料整体蠕变性能的细观力学方法,同时考虑了基体蠕变和界面扩散蠕变两种蠕变机理,导出了具有张量形式并满足不可压缩性的界面扩散蠕变应变表达式．采用Mori-Tanaka法和自洽法二者结果的平均以便更准确地计算纤维中的应力,揭示了两种蠕变机理相互影响的竞争关系．研究了恒定双轴荷载下的总体蠕变和固定位移约束下的应力松弛这两种常见蠕变问题,探究了基体蠕变与界面扩散蠕变两种蠕变机理在总蠕变中发挥的作用,考察了不同加载条件和不同纤维体积分数对复合材料整体蠕变行为的影响．     

A coupled electromechanical analysis of a piezoelectric layer bonded to an elastic substrate: Part I,development of governing equations

This two-part contribution presents a novel and efficient method to analyze the two-dimensional (2-D) electromechanical fields of a piezoelectric layer bonded to an elastic substrate, which takes into account the fully coupled electromechanical behavior. In Part I, Hellinger–Reissner variational principle for elasticity is extended to electromechanical problems of the bimaterial, and is utilized to obtain the governing equations for the problems concerned. The 2-D electromechanical field quantities in the piezoelectric layer are expanded in the thickness-coordinate with seven one-dimensional (1-D) unknown functions. Such an expansion satisfies exactly the mechanical equilibrium equations, Gauss law, the constitutive equations, two of the three displacement–strain relations as well as one of the two electric field-electric potential relations. For the substrate the fundamental solutions of a half-plane subjected to a vertical or horizontal concentrated force on the surface are used. Two differential equations and two singular integro-differential equations of four unknown functions, the axial force, N, the moment, M, the average and the first moment of electric displacement, D₀ and D₁, as well as the associated boundary conditions have been derived rigorously from the stationary conditions of Hellinger–Reissner variational functional. In contrast to the thin film/substrate theory that ignores the interfacial normal stress the present one can predict both the interfacial shear and normal stresses, the latter one is believed to control the delamination initiation.     

Film thickness measurement with an ultrasonic transducer

A technique for measuring condensate film thickness using an ultrasonic transducer is described. In the experiment, the condensate film thickness with R-113 and FC-72 (a fluorinert compound developed by the 3M Company) condensing on the horizontal lower surface of a rectangular duct was measured at several locations. From the measured values a power law relation between the condensate film thickness and the axial distance from the leading edge of the condensing surface was derived by regression analysis. Assuming a linear temperature profile in the condensate film, local and average heat transfer coefficients were computed from the condensate film thickness. The average heat transfer coefficients were compared with the values obtained by measuring the heat transfer rate to the coolant. The two values were within ±12% of each other. As yet there is no satisfactory analytical model to predict the local heat transfer coefficient even in the annular condensation regime. One of the main difficulties in modeling the condensation is the lack of a suitable model to predict the interfacial shear stress. With the measurement of the film thickness it is possible to determine the interfacial shear stress. It is hoped that the shear stresses so determined will lead to the development of a satisfactory model for interfacial shear stress with condensation.     

含裂粘补结构胶层剪应力的计算

提出了一种计算含裂粘补结构胶层剪应力的新方法。采用有限元计算补片的影响系数，使用已有解析解计算蒙皮的影响系数。利用位移协调条件，列出有限多组联立方程。解出含裂粘补结构胶层剪应力。计算结果较好的反映了含裂粘补结构胶层剪应力的实际变化趋势，该方法同样可应用于其他结构如高梯度基膜复合材料剪应力的计算。     

Mixed-mode failure of thin films using laser-generated shear waves

A new test method is developed for studying mixed-mode interfacial failure of thin films using laser generated stress waves. Guided by recent parametric studies of laser-induced tensile spallation, we successfully extend this technique to achieve mixed-mode loading conditions. By allowing an initial longitudinal wave to mode convert at an oblique surface, a high amplitude shear wave is generated in a fused silica substrate and propagated toward the thin-film surface. A shear wave is obtained with amplitude large enough to fail an Al film/fused silica interface and the corresponding shear stress calculated from high-speed interferometric displacement measurements. Examination of the interfaces failed under mixed-mode conditions reveals significant wrinkling and tearing of the film, in great contrast to blister patterns observed in similar Al films failed under tensile loading.     

Effect of geometrical and material parameters in nanoindentation of layered materials with an interphase

Indentation models for thin layer-substrate geometry with an interphase have been developed. The interphase can be modeled either as a nonhomogeneous layer or as a homogeneous layer. Between the two models of the interphase, contact depth and critical interfacial stresses are compared to find the effect of indentation area, film and substrate Young’s moduli, and the interphase and film thicknesses. Although contact depth is found not to be sensitive to the type of interphase model used, critical interfacial stresses are significantly different (up to 15%) for film to substrate elastic Young’s moduli ratios of more than 25. A formal sensitivity analysis based on design of experiments shows that on critical interfacial stresses, interphase to film thickness ratio and film to substrate Young’s moduli ratio has the most impact, while type of elastic moduli variation in the interphase and indentor width to film thickness ratio has the least impact.     

Tensile and mixed-mode strength of a thin film-substrate interface under laser induced pulse loading

Laser induced stress waves are used to characterize intrinsic interfacial strength of thin films under both tensile and mixed-mode conditions. A short-duration compressive pulse induced by pulsed-laser ablation of a sacrificial layer on one side of a substrate is allowed to impinge upon a thin test film on the opposite surface. Laser-interferometric measurements of test film displacement enable calculation of the stresses generated at the interface. The tensile stress at the onset of failure is taken to be the intrinsic tensile strength of the interface. Fused-silica substrates, with their negative nonlinear elasticity, cause the compressive stress wave generated by the pulse laser to evolve a decompression shock, critical for generation of the fast fall times needed for significant loading of surface film interfaces. By allowing the stress pulse to mode convert as it reflects from an oblique surface, a high amplitude shear wave with rapid fall time is generated and used to realize mixed-mode loading of thin film interfaces. We report intrinsic strengths of an aluminum/fused silica interface under both tensile and mixed-mode conditions. The failure mechanism under mixed-mode loading differs significantly from that observed under pure tensile loading, resulting in a higher interfacial strength for the mixed-mode case. Inferred strengths are found to be independent, as they should be, of experimental parameters.     

Instabilities in free-surface electroosmotic flows

Instability of a thin electrolyte film undergoing a direct current electroosmotic flow has been investigated. The film with a compliant electrolyte–air interface is flowing over a rigid charged substrate. Unlike previous studies, inclusion of the Maxwell stresses in the formulation shows the presence of a new finite wavenumber shear-flow mode of instability, alongside the more frequently observed long-wave interfacial mode. The shear mode is found to be the dominant mode of instability when the electrolyte–solid and electrolyte–air interfaces are of opposite charge or of same charge but have very large zeta-potential at the electrolyte–air interface. The conditions for mode-switch (interfacial to shear) and the direction of the travelling waves are discussed through stability diagrams. Interestingly, the analysis shows that when the interfaces are of nearly same zeta potential, the ‘free’ electrolyte–air interface behaves more like a ‘stationary’ wall because of the ion transport in the reverse direction of the flow.     

A new zig-zag theory and C<Superscript>0</Superscript> plate bending element for composite and sandwich plates

A higher-order zig-zag theory for laminated composite and sandwich structures is proposed. The proposed theory satisfies the interlaminar continuity conditions and free surface conditions of transverse shear stresses. Moreover, the number of unknown variables involved in present model is independent of the number of layers. Compared to the zig-zag theory available in literature, the merit of present theory is that the first derivatives of transverse displacement have been taken out from the in-plane displacement fields, so that the C⁰ interpolation functions is only required during its finite element implementation. To obtain accurately transverse shear stresses by integrating three-dimensional equilibrium equations within one element, a six-node triangular element is employed to model the present zig-zag theory. Numerical results show that the present zig-zag theory can predict more accurate in-plane displacements and stresses in comparison with other zig-zag theories. Moreover, it is convenient to obtain transverse shear stresses by integration of equilibrium equations, as the C⁰ shape functions is only used when implemented in a finite element.