• 专利附录
  • 专利说明
  • 权利要求
  • 类似技术
  • 同族专利
  • 引用文献
  • 法律状态
  • 优先权
    • 专利摘要:
    The present invention discloses a multidimensional hybrid numerical simulation method for an asphalt mix. Aggregates of an asphalt mix are processed by a discrete element method (DEM), and a finite difference method is used to implement a continuous medium simulation of an asphalt binder in the mix. asphalt. The present invention can accurately restore a cavity structure as well as the actual shapes and sizes of aggregates and a binder in the asphalt mix, and can characterize distribution characteristics thereof.
    公开号:BE1027548B1
    申请号:E20205349
    申请日:2020-05-19
    公开日:2021-05-05
    发明作者:Huanan Yu；Shunjun Li；Guoping Qian；Xiangbing Gong；Xuan Zhu
    申请人:Univ Changsha Sci & Tech；
    IPC主号:
    专利说明:

    TECHNICAL FIELD The present invention relates to a digital simulation of an asphalt mixture and, in particular, to a multidimensional hybrid digital simulation method for an asphalt mixture.
    Background Currently, most new roads in the world use asphalt concrete pavements.
    However, in an asphalt mix, aggregates form inconspicuously with uneven shapes, and an asphalt mortar material exhibits some continuity.
    Therefore, it is very difficult to simulate the aggregate and the asphalt mortar during a digital simulation of the asphalt mix.
    In most previous studies, a Discrete Element Method (DEM) is used to perform a simulated calculation of an asphalt mix.
    A DEM has several advantages for a simulated calculation when a contact network, a distribution of cavities, etc. of aggregates are carefully considered.
    However, when an asphalt binder and large dynamic strain are taken into consideration, it is still difficult to use DEM to perform a simulated calculation.
    Compared to a finite element method (FEM) or other methods, a finite difference method (FDM) has the greatest advantage of a relatively large plastic strain analysis capacity.
    In addition, it can flexibly handle any constitutive model without introducing element stress into a plastic surface.
    Therefore, it is more efficient and accurate to simulate high temperature plastic flow of an asphalt mixture.
    However, FDM can only deal with stress, strain and displacement of the asphalt mix at a macro level, but cannot analyze cavity characteristics and cannot achieve contacting and interlocking between aggregates. at a micro level.
    Summary To accurately perform a three-dimensional numerical simulation of an asphalt mix, the present invention provides a three-dimensional modeling and simulation method for an asphalt mix by combining a discrete element method (DEM) and a difference method. finished (FDM). Aggregates of an asphalt mix are processed by DEM, and FDM is used to implement a continuous medium simulation of an asphalt binder in the asphalt mix.
    A simulation of deformation, shrinkage, cracking, etc. of a multiphase material in the asphalt mix under different load is implemented considering the influence of the combination of DEM and FDM on the mechanical properties of the asphalt mix, such as strength and elasticity module.
    To solve the problem of simulating aggregates and a binder in an asphalt mixture, the present invention provides a technical solution: a multidimensional hybrid numerical simulation method for an asphalt mixture.
    The method comprises the following steps: (1) scanning a test sample of an asphalt mixture by computer-aided industrial tomography (CT), and performing post-processing, in order to obtain three-dimensional coordinates of the aggregates, d a binder, and a cavity structure of the asphalt mixture; (2) construct a three-dimensional model of the asphalt mix, assign a three-dimensional DEM attribute to an aggregate shape of the mix, and perform a PE 2020/5368 asphalt mortar continuity simulation using an FDM; (3) consider the influence of the aggregate shape on the mechanical properties of the asphalt mix, relate a microscopic characteristic of the aggregates of the asphalt mix and a stress field using a three-dimensional DEM; (4) consider the influence of asphalt mortar on the mechanical properties of the asphalt mix, relate a characteristic of the asphalt mortar of the asphalt mix and a stress field using FDM; (5) defining an aggregate parameter, an asphalt binder parameter, a displacement boundary, a model stress condition and a load condition of a computational model; and (6) calculate the deformation and errors of the three-dimensional DEM and a continuous element of the asphalt mix by combining the DEM and FDM, in order to implement numerical simulation and modeling of displacement and migration aggregate as well as cracking and deformation of the asphalt mortar in the asphalt mix.
    Advantageous Effects: Compared with the prior art, the multidimensional hybrid numerical simulation method for an asphalt mix in the present invention has the following advantages: (1) the present invention can accurately restore a cavity structure as well as the shapes and actual sizes of aggregates and a binder in an asphalt mixture, and can characterize the distribution characteristics thereof; (2) the present invention overcomes the limitations of the analysis methods based on DEM and FDM respectively in a digital simulation of the asphalt mix; (3) in the present invention, the actual shape of the aggregates can be taken into consideration, a microscopic phenomenon of the aggregates can be analyzed for further processing, and processing can be performed on the basis of the large macroscopic dynamic deformation and a boundary condition of a continuous medium during compaction; and (4) the method provided in the present invention exhibits high precision. Brief Description of the Drawings Figure 1 is a flowchart for implementing a multidimensional hybrid numerical simulation method for an asphalt mix; Figure 2 is a flowchart of a method for establishing an information exchange boundary for a combined calculation; and Figure 3 is a flowchart for implementing a combined calculation for aggregate and asphalt binder. Detailed Description A work flow diagram of a method according to the present invention is shown in Figures 1 to 3.
    Figure 1 is a flowchart for implementing a multidimensional hybrid numerical simulation method for an asphalt mix.
    (1) perform an industrial CT scan on the test sample of the asphalt mix to obtain an erroneous CT image of the asphalt mix; (2) perform processing and three-dimensional reconstruction of the erroneous CT image of the asphalt mixture to obtain three-dimensional pixel coordinates of the aggregates, binder and cavity structure of the asphalt mixture in the image Erroneous inspection; (3) construct a three-dimensional model of the asphalt mix, establish a discontinuous aggregate model in a three-dimensional DEM, and establish a continuous model of the asphalt binder in an FDM;
    (4) input an aggregate parameter, an asphalt binder PE 2020/5368 parameter, a displacement boundary, a model stress condition and a computational model load condition; (5) determining a contact boundary between the aggregates and the asphalt mortar, and defining the contact boundary as the information exchange boundary for the calculation by combining the DEM and the FDM; (6) calculate the deformation and errors of the three-dimensional DEM and a continuous element of the asphalt mix by combining the DEM and the FDM, in order to implement a numerical simulation of the movement and migration of aggregates as well as cracking and warping of the asphalt mortar in the asphalt mix; and (7) generating a result of the digital simulation.
    Figure 2 is a flowchart of a method for establishing an information exchange boundary for a combined computation. A method for establishing the information exchange boundary in step (5) is described in detail from step (21) to step (23): (21) traversing discrete aggregates and the binder d 'asphalt in a computational domain of the complete model; (22) determining a contact boundary between the discrete aggregates and the asphalt binder; and (23) partition spatial grids of a contact surface of an element of finite differences on the contact boundary according to spatial grids of aggregates, and define partitioned spatial grids as the information exchange boundary for the computation by combining a DEM and an FDM, in order to implement an exchange of information between aggregate particles and the element of finite differences.
    Figure 3 is a flowchart for implementing a combined calculation for aggregates and asphalt binder.
    A calculation simulation method combined with step (6) is described in detail from step (31) to step (36): (31) establishing data communication between three-dimensional DEM calculation software and FDM calculation software according to the information exchange frontier determined in step (5); (32) start a large deformation mode in the software of FDM, in order to adapt it to the large dynamic deformation of the asphalt binder; (33) calculate a cycle in an FDM using a computational equation, transcribe both a velocity of a boundary node and updated position coordinates of it into a network, and send the data to a DEM model through a data interface connection; (34) once a DEM has received the boundary node velocity along with its updated position coordinates, recalculate the resulting force and torque according to an equivalent force system, and then return the data to a finite difference model; (35) in each time step, determining whether to continue an iteration based on a status of crack development or if a specific iteration condition is met; and if iteration is required, returning to step (32); or otherwise, generate a result of the iteration; and (36) returning a simulated result.
    权利要求:
    Claims (3)
    [1]
    A multidimensional hybrid numerical simulation method for an asphalt mixture, comprising the following steps: (1) scanning a test sample of an asphalt mixture by industrial computer-aided tomography (CT), and performing a post -treatment, in order to obtain three-dimensional coordinates of the aggregates, of a binder and of a structure of cavities of the asphalt mixture; (2) construct a three-dimensional model of the asphalt mix, assign a three-dimensional discrete element method (DEM) attribute to an aggregate shape of the mix, and perform a simulation of asphalt mortar continuity using a method of finite differences (FDM); (3) consider the influence of the aggregate shape on the mechanical properties of the asphalt mix, relate a microscopic characteristic of the aggregates of the asphalt mix and a stress field using a three-dimensional DEM; (4) consider the influence of asphalt mortar on the mechanical properties of the asphalt mix, relate a characteristic of the asphalt mortar of the asphalt mix and a stress field using FDM; (5) defining an aggregate parameter, an asphalt binder parameter, a displacement boundary, a model stress condition and a load condition of a computational model; and (6) calculate the strain and errors of the three-dimensional DEM and a continuous element of the asphalt mix by combining the DEM and FDM, in order to implement numerical simulation and modeling of displacement and migration aggregate as well as cracking and deformation of the asphalt mortar in the asphalt mix.
    [2]
    A multidimensional hybrid numerical simulation method for an asphalt mixture according to claim 1, comprising the following specific steps: (1) performing an industrial CT scan on the test sample of the asphalt mixture in order to obtain a test sample. erroneous CT image of asphalt mix; (2) Perform three-dimensional processing and reconstruction of the erroneous CT image of the asphalt mixture in order to obtain three-dimensional coordinates of the aggregate, binder and cavity structure of the asphalt mix in the erroneous CT image ; (3) construct the three-dimensional model of the asphalt mix, establish a discontinuous aggregate model in the three-dimensional DEM, and establish a continuous model of the asphalt binder in the FDM; (4) input the aggregate parameter, the asphalt binder parameter, the displacement boundary, the model stress condition and the load condition of the computational model; (5) determine a {contact boundary between the aggregates and the asphalt mortar, and define the contact boundary as the information exchange boundary for the calculation by combining the DEM and the FDM; (6) calculate the deformation and errors of the three-dimensional DEM and the continuous element of the asphalt mix by combining the DEM and the FDM, in order to implement a numerical simulation of the movement and migration of aggregates as well as cracking and deformation of the asphalt mortar in the asphalt mix; and (7) generating a result of the digital simulation.
    [3]
    The multidimensional hybrid digital simulation method for an asphalt mix according to claim 2, wherein a method for establishing the information exchange boundary in step (5) is described in detail in step (21). ) in step (23): (21) traversing discrete aggregates and the asphalt binder in a computational domain of the complete model; (22) determining a contact boundary between the discrete aggregates and the asphalt binder; and (23) partition spatial grids of a contact surface of an element of finite differences on the contact boundary according to spatial grids of aggregates, and define partitioned spatial grids as the information exchange boundary for computation by combining DEM and FDM, in order to implement an information exchange between aggregate particles and the finite difference element.
    4. A multidimensional hybrid numerical simulation method for an asphalt mix according to claim 2, wherein a computational simulation method combined with step (6) is described in detail from step (31) to step. (36): (31) establishing data communication between three-dimensional DEM calculation software and FDM calculation software according to the information exchange boundary determined in step (5); (32) start a large deformation mode in the software of FDM, in order to adapt it to the large dynamic deformation of the asphalt binder; (33) calculate a cycle in the FDM using a computational equation, transcribe both a velocity of a boundary node and updated position coordinates of it into a network, and send the data to a DEM model through a data interface connection;
    (34) Once the DEM has received the speed of the boundary node along with its updated position coordinates, recalculate the resulting force as well as a torque according to an equivalent force system, and then return the data to a finite difference model;
    (35) in each time step, determining whether to continue iteration based on a crack development status or whether a specific iteration condition is met; and if iteration is required,
    return to step (32); or otherwise, generate a result of the iteration; and
    (36) return a simulated result.
    类似技术:
    公开号 | 公开日 | 专利标题
    Ye et al.2003|Three-dimensional microstructure analysis of numerically simulated cementitious materials
    US7451637B2|2008-11-18|Method of evaluating contact characteristics, and computer product
    Sadowski et al.2015|ANN modeling of pull-off adhesion of concrete layers
    CN107038524B|2020-03-10|Comprehensive evaluation method for construction quality of roller compacted concrete dam considering parameter uncertainty
    Kim et al.2016|Two-dimensional virtual microstructure generation of particle-reinforced composites
    CN108153989B|2020-04-21|Concrete dam compaction quality rapid prediction method considering parameter uncertainty influence
    Baldo et al.2018|Analysis of the mechanical behaviour of asphalt concretes using artificial neural networks
    CN102162577A|2011-08-24|Pipeline defect surface integrity detection device and detection method
    BE1027548B1|2021-05-05|MULTIDIMENSIONAL HYBRID DIGITAL SIMULATION PROCESS FOR AN ASPHALT MIXTURE
    CN109033537B|2020-12-01|Calculation method and system for numerical simulation in rock-fill concrete pouring process
    Majidi et al.2018|On the use of the extended finite element and incremental methods in brittle fracture assessment of key-hole notched polystyrene specimens under mixed mode I/II loading with negative mode I contributions
    CN104021199B|2017-04-19|Function module detecting method based on node domination capacity similarity
    Sadok et al.2016|Genetic programming for granular compactness modelling
    Han et al.2011|Technical comparisons of simulation-based productivity prediction methodologies by means of estimation tools focusing on conventional earthmovings
    Leonardi et al.2018|Road degradation survey through images by drone
    Liu et al.2021|An improved differential box counting method to measure fractal dimensions for pavement surface skid resistance evaluation
    CN104809756B|2018-07-06|Asphalt mixture gap spatial construction method based on X-ray CT images
    Ma et al.2021|Experimental study on the influence of height and dip angle of asperity on the mechanical properties of rock joints
    Broekman et al.2020|Digital twinning of asphalt pavement surfacings using Visual Simultaneous Localization and Mapping
    Hu et al.2015|Slope excavation quality assessment and excavated volume calculation in hydraulic projects based on laser scanning technology
    CN110864610A|2020-03-06|Method for testing three-dimensional roughness of reconstructed concrete surface based on 3D scanning
    Zeghal2008|Modeling the creep compliance of asphalt concrete using the artificial neural network technique
    Dogariu et al.2016|Model preparation for structural FEA on main components of an internal combustion engine
    CN112164133A|2021-01-01|Pavement structure interlayer interface extraction and analysis method based on reverse modeling
    Gong et al.2016|A MATLAB program for Fourier analysis and evaluation on mineral aggregate morphological features
    同族专利:
    公开号 | 公开日
    US20210090326A1|2021-03-25|
    ZA202002980B|2021-05-26|
    CN110750819A|2020-02-04|
    BE1027548A1|2021-03-26|
    引用文献:
    公开号 | 申请日 | 公开日 | 申请人 | 专利标题
    法律状态:
    2021-05-07| FG| Patent granted|Effective date: 20210505 |
    优先权:
    申请号 | 申请日 | 专利标题
    CN201910884357.3A|CN110750819A|2019-09-19|2019-09-19|Discrete element and finite difference method coupled asphalt mixture simulation modeling method|
    [返回顶部]