非均质材料力学性能与失效分析的多尺度有限元法研究

来源 :大连理工大学 | 被引量 : 19次 | 上传用户:laobo999
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
自然界中存在的、以及工程应用中的大部分材料都具有多尺度特征,且当考察尺度小到一定程度后,都将表现为非均质性。对这类材料而言,其微观非均质特征对宏观行为影响非常大,实际上,材料在宏观尺度所表现出来的大部分物理力学现象大多是由微观非均质结构决定的。以航空航天工业中广泛使用的复合材料为例,其力学性能与各组分相的材料属性、尺寸、形状、空间分布及界面连接特性等细观特征密切相关。另一方面,材料的宏观失效通常是由微观缺陷、微观裂纹等引起并发展所致。所以,研究微观现象和宏观行为的相互影响已成为工程应用中的一个迫切需要解决的问题。当采用传统数值方法(如有限元法)直接求解这类多尺度特征问题时,常常因为需要大量的计算机存储和计算时间而遇到求解困难。多尺度计算方法作为求解这类问题的一种有效方法,已成为近年来的一个研究热点。本文首先提出了非均质材料弹性力学性能分析的扩展多尺度有限元法(EMsFEM)。基于多尺度有限元法(MsFEM)的思想,数值构造了能反映粗网格单元内部小尺度非均质信息的矢量场多尺度基函数。和标量场问题不同的是,在考虑固体变形的矢量场问题中材料内部会发生体积膨胀/收缩效应(泊松效应),从而导致材料各个方向的变形是相互耦合的。为了处理这种耦合关系,我们提出了一种新型的多尺度基函数构造方式,并在多尺度基函数中添加了耦合附加项。对不同构造方法下的多尺度基函数性能进行比较,结果表明,新型多尺度基函数能正确地反映粗网格单元内部的耦合变形,在相同的边界条件下可得到和常规精细有限元解完全一致的微观变形场,从而很大程度地提高了多尺度方法的计算精度。比较了不同边界条件下多尺度基函数的构造对多尺度方法精度的影响,并分析了尺度效应和误差产生的原因。所发展的方法实施方便,具有较好的数值精度。与传统均匀化方法相比,本文方法无尺度分离和周期性假设的限制,且易于进行降尺度计算以获得微观尺度的真实应力应变信息。其次,基于扩展多尺度有限元法基本思想,提出了一类新的平面矩形单元和空间长方体单元(广义单元)。通过弹性力学伽辽金方程,理论推导出单元内部位移场和节点值之间的函数关系,并应用于新单元形函数的构造中。新单元能更合理地考虑单元内部的耦合变形,在不增加计算自由度的情况下得到比传统单元更精确的结果,为构造精度更高的新型有限单元提供新的思路。进一步,发展了非均质材料小变形弹塑性分析的扩展多尺度有限元法。对于弹塑性材料的多尺度分析,随着微观物理变量的演化,微观尺度节点上将产生不平衡节点力。为了处理这些不平衡节点力产生的响应,我们提出了一种位移分解技术,将微观不平衡节点力响应等效成宏观等效力和微观扰动力的共同作用。其中宏观等效力用来驱动结构的宏观位移响应,而微观扰动力用来获得局部的微观位移扰动响应。在位移分解技术的基础上,进一步提出了两尺度协同迭代算法。在两尺度计算框架下,对宏微观非线性方程组进行协同增量迭代计算。最后,对方法的计算复杂度进行分析并与传统有限元法进行比较。数值结果表明,所发展的方法在模拟非均质材料弹塑性问题时具有很好的精度,且能显著地减少计算机存储需求和计算时间。此外,发展了周期性桁架材料应变局部化分析的自适应多尺度方法(AMM)。局部化现象是引起工程结构破坏的主要原因之一,传统基于体积平均思想的均匀化方法由于受到局部周期性假设和尺度分离假设的限制而不能很好地用于模拟这类问题。在计算过程中,AMM自适应地将整体计算域分解为两个部分,即细网格计算域和粗网格计算域,EMsFEM提供了多尺度框架。这样,变形剧烈的局部化域可采用有限元法在细尺度网格上直接求解,保证了算法的精度;而对于其它变形较缓和的域,则仍采用EMsFEM在粗网格尺度上进行非线性计算。同时,尺度连接通过建立在尺度界面上的主从节点约束方程完成。提出一种基于宏观单元位移梯度的网格识别器来自动更新细网格计算域。数值结果表明AMM具有较高的计算效率,能准确预测局部化问题的剪切带位置,且避免了宏观单元网格依赖性问题。最后,发展了一种模拟准静态裂纹扩展问题的自适应多尺度方法。采用扩展有限元法(XFEM)和水平集法(LSM)的组合对细尺度模型的不连续问题(包括强不连续和弱不连续问题)进行模拟,其中,XFEM用来计算不连续模型的应力应变场,而LSM用来对固定不连续界面(材料界面)和移动不连续界面(裂纹面)进行几何描述。这样,有限元网格可独立于不连续界面,模拟裂纹扩展问题时不需要对有限元网格重新剖分。而对于不需要精细分析的其它子域,采用EMsFEM在粗尺度网格上进行求解。提出一种基于距离参数的自适应网格识别技术来对求解域进行自适应分解,其距离参数可由水平集函数的演化方程自动获取。方法具有很好的精度,能准确地预测裂纹扩展路径,有效地减少裂纹扩展模拟所需的计算自由度和计算时间。
其他文献
目的讨论并分析在小儿哮喘护理中应用健康教育的临床效果。方法选取我院2012年10月至2014年9月期间收治的60例小儿哮喘患者为研究对象,采用随机性分组的方式将患儿平均分为观
总结昆明冶金高等专科学校人才培养的实践经验,提出:要有效培养高技能人才,必须营造有利于企业高技能人才成长的良好氛围,实施高技能人才培养工程,加强校企联合,建立一批高技
目的:从制川乌与白芍组分配伍层面上研究其抗炎、镇痛、解毒、以及治疗类风湿性关节炎的作用,确定其增效的最佳配伍组分和最佳配伍比例,探讨制川乌与白芍组分配伍治疗类风湿
目的:探讨多囊卵巢综合征(polycystic ovary syndrome,PCOS)患者血清抗苗勒管激素(anti-Müllerian hormone,AMH)水平与肥胖、胰岛素抵抗(insulin resistance,IR)程度的相关
实木椅结构形式多样,以往进行家具结构设计并没有一种科学而精确的设计分析方法,文章通过归类整理实木椅类家具结构部件外观设计要素,分析椅类家具使用载荷情况,利用ANSYS有
焊接机器人作为工业机器人应用的一个重要领域,对提高企业工作效率、提升产品质量、降低企业生产成本等方面有着非常重要的意义。本文主要针对焊接机器人的示教盒进行研究。
目的研究新疆医科大学临床医学专业实践教学及其体系,总结实践教学体系的结构特征,分析其不足之处并提出完善体系的建议。方法以新疆医科大学临床医学专业2000年和2005年教学
近些年来,我国各领域的发展逐渐步入互联网时代,国家也开始大力倡导"互联网+"产业的融合发展。网络信息技术是支撑互联网时代稳定发展的重要因素,而网络信息技术自身具备较大
<正>大部分大陆碰撞造山带是由俯冲碰撞的大陆板块和夹杂其中的地体组成,构成了"多地体"岩石圈结构。这些地体之间的缝合带作为构造软弱带、变形易发地,显著影响着造山带的构
会议
【案情简介】李某于2012年12月进入某连锁快餐机构担任店长,2016年8月进人该连锁机构直营事业部任营运助理,其劳动合同期限为2012年12月7日至2018年12月6日。李某于2017年2月