CN113887094A - 运用多尺度离散实体元解决连续体结构力学问题的仿真方法 - Google Patents
运用多尺度离散实体元解决连续体结构力学问题的仿真方法 Download PDFInfo
- Publication number
- CN113887094A CN113887094A CN202111037270.6A CN202111037270A CN113887094A CN 113887094 A CN113887094 A CN 113887094A CN 202111037270 A CN202111037270 A CN 202111037270A CN 113887094 A CN113887094 A CN 113887094A
- Authority
- CN
- China
- Prior art keywords
- ball
- spring
- area
- elements
- small
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000004088 simulation Methods 0.000 title claims abstract description 27
- 238000004364 calculation method Methods 0.000 claims abstract description 28
- 230000007704 transition Effects 0.000 claims abstract description 26
- 239000007787 solid Substances 0.000 claims abstract description 20
- 239000002245 particle Substances 0.000 claims description 16
- 238000006073 displacement reaction Methods 0.000 claims description 15
- 239000000463 material Substances 0.000 claims description 11
- 230000001133 acceleration Effects 0.000 claims description 8
- 238000013016 damping Methods 0.000 claims description 8
- 239000008188 pellet Substances 0.000 claims description 3
- OIGNJSKKLXVSLS-VWUMJDOOSA-N prednisolone Chemical compound O=C1C=C[C@]2(C)[C@H]3[C@@H](O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 OIGNJSKKLXVSLS-VWUMJDOOSA-N 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 4
- 238000011439 discrete element method Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000010008 shearing Methods 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000000712 assembly Effects 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000011089 mechanical engineering Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- 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)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种多尺度离散实体元解决连续体结构力学问题的仿真方法。多尺度离散实体元模型包括大球区域、小球区域和过渡区域,球元之间通过弹簧系统连接,弹簧系统包括轴向弹簧和两个方向的切向弹簧。基于连续介质力学理论以能量等效为原则,推导了细观模型参数(弹簧刚度)和宏观弹性常数(弹性模量和泊松比)之间的解析关系。本发明可以有效减少球元和弹簧的数量,提高计算效率,使离散实体元法提供具有更多可能性,适用于结构或构件的大变形、损伤断裂和倒塌破坏等复杂力学仿真。
Description
技术领域
本发明涉及一种多尺度离散实体单元法力学仿真方法,用于解决结构或构件的大变形、损伤断裂和倒塌破坏等强非线性力学问题,属于离散元分析技术领域。
背景技术
离散单元法是美国学者Cundall于1971年提出来的一种非连续性数值计算方法[1-2],最初被用于分析岩石边坡运动的分析计算。([1]Cundall P A,Strack O D L.Adiscrete numerical model for granular assemblies[J].geotechnique,1979,29(1):47-65.[2]Cundall P A,Hart R D.Numerical modelling of discontinua[J].Engineering computations,1992,9(2):101-113.)离散单元法的基本原理为:把研究对象划分成一个个离散的单元,在结构受力变形、运动过程中,单元可以分离,即一个单元与其邻近单元可以接触,也可以分开[3-4]。([3]邢继波,王泳嘉.离散元法的改进及其在颗粒介质研究中的应用[J].岩土工程学报,1990,12(5):51-57.[4]王泳嘉,邢纪波.离散单元法及其在岩土力学中的应用.东北工学院出版社[M].1991.)传统离散单元法主要研究对象为散粒体,对连续体结构不适用。近年来开始有学者尝试将离散单元法用于钢结构倒塌、稳定、冲击破坏问题的数值计算。叶继红课题组[5-7]提出了杆系结构的离散元模型,推导了模型中弹簧的刚度系数,建立了适用于杆系结构的离散单元法塑性铰模型、考虑截面塑性开展的纤维模型、断裂准则和单元碰撞模型。([5]齐念,叶继红.基于离散元法的杆系结构几何非线性大变形分析[J].东南大学学报(自然科学版),2013(5):917-922.[6]齐念.DEM/FEM耦合计算方法研究及其在网壳倒塌破坏模拟中的应用[D].南京:东南大学,2016.[7]覃亚男.基于离散单元法的简单结构静动力响应数值模拟研究[D].南京:东南大学,2015.)。
发明人所在的冯若强课题组[8-10]提出了离散实体单元法,改进了球元间接触弹簧的设置规则,在球元立方体排列模型的棱边和面对角线上采用弹簧进行连接,球元间弹簧包括法向弹簧和两个方向切线弹簧,使其能够准确地体现连续体结构的力学性能,同时可以考虑材料的泊松比效应。([8]Zhu B,Feng R,Wang X.3D discrete solid-elementmethod for elastoplastic problems of continuity[J].Journal of engineeringmechanics,2018,144(7):04018051.[9]Zhu B,Feng R.Discrete solid element modelapplied to plasticity and dynamic crack propagation in continuous medium[J].Computational Particle Mechanics,2019:1-17.[10]Zhu B,Feng R.Investigation ofa boundary simulation of continuity using the discrete solid element method[J].Advances in Mechanical Engineering,2019,11(1):1687814018822397.)
综上所述,现有的离散实体单元法已可以对简单的连续体结构进行力学分析,但在解决受力复杂的复杂结构还存在球元数目大、计算效率低等问题。因此,迫切需要多尺度离散元计算方法以高效模拟受力复杂的复杂结构力学问题。
发明内容
发明目的:针对传统离散单元法不能很好模拟连续介质,现有的离散实体单元法已可以对简单的连续体结构进行力学分析,但在解决受力复杂的复杂结构还存在球元数目大、计算效率低等问题。本发明提供了一种多尺度离散实体单元法,基于能量等效原则,构建了多尺度离散实体元模型,推导了细观模型参数(弹簧刚度)与宏观弹性常数(弹性模量和泊松比)之间的解析关系,从而解决受力复杂的复杂结构采用传统离散实体元法计算中球元数目大、计算效率低等问题。
技术方案:为了实现上述发明目的,本发明提出的多尺度离散实体单元法力学仿真方法,步骤如下:
步骤1,根据仿真需求,将结构或构件离散为多尺度离散元模型,建立模型的球元为规则排列,材料的变形完全由球元间的接触弹簧来存储和表示;
多尺度离散实体元模型包括大球区域、小球区域和过渡区域,球元之间通过弹簧系统连接,弹簧系统包括轴向弹簧和两个方向的切向弹簧。
步骤2,t=0,t≥0时刻,对各颗粒单元的外力、内力、位移、速度和加速度赋初值。
步骤3,对所有颗粒单元:计算t时刻的外力及阻尼力,应用牛顿第二定律求解运动控制方程,得到t+Δt时刻的颗粒单元的位置与速度。
步骤5,根据弹塑性接触本构方程,计算Δt时步内颗粒单元的接触力增量ΔF,并更新当前时刻的颗粒单元间的接触力。
步骤6,根据平衡方程,计算各颗粒单元的内力、阻尼力以及外力,为t+2Δt时刻的计算准备,如此循环往复,直至模型达到稳定状态。
具体来说,包括以下具体步骤:
(1)初步估计仿真对象受力复杂的区域
根据力学概念、经验估计仿真对象受力复杂的区域,确定大球、小球球元的半径,将受力复杂区域划分为小球区域,将受力较简单的区域划分为大球区域,小球区域和大球区域之间为过渡区域。
(2)将仿真对象离散为多尺度离散元模型
根据步骤1确定的大球区域、小球区域和过渡区域,建立仿真对象的多尺度离散元模型,大球区域、小球区域和过渡区域球元按各自区域的规律排列,材料的变形完全由单元的接触弹簧来存储和表示,输入模型参数以及外荷载信息,确定边界条件。
大球区域弹簧刚度可表示为:
其中,为大球区域棱弹簧法向弹簧刚度,为大球区域面对角弹簧法向弹簧刚度, 分别为大球区域棱弹簧两个方向切向弹簧刚度,分别为大球区域面对角弹簧两个方向切向弹簧刚度,R为大球半径,E为弹性模量,μ为材料的泊松比。
小球区域弹簧刚度可表示为:
其中,为小球区域棱弹簧法向弹簧刚度,为小球区域面对角弹簧法向弹簧刚度, 分别为小球区域棱弹簧两个方向切向弹簧刚度,分别为小球区域面对角弹簧两个方向切向弹簧刚度,r为小球半径,E为弹性模量,μ为材料的泊松比。
过渡区域弹簧刚度可表示为:
其中,分别为过渡区域与大球区域、小球区域相连部分棱弹簧轴向刚度;分别为过渡区域与大球区域、小球区域相连部分面对角弹簧轴向刚度;分别为过渡区域棱弹簧和面对角弹簧轴向刚度;分别为过渡区域与大球区域、小球区域相连部分棱弹簧两个方向切向刚度;分别为过渡区域与大球区域、小球区域相连部分面对角弹簧两个方向切向刚度;分别为过渡区域棱弹簧、面对角弹簧两个方向切向刚度;R为小球半径,r为小球半径,E为弹性模量,μ为材料的泊松比。
(3)t=0时刻,对球元的速度、加速度等赋初值,对弹簧系统的内力赋初值;
(4)对所有球元,由t时刻的外力及阻尼力,应用牛顿第二定律,求得t+Δt时刻球元的位置与速度;
在公式(1)的基础上,运用中心有限差分算法确定球元球心的加速度,进而确定球元球心的位移。计算球元i在t+△t时刻的位移如公式(2)所示。
(5)计算时步Δt内,弹簧法向和切向接触力增量与法向和切向位移增量的关系如公式(3)所示。采用坐标转换矩阵将弹簧的法向和切向接触力增量转换为整体坐标下接触力增量,如公式(4)所示。采用增量的形式计算球元间的接触力,t+Δt时刻与t时刻整体坐标下弹簧接触力的关系如公式(5)所示。
(6)实体单元法计算模型中,一个球元与若干个球元通过弹簧连接,弹簧接触力和外力之和即为该球元在t+Δt时刻的不平衡力,这个不平衡力将用于求解下一计算时步t+2Δt球元的运动计算,进而求得球元的加速度、速度和位移。
有益效果:
(1)本发明基于连续介质力学理论以能量等效为原则,得到细观模型参数(弹簧刚度)和宏观弹性常数(弹性模量和泊松比)之间的解析关系。有效减少球元和弹簧的数量,提高计算效率,使离散实体元法提供具有更多可能性,适用于结构或构件的大变形、损伤断裂和倒塌破坏等复杂力学仿真。
(2)在多尺度离散实体单元法中,无需组集刚度矩阵,进行迭代求解,计算流程简单清晰。结构越复杂,单元数量就越多,相应的计算量也越大,但只是循环重复计算,不会给计算带来实质性困难。
(3)传统方法相比,在大型结构和复杂力学行为的模拟中离散实体单元法将表现出更大的优势。构建的三维离散实体元模型能够体现泊松比效应,实现不同尺度球元的连接。
附图说明
图1为本发明的流程图;
图2为多尺度离散实体单元法物理模型示意;
1-大球区域面弹簧、2-大球区域棱弹簧、3-过渡区域棱弹簧、4-小球区域面弹簧、5-小球区域棱弹簧、6-过渡区域棱弹簧、7-过渡区域面弹簧;
A-大球区域、B-过渡区域、C-小球区域;
a-端部固定约束、b-端部大球固定约束;
图3为多尺度离散实体单元法物理模型分区示意;
图4为离散实体元模型;
图5为多尺度离散实体元模型;
图6为离散实体元模型变形图;
图7为多尺度离散实体元模型变形图;
图8为离散实体元(DSEM)、多尺度离散实体元(MDSEM)与有限元模型(FEM)变形对比图。
具体实施方式
以下将结合实施例具体说明本发明的技术方案:
参照图1所示,本发明运用多尺度离散实体元解决连续体结构力学问题的仿真方法,步骤如下:
(1)初步估计仿真对象受力复杂的区域,根据力学概念、经验估计仿真对象受力复杂的区域,确定大球、小球球元的半径,将受力复杂区域划分为小球区域,将受力较简单的区域划分为大球区域,小球区域和大球区域之间为过渡区域。
(2)将仿真对象离散为多尺度离散元模型,建立仿真对象的多尺度离散元模型,大球区域、小球区域和过渡区域球元按各自区域的规律排列,材料的变形完全由单元的接触弹簧来存储和表示,输入模型参数以及外荷载信息,确定边界条件,参照图2、3所示。
(3)t=0时刻,对球元的速度、加速度等赋初值,对弹簧系统的内力赋初值;
(4)对所有球元,由t时刻的外力及阻尼力,应用牛顿第二定律,求得t+Δt时刻球元的位置与速度;
(5)计算时步Δt内法向和切向位移增量,采用坐标转换矩阵将弹簧的法向和切向接触力增量转换为整体坐标下接触力增量,计算球元间的接触力。
(6)实体单元法计算模型中,一个球元与若干个球元通过弹簧连接,弹簧接触力和外力之和即为该球元在t+Δt时刻的不平衡力,这个不平衡力将用于求解下一计算时步t+2Δt球元的运动计算,进而求得球元的加速度、速度和位移。
实施例1
下面用算例验证多尺度离散实体元在模拟构件大变形的优越性。
悬臂梁长160mm,截面尺寸为30mm×30mm,弹性模量为2.06×105MPa,泊松比为0.2,密度7850kg/m3,梁端部施加剪切力1000MPa。
离散实体元模型球元半径为2.5mm,球元数量1617,连接弹簧数量12092×3,如图4所示。多尺度离散实体元模型大球区域长80mm,大球直径10mm,小球区域长70mm,小球直径5mm,过渡区域长10mm,球元数量879,连接弹簧数量6332×3,如图5所示。本实施例考察了悬臂梁端部在弯矩M作用下Y向位移的位移变化情况。
图6和图7为离散实体元模型和多尺度离散实体元模型剪力作用下变形。离散实体元、多尺度离散实体元与有限元计算结果对比如表1所示。
图8为离散实体元、多尺度离散实体元与有限元模型变形对比图。比较离散实体元模型,多尺度离散实体元模型的球元数量和连接弹簧数量下降约50%,计算效率提高50%,但计算精度下降约2%,表明多尺度离散实体元在满足精度的前提下,计算效率有显著提高。
表1悬臂梁受剪力(y=0,z=0)直线上各点Y向位移计算结果
Claims (5)
1.运用多尺度离散实体元解决连续体结构力学问题的仿真方法,其特征在于,步骤如下:
步骤1,根据仿真需求,将结构或构件离散为多尺度离散元模型,建立模型的球元为规则排列,材料的变形完全由球元间的接触弹簧来存储和表示;
步骤2,t=0,t≥0时刻,对各颗粒单元的外力、内力、位移、速度和加速度赋初值;
步骤3,对所有颗粒单元:计算t时刻的外力及阻尼力,应用牛顿第二定律求解运动控制方程,得到t+Δt时刻的颗粒单元的位置与速度;
步骤5,根据弹塑性接触本构方程,计算Δt时步内颗粒单元的接触力增量ΔF,并更新当前时刻的颗粒单元间的接触力;
步骤6,根据平衡方程,计算各颗粒单元的内力、阻尼力以及外力,为t+2Δt时刻的计算准备,如此循环往复,直至模型达到稳定状态。
2.根据权利要求1所述的运用多尺度离散实体元解决连续体结构力学问题的仿真方法,其特征在于,步骤1所述的多尺度离散实体元模型包括大球区域、小球区域和过渡区域,球元之间通过弹簧系统连接,弹簧系统包括轴向弹簧和两个方向的切向弹簧;
大球球元半径R,小球球元半径r,大球、小球半径之比为2:1;
过渡区域包括大球、小球及连接大球和小球的弹簧系统,过渡区域的一面球元排列与大球区域球元排列相同,大球与相邻大球球心间距为2R,该面与大球区域相连;
过渡区域的一面与小球区域排列相同,小球与相邻小球球心间距为2r,该面与小球区域相连;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111037270.6A CN113887094A (zh) | 2021-09-06 | 2021-09-06 | 运用多尺度离散实体元解决连续体结构力学问题的仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111037270.6A CN113887094A (zh) | 2021-09-06 | 2021-09-06 | 运用多尺度离散实体元解决连续体结构力学问题的仿真方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113887094A true CN113887094A (zh) | 2022-01-04 |
Family
ID=79008198
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111037270.6A Pending CN113887094A (zh) | 2021-09-06 | 2021-09-06 | 运用多尺度离散实体元解决连续体结构力学问题的仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113887094A (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107391788A (zh) * | 2017-06-09 | 2017-11-24 | 东南大学 | 运用三维离散实体解决连续介质构件非线性力学问题的方法 |
CN112329290A (zh) * | 2020-10-28 | 2021-02-05 | 河海大学 | 可用于施工过程模拟的有限元离散元多尺度耦合计算方法 |
-
2021
- 2021-09-06 CN CN202111037270.6A patent/CN113887094A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107391788A (zh) * | 2017-06-09 | 2017-11-24 | 东南大学 | 运用三维离散实体解决连续介质构件非线性力学问题的方法 |
CN112329290A (zh) * | 2020-10-28 | 2021-02-05 | 河海大学 | 可用于施工过程模拟的有限元离散元多尺度耦合计算方法 |
Non-Patent Citations (1)
Title |
---|
冯若强;朱宝琛;胡椿昌;: "基于离散单元法的单层网壳结构冲击破坏数值模拟研究", 土木工程学报, no. 05, 15 May 2019 (2019-05-15) * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Analysis of three-dimensional curved beams using isogeometric approach | |
Luo et al. | Model order reduction for dynamic simulation of a flexible multibody system via absolute nodal coordinate formulation | |
CN108986220B (zh) | 一种加速有限元求解物体网格模型弹性变形的方法 | |
CN107391788B (zh) | 运用三维离散实体解决连续介质构件非线性力学问题的方法 | |
CN112949065B (zh) | 模拟层状岩体力学行为的双尺度方法、装置、存储介质及设备 | |
CN113361156B (zh) | 一种基于sph方法的飞机水箱投汲水数值模拟方法 | |
CN108984829B (zh) | 堆石混凝土堆石体堆积过程的计算方法和系统 | |
Sevim et al. | Water length and height effects on the earthquake behavior of arch dam-reservoir-foundation systems | |
CN106021644A (zh) | 确定混合维模型界面约束方程系数的方法 | |
CN109740182A (zh) | 一种基于再生核粒子的无网格物理变形仿真方法 | |
CN109657408A (zh) | 一种再生核粒子算法实现结构线性静力学仿真方法 | |
CN104281730A (zh) | 一种大转动变形的板壳结构动响应的有限元分析方法 | |
Yu et al. | Integration of geometry and analysis for vehicle system applications: continuum-based leaf spring and tire assembly | |
Wang et al. | A GPU-based parallel algorithm for 2D large deformation contact problems using the finite particle method | |
Stricklin et al. | Computation and solution procedures for non-linear analysis by combined finite element—finite difference methods | |
Hu et al. | Dynamic fracture simulation of flexible multibody systems via coupled finite elements of ANCF and particles of SPH | |
Park et al. | Multidisciplinary design optimization of a structurally nonlinear aircraft wing via parametric modeling | |
CN113887094A (zh) | 运用多尺度离散实体元解决连续体结构力学问题的仿真方法 | |
McDaniel et al. | Efficient mesh deformation for computational stability and control analyses on unstructured viscous meshes | |
CN114692446A (zh) | 一种矢量转动球坐标参数化的非线性壳有限元建模方法 | |
Guo et al. | Contact model and elastic deformation analysis of plate structure based on the discrete element method | |
Raju et al. | Nonlinear structural analysis using integrated force method | |
CN103745038A (zh) | 摩天轮结构有限元模型的建模方法 | |
CN114638143A (zh) | 一种适用于模拟水弹性问题的耦合数值计算方法 | |
Ozkul | A finite element formulation for dynamic analysis of shells of general shape by using the Wilson-θ method |
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 |