CN110750819A - 一种离散元与有限差分法耦合的沥青混合料模拟建模方法 - Google Patents

一种离散元与有限差分法耦合的沥青混合料模拟建模方法 Download PDF

Info

Publication number
CN110750819A
CN110750819A CN201910884357.3A CN201910884357A CN110750819A CN 110750819 A CN110750819 A CN 110750819A CN 201910884357 A CN201910884357 A CN 201910884357A CN 110750819 A CN110750819 A CN 110750819A
Authority
CN
China
Prior art keywords
asphalt mixture
finite difference
asphalt
aggregate
difference method
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
Application number
CN201910884357.3A
Other languages
English (en)
Inventor
于华南
李顺俊
钱国平
龚湘兵
祝轩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changsha University of Science and Technology
Original Assignee
Changsha University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Changsha University of Science and Technology filed Critical Changsha University of Science and Technology
Priority to CN201910884357.3A priority Critical patent/CN110750819A/zh
Publication of CN110750819A publication Critical patent/CN110750819A/zh
Priority to US16/862,878 priority patent/US20210090326A1/en
Priority to BE20205349A priority patent/BE1027548B1/fr
Priority to ZA2020/02980A priority patent/ZA202002980B/en
Pending legal-status Critical Current

Links

Images

Classifications

    • 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
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/26—Composites
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2200/00—Indexing scheme for image data processing or generation, in general
    • G06T2200/04—Indexing scheme for image data processing or generation, in general involving 3D image data

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种离散元与有限差分法耦合的沥青混合料模拟建模方法，通过沥青混合料骨料离散元处理，采用有限差分方法实现沥青混合料中的沥青胶结料连续介质模拟；通过考虑离散元和有限差分法耦合作用对沥青混合料强度、模量等力学性质的影响，以实现沥青混合料多相材料在不同荷载作用下变形、收缩和开裂等过程的模拟。本发明可以精确还原沥青混合料中集料和胶结料的真实形状和大小，以及空隙结构，并表征其分布特征，克服了离散元和有限差分法在沥青混合料数值模拟上各自分析方法的局限性，适用于分析沥青混合料多相材料不同条件下黏塑性大变形问题。

Description

一种离散元与有限差分法耦合的沥青混合料模拟建模方法
技术领域
本发明涉及沥青混合料数值模拟，尤其是三维离散元和有限差分法耦合建模和模拟方法。
背景技术
目前世界上新建道路绝大部分是沥青混凝土路面。而沥青混合料骨料形状的不均匀性和离散性，以及沥青胶浆材料的连续性，使得在沥青混合料的数值模拟中对骨料和沥青胶浆的模拟非常困难。之前研究大多采用离散元对沥青混合料进行模拟计算，虽然离散元在细致考虑集料的接触网络和空隙分布等方面具有优越性，但在考虑沥青胶结料和动态大变形的时候仍存在较大难度。而有限差分法相对于与有限元等其他方法，最大的优点是其处理塑性变形分析能力较强，且可以灵活处理任意本构模型，无需将单元应力带到屈服界面，对于模拟沥青混合料高温塑性流动问题更加高效和准确；然而有限差分法只能处理沥青混合料宏观层面的应力应变和位移现象，无法从细观上分析空隙特性与集料之间的接触互锁。
发明内容
为准确进行沥青混合料三维数值模拟，本发明提出一种离散元和有限差分法耦合条件下的三维建模模拟方法。通过沥青混合料骨料离散元处理，采用有限差分方法实现沥青混合料中的沥青胶结料连续介质模拟。通过考虑离散元和有限差分法耦合作用对沥青混合料强度、模量等力学性质的影响，以实现沥青混合料多相材料在不同荷载作用下变形、收缩和开裂等过程的模拟。
本发明为解决沥青混合料骨料和胶结料模拟问题提出的技术方案是，一种离散元与有限差分法耦合的沥青混合料模拟方法，其步骤包括：(1)通过工业CT对沥青混合料试件进行扫描及后处理，获得沥青混合料骨料、胶结料和空隙结构的三维坐标值；(2)构建沥青混合料的三维模型，赋予混合料骨料形状的三维离散单元属性，采用有限差分法实现沥青砂浆的连续性模拟；(3)考虑骨料形状对沥青混合料力学性质的影响，通过三维离散元建立沥青混合料骨料的细观特征和应力场的耦合；(4)考虑沥青砂浆对沥青混合料力学性质的影响，通过有限差分法建立沥青混合料沥青砂浆特征和应力场的耦合；(5)设定模型荷载位置、荷载大小、荷载方向以及作用时间；(6)利用离散元与有限差分法耦合计算沥青混合料的三维离散和连续单元的变形与破坏，实现沥青混合料骨料运动与迁移、沥青胶浆开裂变形过程的数值模拟与建模。
有益效果：与现有技术相比，本发明的一种离散元与有限差分法耦合的沥青混合料模拟建模方法的优点在于：(1)可以精确还原沥青混合料中集料和胶结料的真实形状和大小，以及空隙结构，并表征其分布特征；(2)克服了离散元和有限差分法在沥青混合料数值模拟上各自分析方法的局限性；(3)可同时考虑集料的真实形状、分析和处理细观现象，并处理压实过程中宏观动态大变形和连续介质的边界条件；(4)本发明所提出的方法精确性好。
附图说明
图1为离散元与有限差分法耦合的沥青混合料模拟建模方法的实现流程图。
图2为耦合计算信息交互边界的建立过程流程图。
图3为骨料与沥青胶结料耦合计算实现流程图。
具体实施方式
本发明方法工作流程如图1-图3所示。
图1为用于沥青混合料的三维离散元与有限差分耦合模拟建模方法实现流程图。
(1)采用工业CT扫描对沥青混合料试件进行扫描，得到沥青混合料CT断层图片；
(2)对沥青混合料CT断层图片进行处理和三维重构，获得沥青混合料的骨料、胶结料和空隙结构中像素的三维坐标值；
(3)构建沥青混合料的三维模型；在三维离散元中建立骨料非连续模型，并在有限差分法中建立沥青胶结料连续模型；
(4)输入计算模型的骨料参数、沥青胶结料参数、位移边界、模型约束条件、荷载条件；
(5)确定骨料和沥青胶浆接触边界，并将其设置为离散元与有限差分法耦合计算的信息交互边界；
(6)利用离散元与有限差分法耦合计算沥青混合料的三维离散和连续单元的变形与破坏，实现沥青混合料骨料运动与迁移、沥青胶浆开裂变形过程的数值模拟；
(7)输出数值模拟结果。
图2为耦合计算信息交互边界的建立过程流程图。步骤(21)-(23)详细说明步骤(5)中信息交互边界的建立过程。
(21)遍历全模型计算域内的离散骨料和沥青胶结料；
(22)确定离散骨料和沥青胶结料的接触边界；
(23)根据骨料空间网格对接触边界上的有限差分单元接触面的空间网格进行网格划分，并将其设置为离散元与有限差分法耦合计算的信息交互边界，以便骨料颗粒与有限差分单元进行信息交互。
图3为骨料与沥青胶结料耦合计算实现流程图。
步骤(31)-(36)详细说明步骤(6)中耦合计算的模拟过程。
(31)根据步骤(5)确定的信息交互边界，建立三维离散元计算软件和有限差分法计算软件的数据通信；
(32)在有限差分法软件中打开大应变模式，以适应沥青胶结料的动态大变形；
(33)在有限差分法中使用计算方程式计算一个循环，并将边界节点速度及其更新后的位置坐标一起写入阵列，再这些数据通过数据接口连接发送到离散元模型；
(34)离散元接收到边界节点速度及其更新后的位置坐标后，根据等效力系统重新计算合力和力矩，再将这些数据反馈给有限差分法模型；
(35)在每个时间步中，通过裂隙发展情况或是否满足设定迭代条件判断是否继续迭代，如果需要迭代，进入步骤(32)，否则，输出迭代结果；
(36)返回模拟结果。

Claims (4)

1.一种离散元和有限差分法耦合的沥青混合料模拟建模方法，其特征是，包括以下步骤：
(1)通过工业CT对沥青混合料试件进行扫描及后处理，获得沥青混合料骨料、胶结料和空隙结构的三维坐标值；
(2)构建沥青混合料的三维模型，赋予混合料骨料形状的三维离散单元属性，采用有限差分法实现沥青砂浆的连续性模拟；
(3)考虑骨料形状对沥青混合料力学性质的影响，通过三维离散元建立沥青混合料骨料的细观特征和应力场的耦合；
(4)考虑沥青砂浆对沥青混合料力学性质的影响，通过有限差分法建立沥青混合料沥青砂浆特征和应力场的耦合；
(5)设定计算模型的骨料参数、沥青胶结料参数、位移边界、模型约束条件、荷载条件；
(6)利用离散元与有限差分法耦合计算沥青混合料的三维离散和连续单元的变形与破坏，实现沥青混合料骨料运动与迁移、沥青胶浆开裂变形过程的数值模拟与建模。
2.根据权利要求1所述离散元和有限差分法耦合的沥青混合料模拟建模方法，其特征是，具体步骤为：
(1)采用工业CT扫描对沥青混合料试件进行扫描，获得沥青混合料CT断层图片；
(2)对沥青混合料CT断层图片进行处理和三维重构，得到断层图片中沥青混合料骨料、胶结料和空隙结构的三维坐标值；
(3)构建沥青混合料的三维模型；在三维离散元中建立骨料非连续模型，并在有限差分法中建立沥青胶结料连续模型；
(4)输入计算模型的骨料参数、沥青胶结料参数、位移边界、模型约束条件、荷载条件；
(5)确定骨料和沥青胶浆接触边界，并将其设置为离散元与有限差分法耦合计算的信息交互边界；
(6)利用离散元与有限差分法耦合计算沥青混合料的三维离散和连续单元的变形与破坏，实现沥青混合料骨料运动与迁移、沥青胶浆开裂变形过程的数值模拟；
(7)输出数值模拟结果。
3.根据权利要求2所述离散元和有限差分法耦合的沥青混合料模拟建模方法，其特征是，步骤(21)-(23)详细说明步骤(5)中信息交互边界的建立过程；
(21)遍历全模型计算域内的离散骨料和沥青胶结料；
(22)确定离散骨料和沥青胶结料的接触边界；
(23)根据骨料空间网格对接触边界上的有限差分单元接触面的空间网格进行网格划分，并将其设置为离散元与有限差分法耦合计算的信息交互边界，以便骨料颗粒与有限差分单元进行信息交互。
4.根据权利要求2所述离散元和有限差分法耦合的沥青混合料模拟建模方法，其特征是，步骤(31)-(36)详细说明步骤(6)中耦合计算的模拟过程；
(31)根据步骤(5)确定的信息交互边界，建立三维离散元计算软件和有限差分法计算软件的数据通信；
(32)在有限差分法软件中打开大应变模式，以适应沥青胶结料的动态大变形；
(33)在有限差分法中使用计算方程式计算一个循环，并将边界节点速度及其更新后的位置坐标一起写入阵列，并将这些数据通过数据接口连接发送到离散元模型；
(34)离散元接收到边界节点速度和及其更新后的位置坐标后，根据等效力系统重新计算合力和力矩，再将这些数据反馈给有限差分模型；
(35)在每个时间步中，通过裂隙发展情况或是否满足设定迭代条件判断是否继续迭代，如果需要迭代，进入步骤(32)，否则，输出迭代结果；
(36)返回模拟结果。
CN201910884357.3A 2019-09-19 2019-09-19 一种离散元与有限差分法耦合的沥青混合料模拟建模方法 Pending CN110750819A (zh)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN201910884357.3A CN110750819A (zh) 2019-09-19 2019-09-19 一种离散元与有限差分法耦合的沥青混合料模拟建模方法
US16/862,878 US20210090326A1 (en) 2019-09-19 2020-04-30 Modeling method for asphalt mixture by coupling discrete element method and finite difference method
BE20205349A BE1027548B1 (fr) 2019-09-19 2020-05-19 Procede de simulation numerique hybride multidimensionnel pour un melange d'asphalte
ZA2020/02980A ZA202002980B (en) 2019-09-19 2020-05-21 Numerical simulation method for asphalt mixture based on discrete element method-finite difference method coupling

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910884357.3A CN110750819A (zh) 2019-09-19 2019-09-19 一种离散元与有限差分法耦合的沥青混合料模拟建模方法

Publications (1)

Publication Number Publication Date
CN110750819A true CN110750819A (zh) 2020-02-04

Family

ID=69276760

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910884357.3A Pending CN110750819A (zh) 2019-09-19 2019-09-19 一种离散元与有限差分法耦合的沥青混合料模拟建模方法

Country Status (4)

Country Link
US (1) US20210090326A1 (zh)
CN (1) CN110750819A (zh)
BE (1) BE1027548B1 (zh)
ZA (1) ZA202002980B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112014216A (zh) * 2020-09-04 2020-12-01 长沙理工大学 一种沥青混合料三维应变破坏准则模型的计算方法
CN112084683A (zh) * 2020-07-22 2020-12-15 国网江西省电力有限公司电力科学研究院 一种输电线路用铝合金材料金具的塑性变形分析方法
CN112417707A (zh) * 2020-12-10 2021-02-26 哈尔滨工业大学 一种构建沥青混合料三维空隙模型的方法

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113434933B (zh) * 2021-06-25 2023-05-16 长安大学 一种微胶囊在沥青混合料拌和过程中破裂的仿真分析方法
CN113927738B (zh) * 2021-10-15 2023-02-17 中国重汽集团青岛重工有限公司 一种搅拌罐螺旋叶片的优化设计方法及成型方法
KR102686034B1 (ko) 2022-01-10 2024-07-17 한국원자력연구원 삼차원 불연속체 기반 수리역학 복합거동 해석 방법 및 시스템
CN114861485B (zh) * 2022-04-14 2024-08-23 大连理工大学 基于形函数插值的沥青混合料三维数值模型跨尺度分析法
CN114757948B (zh) * 2022-06-14 2022-09-06 福建南方路面机械股份有限公司 一种基于深度学习的再生骨料砂浆含量检测方法及装置
CN115060749B (zh) * 2022-06-21 2024-08-27 西安建筑科技大学 一种考虑胶结料随机非均匀性的沥青混合料开裂行为模拟方法
CN116305451B (zh) * 2023-03-02 2024-01-09 中国地质大学（北京） 连续-非连续地质模型建立方法及装置
CN117454698B (zh) * 2023-10-27 2024-05-07 华中科技大学 一种混凝土中纤维拉拔性能的数值模拟方法
CN118016213B (zh) * 2024-02-02 2024-07-16 哈尔滨工业大学 一种基于真实质量分布的沥青混合料离散元模型构建方法
CN118553362B (zh) * 2024-07-30 2024-09-24 内蒙古工业大学 沥青混合料的低温开裂数据的预测方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106446402A (zh) * 2016-09-22 2017-02-22 南京大学（苏州）高新技术研究院 一种土体失水开裂多场耦合离散元快速模拟建模方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106446402A (zh) * 2016-09-22 2017-02-22 南京大学（苏州）高新技术研究院 一种土体失水开裂多场耦合离散元快速模拟建模方法

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
LI MINGGUANG ET AL.: "A multiscale coupling approach between discrete element method and finite difference method for dynamic analysis", 《INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING》 *
ZHANPING YOU ET AL: "Three-Dimensional Discrete Element Models for Asphalt Mixtures", 《JOURNAL OF ENGINEERING MECHANICS》 *
万蕾等: "沥青混合料三维离散元几何模型重构研究", 《低温建筑技术》 *
张楚汉: "论岩石、混凝土离散-接触-断裂分析", 《岩石力学与工程学报》 *
张泉等: "含碎石层沥青加铺结构外部耦合计算方法", 《岩土力学》 *
张泉等: "碎石化混凝土路面沥青加铺层结构的数值分析", 《公路交通科技》 *
陈俊等: "基于三维离散元法的沥青混合料断裂过程模拟", 《华南理工大学学报(自然科学版)》 *
高亮等: "离散元-有限差分耦合法在铁路有砟道床研究中的应用", 《北京交通大学学报》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112084683A (zh) * 2020-07-22 2020-12-15 国网江西省电力有限公司电力科学研究院 一种输电线路用铝合金材料金具的塑性变形分析方法
CN112084683B (zh) * 2020-07-22 2022-07-08 国网江西省电力有限公司电力科学研究院 一种输电线路用铝合金材料金具的塑性变形分析方法
CN112014216A (zh) * 2020-09-04 2020-12-01 长沙理工大学 一种沥青混合料三维应变破坏准则模型的计算方法
CN112014216B (zh) * 2020-09-04 2024-04-05 长沙理工大学 一种沥青混合料三维应变破坏准则模型的计算方法
CN112417707A (zh) * 2020-12-10 2021-02-26 哈尔滨工业大学 一种构建沥青混合料三维空隙模型的方法
CN112417707B (zh) * 2020-12-10 2022-09-09 哈尔滨工业大学 一种构建沥青混合料三维空隙模型的方法

Also Published As

Publication number Publication date
BE1027548A1 (fr) 2021-03-26
US20210090326A1 (en) 2021-03-25
ZA202002980B (en) 2021-05-26
BE1027548B1 (fr) 2021-05-05

Similar Documents

Publication Publication Date Title
CN110750819A (zh) 一种离散元与有限差分法耦合的沥青混合料模拟建模方法
Meng et al. A numerical homogenization study of the elastic property of a soil-rock mixture using random mesostructure generation
Mazzucco et al. Elastoplastic-damaged meso-scale modelling of concrete with recycled aggregates
Shuguang et al. Method of meshing ITZ structure in 3D meso-level finite element analysis for concrete
Asahina et al. Elastically-homogeneous lattice models of damage in geomaterials
Shieh-Beygi et al. Numerical analysis of structural masonry: mesoscale approach
CN209368648U (zh) 水泥稳定rap混合料基层沥青路面反射裂缝试验及评价装置
CN102592029B (zh) 湿陷性黄土路基工后沉降的分析预测方法
Chen et al. Centrifuge modeling testing and multiscale analysis of cemented sand and gravel (CSG) dams
CN114969951B (zh) 加筋土结构的数值计算方法、装置、存储介质及电子设备
CN114186434B (zh) 基于逐级划分的沥青混合料三维细观结构模型构建方法
Shen et al. A mesoscale bond model for discrete element modeling of irregular cemented granular materials
Li et al. A coupled MPM-DEM method for modelling soil-rock mixtures
Meng et al. Rigid-block DEM modeling of mesoscale fracture behavior of concrete with random aggregates
Tsitova et al. Identification of microscale fracture models for mortar with in-situ tests
Mohajerani et al. A novel peridynamics modelling of cemented granular materials
Shi et al. Deformation coordination analysis of CC-RCC combined dam structures under dynamic loads
CN112131805A (zh) 基于离散元方法的双层多孔沥青路面排水特性分析方法
Ling et al. Predicting the mechanical properties of semi-flexible pavement material with micromechanical modeling
Ghadimi et al. Effects of geometrical parameters on numerical modeling of pavement granular material
CN108460838A (zh) 三维可视化技术与数值模拟技术融合的实现方法与系统
Gherbi et al. Study on two way reinforced concrete slab using ANSYS with different boundary conditions and loading
Qian et al. Aggregate fatigue failure on macro texture polishing of asphalt pavement
CN113408171B (zh) 超高性能混凝土的力学性能预测方法
CN109596423A (zh) 一种基于单轴数值试验的沥青混合料性能测试方法及系统

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