CN113343510A - 二维随机裂隙网格生成方法 - Google Patents

二维随机裂隙网格生成方法 Download PDF

Info

Publication number
CN113343510A
CN113343510A CN202110490762.4A CN202110490762A CN113343510A CN 113343510 A CN113343510 A CN 113343510A CN 202110490762 A CN202110490762 A CN 202110490762A CN 113343510 A CN113343510 A CN 113343510A
Authority
CN
China
Prior art keywords
fracture
crack
grid
simulation
nodes
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.)
Granted
Application number
CN202110490762.4A
Other languages
English (en)
Other versions
CN113343510B (zh
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.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
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 Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN202110490762.4A priority Critical patent/CN113343510B/zh
Publication of CN113343510A publication Critical patent/CN113343510A/zh
Application granted granted Critical
Publication of CN113343510B publication Critical patent/CN113343510B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computing Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供一种二维随机裂隙网格生成方法，能够有效模拟真实裂隙形态，该方法包括：步骤1.根据模拟对象的实际情况，确定模拟区域的长和宽；步骤2.根据模拟对象的实际情况，确定裂隙宽度，并利用正三角形网格划分整个模拟区域，设置模拟区域长边上网格数量；步骤3.根据模拟对象的实际情况，确定最大迭代步数、以及裂隙条数；最大迭代步数应根据实际裂隙长度进行确定，每进行一次迭代计算，裂隙会向两端发展一个网格边长的长度，计算最大裂隙长度；步骤4.利用Matlab编程在模拟区域产生初始断裂点；步骤5.在初始断裂点的基础上进行迭代计算，设置裂隙向不同方向发展的概率模拟裂隙的发育过程；步骤6.生成裂隙模型。

Description

二维随机裂隙网格生成方法
技术领域
本发明属于材料随机裂隙生成技术领域，具体涉及一种二维随机裂隙网格生成方法。
技术背景
裂隙广泛存在于各种材料中，裂隙会显著降低材料的强度和刚度并增加其渗透性，可能会导致严重的工程安全问题。为了预测裂隙对于工程安全及破坏模式的影响，工程中常使用数值模拟方法进行设计，传统的确定性裂隙建模方法是在数值模型的关键位置设置几条直线或平面裂隙结构，但是该方法无法准确描述真实裂隙的形态和位置。离散元方法、CT扫描图像处理等方法是当前常用的随机裂隙生成方法，但是受限于计算能力、样品尺寸等问题，无法进行大尺度场地上随机裂隙的模拟。与此同时，由于裂隙会影响材料的失效模式，如果不考虑裂隙的存在或仅将材料强度进行整体折减，往往会低估裂隙对于结构体安全的影响。因此，如何模拟工程中真实裂隙形态并进行数值计算是当前亟待解决的问题。
另外，裂隙的发展不是完全随机的，由于裂隙尖端会产生应力集中现象，裂隙将沿裂隙尖端继续发展。并且因为材料性质或受力不均匀，裂隙在发育过程中会在一定范围内偏转，从而导致裂隙形态并非是直线。因此，人为假定的直线裂隙与真实裂隙形态相差较大。而解决这些问题的基础和关键是寻找一种高效的随机裂隙网络生成方法，生成的裂隙模型应当可以应用于复杂的数值计算中，可为工程设计提供一定的参考。
发明内容
本发明是为了解决上述问题而进行的，目的在于提供一种能够模拟真实裂隙形态，并具有较高计算效率的二维随机裂隙网格生成方法。
本发明为了实现上述目的，采用了以下方案：
如图1所示，本发明提供一种二维随机裂隙网格生成方法，其特征在于，包括以下步骤：
步骤1.根据模拟对象的实际情况，确定模拟区域的长La和宽Wa；
步骤2.根据模拟对象的实际情况，确定裂隙宽度wc，并利用正三角形网格划分整个模拟区域，设置模拟区域长边上网格数量ngb；
步骤3.根据模拟对象的实际情况，确定最大迭代步数nt、以及裂隙条数fb；
最大迭代步数应根据实际裂隙长度进行确定，每进行一次迭代计算，裂隙会向两端发展一个网格边长的长度，采用下式1计算最大裂隙长度lc,max：
lc，max＝(2nt+1)a (1)
a为正三角形网格的边长，计算公式为：
Figure BDA0003051933100000021
当裂隙发展至模拟区域边界处或其迭代步数达到nt时，停止继续发展；
步骤4.利用Matlab编程在模拟区域产生初始断裂点；包括以下子步骤：
步骤4-1.将网格所有节点的X,Y坐标分别存储进m×n的矩阵中，其中n为节点的最大列数，表示为
n＝ngb+l (3)
由于模拟区域通过正三角形网格划分，当划分网格总行数为奇数时，奇数行节点的列数为ngb，偶数行节点的列数为ngb+1；当划分网格总行数为偶数时，奇数行节点的列数为ngb+1，偶数行节点的列数为ngb；当节点列数为ngb时，则矩阵中该行的第n列储存为0；
m为节点的行数，表示为
Figure BDA0003051933100000022
步骤4-2.对模拟区域中所有节点进行编号，使用randperm函数在节点中随机抽取fb个随机数，为了防止裂隙超出模拟区域的范围，边界节点在抽样的过程中应被过滤掉，抽取到的节点作为裂隙的初始断裂点；如图2所示，在正三角形网格中，裂隙具有三种基本断裂模式：水平发展的平裂隙断裂模式，左上右下发展的左上右下斜裂隙断裂模式，右上左下发展的右上左下斜裂隙断裂模式；对于平裂隙，破裂点处的节点分裂成上下两个节点；对于斜裂隙，破裂点处的节点分裂为左右两个节点，每个节点的偏移量均为wc/2；
步骤5.在初始断裂点的基础上进行迭代计算，设置裂隙向不同方向发展的概率(P1、P2、P3)模拟裂隙的发育过程；
由于裂隙尖端处存在应力集中，垂直裂隙的方向上会产生应力释放，因此裂隙会沿尖端方向继续开展，如图3所示，在正三角形网格中裂隙只能是沿着三个方向发展：平行于裂隙尖端的方向和与该方向夹角为60°的另外两个方向进行发展；
步骤6.生成裂隙模型。
本发明提供的二维随机裂隙网格生成方法，还可以具有以下特征：模拟对象为场地或试样，在步骤1中，模拟对象应为矩形，该矩形的上下边长为La，左右边宽为Wa，该参数的选择应根据模拟对象的实际尺寸进行确定。
优选地，本发明提供的二维随机裂隙网格生成方法，还可以具有以下特征：模拟对象为场地或试样，在步骤1中，模拟对象应为矩形，该矩形的上下边长为La，左右边宽为Wa，该参数的选择应根据模拟对象的实际尺寸进行确定。
优选地，本发明提供的二维随机裂隙网格生成方法，还可以具有以下特征：在步骤2中，ngb的取值应基于裂隙宽度wc而确定，裂隙宽度wc计算公式为：
wc＝2δa (5)
式中，δ为裂隙宽度系数，取值范围为(0,1]；假设裂隙的宽度恒定不变，对应实际场地中裂隙的平均宽度；
模型长边上的网格数ngb的取值应保证裂隙宽度wc与模拟对象的平均裂隙宽度一致，ngb确定后，网格总数ng通过下式计算：
Figure BDA0003051933100000031
优选地，本发明提供的二维随机裂隙网格生成方法，还可以具有以下特征：在步骤3中，对于裂隙发展至模拟区域边界处或其迭代步数达到nt时停止继续发展的设定通过if语句进行实现：
ifx>＝La||x<＝0||y>＝Wa||y<＝0||nc>nt；
break(跳出循环)；
end；
式中，x,y分别为裂隙下一步发展的坐标；nc为当前的迭代步数；由于模拟的裂隙长度不会超过最大裂隙长度lc,max，故nt的取值可以模拟对象中最长裂隙的长度进行确定。
优选地，本发明提供的二维随机裂隙网格生成方法，还可以具有以下特征：在步骤3中，裂隙条数fb的选择应根据模拟对象裂隙的数量确定，如果未知模拟对象的裂隙数量，可以根据模拟的工况进行参数确定，若研究的是材料中的一条或多条主裂隙，例如边坡的软弱带，岩石节理，水力致裂等，建议fb取值小于10；若研究的是裂隙网络，例如土壤表层干缩裂隙形成的裂隙网，建议fb取值大于10，以便形成裂隙网络。
优选地，本发明提供的二维随机裂隙网格生成方法，还可以具有以下特征：在步骤5中，每条裂隙的发展是独立的，相互之间不会产生影响，在裂隙交界处，后发育到此处的裂隙不会对已存在的裂隙上的节点位置产生影响；该步骤主要通过设置虚拟网格和实际网格进行实现，其中虚拟网格与裂隙产生前的实际网格完全相同，但是在裂隙发展过程中虚拟网格的节点位置不会发生改变；通过使用if语句判断破裂点处的实际网格与虚拟网格节点坐标是否一致，如果一致则无裂隙交汇；如果不一致则发生裂隙交汇，此时两条裂隙重叠部分中，新产生裂隙导致的节点变化不会反映到实际网格中。
优选地，本发明提供的二维随机裂隙网格生成方法，还可以具有以下特征：在步骤5中，裂隙向不同方向发展的概率应根据实际情况进行设置，由于裂隙尖端存在的应力集中现象，裂隙沿裂隙尖端方向发展的概率P1应大于向另外两个方向发展的概率P2、P3，且应符合：
P1+P2+P3＝100％ (7)
通常情况下，P2＝P3；P1、P2、P3的具体取值应按照实际裂隙的曲直进行确定；在确定每个方向的发展概率后，使用unifrnd(0,1)函数实现特定概率下裂隙的随机发展。
优选地，本发明提供的二维随机裂隙网格生成方法，还可以具有以下特征：若裂隙为岩石节理类的较直裂隙，则P1建议大于60％；若裂隙为水力压裂、水泥裂缝类的较弯曲裂隙，则P1建议取值在40％至60％之间。
发明的作用与效果
(1)本发明所提供的二维随机裂隙生成方法能够模拟真实裂隙形态，并具有较高的计算效率，可以模拟不同材料在任意尺度下的裂隙生成过程，特别适用于大尺度下裂隙高效建模；
(2)本发明方法考虑了应力集中对裂隙发展方向的影响，能够获得最为接近真实裂隙形态的条带状或网状裂隙模型，准确模拟实际工程中的真实裂隙形态结构；
(3)本发明方法生成的裂隙模型经过简单处理后可以导入数值计算软件，可以根据实际工况进行复杂的数值计算；
(4)本发明方法采用的参数基本都可根据实际场地的实测、经验或统计规律方便获得，易于实施。
附图说明
图1为本发明涉及的二维随机裂隙网络生成方法的流程图；
图2为本发明涉及的裂隙开裂的三种基本模式示意图；
图3为本发明涉及的裂隙发展方向示意图；
图4为本发明实施例中涉及的模拟裂隙的结果图像。
具体实施方式
以下结合附图对本发明涉及的二维随机裂隙网格生成方法的具体实施方案进行详细地说明。
<实施例>
如图1所示，本实施例所提供的二维随机裂隙网格生成方法包括以下步骤：
步骤1：打开Matlab并添加变量，设置模拟区域的长La和宽Wa均为200，最大迭代步数nt分别为25、50、100，裂隙宽度系数δ为0.6，以及裂隙条数fb分别为10、20、30，设置模拟区域左下角的坐标为(0,0)；
步骤2：利用正三角形网格划分整个模拟区域，设置模拟区域长边上网格数量ngb为150，即网格边长a为4/3，共150列260行，网格数总计为3900个，进而可以得到正三角形网格每个节点的坐标，以矩阵的形式储存节点坐标；
步骤3：利用Matlab编程在模拟区域产生初始断裂点，使用randperm函数在所有网格线中随机寻找并基于图2所示三种基本断裂方式生成fb个初始断裂点，对节点坐标进行更新；
步骤4：在fb个初始断裂点的基础上进行迭代计算，模拟裂隙的发育过程。如图3所示，首先应根据初始断裂点及裂隙尖端坐标判断初始断裂点形态，再根据裂隙形态确定本迭代步中裂隙左右两侧可能发展的三个方向。本实施例中P1为60％，P2和P3均为20％。确定裂隙向各个方向发展的概率后，可以使用unifrnd(0,1)函数实现特定概率下的裂隙随机发展，每进行一次迭代计算，都应根据当前的裂隙形态对节点坐标进行更新。当迭代步数达到nt或裂隙发展至模拟区域边缘时，该条裂隙将停止发育。模拟结果如图4所示，与真实裂隙发展情况非常相符，表明了本方法能够准确地模拟出实际工程中的真实裂隙形态结构。
以上实施例仅仅是对本发明技术方案所做的举例说明。本发明所涉及的二维随机裂隙网格生成方法并不仅仅限定于在以上实施例中所描述的内容，而是以权利要求所限定的范围为准。本发明所属领域技术人员在该实施例的基础上所做的任何修改或补充或等效替换，都在本发明的权利要求所要求保护的范围内。

Claims (8)

1.一种二维随机裂隙网格生成方法，其特征在于，包括以下步骤：
步骤1.根据模拟对象的实际情况，确定模拟区域的长La和宽Wa；
步骤2.根据模拟对象的实际情况，确定裂隙宽度wc，并利用正三角形网格划分整个模拟区域，设置模拟区域长边上网格数量ngb；
步骤3.根据模拟对象的实际情况，确定最大迭代步数nt、以及裂隙条数fb；
最大迭代步数应根据实际裂隙长度进行确定，每进行一次迭代计算，裂隙会向两端发展一个网格边长的长度，采用下式1计算最大裂隙长度lc,max：
lc，max＝(2nt+1)a (1)
a为正三角形网格的边长，计算公式为：
Figure FDA0003051933090000011
当裂隙发展至模拟区域边界处或其迭代步数达到nt时，停止继续发展；
步骤4.利用Matlab编程在模拟区域产生初始断裂点；包括以下子步骤：
步骤4-1.将网格所有节点的X,Y坐标分别存储进m×n的矩阵中，其中n为节点的最大列数，表示为
n＝ngb+1 (3)
由于模拟区域通过正三角形网格划分，当划分网格总行数为奇数时，奇数行节点的列数为ngb，偶数行节点的列数为ngb+1；当划分网格总行数为偶数时，奇数行节点的列数为ngb+1，偶数行节点的列数为ngb；当节点列数为ngb时，则矩阵中该行的第n列储存为0；
m为节点的行数，表示为
Figure FDA0003051933090000012
步骤4-2.对模拟区域中所有节点进行编号，使用randperm函数在节点中随机抽取fb个随机数，边界节点在抽样的过程中应被过滤掉，抽取到的节点作为裂隙的初始断裂点；在正三角形网格中，裂隙具有三种基本断裂模式：水平发展的平裂隙断裂模式，左上右下发展的左上右下斜裂隙断裂模式，右上左下发展的右上左下斜裂隙断裂模式；对于平裂隙，破裂点处的节点分裂成上下两个节点；对于斜裂隙，破裂点处的节点分裂为左右两个节点，每个节点的偏移量均为wc/2；
步骤5.在初始断裂点的基础上进行迭代计算，设置裂隙向不同方向发展的概率模拟裂隙的发育过程；
由于裂隙尖端处存在应力集中，垂直裂隙的方向上会产生应力释放，因此裂隙会沿尖端方向继续开展，在正三角形网格中裂隙是沿着平行于裂隙尖端的方向和与该方向夹角为60°的另外两个方向进行发展；
步骤6.生成裂隙模型。
2.根据权利要求1所述的二维随机裂隙网格生成方法，其特征在于：
其中，模拟对象为场地或试样，在步骤1中，模拟对象应为矩形，该矩形的上下边长为La，左右边宽为Wa，该参数的选择应根据模拟对象的实际尺寸进行确定。
3.根据权利要求1所述的二维随机裂隙网格生成方法，其特征在于：
其中，在步骤2中，ngb的取值应基于裂隙宽度wc而确定，裂隙宽度wc计算公式为：
wc＝2δa (5)
式中，δ为裂隙宽度系数，取值范围为(0,1]；
模型长边上的网格数ngb的取值应保证裂隙宽度wc与模拟对象的平均裂隙宽度一致，ngb确定后，网格总数ng通过下式计算：
Figure FDA0003051933090000021
4.根据权利要求1所述的二维随机裂隙网格生成方法，其特征在于：
其中，在步骤3中，对于裂隙发展至模拟区域边界处或其迭代步数达到nt时停止继续发展的设定通过if语句进行实现：
ifx>＝La||x<＝0||y>＝Wa||y<＝0||nc>nt；
break(跳出循环)；
end；
式中，x,y分别为裂隙下一步发展的坐标；nc为当前的迭代步数；nt的取值依据模拟对象中最长裂隙的长度进行确定。
5.根据权利要求1所述的二维随机裂隙网格生成方法，其特征在于：
其中，在步骤3中，裂隙条数fb的选择应根据模拟对象裂隙的数量确定，如果未知模拟对象的裂隙数量，可以根据模拟的工况进行参数确定，若研究的是材料中的一条或多条主裂隙，建议fb取值小于10；若研究的是裂隙网络，建议fb取值大于10，以便形成裂隙网络。
6.根据权利要求1所述的二维随机裂隙网格生成方法，其特征在于：
其中，在步骤5中，每条裂隙的发展是独立的，相互之间不会产生影响，在裂隙交界处，后发育到此处的裂隙不会对已存在的裂隙上的节点位置产生影响；该步骤主要通过设置虚拟网格和实际网格进行实现，其中虚拟网格与裂隙产生前的实际网格完全相同，但是在裂隙发展过程中虚拟网格的节点位置不会发生改变；通过使用if语句判断破裂点处的实际网格与虚拟网格节点坐标是否一致，如果一致则无裂隙交汇；如果不一致则发生裂隙交汇，此时两条裂隙重叠部分中，新产生裂隙导致的节点变化不会反映到实际网格中。
7.根据权利要求1所述的二维随机裂隙网格生成方法，其特征在于：
其中，在步骤5中，裂隙向不同方向发展的概率应根据实际情况进行设置，由于裂隙尖端存在的应力集中现象，裂隙沿裂隙尖端方向发展的概率P1应大于向另外两个方向发展的概率P2、P3，且应符合：
P1+P2+P3＝100％ (7)
通常情况下，P2＝P3；P1、P2、P3的具体取值应按照实际裂隙的曲直进行确定；在确定每个方向的发展概率后，使用unifrnd(0,1)函数实现特定概率下裂隙的随机发展。
8.根据权利要求1所述的二维随机裂隙网格生成方法，其特征在于：
其中，若裂隙为岩石节理类的较直裂隙，则P1建议大于60％；若裂隙为水力压裂、水泥裂缝类的较弯曲裂隙，则P1建议取值在40％至60％之间。
CN202110490762.4A 2021-05-06 2021-05-06 二维随机裂隙网格生成方法 Active CN113343510B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110490762.4A CN113343510B (zh) 2021-05-06 2021-05-06 二维随机裂隙网格生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110490762.4A CN113343510B (zh) 2021-05-06 2021-05-06 二维随机裂隙网格生成方法

Publications (2)

Publication Number Publication Date
CN113343510A true CN113343510A (zh) 2021-09-03
CN113343510B CN113343510B (zh) 2022-05-17

Family

ID=77469595

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110490762.4A Active CN113343510B (zh) 2021-05-06 2021-05-06 二维随机裂隙网格生成方法

Country Status (1)

Country Link
CN (1) CN113343510B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113946984A (zh) * 2021-12-20 2022-01-18 北京科技大学 一种三维随机裂隙模型建立方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100312529A1 (en) * 2009-06-05 2010-12-09 Schlumberger Technology Corporation Fracture network characterization method
CN105787220A (zh) * 2016-04-22 2016-07-20 山东科技大学 一种煤层高压注水致裂-渗流数值模拟方法
CN108399290A (zh) * 2018-02-09 2018-08-14 河海大学 基于裂隙网络连通性的裂隙流的modflow模拟方法
WO2019238451A1 (en) * 2018-06-13 2019-12-19 Danmarks Tekniske Universitet A method and a system for modelling and simulating a fractured geological structure
CN110851972A (zh) * 2019-11-04 2020-02-28 武汉大学 基于高斯随机场的岩土体结构随机裂隙模拟方法及装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100312529A1 (en) * 2009-06-05 2010-12-09 Schlumberger Technology Corporation Fracture network characterization method
CN105787220A (zh) * 2016-04-22 2016-07-20 山东科技大学 一种煤层高压注水致裂-渗流数值模拟方法
CN108399290A (zh) * 2018-02-09 2018-08-14 河海大学 基于裂隙网络连通性的裂隙流的modflow模拟方法
WO2019238451A1 (en) * 2018-06-13 2019-12-19 Danmarks Tekniske Universitet A method and a system for modelling and simulating a fractured geological structure
CN110851972A (zh) * 2019-11-04 2020-02-28 武汉大学 基于高斯随机场的岩土体结构随机裂隙模拟方法及装置

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113946984A (zh) * 2021-12-20 2022-01-18 北京科技大学 一种三维随机裂隙模型建立方法

Also Published As

Publication number Publication date
CN113343510B (zh) 2022-05-17

Similar Documents

Publication Publication Date Title
Chen et al. The enhanced extended finite element method for the propagation of complex branched cracks
CN109885864B (zh) 一种三维钢桥塔涡激振动计算方法
CN108399290B (zh) 基于裂隙网络连通性的裂隙流的modflow模拟方法
CN104574472B (zh) 基于嵌入网格的固体碎裂模拟和动画方法
CN112699623A (zh) 基于非结构网格规则化重构技术的高精度热流计算方法
CN109165404B (zh) 一种扩展光滑无网格伽辽金法
CN113343510B (zh) 二维随机裂隙网格生成方法
CN111460568B (zh) 一种混凝土重力坝运行期裂缝扩展判别方法
CN110162849B (zh) 一种混杂纤维混凝土的建模方法
CN108829954A (zh) 一种基于结构模型剖面特征的钢筋配置方法
Sun et al. Research of large scale mechanical structure crack growth method based on finite element parametric submodel
CN113343423B (zh) 基于强度空间变异性的随机裂隙网络生成方法
CN113158425B (zh) 一种基于边界元方法的全三维裂缝相交过程模拟方法
CN106844963A (zh) 模拟开挖至运行全过程的拱坝三维网格模型自动剖分方法
Saucedo-Mora et al. Method for the explicit insertion of microstructure in Cellular Automata Finite Element (CAFE) models based on an irregular tetrahedral Finite Element mesh: Application in a multi-scale Finite Element Microstructure MEshfree framework (FEMME)
CN111159794A (zh) 多裂隙类岩石试样力学性质的几何损伤流变分析方法
CN104484527B (zh) 一种离散结构拓扑优化过程中均布载荷自动动态修改方法
CN109241579B (zh) 一种考虑流体附加质量多维空间瞬变特性的快速建模方法
CN101354710B (zh) 一种线段求交的方法及装置
CN111222215A (zh) 节理岩体力学性质的几何损伤流变模型分析方法
CN106934729A (zh) 建筑检测鉴定方法及装置
Yu et al. An h-adaptive numerical manifold method for solid mechanics problems
Yuan et al. Improved random aggregate model for numerical simulations of concrete engineering simulations of concrete engineering
CN106874611A (zh) 一种基于超体积迭代策略的含区间参数结构响应区间的分析方法
CN107577899B (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
GR01 Patent grant
GR01 Patent grant