CN108334676A - 一种基于python再生混凝土三维随机球形骨料模型的构建方法 - Google Patents
一种基于python再生混凝土三维随机球形骨料模型的构建方法 Download PDFInfo
- Publication number
- CN108334676A CN108334676A CN201810051411.1A CN201810051411A CN108334676A CN 108334676 A CN108334676 A CN 108334676A CN 201810051411 A CN201810051411 A CN 201810051411A CN 108334676 A CN108334676 A CN 108334676A
- Authority
- CN
- China
- Prior art keywords
- aggregate
- random
- model
- parameter
- size
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于python再生混凝土三维随机球形骨料模型的构建方法。定义数组,用于存储直径、坐标、体积参数;输入基本参数,包括骨料含量、模型尺寸、骨料细观尺寸;内置粒径分布函数,内置随机函数;从粒径d1至d2中随机生成一个浮点数,获得所需范围内的随机数;将骨料粒径进行分级;将投放空间划分为均匀的小格子,按照骨料粒径从大到小的方式进行投放,投放位置随机生成,通过循环投放得到满足要求的骨料模型;对所得球形骨料进行面平移,得到外界面区、老硬化砂浆、内界面区;该模型选用ABAQUS软件自带的损伤塑性模型,各参数值存储于mat数组中,保存模型,为再生混凝土的数值模拟研究提供了一种更加快捷高效的方法。
Description
技术领域
本发明涉及建筑技术领域,尤其涉一种基于python再生混凝土三维随机球形骨料模型的构建方法。
背景技术
再生混凝土是指利用再生粗骨料部分或者全部代替天然骨料配置而成的混凝土,再生混凝土技术的开发和利用是建筑废弃物资源化利用的关键。由于再生骨料细观结构的复杂性和试验条件的差异性,单纯依赖试验的研究方法,很难对再生混凝土的研究有进一步的突破。选取合理计算模型进行数值模拟,能够避开试验条件的限制和人为因素的影响,从而实现将再生混凝土的细观结构和宏观力学性能联系起来。
目前,国内外诸多学者在再生混凝土数值模拟方面进行了研究。肖建庄,李文贵等通过对模型再生混凝土的数值模拟,研究了再生混凝土内部的应力分布特征和各项介质对再生混凝土试件力学性能的影响。党娜娜,彭一江等基于蒙特卡罗法建立二维再生混凝土随机骨料模型,研究了再生混凝土的细观损伤机理。许岳周提出颗粒弹簧模型,研究了不同级配下再生混凝土的抗压强度。现阶段的研究取得了一定的成果,但还存在诸多不足之处。首先,固定颗粒数目的计算模型不能体现出骨料的随机性;其次,再生骨料细观结构复杂,在建立细观模型时未能考虑其多界面的特性;最后,大多数的再生混凝土骨料模型是建立在Matlab,Fortran,C++等语言的基础上的,由于这些语言不具备定义材料属性、划分网格等功能,需要建立好骨料模型后导入有限元软件进行后续步骤,增加了许多工作量。
发明内容
本发明的目的在于提供一种基于python再生混凝土三维随机球形骨料模型的构建方法,以解决上述技术问题。
为实现上述目的本发明采用以下技术方案:
一种基于python再生混凝土三维随机球形骨料模型的构建方法,包括如下步骤:
1)、定义数组,用于存储直径、坐标、体积参数;输入基本参数,包括骨料含量、模型尺寸、骨料细观尺寸;
2)、内置粒径分布函数,内置随机函数;从粒径d1至d2中随机生成一个浮点数,获得所需范围内的随机数;
3)、将骨料粒径进行分级,当每级累积骨料体积总和小于该级骨料体积总和时,调取随机粒径函数,继续增加骨料,对每级骨料分别进行计算,获得各级骨料粒径值,并将值计入粒径库;
4)将投放空间划分为均匀的小格子,为了提高投放效率,把骨料位置投放在格子的中心上,按照骨料粒径从大到小的方式进行投放,投放位置随机生成,当骨料投放位置与已投放骨料位置重合时,重新进行投放,相邻骨料间需满足最小容差,通过循环投放得到满足要求的骨料模型;
5)、对所得球形骨料进行面平移,得到外界面区、老硬化砂浆、内界面区;通过面平移实现分区,得到外界面区,通过面平移实现分区,得到老硬化砂浆,通过面平移实现分区,得到内界面区;
6)、编写定义材料属性、创建分析步、划分网格、定义单元类型、定义约束、输出脚本程序,该模型选用ABAQUS软件自带的损伤塑性模型,需定义塑性参数和损伤参数,塑性参数包括屈服函数、流动法则、粘性相关参数,损伤参数包括压缩损伤参数和拉伸损伤参数,各参数值存储于mat数组中,保存模型。
与现有技术相比,本发明具有以下优点:本发明建立的再生骨料计算模型能满足骨料生成随机性要求;本发明构建的含双界面、五相材料组成的再生骨料细观计算模型与真实骨料细观结构更为符合。
本发明将再生骨料模型的生成和数值模拟模块整合在一起只需修改相关参数就可进行相应计算,为再生混凝土的数值模拟研究提供了一种更加快捷高效的方法。
附图说明
图1为本发明再生骨料投放流程方法步骤图。
图2是本发明不同骨料含量的计算模型图。
图3取骨料含量为50%的计算模型进行计算得到应力云图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步详细阐述。
一种基于python再生混凝土三维随机球形骨料模型的构建方法,包括如下步骤:
1)、定义数组,用于存储直径、坐标、体积参数;输入基本参数,包括骨料含量、模型尺寸、骨料细观尺寸;
2)、内置粒径分布函数,内置随机函数;从粒径d1至d2中随机生成一个浮点数,获得所需范围内的随机数;
3)、将骨料粒径进行分级,当每级累积骨料体积总和小于该级骨料体积总和时,调取随机粒径函数,继续增加骨料,对每级骨料分别进行计算,获得各级骨料粒径值,并将值计入粒径库;
4)将投放空间划分为均匀的小格子,为了提高投放效率,把骨料位置投放在格子的中心上,按照骨料粒径从大到小的方式进行投放,投放位置随机生成,当骨料投放位置与已投放骨料位置重合时,重新进行投放,相邻骨料间需满足最小容差,通过循环投放得到满足要求的骨料模型,如图2所示;
5)、对所得球形骨料进行面平移,得到外界面区、老硬化砂浆、内界面区;通过面平移实现分区,得到外界面区,通过面平移实现分区,得到老硬化砂浆,通过面平移实现分区,得到内界面区;
6)、编写定义材料属性、创建分析步、划分网格、定义单元类型、定义约束、输出脚本程序,该模型选用ABAQUS软件自带的损伤塑性模型,需定义塑性参数和损伤参数,塑性参数包括屈服函数、流动法则、粘性相关参数,损伤参数包括压缩损伤参数和拉伸损伤参数,各参数值存储于mat数组中,保存模型。
例如再生骨料最大粒径为20mm,最小粒径5mm。将骨料粒径分为三级,分别为20-15mm、15-10mm、10-5mm。投放区域为100mm×100mm×100mm 的试件空间,骨料含量分别为45%、50%、55%,计算模型如图3所示,其中 (a)为骨料含量为45%的计算模型,(b)为骨料含量为50%的计算模型, (c)为骨料含量为55%的计算模型,(d)为骨料含量为60%的计算模型。
选取骨料含量为50%的计算模型进行计算得到应力云图,图3所示。(a) 为s11应力云图,(b)为s22应力云图,(c)为s33应力云图。
本发明的第一个目的是解决再生骨料计算模型中骨料生成的随机性问题。
实现本发明第一个目的得技术方案是:应用蒙特卡洛法使得再生骨料随机产生。通过构造随机模型使得某一随机变量的数学期望等于所求的解,最基本的随机变量在区间[0,1]上服从均匀分布。若产生[0,1]的均匀随机变量x,则经变换X=a+(b-a)x后,就可求得在区间[a,b]上的均匀分布的随机变量X。对于骨料颗粒的随机位置坐标(Xi,Yi,Zi)、骨料随机半径Ri等参数,均可用均匀随机变量经过变换得到。
本发明的第二个目的是构建出符合真实骨料细观结构的再生骨料计算模型。
实现本发明第二个目的的技术方案是:通过查阅再生骨料细观结构方面的试验资料,建立含双界面、五相材料组成的的三维球形再生骨料计算模型。再生粗骨料来源于原生混凝土,本身性能就有很多差异,其细观结构较普通混凝土具有更高的复杂性和随机性。由于再生骨料上附着砂浆(即老硬化砂浆)的存在,使得再生混凝土内部存在较多的过渡界面区。在细观层次上,将再生混凝土看作由新硬化砂浆、外界面区、老硬化砂浆、内界面区、天然骨料五相材料组成。结合相关研究成果,将再生骨料简化为球形,内界面区、老硬化砂浆、外界面区均以球状包裹于天然骨料外侧。老硬化砂浆厚度取为 1mm,内外界面区厚度取为0.5mm。
本发明的第三个目的是提高模型生成效率,并将再生骨料模型的生成和数值模拟模块整合在一起,简化计算过程。
实现本发明第三个目的的技术方案是:Python是一种面向对象并具有动态语义特征的编程语言,具有高级的数据类型和简单有效的程序设计方法。 Python具有面向对象、可移植性、可拓展性、免费开源等特点,所以特别适合作为脚本嵌入于应用程序中。ABAQUS为二次开发用户提供了Python语言接口,通过集成Python脚本向二次开发用户提供丰富的库函数。可直接操纵ABAQUS的内核,实现建模、定义材料属性、划分网格、后处理等分析功能。
以上所述为本发明较佳实施例,对于本领域的普通技术人员而言,根据本发明的教导,在不脱离本发明的原理与精神的情况下,对实施方式所进行的改变、修改、替换和变型仍落入本发明的保护范围之内。
Claims (1)
1.一种基于python再生混凝土三维随机球形骨料模型的构建方法,其特征在于,包括如下步骤:
1)、定义数组,用于存储直径、坐标、体积参数;输入基本参数,包括骨料含量、模型尺寸、骨料细观尺寸;
2)、内置粒径分布函数,内置随机函数;从粒径d1至d2中随机生成一个浮点数,获得所需范围内的随机数;
3)、将骨料粒径进行分级,当每级累积骨料体积总和小于该级骨料体积总和时,调取随机粒径函数,继续增加骨料,对每级骨料分别进行计算,获得各级骨料粒径值,并将值计入粒径库;
4)将投放空间划分为均匀的小格子,为了提高投放效率,把骨料位置投放在格子的中心上,按照骨料粒径从大到小的方式进行投放,投放位置随机生成,当骨料投放位置与已投放骨料位置重合时,重新进行投放,相邻骨料间需满足最小容差,通过循环投放得到满足要求的骨料模型;
5)、对所得球形骨料进行面平移,得到外界面区、老硬化砂浆、内界面区;通过面平移实现分区,得到外界面区,通过面平移实现分区,得到老硬化砂浆,通过面平移实现分区,得到内界面区;
6)、该模型选用ABAQUS软件自带的损伤塑性模型,需定义塑性参数和损伤参数,塑性参数包括屈服函数、流动法则、粘性相关参数,损伤参数包括压缩损伤参数和拉伸损伤参数,各参数值存储于mat数组中,保存模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810051411.1A CN108334676A (zh) | 2018-01-19 | 2018-01-19 | 一种基于python再生混凝土三维随机球形骨料模型的构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810051411.1A CN108334676A (zh) | 2018-01-19 | 2018-01-19 | 一种基于python再生混凝土三维随机球形骨料模型的构建方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108334676A true CN108334676A (zh) | 2018-07-27 |
Family
ID=62926354
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810051411.1A Pending CN108334676A (zh) | 2018-01-19 | 2018-01-19 | 一种基于python再生混凝土三维随机球形骨料模型的构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108334676A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109766665A (zh) * | 2019-02-02 | 2019-05-17 | 中交公路长大桥建设国家工程研究中心有限公司 | 一种基于有限元分析的轻骨料混凝土配合比设计方法 |
CN109979009A (zh) * | 2019-03-04 | 2019-07-05 | 天津大学 | 一种具有功能梯度的胶结颗粒材料三维细观数字模型重构方法 |
CN110210178A (zh) * | 2019-06-26 | 2019-09-06 | 西安理工大学 | 一种基于Python再生混凝土三维随机球形骨料模型的构建方法 |
CN110427732A (zh) * | 2019-09-11 | 2019-11-08 | 河海大学 | 一种双层非均质土层结构模型构建方法 |
CN112632814A (zh) * | 2020-12-04 | 2021-04-09 | 西安理工大学 | 一种基于椭球体随机骨料的再生混凝土三维模型构建方法 |
CN113092248A (zh) * | 2021-03-23 | 2021-07-09 | 武汉大学 | 一种基于耦合格构模型的长龄期混凝土开裂模拟方法 |
CN113177315A (zh) * | 2021-04-29 | 2021-07-27 | 东南大学 | 一种研究微生物自修复混凝土中修复剂颗粒分布特征的方法 |
CN113204809A (zh) * | 2021-05-25 | 2021-08-03 | 湖北工业大学 | 一种骨料投放过程中考虑骨料均匀性的随机骨料投放方法 |
CN113268909A (zh) * | 2021-05-26 | 2021-08-17 | 太原理工大学 | 一种随机分布属性的细观混凝土损伤演化分析方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20160060793A (ko) * | 2014-10-21 | 2016-05-31 | 군산대학교산학협력단 | 능동형 정보수집 스크립트 및 수치모형을 이용한 태풍해일 자동 예측 방법 |
US9589038B1 (en) * | 2014-03-28 | 2017-03-07 | Amazon Technologies, Inc. | Attribute tracking, profiling, and recognition |
CN106503405A (zh) * | 2016-12-02 | 2017-03-15 | 东南大学 | 一种基于综合模型的自然火灾下梁柱式木结构框架整体抗火性能的分析系统及分析方法 |
CN106650098A (zh) * | 2016-12-22 | 2017-05-10 | 中铁二院工程集团有限责任公司 | 一种混凝土徐变应变计算方法 |
-
2018
- 2018-01-19 CN CN201810051411.1A patent/CN108334676A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9589038B1 (en) * | 2014-03-28 | 2017-03-07 | Amazon Technologies, Inc. | Attribute tracking, profiling, and recognition |
KR20160060793A (ko) * | 2014-10-21 | 2016-05-31 | 군산대학교산학협력단 | 능동형 정보수집 스크립트 및 수치모형을 이용한 태풍해일 자동 예측 방법 |
CN106503405A (zh) * | 2016-12-02 | 2017-03-15 | 东南大学 | 一种基于综合模型的自然火灾下梁柱式木结构框架整体抗火性能的分析系统及分析方法 |
CN106650098A (zh) * | 2016-12-22 | 2017-05-10 | 中铁二院工程集团有限责任公司 | 一种混凝土徐变应变计算方法 |
Non-Patent Citations (1)
Title |
---|
尤作磊: "基于细观尺度的水泥混凝土受荷开裂模拟", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109766665A (zh) * | 2019-02-02 | 2019-05-17 | 中交公路长大桥建设国家工程研究中心有限公司 | 一种基于有限元分析的轻骨料混凝土配合比设计方法 |
CN109766665B (zh) * | 2019-02-02 | 2023-04-07 | 中交公路长大桥建设国家工程研究中心有限公司 | 一种基于有限元分析的轻骨料混凝土配合比设计方法 |
CN109979009B (zh) * | 2019-03-04 | 2023-03-24 | 天津大学 | 一种具有功能梯度的胶结颗粒材料三维细观数字模型重构方法 |
CN109979009A (zh) * | 2019-03-04 | 2019-07-05 | 天津大学 | 一种具有功能梯度的胶结颗粒材料三维细观数字模型重构方法 |
CN110210178A (zh) * | 2019-06-26 | 2019-09-06 | 西安理工大学 | 一种基于Python再生混凝土三维随机球形骨料模型的构建方法 |
CN110427732A (zh) * | 2019-09-11 | 2019-11-08 | 河海大学 | 一种双层非均质土层结构模型构建方法 |
CN110427732B (zh) * | 2019-09-11 | 2023-04-18 | 河海大学 | 一种双层非均质土层结构模型构建方法 |
CN112632814A (zh) * | 2020-12-04 | 2021-04-09 | 西安理工大学 | 一种基于椭球体随机骨料的再生混凝土三维模型构建方法 |
CN113092248A (zh) * | 2021-03-23 | 2021-07-09 | 武汉大学 | 一种基于耦合格构模型的长龄期混凝土开裂模拟方法 |
CN113177315A (zh) * | 2021-04-29 | 2021-07-27 | 东南大学 | 一种研究微生物自修复混凝土中修复剂颗粒分布特征的方法 |
CN113204809B (zh) * | 2021-05-25 | 2022-07-05 | 湖北工业大学 | 一种骨料投放过程中考虑骨料均匀性的随机骨料投放方法 |
CN113204809A (zh) * | 2021-05-25 | 2021-08-03 | 湖北工业大学 | 一种骨料投放过程中考虑骨料均匀性的随机骨料投放方法 |
CN113268909A (zh) * | 2021-05-26 | 2021-08-17 | 太原理工大学 | 一种随机分布属性的细观混凝土损伤演化分析方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108334676A (zh) | 一种基于python再生混凝土三维随机球形骨料模型的构建方法 | |
CN110210178A (zh) | 一种基于Python再生混凝土三维随机球形骨料模型的构建方法 | |
Ranjbar et al. | Structural-acoustic optimization of a rectangular plate: A tabu search approach | |
CN102628861B (zh) | 一种大体积混凝土温度开裂数值模拟方法 | |
JP6784780B2 (ja) | 大規模再生可能エネルギーのデータについて確率モデルを構築する方法 | |
CN102207987A (zh) | 基于OpenCL的GPU加速三维时域有限差分电磁场仿真的方法 | |
CN104077428A (zh) | 一种服务产业联盟的远程有限元分析方法 | |
Eleftheriadis et al. | Automated specification of steel reinforcement to support the optimisation of RC floors | |
CN103729506B (zh) | 一种复杂模型完全六面体建模及几何体重塑加密方法 | |
CN109446731A (zh) | 一种基于abaqus的岩土工程数值模拟方法 | |
Dunant | Experimental and modelling study of the alkali-silica-reaction in concrete | |
Meng et al. | A discrete numerical method for brittle rocks using mathematical programming | |
CN106502771A (zh) | 基于kernel函数的时间开销模型构建方法及系统 | |
Fu et al. | CPS information security risk evaluation system based on Petri net | |
CN103425898B (zh) | 一种水泥净浆微观数字模型生成算法 | |
Kamyab et al. | Size optimization of nonlinear scallop domes by an enhanced particle swarm algorithm | |
CN108038329A (zh) | 一种微胶囊自愈合材料损伤和愈合机理的数值仿真方法 | |
Yang et al. | Isogeometric double-objective shape optimization of free-form surface structures with Kirchhoff–Love shell theory | |
CN102750428B (zh) | 供水管网节点水力计算方程并行算法分区域构造方法 | |
Badel et al. | A note on integration schemes for the microplane model of the mechanical behaviour of concrete | |
CN114840891A (zh) | 一种任意形状试件的三维混凝土细观模型建模方法 | |
CN110675088B (zh) | 一种数字孪生系统复杂任务高效划分方法 | |
CN106354913A (zh) | 一种用于斜齿轮传动系统参数化设计平台及使用方法 | |
Yang et al. | [Retracted] Stress Analysis of Concrete Materials Based on Finite Element Analysis | |
Krol et al. | Solving PDEs in modern multiphysics simulation software |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180727 |