CN107391788A - 运用三维离散实体解决连续介质构件非线性力学问题的方法 - Google Patents
运用三维离散实体解决连续介质构件非线性力学问题的方法 Download PDFInfo
- Publication number
- CN107391788A CN107391788A CN201710436025.XA CN201710436025A CN107391788A CN 107391788 A CN107391788 A CN 107391788A CN 201710436025 A CN201710436025 A CN 201710436025A CN 107391788 A CN107391788 A CN 107391788A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- contact
- continuous media
- dem
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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
-
- 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/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种运用三维离散实体解决连续介质构件非线性力学问题的方法,针对传统的离散单元法计算连续介质问题的局限性,基于连续介质力学理论,以能量等效为原则,构建了两种三维离散实体离散元模型,通过增加颗粒单元之间的接触弹簧,体现了材料泊松比效应,并推导出了细观模型参数(弹簧刚度)和宏观弹性常数(弹性模量和泊松比)之间的解析关系,对在PFC软件中,能够模拟大泊松比的连续介质是一个很大的突破。另外在颗粒单元间的接触本构方程中增加了塑性部分。本发明能够有效地实现采用离散实体单元法对连续介质弹塑性问题的模拟计算分析,适用于结构或构件的大变形、损伤断裂和倒塌破坏等非线性力学问题。
Description
技术领域
本发明涉及一种三维离散元计算模拟方法,具体涉及采用一种模拟连续介质的三维离散实体单元法,用于解决结构或构件的大变形、损伤断裂和倒塌破坏等非线性力学问题,属于有限元分析技术领域。
背景技术
传统有限元方法在模拟钢结构,例如装配式节点中会存在一些难点:如螺栓和节点的面面接触需要设置接触对、在实际加载的过程中螺栓处容易产生应力集中、模拟时网格划分需加密、杆件大变形时导致不收敛以及断裂等。因此需寻找新的数值分析方法。
由于离散元法不要求满足位移连续和变形协调条件,所以对于结构或构件的大变形、损伤断裂和倒塌破坏等非线性力学问题,离散元法具有明显的算法优越性。而传统的离散元法计算连续介质问题时计算精度远远不如有限元法,并且反映单元作用力关系的参数往往需要通过实验才能确定,严重制约了离散元法的应用。因此,迫切需要一种新的离散元计算方法能有效地模拟连续介质力学问题。
发明内容
技术问题:针对传统离散单元方法在模拟连续介质的一些问题,提供了一种模拟连续介质的三维离散实体单元法,基于连续介质力学理论,以能量等效为原则,构建了两种三维离散实体离散元模型,体现了材料泊松比效应,并推导出了细观模型参数(弹簧刚度)和宏观弹性常数(弹性模量和泊松比)之间的解析关系,对在PFC软件中,能够模拟大泊松比的连续介质是一个很大的突破。另外在颗粒单元间的接触本构方程中增加了塑性部分,从而有效地实现采用离散实体单元法对连续介质弹塑性力学问题的模拟计算分析。
技术方案:本发明提出的模拟连续介质的三维离散实体单元法,包括以下具体步骤:
(1)对结构或构件进行离散实体建模
采用体心立方模型(简称BCC模型——Body Centered Cubic)或面心立方模型(简称FCC模型——Face Centered Cubic)对结构或构件进行建模,建立模型的颗粒单元为规则排列且单元尺寸相同,材料的变形完全由单元的接触弹簧来存储和表示,输入模型参数以及外力信息,并确定边界条件;
(2)t=0时刻,对各颗粒单元的内力、速度等赋初值;
(3)对所有单元:计算t(t≥0)时刻的外力及阻尼力,应用牛顿第二定律求解运动控制方程,得到t+Δt时刻的颗粒单元的位置与速度;
任取一个单元α,设有n个单元与其相邻,作用在单元α上的外力为Fext,外力矩为Mext。根据牛顿第二定律,其运动控制方程为
其中,m,J分别为单元α的质量和转动惯量;r,ω分别为单元α的位置矢量和角速度矢量;分别为相邻单元j对单元α产生的接触力和接触力矩;t为时间;
(4)计算位移增量ΔU,并计算试探接触力增量从而形成t+Δt的接触力S,根据屈服准则判断当前的颗粒单元的弹塑性状态;
对于颗粒-颗粒接触,接触面的单位法向量可定义为:
其中,代表颗粒的形心位置矢量;ni代表颗粒单元之间的接触法向单位矢量;d代表颗粒单元A、B球心之间的距离,其公式为:
球体颗粒-球体颗粒接触时,法向的相对接触位移Un为:
Un=R[A]+R[B]-d
(4)
其中,R[A]、R[B]分别为颗粒单元A、B的半径;
(5)根据弹塑性接触本构方程,计算Δt时步内颗粒单元的接触力增量ΔF,并更新当前时刻的颗粒单元间的接触力;
法向接触力矢量Fi n为:
Fi n=KnUnni (5)
而切向接触力矢量Fi s以增量的形式进行计算。球体颗粒之间的初始切向接触力为零,随后由于颗粒之间出现了切向位移增量,从而产生附加接触力,这些附加的接触力经过叠加运算,不断的更新球体之间的接触力的大小,即:
Fi s←Fi s+ΔFi s (6)
切向接触力增量矢量为:
ΔFi s=-KsΔUis (7)
其中,Kn、Ks分别为接触的法向和切向刚度,ΔUi s为切向的相对位移增量;
(6)根据平衡方程,计算各颗粒单元的内力、阻尼力以及外力,为下一个t+Δt时刻的计算作准备,如此循环往复,直至模型达到稳定状态;
具体地,所述的BCC模型由八个处于立方体角上的颗粒和一个处于立方体中心的颗粒构成,所述的颗粒半径r均相等,所述的角上八个颗粒与中心颗粒相切,并通过弹簧***(包括一个法向弹簧和两个切向弹簧)连接。所述的FCC模型由分布在立方体的八个角上和六个面中心的颗粒构成,所述的颗粒半径r均相等,所述的面中心的颗粒与该面四个角上的颗粒相切,中心颗粒与周围的12个颗粒通过弹簧***(包括一个法向弹簧和两个切向弹簧)连接;
所述的BCC模型的弹簧***,其特征在于:基于能量等效原则,推导出了细观模型参数(弹簧刚度)和宏观弹性常数(弹性模量和泊松比)之间的解析关系:
法向弹簧刚度:
切向弹簧刚度:
具体地,所述的FCC模型的弹簧***,基于能量等效原则,推导出了细观模型参数(弹簧刚度)和宏观弹性常数(弹性模量和泊松比)之间的解析关系:
法向弹簧刚度:
切向弹簧刚度:
具体地,屈服法则是根据第四强度理论,只要材料的畸变能密度达到材料限值,材料便发生屈服,因此离散实体单元法的屈服方程可表示为:
其中,是单位体积的剪切变形能,是材料的屈服应力,G为剪切模量;
具体地,当弹簧接触内力满足屈服条件时,应建立弹塑性接触本构方程。接触位移包括弹性和塑性两部分,需建立离散实体单元法流动准则,进而推导弹塑性接触本构方程:
其中,ΔFn、ΔFs、ΔFt分别为法向和切向接触力,ΔUn、ΔUs、ΔUt分别为法向和切向位移增量,dλ为比例系数,Ke为弹性接触刚度矩阵。
有益效果:在离散实体单元法中,求解时需遍历结构中所有单元,并计算每个单元的接触内力和运动方程式,但无需组集所谓的刚度矩阵和迭代求解,计算流程比较简单清晰。结构越复杂,单元数量就越多,相应的计算量也越大,但只是循环重复计算而已,并不会给分析带来实质性困难。计算效率与传统方法相比,在大型结构和复杂力学行为的模拟中离散实体单元法将表现出更大的优势。
结构的弹塑性行为仅与内力的求解有关,除需采用屈服准则判断单元受力状态、以及采用弹塑性接触本构模型计算单元接触力外,其它的求解流程都不改变。
构建的三维离散实体元模型能够体现泊松比效应,对在PFC3D软件中,能够模拟大泊松比的连续介质是一个突破。
下面结合说明书附图和具体算例对本发明进行进一步的说明,验证了本发明的使用方式和有效性,其具体的数据及内容不在本发明保护范围之内。
附图说明
图1为本发明的流程图;
图2为BCC模型;
图3为FCC模型;
图4为受弯矩M作用的悬臂梁示意图;
图5为受不同端弯矩作用的悬臂梁变形图(离散元结果);
图6为悬臂端的荷载-水平位移曲线;
图7为悬臂端的荷载-竖向位移曲线。
具体实施方式
在线弹性力学中,由于假设结构位移很小,荷载与位移为线性关系,基于变形前位置给出的平衡方程求解值即可满足工程上的精度要求。但当结构大变形时,大变形引起二阶或更高阶效应,荷载与位移为非线性关系,结构刚度变化明显,平衡方程必须建立在变形后的状态下。离散单元法中各单元之间没有变形协调条件的约束,采用动态松弛法求解,单元位置随计算时步的增加而更新,并用于计算下一时步的内力,由此,离散单元法能有效考虑结构大变形引起的高阶效应。
下面用算例验证三维离散元实体模型在模拟大变形构件中的优越性。
算例:如图4所示,一尺寸为1m×1m×10m的悬臂梁左端部固定,悬臂端受弯矩M作用。材料参数为:弹性模量E=2.06×1011N/m2;泊松比μ=0.3;密度ρ=7850kg/m3。本算例考察了悬臂端在弯矩M作用下水平和竖向的位移变化情况。
当悬臂梁受横向端部弯矩M作用时,无因次量纲ML/(EI)与梁变形的关系如表1所示,PFC计算结果如图5所示,可以看出:结构承受静力荷载过程中,随着弯矩的逐渐增加,结构出现大变形和大转角,并且与表1所列的结果相吻合,验证了BCC模型的有效性。
表1悬臂梁端弯矩与变形的关系
图6和图7分别为悬臂端的荷载-水平位移曲线和悬臂端的荷载-竖向位移曲线,可见PFC结果与ANSYS结果很好地吻合,三维离散实体单元法在分析结构大变形问题时,不需要修改计算模型,直接可得到准确的几何非线性解,从而验证了本发明的有效性。
以上所述仅为本发明的优选实施方式而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化,在本发明的原理和技术思想的范围内,对这些实施方式进行多种变化、修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种运用三维离散实体解决连续介质构件非线性力学问题的方法,其特征在于包括以下步骤:
步骤1,对结构或构件进行离散实体建模:采用体心立方模型或面心立方模型对结构或构件进行建模,建立模型的颗粒单元为规则排列且单元尺寸相同,材料的变形完全由单元的接触弹簧来存储和表示;
步骤2,t=0时刻,对各颗粒单元的内力、速度等赋初值;
步骤3,对所有颗粒单元:计算t(t≥0)时刻的外力及阻尼力,应用牛顿第二定律求解运动控制方程,得到t+Δt时刻的颗粒单元的位置与速度;
步骤4,计算位移增量ΔU,并计算试探接触力增量从而形成t+Δt的接触力S,根据屈服准则判断当前的颗粒单元的弹塑性状态;
步骤5,根据弹塑性接触本构方程,计算Δt时步内颗粒单元的接触力增量ΔS,并更新当前时刻的颗粒单元间的接触力;
步骤6,根据平衡方程,计算各颗粒单元的内力、阻尼力以及外力,为下一个t+Δt时刻的计算作准备,如此循环往复,直至模型达到稳定状态。
2.根据权利要求书1所述的运用三维离散实体解决连续介质构件非线性力学问题的方法,其特征在于:所述体心立方模型由八个处于立方体角上的颗粒和一个处于立方体中心的颗粒构成,所述的颗粒半径r均相等,所述的角上八个颗粒与中心颗粒相切,并通过弹簧***连接,所述弹簧***包括一个法向弹簧和两个切向弹簧;
所述面心立方模型由分布在立方体的八个角上和六个面中心的颗粒构成,所述的颗粒半径r均相等,所述的面中心的颗粒与该面四个角上的颗粒相切,中心颗粒与周围的12个颗粒通过弹簧***连接。
3.根据权利要求书1所述的运用三维离散实体解决连续介质构件非线性力学问题的方法,其特征在于:所述体心立方模型细观模型参数与宏观弹性常数之间的解析关系为:
法向弹簧刚度:
切向弹簧刚度:
其中,kn为法向弹簧刚度,ks、kt分别为两个切向弹簧刚度,V为一个颗粒单元所占据的平均体积,r为颗粒半径,E为弹性模量,μ为材料的泊松比;所述细观模型参数为弹簧刚度,所述宏观弹性常数为弹性模量和泊松比。
4.根据权利要求书1所述的运用三维离散实体解决连续介质构件非线性力学问题的方法,其特征在于:所述面心立方模型细观模型参数与和宏观弹性常数之间的解析关系为:
法向弹簧刚度:
切向弹簧刚度:
其中,kn为法向弹簧刚度,ks、kt分别为两个切向弹簧刚度,V为一个颗粒单元所占据的平均体积,r为颗粒半径,E为弹性模量,μ为材料的泊松比;所述细观模型参数为弹簧刚度,所述宏观弹性常数为弹性模量和泊松比。
5.根据权利要求书1所述的运用三维离散实体解决连续介质构件非线性力学问题的方法,其特征在于:根据第四强度理论,离散实体单元法的屈服方程表示为:
<mrow>
<mi>f</mi>
<mo>=</mo>
<msubsup>
<mi>U</mi>
<mi>&phi;</mi>
<mi>C</mi>
</msubsup>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>&sigma;</mi>
<mi>s</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mn>6</mn>
<mi>G</mi>
</mrow>
</mfrac>
</mrow>
其中,是单位体积的剪切变形能,是材料的屈服应力,G为剪切模量。
6.根据权利要求书1所述的运用三维离散实体解决连续介质构件非线性力学问题的方法,当弹簧接触内力满足屈服条件时,建立弹塑性接触本构方程;
接触位移包括弹性和塑性两部分,建立离散实体单元法流动准则,推导弹塑性接触本构方程:
<mrow>
<msub>
<mi>&Delta;F</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<msub>
<mi>K</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&Delta;U</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<mi>d</mi>
<mi>&lambda;</mi>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>f</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>F</mi>
<mi>n</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
<msubsup>
<mi>&Delta;F</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
<mn>1</mn>
</msubsup>
<mo>=</mo>
<msub>
<mi>K</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&Delta;U</mi>
<mi>s</mi>
</msub>
<mo>-</mo>
<mi>d</mi>
<mi>&lambda;</mi>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>f</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>F</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
<msub>
<mi>&Delta;F</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<msub>
<mi>K</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&Delta;U</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<mi>d</mi>
<mi>&lambda;</mi>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>f</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>F</mi>
<mi>t</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
其中,ΔFn、ΔFs、ΔFt分别为法向和切向接触力,ΔUn、ΔUs、ΔUt分别为法向和切向位移增量,dλ为比例系数,Ke为弹性接触刚度矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710436025.XA CN107391788B (zh) | 2017-06-09 | 2017-06-09 | 运用三维离散实体解决连续介质构件非线性力学问题的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710436025.XA CN107391788B (zh) | 2017-06-09 | 2017-06-09 | 运用三维离散实体解决连续介质构件非线性力学问题的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107391788A true CN107391788A (zh) | 2017-11-24 |
CN107391788B CN107391788B (zh) | 2020-10-02 |
Family
ID=60332400
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710436025.XA Active CN107391788B (zh) | 2017-06-09 | 2017-06-09 | 运用三维离散实体解决连续介质构件非线性力学问题的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107391788B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109492285A (zh) * | 2018-10-30 | 2019-03-19 | 河海大学 | 一种可变形三维任意圆化凸多面体块体离散单元法 |
CN109800462A (zh) * | 2018-12-20 | 2019-05-24 | 江苏大学 | 一种空心杆状柔性体动力学特性的建模与仿真方法 |
CN109815599A (zh) * | 2019-01-28 | 2019-05-28 | 南京大学 | 一种离散元材料自动训练方法 |
CN113108002A (zh) * | 2021-04-23 | 2021-07-13 | 天津大学 | 一种具有颗粒阻尼特性的立体形式声子晶体隔振降噪装置 |
CN113297749A (zh) * | 2021-06-07 | 2021-08-24 | 吉林大学 | 一种基于弹塑性的连接力学模型 |
CN113887094A (zh) * | 2021-09-06 | 2022-01-04 | 东南大学 | 运用多尺度离散实体元解决连续体结构力学问题的仿真方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU932358A1 (ru) * | 1980-11-05 | 1982-05-30 | Калининский Ордена Трудового Красного Знамени Политехнический Институт | Устройство дл упруго-пластического деформировани длинномерных образцов при сжатии |
CN101667220A (zh) * | 2009-09-22 | 2010-03-10 | 同济大学 | 盾构隧道接头构造的模拟方法 |
CN102385663A (zh) * | 2011-08-22 | 2012-03-21 | 天津大学 | 一种高层钢框架结构的地震倒塌分析方法 |
CN105659887B (zh) * | 2011-06-15 | 2014-01-22 | 北京空间飞行器总体设计部 | 一种月壤力学性能模拟*** |
CN104091009A (zh) * | 2014-07-01 | 2014-10-08 | 东南大学 | 颗粒流与有限差分法耦合计算方法 |
CN106646597A (zh) * | 2016-12-14 | 2017-05-10 | 中国石油大学(北京) | 基于弹簧网络模型的正演模拟方法及装置 |
-
2017
- 2017-06-09 CN CN201710436025.XA patent/CN107391788B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU932358A1 (ru) * | 1980-11-05 | 1982-05-30 | Калининский Ордена Трудового Красного Знамени Политехнический Институт | Устройство дл упруго-пластического деформировани длинномерных образцов при сжатии |
CN101667220A (zh) * | 2009-09-22 | 2010-03-10 | 同济大学 | 盾构隧道接头构造的模拟方法 |
CN105659887B (zh) * | 2011-06-15 | 2014-01-22 | 北京空间飞行器总体设计部 | 一种月壤力学性能模拟*** |
CN102385663A (zh) * | 2011-08-22 | 2012-03-21 | 天津大学 | 一种高层钢框架结构的地震倒塌分析方法 |
CN104091009A (zh) * | 2014-07-01 | 2014-10-08 | 东南大学 | 颗粒流与有限差分法耦合计算方法 |
CN106646597A (zh) * | 2016-12-14 | 2017-05-10 | 中国石油大学(北京) | 基于弹簧网络模型的正演模拟方法及装置 |
Non-Patent Citations (6)
Title |
---|
B.D. LE, F. DAU ET AL: "Modeling damages and cracks growth in composite with a 3D discrete element method", 《COMPOSITES PART B: ENGINEERING》 * |
CHIAO-YIN LU ET AL: "Forecasting landslide hazard by the 3D discrete element method: A case study of the unstable slope in the Lushan hot spring district, central Taiwan", 《ENGINEERING GEOLOGY》 * |
冯春等: "基于连续介质力学的块体单元离散弹簧法研究", 《岩石力学与工程学报》 * |
刘凯欣等: "离散元法在求解三维冲击动力学问题中的应用", 《固体力学学报》 * |
张洪武等: "基于非线性接触本构的颗粒材料离散元数值模拟", 《岩土工程学报》 * |
魏怀鹏等: "基于连续介质模型的离散元方法中弹簧性质研究", 《岩石力学与工程学报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109492285A (zh) * | 2018-10-30 | 2019-03-19 | 河海大学 | 一种可变形三维任意圆化凸多面体块体离散单元法 |
CN109492285B (zh) * | 2018-10-30 | 2022-08-26 | 河海大学 | 一种可变形三维任意圆化凸多面体块体离散单元法 |
CN109800462A (zh) * | 2018-12-20 | 2019-05-24 | 江苏大学 | 一种空心杆状柔性体动力学特性的建模与仿真方法 |
CN109815599A (zh) * | 2019-01-28 | 2019-05-28 | 南京大学 | 一种离散元材料自动训练方法 |
CN109815599B (zh) * | 2019-01-28 | 2023-03-24 | 南京大学 | 一种离散元材料自动训练方法 |
CN113108002A (zh) * | 2021-04-23 | 2021-07-13 | 天津大学 | 一种具有颗粒阻尼特性的立体形式声子晶体隔振降噪装置 |
CN113108002B (zh) * | 2021-04-23 | 2022-12-13 | 天津大学 | 一种具有颗粒阻尼特性的立体形式声子晶体隔振降噪装置 |
CN113297749A (zh) * | 2021-06-07 | 2021-08-24 | 吉林大学 | 一种基于弹塑性的连接力学模型 |
CN113297749B (zh) * | 2021-06-07 | 2023-02-03 | 吉林大学 | 一种基于弹塑性连接力学模型的玉米脱粒过程仿真方法 |
CN113887094A (zh) * | 2021-09-06 | 2022-01-04 | 东南大学 | 运用多尺度离散实体元解决连续体结构力学问题的仿真方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107391788B (zh) | 2020-10-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107391788A (zh) | 运用三维离散实体解决连续介质构件非线性力学问题的方法 | |
Caliri Jr et al. | A review on plate and shell theories for laminated and sandwich structures highlighting the Finite Element Method | |
Carrera et al. | Beam structures: classical and advanced theories | |
Tran et al. | Dynamic analysis of sandwich auxetic honeycomb plates subjected to moving oscillator load on elastic foundation | |
Han et al. | Efficient optimization of transient dynamic problems in MEMS devices using model order reduction | |
Yakoub et al. | Three dimensional absolute nodal coordinate formulation for beam elements: implementation and applications | |
Yang et al. | Unstructured dynamic meshes with higher-order time integration schemes for the unsteady Navier-Stokes equations | |
Cárdenas et al. | Numerical validation of a finite element thin‐walled beam model of a composite wind turbine blade | |
Phung‐Van et al. | A cell‐based smoothed discrete shear gap method (CS‐FEM‐DSG3) based on the C0‐type higher‐order shear deformation theory for dynamic responses of Mindlin plates on viscoelastic foundations subjected to a moving sprung vehicle | |
CN101887474B (zh) | 基于有限元法与广义傅里叶级数法的结构振动分析方法 | |
Zhuang et al. | Modal and aeroelastic analysis of trapezoidal corrugated-core sandwich panels in supersonic flow | |
CN109740182A (zh) | 一种基于再生核粒子的无网格物理变形仿真方法 | |
Thai et al. | Finite-element modeling for static bending analysis of rotating two-layer FGM beams with shear connectors resting on imperfect elastic foundations | |
Yu et al. | Integration of geometry and analysis for vehicle system applications: continuum-based leaf spring and tire assembly | |
CN104504189B (zh) | 随机激励下大规模结构设计方法 | |
CN112001109A (zh) | 再生核粒子算法实现结构冲击动力学仿真方法 | |
Orzechowski et al. | Integration of the equations of motion of multibody systems using absolute nodal coordinate formulation | |
Tai et al. | A hierarchic high-order Timoshenko beam finite element | |
Zhang et al. | Elastodynamic analysis of regular polygonal CNT-reinforced composite plates via FSDT element-free method | |
CN106354954A (zh) | 一种基于叠层基函数的三维力学模态仿真模拟方法 | |
Wang et al. | The spectral difference method for 2d Euler equations on unstructured grids | |
Tran et al. | Methods of fluid–structure coupling in frequency and time domains using linearized aerodynamics for turbomachinery | |
Nguyen et al. | A refined quasi-3D isogeometric nonlinear model of functionally graded triply periodic minimal surface plates | |
Ngo-Cong et al. | A numerical procedure based on 1D-IRBFN and local MLS-1D-IRBFN methods for fluid-structure interaction analysis | |
VeisiAra et al. | Computational analysis of the nonlinear vibrational behavior of perforated plates with initial imperfection using NURBS-based isogeometric approach |
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 |