An automatic full 3D mesh generator for the non-linear 3D homogenized discrete FE analysis of masonry double curvature structures is presented. The method can reproduce the main features exhibited by masonry, preserving the orthotropy and accounting for the influence of membrane loads on the flexural behavior. Previous homogenization procedures allow reducing the number of variables, leading to non-linear analyses with a low computational burden. In the approach proposed, masonry is modeled through an assemblage of elastic cells joined by non-linear interfaces with finite thickness. The latter are modeled as brick elements with softening replicated by a damage constitutive relationship (known in Abaqus - the commercial FE code used - as Concrete Damage Plasticity). The discrete mesh is obtained automatically by means of an ad-hoc routine which transforms a standard 2D homogeneous FE mesh into a 3D orthotropic elastic model with non-linear interfaces connecting contiguous elastic elements (cells). The proposed method is applied to a hemispherical dome and a cloister vault, with abundant experimental and numerical data available in the literature to compare with. Sensitivity analyses are conducted in order to have an insight into the reliability of the method, which proved excellent predictivity capabilities, as far as both the peak load and the activated failure mechanism are concerned.

提出了一种自动全3D网格生成器，用于砌体双曲率结构的非线性3D均质化离散有限元分析。该方法可以重现砌体所表现出的主要特征，保留正交性并考虑膜荷载对挠曲行为的影响。以前的均质化程序允许减少变量的数量，从而导致非线性分析且计算负担低。在提出的方法中，砌筑是通过由有限厚度的非线性界面连接的弹性单元的组合建模的。后者被建模为具有破坏本构关系的软化砖元素（在Abaqus中使用-商业FE代码-称为混凝土破坏可塑性）。离散网格是通过ad-hoc例程自动获得的，该例程将标准的2D均匀有限元网格转换为3D正交各向异性弹性模型，该模型具有连接相邻弹性元件（单元）的非线性界面。该方法应用于半球形穹顶和回廊，具有丰富的实验和数值数据可供参考。进行灵敏度分析是为了深入了解该方法的可靠性，就峰值负载和激活的失效机理而言，该方法具有出色的预测能力。该方法应用于半球形穹顶和回廊，具有丰富的实验和数值数据可供参考。进行灵敏度分析是为了深入了解该方法的可靠性，就峰值负载和激活的失效机制而言，该方法具有出色的预测能力。该方法应用于半球形穹顶和回廊，具有丰富的实验和数值数据可供参考。进行灵敏度分析是为了深入了解该方法的可靠性，就峰值负载和激活的失效机理而言，该方法具有出色的预测能力。