首页/文章/ 详情

虚单元法(Virtual Element Method)前沿研究分享

10月前浏览1034

文一:

 

Timoshenko 梁的一种无锁定虚元方法

摘要:

从Beirão da Veiga等人(2013)对弹性固体的贡献开始,虚拟单元法(VEM)为结构和固体力学问题的离散化提供了新的方法。有趣的是,虚拟单元法还允许重新审视与有限元具有相同形状的不同单元的构造。这甚至适用于桁架和梁等一维结构。在这里,我们研究了Timoshenko梁的虚拟单元发展,这令人惊讶地导致了锁定自由的、在线性范围内甚至精确的Timoshenk梁单元的直接公式。这些单元可以很容易地合并到经典的有限元代码中,因为它们与梁的有限元具有相同数量的未知数。该公式允许使用Reissner(1972)中提供的公式来计算经历大挠度和旋转的非线性结构问题。

 

图:梁在恒定线荷载作用下与相关单元自由度的关系。

 

图:线性变化线荷载下的梁。

 

图:(1st/2nd和(2nd/3rd)插值的两种不同虚元公式的收敛行为。

 

图:(a) 具有垂直荷载的框架和(b)离散化。

 

图:负载偏转,包括咬合。

文二:

 

有限应变应用中的无稳定虚拟单元法

摘要:

本文针对二维可压缩超弹性材料,提出了一种新的无稳定度较高的虚拟单元方法。与最传统的虚拟单元公式不同,该方法不需要任何稳定性。主要思想是修改虚拟单元空间,以允许计算高阶多项式𝐿2梯度的投影。在此基础上,可以直接得到刚度矩阵,极大地简化了分析过程,尤其是对于非线性问题。考虑了超弹性材料,解决了一些基准非线性问题,验证了无稳定虚拟单元法的性能和准确性。

 

图:几何非线性分析中的参考配置和当前配置。

 

图:一阶和二阶虚拟单元。

 

图:多边形单元的三角剖分。

 

图:具有不同网格的库克膜问题的初始配置。

 

图:具有非凸单元的不同网格垂直位移的变形形状和轮廓曲线。

 

图:具有不同大小的非结构化多边形网格。

 

图:通过不同网格的无稳定VEM获得的数值解。

文三:

 

不可压缩和不可扩展问题的混合虚拟单元公式

摘要:

在弹性材料的数值模拟过程中,锁定效应可能是一个主要问题,尤其是对于大应变。这些影响源于体积约束,例如底层材料类别的不可压缩性或各向异性影响。一种特定的解决方案策略是采用混合配方,该配方提供针对手头特定锁定现象量身定制的解决方案。虽然这种解决方案策略通常是强大的,但其一个缺点是,所提供的克服锁定的策略往往局限于某些特定拓扑类型的有限元,从而丧失了通用性。相反,虚拟单元法(VEM)从定义上受益于在元素级别允许任意元素形状和节点数量。在这篇文章中,我们为低阶模拟的虚拟单元法制定了多种处理超弹性材料锁定现象的方法。在VEM中实现多场混合原理的关键因素是每个场只考虑一个常量变量和整个虚拟单元上的一个相应拉格朗日乘子。因此,稳定贡献利用混合公式,但与虚拟单元的投影部分共享单元常量变量。这种相当简单的实现策略的直接结果是将强大的混合公式与能够处理一般单元形状的计算方法相结合。所提出的公式在计算力学的标准示例以及也使用非结构化网格的特定计算工程应用中关于结构化网格进行了测试。

 

图:二维冲压问题:原理图设置;Mooney-Rivlin型高度不可压缩各向同性材料行为的b VEM Q1 CoJP预测。

 

图:二维和三维库克膜问题的示意图设置

 

图:具有横向各向异性Neo-Hookean材料定律的二维Cook膜问题的变形构形。

 

图:具有横向各向异性Mooney-Rivlin材料定律的三维Cook膜问题的变形构形。

 

图:具有100个晶粒的结晶微观结构的网格。

 

图:负载下的二维双夹紧补片问题。VEM Q1 TrJP预测的结果。

 

图:多尺度方案的示意图:微观计算问题的局部化和均匀化。

文四:

 

计算各向异性晶体塑性的虚拟单元

摘要:

在这篇文章中,将具有线性模拟的虚拟单元法(VEM)应用于微结构环境中的计算晶体塑性框架。此外,基于张量变形测度和结构张量的不变公式,对底层晶体结构的立方弹性各向异性提出了一个简单的各向异性能量贡献。各向异性弹性公式在小变形的极限下恢复了立方体材料的弹性张量结构。作者提出了一个新的稳定退化公式,该公式完全基于问题的耗散响应。代表性的例子说明了VEM在晶体塑性框架中锁定现象方面的稳健性和性能,当与经典有限元方法的解进行对比时。进一步的例子研究了VEM在晶体塑性框架内的性能和电流限制,当应用于结构化元件拓扑和柔性元件拓扑的异质微结构时。

 

图:示意图。

 

图:由均匀化杨氏模量的Reuss估计加权的定向杨氏模量的图形表示。

 

图:三维库克膜问题的几何和加载设置图解。

 

图:(a)、(b)的等值线图:纯弹性响应的等效von Mises应力分布;(c) ,(d):弹塑性响应的累积塑性剪切分布。

 

图:(a)-(c)旋转角度预测图解;(d) –(f)旋转角绝对值的差异,由VEM和FEM在正方形2D域的纯剪切变形后预测,由四边形单元离散。

 

图:具有300个晶粒的有限条域的图解,通过(a)VEM-VO方法进行网格划分;(b) FEM-O2方法;(c) 单个晶粒的方向颜色图例。单个颗粒的颜色是针对旋转角度获得的三种颜色的混合。

 

图:悬臂梁域图示为(a)示意图;(b) 粗网格,将单个六面体元素视为单个晶粒;(c) 精细网格域。

 

图:含软夹杂物的晶体试样及其网格。

文五:

 

非协调网格三维接触问题的虚拟单元法

摘要:

虚拟单元法非常适合于不一致网格的问题,因为可以在虚拟单元边界处插入节点以获得一致的网格结构,而无需改变单元的模拟。接触条件可以用于不同的执行策略。这种特征也可以用于三维接触问题。在这项工作中,我们介绍了一个新的框架,包括重新网格策略,它可以产生一致的网格。将显示3D中的各种数值示例,包括接触贴片测试和赫兹接触等标准情况。

 

图:联系边值问题。

 

图:参数空间中的三角形(左)和三维欧氏空间(右)。

 

图:节点到曲面的离散化。

 

图:接触片在接触界面测试网格的边值问题。

 

图:节点到节点和节点到曲面接触的修补程序测试。

 

图:在接触界面上与不匹配的三角形和四面体进行节点到节点的接触。

来源:STEM与计算机方法
非线性通用理论自动驾驶材料多尺度数字孪生人工智能曲面
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-11-01
最近编辑:10月前
江野
博士 等春风得意,等时间嘉许。
获赞 58粉丝 96文章 153课程 0
点赞
收藏
作者推荐

高分子材料(Polymer materials)的数值模拟前沿研究分享

文一: 聚合物复合材料断裂的梯度相间增强相场模拟方法摘要:与微米级夹杂物相比,纳米级填料颗粒在改善聚合物的弹性和断裂行为方面的优势已通过几项实验观察得到证实。然而,对这种非均质材料中裂纹扩展的精确建模涉及复杂裂纹拓扑结构的解析,同时考虑多个裂纹的聚结和分支。这种复杂性使得传统的尖锐裂纹建模方法,例如基于单位分割思想的方法,适用性有限,尤其是在3D中,因为这些方法涉及对演化裂纹表面的显式跟踪。因此,相场断裂方法已经成为传统尖锐裂纹模型的一种有吸引力的替代方法,尤其是当需要处理复杂的裂纹拓扑结构时。此外,标准的两阶段连续谱模型由于缺乏必要的长度尺度,无法捕捉到前面提到的较小即较强的尺寸效应。为了弥补这一缺陷,在有限应变超弹性的背景下,我们在先前的工作中引入了标准一阶连续体模型的基于梯度界面的增强。本贡献将梯度界面的思想扩展到相场断裂方法。在此,考虑了如实验所观察到的具有连续变化或分级材料性质的填料颗粒周围的界面区域。复合材料的断裂行为可以通过用于确定界面区域内的弹性和断裂特性的分级程度来控制。与标准相场断裂模型相比,分级相间参数的最佳组合产生了更硬的宏观响应,如在更硬相间的情况下所观察到的。通过广泛的数值实验,描述了所引入的技术对广泛的实验观察到的断裂行为进行建模的适当性,这取决于填料和基体相之间的粘附程度。 图:连续体的尖端裂纹示意图和弥散相场裂缝模型的示意图。 图:梯度相界概念。 图:建模示意图。 图:标准算例的断裂相场建模示意图。 图:在低相间断裂能的情况下,分级相间弹性对复合材料断裂响应的影响。文二: 模拟聚合物复合材料尺寸效应的增强计算均匀化技术摘要:几项实验研究证实,纳米夹杂物是比微米夹杂物更有效的增强聚合物的增强材料。标准的一阶计算均匀化方案由于缺乏必要的长度尺度来模拟这种尺度效应而存在不足,需要对标准方案进行改进。本文基于界面能量学的思想,对一个这样的扩展进行了深入的评估。系统的数值试验和分析表明了上述方法在复合材料力学行为建模中的局限性,即填充材料比基体更加刚性。提出了一种基于连续分级中间相思想的方法。通过具有代表性的数值算例对该方法进行了综合评价,表明该方法适用于不同填料-基体刚度组合的纳米复合材料的建模。 图:计算均匀化过程示意图,描绘了两个相关尺度下的BVP。 图:RVE 由基体材料(蓝色显示)和圆形夹杂物(红色显示)组成,具有(a)含能界面或(b)厚度为 t 的梯度界面(绿色显示)。 图:二氧化硅填料颗粒在环氧树脂基体中的相间性质随指数 n 值的变化。 图:IECH尺寸效应缩放:环氧树脂基体中橡胶颗粒的宏观Piola应力随填料颗粒尺寸的变化。 图:在单轴拉伸载荷作用下,橡胶颗粒在环氧树脂基体中的粒径和界面强度发生变化时,相应的柯西应力分量 σ11的节点值的云图。 图:在选定的颗粒半径下,RVE 在剪切载荷作用下的变形形状。为了显示填料和基体的相对位移,将变形形状叠加在未变形形状上。 图:相间厚度对容积加载尺寸效应标度的影响。 图:三维RVE的结构化六边形网格,a=100,r=25(和t=5),分别用于IECH和GICH。 图:三维强度尺度参数对两种方法的影响。本文研究了在体积载荷作用下,环氧树脂基体中二氧化硅颗粒在 a = 100,r = 25(和 t = 5)的固定微结构构型。文三: 扩展一个通用的快速粗粒度分子动力学模型来研究接枝聚合物纳米复合材料的力学行为摘要:聚合物纳米复合材料因其高通用性和良好的力学性能以及低密度而成为工程应用的一类重要材料。通过将聚合物链直接连接到纳米填料上,即所谓的接枝,实现了基质和填料之间更好的负载转移,此外,还获得了填料更好的分散性。两者都能增强力学性能。由于纳米尺度的实验研究极具挑战性,因此需要补充的数值研究来揭示聚合物纳米复合材料的力学行为。为此,分子动力学非常适合,因为它可以捕捉微观结构,但在数值上计算成本也很高。因此,这一贡献为研究接枝聚合物纳米复合材料的力学行为提供了一个快速粗粒分子动力学模型。为此,我们通过接枝键扩展了现有的模型,这使我们能够直接比较未处理和接枝填料的效果。特别是,我们研究了填料含量、接枝度和填料尺寸对聚合物(接枝)纳米复合材料刚度和强度的影响。我们得出的结论是,与未处理的填料相比,接枝键对刚度的影响很小,而强度显著提高,这与文献一致。所提出的聚合物接枝纳米复合材料的分子动力学模型为进一步研究,特别是关键的基质-填料界面提供了基础。此外,这一贡献将分子动力学的见解转化为力学性能,从而弥合了与工程规模的差距,从而代表着朝着开发聚合物(接枝)纳米复合材料的全部潜力迈出了一步。 图:嵌入热塑性基体的无机填料: a)未经处理填料的聚合物纳米复合材料(PNC) ; b)接枝填料的聚合物接枝纳米复合材料(PGN) ,即聚合物链与填料表面之间的共价键。 图:单轴拉伸模拟的刚度和强度。 图:填料含量的影响。 图:接枝键在聚合物链中的位置。文四: 非线性光弹性的数值模拟摘要:当分子光开关,如偶氮苯或降冰片二烯,嵌入到足够柔软的聚合物基体中时,所得化合物可以经历由特定波长的光引起的机械变形。这些光敏化合物具有作为软激励器应用的潜力,而不需要硬接线电子器件或单独的能源。这种特性在微型机器人的设计中 特别有吸引力,但其他应用,如高速数据传输或将光子能量转换为机械响应,也有很大的前景。尽管有这些几乎未来主义的可能性,光敏聚合物在工业应用中还没有得到足够的关注。提高这类软智能材料接受度的一个重要因素是结合数值模拟方法制定严格的本构建模方法。因此,在这篇文章中,我们提出了一种光机械建模方法,与之前发表的基本原理不同。我们简要介绍了必要的本构方程,这些方程随后与各自的平衡定律结合使用到有限元实现中。最后,通过一个简单的二维基准点示例说明了数值求解方法的能力,并随后将其扩展到更复杂的三维问题。 图:光致变形下未变形体和变形体的示意图。 图:(左)在可变强度光源照明下的材质样本草图。(右)施加10 mV电势差时有限元模型的变形。颜色映射和箭头指的是均匀化电子序参数的大小和方向。 图:(左) X 方向(虚线) ,Y 方向(虚线)和不同穿透深度的大小(实线)材料样品顶边中心点的位移。(右)变形材料样本。颜色是指在左边的图中选择的颜色。 图:(左)花瓣的一半与140个六面体单元相啮合的模型。(中)闭合的花,由八个完整的花瓣组成,以2 V的施加电势排列成花朵。(右)开放的花,以0.2 V的施加电位开放。彩色图是指位移大小。 图:(左)对于各种材质厚度值,花瓣尖端在Y方向(虚线)、Z方向(短划线)和大小(实线)上的位移。(右)变形的花瓣。着色是指变形。文五: 柔性管状金属-聚合物粘接接头的数值研究摘要:粘结管状接头已在各种工程和工业应用中得到应用,为连接管状基材提供了一种高效耐用的方法,同时确保了良好的密封性能。在这项工作中,提出了一种新的设计方法,用于在复杂载荷和被粘物大变形下承受大扭转旋转的粘结柔性管状金属-聚合物接头。为此,开发了两个不同的有限元模型。一种是针对较小的旋转值优化接头性能,另一种是使用内聚区建模并引入聚合物管状粘附物中的屈曲效应来预测整体接头强度。为了评估超弹性特性对接头性能和强度的影响,还考虑了不同的聚合物材料。三种不同接头布局的实验结果与数值结果吻合良好。结果表明,对于较小的旋转值,可以使用二维模型来准确预测关节行为。尽管如此,对于聚合物材料的高旋转值和大变形,复杂的3D模型(包括屈曲的影响)对于预测接头强度是必要的。 图:纯I型和II型的线性软化牵引分离定律。 图:管接头的布局、尺寸和胶粘剂厚度(以毫米为单位) : (a)-布局 X,(b)-布局 Y 和(c)-布局 Z。 图:用于制造管状接头的模具。 图:扭转试验装置。 图:不同关节布局的扭矩-扭转曲线实验研究。 图:通过解析解和本模型得到粘结层中的剪应力分布,采用粗网格和细网格。 图:变形形状,最大旋转角度沿轴线 = 180 ° 。 图:三维模型用于预测接头的最终旋转,网格表示。 图:对于布局 X、布局 Y 和布局 Z,前三个折叠模式形状。 图:不同布局的数值和实验扭矩-扭转曲线。 图:在旋转800 ° 时 X 排样的变形形状和粘接层中的 SDEG 损伤变量。来源:STEM与计算机方法

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈