超弹性薄膜与可压缩基底双层结构表面失稳分析

当上层超弹性硬质薄膜和下层可膨胀基底构成的双层结构受压时,薄膜的自由表面可通过形成褶皱降低系统能量.研究表明,上下两层的模量比不同时,上层弹性硬质薄膜将表现出不同的表面失稳模式.本文提出了一种新颖的方法可有效抑制双层软材料的表面失稳,即改变基底材料的泊松比,这种方法同时适用于不具有应变硬化的软材料.首先基于Neo-Hookean模型发展了小变形条件下双层结构表面失稳的理论模型,通过半解析的方法得到了表面失稳的临界应变;然后通过有限元计算与模拟,进一步验证了负泊松比基底可延缓表面失稳.结果表明:(1)当双层结构基底泊松比为正且趋于0.5(不可压缩)时,双层结构在较小的压缩应变下出现表面失稳;(2)当基底的泊松比为负且趋于-1时,可被压缩至46%而不出现表面失稳,即可膨胀基底能有效抑制薄膜的表面失稳.本文发展的方法及主要结果可为延展性电子器件的设计提供指导.     

基于数字散斑相关实验测量的材料构型力的计算方法

材料构型力学主要研究材料中的缺陷(夹杂、空穴、位错、裂纹、塑性区等)的构型(形状、尺寸和位置)改变时,所引起的系统自由能的变化。本研究将基于数字散斑相关技术,实验测量材料试件的位移场分布,随后通过材料构型力的定义式,计算求得弹塑性材料中缺陷构型力的分布。其方法概括如下:位移场通过数字图像相关技术测得;应变及位移梯度场利用三次样条拟合获得;线弹性材料应力通过简单线弹性本构方程获取,而塑性材料的表面应力场通过Ramberg-Osgood本构方程计算求得;弹塑性应变能密度分布则由应力-应变曲线数值积分获得。该方法对普通弹性材料或者弹塑性材料均适用,可以用于各种不同的缺陷及缺陷群的材料构型力测量。     

弯曲波彩虹捕获效应及其在能量俘获中的应用

弹性波在色散关系经过设计的梯度结构中传播时会产生空间分频现象和波场能量增强现象,即不同频率的弹性波会在结构的不同位置停止向前传播并发生能量聚集,这就是弹性波彩虹捕获效应.其相关研究成果可以促进结构健康监测、振动控制以及能量俘获等领域的发展.本文通过所设计的梯度结构梁,系统地研究了弯曲波彩虹捕获效应及其在压电能量俘获中的应用.首先,利用传递矩阵法获得了梯度结构梁元胞能带结构的解析解,进而分析了弯曲波彩虹捕获效应的产生机理:不同频率的弯曲波会在不同元胞附近群速度减小到零,从而停止向前传播并发生反射;入射波和反射波的叠加,以及群速度减小带来的能量聚集,会显著增强反射处的波场能量.其次,通过有限元仿真和实验验证了弯曲波彩虹捕获效应的空间分频现象和波场能量增强现象.最后,通过有限元多物理场耦合仿真和实验,研究了粘贴PVDF压电薄膜的梯度结构梁相对于均匀梁的弯曲波能量俘获效果及其随入射波频率的变化规律.结果表明,在弯曲波彩虹捕获效应发生频带内,粘贴在梯度结构梁上的PVDF压电薄膜的输出电压约为粘贴在均匀梁相应位置处的PVDF压电薄膜的输出电压的2倍.     

压电半平面梯形夹杂的梯度特征应变问题

针对边界处自由和绝缘以及固定和绝缘两种不同的条件,分别计算分析了均匀和梯度特征应变下梯形夹杂内部和外部诱发产生的弹性场和电场,并且讨论了梯形夹杂角点处的奇异性．最后,计算了平均梯度特征应变为零时梯形夹杂内部产生的平均弹性场和电场．所得结果揭示了基体不同边界条件对诱导场的影响．     

单畴外延铁电薄膜的相结构与稳定性

本文采用动态金茨堡-朗道（DGL）方程研究了薄膜厚度与错配应变对 取向单畴外延PbTiO3（PTO）铁电薄膜相结构与稳定性的影响。结合平面内松弛应变（等效应变）、表面效应与退极化场等机电耦合边界条件,通过数值求解DGL方程获得外延单畴铁电薄膜错配应变-厚度相图和错配应变-温度相图。数值分析结果显示,由于生成的界面位错松弛了薄膜内错配应变,在理论高应变区相图与传统分析结果有较大差别,文中发现在更广的理论错配拉应变区出现稳定的四方相（c相）结构和单斜相（r相）结构。结果也显示,随着薄膜厚度的减小,表面效应与退极化效应会把顺电相扩展到更低温度区域,从而压缩稳定的铁电相存在的温度区域。     

虚位移场方法在石墨材料力学参数测量中的应用

光测实验技术在现代力学研究中得到了广泛的应用。对于材料力学参数如杨氏模量和泊松比的测量,可利用典型加载试验如拉伸试验、弯曲试验并结合光测方法(如云纹和数字图像相关技术)得到位移值,利用载荷信息和应变场信息通过计算获得相关的力学参数。本文利用虚位移场方法测量石墨材料的力学参数。结合石墨材料的三点弯曲实验,由数字图像相关法测量得到试件表面的非均匀变形场。通过选择两组不同的虚位移场,可以反算出材料的力学参数:杨氏模量和泊松比。结果表明这种方法可以有效测量石墨材料的弹性参数。该方法可望在材料力学行为检测中得到推广应用。     

一种新的宏、细观变形全场测量方法

基于光学干涉和衍射理论,提出了一种新的用于宏、细观面内变形测量的光学方法。利用三束准直相干光之间的干涉,分别得到面内位移U V和U-V的两幅耦合条纹图。借助于FFT相移技术,直接计算出全场的U和V分布。本实验系统中包含有长距离显微镜,三光束系统,数字图像采集和处理系统等。它具有光路简单,避开条纹级次的确定,直接计算得到U、V场等优点。将它用于热塑性复合材料 56层板AS4/PEEK[0/±45/90]7S三点弯曲梁位移场的宏、细观测量。得到了层间的变形规律。本文还利用三维各向异性弹性有限元法对AS4/PEEK[0/±45/90]7S三点弯曲梁进行了数值计算,得到了自由表面的位移场和应变场的数值解。定量分析了层间的变形传递规律,有助于进一步分析AS4/PEEK[0/±45/90]7S层合板宏、细观变形。     

双层薄膜与弹性梯度基底三层结构表面失稳分析

硬薄膜/软基底结构的表面失稳问题一直是柔性电子器件的难题,基于此,本文考虑了双层结构与弹性梯度基底间的界面剪切力,建立了双层薄膜/弹性梯度基底模型;利用位移协调条件,理论推导得到了双层薄膜/弹性梯度基底结构的临界应变和失稳波长的表达式并通过有限元仿真,验证了本研究解析解的有效性。在此基础上,应用此解析解进一步研究了弹性梯度基底的材料、双层薄膜结构厚度比等参数对临界应变和波长的影响。结果表明：减小器件层的厚度或者增加封装层的厚度,可以提高双层膜/弹性梯度基底结构的稳定性;当弹性梯度材料基底表面“较软”或器件层“较硬”时,器件层与基底界面的剪切力的影响较大,可以提升三层膜/基结构抵抗界面破坏的能力。本研究成果将为硬薄膜/弹性梯度基底结构的柔性电子器件的制备提供理论支撑。     

应力作用下NiTi形状记忆合金微结构演化的相场模拟及其本征应变率敏感性

基于Ginzburg-Landau动力学控制方程建立了NiTi形状记忆合金非等温相场模型,实现了对NiTi合金内应力诱导马氏体相变的数值模拟。同时将晶界能密度引入系统局部自由能密度,从而考虑多晶系统中晶界的重要作用。数值计算了单晶和多晶NiTi形状记忆合金在单轴机械载荷作用下微结构的动态演化过程和宏观力学行为,并重点研究了晶粒尺寸为60 nm的NiTi纳米多晶在低应变率下（0.000 5～15 s-1）力学行为的本征应变率敏感性。研究结果表明,单晶NiTi合金系统高温拉伸-卸载过程中马氏体相变均匀发生,未形成奥氏体-马氏体界面。而纳米多晶系统在加载阶段出现了马氏体带的形成-扩展现象,在卸载阶段出现了马氏体带的收缩-消失现象。相同外载作用过程中,NiTi单晶系统的宏观应力-应变曲线具有更大的滞回环面积,拥有更优的超弹性变形能力。计算结果显示,在中低应变率下纳米晶NiTi形状记忆合金应力-应变关系表现出较明显的应变率相关性,应变率升高导致材料相变应力提升。这一应变率相关性主要源于相场模型中外加载荷速率与马氏体空间演化速度的相互竞争关系。     

嵌入薄膜传感器悬臂梁三向力测量技术研究

以薄膜传感器悬臂梁作为等效模型,通过传感器的应变效应对三向力测量技术进行了研究。为提高薄膜传感器的应变输出响应,对悬臂梁上布放薄膜传感器的位置加设弹性结构,研究了三向力测力模型输出电压与传感器所在位置应变的关系;分析了受力位置对测力模型输出响应的影响关系,结合实验验证了其工作原理、测力模型应变输出响应与可控尺寸参数的关系。研究表明:该测力模型可实现三向力测量,各个方向最大测量误差均在9%以内,悬臂梁宽度方向x和高度方向z的交叉干涉误差分别为2.84%和3.37%;当悬臂梁自由端受力位置发生变化时,测力模型输出响应只在梁长度方向y上发生变化。     

Stress-driven interfacial instability of a bilayer structure: The interaction between free surface and interface

The morphological evolution of strained films is of technological importance to microelectronics and nanotechnology. The morphological instability of a bilayer system is analyzed, consisting of an elastic film and an elastic substrate with a misfit strain on the coherent interface. A kinetic model is derived by considering the morphological fluctuations of different perturbation amplitudes along both the free surface and the interface and the coupling effect between the film and the substrate. The couplings include the misfit strain, surface/interface energy, and surface/interface diffusion, which determine the morphological instability of the system. A quadratic dispersion relationship is established for the growth rate of the longitudinal surface and interfacial perturbations along the free surface and the interface, respectively. The propagation of the surface perturbations is revealed from the free surface to the interface, and the characteristic frequencies are identified for the initiation of the morphological instability.     

Equilibrium Configurations of Epitaxially Strained Crystalline Films: Existence and Regularity Results

Strained epitaxial films grown on a relatively thick substrate are considered in the context of plane linear elasticity. The total free energy of the system is assumed to be the sum of the energy of the free surface of the film and the strain energy. Because of the lattice mismatch between film and substrate, flat configurations are in general energetically unfavorable and a corrugated or islanded morphology is the preferred growth mode of the strained film. After specifying the functional setup in which the existence problem can be properly framed, a study of the qualitative properties of the solutions is undertaken. New regularity results for volume-constrained local minimizers of the total free energy are established, leading, as a byproduct, to a rigorous proof of the zero-contact-angle condition between islands and wetting layers.     

A theoretical analysis of the breakdown of electrostrictive oxide film on metal

Oxide films that form to protect (passivate) metal substrates from corrosive environments can be severely damaged when they are subjected to sufficient levels of electric potential. A continuum mechanics model is presented that captures the intimate electromechanical coupling of the environment and the film responsible for either growth or dissolution of the oxide. Analytical solutions, obtained for a finite-thick film experiencing a uniform electric field, illustrate the existence of a critical combination of electric field strength, initial film thickness and shape, beyond which the passivating oxide can become thin enough to undergo dielectric breakdown, or the substrate can become exposed to the corrosive environment. An experimental procedure is proposed to measure combinations of material properties required by the theoretical model to predict the lifetime of the oxide or to avoid the critical state. Illustrative numerical examples are provided to describe the morphological evolution of oxide films with a periodically wavy surface.     

Kinetic Monte Carlo simulation of strained heteroepitaxial growth with intermixing

An efficient method for the simulation of strained heteroepitaxial growth with intermixing using kinetic Monte Carlo is presented. The model used is based on a solid-on-solid bond counting formulation in which elastic effects are incorporated using a ball and spring model. While idealized, this model nevertheless captures many aspects of heteroepitaxial growth, including nucleation, surface diffusion, and long-range effects due to elastic interaction. The algorithm combines a fast evaluation of the elastic displacement field with an efficient implementation of a rejection-reduced kinetic Monte Carlo based on using upper bounds for the rates. The former is achieved by using a multigrid method for global updates of the displacement field and an expanding box method for local updates. The simulations show the importance of intermixing on the growth of a strained film. Further, the method is used to simulate the growth of self-assembled stacked quantum dots.     

The stability of elastically strained nanorings and the formation of quantum dot molecules

Self-assembled nanorings have recently been identified in a number of heteroepitaxially strained material systems. Under some circumstances these rings have been observed to break up into ring-shaped quantum dot molecules. A general non-linear model for the elastic strain energy of non-axisymmetric epitaxially strained nanostructures beyond the small slope assumption is developed. This model is then used to investigate the stability of strained nanorings evolving via surface diffusion subject to perturbations around their circumference. An expression for the fastest growing mode is determined and related to experimental observations. The model predicts a region of stability for rings below a critical radius, and also a region for larger rings which have a proportionally small thickness. The predictions of the model are shown to be consistent with the available results. For the heteroepitaxial InP on In_0.5Ga_0.5P system investigated by Jevasuwan et al. (2013), the nanorings are found to be stable below a certain critical size. This is in good quantitative agreement with the model predictions. At larger sizes, the rings are unstable. The number of dots in the resulting quantum dot molecule is similar to the mode number for the fastest growing mode. Second order terms show that the number of dots is expected to reduce as the height of the ring increases in proportion to its thickness. The strained In_0.4Ga_0.6As on GaAs nanorings of Hanke et al. (2007) are always stable and this is in accordance with the findings of the analysis. The Au nanorings of Ruffino et al. (2011) are stable as well, even as they expand during annealing. This observation is also shown to be consistent with the proposed model, which is expected to be useful in the design and tailoring of heteroepitaxial systems for the self-organisation of quantum dot molecules.     

饱和蒸发和热毛细力作用下低雷诺数液膜在波纹竖壁上的流动

考虑表面蒸发压力和热毛细力作用情况下,对饱和蒸发状态下低雷诺数自由降落液膜在小波幅正弦型波纹壁面上的流动进行理论分析。对控制微分方程及边界条件进行量纲一化并引入流函数,对微分方程及边界条件进行摄动展开,得到了这种情况下液膜流动的简化分析模型,求出了近似解析解。讨论了壁面波纹、表面张力、蒸发压力、热毛细力对液膜流动的影响。研究表明:液膜的波动幅度随蒸发强度和热毛细力的增大而增大;液膜波动与壁面波纹的相位差随蒸发强度增大而增大,随热毛细力增大而减小。     

Strain Relaxation in an Alloy Film with a Rough Free Surface

We study strain relief by surface roughness and composition variation in a stressed alloy film. Instead of using common perturbation techniques, we derive a rigorous relaxation formula based on the energy approach in the case of slightly undulating surface and fluctuating composition. We do not require any a priori assumption of elastic isotropy or identical material properties between film and substrate in deriving our result. We show that the change of elastic energy is negative, giving rise to energy relief due to the presence of free surface. We apply our result to the study of compositional and morphological instabilities of a stressed thin layer with a free surface. The critical wave number of instability is determined by the competition between the destabilizing influence of elastic strain energy and the stabilizing influence of chemical and surface energies.     

Post-buckling evolution of wavy patterns in trapezoidal film/substrate bilayers

This paper investigates quantitatively the post-buckling response of irregular wrinkles in a trapezoidal film/substrate bilayer. The geometric gradient can change the wrinkling profile to create ribbed and graded structural patterns with variations of wave direction, amplitude and wavelength. The tapered angle and edge dimension are examined numerically using a nonlinear shell/solid coupled finite element model that incorporates a path-following continuation technique, which explores their influences on secondary bifurcations, localization and surface mode transition. For instance, the competition between plate-like and beam-like post-buckling behavior is discussed. An analysis of graded amplitude is also provided based on Fourier envelope equations of beam/foundation models, which gives an insightful understanding of these fading wrinkles that differ from the ones usually observed in rectangular geometric cases. The results of this work can be used to guide the design of geometrically gradient film/substrate systems to achieve desired wavy instability patterns.     

Thermo-mechanical solution of film/substrate systems under local thermal load and application to laser lift-off of GaN/sapphire structures

Film/substrate structures may undergo a localized thermal load, which can induce stresses, deformation and defects. In this paper, we present the solutions of temperature and stresses in a film/substrate structure under a local thermal load on the film surface. Then, the generalized Stoney formula, which connects the curvature of deformation and the stress field is obtained. The present solution takes into account the non-uniformity of the temperature field both in the width and thickness directions of the film. The thermo-mechanical solution is applied to the analysis of the temperature distribution, stresses, and damage of a GaN/sapphire system during the laser lift-off (LLO) process. It is shown that the laser with the Gaussian distribution of energy density causes much smaller tensile stresses at the edge of the heated area in the film than the laser with the uniform distribution of energy density, and thus can avoid damage to the GaN films separated from the substrate.     

Experimental study of filmwise condensation considering wave initiation

The problem of wavy filmwise condensation was recently investigated using both analytical and numerical concepts. The analytical solution assumed the film profile as a simple sinwave moving with a constant phase velocity and the problem was reduced to a closed integral form. While the numerical solution considers variations in the wave length and amplitude of the proposed sin wavy free surface. In order to find out the validity of the suggested theoretical models, an experimental test rig was installed, where the experimental parameters were controlled to generate data in the range of 40 <Re < 200. The results of this work demonstrate the possibility of accurate simulation of long vertical condenser using laboratory set of a reasonable dimensions. This simulation yields better results when the thermal and hydrodynamical average quantities at the entrance are properly corrected. Good agreement was found by comparing the generated experimental data, together with others data, with the analytical and numerical models. Imperical formulas for heat transfer were considered and differences between the correlation coefficients was found to be mainly due to unproper definition of the validity ranges of the used dimensionless groups. More realistic definition for the range of validity was introduced with these imperical formulas.