非均匀脆性材料数值试验多重随机等效空间桁架模型.docxVIP

-1-

非均匀脆性材料数值试验多重随机等效空间桁架模型

远方，刁可

天津大学建筑工程学院，天津（300072）

摘要：本文提出的数值试验多重随机等效空间桁架模型理论是二维多重随机等效平面桁架模型在三维空间上的推广，用以模拟混凝土、岩石等非均匀脆性材料的开裂过程，研究分析材料的变形性质和力学性能。文中给出了空间随机杆件方向余弦常量，推导了等效空间桁架杆件截面面积计算公式，并通过数值试验证明了方法的有效性，为非均匀脆性材料力学性质研究提供了一种新的手段。

关键词：多重随机等效空间桁架模型；数值仿真试验；非均匀脆性材料；自动格子机

1．引言

本文提出的数值试验多重随机等效空间桁架模型是多重随机等效平面桁架模型[1]在三维空间上的推广。

以混凝土、岩石为代表的非均匀脆性材料在土木工程中有着非常广泛的应用，但是它们在荷载作用下的破坏机理却极为复杂，如何建立合理、有效、可靠和使用方便的数理模型，描述这类材料的受荷破坏过程，至今仍是科研人员为之奋斗而未克服的难题。

非均匀脆性材料的一个重要特点是：细观结构上内部充满空隙、孔洞、裂纹、杂质等各种缺陷，物质分布极不均匀；另一个重要特点是：受力破坏过程中，材料内部伴随大量微损伤和微破坏。一般的破坏过程是：随着应力增加，微损伤不断产生，积累；微裂缝相互连接，贯穿，形成大的裂缝；同时，材料刚度逐渐减小，直到完全丧失承载能力。

传统的弹塑性力学以及断裂力学、损伤力学都以连续介质模型为基础。显然，以此分析非均匀脆性材料不符合假设前提，除非将研究阶段限定在材料变形的近似线弹性阶段，否则必须考虑材料的非均匀性影响。为了研究分析非均匀脆性材料变形进入非线性阶段和承载力下降段的力学特征，近年来细观力学蓬勃发展起来[2]。细观力学以各种计算模型[3][4][5][6]为基础，用随机概率方法刻画材料在分布和强度上的不均匀性，通过细观层次上的本构关系描述宏观层次上的力学性能。细观力学在一定程度上取得了成功，科研工程人员可以以此为工具比较准确地分析研究非均匀脆性材料的各种力学特性。然而细观力学对材料性质的描述是建立在细观尺度上的，对于那些相对于细观尺度比较大的形体，单元剖分数量必然很大，需要占用相当大的计算资源，因而，限制了细观力学的广泛应用。

2000年BrianTatting和ZaferGurdal将自动格子机原理用于结构仿真分析，提出了自动格子机二维结构分析方法[7][8]。文中作者应用有限元理论对材料进行了几何剖分，然后将剖分单元质点化，用质点代替剖分单元，质点间用杆件相连，形成平面桁架，并以此平面桁架模拟连续材料介质的力学性能，见图1（a）。显然，如果一个剖分单元以及相邻剖分单元间的物理联系仅用一个质点和一根杆件来代表，不足以描述材料内部物质间的相互关系，除非剖分单元足够小，小到细观尺寸水平，如此已与细观力学计算模型无异。为了能够在细观尺度与宏观尺度之间的任意尺度上构建反映材料特性的计算模型，本文采用多个质点代表一个剖分单元，并用多个杆件代表相邻单元间的物理联系，见图1（b）；通过增加相邻单元间的物理联系，而不是增加剖分单元数量的方法，给出了描述非均匀脆性材料的结构组成和力学性能的一种新方法。

-2-

剖分单元单元间杆件

单元质点

(a)

(b)

图1单质点单杆件与多质点多杆件的对比

文中第二、三部分阐述了多重随机等效空间桁架模型的构造方法和杆件截面计算公式以及使用该模型进行数值仿真的限制条件。第四部分以混凝土材料为背景说明了模型的使用情况。结论部分对全文作了总结，并提出了进一步需要开展的研究工作。

2．多重随机等效空间桁架模型的几何构成

由基本等效空间桁架结构构造多重随机等效空间桁架结构。

2.1基本等效空间桁架结构

设有一定体积的混凝土，将混凝土按立方体网格剖分，由每一立方体的中心出发，分别与其相邻的26个立方体中心相连，将连线看作空间二力杆件，形成基本等效空间桁架结构，见图2。在一个剖分中，分别取由中心点出发向外辐射杆件的一半，形成等效空间桁架单元，见图3。剖分中杆件的汇聚点称为等效空间桁架结构节点。

立方体单元等效桁架单元图3微元体与桁架单元等效关系图

立方体单元

等效桁架单元

图3微元体与桁架单元等效关系图

图2混凝土按立方体剖分，形成基本等效桁架结构

2.2等效空间桁架结构节点的随机扰动

为了反映材料的不均匀性，将节点在各自所在立方体范围内作随机扰动，形成随机等效空间桁架结构，见图4。其中边界节点的处理规则为：顶点节点保持不动，棱线上节点在棱线上作随机扰动，各侧面上节点在面内作随机扰动。

2.3等效空间桁架结构的多重叠加

参照文献[1]介绍的幻化概念，对随机等效空间桁架