周期性点阵桁架材料力学性能分析的一种新的多尺度计算方法.pdfVIP

第32卷第2期 固体力学学报 V01．32No．2 CHINESE OFSOLID 2011年4月 JOURNAL MECHANICS 2011 April 周期性点阵桁架材料力学性能分析 的一种新的多尺度计算方法’ 张洪武“ 吴敬凯付振东 (大连理工大学工业装备结构分析国家重点实验室，运载7-程与力学学部，工程力学系，大连，116024) 摘要论文提出一种周期性点阵桁架材料力学性能分析的新的多尺度方法．方法的主要思想是通过数值构 造能反映周期性桁架材料单胞内部非均质性的多尺度基函数，从而在大尺度上求得单胞的等效刚度阵，大大减小 了模型计算量．通过引入基函数的耦合附加项，以考虑多维矢量场问题不同方向问的耦合作用．数值研究表明，采 用线性边界条件构造基函数有时会产生较强的边界效应，而超样本技术的振荡边界和周期性边界条件能很大程度 地减少单胞边界强制变形产生的误差．特别是对于非均质特征尺度跟宏观单元尺度相近的单胞，论文提出的周期 性边界条件具有很好的效果．所提出的方法的优点是，在所构造的多尺度有限元法的基础上，能较容易地进行降尺 度计算，较准确地得到单胞内部小尺度上的应力应变信息，为材料强度分析的多尺度计算打下基础． 关键词 多尺度计算方法，桁架材料，降尺度计算，基函数 度上对问题精确求解而直接求得问题的大尺度解， 0 引言 大大降低了计算量． 尽管多尺度有限元方法在求解标量场问题中已 具有周期性微观结构的点阵桁架材料由于质量 经得到成功应用，但是对于求解固体变形的矢量场 轻、比强度、比刚度高等优点，在土木工程和航空航 问题的成果国内外还少见报道．其难点之一是如何 天材料结构等领域的运用越来越广泛．通常的有限 利用该方法处理固体力学问题中的体积膨胀／收缩 元技术可用于求解这种周期点阵结构，但随着单元 及泊松效应等问题，以便获得较高精度的结果． 数量的增长，计算成本将急剧增加，甚至不可行．解 本文针对周期性桁架材料的力学性能分析问 决问题的方法之一是通过均匀化方法将原结构连续 题，提出了扩展多尺度有限元法(EMsFEM)．通过 体化，由此可采用较大尺寸的连续体有限单元离散 引入基函数的耦合附加项，以考虑多维矢量场问题 结构，大大提高计算效率．目前关于桁架材料结构的 不同方向间变形的耦合作用．同时建议采取三类边 等效性能研究已经开展了大量的工作[1‘7]．如文献界条件进行数值基函数构造．数值结果发现，基函数 E3]总结了一些常用的连续体化方法，阐述了它们的 构造时所施加的边界条件对多尺度计算结果的精度 一些优点，并证明这些等效方法的适用性．文献[7] 具有较大的影响，线性边界条件下有时会产生较强 比较了Dirichlet型和Neumann型边界条件下的代的边界效应，而由超样本技术得到的振荡边界条件 表体元法及均匀化法对具有周期性微观结构的点阵 以及周期边界条件能有效地降低单胞边界强制变形 类桁架材料等效性能的预测结果，同时对尺度效应 所引起的误差，使结果更加精确．特别地，对于非均 进行了讨论．近年来国内外学者针对渗流等标量场 质特征尺度跟宏观单元尺度相近的单胞，本文提出 问题求解提出了多尺度有限元方法MsFEM(Mul—的周期性边界条件能获得很好的结果．同其它均匀 tiscaleFiniteElement Method)[8。0。，其基本思想是化方法一样，本文方法能显著地减少计算量和节约 通过数值构造满足局部特