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

非均匀脆性材料数值试验多重随机等效空间 桁架模型 远方，刁可 天津大学建筑工程学院，天津（300072 ） E-mail：yf_CAD@126.com 摘 要：本文提出的数值试验多重随机等效空间桁架模型理论是二维多重随机等效平面桁架 模型在三维空间上的推广，用以模拟混凝土、岩石等非均匀脆性材料的开裂过程，研究分析 材料的变形性质和力学性能。文中给出了空间随机杆件方向余弦常量，推导了等效空间桁架 杆件截面面积计算公式，并通过数值试验证明了方法的有效性，为非均匀脆性材料力学性质 研究提供了一种新的手段。 关键词：多重随机等效空间桁架模型；数值仿真试验；非均匀脆性材料；自动格子机 1．引言 [1] 本文提出的数值试验多重随机等效空间桁架模型是多重随机等效平面桁架模型 在三 维空间上的推广。 以混凝土、岩石为代表的非均匀脆性材料在土木工程中有着非常广泛的应用，但是它们 在荷载作用下的破坏机理却极为复杂，如何建立合理、有效、可靠和使用方便的数理模型， 描述这类材料的受荷破坏过程，至今仍是科研人员为之奋斗而未克服的难题。 非均匀脆性材料的一个重要特点是：细观结构上内部充满空隙、孔洞、裂纹、杂质等各 种缺陷，物质分布极不均匀；另一个重要特点是：受力破坏过程中，材料内部伴随大量微损 伤和微破坏。一般的破坏过程是：随着应力增加，微损伤不断产生，积累；微裂缝相互连接， 贯穿，形成大的裂缝；同时，材料刚度逐渐减小，直到完全丧失承载能力。 传统的弹塑性力学以及断裂力学、损伤力学都以连续介质模型为基础。显然，以此分析 非均匀脆性材料不符合假设前提，除非将研究阶段限定在材料变形的近似线弹性阶段，否则 必须考虑材料的非均匀性影响。为了研究分析非均匀脆性材料变形进入非线性阶段和承载力 下降段的力学特征，近年来细观力学蓬勃发展起来[2] [3][4][5][6] 。细观力学以各种计算模型 为基 础，用随机概率方法刻画材料在分布和强度上的不均匀性，通过细观层次上的本构关系描述 宏观层次上的力学性能。细观力学在一定程度上取得了成功，科研工程人员可以以此为工具 比较准确地分析研究非均匀脆性材料的各种力学特性。然而细观力学对材料性质的描述是建 立在细观尺度上的，对于那些相对于细观尺度比较大的形体，单元剖分数量必然很大，需要 占用相当大的计算资源，因而，限制了细观力学的广泛应用。 2000年Brian Tatting 和Zafer Gurdal将自动格子机原理用于结构仿真分析，提出了自动格 子机二维结构分析方法[7][8] 。文中作者应用有限元理论对材料进行了几何剖分，然后将剖分 单元质点化，用质点代替剖分单元，质点间用杆件相连，形成平面桁架，并以此平面桁架模 拟连续材料介质的力学性能，见图1 （a ）。显然，如果一个剖分单元以及相邻剖分单元间的 物理联系仅用一个质点和一根杆件来代表，不足以描述材料内部物质间的相互关系，除非剖 分单元足够小，小到细观尺寸水平，如此已与细观力学计算模型无异。为了能够在细观尺度 与宏观尺度之间的任意尺度上构建反映材料特性的计算模型，本文采用多个质点代表一个剖 分单元，并用多个杆件代表相邻单元间的物理联系，见图1 （b ）；通过增加相邻单元间的物 理联系，而不是增加剖分单元数量的方法，给出了描述非均匀脆性材料的结构组成和力学性 能的一种新方法。 -1- 剖分单元 单元间杆件 单元质点 (a) (b) 图 1 单质点单