首页/文章/ 详情

随机有限元分析的前沿研究分享

9月前浏览584

文一:

 

基于Karhunen–Loève展开和修正摄动随机有限元法的结构静力随机分析新计算方案

摘要:

由于不确定性,确定性分析不能充分反映结构的性能。随机分析可以考虑多种不确定性因素的影响,提高分析结果的可信度。针对具有低水平不确定性的结构,提出了一种新的随机计算格式,该格式具有Karhunen–Loève(K–L)展开和修正微扰随机有限元法(MPSFEM)的特点,简称KL-MPSM。材料参数被视为随机场,并通过K–L展开进行离散。将获得的随机变量代入MPSFEM,以获得结构响应的前两阶矩(均值和方差)的估计值。引入JC方法,利用二阶估计计算结构的可靠性指标和失效概率。以深梁和平面框架结构为算例,验证了KL-MPSM的可行性,并研究了一些随机场的性质。结果表明,KL-MPSM具有良好的精度、效率和编程优势。因此,KL-MPSM非常适合于具有低水平不确定性的结构的静态随机分析。

 

图:非正态分布随机变量的等价归一化。

 

图:悬臂深梁的几何形状。

 

图:前五阶本征值与本征函数的叠加。

 

图:MCM和KL-MPSM在不同变异系数下计算的A点垂直位移的平均值和标准偏差。

 

图:平面多层框架的几何形状。

文二:

 

考虑土壤参数空间变异的土石坝随机有限元分析

摘要:

土壤参数的空间变异性显著影响土石坝数值分析的准确性。本研究建立了一种考虑土壤参数空间变异性的随机有限元方法。根据现场监测值确定土壤干密度的自相关距离和分布函数。利用三轴试验结果,建立了干密度与邓肯E-B模型参数之间的回归函数。根据干密度分布函数和这些回归函数,确定了E-B模型参数的累积分布函数,并采用异概率变换方法建立了土壤参数随机场。以观音岩大坝为例,讨论了随机有限元的收敛性、随机场单元大小的影响以及E-B模型参数的敏感性。观音岩组合坝随机有限元分析表明,考虑土体参数空间变异性的最大沉降平均值与反演参数计算结果接近。当随机有限元分析计算超过400倍时,最大沉降的平均值开始在该区域收敛。土壤参数的随机场元大小对随机有限元计算的影响小于对有限元计算。Kb对随机有限元计算的影响最大,Rf的影响最小。

 

图:最佳拟合曲线,邓肯E-B模型参数与堆石料干密度之间的最佳拟合曲线±一个和两个标准差(SD)。

 

图:考虑土石坝土壤参数空间变异性的随机有限元分析流程图。

 

图:观音岩大坝土壤干密度随机场样。

 

图:波动堆石料干密度的尺度: (a)垂直方向; (b)水平方向。

 

图:邓肯E-B模型参数与观音岩复合坝堆石干密度的回归函数。

 

图:使用(a)检查参数计算的沉降量;(b) 反演参数;和(c)一组随机田间土壤参数。

 

图:使用(a)检查参数计算的应力;(b) 反演参数;和(c)一组随机田间土壤参数。

 

图:Kb和Rf随机场的最大沉降分布。

文三:

 

大型混凝土结构在热力-水力-机械载荷作用下的可用性随机有限元分析——以核安全壳建筑为例

摘要:

这项工作提出了一种全局随机有限元方法(SFEM)来模拟混凝土老化不确定性对大型钢筋和预应力结构的可用性和耐久性的影响。由于它们的建模需要具有大量参数的强非线性、耦合和昂贵的计算,因此需要定义适应和有效的概率策略,以在合理成本和物理上可容许的代表性内进行随机分析。在这一贡献中,这是通过四个步骤实现的:(a)基于交错的热电-机械+泄漏(THM-L)模型定义一个公认的物理框架;(b) 使用基于方差的分层和局部灵敏度分析(HLSA)将不确定性传播的随机输入限制为最具影响力的输入;(c) 使用多项式混沌展开(PCE)构建THM-L响应元模型;(d) 将粗蒙特卡罗方法(CMCM)应用于所开发的元模型,对可用性准则进行可靠性分析。为了验证和证明在复杂的工业框架内的可实现性,将这种全球方法应用于核反应堆1:3比例的实验性安全壳厂房。最终,研究表明,仅在几十天的计算时间内,就可以获得物理容许的总干空气泄漏率(表示核安全壳结构的性能)及其随时间演变的完整概率分析。这样的结果可以在大型结构的设计、维护和风险评估的决策过程中提供见解和帮助。对于核安全壳建筑(NCB),一个直接的应用是基于运行负荷下泄漏率定义的标准来评估寿命延长。

 

图:应用于同时承受THM荷载的大型钢筋和预应力混凝土结构的全局THM-L建模策略方案。

 

图:全局随机THM-L建模策略的方案。

 

图:HLSA方案应用于交错THM-L计算。

 
 

图:实验1:3比例的双壁NCB–VeRCoRS实体模型(a)内壁(b)外壁(c)内圆顶(d)外圆顶(e)全实体模型(f)内壁结构部件示意图。

 

图:NCB中预应力电缆的类型(a)限制在墙壁高度的水平电缆(b)限制在墙高度的垂直电缆(c)穿过墙壁和圆顶的伽马电缆(d)圆顶电缆。

 

图:将安全壳厂房分解为具有代表性的结构体积。

 

图:最具影响力的输入对(a)角撑板早期峰值温度的全局方差(在局部意义上——第3.1节)的贡献;(b) 每个角撑板的RSV的裂缝数量;(c) 墙在1年和6年时的平均饱和率;(d) 第一次和最后一次加压试验期间通过角撑板裂缝的泄漏率;(e) 第一次和最后一次加压试验期间通过VeRCoRs质量的泄漏率,以及(f)第一次和最近一次加压测试期间通过VeRCoRs模型的总泄漏率。

 

图:总漏气率和部分漏气率的概率描述。(a)平均值随时间的演变;(b) 变异系数和(c)CDF。

文四:

 

预测人腰椎间盘生物力学反应的随机有限元模型分析

摘要:

背景:由于椎间盘复杂的结构和材料特征,很难对其进行准确的生物力学研究。

目的:通过使用随机有限元(FE)模型考虑不同的非线性结构和材料特性,研究IVD生物力学响应的概率分布。

方法:重建L3-4椎间盘的有限元模型,包括髓核、环形基质和纤维。使用蒙特卡罗方法随机生成500组圆盘的非线性材料特性和纤维取向,这些特性和取向被实现到有限元模型中。有限元模型在七种载荷条件下进行了分析:分别为500 N的压缩力、模拟屈曲、伸展、左右侧向弯曲和左右轴向旋转的7.5Nm力矩。分析了椎间盘的运动范围(ROM)、椎间盘内压力(IDP)、纤维应力和基质应变的分布。

结果:在压缩载荷作用下,位移在0.29mm到0.76mm之间变化。在7.5Nm力矩下,ROM在主旋转中在3.0和6.0之间变化。在所有负载条件下,IDP变化在0.3MPa以内。在屈曲和伸展力矩下分别观察到最大纤维应力(3.22±0.64MPa)和基质应变(0.27±0.12%)。

结论:用于构建有限元模型的结构和材料参数会显著影响IVD生物力学。随机有限元模型包括结构和材料参数的概率分布,为分析IVD生物力学反应的统计范围提供了一种有用的方法。

 

图:IVD的有限元模型和纤维环截面示意图。

 

图:前外侧(AO)、后外侧(PO)、前内侧(AI)和后内侧(PI)纤维环纤维位置处纤维的相应响应与实验测量的单层拉伸试验数据的比较。

 

图:随机有限元模型在压缩过程中的力-位移响应。深绿线、红线和黄线之间的区域分别代表90%、70%和50%的病例概率。

 

图:屈曲、伸展、侧向弯曲和轴向旋转的随机有限元模型的椎间盘位移/旋转累积概率(mm/˚)和压缩时椎间盘位移的箱形图(a)和累积概率(b)。

 

图:随机有限元模型在压缩、屈曲、伸展、侧向弯曲和轴向旋转中的椎间盘内压力的(a)箱形图和(b)累积概率(b)。

 

图:随机有限元模型在压缩、弯曲、拉伸、侧向弯曲和轴向旋转中的最大纤维应力(a,b)和最大应变(c,d)的盒图和累积概率。

 

图:将所开发的模型在不同荷载下的响应与实验和有限元在相同荷载条件下的响应进行比较。

文五:

 

磁性填料聚合物复合材料ーー形态角色塑造及力学性能的实验和随机有限元分析

摘要:

含有磁性填料的聚合物复合材料是一种很有前途的材料,用于各种应用,如储能和医疗领域。为了便于各个部件的工程设计,全面了解这种不均匀和潜在的高度各向异性材料的力学行为是很重要的。因此,作者通过压缩成型制造了磁性复合材料。用商业级增稠剂对环氧聚合物基体进行改性。添加各向同性磁性颗粒作为功能填料。使用显微镜技术对微观结构形态,特别是填料的分布、分散和排列进行了表征。采用随机有限元分析方法对复合材料的力学性能进行了实验表征和研究。采用四种情况进行建模以预测弹性模量:完全随机分布、随机排列分布、所谓的“粗糙”界面接触和粘结界面接触。对实验结果和SFEA建模结果进行了比较和讨论。

 

图:(a) SFEA框架的算法,(b)表示蒙特卡罗模拟(MCS)模块的算法的示意图,以及(c)说明随机数生成器(RNG)模块的示意图。

 

图:用于模型生成和分析子过程的有限元分析平台示意图。

 

图:SEM成像磁粉形态的简化数据:(a)磁粉的SEM图像,(b)等效圆盘直径的直方图,以及(c)圆盘厚度的直方图。

 

图:X射线显微镜的示例图像。

 

图:拉伸试样照片:(a)顶面,(b)底面,(c,d)厚度方向的侧视图。

 

图:循环和拉伸力学试验数据摘要: (a)弹性模量,和(b)抗拉强度。误差线表示与平均值的标准偏差。

 

图:RVE内的圆盘颗粒,用于(a)完全随机分布和(b)随机排列的分布情况,去除基质,以及(c)基质嵌入颗粒的剖面图。

来源:STEM与计算机方法
复合材料非线性建筑理论自动驾驶材料储能数字孪生试验人工智能
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-11-01
最近编辑:9月前
江野
博士 等春风得意,等时间嘉许。
获赞 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
联系我们
帮助与反馈