[发明专利]含界面结构复合材料界面相的力学性能三维仿真评估方法在审
| 申请号: | 202111290207.3 | 申请日: | 2021-11-02 |
| 公开(公告)号: | CN114329658A | 公开(公告)日: | 2022-04-12 |
| 发明(设计)人: | 关康;龚振远;饶平根;曾庆丰;刘建涛 | 申请(专利权)人: | 华南理工大学 |
| 主分类号: | G06F30/10 | 分类号: | G06F30/10;G06F30/20;G06T17/00;G06F113/26;G06F119/14 |
| 代理公司: | 深圳紫晴专利代理事务所(普通合伙) 44646 | 代理人: | 郭清秀 |
| 地址: | 510630 广*** | 国省代码: | 广东;44 |
| 权利要求书: | 查看更多 | 说明书: | 查看更多 |
| 摘要: | |||
| 搜索关键词: | 界面 结构 复合材料 力学性能 三维 仿真 评估 方法 | ||
1.含界面结构复合材料界面相的力学性能三维仿真评估方法,其特征在于:具体检测评估如下:
S1:建立具有三维几何结构界面相复合材料代表性体积元模型;
S2:为模型布置全局网格种子点,划分网格,建立覆盖全模型零厚度内聚力单元;
通过软件Abaqus的全局种子点划分功能,为模型布置合适的全局种子点,并划分网格;为实现任意裂纹扩展,基于Abaqus二次开发接口,将已划分好网格的模型生成覆盖全模型的零厚度内聚力单元;内聚力模型是以应力-位移曲线的概念提出的,常见的内聚力单元模型的应力-位移曲线为双线性模型;裂纹尖端前部介于裂纹张开区域和未发生损伤的线弹性固体区域之间,该区域的材料因局部损伤和断裂而具有非线性特点;在内聚力单元损伤起始和演化前,单元服从线弹性关系;弹性关系用界面处的名义应力与名义应变的弹性本构矩阵的形式表示;
S3:输入材料力学参数并施加载荷,赋予正确的网格类型,并采用准静态分析进行模拟拉伸测试;
利用Abaqus软件进行仿真分析,输入合理的复合材料力学参数,施加载荷;为内聚力单元与实体单元赋予正确的网格类型;模拟三维模型拉伸断裂测试,采用准静态分析,使用显式分析步;显式动力学的基本运动方程定义为:
其中M时质量矩阵,是结构的加速度,I是结构中的内力,P为所施加的外力;该公式表述的是牛顿第二运动定律;在显式动力学中,I不仅包括了结构的内力,也包括了源于运动如阻尼和结构变形的贡献;属于无阻尼系统,即I=Ku,因此式写为:
其中K为材料的刚度矩阵,u为节点位移矢量;
S4:界面相力学性能评估:采用SDEG值评估界面相在复合材料中的力学行为影响;在ABAQUS后处理模块中,获得界面内聚力单元的SDEG云图,通过Abaqus软件自带的数据输出功能,输出界面相内聚力单元的SDEG值,并作出平均化处理,用来表征界面相的损伤情况,定量分析界面相对裂纹偏转的影响。
2.根据权利要求1所述的含界面结构复合材料界面相的力学性能三维仿真评估方法,其特征在于:所述步骤S1中建立含有不同几何特征的界面相结构的三维复合材料代表性体积元几何模型,从外侧到最内侧依次是:基体、中间相和纤维。
3.根据权利要求1所述的含界面结构复合材料界面相的力学性能三维仿真评估方法,其特征在于:所述步骤S2在内聚力单元中,名义应力矢量t由三个分量组成:tn,ts,tt;相对应的位移量记为δn,δs,δt;内聚力单元的原始厚度记为T0,名义应变定义为:
线弹性本构方程矩阵形式写成:
其中,tn为法向应力,ts、tt为两个切向应力,δn为法向应变,δs、δt为两个切向应变,Kij为刚度矩阵。
4.根据权利要求3所述的含界面结构复合材料界面相的力学性能三维仿真评估方法,其特征在于:所述步骤S2中的内聚力单元当应力或者应变满足损伤起始准则的时候,损伤起始发生,采用二次名义应力准则,假设损伤起始于名义应力的二次函数达到1,其定义为:
其中,分别代表当变形是纯法向、纯剪切1方向和纯剪切2方向时名义应力的峰值。
5.根据权利要求4所述的含界面结构复合材料界面相的力学性能三维仿真评估方法,其特征在于:所述步骤S2中当内聚力单元的应力状态满足损伤起始准则后,单元就将发生损伤,然后进入损伤演化阶段;采用线性衰减演化;针对内聚力单元的演化,损伤变量D,将用来表示材料中整体损伤和各种机制的联合作用;D的初始值为0,当内聚力单元发生损伤,进入损伤演化阶段,损伤变量D随着载荷的加载从0到1单调增加;当D达到1时,单元完全失效,失去承载能力;于是损伤演化过程中,应力表达式为:
其中,对于法向方向压缩载荷的情况,不考虑损伤演化过程,即此时D=0。
6.根据权利要求1所述的含界面结构复合材料界面相的力学性能三维仿真评估方法,其特征在于:所述步骤S3中能量平衡方程有助于评估准静态分析的合理性,在Abaqus/Explicit中,能量平衡写为:
E1+EV+EFD+EKE-EW=ETOT=cons tan t
其中E1为内能;EV为粘性机制耗散的能量;EFD为摩擦耗散的能量;EKE为动能;EW为外力功;ETOT为系统的总能量;
为确保获得合理的准静态分析过程,EKE只能占据E1的1%-5%。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于华南理工大学,未经华南理工大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202111290207.3/1.html,转载请声明来源钻瓜专利网。
- 上一篇:小区切换方法、装置及计算机可读存储介质
- 下一篇:一种立式气化器





