CN108256281A - 一种考虑搭接界面形貌及搭接物梯度性质的强度预测方法 - Google Patents
一种考虑搭接界面形貌及搭接物梯度性质的强度预测方法 Download PDFInfo
- Publication number
- CN108256281A CN108256281A CN201810250549.4A CN201810250549A CN108256281A CN 108256281 A CN108256281 A CN 108256281A CN 201810250549 A CN201810250549 A CN 201810250549A CN 108256281 A CN108256281 A CN 108256281A
- Authority
- CN
- China
- Prior art keywords
- parameter
- model
- overlap joint
- glue
- field variable
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及一种考虑搭接界面形貌及搭接物梯度性质的强度预测方法。首先得到拟合后的不平的搭接面形貌后,给定坐标系,算出粘结层厚度的表达式。然后用ABAQUS进行建模,根据模型坐标系确定被粘物梯度参数和胶层参数。搭接物梯度材料参数和粗糙界面的粘结层的参数赋予模型需要通过ABAQUS中的用户子程序USDFLD来实现。材料参数全部赋予模型后创建分析步,添加相应的约束和边界条件,创建工况计算后就可以提取最粘结体系的强度。
Description
技术领域
本发明属于粘结搭接领域,尤其涉及一种考虑搭接界面形貌及搭接物梯度性质的强度预测方法。
背景技术
随着胶黏剂的快速发展,粘结结构广泛运用于机械、建筑、航空航天、医学等领域。与传统的钉子连接、铆接、螺栓连接等传统方法相比,胶接技术有着不受搭接物结合形状的限制,应力集中小,胶接结构质量轻等,连接密封性好,成本较低等特点,使之越来越受到青睐。
在具体应用中,胶接结构的破坏往往都发生在粘结界面上,所以对界面强度的研究具有重要意义。而对于粘接结构强度的研究,早期的学者通过各种专业力学试验设备进行研究,需要耗费大量的人力物力,甚至现在有些实验条件还没法达到。所以研究者们开始重视理论模型的建立,但是往往由于构建形状,载荷,边界条件的复杂性,很多时候想找到一个解析解很困难。现在计算机技术和数值计算方法的飞速发展,使得粘结结构的力学性能研究进入崭新空间。
粘结结构的强度受到粘结层参数、被粘物材料参数,以及粘结构件的几何参数所影响。现在对粘结搭接结构的的研究几乎考虑的搭接面都是平滑的,但是在实际中的搭接面很可能是凹凸不平的。搭接强度对胶层厚度很敏感,搭接粗糙不平导致搭接面的胶层厚度是梯度变化的,此时用平滑界面搭接面的强度来预测不平滑的情况就很不准确。很多强度预测模型考虑的搭接物都是均质的或者是复合材料,但是在实际应用中,搭接物材料性质有时是梯度变化的,比如有些航空航天材料的梯度涂层,生物医学领域里面的牙齿骨骼。
因此,一种考虑搭接界面形貌及搭接物梯度性质的强度预测方法亟需建立。
发明内容
本发明提供一种考虑搭接界面形貌及搭接物梯度性质的强度预测方法,弥补对于非平滑搭接界面和梯度搭接物的强度预测的空白。
为达到上诉目的,本发明采用的技术方案为:
一种考虑搭接界面形貌及搭接物梯度性质的强度预测方法,其特征在于,包括以下步骤:
(1)用直线和曲线近似拟合搭接界面的形貌,给定一个坐标系,计算出拟合过的上下搭接界面形状的函数f1(x,y)和f2(x,y);
(2)按照给定的坐标系,用有限元对近似的搭接模型进行建模;
(3)在ABAQUS软件中通用的材料赋值选项中,只有针对均质材料或者各项异性材料,没有材料梯度参数的赋值选项,所以我们对梯度参数的赋值需要调用ABAQUS的用户子程序USDFLD,目的是把梯度参数以场的形式赋予模型。首先需要编写一个Fortran文件,这文件作用是定义场变量为X/Y/Z坐标,场变量可以是一维二维三维。对于一维场变量就只需定义一个场变量1,如果梯度参数是随X轴变化,那就只需要定义场变量1是X轴坐标,如果梯度参数是随Y轴或者Z轴变化,那就定义场变量1是Y轴坐标或者Z轴坐标;如果二维场变量家具需要定义两个场变量,例如梯度参数在X0Y平面内变化就需要定义场变量1和场变量2分别为X坐标和Y坐标;如果是三维场变量就需要三个场变量分别为X、Y、Z坐标。所以首先根据模型参数的变化来定义相应的场变量,然后通过ABAQUS的用户子程序USDFLD的接口提交编写好的Fortran文件就实现好了场变量的定义,最后根据已经定义好的场变量和模型的坐标系位置算出把每个网格节点的参数赋予模型,就完成了梯度材料参数的设置;
(3.1)梯度材料搭接物的弹性模量或者泊松比大多数情况下都是一位变化或者二维变化,梯度变化规律为E(x,y),υ(x,y);
(3.2)胶层材料参数的确定。在实际粘结体系工程计算中,内聚力模型能够描述界面的开裂过程,其中双线性内聚力模型受到了很多研究者的青睐。界面分离位移大到一定程度后,界面层就开始损伤,需要选择相应的初始损伤准则和断裂准则。我们采用了最常用的二次名义应力损伤判据和判断界面是否断裂的幂指数准则。在有限元模拟中,需要输入胶层的内聚力参数初始法向和切向的初始刚度、总断裂能。胶层的初始法向和切向的初始刚度、总断裂能、分离强度的表达式分别为:
界面法向和切向初始刚度:
总断裂能
分离强度
其中E和G分别为杨氏模量和剪切模量;Γ0为本征断裂能;U可以看成为塑性区内单位体积胶黏剂所耗散的塑性能;为裂纹尖端塑性区估计值;σf为宏观分离强度。
从上述内容可以看出胶层的初始法向和切向的初始刚度、总断裂能、分离强度对胶层的厚度t都有依赖。而对于粗糙的搭接面,胶层的厚度 t(x,y)=f1(x,y)-f2(x,y)为梯度变化的,从而导致胶的性质呈现梯度变化。把梯度的胶层厚度带入上面初始刚度、总断裂能、分离强度表达式,就可以得到粗糙搭接界面的粘结层的内聚力参数;
(4)模型的材料参数赋值完后,就可以创建分析步;然后根据实际情况,给模型添加边界条件,外载等条件;对模型进行单元的划分,被粘物采用四节点平面应变单元。而粘结层采用单层四节点内聚力单元,最后创建工况就可以开始计算;
(5)强度的提取,工况计算完成后,提取出模型的力位移曲线,力位移曲线的最高点所对应的FP就是作为衡量粘结体系总体强度的参考值。
本发明的有益效果是:考虑了搭接界面的形貌特征和搭接材料的梯度特性,填补了考虑搭接界面形貌和搭接物梯度特性的强度预测的空白,为搭接界面的优化设计以及搭接物梯度特性的设计提供有效依据。
附图说明
图1为双线性内聚力模型的T-S曲线图。
图2是实例1模型示意图。
图3是实例2模型示意图。
图4是实例3模型示意图。
图5是3个实施例的力位移载荷图。
具体实施方式
实施例1
下面结合附图和具体实施例对本发明做进一步说明。下面所描述的只是发明的一部分实施例。对于本领域的普通技术人员来说,不需要付出创作性的劳动,只需要通过改变具体搭接界面的形状以及搭接物梯度特性就能获得相应的粘结体系的强度。
1)图2是一个对实际中不平滑搭接面近似拟合后的单搭接模型。搭接面的上表面是4段等长的线段组成,每条线段的斜率k绝对值为0.4/12.5,搭接面的下表面是平滑的直线,上下搭接物的长L为200mm,厚度h为5mm,搭接长度l为50mm。把直角坐标系原点建在搭接面左下角,则可计算出上下搭接面的形状函数为:
下搭接面的形状函数为:f2(x)=0(0≤x≤50)。
2)然后用ABAQUS有限元商业软件,按照图2模型尺寸和坐标系进行二维建模。
3)接下来我们对模型进行材料属性定义。本实例中被粘物的梯度参数和粗糙搭接界面粘结层的参数都只是随X轴一维变化,首先需要编写一个 Fortran文件定义一个场变量1为X坐标。
3.1)如图2所示,本实例中下搭接物中含搭接区域的这一截的杨氏模量沿着X方向梯度变化,梯度变化规律为泊松比为 v1=0.29;上搭接物和下搭接物其余部分都是均质材料杨氏模量E1(x)=209GPa,泊松比v1=0.29。梯度部分材料参数根据梯度层网格的大小和坐标系位置,利用 MATLAB软件算出每个节点的参数,把这些参数赋予模型的相应梯度材料部分就完成赋值。
3.2)胶层材料参数的确定。双线性内聚力的张力和开裂位移的关系曲线如图1所示。图1中(a)为法向内聚应力和法向分离量的关系;(b)为切向内聚应力和切向分离量的关系。图中σ和τ表示拉伸和剪切方向的内聚应力值,而σu和τu分别是纯Ⅰ型内聚力模型(拉伸)和纯Ⅱ型内聚力模型(剪切)的内聚应力的峰值,称作分离强度。δc和δm分别代表临界位移和完全失效位移.k是内聚力模型中的初始刚度值,下标"Ⅰ"和"Ⅱ"分别表示拉伸和剪切方向的对应值。胶层的初始法向和切向的初始刚度、总断裂能、分离强度对胶层的厚度t 都有依赖。而对于粗糙的搭接面,胶层的厚度t(x,y)=f1(x,y)-f2(x,y)为梯度变化的,从而导致胶的性质呈现梯度变化。把梯度的胶层厚度带入上面初始刚度、总断裂能、分离强度表达式,就可以得到粗糙搭接界面的粘结层的内聚力参数。
本实例中粘结层采用环氧树脂胶黏剂Hysol EA9361,他的材料属性如表1 所示:
表1
接下来需要确定胶层的内聚力参数初始法向和切向的初始刚度、总断裂能、分离强度。胶层的初始法向和切向的初始刚度、总断裂能、分离强度的表达式分别为:
在本实例中胶层的厚度t(x)=f1(x)-f2(x),即
把梯度的胶层厚度表达式带入上面初始刚度、总断裂能、分离强度表达式,最后表中胶层的材料参数代入就可以得到本实例的粘结层的随坐标x变化的内聚力参数。得到粘结层内聚力参数后,和上面搭接物梯度参数赋值一样,根据粘结层网格的大小和坐标系位置,利用MATLAB软件算出每个节点的参数,把这些参数赋予模型的相应粘结层部分就完成赋值。
4)模型的材料参数赋值完后,就可以创建分析步;接下来按照图1所给出边界条件,给模型施加相应的边界条件和外载荷;然后对模型进行单元的划分,被粘物采用四节点平面应变单元,而粘结层采用单层四节点内聚力单元,最后创建工况就可以开始计算。
5)强度的提取。工况计算完成后,输出上面参考点处的U1和RT1,就可以合成本粘结体系的力位移曲线,曲线最高点的所对应的力,即最大拉脱离,就可以代表这个粘结体系的强度。
为了反应搭接界面和搭接物梯度特性会对粘结体系的强度产生影响,突出本发明的意义,发明人还做了两个实例。
实施例2,如图3所示把实例1中的被粘物梯度特性去掉,弹性模量都取 209GPa,其他条件不变。
实实例3,如图4所示把实例2中搭接界面的不平滑去掉,改成平滑界面,其他条件不变。
图5是三个实施例的力位移曲线图,从图中可以看出,同时考虑搭接面形貌和被粘物梯度性质的最大拉脱离最大,其次是只考虑不光滑搭接面的最大拉脱力,最小的是光滑搭接面和均质搭接物的最大拉脱力,说明搭接面的形貌和搭接物的梯度特性都能对粘结体系的强度有很大影响,从而就更直接体现本发明的意义。
Claims (2)
1.一种考虑搭接界面形貌及搭接物梯度性质的强度预测方法,其特征在于,包括以下步骤:
步骤1)用直线和曲线近似拟合搭接界面的形貌,给定一个坐标系,计算出拟合过的上下搭接界面形状的函数f1(x,y)和f2(x,y);
步骤2)按照给定的坐标系,用有限元对近似的搭接模型进行建模;
步骤3)在ABAQUS软件中把梯度参数以场的形式赋予模型,根据模型参数的变化来定义相应的场变量,然后通过ABAQUS的用户子程序USDFLD的接口提交编写好的Fortran文件实现场变量的定义,最后根据已经定义好的场变量和模型的坐标系位置算出把每个网格节点的参数赋予模型,完成梯度材料参数的设置;
步骤3.1)梯度材料搭接物的弹性模量或者泊松比大多数情况下都是一维变化或者二维变化,梯度变化规律为E(x,y),υ(x,y);
步骤3.2)粗糙搭接面胶层材料参数的确定,以双线性内聚力模型描述界面的开裂过程,即张力和开裂位移的关系:
当界面分离位移大到一定程度后,界面层就开始损伤,需要选择相应的初始损伤准则和断裂准则;采用最常用的二次名义应力损伤判据和判断界面是否断裂的幂指数准则,在有限元模拟中,输入胶层的内聚力参数初始法向和切向的初始刚度、总断裂能、分离强度,胶层的初始法向和切向的初始刚度、总断裂能、分离强度的表达式分别为:
界面法向初始刚度kⅠ和切向初始刚度kⅡ:
总断裂能
分离强度
其中E和G分别为杨氏模量和剪切模量;Γ0为本征断裂能;U为塑性区内单位体积胶黏剂所耗散的塑性能;为裂纹尖端塑性区估计值;σf为宏观分离强度;
胶层的初始法向和切向的初始刚度、总断裂能、分离强度对胶层的厚度t都有依赖,而对于非平滑的搭接面,胶层的厚度t(x,y)=f1(x,y)-f2(x,y)为梯度变化的,从而导致胶的性质呈现梯度变化;把梯度的胶层厚度带入上面初始刚度、总断裂能、分离强度表达式,就可以得到粗糙搭接界面的胶层的材料参数;
步骤4)模型的材料参数赋值完后创建分析步;然后根据实际情况,给模型添加边界条件、外载条件对模型进行单元的划分,被粘物采用四节点平面应变单元,而粘结层采用单层四节点内聚力单元,最后创建工况开始计算;
步骤5)强度的提取,工况计算完成后,提取出模型的力位移曲线,力位移曲线的最高点所对应的FP作为衡量粘结体系总体强度的参考值。
2.根据权利要求1所述的一种考虑搭接界面形貌及搭接物梯度性质的强度预测方法,其特征在于,步骤3)中场变量的定义为:
首先编写一个Fortran文件,定义场变量为X/Y/Z坐标,场变量是一维二维或三维;对于一维场变量就只需定义一个场变量1,如果梯度参数是随X轴变化,那就只需要定义场变量1是X轴坐标,如果梯度参数是随Y轴或者Z轴变化,那就定义场变量1是Y轴坐标或者Z轴坐标;如果是二维场变量需要定义两个场变量,定义场变量1和场变量2分别为X坐标和Y坐标;如果是三维场变量就需要三个场变量分别为X、Y、Z坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810250549.4A CN108256281B (zh) | 2018-03-26 | 2018-03-26 | 一种考虑搭接界面形貌及搭接物梯度性质的强度预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810250549.4A CN108256281B (zh) | 2018-03-26 | 2018-03-26 | 一种考虑搭接界面形貌及搭接物梯度性质的强度预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108256281A true CN108256281A (zh) | 2018-07-06 |
CN108256281B CN108256281B (zh) | 2020-03-10 |
Family
ID=62746378
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810250549.4A Active CN108256281B (zh) | 2018-03-26 | 2018-03-26 | 一种考虑搭接界面形貌及搭接物梯度性质的强度预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108256281B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109101701A (zh) * | 2018-07-20 | 2018-12-28 | 昆明理工大学 | 一种自冲铆接头钉脚张开度的预测方法 |
CN109918838A (zh) * | 2019-03-26 | 2019-06-21 | 孙凌玉 | 一种预测含胶接大型复杂结构服役性能的方法及系统 |
CN114235576A (zh) * | 2021-12-16 | 2022-03-25 | 中南大学 | 一种利用拉伸单剪方法定性分析多层异质梯度材料最弱界面的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7203628B1 (en) * | 1999-08-23 | 2007-04-10 | St Ville James A | Manufacturing system and method |
CN101567025A (zh) * | 2009-05-31 | 2009-10-28 | 湘潭大学 | 用于涡轮叶片热障涂层系统破坏过程的有限元建模方法 |
CN103605865A (zh) * | 2013-11-28 | 2014-02-26 | 北京航空航天大学 | 一种基于内聚力模型的复合材料多向铺层板分层扩展行为模拟方法 |
CN103955604A (zh) * | 2014-04-11 | 2014-07-30 | 南京航空航天大学 | 一种含裂纹金属梯度材料剩余强度预测方法 |
CN105740566A (zh) * | 2016-02-18 | 2016-07-06 | 浙江大学 | 一种预测层状复合材料层内损伤和层间分层的有限元方法 |
-
2018
- 2018-03-26 CN CN201810250549.4A patent/CN108256281B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7203628B1 (en) * | 1999-08-23 | 2007-04-10 | St Ville James A | Manufacturing system and method |
CN101567025A (zh) * | 2009-05-31 | 2009-10-28 | 湘潭大学 | 用于涡轮叶片热障涂层系统破坏过程的有限元建模方法 |
CN103605865A (zh) * | 2013-11-28 | 2014-02-26 | 北京航空航天大学 | 一种基于内聚力模型的复合材料多向铺层板分层扩展行为模拟方法 |
CN103955604A (zh) * | 2014-04-11 | 2014-07-30 | 南京航空航天大学 | 一种含裂纹金属梯度材料剩余强度预测方法 |
CN105740566A (zh) * | 2016-02-18 | 2016-07-06 | 浙江大学 | 一种预测层状复合材料层内损伤和层间分层的有限元方法 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109101701A (zh) * | 2018-07-20 | 2018-12-28 | 昆明理工大学 | 一种自冲铆接头钉脚张开度的预测方法 |
CN109101701B (zh) * | 2018-07-20 | 2022-04-12 | 昆明理工大学 | 一种自冲铆接头钉脚张开度的预测方法 |
CN109918838A (zh) * | 2019-03-26 | 2019-06-21 | 孙凌玉 | 一种预测含胶接大型复杂结构服役性能的方法及系统 |
CN114235576A (zh) * | 2021-12-16 | 2022-03-25 | 中南大学 | 一种利用拉伸单剪方法定性分析多层异质梯度材料最弱界面的方法 |
CN114235576B (zh) * | 2021-12-16 | 2023-11-10 | 中南大学 | 一种利用拉伸单剪方法定性分析多层异质梯度材料最弱界面的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108256281B (zh) | 2020-03-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108256281A (zh) | 一种考虑搭接界面形貌及搭接物梯度性质的强度预测方法 | |
Guo et al. | A layerwise theory for laminated composites in the framework of isogeometric analysis | |
Sadeghi et al. | Failure load prediction of adhesively bonded single lap joints by using various FEM techniques | |
US20090192766A1 (en) | Method for simulating the behavior of a bonded joint of two parts | |
Stein et al. | Analytical models for functionally graded adhesive single lap joints: A comparative study | |
Sarrado et al. | Finite-thickness cohesive elements for modeling thick adhesives | |
CN103345545A (zh) | 一种基于三角包络的复合材料π形非平面胶接连接强度预测方法 | |
Hua et al. | Three-dimensional modeling of carbon/epoxy to titanium single-lap joints with variable adhesive recess length | |
Marzi et al. | A rate-dependent, elasto-plastic cohesive zone mixed-mode model for crash analysis of adhesively bonded joints | |
Mortensen et al. | Coupling effects in adhesive bonded joints | |
Cognard et al. | Influence of the geometry of coaxial adhesive joints on the transmitted load under tensile and compression loads | |
Salve et al. | Implementation of cohesive zone in ABAQUS to investigate fracture problems | |
Liao et al. | Numerical analysis on load-bearing capacity and damage of double scarf adhesive joints subjected to combined loadings of tension and bending | |
Stapleton et al. | Adaptive shape functions and internal mesh adaptation for modeling progressive failure in adhesively bonded joints | |
Kim et al. | Comparative study of failure models for prediction of mixed-mode failure characteristics in composite adhesively bonded joint with brittle/Quai-brittle adhesive using finite element analysis | |
CN110889253B (zh) | 复合材料层合板等效方法 | |
Gleich et al. | Shear and peel stress analysis of an adhesively bonded scarf joint | |
Pickett et al. | The analysis of elastic adhesive stresses in bonded lap joints in FRP structures | |
Das et al. | Three-dimensional nonlinear analyses of scarf repair in composite laminates and sandwich panels | |
Panigrahi et al. | Adhesion failure propagation in adhesively-bonded single-lap laminated FRP composite joints | |
Barros et al. | Numerical simulations of crack propagation tests in adhesive bonded joints | |
Muñoz-Reja et al. | Convergence of the BEM Solution Applied to the CCFFM for LEBIM | |
Behera et al. | 3-D interfacial stress analysis of adhesively bonded curved laminated FRP composite single lap joint | |
Larsson et al. | Geometrically non‐linear damage interface based on regularized strong discontinuity | |
Barros et al. | Use of cohesive zone models to design automotive bonded joints |
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 |