纤维混凝土纤维桥联数值分析模型的构建方法

1.一种纤维混凝土纤维桥联数值分析模型的构建方法，其特征在于，包括以下步骤：步骤S1、建立纤维与基体界面的非线性滑移本构模型，确定界面剪切应力与滑移量之间的关系；步骤S2、生成三维纤维方向角和埋置长度的随机分布场；步骤S3、将微观界面的非线性滑移本构模型与宏观相场法结合，通过粘聚力单元表征界面行为，并基于纤维方向角和埋置长度的随机分布场，利用异构数据传递协议将界面脱粘行为映射至宏观断裂面；步骤S4、在多尺度分析中，基于后验误差估计方法，采用非结构化六面体网格局部细分技术，对裂缝尖端区域的网格尺寸进行动态调整；步骤S1中，非线性滑移本构模型为分段双线性软化模型，其表达式为： ；其中， 表示界面剪切应力， 表示滑移量， 表示峰值剪切应力， 表示达到峰值剪切应力时的滑移量， 表示软化阶段的残余滑移量；步骤S3中，异构数据传递协议通过如下步骤实现：步骤S31、在微观尺度上，获取粘聚力单元的脱粘位移场和应力场；步骤S32、利用插值算法将微观尺度的脱粘位移场和应力场映射到宏观断裂面上；步骤S33、根据映射结果，更新宏观断裂面的力学状态，包括位移和应力；步骤S34、将更新后的宏观断裂面力学状态反馈至微观尺度，实现微观与宏观尺度的耦合；异构数据传递协议的具体表达式为： ；其中， 表示宏观断裂面上点 处的位移， 表示微观界面脱粘处第 个粘聚力单元的位移， 表示插值函数， 表示微观处的位移梯度， 表示第 个粘聚力单元的坐标。 2.根据权利要求1所述的一种纤维混凝土纤维桥联数值分析模型的构建方法，其特征在于，步骤S2中，三维纤维方向角 按如下概率密度分布 ： ；其中， 表示极角， 表示方位角， 表示各向异性系数， 表示第二阶勒让德多项式，表达式为： 。 3.根据权利要求2所述的一种纤维混凝土纤维桥联数值分析模型的构建方法，其特征在于，步骤S2中，通过修正因子 修正纤维方向角分布函数，修正因子 的表达式为： ；其中， 表示单位体积内纤维数量， 表示纤维直径， 表示混凝土密度； 表示纤维埋置长度；修正后的纤维方向角分布函数为： 。 4.根据权利要求3所述的一种纤维混凝土纤维桥联数值分析模型的构建方法，其特征在于，步骤S2中，纤维埋置长度遵循均值为 、标准差为 的对数正态分布： ；其中， 表示纤维埋置长度为 时的概率密度。 5.根据权利要求4所述的一种纤维混凝土纤维桥联数值分析模型的构建方法，其特征在于，步骤S3中，粘聚力单元的力学参数通过如下优化公式确定： ； ；其中， 表示断裂能， 表示拉拔试验中的峰值荷载， 表示临界脱粘位移， 表示试件宽度， 表示峰值粘结应力， 表示纤维弹性模量；同时，引入数值反分析修正系数 ，通过对比数值模拟与试验的荷载-位移曲线，迭代优化 使模拟结果与试验结果的均方根误差RMS低于5%。 6.根据权利要求1所述的一种纤维混凝土纤维桥联数值分析模型的构建方法，其特征在于，步骤S4中，非结构化六面体网格局部细分技术包括以下步骤：步骤S41、确定需要细分的裂缝尖端区域；步骤S42、基于后验误差估计结果，计算该区域所需的网格密度；步骤S43、采用非结构化六面体网格生成算法，对目标区域进行网格细分；步骤S44、对细分后的网格进行质量检查和优化，确保网格质量满足数值分析要求。 7.根据权利要求6所述的一种纤维混凝土纤维桥联数值分析模型的构建方法，其特征在于，后验误差估计采用如下能量范数 形式： ；其中， 表示数值计算得到的应力场， 表示基于试验或高精度基准解的应力场， 表示分析区域；当 超过预设阈值 时，触发非结构化六面体网格在裂缝尖端区域的细分，细分准则为： ；其中， 和 分别表示细分后和细分前的网格尺寸， 表示局部误差估计值， 表示收敛阶数，取值范围为[1,3]。