CN103439177B - 轴心受力状态下水泥净浆早期宏观弹性模量计算方法 - Google Patents
轴心受力状态下水泥净浆早期宏观弹性模量计算方法 Download PDFInfo
- Publication number
- CN103439177B CN103439177B CN201310401803.3A CN201310401803A CN103439177B CN 103439177 B CN103439177 B CN 103439177B CN 201310401803 A CN201310401803 A CN 201310401803A CN 103439177 B CN103439177 B CN 103439177B
- Authority
- CN
- China
- Prior art keywords
- cement paste
- elastic modulus
- cement
- micromodel
- early days
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提供了一种轴心受力状态下水泥净浆早期宏观弹性模量计算方法,在已知水泥净浆微观模型和各类水化物材料性能的基础上,为计算轴心受力状态下水泥净浆早期宏观弹性模量提供一种概念直观,易于掌握使用有限元分析方法。基于已知水泥净浆微观模型和各类水化物材料性能的基础上,借助有限单元仿真技术,计算轴心受力状态下微观模型各单元平均应变,并结合弹性模量基本定义计算水泥浆体微观模型的平均弹性模量,是水泥基材料微观力学理论和有限元技术的结合创新和扩展应用。本发明可得到水泥净浆早期宏观弹性模量发展规律并为分析水泥水化产物组成、分布以及孔隙率、孔径分布对水泥净浆宏观弹性模量的影响奠定基础。本发明算法简单、概念清晰、易于使用者掌握和扩展。
Description
技术领域
本发明属于水泥基材料早期力学性能仿真技术领域,具体为一种基于水泥净浆微观模型和有限元仿真技术的早期宏观弹性模量计算方法。
背景技术
研究表明水泥基材料的微观结构对其宏观行为有着决定性作用,对于某一种水泥基材料,其强度、弹性模量、导热系数等宏观性能的发展与其内部微观结构密切联系。因此只有从本质上认识水泥基材料,理解水泥水化过程中材料内部微观结构的组成、分布及其对宏观性能的作用等一套完整的知识体系,才能控制水泥基材料的性能,进而开发高性能材料、绿色环保材料及基于性能设计的材料,为工程应用服务。
将微观力学引入水泥浆体、水泥砂浆和混凝土等水泥基复合材料的研究是对水泥基材料性能研究的重大突破。微观力学(micromechanics)理论基于各相几何和物理特性来研究复合材料的物理性能。其中,预测弹性模量的微观力学模型常采用解析法,包括Mori-Tanaka理论、Hashin-Shtrikman上下限理论、自洽法和改进直接法。解析法的应用使从结构组成出发预测多相多尺度水泥基材料的宏观弹性模量成为现实。但目前另一种常用分析方法——有限元理论在水泥基材料微观结构分析中应用较少。而有限元理论具有解析法所不具备的概念直观,易于掌握使用的优点。
因此探讨有限元理论在水泥基材料微观力学领域的应用,通过有限单元法预测水泥基材料宏观力学性能,是一种新方法的尝试。对于开展水泥基材料早期宏观性能发展机理研究具有十分重要的现实意义。
发明内容
本发明的目的在于提供一种轴心受力状态下水泥净浆早期宏观弹性模量计算方法,在已知水泥净浆微观模型和各类水化物材料性能的基础上,为计算轴心受力状态下水泥净浆早期宏观弹性模量提供一种概念直观,易于掌握使用有限元分析方法。
为实现上述目的,本发明采用以下技术方案:
一种轴心受力状态下水泥净浆早期宏观弹性模量计算方法,基于已知水泥净浆微观模型和各类水化物材料性能的基础上,借助有限单元仿真技术,计算轴心受力状态下微观模型各单元平均应变,并结合弹性模量基本定义计算水泥浆体微观模型的平均弹性模量。
具体在水泥水化各阶段中,提取含有水泥净浆内部微观结构组成和分布信息的水泥净浆微观模型,于模型顶面施加微小均布荷载,通过有限元基本理论,计算模型各单元应变值,得到模型的平均应变;最后根据弹性模量的基本定义求解各水化阶段水泥净浆弹性模量,分析水泥净浆早期宏观弹性模量发展规律。
水泥净浆微观模型包括低密度C-S-H单元、高密度C-S-H单元、CH单元、水化硫铝酸钙单元、毛细孔单元、自由水单元以及未水化水泥颗粒单元,全面考虑水泥净浆水化过程中微观结构的组成和分布。
以下针对本发明方法做进一步说明,具体内容如下:
一、水泥净浆微观结构
水泥净浆微观结构组成可以分为:水化产物、孔、水和未水化的水泥颗粒。
(1)水化产物:主要包括水化硅酸钙凝胶C-S-H、CH和水化硫铝酸钙。
其中水化硅酸钙凝胶C-S-H为最主要的水化产物,占水化水泥浆体体积的50~60%,是水泥浆体中决定硬化水泥浆体物理结构和性能的主要成分。根据C-S-H中凝胶孔隙率大小,C-S-H可分为低密度C-S-H(LD C-S-H,凝胶孔隙率为37%)和高密度C-S-H(HD C-S-H,凝胶孔隙率为24%),低密度C-S-H和高密度C-S-H的弹性模量略有不同。
氢氧化钙CH是水泥净浆的第二大水化产物,最新的研究表明氢氧化钙占水化水泥浆体体积的15%。
水化硫铝酸钙是水泥浆体的次要组成成分,占水化水泥浆体体积的10~15%。
(2)孔:孔隙是水泥浆体的主要组成成分之一,其微观结构对水泥浆体的物理力学性能产生重要的影响。孔隙分类两类:毛细孔和凝胶孔,毛细孔对水泥浆体的强度、渗透性有一定影响。凝胶孔被认为是C-S-H的一部分,位于C-S-H所占的体积内部。通常认为孔隙率是决定强度的最为重要的单因素。
(3)一般情况下水泥浆体中的水分为三种状态:化学结合水、物理吸附水和自由水。其中化学结合水是水化反应形成的水化产物的一部分;物理吸附水是指存在于凝胶中的物理水,和化学结合水一样,都属于水化硅酸钙凝胶C-S-H的一部分,对浆体结构和性能无直接影响;自由水是水泥浆体中尚未参与水化反应或反应后多余的水,宏观性能的产生、发展和变化具有重要的作用。
(4)残留未水化水泥颗粒:水泥颗粒的粒径普遍在0.1~100μm范围内,平均粒径约为10~15μm,其中较大的水泥颗粒很难完全水化,因此在水化过程中水泥浆体内不可避免地含有一定量残留的未水化水泥颗粒。
基于水泥净浆微观结构的组成,本发明涉及水泥净浆微观有限元模型中包括低密度C-S-H单元、高密度C-S-H单元、CH单元、水化硫铝酸钙单元、毛细孔单元、自由水单元以及未水化水泥颗粒单元。
二、有限元基本理论
本发明涉及有限元模型中单元采用空间8节点六面体单元,如附图1所示。
有限单元法根据虚功原理建立的非线性有限元平衡方程为:
[K]{δ}={P} (1)
式中,[K]为结构整体刚度矩阵;{δ}为节点位移列阵;{P}为节点荷载列阵。
求解方程(1)可得模型各节点位移矩阵{δ},然后进行单元分析,即根据单元的节点位移列阵{δ}e确定单元的位移分量列阵{f}、应变分量列阵{ε}。
空间8节点等参单元的位移模式和坐标变换式分别为
以及
式中,ui、vi、wi以及xi、yi、zi(i=1,2,…,8)分别为8个节点的节点位移分量与节点整体坐标分量,Ni为形函数,采用单元的局部坐标ξ、η、ζ表示,即得
将式(4)代入式(2)得
根据弹性力学的几何方程,应变的计算公式如下
{ε}=[B]{δ}e (6)
式中,[B]为应变矩阵。
三、轴心受力构件弹性模量计算
在宏观层次上,水泥净浆视为均质材料,在轴心受力状态下,根据材料力学中弹性模量的定义:
E=σ/ε (7)
式中,σ为试件的平均应力,MPa;ε为试件的平均应变。
本发明提及的水泥净浆有限元微观模型由低密度C-S-H单元、高密度C-S-H单元、CH单元、水化硫铝酸钙单元、毛细孔单元、自由水单元以及未水化水泥颗粒单元组成,各类单元弹性模量不同;因此在轴心荷载作用下,各单元的应变不同;
为计算整个试件的平均应变,取试件1/4高度处截面、1/2高度处截面、3/4高度处截面以及顶面各单元的应变进行加权平均,得到各截面处平均应变
式中,ε1/4、ε1/2、ε3/4和ε1分别为1/4高度处截面、1/2高度处截面、3/4高度处截面以及顶面的平均应变;和分别为1/4高度处截面、1/2高度处截面、3/4高度处截面以及顶面第i单元的平均应变;n1/4、n1/2、n3/4和n1分别为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)
对以上弹性进行再次加权平均,可得试件整体弹性模量。
本发明的有益效果是:本发明在已知水泥净浆微观模型和各类水化物材料性能的基础上,为计算轴心受力状态下水泥净浆早期宏观弹性模量提供一种概念直观,易于掌握使用有限元分析方法。
附图说明
图1为空间8节点六面体单元示意图;
图2为水泥净浆微观数字模型生成流程图。
具体实施方式
本发明通过探讨有限单元法在水泥基材料早期宏观性能预测领域的研究,提出了一种轴心受力状态下水泥净浆早期宏观弹性模量计算方法。本发明主要步骤流程图如图2所示,以下具体结合附图说明本发明的具体实施方法:
步骤1,输入计算起始、终止时间和计算步长,计算总次数n=(终止时间-起始时间)/计算步长;
步骤2,输入第一计算步当前时刻的水泥净浆微观模型,模型中已将各单元划分成低密度C-S-H单元、高密度C-S-H单元、CH单元、水化硫铝酸钙单元、毛细孔单元、自由水单元以及未水化水泥颗粒单元;
步骤3,输入低密度C-S-H单元、高密度C-S-H单元、CH单元、水化硫铝酸钙单元自由水单元以及未水化水泥颗粒单元的弹性模量,将毛细孔单元单元弹性模量设为1;
步骤4,模型顶面施加一微小均布荷载,底部施加竖向约束;
步骤5,通过常用有限元分析软件进行位移场、应力场求解,得到模型各单元的平均应变值;
步骤6,输出试件1/4高度处截面各单元的应变值,根据式(8)计算试件1/4高度处截面的截面平均应变;
步骤7,对试件1/4高度处截面各单元循环,验证各单元应变值的有效性:若某单元的[(应变—截面平均应变)/截面平均应变]小于±30%,认为数据有效,反之将数据剔除;
步骤8,对有效的单元应变值按式(8)再次进行加权平均,得到试件1/4高度处截面的最终截面平均应变;
步骤9,按式(9),计算1/4高度处截面的平均弹性模量;
步骤10,重复步骤6~8,得到1/2高度处截面、3/4高度处截面以及顶面的平均弹性模量;
步骤11,按式(10)计算模型平均弹性模量;
步骤12,重复步骤2~11,计算第2计算步直至第n计算步内模型的平均弹性模量;
步骤13,以各计算步对应的时刻为横坐标,以计算得到的各计算步内模型的平均弹性模量为纵坐标,绘制水泥净浆早期宏观弹性模量发展曲线。
Claims (1)
1.一种轴心受力状态下水泥净浆早期宏观弹性模量计算方法,其特征在于:基于已知水泥净浆微观模型和各类水化物材料性能的基础上,借助有限单元仿真技术,在水泥净浆有限元微观模型上选取试件1/4高度处截面、1/2高度处截面、3/4高度处截面以及顶面各单元的应变进行加权平均,并通过均布外荷载q获取各个截面上的弹性模量,然后加权平均获得的试件整体的弹性模量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310401803.3A CN103439177B (zh) | 2013-09-06 | 2013-09-06 | 轴心受力状态下水泥净浆早期宏观弹性模量计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310401803.3A CN103439177B (zh) | 2013-09-06 | 2013-09-06 | 轴心受力状态下水泥净浆早期宏观弹性模量计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103439177A CN103439177A (zh) | 2013-12-11 |
CN103439177B true CN103439177B (zh) | 2016-12-07 |
Family
ID=49692877
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310401803.3A Active CN103439177B (zh) | 2013-09-06 | 2013-09-06 | 轴心受力状态下水泥净浆早期宏观弹性模量计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103439177B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103983547B (zh) * | 2014-05-23 | 2016-06-29 | 重庆大学 | 一种细长金属管粉料夯实过程粉料颗粒尺度的划分方法 |
CN104062224B (zh) * | 2014-06-19 | 2016-05-11 | 合肥工业大学 | 一种颗粒污泥弹性模量的测量装置及其测量方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1069411C (zh) * | 1996-04-26 | 2001-08-08 | 柳州铁路局工程处第一建筑工程段 | 评定混凝土构件质量的弹性模量法 |
US6301970B1 (en) * | 1998-08-06 | 2001-10-16 | The United States Of America The Secretary Of The Navy | Cumulative damage model for structural analysis of filed polymeric materials |
CN1696646A (zh) * | 2005-05-23 | 2005-11-16 | 喻勇 | 采用立方体试样测试脆性材料抗拉强度的改进方法 |
US7945432B2 (en) * | 2008-11-26 | 2011-05-17 | Livermore Software Technology Corporation | Spot weld failure determination method in a finite element analysis |
CN101923024B (zh) * | 2010-08-04 | 2012-01-25 | 重庆大学 | 一种确定带裂钢筋混凝土梁弹性模量及抗弯刚度的方法 |
CN102116716A (zh) * | 2010-08-10 | 2011-07-06 | 中国水利水电第十四工程局有限公司 | 一种测定塑性混凝土弹性模量的试验方法 |
CN102636634B (zh) * | 2012-01-20 | 2014-08-27 | 中国寰球工程公司 | 一种混凝土构件温度效应确定方法 |
CN103226076A (zh) * | 2013-04-01 | 2013-07-31 | 河海大学 | 一种水泥基材料早期弹性模量测量方法及测量装置 |
-
2013
- 2013-09-06 CN CN201310401803.3A patent/CN103439177B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN103439177A (zh) | 2013-12-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Fu et al. | Assessment of the elastic properties of amorphous calcium silicates hydrates (I) and (II) structures by molecular dynamics simulation | |
Benedetti et al. | A grain boundary formulation for crystal plasticity | |
Hunter et al. | Large-scale 3D phase field dislocation dynamics simulations on high-performance architectures | |
Shen et al. | Predicting structure and energy of dislocations and grain boundaries | |
Kamali-Bernard et al. | Effect of tensile cracking on diffusivity of mortar: 3D numerical modelling | |
Li et al. | Modeling the apparent and intrinsic viscoelastic relaxation of hydrating cement paste | |
WO2021218070A1 (zh) | 基于rev全区域覆盖的岩体工程跨尺度模拟计算方法 | |
Tian et al. | Dynamic crushing behavior and energy absorption of hybrid auxetic metamaterial inspired by Islamic motif art | |
CN103439177B (zh) | 轴心受力状态下水泥净浆早期宏观弹性模量计算方法 | |
Ma et al. | Refined 18-DOF triangular hybrid stress element for couple stress theory | |
Pourbehi et al. | Analysis of combined action of seismic loads and alkali-silica reaction in concrete dams considering the key chemical-physical-mechanical factors and fluid-structure interaction | |
Eghbalian et al. | Micromechanical approach to swelling behavior of capillary‐porous media with coupled physics | |
Pan et al. | Chloride binding to calcium silicate hydrates (CSH) in cement paste: A molecular dynamics analysis | |
Bernard et al. | Multiscale modelling approach to determine the specific heat of cementitious materials | |
Zhang et al. | Parametric variational principle for bi-modulus materials and its application to nacreous bio-composites | |
Cusatis et al. | Multiscale modeling of Alkali Silica reaction degradation of concrete | |
Yang et al. | Evolution of properties in hydration of cements—A numerical study | |
Zhao et al. | Discontinuum based micromechanics modelling methods | |
Dorival et al. | A lattice model to simulate ice-structure interaction | |
Liu et al. | Multiphysics for early stage cement hydration: theoretical framework | |
Schrader et al. | Hybrid computing models for large-scale heterogeneous 3d microstructures | |
Wang et al. | An efficient peridynamic method and its MPI parallelization for simulating the continuous icebreaking process | |
Abali | Asymptotic homogenization and parameter determination in metamaterials | |
Li et al. | Shrinkage and creep caused by dissolution | |
Li et al. | Modeling Creep and Relaxation Caused by Phase Dissolution |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |