CN108716959B - 有效预测压电薄膜与梯度非均匀基底界面应力分布的方法 - Google Patents
有效预测压电薄膜与梯度非均匀基底界面应力分布的方法 Download PDFInfo
- Publication number
- CN108716959B CN108716959B CN201810310838.9A CN201810310838A CN108716959B CN 108716959 B CN108716959 B CN 108716959B CN 201810310838 A CN201810310838 A CN 201810310838A CN 108716959 B CN108716959 B CN 108716959B
- Authority
- CN
- China
- Prior art keywords
- interface
- film
- substrate
- piezoelectric film
- stress
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L5/00—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
- G01L5/0047—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes measuring forces due to residual stresses
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及一种有效预测压电薄膜与非均匀基底间界面应力分布的方法,首先利用实验测量梯度非均匀基底的几何和材料参数,获得压电薄膜的电弹参数,利用力学模型获得薄膜/非均匀基底系统接触力学问题的控制方程,利用边界配置法求解控制方程,即可得到薄膜/非均匀基底系统界面应力分布的信息,为薄膜/非均匀基底系统界面的优化设计以及基底梯度特性设计提供有效依据。
Description
技术领域
本发明属于结构无损检测领域,尤其涉及一种有效预测压电薄膜与梯度非均匀基底界面应力分布的方法。
技术背景
作为新材料的典型代表,梯度非均匀材料在机械、土木工程、航空航天、医学等领域得到了广泛的应用。与传统的均匀材料相比,梯度非均匀材料能够将异质材料的优势结合,并极大地减少由于材料适配引起的应力集中,并提高结构的表面抗磨损性能等,因而受到了各国科学界和工程界的高度重视。然而由于梯度非均匀材料的制备工艺、受载损伤等因素,使得该类材料不可避免地产生缺陷和破坏。因此,探索和发展梯度非均匀材料的无损检测技术成为重要的研究方向。
由于压电薄膜、柔性压电器件等都具有优异的电弹性质,将其植入或者黏贴在结构内部或者表面,可以实现结构的无损检测和实时控制,这无疑为梯度非均匀材料无损检测的发展提供了新的启示。薄膜/基底系统界面应力传递作为无损检测的理论基础,势必受压电薄膜的电弹性质、梯度非均匀材料的梯度变化规律、薄膜与梯度基底的粘结情况的影响。因此,一种预测压电薄膜与梯度非均匀基底界面应力响应的方法亟需建立。
发明内容
为了解决上述技术问题,本发明提供一种预测压电薄膜与梯度非均匀基底界面应力分布预测方法。
为达到上述目的,本发明采用的技术方案为:
一种有效预测压电薄膜与梯度非均匀基底界面应力分布的方法,其特征是,步骤如下:
步骤1)获取压电薄膜/梯度非均匀基底系统材料参数的表征
步骤1.1)非均匀基底弹性参数和梯度变化规律的测量和表征
采用接触压入技术获得梯度非均匀基底沿厚度方向的弹性模量分布,并用指数函数进行拟合,拟合公式为:
μ(z)=μ1exp(γz)
其中μ1为非均匀基底表面的剪切模量,z为从非均匀基底从表面到内部的厚度值,μ(z)为基底厚度为z处的剪切模量值,γ表征梯度基底的不均匀度;
步骤1.2)获取压电薄膜材料参数
压电薄膜通常为横观各向同性材料,它的弹性参数、压电参数、介电常数可通过查询文献得到,或者将制备好的压电薄膜采用已经成熟直接测量法,如悬臂法和间接测量法,如复合谐振法方法获得;
1.3)压电薄膜与非均匀基底的界面力学分析
压电薄膜与梯度非均匀基底间的界面根据结合性能定义为三种形式,包括应力连续的界面,有限厚度的胶层,以及界面处存在界面脱黏的情况,界面性质取决于薄膜与基底间的结合性能;
1.31)对于应力连续的界面,界面关系为
其中和σxz(x)分别为薄膜和基底界面上的切应力,tl为薄膜左端的x坐标,为膜内正应力,hf为压电薄膜的厚度;
1.32)对于有限厚度的胶层界面,界面关系为
其中μb和分别为胶层的剪切模量和x方向正应变,其中等于胶层上下表面的位移差与其剪切模量的比值,即
1.33)对于存在界面脱黏的情况,薄膜内部的正应力与切应力满足如下关系:
其中σd为待定常数,即脱黏区域的膜内正应力相等;dl和dr为脱黏区域两端的x坐标,tl和tr为薄膜两端的x坐标;
2)考虑具体的载荷类型,建立力学模型得到薄膜/梯度非均匀基底系统的控制方程
2.1)对压电薄膜进行受力分析,利用界面切应力表示压电薄膜的横向正应变,压电薄膜应力和应变关系为
其中Ez为施加的电场强度,Ef为等效弹性参数,ef为等效压电材料参数;
2.2)考虑基底的梯度属性,给出表面切应力下非均匀基底表面的x方向的位移,梯度非均匀基底表面的x方向的位移为积分形式的表达式:
其中r为积分变量,G为广义表达式,a为薄膜长度的一半,为待求函数;
如果薄膜和基底间为连续应力界面情况,薄膜/基底系统只承受电场作用,则界面相容方程则为
2.3)根据载荷和界面特性,给出非均匀基底与压电薄膜应变的相容关系,即压电薄膜的横向正应变和非均匀基底表面位移的关系,从而得到表征压电薄膜/非均匀基底系统界面力学特性的控制方程
薄膜/基底系统的控制方程为具有奇异核的积分方程,其中未知函数为界面的切应力分布r为积分变量,;
3)利用边界配置法将控制方程离散化,将积分方程离散成为线性代数方程组;
4)求解线性代数方程,从而得到界面应力分布信息。
本发明的有益效果是,提出了一种压电薄膜/非均匀基底系统界面应力预测的方法,填补了薄膜/基底结构界面应力预测方面的空白,为薄膜/非均匀基底系统界面的优化设计以及基底梯度特性设计提供有效依据。
附图说明
图1为完美黏附的应力连续界面压电薄膜/非均匀基底系统。
图2是具有胶层的压电薄膜/梯度非均匀基底系统。
图3是存在界面脱胶的薄膜/基底系统
图4是电场作用下,完美黏附的应力连续界面压电薄膜/非均匀基底界面的切应力分布。
图5是是电场作用下,完美黏附的应力连续界面压电薄膜/非均匀基底界面的横向正应力分布。
具体实施方式
1)压电薄膜/梯度非均匀基底系统材料参数的表征
1.1)非均匀基底弹性参数和梯度变化规律的测量和表征
采用接触压入技术获得梯度非均匀基底沿厚度方向的弹性模量分布,并用指数函数进行拟合,拟合公式为:
μ(z)=μ1exp(γz)
其中μ1为非均匀基底表面的剪切模量,z为从非均匀基底从表面到内部的厚度值,μ(z)为基底厚度为z处的剪切模量值,γ表征梯度基底的不均匀度。
1.2)通过查询文献,得到压电薄膜的弹性参数、压电参数、介电常数等。
1.3)压电薄膜与非均匀基底的界面力学分析
本发明中压电薄膜与梯度非均匀基底间的界面既可以处理成应力连续的界面,亦可以将薄膜与基底间的界面考虑成有限厚度的胶层,并可以分析界面处存在界面脱黏的情况;界面性质的选择取决于薄膜与基底间的结合性能和设计要求。
i)对于应力连续的界面(如图1所示),可考虑
其中和σxz(x)分别为薄膜和基底界面上的切应力,tl和tr为薄膜两端的x坐标,为膜内正应力,hf为薄膜厚度;
ii)对于有限厚度的胶层界面(如图2所示),界面关系为
其中μb和分别为胶层的剪切模量和x方向正应变,其中等于胶层上下表面的位移差与其剪切模量的比值,即
iii)对于存在界面脱黏的情况(如图3所示),薄膜内部的正应力与切应力满足如下关系:
其中hf为压电薄膜的厚度,σd为待定常数,即脱黏区域的膜内正应力相等,dl和dr为脱黏区域两端的x坐标;
2.1)压电薄膜应力和应变关系为
其中Ez为施加的电场强度,Ef为等效弹性参数,ef为等效压电材料参数。
具体参考(Jin,C.,Wang,X.,2011.Analytical modelling of theelectromechanical behaviour of surface-bonded piezoelectric actuatorsincluding the adhesive layer.Eng.Fract.Mech.78,2547-2562.)的论文。
得到梯度非均匀基底表面的x方向的位移为积分形式的表达式:
其中r为积分变量,G为广义表达式,a为薄膜长度的一半。
具体参考(Chen,P.,Chen,S.,2013.Partial slip contact between a rigidpunch with an arbitrary tip-shape and an elastic graded solid with a finitethickness.Mech.Mater.59,24-35.)的论文。
如果薄膜和基底间为完美应力连续界面情况,薄膜/基底系统只承受电场作用,则界面相容方程则为
即可得到控制方程
3)将控制方程离散化,将积分方程离散成为线性代数方程组。
选择坐标值的离散点x=xi(i=0,...,N),其中xi为取值范围内的第i个取值;
根据薄膜/基底结构的界面力学属性,可将界面切应力离散表达为其中An(n=0,...,N)为待定常数,Tn(·)为第n次的多项式展开项,为求和符号。
将坐标值和界面切应力的离散表达式带入控制方程,则包含未知函数的积分方程变成仅含有未知常数An的代数方程组。
[K]{a}={b}
其中[K]为已知参数矩阵,{a}={A0,A2,...,AN}T为未知向量,{b}为已知参数向量。
4)求解线性代数方程,从而得到界面切应力和横向正应力分布信息,如图4、5所示。
4.1)线性代数方程[K]{a}={b}的求解,既可选择Matlab或者Fortran编辑器等商用软件提供的命令直接计算,也可以根据高斯迭代、三角分解法等方法编制数值计算程序进行计算,线性代数方程组[K]{a}={b}求解后的值即为待定常数向量{a}={A0,A2,...,AN}T;
4.2)进一步利用和界面应力的关系,即可得切应力和横向正应力。
Claims (6)
1.一种有效预测压电薄膜与梯度非均匀基底界面应力分布的方法,其特征是,步骤如下:
步骤1)获取压电薄膜/梯度非均匀基底系统材料参数的表征
步骤1.1)非均匀基底弹性参数和梯度变化规律的测量和表征
采用接触压入技术获得梯度非均匀基底沿厚度方向的弹性模量分布,并用指数函数进行拟合,拟合公式为:
μ(z)=μ1exp(γz)
其中μ1为非均匀基底表面的剪切模量,z为非均匀基底从表面到内部的厚度值,μ(z)为基底厚度为z处的剪切模量值,γ表征梯度基底的不均匀度;
步骤1.2)获取压电薄膜材料参数
压电薄膜的材料参数通过查询文献得到,或者将制备好的压电薄膜采用已经成熟的直接测量法和间接测量法获得;
1.3)压电薄膜与非均匀基底的界面力学分析
压电薄膜与梯度非均匀基底间的界面根据结合性能定义为三种形式,包括应力连续的界面,有限厚度的胶层,以及界面处存在界面脱黏的情况,界面性质取决于薄膜与基底间的结合性能;
2)考虑具体的载荷类型,建立力学模型得到薄膜/梯度非均匀基底系统的控制方程;
2.1)对压电薄膜进行受力分析,利用界面切应力表示压电薄膜的横向正应变;
2.2)考虑基底的梯度属性,给出表面切应力下非均匀基底表面的x方向的位移,如果薄膜和基底间为连续应力界面情况,薄膜/基底系统只承受电场作用;
2.3)根据载荷和界面特性,给出非均匀基底与压电薄膜应变的相容关系,即压电薄膜的横向正应变和非均匀基底表面位移的关系,从而得到表征压电薄膜/非均匀基底系统界面力学特性的控制方程,薄膜/基底系统的控制方程为具有奇异核的积分方程,其中未知函数为界面的切应力分布;
3)利用边界配置法将控制方程离散化,将积分方程离散成为线性代数方程组;
4)求解线性代数方程,从而得到界面应力分布信息。
2.根据权利要求1所述的一种有效预测压电薄膜与梯度非均匀基底界面应力分布的方法,其特征是,步骤1.2)中所述的压电薄膜材料参数至少包括,弹性参数、压电参数和介电常数。
3.根据权利要求1所述的一种有效预测压电薄膜与梯度非均匀基底界面应力分布的方法,其特征是,步骤1.3)压电薄膜与非均匀基底的界面力学分析,参照界面的不同形式,具体步骤为:
1.31)对于应力连续的界面,界面关系为
其中和σxz(x)分别为薄膜和基底界面上的切应力,tl为薄膜左端的x坐标,为膜内正应力,hf为压电薄膜的厚度;
1.32)对于有限厚度的胶层界面,界面关系为
其中μb和分别为胶层的剪切模量和x方向正应变,其中等于胶层上下表面的位移差与其剪切模量的比值,即
1.33)对于存在界面脱黏的情况,薄膜内部的正应力与切应力满足如下关系:
其中hf为压电薄膜的厚度;σd为待定常数,它等于脱黏区域的膜内正应力;dl和dr为脱黏区域两端的x坐标,tl和tr为薄膜两端的x坐标。
4.根据权利要求1所述的一种有效预测压电薄膜与梯度非均匀基底界面应力分布的方法,其特征是,步骤2.1)中,所述的压电薄膜切应力和横向正应变关系为
其中Ez为施加的电场强度,Ef为等效弹性参数,ef为等效压电材料参数。
5.根据权利要求1所述的一种有效预测压电薄膜与梯度非均匀基底界面应力分布的方法,其特征是,步骤2.2)中,所述的表面切应力下梯度非均匀基底表面的x方向的位移为:
其中r为积分变量,G为广义表达式,a为薄膜长度的一半,为待求函数;
如果薄膜和基底间为连续应力界面情况,薄膜/基底系统只承受电场作用,则界面相容方程则为
6.根据权利要求1所述的一种有效预测压电薄膜与梯度非均匀基底界面应力分布的方法,其特征是,步骤2.3)中,所述表征压电薄膜/非均匀基底系统界面力学特性的控制方程为
该方程为具有奇异核的积分方程,其中未知函数为界面的切应力分布r为积分变量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810310838.9A CN108716959B (zh) | 2018-04-09 | 2018-04-09 | 有效预测压电薄膜与梯度非均匀基底界面应力分布的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810310838.9A CN108716959B (zh) | 2018-04-09 | 2018-04-09 | 有效预测压电薄膜与梯度非均匀基底界面应力分布的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108716959A CN108716959A (zh) | 2018-10-30 |
CN108716959B true CN108716959B (zh) | 2019-12-06 |
Family
ID=63898810
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810310838.9A Active CN108716959B (zh) | 2018-04-09 | 2018-04-09 | 有效预测压电薄膜与梯度非均匀基底界面应力分布的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108716959B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110705076A (zh) * | 2019-09-25 | 2020-01-17 | 哈尔滨理工大学 | 一种求解具有任意属性的功能梯度压电材料断裂问题的方法 |
CN111797503B (zh) * | 2020-06-05 | 2023-10-20 | 中国飞机强度研究所 | 一种控制压电薄膜连接器脱离屈曲的方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2812082B1 (fr) * | 2000-07-20 | 2002-11-29 | Vai Clecim | Rouleau de mesure de planeite |
JP2004317403A (ja) * | 2003-04-18 | 2004-11-11 | Alps Electric Co Ltd | 面圧分布センサ |
JP2012088827A (ja) * | 2010-10-18 | 2012-05-10 | Mitsubishi Heavy Ind Ltd | 負荷イナーシャ推定方法及び制御パラメータ調整方法 |
CN103245437B (zh) * | 2012-02-13 | 2017-02-08 | 付康 | 一种确定非线性薄膜应力的系统与方法 |
CN104280155B (zh) * | 2013-07-14 | 2019-03-12 | 付康 | 一种确定多层膜薄膜应力的系统与方法 |
CN204115936U (zh) * | 2014-09-29 | 2015-01-21 | 西北农林科技大学 | 一种基于多点pvdf压电薄膜传感器的雨滴打击动能测试系统 |
CN106248280B (zh) * | 2016-08-22 | 2018-03-20 | 东南大学 | 一种导电薄膜材料残余应力的在线测量方法及测量装置 |
-
2018
- 2018-04-09 CN CN201810310838.9A patent/CN108716959B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108716959A (zh) | 2018-10-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ghosh et al. | Interfacial debonding analysis in multiple fiber reinforced composites | |
CN104112042B (zh) | 一种基于分子动力学的金刚石涂层膜基界面结合强度测量方法 | |
Van Do et al. | Thermal buckling analysis of cracked functionally graded plates | |
CN108716959B (zh) | 有效预测压电薄膜与梯度非均匀基底界面应力分布的方法 | |
Boll et al. | Miniaturized flexible interdigital sensor for in situ dielectric cure monitoring of composite materials | |
Tay et al. | Mesh design in finite element analysis of post-buckled delamination in composite laminates | |
CN110082029B (zh) | 一种pvdf压电薄膜传感器动态力学特性标定方法 | |
CN106096163A (zh) | 一种碳纤维复合材料切削模型的建立方法 | |
JP3809374B2 (ja) | 応力−ひずみ関係シミュレート方法および除荷過程における降伏点を求める方法 | |
Cavallo et al. | “Soft Si”: Effective stiffness of supported crystalline nanomembranes | |
Liprandi et al. | A theoretical-numerical model for the peeling of elastic membranes | |
Intarit et al. | Elastic layer under axisymmetric indentation and surface energy effects | |
Mekonnen et al. | Effects of size and location of initial delamination on post-buckling and delamination propagation behavior of laminated composites | |
Xu et al. | Surface tension measurement from the indentation of clamped thin films | |
Crocombe et al. | A unified approach for predicting the strength of cracked and non-cracked adhesive joints | |
Holcomb et al. | Inelastic constitutive properties and shear localization in Tennessee marble | |
Zhao | Three-parameter approach for elastic–plastic stress field of an embedded elliptical crack | |
Wan et al. | Evaluation of interfacial friction characteristics in multi-mode ultrasonic vibration assisted micro-extrusion process | |
Jiang et al. | A dynamic model of laminated material extrusion additive manufacturing plate with the property of orthogonal anisotropy | |
Attia et al. | Numerical analysis of frictionless nano-contact problems based on surface elasticity | |
Kam et al. | Material characterization of laminated composite materials using a three-point-bending technique | |
Jeong et al. | Effective mechanical properties of micro/nano-scaleporous materials considering surface effects | |
CN111141602B (zh) | 利用间接拉伸试验同步确定沥青混合料拉、压模量的方法 | |
Sapsathiarn et al. | Electro-mechanical load transfer from a fiber in a 1–3 piezocomposite with an imperfect interface | |
Heinrich et al. | On the identification of cohesive parameters for printed metal-polymer interfaces |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |