界面内聚力模型及有限元法(20200825195115)
FEM_有限元法 PPT
❖ 应用范围
广泛地被应用于各种结构工程 成功地用来解决其他工程领域中的问题
➢热传导、渗流、流体力学、空气动力学、土壤力学、 机械零件强度分析、电磁工程问题等等
有限元法
❖ Finite Element Method的缩写，有限单元法，其实际应用 中往往被称为有限元分析（FEA），是一个数值方法解偏微 分方程。FEM是一种高效能、常用的计算方法，它将连续体 离散化为若干个有限大小的单元体的集合，以求解连续体力 学问题。有限元法在早期是以变分原理为基础发展起来的， 所以它广泛地应用于以拉普拉斯方程和泊松方程所描述的各 类物理场中（这类场与泛函的极值问题有着紧密的联系）。 自从1969年以来，某些学者在流体力学中应用加权余数法 中的迦辽金法(Galerkin)或最小二乘法等同样获得了有限元 方程，因而有限元法可应用于以任何微分方程所描述的各类 物理场中，而不再要求这类物理场和泛函的极值问题有所联 系．
❖ FEM是应用于现代复杂机械结构优化设计的非常重要的计算 机辅助分析方法。FEM早期主要应用于航空航天制造、船舶 工业及高端军事领域 。
方法运用的基本步骤
❖ 基本思想：由解给定的泊松方程化为求解泛函的极值问题。
❖ 步骤1：剖分 ❖ 将待解区域进行分割，离散成有限个元素的集合．元素பைடு நூலகம்单
元）的形状原则上是任意的．二维问题一般采用三角形单元 或矩形单元，三维空间可采用四面体或多面体等．每个单元 的顶点称为节点（或结点）。
有限元法介绍
通俗地说，有限元法就是一种计算机模拟技术，使人们能够在计算机上用软件模拟一个工程问题的发生过程而无需把东西真的做出来。
这项技术带来的好处就是，在图纸设计阶段就能够让人们在计算机上观察到设计出的产品将来在使用中可能会出现什么问题，不用把样机做出来在实验中检验会出现什么问题，可以有效降低产品开发的成本，缩短产品设计的周期。
有限元法也叫有限单元法（finite element m ethod, FEM），是随着电子计算机的发展而迅速发展起来的一种弹性力学问题的数值求解方法。
五十年代初，它首先应用于连续体力学领域—飞机结构静、动态特性分析中，用以求得结构的变形、应力、固有频率以及振型。
由于这种方法的有效性，有限单元法的应用已从线性问题扩展到非线性问题，分析的对象从弹性材料扩展到塑性、粘弹性、粘塑性和复合材料，从连续体扩展到非连续体。
有限元法最初的思想是把一个大的结构划分为有限个称为单元的小区域，在每一个小区域里，假定结构的变形和应力都是简单的，小区域内的变形和应力都容易通过计算机求解出来，进而可以获得整个结构的变形和应力。
事实上，当划分的区域足够小，每个区域内的变形和应力总是趋于简单，计算的结果也就越接近真实情况。
理论上可以证明，当单元数目足够多时，有限单元解将收敛于问题的精确解，但是计算量相应增大。
为此，实际工作中总是要在计算量和计算精度之间找到一个平衡点。
有限元法中的相邻的小区域通过边界上的结点联接起来，可以用一个简单的插值函数描述每个小区域内的变形和应力，求解过程只需要计算出结点处的应力或者变形，非结点处的应力或者变形是通过函数插值获得的，换句话说，有限元法并不求解区域内任意一点的变形或者应力。
大多数有限元程序都是以结点位移作为基本变量，求出结点位移后再计算单元内的应力，这种方法称为位移法。
有限元法本质上是一种微分方程的数值求解方法，认识到这一点以后，从70年代开始，有限元法的应用领域逐渐从固体力学领域扩展到其它需要求解微分方程的领域，如流体力学、传热学、电磁学、声学等。
两相介质动力反应计算分析的时域显式有限元方法
目录
• 引言 • 两相介质动力反应理论基础 • 时域显式有限元方法概述 • 两相介质动力反应计算的时域显式有
限元方法 • 方法改进与优化 • 结论与展望
01
引言
研究背景与意义
两相介质动力反应在石油、化工、能 源等领域具有广泛应用，准确预测其 反应过程对于提高生产效率和安全性 具有重要意义。
通过与实验数据的对比，验证了该方法的有效性和准确性，为两相介质动力反应的数值模拟提供了新的 工具。
该方法在石油、化工、水力等领域具有广泛的应用前景，能够为多相流反应器的设计、优化和安全运行 提供理论支持。
未来研究方向
结合高性能计算技术，实现大规模两相介质动 力反应的实时模拟，为工程实际提供更为准确
两相介质动力反应方程
方程形式
采用拉格朗日坐标或欧拉坐标描述， 涉及到相间作用力、相对速度、界面 条件等复杂因素。
求解方法
采用数值方法进行求解，如有限元法 、有限差分法等。
初始条件和边界条件
初始条件
描述介质在初始时刻的状态，如初始压力、速度、温度等。
边界条件
描述介质在边界上的行为，如给定压力、速度、热流量等。
网格更新与调整
根据计算结果和误差分析，动态更新和调整网格， 以适应不同阶段和区域的问题求解需求。
3
网格技术与有限元方法的结合
将自适应网格技术与有限元方法有机结合起来， 形成高效的两相介质动力反应计算分析方法。
06
结论与展望
研究成果总结
提出了一种基于时域显式有限元方法的两相介质动力反应计算分析方法，该方法能够准确模拟两相介 质的动态行为，并考虑了流体与固体之间的耦合作用。
求解离散化方程
基于微观图像及内聚力模型的复合材料裂纹扩展模拟
基于微观图像及内聚力模型的复合材料裂纹扩展模拟秦丽媛;孟松鹤;李金平;金华【摘要】为研究材料微观结构及晶界强度对材料力学性能的影响,在晶界处引入内聚力单元模型,模拟晶间破坏过程.以ZrB2-SiC复合材料为研究对象,将其扫描的微观结构图片进行矢量化处理,并导入ABAQUS有限元软件中建立模型,同时在其晶界处,设置内聚力单元模拟晶界破坏过程.通过改变ZrB2与SiC相界面强度,得到了晶界及材料不均匀对材料应力分布及裂纹扩展的影响.结果表明,由于晶界的存在,材料内部出现应力分布不均匀现象并产生应力集中.随着晶界强度的改变,裂纹起始位置及扩展方向发生改变,且裂纹沿低强度的界面进行扩展.随着ZrB2-SiC界面强度增大,材料的强度提高,拉伸模量不变.%The cohesive zone model is used to simulate the intergranular fracture propagation to get the effects of the microstructure and the interfacial strength on the mechanical properties of the composites.The microstructure morphology of ZrB2-SiC ceramics illustrated by SEM was imported into ABAQUS finite element software,with the cohesive elements embedded along grain boundaries,to predict the intergranular fracture propagation.The effects of grain boundary and material inhomogeneity on the stress distribution and crack propagation were investigated by changing the ZrB2-SiC interfacial strength.The result shows an uneven stress distribution and stress concentration are observed due to the presence of grain boundaries.The locations of the crack initiation and propagation direction vary with the interfacial strengths and the crack propagates along the grain boundaries with low strength.Thestrength of polycrystals increases as ZrB2-SiC interfacial strength increases, and the elastic modulus remains substantially unchanged.【期刊名称】《固体火箭技术》【年(卷),期】2017(040)004【总页数】5页(P501-505)【关键词】微观结构;内聚力单元法;晶界强度;ZrB2-SiC陶瓷【作者】秦丽媛;孟松鹤;李金平;金华【作者单位】哈尔滨工业大学复合材料与结构研究所,哈尔滨 150008;哈尔滨工业大学复合材料与结构研究所,哈尔滨 150008;哈尔滨工业大学复合材料与结构研究所,哈尔滨 150008;哈尔滨工业大学复合材料与结构研究所,哈尔滨 150008【正文语种】中文【中图分类】TB332复合材料是两种或两种以上的组分组成的多相材料，且各相之间具有明显的界面。
基于内聚力理论的二维二次界面单元在ABAQUS中的UEL程序实现
基于内聚力理论的二维二次界面单元在ABAQUS中的UEL程序实现刘敏; 李旭【期刊名称】《《计算力学学报》》【年(卷),期】2019(036)005【总页数】6页(P693-698)【关键词】内聚力模型; 界面单元; 有限元法; ABAQUS; UEL【作者】刘敏; 李旭【作者单位】武汉理工大学理学院工程结构与力学系武汉430070; 武汉理工大学理学院新材料力学理论与应用湖北省重点实验室武汉430070【正文语种】中文【中图分类】O341; TB1241 引言内聚力模型的概念最早由Dugdale等[1,2]提出，该模型假设在实际裂纹尖端的前部存在一个内聚力区,如图1所示，该区域由两个相邻(不受力时无厚度)的上下表面组成，通过界面粘结力(内聚力)连接在一起。
当存在外力作用时，上下表面将逐渐发生分离，造成界面开裂或裂纹扩展，且此过程中的界面粘结力大小由界面间的相对分离位移来确定。
内聚力模型立足于相邻界面间的粘结应力-分离位移关系，也称内聚力法则，经历了数十年的研究发展，实践证明该模型能够较有效地表征粘结界面受载荷作用时的损伤失效过程[3]。
将有限元法与内聚力模型相结合，可定义得到内聚力界面有限单元[4，5]，进而较方便地实现各种材料和结构内部界面破坏或裂纹扩展过程的数值分析/求解[6-8]。
在目前主流的有限元软件(ABAQUS，ANSYS及MARC等)中，都能够使用内聚力界面单元进行相关的模拟计算。
当然，现阶段通用有限元软件在这方面的功能还比较有限，以ABAQUS为例，其单元库中仅有位移函数线性插值的一次界面单元，内聚力模型也只提供双线性和指数型两种法则[9]。
然而，当考察的问题几何形状比较复杂时，如颗粒随机分布的橡胶高弹体复合材料，其内部结构呈现出高度非均匀性，一般需采用三角形或四面体单元建模，此时线性单元的精度往往有所不足，无论是表征橡胶材料自身或是界面层的力学性能，都最好选择二次单元[10,11]，这种情况下商用有限元软件便无法满足相应的分析需求。
有限元法概述
（2）MSC/NASTRAN。 MSC/NASTRAN是在原NAST RAN基础上进行大量改进后的系统软件，主要包括MS C.Patran并行框架式有限元前后处理及分析系统、 MS C.GS-Mesher快速有限元网格、 MSC.MARC非线性有 限元软件等。其中MSC.MARC具有较强的结构分析能
.
5.在产品制造或工程施工前预先发现潜在的问题； 6. 模拟各种试验方案，减少试验时间和经费； 7. 进行机械事故分析，查找事故原因。
轴承强度分析
.
汽车碰撞实验
.
刹车制动时地盘的应力分析
.
钢板精轧机热轧制分析
.
三维椭圆封头开孔补强
.
水轮机叶轮的受力分析模拟
.
人体股骨端受力分析
.
半导体芯片温度场的数值仿真
知量时称为混合法。 位移法易于实现计算自动化，所以，在有限单元法
中位移法应用范围最广。
.
2、有限元法的发展
有限单元法基本思想的提出，可以追溯到Courantl在1 943年的工作，他第一次尝试应用定义在三角形区域上的 分片连续函数和最小位能原理相结合，来求解St·Venant 扭转问题。相继一些应用数学家、物理学家和工程师由于 各种原因都涉足过有限单元的概念。
.
4、有限元的特点
(1) 概念清楚,容易理解。可以在不同的专业背景和水平 上建立起对该方法的理解。从使用的观点来讲,每个人的 理论基础不同,理解的深度也可以不同,既可以通过直观的 物理意义来学习,也可以从严格的力学概念和数学概念推 导。
沥青混凝土的细观开裂模拟方法研究
本研究的目的是建立一个更有效、更精确的细观沥青混凝土模型研究其开裂 问题。本文采用两种数值建模方法分析沥青混凝土的开裂问题。第一种是随机骨 料分布及嵌入粘结单元有限元模型。本方法将沥青混凝土分为四个重要组成部分： 骨料、沥青砂浆、界面过渡区和初始缺陷。根据由 Full 和 Thompson 开发的级配 曲线，将骨料颗粒模拟成大小不同的随机分布多边形。采用了一些有效的方法来 提高骨料投递的成功率。对于初始缺陷，本研究仅考虑空隙。在沥青砂浆的初始 网格内，沿骨料与沥青砂浆的界面插入零厚粘性单元，模拟混凝土的开裂过程。 修正了以往的开裂初始准则和牵引分离规律，并用以描述粘结构件的破坏行为。 基于 python 语言，开发了骨料投递、内聚力单元插入和修正的构造本构模型的 用户自定义程序，并将其嵌入到商业有限元软件包 abaqus 中。通过与试验结果 的比较，验证了所提出的细观模型的有效性和准确性，并研究了细观结构对沥青 混凝土的宏观性能的影响。第二种方法采用基于 Voronoi 多边形的细观刚体弹簧 法进行数值建模，在第一种方法的基础上，采用随机骨料投递技术随机投放圆形 骨料，然后根据骨料形心点生成 Voronoi 网格，过 Voronoi 单元网格边界的几何 形状寻找骨料之间的相互作用关系，骨料间的粘结材料被浓缩为骨料间的界面弹 簧，其刚度由粘结材料的厚度定义。通过二维和三维数值单轴压缩实验进行验证 其有效性。 关键词：细观沥青混凝土；随机骨料投递算法；粘结单元嵌入算法；Voronoi 多 边形；细观刚体弹簧法
i
ABSTRACT
The purpose of this study is to establish a more effective and accurate mesoscopic asphalt mixture model to study its cracking problem.In this paper, two numerical modeling methods are used to analyze the cracking of asphalt mixture.The first is the random aggregate distribution and the embedded bond element finite element model.This method divides asphalt mixture into four important components: aggregate, asphalt mortar, interface transition zone and initial defect.Aggregate particles were simulated as randomly distributed polygons of different sizes according to the grading curves developed by Full and Thompson.Some effective methods are adopted to improve the success rate of aggregate delivery. For the initial defects, only the void was considered in this study. In the initial grid of asphalt mortar, zero-thickness viscous unit is inserted along the interface between aggregate and asphalt mortar to simulate the cracking process of concrete.The former cracking initial criterion and traction separation rule are modified to describe the failure behavior of bonded components. Based on the python language, a user-defined program for aggregate delivery, cohesive force unit insertion and modification of constitutive model construction was developed and embedded in abaqus, a commercial finite element software package. Compared with the experimental results, the validity and accuracy of the proposed meso-structure model are verified, and the effect of meso-structure on the macro performance of asphalt mixture is studied. The second method based on Voronoi polygon of the mesoscopic numerical model based on the rigid spring method and random on the circular aggregate by random aggregate delivering technology, then Voronoi grids were generated according to the aggregate centroid and a Voronoi unit grid boundary geometry for the interaction relationship between aggregate and bond between aggregate material is concentrated to aggregate the interface between the spring, defined by the thickness of the bonding material, The validity of the method is verified by two-dimensional and three-dimensional numerical uniaxial compression experiments.
内聚力模型在复合材料界面层中应用
内聚力模型在复合材料界面层中应用加工制造复合材料时，界面最难控制且大多是最为薄弱的部分，在承受载荷时，界面层处也往往最先出现破坏，从而大大降低了复合材料的强度。
因此在进行复合材料细观模型分析时，界面层不能够被忽略。
界面的强度对整个单胞模型强度的影响非常大，本文通过纤维顶出法计算出复合材料界面层内聚力模型的本构参数，然后分别对代表体积元RVE模型在纵向拉伸载荷和横向拉伸载荷下做有限元计算，分析界面损伤和基体塑性变形的数值关系，结果表明，复合材料的纵向拉伸强度随界面强度增加而增加，但到一定程度后影响就变小了，代替的是基体拉伸强度起主要作用。
不论界面强度值增加到多大，所得到的复合材料的横向拉伸强度都永远低于基体的拉伸强度980Mpa，这充分的说明了复合材料的强度极限是由最最薄弱的部分决定的。
本文的研究结论可以为复合材料的制备提供理论参照。
标签：内聚力模型;界面层;失效模式0 引言在对TiC/TC4复合材料应用内聚力模型进行细观力学模型计算之前，必须先确定界面层的cohesive本构关系。
内聚力本构模型的参数是由实验测得的，但大多数情况下，一些材料的基本力学性能参数是没有的，此时想要继续研究，就需要做一些理论推导与假设。
本文就是由纤维顶出法测得界面结合强度与界面反应厚度的关系式，然后经过一系列计算，就可以确定界面层内聚力模型本构的基本参数，进而展开后续的分析计算研究。
1 理论分析1.1 参数的计算采用纤维顶出法试验测得SiC/TC4复合材料的界面结合强度可以用下式表示：，其中是纤维与基体之间界面的剪切强度值，H是界面层的厚度。
经有限元计算分析，复合材料的界面结合强度与载荷大小的关系为，其中是加载时的最大载荷。
其中K是内聚力本构模型的刚度值，h是薄片试样的厚度，E是弹性模量，是应力，是应变，是位移。
这里假设界面的结合强度等于本构关系中的最大应力取H=1um，E=330Gpa，=6.2N，代入计算，可以求最大应力为102.5Mpa，试件的厚度为0.21mm，位移为0.065um，本构关系=0时对应的位移取0.2um，计算得到断裂韧性为10.25。