CN109657408A - 一种再生核粒子算法实现结构线性静力学仿真方法 - Google Patents
一种再生核粒子算法实现结构线性静力学仿真方法 Download PDFInfo
- Publication number
- CN109657408A CN109657408A CN201910023402.6A CN201910023402A CN109657408A CN 109657408 A CN109657408 A CN 109657408A CN 201910023402 A CN201910023402 A CN 201910023402A CN 109657408 A CN109657408 A CN 109657408A
- Authority
- CN
- China
- Prior art keywords
- displacement
- function
- particle
- equation
- discrete
- 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
- 239000002245 particle Substances 0.000 title claims abstract description 123
- 238000000034 method Methods 0.000 title claims abstract description 63
- 230000008929 regeneration Effects 0.000 title claims abstract description 43
- 238000011069 regeneration method Methods 0.000 title claims abstract description 43
- 230000003068 static effect Effects 0.000 title claims abstract description 35
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 31
- 238000004088 simulation Methods 0.000 title claims abstract description 27
- 238000006073 displacement reaction Methods 0.000 claims abstract description 80
- 239000011159 matrix material Substances 0.000 claims abstract description 33
- 238000012937 correction Methods 0.000 claims abstract description 10
- 238000010276 construction Methods 0.000 claims abstract description 7
- 238000012545 processing Methods 0.000 claims abstract description 7
- 230000009466 transformation Effects 0.000 claims abstract description 7
- 230000001133 acceleration Effects 0.000 claims description 12
- 238000013016 damping Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 7
- 238000013461 design Methods 0.000 abstract description 2
- 238000005457 optimization Methods 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 5
- 238000004590 computer program Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000002203 pretreatment Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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]
Abstract
本发明公开了一种再生核粒子算法实现结构线性静力学仿真方法，包括以下步骤：按照粒子密度对仿真几何模型进行粒子离散，得到离散颗粒；选用修正函数和窗函数构造形函数；根据影响域内离散颗粒的位移和形函数构造位移场函数；利用边界转化方法进行边界处理；联立动量方程和几何方程、物理方程、边界条件构造等效积分弱形式方程；将全场位移场函数带入等效积分弱形式方程求解，得到等效内力和外力矩阵；利用中心差分法进行代数求解，得到不同时间步的形变和应力。通过对比有限元方法，该方法能够大大简化复杂模型的前处理过程，提升了结构线性静力学仿真的求解精度和求解效率，有利于工程实际问题结构的设计和优化。
Description
技术领域
本发明涉及计算机辅助工程中的仿真计算领域，尤其涉及一种再生核粒子算法实现结构线性静力学仿真方法。
背景技术
结构静力学仿真是利用计算机程序对工程中的实际结构进行力学性能计算的过程，利用求解得到的结果可以对结构进行强度校核和刚度校核。当在其实际工况作用下，结构的位移与载荷的关系成线性变化时，则被认为是线性问题。目前，用于结构线性静力学仿真的计算机程序都是基于有限元法。有限元法在早期是以变分原理为基础发展起来的，所以广泛地应用于以拉普拉斯方程和泊松方程所描述的各类物理场中。其基本思想是：将无限自由度的结构近似离散成有限个自由度，通过构建自由度之间的关系函数，求解泛函的极值。
尽管有限元法在结构线性静力学仿真中的应用极其广泛，而且也非常成熟，但是由其自身算法带来的局限性，技术人员一直也在需求突破。一方面，有限元法的前提是利用单元对结构进行离散，即在程序求解前，需要生成单元，对于复杂的几何模型，这个过程需要花费工程师大量的时间和精力，很大程度上限制了工程师对结构设计方案进行快速验证迭代；另一方面，有限元法对网格尺寸和质量特别敏感和依赖，结构线性静力学仿真的结果精度很大程度依赖于结构网格处理的好坏，而求解效率又会随着网格尺寸的细化而大幅降低。
近年来，再生核粒子算法作为先进的无网格算法被引入到结构线性静力学仿真计算中，其具有基于网格的有限元法所不具备的优点。其在模型前处理时，无需如同有限元法那样划分网格，因此可以对任意复杂的几何模型进行快速前处理，可以大大节省复杂几何模型进行前处理的时间。没有网格尺寸和质量的牵制，结构线性静力学仿真的求解精度和效率也到了最大程度的保障。基于这么多的核心优势，目前再生核粒子算法成为了最前沿的研究方向，受到技术人员的广泛关注。
发明内容
本发明的目的是提供一种再生核粒子算法实现结构线性静力学仿真方法，可以克服传统方法中对网格的依赖，简化了复杂模型的前处理过程，提升了结构线性静力学仿真的求解精度和求解效率。
为实现上述目的，本发明采用以下技术方案：
一种再生核粒子算法实现结构线性静力学仿真方法，包括如下步骤：
(1)按照给定的粒子密度对仿真几何模型进行粒子离散，得到离散颗粒，确定所述离散颗粒的初始位置、影响域Ω和影响域半径ρ；
(2)基于再生核粒子算法，选用修正函数和窗函数的类型，构造所述影响域内所述离散颗粒的形函数；
(3)根据所述离散颗粒的位置和所述影响域半径ρ，利用所述影响域Ω内的所述离散颗粒的位移和所述形函数构造位移场函数；
(4)利用边界转化方法进行边界处理，消除边界不一致性引起的误差；
(5)联立动量方程和几何方程、物理方程、边界条件构造等效积分弱形式方程；
(6)构造位移矩阵，将所述位移场函数带入所述等效积分弱形式方程求解，得到等效内力和外力矩阵；
(7)利用中心差分法对等效积分弱形式方程进行代数求解，得到不同时间步的形变和应力。
优选地，步骤(1)中，所述离散颗粒的几何形状包括、但不限于圆形、球形、矩形、立方体。
优选地，步骤(1)中，所述影响域Ω由每个所述离散颗粒指定。
优选地，步骤(1)中，所述影响域半径ρ对于所有所述离散颗粒为恒定值或者对于每个单独的所述离散颗粒为不同的。
优选地，步骤(2)中，基于再生核粒子算法，所述形函数Φρ(x,x-y)为：
Φρ(x；x-y)＝Cρ(x；x-y)kρ(x-y)
式中，Cρ(x,x-y)为满足再生条件的修正函数，kρ(x-y)为满足再生条件的窗函数；
所述修正函数Cρ(x,x-y)是由多项式基函数的线性组合表示而成，如下：
Cρ(x；x-y)＝C0(x)+C1(x)(x-y)+C2(x)(x-y)2+…+CN(x)(x-y)N＝HT(x-y)C(x)
其中，N是满足再生条件的次数，HT(x-y)是N次多项式基函数的向量；C(x)为系数向量，通过施加N次再生条件来确定；Ci(x)为系数(i＝1,2,3……N)；
HT(x-y)＝[1,x-y,(x-y)2,…,(x-y)N]
C(x)＝[C0(x),C1(x),…,CN(x)]
所述窗函数κρ(x-y)是二维窗函数：
其中，为中间变量。
更优选地，所述窗函数选用三次样条函数，如下：
式中，
更优选地，所述窗函数包括、但不限于三次样条函数、B样条函数、高斯函数。
优选地，步骤(3)中，所述位移场函数uR(x)为：
式中，u(y)表示所述影响域Ω内的所述离散颗粒的位移；
将影响域Ω用NP个所述离散粒子离散，获得所述全场位移场函数的离散近似式：
式中，NP表示所述影响域Ω内的所述离散粒子的总点数，uI对应所述离散粒子I处的位移,ΔVI为离散粒子的体积权重，NI(x)为所述离散粒子I处的形函数，I＝1，2，……NP。
优选地，所述步骤(4)中，所述边界转化方法将所述位移场函数uR(x)分为两部分
式中，b表示位移边界上的离散颗粒，nb表示除去位移边界上的离散颗粒剩下的所有离散颗粒；位移边界上的离散颗粒总数为Nb，Nnb＝NP-Nb。Nb(x)为边界插值点上的形函数，Nnb(x)为非边界插值点的形函数，ub为边界插值点位移，unb为非边界插值点位移；
通过强制位移边界条件消除误差，假定
ub＝(Db)-1g
式中：g表示边界位移，Db通过下式构造
优选地，所述步骤(5)中，
动量方程：
几何方程：εij＝Luij (2)
物理方程：σij＝Dεij (3)
边界条件：
式中：σij表示应力分量，εij表示应变分量，ρ为密度，bi为体积力，L为计算微分算子，D为弹性矩阵，nj表示边界法向量，Ti为边界载荷力，为边界位移；uij为位移，Sσ为载荷边界域，Su为位移边界域，νi为边界位移；
把式(2)、式(3)、式(4)带入式(1)，联立上述方程得到动量方程的等效积分弱形式方程：
上式中，δui,j为i方向位移对xj坐标求偏导，V为求解体积域，s为边界域，A为面积域，uij为ij方向的位移，ui为i方向位移，为i方向的加速度，为Sσ边界域上的外力。
优选地，所述步骤(6)中，将重构的位移场函数uR(x)带入等效积分弱形式方程，并将该方程写成矩阵形式得到下式：
Mü+fint＝fevt
式中，M为质量矩阵，为加速度场，fint为等效内力矩阵，fext为等效外力矩阵。
优选地，所述步骤(7)中，所述中心差分法将时间平均分为n个时间间隔，0，Δt，…t时刻的位移、速度、加速度已知，求解t+Δt时刻对应的解；
t时刻的加速度和速度利用中心差分法表示为
式中，ut-Δt(x)为离散颗粒在t-Δt时刻的位移，ut(x)为离散颗粒在t时刻的位移，ut+Δt(x)为离散颗粒在t+Δt时刻的位移；
t+Δt时刻的位移解答可以根据时间t时刻的运动方程得到，即为下式
式中，M、C和K分别代表影响域的质量矩阵、阻尼矩阵和刚度矩阵；Ft是t时刻外力和；
由此可得到中心差分法的递推公式
若t与t-Δt时刻的位移已知，则可以进一步求解出t+Δt时刻的位移。
利用步骤(5)中的几何方程，将不同时刻位移带入计算，得到形变，即应变，再利用步骤(5)中的物理方程求解，将应变带入计算得到应力。
本发明提供了一种实现结构线性静力仿真的全新的求解方法，与现有技术相比，有益效果在于：
1.简化了复杂模型的前处理过程。现在广泛应用的有限元法，需要在前处理时对结构进行离散化，将几何模型划分成有限个相互连接、互不重合的单元，对于复杂的几何模型，边界几何复杂，需要表征的小特征特别多，这些都使得划分网格的过程耗时耗力；而再生核粒子算法本质上是无网格法，无需划分网格，使用粒子对结构进行离散，只需给定粒子的密度，大大简化了前处理的过程。
2.提升了结构线性静力学仿真的求解精度和求解效率。现广泛应用的有限元法，在每一单元上通过多项式插值构造未知场函数的近似解，因此程序求解的精度和效率会很大程度上依赖于网格的尺寸和质量；而再生核粒子算法，是一种无网格法，回避了网格尺寸和质量对求解精度的影响，另外，再生核粒子算法依据离散颗粒的初始位置，结合选用的修正函数和窗函数，构造形函数，进而根据支持域的半径，生成位移场函数，其构造的形函数和常函数，相比较于有限元法，阶次更高，因此求解精度也更高。
附图说明
构成本申请的一部分附图用来提供对本申请的进一步理解，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。在附图中：
图1是本发明的一种再生核粒子实现结构线性静力仿真方法的流程图。
图2是本发明的一种再生核粒子实现结构线性静力仿真方法的示意性的离散颗粒的影响域二维示图。
图3是本发明中再生核粒子算法与有限元法实现结构线性静力学仿真的求解精度比较图。
具体实施方式
本发明提供一种再生核粒子实现结构线性静力仿真方法，为使本发明的目的、技术方案及效果更加清楚、明确，以下参照附图并举实例对本发明进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
图1是本发明一种再生核粒子实现结构线性静力仿真方法的流程图。
结合图1，本发明所述的一种再生核粒子实现结构线性静力仿真方法，包括如下步骤：
步骤1：按照给定的粒子密度对仿真几何模型进行粒子离散，得到离散颗粒，确定所述离散颗粒的初始位置、影响域Ω和影响域半径ρ。
其中，所述离散颗粒可以为圆形或球形，也可以为其他不同的几何形状，例如，矩形、立方体等。粒子离散后的所述几何模型由多个所述离散颗粒代表。影响域是指一个离散颗粒的影响区域，所述影响域由每个离散颗粒指定，这通过为每个离散颗粒指定影响域半径而完成。
图2是本发明本发明的一种再生核粒子实现结构线性静力仿真方法的示意性的离散颗粒的影响域二维示图。如图2所示，示意性离散颗粒100(黑色实心圆)指定影响域101(虚线圆区域)。设置在所述影响域101内的离散颗粒为实线圆103，而设置在影响域101外的离散颗粒为虚线圆102。每个离散颗粒具有它自身的影响域。所述影响域半径ρ104用来控制所述离散颗粒的影响域的尺寸，对于所有离散颗粒为恒定值或者对于每个单独的所述离散颗粒为不同的。
步骤2：基于再生核粒子算法，选用修正函数和窗函数的类型，构造所述影响域内所述离散颗粒的形函数。
基于再生核粒子算法，所述形函数Φρ(x,x-y)为：
Φρ(x；x-y)＝Cρ(x；x-y)kρ(x-y)
式中，Cρ(x,x-y)为满足再生条件的修正函数或改进函数，它决定近似的完备性与一致性；kρ(x-y)为满足再生条件的窗函数，它是用来表征影响域内离散颗粒的光滑性和紧支性的函数，用于控制计算的连续性。
所述修正函数Cρ(x,x-y)是由多项式基函数的线性组合表示而成，如下：
Cρ(x；x-y)＝C0(x)+C1(x)(x-y)+C2(x)(x-y)2+…+CN(x)(x-y)N＝HT(x-y)C(x)
其中，N是满足再生条件的次数，HT(x-y)是N次多项式基函数的向量；C(x)为系数向量，通过施加N次再生条件来确定；Ci(x)为系数(i＝1,2,3……N)；
HT(x-y)＝[1,x-y,(x-y)2,…,(x-y)N]
C(x)＝[C0(x),C1(x),…,CN(x)]
所述窗函数κρ(x-y)是二维窗函数：
其中，为中间变量。
所述窗函数可以是、但不限于三次样条函数、B样条函数、高斯函数。
由于样条型核函数具有紧支撑特性，因此在大多数数值计算中经常采用该类型的核函数以节省计算时间。本实例中，所述窗函数选用三次样条函数，具体如下：
式中，
步骤3：根据所述影响域内的所述离散颗粒的位移和所述形函数构造位移场函数。
再生核粒子算法通过构造形函数建立位移场函数和离散颗粒处函数的关系，从而得到位移场函数，进一步进行求解。
所述位移场函数uR(x)具体构造如下所示：
式中，u(y)表示所述影响域内的所述离散颗粒的位移；
将影响域Ω用NP个所述离散粒子离散，获得所述全场位移场函数的离散近似式：
式中，NP表示所述影响域Ω内的所述离散粒子的总点数，uI对应所述离散粒子I处的位移,ΔVI为离散粒子的体积权重，NI(x)为所述离散粒子I处的形函数，I＝1，2，……NP。
步骤4：利用边界转化方法进行边界处理，消除边界不一致性引起的误差。
无网格法均不能自然满足位移边界条件，因此在计算过程中必须采用适当的方法来强制满足边界条件，常用的方法有：拉格朗日乘子法、再生核粒子插值技术、修正变分原理法、坐标变换法等。这些方法大多数致力于构建具有Kronecker dalta特性的近似函数，各有其特点，适用于不同场合。例如拉格朗日乘子法其系数矩阵即不正定又不呈带状，所以运算量很大，但这一方法的精确度很高，适用于求解小规模二维问题；修正变原理法的系数矩阵虽然保持正定和带状性，但精度不高。但总的说来，和其他无网络方法一样，还没有一种简单而又完美的方法来处理再生核粒子的边界条件。因此，边界条件的处理是目前无网格方法研究的热点之一，直接关系到无网格方法的实用性。
本实施例中，采用一种对再生核粒子方法适用的一致方法。所述一致方法将所述位移场函数uR(x)分为两部分
式中，b表示位移边界上的离散颗粒，nb表示除去位移边界上的离散颗粒剩下的所有离散颗粒；位移边界上的离散颗粒总数为Nb，Nnb＝NP-Nb。Nb(x)为边界插值点上的形函数，Nnb(x)为非边界插值点的形函数，ub为边界插值点位移，unb为非边界插值点位移；
通过强制位移边界条件消除误差，假定
ub＝(Db)-1g
式中：g表示边界位移，Db通过下式构造
步骤5：联立动量方程和几何方程、物理方程、边界条件构造等效积分弱形式方程。
动量方程：
几何方程：εij＝Luij (2)
物理方程：σij＝Dεij (3)
边界条件：
式中：σij表示应力分量，εij表示应变分量，ρ为密度，bi为体积力，L为计算微分算子，D为弹性矩阵，nj表示边界法向量，Ti为边界载荷力，为边界位移。
把式(2)、式(3)、式(4)带入式(1)，联立上述方程得到动量方程的等效积分弱形式方程：
上式中，δui,j为i方向位移对xj坐标求偏导，V为求解体积域，s为边界域，A为面积域，uij为ij方向的位移，ui为i方向位移，为i方向的加速度，为Sσ边界域上的外力。
步骤6：构造位移矩阵，将所述位移场函数带入所述等效积分弱形式方程求解，得到等效内力和外力矩阵。
将重构的位移场函数uR(x)带入等效积分弱形式方程，并将该方程写成矩阵形式得到下式：
式中：M为质量矩阵，为加速度场，fint为等效内力矩阵，fext为等效外力矩阵。
步骤7：利用中心差分法对等效积分弱形式方程进行代数求解，得到不同时间步的形变和应力。
所述中心差分法将时间平均分为n个时间间隔，0，Δt，…t时刻的位移、速度、加速度已知，求解t+Δt时刻对应的解；
加速度核速度利用中心差分法表示为
式中，ut-Δt(x)为离散颗粒在t-Δt时刻的位移，ut(x)为离散颗粒在t时刻的位移，ut+Δt(x)为离散颗粒在t+Δt时刻的位移；
t+Δt时刻的位移解答可以根据时间t时刻的运动方程得到，即为下式
式中，M、C和K分别代表影响域的质量矩阵、阻尼矩阵和刚度矩阵；Ft是t时刻外力和；
由此可得到中心差分法的递推公式
若t与t-Δt时刻的位移已知，则可以进一步求解出t+Δt时刻的位移。
需要注意的是，该算法开始计算时，即t＝0时，需要知道u-Δt(x)，为此利用和的计算公式可以得到
其中，u0(x)和可以在初始条件中设置，而可利用初始时运动方程得到，即为
利用步骤5中的几何方程，即式(2)，将不同时刻位移带入计算，得到形变，即应变，在利用步骤5中的物理方程，即式(3)求解，将应变带入计算得到应力。
图3是本发明中再生核粒子算法与有限元法实现结构线性静力学仿真的求解精度比较图。可见，再生核粒子方法大大提升了结构线性静力学仿真的求解精度和求解效率。
以上对本发明的具体实施例进行了详细描述，但其只是作为范例，本发明并不限制于以上描述的具体实施例。对于本领域技术人员而言，任何对本发明进行的等同修改和替代也都在本发明的范畴之中。因此，在不脱离本发明的精神和范围下所作的均等变换和修改，都应涵盖在本发明的范围内。
Claims (10)
1.一种再生核粒子算法实现结构线性静力学仿真方法，其特征在于：包括如下步骤：
(1)按照给定的粒子密度对仿真几何模型进行粒子离散，得到离散颗粒，确定所述离散颗粒的初始位置、影响域Ω和影响域半径ρ；
(2)基于再生核粒子算法，选用修正函数和窗函数的类型，构造所述影响域内所述离散颗粒的形函数；
(3)根据所述离散颗粒的位置和所述影响域半径ρ，利用所述影响域Ω内的所述离散颗粒的位移和所述形函数构造位移场函数；
(4)利用边界转化方法进行边界处理，消除边界不一致性引起的误差；
(5)联立动量方程和几何方程、物理方程、边界条件构造等效积分弱形式方程；
(6)构造位移矩阵，将所述位移场函数带入所述等效积分弱形式方程求解，得到等效内力和外力矩阵；
(7)利用中心差分法对等效积分弱形式方程进行代数求解，得到不同时间步的形变和应力。
2.根据权利要求1所述的一种再生核粒子算法实现结构线性静力学仿真方法，其特征在于：所述离散颗粒的几何形状包括、但不限于圆形、球形、矩形、立方体。
3.根据权利要求1所述的一种再生核粒子算法实现结构线性静力学仿真方法，其特征在于：所述步骤(2)中，基于再生核粒子法，所述形函数Φρ(x,x-y)为：
Φρ(x；x-y)＝Cρ(x；x-y)kρ(x-y)
式中，Cρ(x,x-y)为满足再生条件的修正函数，kρ(x-y)为满足再生条件的窗函数；
所述修正函数Cρ(x,x-y)是由多项式基函数的线性组合表示而成，如下：
Cρ(x；x-y)＝C0(x)+C1(x)(x-y)+C2(x)(x-y)2+…+CN(x)(x-y)N
＝HT(x-y)C(x)
其中，N是满足再生条件的次数，HT(x-y)是N次多项式基函数的向量；C(x)为系数向量，通过施加N次再生条件来确定；Ci(x)为系数(i＝1,2,3……N)；
HT(x-y)＝[1,x-y,(x-y)2,…,(x-y)N]
C(x)＝[C0(x),C1(x),…,CN(x)]
所述窗函数κρ(x-y)是二维窗函数：
其中，为中间变量。
4.根据权利要求3所述的一种再生核粒子算法实现结构线性静力学仿真方法，其特征在于：所述窗函数选用三次样条函数，如下：
式中，
5.根据权利要求3所述的一种再生核粒子算法实现结构线性静力学仿真方法，其特征在于：所述窗函数包括、但不限于三次样条函数、B样条函数、高斯函数。
6.根据权利要求1所述的一种再生核粒子算法实现结构线性静力学仿真方法，其特征在于：所述步骤(3)中，所述位移场函数uR(x)为：
式中，u(y)表示所述影响域Ω内的所述离散颗粒的位移；
将影响域Ω用NP个所述离散粒子离散，获得所述全场位移场函数的离散近似式：
式中，NP表示所述影响域Ω内的所述离散粒子的总点数，uI对应所述离散粒子I处的位移,ΔVI为离散粒子的体积权重，NI(x)为所述离散粒子I处的形函数，I＝1，2，……NP。
7.根据权利要求1所述的一种再生核粒子算法实现结构线性静力学仿真方法，其特征在于：所述步骤(4)中，所述边界转化方法将所述位移场函数uR(x)分为两部分
式中，b表示位移边界上的离散颗粒，nb表示除去位移边界上的离散颗粒剩下的所有离散颗粒；位移边界上的离散颗粒总数为Nb，Nnb＝NP-Nb。Nb(x)为边界插值点上的形函数，Nnb(x)为非边界插值点的形函数，ub为边界插值点位移，unb为非边界插值点位移；
通过强制位移边界条件消除误差，假定
ub＝(Db)-1g
式中：g表示边界位移，Db通过下式构造
8.根据权利要求1所述的一种再生核粒子算法实现结构线性静力学仿真方法，其特征在于：所述步骤(5)中，
动量方程：
几何方程：εij＝Luij (2)
物理方程：σij＝Dεij (3)
边界条件：
式中：σij表示应力分量，εij表示应变分量，ρ为密度，bi为体积力，L为计算微分算子，D为弹性矩阵，nj表示边界法向量，Ti为边界载荷力，为边界位移；
把式(2)、式(3)、式(4)带入式(1)，联立上述方程得到动量方程的等效积分弱形式方程：
上式中，δui,j为i方向位移对xj坐标求偏导，V为求解体积域，s为边界域，A为面积域，uij为ij方向的位移，ui为i方向位移，为i方向的加速度，为Sσ边界域上的外力。
9.根据权利要求1所述的一种再生核粒子算法实现结构线性静力学仿真方法，其特征在于：所述步骤(6)中，将重构的位移场函数uR(x)带入等效积分弱形式方程，并将该方程写成矩阵形式得到下式：
式中，M为质量矩阵，为加速度场，fint为等效内力矩阵，fext为等效外力矩阵。
10.根据权利要求1所述的一种再生核粒子算法实现结构线性静力学仿真方法，其特征在于：所述步骤(7)中，所述中心差分法将时间平均分为n个时间间隔，0，Δt，…t时刻的位移、速度、加速度已知，求解t+Δt时刻对应的解；
t时刻的加速度和速度利用中心差分法表示为
式中，ut-Δt(x)为离散颗粒在t-Δt时刻的位移，ut(x)为离散颗粒在t时刻的位移，ut+Δt(x)为离散颗粒在t+Δt时刻的位移；
t+Δt时刻的位移解答可以根据时间t时刻的运动方程得到，即为下式
式中，M、C和K分别代表影响域的质量矩阵、阻尼矩阵和刚度矩阵；Ft是t时刻外力和；
由此可得到中心差分法的递推公式
若t与t-Δt时刻的位移已知，则可以进一步求解出t+Δt时刻的位移。
利用步骤(5)中的几何方程，将不同时刻位移带入计算，得到形变，即应变，再利用步骤(5)中的物理方程求解，将应变带入计算得到应力。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910023402.6A CN109657408A (zh) | 2019-01-10 | 2019-01-10 | 一种再生核粒子算法实现结构线性静力学仿真方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910023402.6A CN109657408A (zh) | 2019-01-10 | 2019-01-10 | 一种再生核粒子算法实现结构线性静力学仿真方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109657408A true CN109657408A (zh) | 2019-04-19 |
Family
ID=66119105
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910023402.6A Pending CN109657408A (zh) | 2019-01-10 | 2019-01-10 | 一种再生核粒子算法实现结构线性静力学仿真方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109657408A (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111220267A (zh) * | 2020-01-09 | 2020-06-02 | 上海索辰信息科技有限公司 | 航空噪声声场复现方法 |
| CN111259516A (zh) * | 2020-01-08 | 2020-06-09 | 上海索辰信息科技有限公司 | 一种基于再生核粒子算法的复合材料铺覆性分析方法 |
| CN111368471A (zh) * | 2020-03-02 | 2020-07-03 | 上海索辰信息科技有限公司 | 设备热传导分析实现方法 |
| CN114077802A (zh) * | 2021-11-22 | 2022-02-22 | 西安交通大学 | 一种利用形函数插值替代核函数近似的粒子建模方法 |
| CN114510677A (zh) * | 2022-01-19 | 2022-05-17 | 西北核技术研究所 | 基于间断有限元的中子输运方程处理方法、计算机程序产品 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060103648A1 (en) * | 2004-11-16 | 2006-05-18 | Livermore Software Technology Corporation | Method and system for mesh-free analysis of general three-dimensional shell structures |
| CN108345741A (zh) * | 2017-12-13 | 2018-07-31 | 湘潭大学 | 基于无网格rkpm的各向异性材料结构二维热变形和热应力分析方法 |
-
2019
- 2019-01-10 CN CN201910023402.6A patent/CN109657408A/zh active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060103648A1 (en) * | 2004-11-16 | 2006-05-18 | Livermore Software Technology Corporation | Method and system for mesh-free analysis of general three-dimensional shell structures |
| CN108345741A (zh) * | 2017-12-13 | 2018-07-31 | 湘潭大学 | 基于无网格rkpm的各向异性材料结构二维热变形和热应力分析方法 |
Non-Patent Citations (2)
| Title |
|---|
| 廉艳平等: "物质点法的理论和应用", 《力学进展》 * |
| 郭彦波: "高应变率下混凝土构件损伤过程的再生核粒子法分析", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111259516A (zh) * | 2020-01-08 | 2020-06-09 | 上海索辰信息科技有限公司 | 一种基于再生核粒子算法的复合材料铺覆性分析方法 |
| CN111220267A (zh) * | 2020-01-09 | 2020-06-02 | 上海索辰信息科技有限公司 | 航空噪声声场复现方法 |
| CN111368471A (zh) * | 2020-03-02 | 2020-07-03 | 上海索辰信息科技有限公司 | 设备热传导分析实现方法 |
| CN114077802A (zh) * | 2021-11-22 | 2022-02-22 | 西安交通大学 | 一种利用形函数插值替代核函数近似的粒子建模方法 |
| CN114077802B (zh) * | 2021-11-22 | 2024-04-16 | 西安交通大学 | 一种利用形函数插值替代核函数近似的粒子建模方法 |
| CN114510677A (zh) * | 2022-01-19 | 2022-05-17 | 西北核技术研究所 | 基于间断有限元的中子输运方程处理方法、计算机程序产品 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109657408A (zh) | 一种再生核粒子算法实现结构线性静力学仿真方法 | |
| Nakata et al. | A fluid–structure interaction model of insect flight with flexible wings | |
| CN109740182A (zh) | 一种基于再生核粒子的无网格物理变形仿真方法 | |
| JP7431029B2 (ja) | 大規模環境用のマルチインスタンス型シミュレーション | |
| CN112100835B (zh) | 一种适用于复杂流动的高效高精度翼型绕流数值模拟方法 | |
| CN112001109A (zh) | 再生核粒子算法实现结构冲击动力学仿真方法 | |
| Ivan et al. | Multi-dimensional finite-volume scheme for hyperbolic conservation laws on three-dimensional solution-adaptive cubed-sphere grids | |
| CN105760602A (zh) | 一种有限体积加权基本无振荡格式的全流场数值模拟方法 | |
| CN108986220A (zh) | 一种加速有限元求解物体网格模型弹性变形的方法 | |
| CN111079326B (zh) | 二维各向异性网格单元度量张量场光滑化方法 | |
| CN104504189A (zh) | 随机激励下大规模结构设计方法 | |
| CN106354954A (zh) | 一种基于叠层基函数的三维力学模态仿真模拟方法 | |
| CN117436322B (zh) | 基于叶素理论的风力机叶片气动弹性仿真方法和介质 | |
| CN109002630B (zh) | 一种超弹性材料的快速仿真方法 | |
| CN111881629B (zh) | 一种气动热-结构热传导耦合非线性降阶模型方法 | |
| McDaniel et al. | Efficient mesh deformation for computational stability and control analyses on unstructured viscous meshes | |
| CN112434451A (zh) | 一种基于分块并行计算的有限元分析方法 | |
| Oliapuram et al. | Realtime forest animation in wind | |
| Poole et al. | Aerofoil design variable extraction for aerodynamic optimization | |
| CN110083946A (zh) | 一种基于无约束优化模型的多状态模型修正的方法 | |
| CN103902764B (zh) | 基于Householder变换的无约束结构静力分析方法 | |
| CN114611423A (zh) | 一种三维多相可压缩流固耦合的快速计算方法 | |
| CN109118561B (zh) | 一种基于位置的层次化动态模拟方法 | |
| CN107463734A (zh) | 动车组车体模态计算方法及装置 | |
| Sun | A robust simplex cut-cell method for adaptive high-order discretizations of aerodynamics and multi-physics problems |
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 | ||
| CB02 | Change of applicant information |
Address after: 201206 2nd floor, building 13, No. 27, Xinjinqiao Road, China (Shanghai) pilot Free Trade Zone, Pudong New Area, Shanghai 201206 Applicant after: Shanghai suochen Information Technology Co.,Ltd. Address before: 201207 floor 2, building 13, No. 27, Xinjinqiao Road, China (Shanghai) pilot Free Trade Zone, Pudong New Area, Shanghai Applicant before: SHANGHAI SUOCHEN INFORMATION TECHNOLOGY CO.,LTD. |
|
| CB02 | Change of applicant information | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190419 |
|
| WD01 | Invention patent application deemed withdrawn after publication |