CN109033537B - Calculation method and system for numerical simulation of rockfill concrete pouring process - Google Patents

Calculation method and system for numerical simulation of rockfill concrete pouring process Download PDF

Info

Publication number
CN109033537B
CN109033537B CN201810699310.5A CN201810699310A CN109033537B CN 109033537 B CN109033537 B CN 109033537B CN 201810699310 A CN201810699310 A CN 201810699310A CN 109033537 B CN109033537 B CN 109033537B
Authority
CN
China
Prior art keywords
stacking
dimensional model
final
initial state
self
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.)
Expired - Fee Related
Application number
CN201810699310.5A
Other languages
Chinese (zh)
Other versions
CN109033537A (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.)
China Agricultural University
Original Assignee
China Agricultural University
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 China Agricultural University filed Critical China Agricultural University
Priority to CN201810699310.5A priority Critical patent/CN109033537B/en
Publication of CN109033537A publication Critical patent/CN109033537A/en
Application granted granted Critical
Publication of CN109033537B publication Critical patent/CN109033537B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Civil Engineering (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Architecture (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The embodiment of the invention provides a calculation method and a system for numerical simulation in a rockfill concrete pouring process, wherein the calculation method comprises the following steps: calculating and analyzing the stacking process and the final stacking form of the three-dimensional model of the stacking stone initial state according to a finite element/discrete element coupling analysis method FEM/DEM; after the stacking process and the final stacking form of the three-dimensional model of the initial state of the stacking blocks are calculated and analyzed, the pouring process and the final filling state of the self-compacting concrete in the final stacking form of the stacking blocks are calculated and analyzed based on a smooth particle fluid dynamics (SPH) method. The calculation method and the system for numerical simulation of the rockfill concrete pouring process, provided by the embodiment of the invention, can accurately simulate the pouring process of the self-compacting concrete, not only can visually observe the flowing state of the self-compacting concrete in the rockfill pores, but also can know the final compacting state of the rockfill concrete, so that a large amount of labor cost, experimental cost and economic cost are saved, and the problem that the experimental method is difficult to observe on line is avoided.

Description

堆石混凝土浇筑过程数值模拟的计算方法和系统Calculation method and system for numerical simulation of rockfill concrete pouring process

技术领域technical field

本发明实施例涉及堆石混凝土浇筑技术领域,更具体地,涉及一种堆石混凝土浇筑过程数值模拟的计算方法和系统。The embodiments of the present invention relate to the technical field of rockfill concrete pouring, and more particularly, to a calculation method and system for numerical simulation of the pouring process of rockfill concrete.

背景技术Background technique

大体积混凝土在现代工程建设,特别是水利水电工程建设中,占有重要地位。我国每年仅在水利水电工程中所浇筑的大体积混凝土就在一千万立方以上,另外,港口、机场建筑以及重型机器基础等也往往采用大体积混凝土。堆石混凝土通过使用大量块石(粒径大于30cm)来降低水泥用量,从而有效减少水化热和降低CO2排放,具有混凝土结构收缩性小、抗裂抗剪能力提高、施工速度快、水化热低,温控容易、施工质量高、工程造价低等特点,因而具有广阔的应用前景,同时也满足我国大力推广绿色低碳技术的迫切需求。Mass concrete occupies an important position in modern engineering construction, especially in water conservancy and hydropower engineering construction. In my country, the mass concrete poured in water conservancy and hydropower projects alone is more than 10 million cubic meters every year. In addition, port, airport buildings and heavy machinery foundations also often use mass concrete. Rockfill concrete reduces the amount of cement by using a large amount of boulders (particle size greater than 30cm), thereby effectively reducing the heat of hydration and CO2 emissions, and has the advantages of small shrinkage of concrete structure, improved crack resistance and shear resistance, fast construction speed, water It has the characteristics of low chemical heat, easy temperature control, high construction quality and low project cost, so it has broad application prospects, and also meets the urgent needs of China to vigorously promote green and low-carbon technologies.

堆石混凝土在施工过程中,堆石体的堆积过程与最后的堆积形态,对成型后的堆石混凝土的力学性能有着较大的影响,但由于实验造价高且实际工程中堆石体的形态样式复杂,因此很难开展较为精确的实验对实际工程进行分析,因此开展数值计算方法对堆石混凝土堆石体堆积过程进行分析计算是必然的趋势。但到目前为止,如何精确的模拟计算堆石混凝土的浇筑施工过程,尚未有任何可供借鉴的成果和方法。During the construction of rockfill concrete, the accumulation process and final accumulation form of rockfill body have a great influence on the mechanical properties of rockfill concrete after forming. However, due to the high experimental cost and the shape of rockfill body in actual engineering The style is complex, so it is difficult to carry out more accurate experiments to analyze the actual engineering. Therefore, it is an inevitable trend to carry out numerical calculation methods to analyze and calculate the accumulation process of rockfill concrete rockfill. But so far, how to accurately simulate and calculate the pouring construction process of rockfill concrete, there are no results and methods that can be used for reference.

发明内容SUMMARY OF THE INVENTION

针对现有技术存在的问题,本发明实施例提供一种堆石混凝土浇筑过程数值模拟的计算方法和系统。In view of the problems existing in the prior art, the embodiments of the present invention provide a calculation method and system for numerical simulation of the pouring process of rockfill concrete.

本发明实施例提供一种堆石混凝土浇筑过程数值模拟的计算方法,包括:根据有限元/离散元耦合分析方法FEM/DEM,计算分析堆石块初始状态三维模型的堆积过程及最终堆积形态;在计算分析所述堆石块初始状态三维模型的堆积过程及最终堆积形态后,基于光滑粒子流体动力学SPH方法,计算分析自密实混凝土在堆石块的最终堆积形态中的浇筑过程及最终填充状态。The embodiment of the present invention provides a calculation method for numerical simulation of the pouring process of rockfill concrete, including: calculating and analyzing the stacking process and final stacking shape of the three-dimensional model of the initial state of the rockfill block according to the finite element/discrete element coupling analysis method FEM/DEM; After calculating and analyzing the stacking process and final stacking shape of the 3D model of the initial state of the rockfill, based on the smooth particle hydrodynamic SPH method, the pouring process and final filling of the self-compacting concrete in the final stacking shape of the rockfill are calculated and analyzed. state.

本发明实施例提供一种堆石混凝土浇筑过程数值模拟的计算系统,包括:堆积分析模块,用于根据有限元/离散元耦合分析方法FEM/DEM,计算分析堆石块初始状态三维模型的堆积过程及最终堆积形态;浇筑分析模块,用于在计算分析所述堆石块初始状态三维模型的堆积过程及最终堆积形态后,基于光滑粒子流体动力学SPH方法,计算分析自密实混凝土在堆石块的最终堆积形态中的浇筑过程及最终填充状态。The embodiment of the present invention provides a computing system for numerical simulation of rockfill concrete pouring process, including: a stacking analysis module, used for calculating and analyzing the stacking of the three-dimensional model of the initial state of rockfill blocks according to the finite element/discrete element coupling analysis method FEM/DEM Process and final stacking form; the pouring analysis module is used for calculating and analyzing the stacking process and final stacking form of the 3D model of the initial state of the rockfill block, based on the smooth particle hydrodynamic SPH method, to calculate and analyze the self-compacting concrete in the rockfill. The pouring process and final filling state in the final packed form of the block.

本发明实施例提供一种堆石混凝土浇筑过程数值模拟的计算设备,包括:至少一个处理器;以及与所述处理器通信连接的至少一个存储器,其中:所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程序指令能够执行上述计算方法。An embodiment of the present invention provides a computing device for numerical simulation of a pouring process of rockfill concrete, including: at least one processor; and at least one memory connected in communication with the processor, wherein: the memory stores data that can be processed by the processor. The program instructions are executed by the processor, and the processor can execute the above calculation method by calling the program instructions.

本发明实施例提供一种非暂态计算机可读存储介质,所述非暂态计算机可读存储介质存储计算机指令,所述计算机指令使所述计算机执行上述计算方法。An embodiment of the present invention provides a non-transitory computer-readable storage medium, where the non-transitory computer-readable storage medium stores computer instructions, and the computer instructions cause the computer to execute the above calculation method.

本发明实施例提供的堆石混凝土浇筑过程数值模拟的计算方法和系统,通过设置有限元与离散元耦合分析方法进行力学分析,可以完整刻画出堆石块在应力作用下的渐进堆积过程,计算结果更为准确。通过设置SPH方法分析自密实混凝土在堆石块的最终堆积形态中的浇筑过程及最终填充状态,能够准确模拟自密实混凝土的浇筑过程,不仅可以直观地观察到自密实混凝土在堆石孔隙中的流动状态,还可以了解堆石混凝土最终的密实状态。本发明节省了大量的人力成本、实验成本以及经济成本,并避免了实验方法难以在线观测的难题。The calculation method and system for numerical simulation of the pouring process of rockfill concrete provided by the embodiment of the present invention can completely describe the progressive accumulation process of rockfill blocks under the action of stress by setting the finite element and discrete element coupling analysis method for mechanical analysis. The results are more accurate. By setting the SPH method to analyze the pouring process and final filling state of the self-compacting concrete in the final stacking form of the rockfill block, the pouring process of the self-compacting concrete can be accurately simulated, and the self-compacting concrete in the rockfill pores can be observed intuitively. The flow state can also be used to understand the final compaction state of the rockfill concrete. The invention saves a lot of labor cost, experiment cost and economic cost, and avoids the problem that the experimental method is difficult to observe online.

附图说明Description of drawings

为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the accompanying drawings that need to be used in the description of the embodiments or the prior art. Obviously, the accompanying drawings in the following description These are some embodiments of the present invention. For those of ordinary skill in the art, other drawings can also be obtained according to these drawings without creative efforts.

图1为本发明堆石混凝土浇筑过程数值模拟的计算方法实施例的流程图;Fig. 1 is the flow chart of the calculation method embodiment of the numerical simulation of rockfill concrete pouring process of the present invention;

图2为本发明堆石混凝土浇筑过程数值模拟的计算系统实施例的模块图;2 is a block diagram of an embodiment of a computing system for numerical simulation of the pouring process of rockfill concrete according to the present invention;

图3为本发明实施例中的堆石混凝土浇筑过程数值模拟的计算设备的框架示意图。3 is a schematic frame diagram of a computing device for numerical simulation of a pouring process of rockfill concrete in an embodiment of the present invention.

具体实施方式Detailed ways

为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。In order to make the purposes, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments These are some embodiments of the present invention, but not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.

图1为本发明堆石混凝土浇筑过程数值模拟的计算方法实施例的流程图,如图1所示,包括:S101、根据有限元/离散元耦合分析方法FEM/DEM,计算分析堆石块初始状态三维模型的堆积过程及最终堆积形态;S102、在计算分析所述堆石块初始状态三维模型的堆积过程及最终堆积形态后,基于光滑粒子流体动力学SPH方法,计算分析自密实混凝土在堆石块的最终堆积形态中的浇筑过程及最终填充状态。1 is a flowchart of an embodiment of a calculation method for numerical simulation of rockfill concrete pouring process according to the present invention, as shown in FIG. 1 , including: S101 , according to the finite element/discrete element coupling analysis method FEM/DEM, calculate and analyze the initial stage of the rockfill block The stacking process and final stacking form of the state three-dimensional model; S102, after calculating and analyzing the stacking process and final stacking form of the initial state three-dimensional model of the rockfill block, based on the smooth particle hydrodynamics SPH method, calculate and analyze the self-compacting concrete in the stack The pouring process and final filling state in the final stacked form of the stone.

需要说明的是,由于物理模型实验中比尺效应的影像,现有模型尺度的物模实验不能正确反应真实结构的强度,同时原型物模实验又受制于实验条件的限制。通过将模拟结构应力和变形的有限单元法和能够跟踪块体运动的离散单元法结合起来称之为有限元/离散元耦合分析方法FEM/DEM,该方法保留了有限元和离散元各自的优势,能够解决多体运动学问题和破坏力学问题,在计算堆石块的最终堆积形态的过程中,存在块石体个数多,受力复杂、堆石体随受力变化会产生形变的特点,采用有限元与离散元耦合分析方法进行力学分析,可以完整刻画出堆石块在应力作用下的渐进堆积过程,计算结果更为准确。It should be noted that due to the image of the scale effect in the physical model experiment, the physical model experiment at the existing model scale cannot correctly reflect the strength of the real structure, and the prototype physical model experiment is limited by the experimental conditions. By combining the finite element method, which simulates structural stress and deformation, and the discrete element method, which can track the motion of the block, it is called the finite element/discrete element coupled analysis method FEM/DEM, which retains the advantages of each of the finite element and discrete elements. , which can solve the problem of multi-body kinematics and failure mechanics. In the process of calculating the final stacking form of rockfill blocks, there are many rock bodies, complex forces, and rockfill bodies will deform with the change of force. , using the finite element and discrete element coupling analysis method for mechanical analysis, the progressive accumulation process of rockfill blocks under the action of stress can be completely described, and the calculation results are more accurate.

SPH(Smoothed Particle Hydrodynamics)方法是光滑粒子流体动力学方法的缩写,是在近20多年来逐步发展起来的一种无网格方法。该方法的基本思想是将连续的流体(或固体)用相互作用的质点组来描述,各个物质点上承载各种物理量,包括质量、速度等,通过求解质点组的动力学方程和跟踪每个质点的运动轨道,求得整个系统的力学行为。本发明实施例使用SPH方法模拟自密实混凝土在堆石体中的浇筑过程及并计算最终的流态。SPH (Smoothed Particle Hydrodynamics) method is the abbreviation of smooth particle hydrodynamics method, which is a meshless method gradually developed in the past 20 years. The basic idea of this method is to describe a continuous fluid (or solid) with an interacting particle group, and each material point carries various physical quantities, including mass, velocity, etc., by solving the dynamic equation of the particle group and tracking each The motion orbit of the particle is obtained, and the mechanical behavior of the whole system is obtained. In the embodiment of the present invention, the SPH method is used to simulate the pouring process of the self-compacting concrete in the rockfill body and calculate the final flow state.

本发明实施例提供的堆石混凝土浇筑过程数值模拟的计算方法,通过设置有限元与离散元耦合分析方法进行力学分析,可以完整刻画出堆石块在应力作用下的渐进堆积过程,计算结果更为准确。通过设置SPH方法分析自密实混凝土在堆石块的最终堆积形态中的浇筑过程及最终填充状态,能够准确模拟自密实混凝土的浇筑过程,不仅可以直观地观察到自密实混凝土在堆石孔隙中的流动状态,还可以了解堆石混凝土最终的密实状态。本发明节省了大量的人力成本、实验成本以及经济成本,并避免了实验方法难以在线观测的难题。In the calculation method for numerical simulation of the pouring process of rockfill concrete provided by the embodiment of the present invention, by setting the finite element and discrete element coupling analysis method for mechanical analysis, the progressive accumulation process of rockfill blocks under the action of stress can be completely described, and the calculation results are more accurate. to be accurate. By setting the SPH method to analyze the pouring process and final filling state of the self-compacting concrete in the final stacking form of the rockfill block, the pouring process of the self-compacting concrete can be accurately simulated, and the self-compacting concrete in the rockfill pores can be observed intuitively. The flow state can also be used to understand the final compaction state of the rockfill concrete. The invention saves a lot of labor cost, experiment cost and economic cost, and avoids the problem that the experimental method is difficult to observe online.

基于上述实施例,所述根据有限元/离散元耦合分析方法FEM/DEM,计算分析堆石块初始状态三维模型的堆积过程及最终堆积形态,之前还包括:建立堆石块初始状态三维模型,所述堆石块包括若干个块石体。Based on the above embodiment, the method of calculating and analyzing the stacking process and the final stacking form of the three-dimensional model of the initial state of the rockfill block according to the finite element/discrete element coupling analysis method FEM/DEM further includes: establishing a three-dimensional model of the initial state of the rockfill block, The rock pile includes several rock bodies.

基于上述实施例,所述建立堆石块初始状态三维模型,具体包括:使用三维建模软件如GID,根据堆石块的预先确定的参数,建立堆石块初始状态三维模型;或者,对所述块石体外形轮廓进行3D扫描形成输入数据输入至三维建模软件中,建立堆石块初始状态三维模型。Based on the above embodiment, the establishment of a three-dimensional model of the initial state of the rockfill specifically includes: using three-dimensional modeling software such as GID to establish a three-dimensional model of the initial state of the rockfill according to predetermined parameters of the rockfill; The outline of the rock body is 3D scanned to form the input data and input into the 3D modeling software to establish a 3D model of the initial state of the rockfill.

由于每个工程或者实验中需要的石块大小、种类、数量都是不同的是,所以根据堆石块的实际参数建立三维模型,在建立前显然需要采集若干真实参数,例如石块的硬度、大小、密度等等。Since the size, type, and quantity of stones required in each project or experiment are different, a three-dimensional model is established based on the actual parameters of the rockfill. size, density, etc.

本发明实施例在采集块石体的外形轮廓时,具体可以采用3D扫描仪或者轮廓测量仪进行扫描,获得块石体的外形参数,本发明实施例不做具体限定。In the embodiment of the present invention, when collecting the outline of the block stone body, a 3D scanner or a contour measuring instrument may be used for scanning to obtain the shape parameters of the block stone body, which is not specifically limited in the embodiment of the present invention.

本实施例有助于后续进行有限元/离散元耦合分析方法分析堆积过程时获得更准确、真实的结果。This embodiment is helpful for obtaining more accurate and real results when analyzing the stacking process with the finite element/discrete element coupling analysis method in the following.

基于上述实施例,所述计算分析所述堆石块初始状态三维模型的堆积过程及最终堆积形态,以及所述基于光滑粒子流体动力学SPH方法,计算分析自密实混凝土在堆石块的最终堆积形态中的浇筑过程及最终填充状态,之间还包括:建立自密实混凝土模型,所述自密实混凝土模型为Bingham流变模型。Based on the above embodiment, the calculation and analysis of the stacking process and final stacking shape of the three-dimensional model of the initial state of the rockfill block, and the calculation and analysis of the final stacking of the self-compacting concrete in the rockfill block based on the smooth particle hydrodynamic SPH method The pouring process and the final filling state in the form further include: establishing a self-compacting concrete model, where the self-compacting concrete model is a Bingham rheological model.

具体地,自密实混凝土模型为:Specifically, the self-compacting concrete model is:

Figure BDA0001714410050000051
Figure BDA0001714410050000051

其中,τ0为屈服强度,μp为塑性粘度,τ为剪切应力;γ为剪切应变率。Among them, τ 0 is the yield strength, μ p is the plastic viscosity, τ is the shear stress; γ is the shear strain rate.

需要说明的是,宾汉模型(Bingham model):将缓冲器、滑块并联,再与弹簧串联而成的反映岩土黏性、弹性、塑性流变特征及其过程的力学模型。It should be noted that the Bingham model: a mechanical model that reflects the viscous, elastic, plastic rheological characteristics of rock and soil and its process by connecting the buffer and the slider in parallel, and then connecting them in series with the spring.

基于上述实施例,所述基于光滑粒子流体动力学SPH方法,计算分析自密实混凝土在堆石块的最终堆积形态中的浇筑过程及最终填充状态,具体包括:基于SPH方法,逐个时间步长迭代更新所述自密实混凝土模型中的流体粒子的密度、位置和速度;基于所述流体粒子的密度、位置和速度,分析自密实混凝土在堆石块的最终堆积形态中的浇筑过程及最终填充状态。Based on the above embodiment, the calculation and analysis of the pouring process and the final filling state of the self-compacting concrete in the final stacking form of the rockfill block based on the smooth particle hydrodynamic SPH method specifically include: based on the SPH method, iterating over time step by step Update the density, position and velocity of the fluid particles in the self-compacting concrete model; based on the density, position and velocity of the fluid particles, analyze the pouring process and final filling state of the self-compacting concrete in the final stacking form of the rockfill block .

下面对本实施例所需的函数及方程做出进一步地说明,本发明实施例中的SPH离散形式为:The required functions and equations of the present embodiment are further described below, and the discrete form of SPH in the embodiment of the present invention is:

Figure BDA0001714410050000052
Figure BDA0001714410050000052

其中,其中:N为支持域内相邻粒子总的个数;xj为粒子j的位置;mj及ρj分别为粒子j的质量、密度;核函数W(x-xj,h)与粒子间的距离和光滑长度h有关。where: N is the total number of adjacent particles in the support domain; x j is the position of particle j ; m j and ρ j are the mass and density of particle j, respectively; The distance is related to the smooth length h.

本实施例选取核函数为:The kernel function selected in this embodiment is:

Figure BDA0001714410050000061
Figure BDA0001714410050000061

其中,d为空间维数;h为光滑长度。Among them, d is the spatial dimension; h is the smooth length.

Figure BDA0001714410050000062
Figure BDA0001714410050000062

其中,C为归一化常数。where C is a normalization constant.

本实施例选取的状态方程为:The state equation selected in this embodiment is:

Figure BDA0001714410050000063
Figure BDA0001714410050000063

其中,B作为初始压力,用于限制密度的最大改变量,本发明取

Figure BDA0001714410050000064
ρ0是粒子的初始密度;γ为一恒定系数;C为人为音速,出于数值计算稳定性和时间步长考虑,通常取一个比真实音速小得多的值,一般取流场最大流速值的10倍,即
Figure BDA0001714410050000065
H为自由液面高度。Among them, B is used as the initial pressure to limit the maximum change of density, and the present invention takes
Figure BDA0001714410050000064
ρ 0 is the initial density of the particle; γ is a constant coefficient; C is the artificial speed of sound. For the stability of numerical calculation and the time step, a value much smaller than the true speed of sound is usually taken, and the maximum flow velocity value of the flow field is generally taken. 10 times the
Figure BDA0001714410050000065
H is the free surface height.

粒子的位移关系式为:The displacement relation of particles is:

Figure BDA0001714410050000066
Figure BDA0001714410050000066

其中,r为粒子点的位置矢量,u为粒子的速度矢量。Among them, r is the position vector of the particle point, and u is the velocity vector of the particle.

本发明实施例中的边界处理方法为:固壁粒子参与控制方程的计算,但它们位置固定,在本发明中位移设置为0。边界粒子与流体粒子间的相互作用一般可采用Lennard-Jones排斥力方法,即假定当流体粒子距离固壁足够近时,为阻止内部流体粒子穿越固壁边界后飞散而引起计算崩溃,边界粒子对行近的流体粒子施加一个中心排斥力,排斥力大小由相互距离、投影速度、高度位置等确定。自由表面采用Koshizuka等提出的算法实现Dirichlet条件,界面处粒子压力赋值为0或一定的外压力值。The boundary processing method in the embodiment of the present invention is as follows: the solid wall particles participate in the calculation of the control equation, but their positions are fixed, and the displacement is set to 0 in the present invention. The Lennard-Jones repulsive force method can generally be used for the interaction between boundary particles and fluid particles, that is, it is assumed that when the fluid particles are close enough to the solid wall, in order to prevent the internal fluid particles from flying through the solid wall boundary and causing the calculation to collapse, the boundary particles will The approaching fluid particles exert a central repulsive force, and the magnitude of the repulsive force is determined by the mutual distance, projection speed, height position, etc. The free surface adopts the algorithm proposed by Koshizuka et al. to realize the Dirichlet condition, and the particle pressure at the interface is assigned as 0 or a certain external pressure value.

基于上述实施例,所述逐个时间步长迭代更新所述自密实混凝土模型中的流体粒子的密度、位置和速度,具体包括:逐个时间步长,通过连续性方程计算所述流体粒子的密度变化,捕捉自由表面粒子并进行自由表面粒子的密度校正,获得所述自密实混凝土模型中每一时刻的密度场;逐个时间步长,通过流体粒子位移方程更新各个流体粒子运动的轨迹,获取任一流体粒子每一时刻的位置;逐个时间步长,通过状态方程确定任一流体粒子的压力,基于每一流体粒子的压力计算外力产生的加速度,基于所述加速度通过动量方程计算所述自密实混凝土模型中每一时刻的速度场。Based on the above embodiment, the iteratively updating the density, position and velocity of the fluid particles in the self-compacting concrete model by time step, specifically includes: calculating the density change of the fluid particles through a continuity equation, time by step. , capture the free surface particles and perform density correction of the free surface particles to obtain the density field at each moment in the self-compacting concrete model; each time step, update the trajectory of each fluid particle motion through the fluid particle displacement equation, and obtain any The position of the fluid particle at each moment; for each time step, the pressure of any fluid particle is determined through the equation of state, the acceleration generated by the external force is calculated based on the pressure of each fluid particle, and the self-compacting concrete is calculated through the momentum equation based on the acceleration. The velocity field at each moment in the model.

需要说明的是,粒子搜索采用树形搜索法或关联链表搜索技术。It should be noted that the particle search adopts the tree search method or the associated linked list search technology.

本发明实施例提供的堆石混凝土浇筑过程数值模拟的计算方法,通过设置逐个时间步长迭代更新所述自密实混凝土模型中的流体粒子的密度、位置和速度,能够准确计算堆石混凝土的浇筑过程。The calculation method for the numerical simulation of the pouring process of rockfill concrete provided by the embodiment of the present invention can accurately calculate the pouring of rockfill concrete by iteratively updating the density, position and velocity of the fluid particles in the self-compacting concrete model by setting time steps. process.

基于上述实施例,所述根据有限元/离散元耦合分析方法FEM/DEM,计算分析堆石块初始状态三维模型的堆积过程及最终堆积形态,具体包括:根据有限元分析法FEM获得所述堆石块初始状态三维模型中每个块石体连续介质力学行为;根据离散元分析法DEM获得所述堆石块初始状态三维模型中块石体间的非连续介质力学行为;将所有块石体的连续介质力学行为和非连续介质力学行为代入堆石块初始状态三维模型的动力平衡方程,获得所述堆积过程及最终堆积形态。Based on the above embodiment, calculating and analyzing the stacking process and final stacking shape of the three-dimensional model of the initial state of rockfill blocks according to the finite element/discrete element coupling analysis method FEM/DEM specifically includes: obtaining the stack according to the finite element analysis method FEM The mechanical behavior of the continuum of each rock mass in the three-dimensional model of the initial state of the rock block; the discontinuous medium mechanical behavior between the rock bodies in the three-dimensional model of the initial state of the rockfill block is obtained according to the discrete element analysis method DEM; The continuum mechanical behavior and discontinuous medium mechanical behavior are substituted into the dynamic balance equation of the three-dimensional model of the initial state of the rockfill block, and the accumulation process and final accumulation form are obtained.

需要说明的是,有限元分析(FEA,Finite Element Analysis)是利用数学近似的方法对真实物理系统(几何和载荷工况)进行模拟。利用简单而又相互作用的元素,就可用有限数量的未知量去逼近无限未知量的真实系统。它将求解域看成是由许多称为有限元的小的互连子域组成,对每一单元假定一个合适的近似解,然后推导求解这个域总的满足条件,从而得到问题的解。It should be noted that Finite Element Analysis (FEA, Finite Element Analysis) uses mathematical approximation to simulate the real physical system (geometry and load conditions). Using simple but interacting elements, a finite number of unknowns can be used to approximate a real system with infinite unknowns. It regards the solution domain as consisting of many small interconnected sub-domains called finite elements, assumes a suitable approximate solution for each element, and then derives the total satisfying conditions for solving this domain, so as to obtain the solution of the problem.

连续介质力学(Continuum mechanics)是物理学当中的一个分支,是处理包括固体和流体的在内的所谓“连续介质”宏观性质的力学。例如,质量守恒、动量和角动量定理、能量守恒等。连续介质力学认为真实流体或固体所占有的空间可以近似地看作连续地无空隙地充满着“质点”(即微观上充分大、宏观上充分小的分子团)。质点所具有的宏观物理量(如质量、速度、压力、温度等)满足一切应该遵循的物理定律,例如质量守恒定律、牛顿运动定律、能量守恒定律、热力学定律以及扩散、粘性及热传导等输运性质。这一假设忽略物质的具体微观结构(对固体的微观结构研究属于凝聚态物理学的范畴),而用一组偏微分方程来表达宏观物理量(如质量,数度,压力等)。Continuum mechanics is a branch of physics that deals with the macroscopic properties of so-called "continuum", including solids and fluids. For example, conservation of mass, theorems of momentum and angular momentum, conservation of energy, etc. In continuum mechanics, the space occupied by real fluids or solids can be approximately regarded as being filled with "particles" continuously without voids (that is, molecular groups that are sufficiently large on the microscopic scale and sufficiently small on the macroscopic scale). The macroscopic physical quantities (such as mass, velocity, pressure, temperature, etc.) of a particle satisfy all physical laws that should be followed, such as the law of conservation of mass, Newton's law of motion, law of conservation of energy, laws of thermodynamics, and transport properties such as diffusion, viscosity, and heat conduction . This assumption ignores the specific microstructure of matter (the study of the microstructure of solids belongs to the category of condensed matter physics), and uses a set of partial differential equations to express macroscopic physical quantities (such as mass, degrees, pressure, etc.).

离散元分析法是专门用来解决非连续介质问题的数值模拟方法。该方法把堆石块体视为由离散的块石体和块石体间的节理面所组成,允许块石体平移、转动和变形,而节理面可被压缩、分离或滑动。因此,堆石块被看作一种不连续的离散介质。其内部可存在大位移、旋转和滑动乃至块体的分离,从而可以较真实地模拟堆石块中的非线性大变形特征。Discrete element analysis is a numerical simulation method specially used to solve the problem of discontinuous medium. This method regards the rockfill block as consisting of discrete blocks and joint surfaces between the blocks, allowing the blocks to translate, rotate and deform, while the joint surfaces can be compressed, separated or slipped. Therefore, the rockfill is regarded as a discontinuous discrete medium. There can be large displacement, rotation and sliding and even block separation inside, so that the nonlinear large deformation characteristics in rockfill blocks can be simulated more realistically.

离散元法的一般求解过程为:将求解空间离散为离散元单元阵,并根据实际问题用合理的连接元件将相邻两单元连接起来;单元间相对位移是基本变量,由力与相对位移的关系可得到两单元间法向和切向的作用力;对单元在各个方向上与其它单元间的作用力以及其它物理场对单元作用所引起的外力求合力和合力矩,根据牛顿运动第二定律可以求得单元的加速度;对其进行时间积分,进而得到单元的速度和位移。从而得到所有单元在任意时刻的速度、加速度、角速度、线位移和转角等物理量。The general solution process of the discrete element method is: discretize the solution space into a discrete element array, and connect two adjacent elements with reasonable connecting elements according to the actual problem; the relative displacement between elements is the basic variable, which is determined by the force and relative displacement. The normal and tangential forces between the two elements can be obtained; the forces between the element and other elements in all directions and the external forces caused by other physical fields acting on the element can be obtained. According to Newton's second law of motion The acceleration of the element can be obtained; it is time-integrated to obtain the velocity and displacement of the element. Thereby, physical quantities such as velocity, acceleration, angular velocity, linear displacement and rotation angle of all units at any time can be obtained.

在上述实施例的基础上,本实施例来进一步地说明根据有限元分析法FEM获得所述堆石块初始状态三维模型中每个块石体连续介质力学行为的具体过程。On the basis of the above embodiment, this embodiment further illustrates the specific process of obtaining the mechanical behavior of the continuum of each rock mass in the three-dimensional model of the initial state of the rockfill block according to the finite element analysis method FEM.

具体地,所述连续介质力学行为包括:网格单元的节点变形引起的内力向量、网格单元的节点的外力向量以及质量矩阵。Specifically, the mechanical behavior of the continuum includes: the internal force vector caused by the deformation of the nodes of the grid element, the external force vector of the nodes of the grid element, and the mass matrix.

与之相对应地,所述根据有限元分析法获得三维模型中每个块石体的连续介质力学行为,具体为:Correspondingly, obtaining the continuum mechanical behavior of each rock body in the three-dimensional model according to the finite element analysis method is specifically:

对堆石块初始状态三维模型进行网格划分,获得由网格单元构成的有限元模型。The three-dimensional model of the initial state of the rockfill block is meshed to obtain a finite element model composed of mesh elements.

根据网格单元体积、标准型函数以及柯西应力张量获得所述内力向量。The internal force vector is obtained from the mesh element volume, the canonical function, and the Cauchy stress tensor.

根据网格单元未变形时的体积、与该网格单元未变形时的体积对应的密度、以及所述标准型函数获得所述质量矩阵。The mass matrix is obtained from the undeformed volume of the mesh element, the density corresponding to the undeformed volume of the mesh element, and the canonical function.

根据块石体的体力、面力以及单元表面积,获得所述外力向量。The external force vector is obtained according to the physical force, the surface force and the surface area of the unit of the rock body.

进一步地,内力向量按下式计算:Further, the internal force vector is calculated as follows:

Figure BDA0001714410050000091
Figure BDA0001714410050000091

其中,V表示单元体积,N为标准形函数,T为柯西应力张量。where V is the unit volume, N is the standard shape function, and T is the Cauchy stress tensor.

节点外力向量按下式计算:The nodal external force vector is calculated as follows:

Figure BDA0001714410050000094
Figure BDA0001714410050000094

其中,b表示体力,t表示面力,S为单元表面积。Among them, b is the physical force, t is the surface force, and S is the unit surface area.

质量矩阵按下式计算:The mass matrix is calculated as:

Figure BDA0001714410050000092
Figure BDA0001714410050000092

其中,V0为单元未变形时的体积,ρ0表示与V0对应的密度。Among them, V 0 is the volume of the unit when it is not deformed, and ρ 0 represents the density corresponding to V 0 .

在上述实施例的基础上,本实施例来进一步地说明根据离散元分析法DEM获得所述堆石块初始状态三维模型中块石体间的非连续介质力学行为的具体过程。On the basis of the above embodiment, this embodiment further illustrates the specific process of obtaining the mechanical behavior of discontinuous medium between the rock bodies in the three-dimensional model of the initial state of the rockfill block according to the discrete element analysis method DEM.

具体地,非连续介质力学行为包括:块石体间的接触力。Specifically, the mechanical behavior of the discontinuous medium includes: the contact force between the boulders.

与之相对应地,所述根据离散元分析法获得三维模型中块石体间的非连续介质力学行为,具体为:Correspondingly, the mechanical behavior of the discontinuous medium between the blocks in the three-dimensional model obtained according to the discrete element analysis method is specifically:

通过非二叉树搜索接触算法(接触判断方法则采用Munjiza和Andrews提出的非二叉树搜索接触算法,即高效的NBS算法,这种算法占用内存小,对于块体单元尺寸相近的情况可以实现单元间的快速接触搜索)获得三维模型中的接触体和目标体。Through the non-binary tree search contact algorithm (the contact judgment method adopts the non-binary tree search contact algorithm proposed by Munjiza and Andrews, that is, the efficient NBS algorithm, this algorithm occupies a small memory, and can achieve fast inter-unit for the case of similar block unit sizes. Contact search) to obtain contact and target bodies in the 3D model.

通过非二叉树搜索接触算法获得三维模型中的接触体和目标体;Obtain the contact body and the target body in the 3D model through the non-binary tree search contact algorithm;

根据公式:

Figure BDA0001714410050000093
获得接触体和目标体之间的法向接触力;其中,βc和βt分别表示接触体和目标体,m和n分别表示用于离散接触体和目标体的有限单元数目,
Figure BDA0001714410050000101
表示接触势,S表示相互嵌入部分的边界,n表示该边界的外法线方向矢量。According to the formula:
Figure BDA0001714410050000093
Obtain the normal contact force between the contact and target; where β c and β t denote the contact and target, respectively, m and n denote the number of finite elements used to discretize the contact and target, respectively,
Figure BDA0001714410050000101
represents the contact potential, S represents the boundary of the mutually embedded part, and n represents the outer normal direction vector of the boundary.

根据公式:

Figure BDA0001714410050000102
获得接触体和目标体之间的切向接触力;其中kt为切向刚度,ηt为切向粘性阻尼系数,dt和vt分别表示切向相对位移和相对速度。According to the formula:
Figure BDA0001714410050000102
The tangential contact force between the contact body and the target body is obtained; where k t is the tangential stiffness, η t is the tangential viscous damping coefficient, and d t and v t represent the tangential relative displacement and relative velocity, respectively.

若判断

Figure BDA0001714410050000103
则根据库伦摩擦定律,重新根据公式
Figure BDA0001714410050000104
计算所述切向接触力。If judged
Figure BDA0001714410050000103
Then according to Coulomb's friction law, according to the formula again
Figure BDA0001714410050000104
Calculate the tangential contact force.

根据法向接触力和切向接触力的矢量和作为所述接触力。The contact force is taken as the vector sum of the normal contact force and the tangential contact force.

在上述实施例的基础上,堆石块初始状态三维模型的动力平衡方程的有限元离散形式为:On the basis of the above embodiment, the finite element discrete form of the dynamic balance equation of the three-dimensional model of the initial state of the rockfill block is:

Figure BDA0001714410050000105
Figure BDA0001714410050000105

x=X+u;x=X+u;

其中,M,C分别表示质量矩阵和阻尼矩阵,fint表示变形引起的内力向量,fc表示接触力,fext表示除接触力外的所有外力向量,

Figure BDA0001714410050000106
Figure BDA0001714410050000107
分别表示有限元节点的加速度和速度,X为有限元节点的初始位置向量,u为有限元节点的位移向量。Among them, M and C represent the mass matrix and damping matrix, respectively, f int represents the internal force vector caused by deformation, f c represents the contact force, f ext represents all the external force vectors except the contact force,
Figure BDA0001714410050000106
and
Figure BDA0001714410050000107
respectively represent the acceleration and velocity of the finite element node, X is the initial position vector of the finite element node, and u is the displacement vector of the finite element node.

作为一个优选实施例,本发明实施例中的堆石混凝土浇筑过程数值模拟的计算方法具体包括:As a preferred embodiment, the calculation method for the numerical simulation of the pouring process of rockfill concrete in the embodiment of the present invention specifically includes:

建立堆石块初始状态三维模型,所述堆石块包括若干个块石体。A three-dimensional model of the initial state of the rockfill block is established, and the rockfill block includes several rock bodies.

根据有限元/离散元耦合分析方法FEM/DEM,计算分析堆石块初始状态三维模型的堆积过程及最终堆积形态。According to the finite element/discrete element coupling analysis method FEM/DEM, the accumulation process and final accumulation form of the 3D model of the initial state of rockfill blocks are calculated and analyzed.

将堆石块初始状态三维模型的最终堆积形态以边界的形式进行输入,并建立自密实混凝土模型,所述自密实混凝土模型为Bingham流变模型。The final stacking form of the three-dimensional model of the initial state of the rockfill block is input in the form of a boundary, and a self-compacting concrete model is established, and the self-compacting concrete model is the Bingham rheological model.

基于光滑粒子流体动力学SPH方法,计算分析自密实混凝土在堆石块的最终堆积形态中的浇筑过程及最终填充状态。Based on the smooth particle hydrodynamic SPH method, the pouring process and final filling state of self-compacting concrete in the final stacking form of rockfill blocks are calculated and analyzed.

本发明实施例使用Paraview软件作为后处理软件,对计算结果进行处理。需要说明的是,Paraview软件本身就是一种数值计算结果处理的成熟软件,主要是分析堆石混凝土浇筑过程的可视化处理,以及速度场,应力场的分析和最终浇筑状态记忆密实度的显示。In the embodiment of the present invention, Paraview software is used as the post-processing software to process the calculation result. It should be noted that Paraview software itself is a mature software for processing numerical calculation results. It mainly analyzes the visual processing of the pouring process of rockfill concrete, as well as the analysis of the velocity field, the stress field and the display of the memory compactness of the final pouring state.

建立堆石块初始状态三维模型,所述堆石块包括若干个块石体,具体包括:使用三维建模软件根据堆石块的预先确定的参数,建立堆石块初始状态三维模型;或者,对所述块石体外形轮廓进行3D扫描形成输入数据输入至三维建模软件中,建立堆石块初始状态三维模型。Establishing a three-dimensional model of the initial state of the rockfill block, where the rockfill block includes several rock bodies, specifically including: using a three-dimensional modeling software to establish a three-dimensional model of the initial state of the rockfill block according to predetermined parameters of the rockfill block; or, 3D scanning is performed on the outline of the rock body to form input data and input into the three-dimensional modeling software to establish a three-dimensional model of the initial state of the rockfill block.

根据有限元/离散元耦合分析方法FEM/DEM,计算分析堆石块初始状态三维模型的堆积过程及最终堆积形态,具体包括:根据有限元分析法FEM获得所述堆石块初始状态三维模型中每个块石体连续介质力学行为;根据离散元分析法DEM获得所述堆石块初始状态三维模型中块石体间的非连续介质力学行为;将所有块石体的连续介质力学行为和非连续介质力学行为代入堆石块初始状态三维模型的动力平衡方程,获得所述堆积过程及最终堆积形态。According to the finite element/discrete element coupling analysis method FEM/DEM, calculate and analyze the accumulation process and final accumulation form of the three-dimensional model of the initial state of the rockfill block, which specifically includes: obtaining the three-dimensional model of the initial state of the rockfill block according to the finite element analysis method FEM. The continuum mechanical behavior of each rock mass; the discontinuous media mechanical behavior between the rock bodies in the three-dimensional model of the initial state of the rockfill block is obtained according to the discrete element analysis method DEM; the continuum mechanical behavior and non-continuous media mechanical behavior of all rock bodies are The mechanical behavior of the continuum is substituted into the dynamic equilibrium equation of the three-dimensional model of the initial state of the rockfill block to obtain the stacking process and final stacking form.

基于光滑粒子流体动力学SPH方法,计算分析自密实混凝土在堆石块的最终堆积形态中的浇筑过程及最终填充状态,具体包括:基于SPH方法,逐个时间步长迭代更新所述自密实混凝土模型中的流体粒子的密度、位置和速度;基于所述流体粒子的密度、位置和速度,分析自密实混凝土在堆石块的最终堆积形态中的浇筑过程及最终填充状态。Based on the SPH method of smooth particle hydrodynamics, calculate and analyze the pouring process and final filling state of self-compacting concrete in the final stacking form of rockfill blocks, including: based on the SPH method, iteratively update the self-compacting concrete model by time step. The density, position and velocity of the fluid particles in the rockfill; based on the density, position and velocity of the fluid particles, the pouring process and final filling state of the self-compacting concrete in the final stacking form of the rockfill block are analyzed.

基于上述实施例,图2为本发明堆石混凝土浇筑过程数值模拟的计算系统实施例的模块图,如图2所示,包括:堆积分析模块201,用于根据有限元/离散元耦合分析方法FEM/DEM,计算分析堆石块初始状态三维模型的堆积过程及最终堆积形态;浇筑分析模块202,用于在计算分析所述堆石块初始状态三维模型的堆积过程及最终堆积形态后,基于光滑粒子流体动力学SPH方法,计算分析自密实混凝土在堆石块的最终堆积形态中的浇筑过程及最终填充状态。Based on the above embodiment, FIG. 2 is a block diagram of an embodiment of a computing system for numerical simulation of the pouring process of rockfill concrete according to the present invention, as shown in FIG. FEM/DEM, calculates and analyzes the stacking process and final stacking shape of the 3D model of the initial state of the rockfill block; the pouring analysis module 202 is used to calculate and analyze the stacking process and final stacking shape of the 3D model of the initial state of the rockfill block, based on The smooth particle hydrodynamic SPH method is used to calculate and analyze the pouring process and final filling state of self-compacting concrete in the final stacking form of rockfill.

本发明实施例的短信分发系统,可用于执行图1所示的堆石混凝土浇筑过程数值模拟的计算方法实施例的技术方案,其实现原理和技术效果类似,此处不再赘述。The short message distribution system of the embodiment of the present invention can be used to implement the technical solution of the calculation method embodiment of the numerical simulation of the pouring process of rockfill concrete shown in FIG.

基于上述实施例,图3为本发明实施例中的堆石混凝土浇筑过程数值模拟的计算设备的框架示意图。请参考图3,本发明实施例提供一种堆石混凝土浇筑过程数值模拟的计算设备,包括:处理器(processor)310、通信接口(Communications Interface)320、存储器(memory)330和总线340,其中,处理器310,通信接口320,存储器330通过总线340完成相互间的通信。处理器310可以调用存储器330中的逻辑指令,以执行如下方法,包括:根据有限元/离散元耦合分析方法FEM/DEM,计算分析堆石块初始状态三维模型的堆积过程及最终堆积形态;在计算分析所述堆石块初始状态三维模型的堆积过程及最终堆积形态后,基于光滑粒子流体动力学SPH方法,计算分析自密实混凝土在堆石块的最终堆积形态中的浇筑过程及最终填充状态。Based on the above embodiment, FIG. 3 is a schematic frame diagram of a computing device for numerical simulation of the pouring process of rockfill concrete in an embodiment of the present invention. Referring to FIG. 3 , an embodiment of the present invention provides a computing device for numerical simulation of a pouring process of rockfill concrete, including: a processor 310 , a communications interface 320 , a memory 330 and a bus 340 , wherein , the processor 310 , the communication interface 320 , and the memory 330 communicate with each other through the bus 340 . The processor 310 can call the logic instructions in the memory 330 to execute the following methods, including: calculating and analyzing the stacking process and final stacking shape of the three-dimensional model of the initial state of the rockfill block according to the finite element/discrete element coupling analysis method FEM/DEM; After calculating and analyzing the stacking process and final stacking shape of the 3D model of the initial state of the rockfill block, based on the smooth particle hydrodynamic SPH method, the pouring process and final filling state of the self-compacting concrete in the final stacking shape of the rockfill block were calculated and analyzed. .

本发明实施例公开一种计算机程序产品,所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序,所述计算机程序包括程序指令,当所述程序指令被计算机执行时,计算机能够执行上述各方法实施例所提供的计算方法,例如包括:根据有限元/离散元耦合分析方法FEM/DEM,计算分析堆石块初始状态三维模型的堆积过程及最终堆积形态;在计算分析所述堆石块初始状态三维模型的堆积过程及最终堆积形态后,基于光滑粒子流体动力学SPH方法,计算分析自密实混凝土在堆石块的最终堆积形态中的浇筑过程及最终填充状态。An embodiment of the present invention discloses a computer program product, where the computer program product includes a computer program stored on a non-transitory computer-readable storage medium, the computer program includes program instructions, and when the program instructions are executed by a computer, The computer can execute the calculation methods provided by the above method embodiments, for example, including: according to the finite element/discrete element coupling analysis method FEM/DEM, calculating and analyzing the stacking process and final stacking shape of the three-dimensional model of the initial state of the rockfill; After the stacking process and final stacking form of the three-dimensional model of the initial state of the rockfill block, based on the smooth particle hydrodynamic SPH method, the pouring process and final filling state of the self-compacting concrete in the final stacking state of the rockfill block are calculated and analyzed.

基于上述实施例,本发明实施例提供一种非暂态计算机可读存储介质,所述非暂态计算机可读存储介质存储计算机指令,所述计算机指令使所述计算机执行上述各方法实施例所提供的计算方法,例如包括:根据有限元/离散元耦合分析方法FEM/DEM,计算分析堆石块初始状态三维模型的堆积过程及最终堆积形态;在计算分析所述堆石块初始状态三维模型的堆积过程及最终堆积形态后,基于光滑粒子流体动力学SPH方法,计算分析自密实混凝土在堆石块的最终堆积形态中的浇筑过程及最终填充状态。Based on the foregoing embodiments, embodiments of the present invention provide a non-transitory computer-readable storage medium, where the non-transitory computer-readable storage medium stores computer instructions, and the computer instructions cause the computer to execute the methods described in the foregoing method embodiments. The provided calculation methods include, for example: calculating and analyzing the stacking process and final stacking shape of the three-dimensional model of the initial state of the rockfill block according to the finite element/discrete element coupling analysis method FEM/DEM; calculating and analyzing the three-dimensional model of the initial state of the rockfill block. After the accumulation process and final accumulation form of rockfill blocks, the pouring process and final filling state of self-compacting concrete in the final accumulation form of rockfill blocks were calculated and analyzed based on the smooth particle hydrodynamic SPH method.

本领域普通技术人员可以理解:实现上述设备实施例或方法实施例仅仅是示意性的,其中所述处理器和所述存储器可以是物理上分离的部件也可以不是物理上分离的,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。Those of ordinary skill in the art can understand that the implementation of the above device embodiments or method embodiments is merely illustrative, wherein the processor and the memory may be physically separated components or may not be physically separated, that is, they may be located in One place, or it can be distributed over multiple network elements. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution in this embodiment. Those of ordinary skill in the art can understand and implement it without creative effort.

通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如U盘、移动硬盘、ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。From the description of the above embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus a necessary general hardware platform, and certainly can also be implemented by hardware. Based on this understanding, the above-mentioned technical solutions can be embodied in the form of software products in essence or the parts that make contributions to the prior art, and the computer software products can be stored in computer-readable storage media, such as U disk, mobile hard disk , ROM/RAM, magnetic disk, optical disk, etc., including several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to perform the methods described in various embodiments or parts of embodiments.

本发明实施例提供的堆石混凝土浇筑过程数值模拟的计算方法和系统,通过设置有限元与离散元耦合分析方法进行力学分析,可以完整刻画出堆石块在应力作用下的渐进堆积过程,计算结果更为准确。通过设置SPH方法分析自密实混凝土在堆石块的最终堆积形态中的浇筑过程及最终填充状态,能够准确模拟自密实混凝土的浇筑过程,不仅可以直观地观察到自密实混凝土在堆石孔隙中的流动状态,还可以了解堆石混凝土最终的密实状态。本发明节省了大量的人力成本、实验成本以及经济成本,并避免了实验方法难以在线观测的难题。通过设置逐个时间步长迭代更新所述自密实混凝土模型中的流体粒子的密度、位置和速度,能够准确计算堆石混凝土的浇筑过程。The calculation method and system for numerical simulation of the pouring process of rockfill concrete provided by the embodiment of the present invention can completely describe the progressive accumulation process of rockfill blocks under the action of stress by setting the finite element and discrete element coupling analysis method for mechanical analysis. The results are more accurate. By setting the SPH method to analyze the pouring process and final filling state of the self-compacting concrete in the final stacking form of the rockfill block, the pouring process of the self-compacting concrete can be accurately simulated, and the self-compacting concrete in the rockfill pores can be observed intuitively. The flow state can also be used to understand the final compaction state of the rockfill concrete. The invention saves a lot of labor cost, experiment cost and economic cost, and avoids the problem that the experimental method is difficult to observe online. By setting iteratively updating the density, position and velocity of fluid particles in the self-compacting concrete model by time step, the pouring process of rockfill concrete can be accurately calculated.

最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, but not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that it can still be The technical solutions described in the foregoing embodiments are modified, or some technical features thereof are equivalently replaced; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims (7)

1. A calculation method for numerical simulation in a rockfill concrete pouring process is characterized by comprising the following steps:
calculating and analyzing the stacking process and the final stacking form of the three-dimensional model of the stacking stone initial state according to a finite element/discrete element coupling analysis method FEM/DEM;
after the stacking process and the final stacking form of the three-dimensional model of the initial state of the stacking blocks are calculated and analyzed, the pouring process and the final filling state of the self-compacting concrete in the final stacking form of the stacking blocks are calculated and analyzed based on a smooth particle fluid dynamics (SPH) method;
iteratively updating the density, the position and the speed of fluid particles in the self-compacting concrete model one by one according to a time step based on an SPH (shortest path H) method;
analyzing the pouring process and the final filling state of the self-compacting concrete in the final stacking form of the stacking blocks based on the density, the position and the speed of the fluid particles;
wherein, the iteratively updating the density, the position and the speed of the fluid particles in the self-compacting concrete model one by one time step specifically comprises:
calculating the density change of the fluid particles one by one according to a continuity equation, capturing free surface particles and correcting the density of the free surface particles to obtain a density field of the self-compacting concrete model at each moment;
updating the motion track of each fluid particle through a fluid particle displacement equation one by one time step to obtain the position of any fluid particle at each moment;
determining the pressure of any fluid particle through a state equation one by one according to the time step, calculating the acceleration generated by the external force based on the pressure of each fluid particle, and calculating the velocity field of each moment in the self-compacting concrete model through a momentum equation based on the acceleration;
the method for calculating and analyzing the stacking process and the final stacking form of the three-dimensional model of the stacking stone initial state according to the finite element/discrete element coupling analysis method FEM/DEM specifically comprises the following steps:
obtaining the continuous medium mechanical behavior of each rock block in the three-dimensional model of the stacking initial state of the rock blocks according to a finite element analysis (FEM);
obtaining the mechanical behavior of a non-continuous medium among the rock masses in the three-dimensional model of the initial state of the rock pile according to a discrete element analysis (DEM);
substituting the continuous medium mechanical behaviors and the non-continuous medium mechanical behaviors of all the rock masses into a dynamic balance equation of the three-dimensional model of the initial state of the rock mass pile to obtain the stacking process and the final stacking form;
the continuous medium mechanical behavior comprises: an internal force vector caused by node deformation of the grid unit, an external force vector of the node of the grid unit and a mass matrix;
the method for obtaining the continuous medium mechanical behavior of each block stone body in the three-dimensional model according to the finite element analysis method specifically comprises the following steps:
carrying out mesh division on the three-dimensional model of the rock-fill block in the initial state to obtain a finite element model consisting of mesh units;
obtaining the internal force vector according to the grid unit volume, the standard type function and the Cauchy stress tensor;
obtaining the quality matrix according to the volume of the grid unit when the grid unit is not deformed, the density corresponding to the volume of the grid unit when the grid unit is not deformed and the standard function;
obtaining the external force vector according to the physical strength, the surface force and the unit surface area of the stone body;
the mechanical behavior of the non-continuous medium comprises: contact force between the stone bodies;
the method for obtaining the mechanical behavior of the non-continuous medium between the rock masses in the three-dimensional model according to the discrete element analysis method specifically comprises the following steps:
obtaining a contact body and a target body in the three-dimensional model through a non-binary tree search contact algorithm;
and taking the vector sum of the normal contact force and the tangential contact force as the contact force.
2. The method according to claim 1, wherein the step of computationally analyzing the stacking process and the final stacking shape of the three-dimensional model of the stacking blocks in the initial state according to the finite element/discrete element coupling analysis method FEM/DEM further comprises:
and establishing a three-dimensional model of the initial state of the rock stacking block, wherein the rock stacking block comprises a plurality of block stone bodies.
3. The method of claim 1, wherein the step of computationally analyzing the stacking process and the final stacking state of the three-dimensional model of the initial state of the stacking blocks, and the step of computationally analyzing the pouring process and the final filling state of the self-compacting concrete in the final stacking state of the stacking blocks based on the Smooth Particle Hydrodynamics (SPH) method further comprises:
and establishing a self-compacting concrete model, wherein the self-compacting concrete model is a Bingham rheological model.
4. The calculation method according to claim 2, wherein the establishing of the three-dimensional model of the stacking block initial state specifically comprises:
establishing a three-dimensional model of the rock stacking block in an initial state by using three-dimensional modeling software according to predetermined parameters of the rock stacking block; or,
and 3D scanning the outline of the rock mass body to form input data, inputting the input data into three-dimensional modeling software, and establishing a three-dimensional model of the rock mass in the initial state.
5. A computing system for numerical simulation of a rockfill concrete placement process, comprising:
the stacking analysis module is used for calculating and analyzing the stacking process and the final stacking form of the three-dimensional model of the stacking stone initial state according to a finite element/discrete element coupling analysis method FEM/DEM;
the pouring analysis module is used for calculating and analyzing the pouring process and the final filling state of the self-compacting concrete in the final stacking form of the stacking blocks based on a smooth particle fluid dynamics (SPH) method after calculating and analyzing the stacking process and the final stacking form of the three-dimensional model of the initial state of the stacking blocks;
iteratively updating the density, the position and the speed of fluid particles in the self-compacting concrete model one by one according to a time step based on an SPH (shortest path H) method;
analyzing the pouring process and the final filling state of the self-compacting concrete in the final stacking form of the stacking blocks based on the density, the position and the speed of the fluid particles;
wherein, the iteratively updating the density, the position and the speed of the fluid particles in the self-compacting concrete model one by one time step specifically comprises:
calculating the density change of the fluid particles one by one according to a continuity equation, capturing free surface particles and correcting the density of the free surface particles to obtain a density field of the self-compacting concrete model at each moment;
updating the motion track of each fluid particle through a fluid particle displacement equation one by one time step to obtain the position of any fluid particle at each moment;
determining the pressure of any fluid particle through a state equation one by one according to the time step, calculating the acceleration generated by the external force based on the pressure of each fluid particle, and calculating the velocity field of each moment in the self-compacting concrete model through a momentum equation based on the acceleration;
the method for calculating and analyzing the stacking process and the final stacking form of the three-dimensional model of the stacking stone initial state according to the finite element/discrete element coupling analysis method FEM/DEM specifically comprises the following steps:
obtaining the continuous medium mechanical behavior of each rock block in the three-dimensional model of the stacking initial state of the rock blocks according to a finite element analysis (FEM);
obtaining the mechanical behavior of a non-continuous medium among the rock masses in the three-dimensional model of the initial state of the rock pile according to a discrete element analysis (DEM);
substituting the continuous medium mechanical behaviors and the non-continuous medium mechanical behaviors of all the rock masses into a dynamic balance equation of the three-dimensional model of the initial state of the rock mass pile to obtain the stacking process and the final stacking form;
the continuous medium mechanical behavior comprises: an internal force vector caused by node deformation of the grid unit, an external force vector of the node of the grid unit and a mass matrix;
the method for obtaining the continuous medium mechanical behavior of each block stone body in the three-dimensional model according to the finite element analysis method specifically comprises the following steps:
carrying out mesh division on the three-dimensional model of the rock-fill block in the initial state to obtain a finite element model consisting of mesh units;
obtaining the internal force vector according to the grid unit volume, the standard type function and the Cauchy stress tensor;
obtaining the quality matrix according to the volume of the grid unit when the grid unit is not deformed, the density corresponding to the volume of the grid unit when the grid unit is not deformed and the standard function;
obtaining the external force vector according to the physical strength, the surface force and the unit surface area of the stone body;
the mechanical behavior of the non-continuous medium comprises: contact force between the stone bodies;
the method for obtaining the mechanical behavior of the non-continuous medium between the rock masses in the three-dimensional model according to the discrete element analysis method specifically comprises the following steps:
obtaining a contact body and a target body in the three-dimensional model through a non-binary tree search contact algorithm;
and taking the vector sum of the normal contact force and the tangential contact force as the contact force.
6. A computing device for numerical simulation of a rockfill concrete placement process, comprising:
at least one processor; and
at least one memory communicatively coupled to the processor, wherein:
the memory stores program instructions executable by the processor, the processor being capable of performing the computing method of any of claims 1 to 4 when invoked by the program instructions.
7. A non-transitory computer-readable storage medium storing computer instructions for causing a computer to perform the computing method of any one of claims 1 to 4.
CN201810699310.5A 2018-06-29 2018-06-29 Calculation method and system for numerical simulation of rockfill concrete pouring process Expired - Fee Related CN109033537B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810699310.5A CN109033537B (en) 2018-06-29 2018-06-29 Calculation method and system for numerical simulation of rockfill concrete pouring process

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810699310.5A CN109033537B (en) 2018-06-29 2018-06-29 Calculation method and system for numerical simulation of rockfill concrete pouring process

Publications (2)

Publication Number Publication Date
CN109033537A CN109033537A (en) 2018-12-18
CN109033537B true CN109033537B (en) 2020-12-01

Family

ID=65520981

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810699310.5A Expired - Fee Related CN109033537B (en) 2018-06-29 2018-06-29 Calculation method and system for numerical simulation of rockfill concrete pouring process

Country Status (1)

Country Link
CN (1) CN109033537B (en)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110750860B (en) * 2019-09-11 2023-05-05 四川轻化工大学 A Whole Process Analysis Method of Soil Slope and Landslide
CN111191373B (en) * 2020-01-03 2023-04-07 长沙有色冶金设计研究院有限公司 High-concentration tailing filling slurry stirring process optimization method
CN111666699B (en) * 2020-04-30 2023-06-02 山东大学 Cross-scale simulation calculation method of rock mass engineering based on REV full area coverage
CN111914321B (en) * 2020-06-09 2024-01-16 西安理工大学 Method for establishing three-phase microscopic model of rock-fill concrete
CN112052587B (en) * 2020-09-02 2023-12-01 中国人民解放军陆军工程大学 Densification method of three-dimensional meso-discrete model of sand cushion
CN112364553B (en) * 2021-01-13 2021-06-25 西南交通大学 Finite element-discrete element coupling evaluation method for surface erosion of high-speed rail ballastless track bed
CN113158413A (en) * 2021-02-05 2021-07-23 天津职业技术师范大学(中国职业培训指导教师进修中心) Diamond circular saw blade sawing system modeling and hard and brittle material sawing processing process numerical simulation method
CN113515878B (en) * 2021-07-07 2023-06-20 重庆交通大学 A Discrete Element 3D Modeling Method of Cumulative Body Based on Block Shape and Fragmentation
CN113338218B (en) * 2021-08-06 2021-10-26 西南交通大学 Multi-scale multi-medium comprehensive inversion method for debris flow flexible protection
CN119397875A (en) * 2024-12-31 2025-02-07 山东交通学院 A simulation method, device, equipment and storage medium for concrete pouring process

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106960070A (en) * 2016-12-28 2017-07-18 山东科技大学 A kind of seepage simulation method that coal body is reconstructed based on finite element discretization member CT

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8494827B2 (en) * 2009-09-25 2013-07-23 Exxonmobil Upstream Research Company Method of predicting natural fractures and damage in a subsurface region

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106960070A (en) * 2016-12-28 2017-07-18 山东科技大学 A kind of seepage simulation method that coal body is reconstructed based on finite element discretization member CT

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
光滑粒子流体动力学算法研究及其在流体力学中的数值模拟应用;王永鹏;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20111215;第2章第2.2节,第3章3.3节,第4章4.2.3节,4.3.2节 *
颗粒材料多尺度分析的连接尺度方法;李锡夔,万柯;《力学学报》;20101029;第42卷(第5期);第4页第2部分,第5页第3部分 *

Also Published As

Publication number Publication date
CN109033537A (en) 2018-12-18

Similar Documents

Publication Publication Date Title
CN109033537B (en) Calculation method and system for numerical simulation of rockfill concrete pouring process
CN108984829B (en) Calculation method and system for accumulation process of rockfill concrete rockfill
Zhang et al. Numerical simulation of a flow-like landslide using the particle finite element method
Jing A review of techniques, advances and outstanding issues in numerical modelling for rock mechanics and rock engineering
Onishi et al. A locking‐free selective smoothed finite element method using tetrahedral and triangular elements with adaptive mesh rezoning for large deformation problems
Dunant et al. A critical comparison of several numerical methods for computing effective properties of highly heterogeneous materials
CN111475937B (en) A simulation method of fluid-solid two-phase fluidized landslide
Jin et al. Two-phase PFEM with stable nodal integration for large deformation hydromechanical coupled geotechnical problems
WO2022041642A1 (en) Unsteady-state seepage simulation method and system for natural gas reservoir
CN106960070A (en) A kind of seepage simulation method that coal body is reconstructed based on finite element discretization member CT
CN109992830B (en) Simulation method of landslide disaster scene based on material point method
González et al. A natural element updated Lagrangian strategy for free-surface fluid dynamics
Zhang et al. Quasi-static collapse of two-dimensional granular columns: insight from continuum modelling
CN111125963A (en) Numerical simulation system and method based on Lagrange integral point finite element
CN117932730B (en) Cross-scale mechanical properties and stability measurement method of deep gravel pile composite foundation
CN108920768A (en) A kind of fluid-structure coupling computation for elastic thin-walled structures
CN105404758A (en) Numerical simulation method of solid continuum deformation based on finite element method
CN117350029A (en) Numerical simulation methods, devices, media and equipment considering the fracture of rock and soil masses
Liang et al. Extended material point method for the three‐dimensional crack problems
CN115099119A (en) Surge simulation method and system based on landslide mass motion model and gridless SPH
CN112818574B (en) Numerical method for simulating debris flow initiation, flow development and re-deposition
Xiong et al. Development of an unresolved CFD-DEM method for interaction simulations between large particles and fluids
CN118762782A (en) A discrete element-based in-situ ground stress equilibrium method for rock and soil analysis
CN114564899B (en) A simulation method, device and storage medium for sand seepage damage
CN116956406A (en) SPH analysis method and system for the whole process of landslide dynamic rupture movement

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
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20201201

Termination date: 20210629