1.一种轴心受力状态下水泥净浆早期宏观弹性模量计算方法，其特征在于：基于已知水泥净浆微观模型和各类水化物材料性能的基础上，借助有限单元仿真技术，计算轴心受力状态下微观模型各单元平均应变，并结合弹性模量基本定义计算水泥浆体微观模型的平均弹性模量。

2.如权利要求1所述的轴心受力状态下水泥净浆早期宏观弹性模量计算方法，其特征在于：具体在水泥水化各阶段中，提取含有水泥净浆内部微观结构组成和分布信息的水泥净浆微观模型，于模型顶面施加微小均布荷载，通过有限元基本理论，计算模型各单元应变值，得到模型的平均应变；最后根据弹性模量的基本定义求解各水化阶段水泥净浆弹性模量，分析水泥净浆早期宏观弹性模量发展规律。

3.如权利要求2所述的轴心受力状态下水泥净浆早期宏观弹性模量计算方法，其特征在于：轴心受力构件弹性模量计算方法如下：在宏观层次上，水泥净浆可视为均质材料，在轴心受力状态下，根据材料力学中弹性模量的定义：E=σ/ε （7）

式中，σ为试件的平均应力，MPa；ε为试件的平均应变。

本发明提及的水泥净浆有限元微观模型由低密度C-S-H单元、高密度C-S-H单元、CH单元、水化硫铝酸钙单元、毛细孔单元、自由水单元以及未水化水泥颗粒单元组成，各类单元弹性模量不同；因此在轴心荷载作用下，各单元的应变不同。

为计算整个试件的平均应变，可取试件1/4高度处截面、1/2高度处截面、3/4高度处截面以及顶面各单元的应变进行加权平均，得到各截面处平均应变式中，ε1、ε1/4、ε1/2和ε3/4分别为1/4高度处截面、1/2高度处截面、3/4高度处截面以及顶面的平均应变； 和 分别为1/4高度处截面、1/2高度处截面、3/4高度处截面以及顶面第i单元的平均应变；n1、n1/4、n1/2和n3/4分别为1/4高度处截面、1/2高度处截面、3/4高度处截面以及顶面的单元个数；

对于轴心受力构件任意截面的平均应力应该等于均布外荷载q，将式（8）代入式（7）可得1/4高度处截面、1/2高度处截面、3/4高度处截面以及顶面的弹性模量，如式（9）所示：E1/4=q/ε1/4 E1/2=q/ε1/2 E3/4=q/ε3/4 E1=q/ε1 (9)对以上弹性进行再次加权平均，可得试件整体弹性模量：

4.如权利要求3所述的轴心受力状态下水泥净浆早期宏观弹性模量计算方法，其特征在于：具体采用以下步骤实现：步骤1，输入计算起始、终止时间和计算步长，计算总次数n=（终止时间-起始时间）/计算步长；

步骤2，输入第一计算步当前时刻的水泥净浆微观模型，模型中已将各单元划分成低密度C-S-H单元、高密度C-S-H单元、CH单元、水化硫铝酸钙单元、毛细孔单元、自由水单元以及未水化水泥颗粒单元；

步骤3，输入低密度C-S-H单元、高密度C-S-H单元、CH单元、水化硫铝酸钙单元自由水单元以及未水化水泥颗粒单元的弹性模量，将毛细孔单元单元弹性模量设为1；

步骤4，模型顶面施加一微小均布荷载，底部施加竖向约束；

步骤5，通过常用有限元分析软件进行位移场、应力场求解，得到模型各单元的平均应变值；

步骤6，输出试件1/4高度处截面各单元的应变值，根据式

计算试件1/4高度处截面的截面平均应变；

步骤7，对试件1/4高度处截面各单元循环，验证各单元应变值的有效性：若某单元的[（应变-截面平均应变）/截面平均应变]小于±30%，认为数据有效，反之将数据剔除；

步骤8，对有效的单元应变值按式（8）再次进行加权平均，得到试件1/4高度处截面的最终截面平均应变；

步骤9，按式

E1/4=q/ε1/4 E1/2=q/ε1/2 E3/4=q/ε3/4 E1=q/ε1 (9)计算1/4高度处截面的平均弹性模量；

步骤10，重复步骤6～8，得到1/2高度处截面、3/4高度处截面以及顶面的平均弹性模量；

步骤11，按式

计算模型平均弹性模量；

步骤12，重复步骤2～11，计算第2计算步直至第n计算步内模型的平均弹性模量；

步骤13，以各计算步对应的时刻为横坐标，以计算得到的各计算步内模型的平均弹性模量为纵坐标，绘制水泥净浆早期宏观弹性模量发展曲线。