CN115081146A - 一种颗粒阻尼动力吸振器耦合数值计算方法 - Google Patents
一种颗粒阻尼动力吸振器耦合数值计算方法 Download PDFInfo
- Publication number
- CN115081146A CN115081146A CN202210845856.3A CN202210845856A CN115081146A CN 115081146 A CN115081146 A CN 115081146A CN 202210845856 A CN202210845856 A CN 202210845856A CN 115081146 A CN115081146 A CN 115081146A
- Authority
- CN
- China
- Prior art keywords
- particle
- dem
- contact
- particles
- mbd
- 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
- 239000002245 particle Substances 0.000 title claims abstract description 121
- 238000013016 damping Methods 0.000 title claims abstract description 50
- 239000006096 absorbing agent Substances 0.000 title claims abstract description 39
- 238000004364 calculation method Methods 0.000 title claims abstract description 30
- 230000008878 coupling Effects 0.000 title claims abstract description 21
- 238000010168 coupling process Methods 0.000 title claims abstract description 21
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 21
- 238000000034 method Methods 0.000 claims abstract description 22
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 8
- 230000010354 integration Effects 0.000 claims abstract description 5
- 230000001133 acceleration Effects 0.000 claims description 20
- 238000006073 displacement reaction Methods 0.000 claims description 7
- 239000000463 material Substances 0.000 claims description 5
- 230000005484 gravity Effects 0.000 claims description 2
- 229910052739 hydrogen Inorganic materials 0.000 claims description 2
- 239000001257 hydrogen Substances 0.000 claims description 2
- 125000004435 hydrogen atom Chemical class [H]* 0.000 claims description 2
- -1 i.e. Substances 0.000 claims description 2
- 238000013519 translation Methods 0.000 claims description 2
- 238000012821 model calculation Methods 0.000 claims 1
- 230000009467 reduction Effects 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012356 Product development Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 150000004292 cyclic ethers Chemical class 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000001737 promoting effect 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/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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明属于减振降噪领域,具体地说,是一种颗粒阻尼动力吸振器耦合数值计算方法,使用MBD和DEM算法分别对每个机械部分和每个颗粒的运动方程进行数值积分析,用直接积分法将二阶运动微分方程转化为一阶微分方程,然后用龙格‑库塔方法求解初值,再用耦合多体动力学和离散元模型计算,本发明采用了多体动力学和离散元耦合的方法,从多体动力学和离散元的原理出发,所建立的颗粒阻尼动力吸振器耦合方法,结果准确,能够充分的反映出结构的运动状态,为颗粒阻尼动力吸振器计算研究提供有效的技术支持。
Description
技术领域
本发明属于减振降噪领域,具体地说,是一种颗粒阻尼动力吸振器耦合数值计算方法。
背景技术
目前,基于调谐质量的强烈振动可以改善颗粒阻尼的阻尼效果,而阻尼效果的改善也将促进动力吸振器的振动控制性能,因此动力吸振器和颗粒阻尼的结合具有一定的研究意义。
耦合的颗粒阻尼动力吸振器的总阻尼效果是颗粒阻尼和动力吸振器两系统相互促进的,而不是动力吸振器和颗粒阻尼的线性叠加。因此,有必要讨论颗粒阻尼动力吸振器的阻尼特性。目前,现有技术中还没有对颗粒阻尼动力吸振耦合的复合阻尼器的仿真学进行研究过,也没有提出相关的算法。
发明内容
为了解决颗粒阻尼动力吸振器的耦合仿真问题,本发明提出一种颗粒阻尼动力吸振器耦合数值计算方法,该计算方法是基于多体动力学和离散元的仿真方法,可以高效准确的对颗粒阻尼动力吸振器进行计算仿真研究,对颗粒阻尼动力吸振器的数值计算具有重要的意义,具有良好的市场应用前景,弥补市场空白。
本发明采用的具体技术方案如下:
一种颗粒阻尼动力吸振器耦合数值计算方法,包括以下步骤:
步骤一:给出机械零件和阻尼颗粒的几何参数和材料特性,并设置了初始条件和边界条件;MBD计算tMBD和DEM计算tDEM的物理时间在计算开始时设置为零;
步骤二:由于MBD计算的时间步长ΔtMBD远大于DEM计算的时间步长ΔtDEM,因此DEM计算首先开始,并遵循DEM循环,DEM循环流程如下:
流程a、通过接触算法计算颗粒之间以及颗粒与机械零件之间的接触位移;
流程b、利用方程式确定接触力:
流程c、计算每个颗粒的线加速度和角加速度;
流程d、更新每个颗粒的位置和速度;
流程e、如果DEM计算tDEM的物理时间小于下一个MBD步骤的物理时间,计算仍然遵循DEM循环,从流程a重复到流程d,否则,调用MBD方案,计算转到步骤三。
步骤三:DEM将机械零件上的等效颗粒力传递给MBD。在DEM循环之间的每个MBD时间步长执行MBD方案。MBD循环流程如下:
流程a、计算广义力。
流程b、通过化简方程求解MBD的控制方程。
流程c、利用RK4方法更新机械零件的位置和速度。
流程d、一旦MBD计算在每个时间步ΔtMBD结束,回到步骤二,计算相应的平移和旋转的机械零件发送回DME循环,然后确定阻尼颗粒的运动行为和等效颗粒部分再次机械零件。
本发明的有益效果:本发明相对于实验来说能明显的缩短产品开发周期、减少实验次数、降低产品的开发成本;动力吸振器和颗粒阻尼器作为低频振动控制的重要措施,动力吸振器能在不移动机械设备、不改变基座结构、不修改船体结构的情况下,添加颗粒阻尼系统吸收主系统振动能量,抑制结构振动,因此研究颗粒阻尼和动力吸振器耦合是抑制低频噪声的一个突破点。
附图说明
图1是本发明颗粒阻尼动力吸振器示意图。
图2a是颗粒与颗粒碰撞原理图。
图2b是颗粒与壁面碰撞原理图。
图3a是颗粒与颗粒碰撞原理图。
图3b是颗粒与壁面碰撞原理图。
图4是颗粒阻尼动力吸振器耦合流程图。
图5是多体动力学和离散元耦合仿真和实验对比图。
具体实施方式
为了加深对本发明的理解,下面将结合附图和实施例对本发明做进一步详细描述,该实施例仅用于解释本发明,并不对本发明的保护范围构成限定。
本发明实施例是对颗粒阻尼动力吸振器进行多体动力学和离散元耦合数值计算,对其方法进行分析,通过本实施例的分析,对本发明做详细的说明,具体实施方式分为以下步骤:
MBD和DEM使用算法,分别对每个机械部分和每个颗粒的运动方程进行数值积分析。用直接积分法将二阶运动微分方程转化为一阶微分方程。然后用龙格-库塔方法求解初值,图4为耦合多体动力学和离散元(MBD-DEM)模型的计算流程图,具体流程如下:
步骤一:给出了机械部件和阻尼颗粒的几何参数和材料特性,并设置了初始条件和边界条件;如表1,MBD计算tMBD、DEM计算tDEM的物理时间,在计算开始时设置为零。
步骤二:由于MBD计算的时间步长(ΔtMBD)远大于DEM计算的时间步长(ΔtDEM),因此DEM计算首先开始,并遵循DEM循环。DEM循环如下:
a)通过接触算法计算颗粒之间以及颗粒与机械零件之间的接触位移(重叠);
b)利用方程式(1)-(13)确定接触力,具体步骤如下:
1、颗粒在离散元(DEM)中,每个颗粒的运动通常由颗粒i表示,遵守牛顿第二运动定律,方程(1)描述质心G的平移运动,方程(2)描述整个质点的旋转运动。
其中m(i)是颗粒i的质量,I(i)是相应的关于质心G的质量惯性矩,是颗粒i的位置矢量,是平移加速度,是角速度,是角加速度。是我在接触点的接触力作用在颗粒j,是重力加速度,上的接触点颗粒的位置向量i在接触j。是颗粒上的联系人数量我和象征“×”表示向量叉乘。
2、在颗粒动力学中,颗粒之间以及颗粒与机械部件之间的接触是一种直接且简单的方法。图2和图3表示出颗粒之间以及颗粒与机械部件之间的接触示意图。颗粒间接触的判据为否则为无接触。同样,颗粒-壁面接触的接触判据由给出,否则为无接触。在上述不等式中,为颗粒j的位置向量,为接触点在壁上的位置向量,ri,rj分别表示颗粒i和j的半径。
3、考虑两个接触的颗粒,即颗粒i和j,mi和mj分别表示颗粒i、j的质量;ri和rj表示这些颗粒各自的半径。分别定义m*和r*为等效质量和半径,表示为m*=mimj/(mi+mj)和r*=rirj/(ri+rj)。两个颗粒具有相同的材料性质,剪切模量为G,泊松比为v。当颗粒i与颗粒j以非零的相对速度斜向碰撞时,颗粒i与颗粒j之间存在接触力。作用在颗粒上的接触力包括法向分量和切向分量表示为:
4、各分量用弹簧力和阻尼力之和表示,表示为:
Fcn=Fcn,s+Fcn,d (4)
Fct=Fct,s+Fct,d (5)
5、在法向上,式(4)中的Fcn,s,Fcn,d分别为法向弹簧力和法向阻尼力。
Fct,s=Fct,s,(n-1)+ΔFct,s (8)
6、ΔFct,s和Fct,s,(n-1)分别为当前和之前时间步长的切向弹簧力,ΔFct,s为切向弹簧力的增量。ξt为切向阻尼比,kt为切向接触刚度,为相对切向速度(颗粒j相对颗粒i)。增量ΔFct,s进一步由ΔFct,s=-ksΔδt给出,Δδt为增量切向位移,表示为其中Δt为时间步长。
切向力受库仑摩擦定律限制,其表达式为:
|Fct|≤μFcn (11)
其中μ是摩擦系数,该接触力模型所需的参数包括剪切模量(G)、泊松比(υ)、颗粒间
的摩擦系数(μb)、颗粒与机械部件之间的摩擦系数(μw)、颗粒间的恢复系数(eb)、颗粒与机械部件之间的恢复系数(ew)。
8、根据牛顿力学方程,含颗粒阻尼动力吸振器方程可以表示为下式:
在方程式(12)中,FCV(t)分别表示颗粒对盒的垂直合力,如式(13)所示。
式中,Fcv,j为颗粒盒第j个接触点接触力的垂直分量,Ncontact为颗粒盒上的总接触数,法向力Fcn一般可以表示为其中K为接触刚度;χ为滞回阻尼因子;a和b是幂指数,右边的第一项和第二项分别表示弹性力分量和阻尼力分量。
c)利用式(1)、(2)计算每个颗粒的线加速度和角加速度;
d)通过显式积分更新每个颗粒的位置和速度。平动速度和角速度用中间间隔的速度值表示为方程式(14)和(15),该位置如式(16)所示。
e)如果DEM(tDEM)的物理时间小于下一个MBD步骤的物理时间(tMBD+ΔtMBD),计算仍然遵循DEM循环,从步骤2(a)重复到步骤2(d)。否则,调用MBD方案,计算转到步骤三。
步骤三:DEM将机械零件上的等效颗粒力传递给MBD。在DEM循环之间的每个MBD时间步长执行MBD方案,MBD循环如下:
a)计算广义力;
b)通过化简方程求解MBD的控制方程;
c)利用(RK4)方法更新机械零件的位置和速度;
d)一旦MBD计算在每个时间步ΔtMBD结束,回到步骤(2),计算相应的平移和旋转的机械零件发送回DME循环,然后确定阻尼颗粒的运动行为和等效颗粒部分再次机械零件;
数值计算具体实施例数据如下表1所示,振动台对基座施加0.12g的加速度,用上述方法,对颗粒阻尼动力吸振器系统进行数值计算,分别对无动力吸振器、有动力吸振器、颗粒阻尼动力吸振器进行仿真和实验的对比,实验与仿真对比如图5所示,可以看出有动力吸振器和无动力吸振器相比加速度的峰值下降了37.5%,加颗粒之后,加速度峰值下降明显,含颗粒动力吸振器和有动力吸振器相比两个峰的加速度峰值分别下降了78%和89%,实验和仿真也吻合较好。
表1动力吸振和颗粒阻尼耦合研究模型参数
以上显示和描述了本发明的基本原理、主要特征及优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (8)
1.一种颗粒阻尼动力吸振器耦合数值计算方法,其特征在于,使用MBD和DEM算法分别对每个机械部分和每个颗粒的运动方程进行数值积分析,用直接积分法将二阶运动微分方程转化为一阶微分方程,然后用龙格-库塔方法求解初值,再用耦合多体动力学和离散元模型计算。
2.根据权利要求1所述的颗粒阻尼动力吸振器耦合数值计算方法,其特征在于,具体包括以下步骤:
步骤一、给出了机械部件和阻尼颗粒的几何参数和材料特性,并设置了初始条件和边界条件,MBD计算tMBD、DEM计算tDEM的物理时间,在计算开始时设置为零;
步骤二、由于MBD计算的时间步长ΔtMBD大于DEM计算的时间步长ΔtDEM,因此DEM计算首先开始,并遵循DEM循环;
步骤三、DEM将机械零件上的等效颗粒力传递给MBD,在DEM循环之间的每个MBD时间步长执行MBD循环方案。
3.根据权利要求2所述的颗粒阻尼动力吸振器耦合数值计算方法,其特征在于,所述步骤二中,DEM循环流程如下:
流程a、通过接触算法计算颗粒之间以及颗粒与机械零件之间的接触位移;
流程b、利用方程式确定接触力:
流程c、计算每个颗粒的线加速度和角加速度;
流程d、更新每个颗粒的位置和速度;
流程e、如果DEM计算tDEM的物理时间小于下一个MBD步骤的物理时间,计算仍然遵循DEM循环,从流程a重复到流程d,否则,调用MBD方案,计算转到步骤三。
4.根据权利要求3所述的颗粒阻尼动力吸振器耦合数值计算方法,其特征在于,所述步骤二中DEM循环的流程b利用方程式确定接触力,具体方法如下:
颗粒在离散元DEM中,每个颗粒的运动由颗粒i表示,遵守牛顿第二运动定律:方程(1)描述质心G的平移运动,方程(2)描述整个质点的旋转运动:
其中m(i)是颗粒i的质量,I(i)是相应的关于质心G的质量惯性矩,是颗粒i的位置矢量,是平移加速度,是角速度,是角加速度,是i在接触点的接触力作用在颗粒j,是重力加速度,上的接触点颗粒的位置向量i在接触j;
考虑两个接触的颗粒,即颗粒i和j,mi和mj分别表示颗粒i、j的质量;ri和rj表示这些颗粒各自的半径;分别定义m*和r*为等效质量和半径,表示为m*=mimj/(mi+mj)和r*=rirj/(ri+rj),两个颗粒具有相同的材料性质,剪切模量为G,泊松比为v,当颗粒i与颗粒j以非零的相对速度斜向碰撞时,颗粒i与颗粒j之间存在接触力,作用在颗粒上的接触力包括法向分量和切向分量表示为:
各分量用弹簧力和阻尼力之和表示,表示为:
Fcn=Fcn,s+Fcn,d (4)
Fct=Fct,s+Fct,d (5)
在法向上,式(4)中的Fcn,s,Fcn,d分别为法向弹簧力和法向阻尼力;
Fct,s=Fct,s,(n-1)+ΔFct,s (8)
ΔFct,s和Fct,s,(n-1)分别为当前和之前时间步长的切向弹簧力,ΔFct,s为切向弹簧力的增量,ξt为切向阻尼比,kt为切向接触刚度,为相对切向速度,增量ΔFct,s进一步由ΔFct,s=-ksΔδt给出,Δδt为增量切向位移,表示为其中Δt为时间步长;
切向力受库仑摩擦定律限制,其表达式为:
|Fct|≤μFcn (11)
其中μ是摩擦系数,该接触力模型所需的参数包括剪切模量(G)、泊松比(v)、颗粒间的摩擦系数(μb)、颗粒与机械部件之间的摩擦系数(μw)、颗粒间的恢复系数(eb)、颗粒与机械部件之间的恢复系数(ew)。
6.根据权利要求5所述的颗粒阻尼动力吸振器耦合数值计算方法,其特征在于,所述步骤二中DEM循环的流程c中利用式(1)、(2)计算每个颗粒的线加速度和角加速度。
8.根据权利要求7所述的颗粒阻尼动力吸振器耦合数值计算方法,其特征在于,所述步骤三中MBD循环方案的流程具体如下:
流程a、计算广义力;
流程b、通过化简方程求解MBD的控制方程;
流程c、利用RK4方法更新机械零件的位置和速度;
流程d、一旦MBD计算在每个时间步ΔtMBD结束,回到步骤二,计算相应的平移和旋转的机械零件发送回DME循环,然后确定阻尼颗粒的运动行为和等效颗粒部分再次机械零件。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210845856.3A CN115081146B (zh) | 2022-07-19 | 2022-07-19 | 一种颗粒阻尼动力吸振器耦合数值计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210845856.3A CN115081146B (zh) | 2022-07-19 | 2022-07-19 | 一种颗粒阻尼动力吸振器耦合数值计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115081146A true CN115081146A (zh) | 2022-09-20 |
CN115081146B CN115081146B (zh) | 2024-07-19 |
Family
ID=83260297
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210845856.3A Active CN115081146B (zh) | 2022-07-19 | 2022-07-19 | 一种颗粒阻尼动力吸振器耦合数值计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115081146B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118133647A (zh) * | 2024-05-06 | 2024-06-04 | 华东交通大学 | 一种垂直分腔颗粒阻尼钢轨吸振器构建方法及系统 |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050183502A1 (en) * | 2004-02-23 | 2005-08-25 | Halliburton Energy Services, Inc. | Motion-responsive coupled masses |
JP2011081530A (ja) * | 2009-10-06 | 2011-04-21 | Univ Of Tokyo | 粒子モデルの離散要素法解析シミュレーション方法、離散要素法解析シミュレーションプログラム、及び離散要素法解析シミュレーション装置 |
US20140358505A1 (en) * | 2013-05-31 | 2014-12-04 | The Board Of Trustees Of The University Of Illinois | Collision impulse derived discrete element contact force determination engine, method, software and system |
CN104376153A (zh) * | 2014-11-03 | 2015-02-25 | 同济大学 | 一种颗粒阻尼器的优化设计方法 |
CN106598912A (zh) * | 2016-10-20 | 2017-04-26 | 浙江工业大学 | 一种基于cfd‑dem耦合模型的磨粒流场分析方法 |
US20170177769A1 (en) * | 2015-12-22 | 2017-06-22 | Dassault Systemes Simulia Corp. | Effectively Solving Structural Dynamics Problems With Modal Damping In Physical Coordinates |
CN110566619A (zh) * | 2019-09-20 | 2019-12-13 | 厦门振为科技有限公司 | 一种粒子阻尼声子晶体结构及其等效理论模型构建方法 |
CN112001017A (zh) * | 2020-08-25 | 2020-11-27 | 青岛理工大学 | 一种高频耗散可控的结构动力响应数值计算方法 |
CN112199883A (zh) * | 2020-07-31 | 2021-01-08 | 中国气象局气象探测中心 | 基于有效碰撞耗能计算的多颗粒阻尼器减振机理分析方法 |
CN112329168A (zh) * | 2020-11-02 | 2021-02-05 | 北京航空航天大学 | 一种二维振动条件下颗粒阻尼器减振效果的3d数值仿真方法 |
CN112464394A (zh) * | 2020-11-02 | 2021-03-09 | 北京航空航天大学 | 一种基于有效碰撞耗能计算的多颗粒阻尼器减振机理分析方法 |
CN113408217A (zh) * | 2021-06-28 | 2021-09-17 | 西北工业大学 | 一种适用于粗颗粒流的气粒混合流动数值计算方法 |
CN114297797A (zh) * | 2021-12-29 | 2022-04-08 | 西安交通大学 | 基于ann的燃机透平阻尼叶片结构等效刚度阻尼分析方法 |
CN114607720A (zh) * | 2022-03-17 | 2022-06-10 | 江苏科技大学 | 一种内置障碍网络的颗粒阻尼器 |
-
2022
- 2022-07-19 CN CN202210845856.3A patent/CN115081146B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050183502A1 (en) * | 2004-02-23 | 2005-08-25 | Halliburton Energy Services, Inc. | Motion-responsive coupled masses |
JP2011081530A (ja) * | 2009-10-06 | 2011-04-21 | Univ Of Tokyo | 粒子モデルの離散要素法解析シミュレーション方法、離散要素法解析シミュレーションプログラム、及び離散要素法解析シミュレーション装置 |
US20140358505A1 (en) * | 2013-05-31 | 2014-12-04 | The Board Of Trustees Of The University Of Illinois | Collision impulse derived discrete element contact force determination engine, method, software and system |
CN104376153A (zh) * | 2014-11-03 | 2015-02-25 | 同济大学 | 一种颗粒阻尼器的优化设计方法 |
US20170177769A1 (en) * | 2015-12-22 | 2017-06-22 | Dassault Systemes Simulia Corp. | Effectively Solving Structural Dynamics Problems With Modal Damping In Physical Coordinates |
CN106598912A (zh) * | 2016-10-20 | 2017-04-26 | 浙江工业大学 | 一种基于cfd‑dem耦合模型的磨粒流场分析方法 |
CN110566619A (zh) * | 2019-09-20 | 2019-12-13 | 厦门振为科技有限公司 | 一种粒子阻尼声子晶体结构及其等效理论模型构建方法 |
CN112199883A (zh) * | 2020-07-31 | 2021-01-08 | 中国气象局气象探测中心 | 基于有效碰撞耗能计算的多颗粒阻尼器减振机理分析方法 |
CN112001017A (zh) * | 2020-08-25 | 2020-11-27 | 青岛理工大学 | 一种高频耗散可控的结构动力响应数值计算方法 |
CN112329168A (zh) * | 2020-11-02 | 2021-02-05 | 北京航空航天大学 | 一种二维振动条件下颗粒阻尼器减振效果的3d数值仿真方法 |
CN112464394A (zh) * | 2020-11-02 | 2021-03-09 | 北京航空航天大学 | 一种基于有效碰撞耗能计算的多颗粒阻尼器减振机理分析方法 |
CN113408217A (zh) * | 2021-06-28 | 2021-09-17 | 西北工业大学 | 一种适用于粗颗粒流的气粒混合流动数值计算方法 |
CN114297797A (zh) * | 2021-12-29 | 2022-04-08 | 西安交通大学 | 基于ann的燃机透平阻尼叶片结构等效刚度阻尼分析方法 |
CN114607720A (zh) * | 2022-03-17 | 2022-06-10 | 江苏科技大学 | 一种内置障碍网络的颗粒阻尼器 |
Non-Patent Citations (11)
Title |
---|
夏兆旺;温华兵;刘献栋;: "颗粒阻尼器结构振动特性耦合算法仿真与试验", 农业机械学报, no. 08, 25 August 2011 (2011-08-25) * |
夏兆旺;魏守贝;温华兵;刘献栋;: "颗粒阻尼技术在制动鼓减振方面的应用研究", 振动工程学报, no. 06, 15 December 2014 (2014-12-15) * |
毛宽民,陈天宁,黄协清: "非阻塞性微颗粒阻尼的散体元模型", 西安交通大学学报, no. 05, 20 May 1999 (1999-05-20) * |
温华兵;左言言;彭子龙;郭夕军;: "船舶结构复合阻尼材料减振性能实验研究", 船舶工程, no. 04, 15 August 2013 (2013-08-15) * |
温华兵;彭子龙;刘林波;: "基于FE-SEA混合法的空心阻振质量阻振性能研究", 振动与冲击, no. 05, 15 March 2015 (2015-03-15) * |
温华兵;彭子龙;夏兆旺;刘林波;: "新型空心结构阻振质量的阻振性能实验", 振动.测试与诊断, no. 02, 15 April 2016 (2016-04-15) * |
温华兵;昝浩;陈宁;杨兴林;: "惯容器对隔振系统动态性能影响研究", 实验力学, no. 04, 19 August 2015 (2015-08-19) * |
王博;于哲舟;袁军;付宏;于建群;: "基于MBD和DEM耦合的新型CAE软件", 吉林大学学报(理学版), no. 02, 26 March 2020 (2020-03-26) * |
胡溧;黄其柏;何雪松;: "非阻塞性颗粒阻尼的离散元法", 华中科技大学学报(自然科学版), no. 07, 15 July 2008 (2008-07-15) * |
谭德昕;刘献栋;单颖春;: "颗粒阻尼减振器对板扭转振动减振的仿真研究", 系统仿真学报, no. 08, 8 August 2011 (2011-08-08) * |
赵月静;秦志英;廖冬梅;: "内置颗粒减振器的机床刀具减振分析", 机床与液压, no. 08, 28 April 2020 (2020-04-28) * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118133647A (zh) * | 2024-05-06 | 2024-06-04 | 华东交通大学 | 一种垂直分腔颗粒阻尼钢轨吸振器构建方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN115081146B (zh) | 2024-07-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115081146A (zh) | 一种颗粒阻尼动力吸振器耦合数值计算方法 | |
Sohoni et al. | Automatic linearization of constrained dynamical models | |
Ray et al. | Hysteretic models for sliding bearings with varying frictional force | |
CN110110364B (zh) | 基于双超卫星平台碰振动力学模型的防碰撞最优控制方法 | |
CN108983734A (zh) | 一种考虑三角形结构下遥操作系统的有限时间控制方法 | |
Song et al. | Second-order sliding mode control with backstepping for aeroelastic systems based on finite-time technique | |
Wang et al. | Improving ride comfort in vehicles with seat vibration isolator embedded with different negative stiffness models | |
Feng et al. | Modeling and robust control of horizontal vibrations for high-speed elevator | |
Ma et al. | A planar shock isolation system with high-static-low-dynamic-stiffness characteristic based on cables | |
Gao et al. | Adaptive vibration control of a flexible structure based on hybrid learning controlled active mass damping | |
Dohrmann | Dynamics of a tire–wheel–suspension assembly | |
Dong et al. | Gray skyhook predictive control of magnetorheological semi-active seat suspension with time delay | |
Gong et al. | Characterization of a novel variable friction connection for semiactive cladding system | |
Wang et al. | Energy‐based finite‐time stabilization and H∞ control of stochastic nonlinear systems | |
Lee | A short note for numerical analysis of dynamic contact considering impact and a very stiff spring-damper constraint on the contact point | |
Park et al. | Should friction losses be included in an electromechanical model of a bioinspired flapping-wing micro aerial vehicle to estimate the flight energetic requirements? | |
Hong et al. | A hybrid PID controller for flexible joint manipulator based on state observer and singular perturbation approach | |
Dias et al. | Dynamics of flexible mechanical systems with contact-impact and plastic deformations | |
Dai et al. | Random displacement and acceleration responses of vehicles with uncertainty | |
Wasfy et al. | Modeling the dynamic frictional contact of tires using an explicit finite element code | |
CN109684766B (zh) | 含转角的大变形柔性梁单元建模方法 | |
Gao et al. | Electrical network optimization based electrically interconnected suspension control for vehicle cabin | |
Zhu et al. | Multiscale analysis including strain rate dependency for transient response of composite laminated shells | |
Weber et al. | Dynamic characterization and isolation performance of triple friction pendulum | |
Gamar et al. | Application of l1 adaptive control for helicopter vibration |
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 |