CN113515868A - 一种磁流变弹性体的链簇模型构建方法及性能估测方法 - Google Patents
一种磁流变弹性体的链簇模型构建方法及性能估测方法 Download PDFInfo
- Publication number
- CN113515868A CN113515868A CN202110889330.0A CN202110889330A CN113515868A CN 113515868 A CN113515868 A CN 113515868A CN 202110889330 A CN202110889330 A CN 202110889330A CN 113515868 A CN113515868 A CN 113515868A
- Authority
- CN
- China
- Prior art keywords
- magnetorheological elastomer
- magnetic
- particle
- particles
- cluster model
- 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
- 229920001971 elastomer Polymers 0.000 title claims abstract description 155
- 239000000806 elastomer Substances 0.000 title claims abstract description 150
- 238000004138 cluster model Methods 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000010276 construction Methods 0.000 title description 2
- 239000002245 particle Substances 0.000 claims abstract description 132
- 239000006249 magnetic particle Substances 0.000 claims abstract description 56
- 230000003993 interaction Effects 0.000 claims abstract description 17
- 230000008859 change Effects 0.000 claims abstract description 11
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 6
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 claims description 6
- 239000010931 gold Substances 0.000 claims description 6
- 229910052737 gold Inorganic materials 0.000 claims description 6
- 230000035699 permeability Effects 0.000 claims description 6
- 239000007788 liquid Substances 0.000 claims description 3
- 229910052757 nitrogen Inorganic materials 0.000 claims description 3
- 238000005507 spraying Methods 0.000 claims description 3
- 230000005415 magnetization Effects 0.000 description 11
- 239000011159 matrix material Substances 0.000 description 11
- 230000000694 effects Effects 0.000 description 8
- 238000013016 damping Methods 0.000 description 6
- 230000004044 response Effects 0.000 description 6
- 238000003860 storage Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 5
- 239000006096 absorbing agent Substances 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000010008 shearing Methods 0.000 description 3
- 230000035939 shock Effects 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical group [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000001816 cooling 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
- 239000013078 crystal Substances 0.000 description 1
- 230000001351 cycling effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000004215 lattice model Methods 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 238000001179 sorption measurement Methods 0.000 description 1
- 239000000758 substrate Substances 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- 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/08—Thermal analysis or thermal optimisation
Abstract
发明公开了一种磁流变弹性体的链簇模型构建方法及性能估测方法，包括以下步骤：对磁流变弹性体样品进行微观结构扫描，获得磁流变弹性体样品中磁性颗粒的粒径均值和粒距均值；按照设定的磁性颗粒体积分数以随机磁性颗粒的粒径和粒距建立磁流变弹性体的链簇模型；将链簇模型中磁性颗粒视为强度相等的磁偶极子，根据两相邻磁偶极子的相互作用能计算总能量密度，再计算得到磁致剪切模量。本发明还公开了一种磁流变弹性体的性能估测方法，由磁流变弹性体中的磁性颗粒的体积分数以及应用中的温度变化代入所述链簇模型，获得所述磁流变弹性体的磁致剪切模量。本发明可解决磁流变弹性体预测性能与实际偏离较大的技术问题。
Description
技术领域
本发明涉及一种磁流变弹性体的模型构建方法及性能估测方法，特别是一种磁流变弹性体的链簇模型构建方法及性能估测方法。
背景技术
磁流变弹性体是近年来发展的一种新型磁流变材料，因可克服磁流变液颗粒沉降严重、稳定性差以及装置需密封等缺点而在高负载减振领域中颇受关注。磁性颗粒和橡胶基体是磁流变弹性体的两种重要组成，组成类型及比例的选取对磁流变弹性体的性能具有重要影响。在磁流变弹性体的制备过程中，磁场的施加使得磁性颗粒沿磁场定向排列并镶嵌于橡胶基体中，制备的磁流变弹性体呈各向异性。由于磁流变弹性体的力学性能可随外加磁场发生快速、可逆、连续的变化，基于其研制的减振装置可实现振动的半主动控制，在减振、降噪等实际工程应用中具有举足轻重的地位。作为磁流变弹性体减振装置的核心元件，磁流变弹性体性能的好坏直接决定磁流变弹性体减振装置的减振性能。
现有的磁流变弹性体微观模型包括磁偶极子模型，以及基于磁偶极子模型发展的耦合场模型、主链吸附模型、晶格模型、多链模型和有限柱模型等。但是上述模型中磁性颗粒以及磁链结构都呈理想分布，与实际偏差较大。此外，在磁流变弹性体的实际应用过程中，电流循环加载会产生大量的热。当磁流变弹性体减振器件的工作时间较长时，热量的大量积累会导致整个器件温度急剧升高，对磁流变弹性体的性能产生影响，进而直接影响磁流变弹性体减振器件的减振性能。目前已有磁流变弹性体微观模型预测的磁流变弹性体性能与实际偏离较大，难以准确获得磁流变弹性体随温度变化的原因，不利于磁流变弹性体的配方设计。
发明内容
针对上述现有技术缺陷，本发明的任务在于提供一种磁流变弹性体的链簇模型构建方法，以及提供一种磁流变弹性体的性能估测方法，解决磁流变弹性体预测性能与实际偏离较大的技术问题。
本发明技术方案如下：一种磁流变弹性体的链簇模型构建方法，包括以下步骤：
步骤一、对磁流变弹性体样品进行微观结构扫描，获得磁流变弹性体样品中磁性颗粒的粒径均值和粒距均值；
步骤二、按照设定的磁性颗粒体积分数以随机磁性颗粒的粒径和粒距建立磁流变弹性体的链簇模型，所述链簇模型的磁性颗粒的粒径均值以及粒距均值与所述步骤一中获得的磁流变弹性体样品中磁性颗粒的粒径均值和粒距均值相等；
步骤三、将链簇模型中磁性颗粒视为强度相等的磁偶极子，两相邻磁偶极子的相互作用能表示为
式中，H和hr分别为两磁偶极子间的距离以及沿磁偶极子连接方向的单位矢量；μ0和μ1分别为真空磁导率和磁性颗粒的相对磁导率，m是磁偶极矩，i表示中心磁偶极子，j表示围绕在中心磁偶极子周边的磁偶极子，总能量Et为单个链簇中所有粒子与围绕在其周边粒子的相互作用能总和，
与一维剪切应变相关的总能量密度表示为
式中，N是单链中的粒子数，V是磁流变弹性体的体积，φ和d分别是颗粒的体积分数和直径，
磁致剪切模量为剪切应力与剪切应变的比值表示为
式中，γ是剪切应变，σ是剪切应力即颗粒间能量密度相对于剪切应变的导数，
进一步地，包括步骤四、通过磁性颗粒粒径和粒距的联合概率密度函数修正磁致剪切模量，修正后的磁致剪切模量表示为
Gm＝∫∫G·f(d,h)dd dh
f(d,h)为磁性颗粒粒径和粒距的联合概率密度函数，d和h分别代表粒径和粒距，μd和μh分别为粒径和粒距的均值，σd和σh分别为粒径和粒距的标准差，ρ为粒径和粒距的相关系数。
进一步地，所述步骤一中对磁流变弹性体样品进行微观结构扫描具体包括步骤：首先将磁流变弹性体样品在液氮中脆断，以保证磁流变弹性体内部结构的完整性；然后，对脆断的磁流变弹性体样品进行表面喷金，以增强样品的导电性；最后，将喷金的磁流变弹性体样品置于电子扫描显微镜中进行结构扫描。
进一步地，所述步骤二建立磁流变弹性体的链簇模型时随机生成若干组(c,d)值并在该尺寸的矩形区域内随机生成粒径为e和粒距为f的粒子，按生成的(c,d)值调整矩形区域尺寸直至矩形区域内粒子体积分数达到所述设定的磁性颗粒体积分数。
进一步地，所述按生成的(c,d)值调整矩形区域尺寸时使单个所述随机生成的粒子不被矩形区域边界所切割。
进一步地，所述总能量Et表示为
一种磁流变弹性体的性能估测方法，包括步骤：根据上述磁流变弹性体的链簇模型构建方法构建磁流变弹性体的链簇模型后由磁流变弹性体中的磁性颗粒的体积分数以及应用中的温度变化代入所述链簇模型，获得所述磁流变弹性体的磁致剪切模量。
本发明与现有技术相比的优点在于：
本发明在传统理想磁流变弹性体微观模型基础上，引入粒径、粒距以及温度影响因子对模型进行修正，使构建的磁流变弹性体链簇模型更准确、更符合实际。本发明可更准确地预测磁流变弹性体的真实性能及影响因素，对磁流变弹性体的配方设计与工艺优化具有重要的指导意义。
附图说明
图1为磁流变弹性体的链簇模型示意图。
图2为磁流变弹性体链簇模型中中心粒子与相邻颗粒作用示意图。
图3为磁性颗粒的体积分数对磁流变弹性体最大磁致模量关系示意图。
图4为基于本发明磁流变弹性体的链簇模型得到的磁流变弹性体的磁响应特性示意图。
图5为基于本发明磁流变弹性体链簇模型得到的磁流变弹性体磁化特性随温度的变化示意图。
图6为磁流变弹性体基体材料的储能模量随温度的变化示意图。
具体实施方式
下面结合实施例对本发明作进一步说明，但不作为对本发明的限定。
磁流变弹性体的链簇模型构建方法以及通过链簇模型构建估算磁流变弹性体性能的方法，包括步骤一、对磁流变弹性体样品进行微观结构扫描，具体扫描方法为：首先将磁流变弹性体样品在液氮中脆断，以保证磁流变弹性体内部结构的完整性；然后，对脆断的磁流变弹性体样品进行表面喷金，以增强样品的导电性；最后，将喷金的磁流变弹性体样品置于电子扫描显微镜中，对仪器进行调试并设定相关参数后，对样品进行结构扫描，获得磁流变弹性体的微观形貌。
步骤二、按照设定的磁性颗粒体积分数以随机磁性颗粒的粒径和粒距建立磁流变弹性体的链簇模型。请结合图1所示，基于磁流变弹性体的微观结构，可将其看作无数个周期性排列的含链簇结构的代表性体积单元(RVE)。为了简化计算过程，在磁流变弹性体链簇模型构建时以其中的一个二维代表性体积单元作为研究对象。选取M条磁链作为一个链簇，以两个无相互作用的链簇作为磁流变弹性体的链簇模型。
磁流变弹性体的链簇模型构建步骤为：首先，基于磁流变弹性体的实际微观结构测试其中磁性颗粒的粒径和粒距大小；然后，分析磁性颗粒的粒径和粒距分布情况，绘制粒径和粒距分布图；再次，对磁性颗粒的粒径和粒距分布进行检验，确定粒径和粒距分布函数，编制程序生成磁流变弹性体的链簇模型。具体是：设定磁性颗粒的平均粒径和粒距分别为m和n，体积分数为V，链簇模型的尺寸为a×b。随机生成若干组(c,d)值并在该尺寸的矩形区域内随机生成粒径为e和粒距为f的粒子。若颗粒的体积分数未达到设定值，则不断调整模型尺寸大小直到达到颗粒的设定体积分数。经过不断调整，获得满足要求的链簇模型尺寸为A×B，进而得到单链中磁性颗粒的数量为N。
步骤三、请结合图2所示，磁流变弹性体链簇模型中任一链簇中中心粒子与相邻颗粒的不同作用方式。磁流变弹性体在磁场条件下的变刚度特性主要源于其中磁性颗粒的磁相互作用。在提出的磁流变弹性体链簇模型中，鉴于颗粒在较大间距下磁相互作用很弱这一特性，只考虑了颗粒与其相邻粒子之间的相互作用。由于磁流变弹性体链簇模型中两个链簇对称，且链簇中边缘行颗粒和中心行颗粒的粒子间相互作用不同，因此选取其中的一个链簇作为研究对象，并依据粒子所处环境对其链簇能量进行计算。对于中心行粒子(i，j1和j2)，中心粒子i与相邻颗粒的作用如图2中虚线框图部分所示，图2(a～d)为中心粒子i(1＜i＜N)与相邻粒子j1～8的不同相互作用方式。
将磁性颗粒视为强度相等的磁偶极子，设其强度为|m|，则两相邻磁偶极子的相互作用能表示为式(1)：
式中，H和hr分别为两磁偶极子间的距离以及沿磁偶极子连接方向的单位矢量；μ0和μ1分别为真空磁导率和磁性颗粒的相对磁导率。m是磁偶极矩，|m|＝JpVi，Jp是极化强度，Vi是粒子的体积，表示为
当施加垂直于铁链的剪应力τ时，产生小的应变γ。考虑到结构的对称性，仅计算了相互作用能Eij8，Eij2，Eij5和Eij6。如果两个偶极子在同一粒子链中，即i,j8(如图2(a)所示)，相互作用能为式(2)：
如果两个偶极子沿着磁场的垂直方向对齐，即i,j2(如图2(b)所示)，相互作用能表示为等式(3):
对于两个偶极子既不平行也不垂直对准磁场的情况，如i,j5和i,j6(如图2(c)和图2(d)所示)，其相互作用能的表达式不是唯一的。对于粒子i,j5，相互作用能可写为式(4):
对于粒子i,j6，相互作用能可写为式(5):
综上，粒子i与相邻颗粒作用的总能量可表示为式(6)：
同理，粒子j1与相邻颗粒作用的总能量可表示为式(7)：
粒子j2与相邻颗粒作用的总能量可表示为式(8)：
粒子i，j1和j2与相邻颗粒作用的总能量可表示为式(9)：
对于链簇中上边缘行和下边缘行的粒子，总能量均为式(10)：
因此，整个链簇的总能量为式(11)：
Et＝(N-2)Ec+E1+EN (11)
与一维剪切应变相关的总能量密度(单位体积能量)表示为式(12)：
式中，V是磁流变弹性体的体积，φ和d分别是颗粒的体积分数和直径，单链中磁性颗粒的数量N通过模型优化获得。剪切应力σ是颗粒间能量密度相对于剪切应变的导数，如式(13)所示：
磁致剪切模量为剪切应力与剪切应变的比值，表示为式(14)：
当磁性颗粒达到饱和状态时，M＝Ms，Jp＝μ0Ms，磁流变弹性体的磁致模量达到最大。
步骤四、鉴于磁性颗粒粒径和粒距的相互影响较大，采用联合概率密度函数来表示二者的概率密度分布，相关系数ρ由相关性检验获得，式(15)为磁性颗粒粒径和粒距的联合概率密度函数表达式：
式中，d和h分别代表粒径和粒距，μd和μh分别为粒径和粒距的均值，σd和σh分别为粒径和粒距的标准差，ρ为粒径和粒距的相关系数。
基于制备的磁流变弹性体中磁性颗粒的粒径和粒距分布函数并结合式(15)，可获得磁性颗粒粒径和粒距的联合概率密度函数表达式。将其带入式(14)，获得修正的磁流变弹性体磁致剪切模量表达式，如式(16)所示：
Gm＝∫∫G·f(d,h)dd dh (16)
式中，磁偶极矩|m|可用磁化强度M表示，如式(17)所示：
|m|＝μ0MVi (17)
由于磁性颗粒的磁化强度M同时是磁场H和温度T的函数，结合式(16)的表达式可以看出，磁场强度H、温度T、颗粒体积分数φ对磁流变弹性体磁致剪切模量有重要的影响。
下面结合实施例验证该模型的合理性。
磁流变弹性体减振器是一种新型智能减振器，其减振效果的好坏取决于核心减振材料磁流变弹性体的磁致性能。制备配方对磁流变弹性体的性能起决定性作用，磁性颗粒的体积分数是影响磁流变弹性体磁致性能的关键因素。请结合图3所示，磁性颗粒的体积分数对磁流变弹性体最大磁致模量影响的数值拟合结果，可以看出，当在原始基础上连续添加5％体积的磁性颗粒时，磁流变弹性体最大磁致剪切模量Gmax的增加百分比逐渐降低，当磁性颗粒的体积分数为30％时，磁流变弹性体的综合性能较好，最大磁致剪切模量可达0.120MPa。
根据上述数值计算结果，在磁流变弹性体制备时，将磁性颗粒体积分数定为30％。对制备的磁流变弹性体采用旋转流变仪进行磁场扫描和温度扫描获得其磁响应特性和温度特性。磁流变弹性体磁致模量Gm和相对磁流变效应(MR effect)由式(13)和(14)计算：
Gm＝Gs-G0 (18)
式中，G0和Gs分别代表零场和磁饱和状态下磁流变弹性体的剪切储能模量。磁流变效应值越大，磁流变弹性体的磁响应特性越佳。请结合图4所示，磁流变弹性体的磁响应特性。图4(a)为电流与磁场的关系，图4(b)为磁流变弹性磁致模量随电流的变化。从图中可以看出，磁流变弹性体的磁致模量Gm最大可达0.118MPa，略低于理论计算值。这可能是由于实际磁流变弹性体中磁性颗粒链并非严格直链引起的。
在磁流变弹性体的应用过程中，温度的变化可能引起磁流变弹性体磁化特性或基体材料性能的变化。为了阐明磁流变弹性体性能随温度变化的根本原因，分别对磁流变弹性体的磁化特性和基体材料性能随温度的变化进行了研究。
结合图5所示，磁流变弹性体磁化特性随温度的变化。图5(a)为磁流变弹性体的最大磁致剪切模量Gmax随温度变化的数值计算结果，可以看出，磁流变弹性体的磁化特性受温度影响很小。随温度的增加，磁流变弹性体的最大磁致剪切模量Gmax仅减小3.31％(由0.121MPa变为0.117MPa)。图5(b)为实验测得的磁流变弹性体的磁温曲线，同样可以看出，磁流变弹性体的磁化特性受温度影响很小。随温度的增加，磁流变弹性体的最大磁化强度仅减小1.93％(196.5KA/m减小为192.7KA/m)。综上所述，数值计算结果和实验测试结果都表明，温度对磁流变弹性体的磁化特性影响较小。
结合图6所示，温度对磁流变弹性体基体材料储能模量的影响。从中可以看出，橡胶基体的储能模量在升温过程中逐渐降低，降低百分比为25.5％。在降温过程中，橡胶基体的储能模量逐渐升高，增加百分比为35.3％。因此，温度对磁流变弹性体基体材料的性能具有重要影响。
由温度对磁流变弹性体磁化特性和橡胶基体储能模量的影响可以得出，磁流变弹性体性能随温度的变化主要源于橡胶基体性能的变化。因此，在实际磁流变弹性体的应用过程中要根据磁流变弹性体的实际工作环境选取满足需求的橡胶基体。
综上，磁流变弹性体的磁响应特性和温度响应特性研究结果都表明，本发明构建的磁流变弹性体的链簇模型合理。该模型不仅可以预测磁流变弹性体在不同条件下的性能，还可指导磁流变弹性体的配方设计和工艺优化。通过调节磁流变弹性体的影响因子，可获得满足不同需求的磁流变弹性体，进而用于不同领域的振动控制。
Claims (7)
1.一种磁流变弹性体的链簇模型构建方法，其特征在于，包括以下步骤：
步骤一、对磁流变弹性体样品进行微观结构扫描，获得磁流变弹性体样品中磁性颗粒的粒径均值和粒距均值；
步骤二、按照设定的磁性颗粒体积分数以随机磁性颗粒的粒径和粒距建立磁流变弹性体的链簇模型，所述链簇模型的磁性颗粒的粒径均值以及粒距均值与所述步骤一中获得的磁流变弹性体样品中磁性颗粒的粒径均值和粒距均值相等；
步骤三、将链簇模型中磁性颗粒视为强度相等的磁偶极子，两相邻磁偶极子的相互作用能表示为
式中，H和hr分别为两磁偶极子间的距离以及沿磁偶极子连接方向的单位矢量；μ0和μ1分别为真空磁导率和磁性颗粒的相对磁导率，m是磁偶极矩，i表示中心磁偶极子，j表示围绕在中心磁偶极子周边的磁偶极子，总能量Et为单个链簇中所有粒子与围绕在其周边粒子的相互作用能总和，
与一维剪切应变相关的总能量密度表示为
式中，N是单链中的粒子数，V是磁流变弹性体的体积，φ和d分别是颗粒的体积分数和直径，
磁致剪切模量为剪切应力与剪切应变的比值表示为
式中，γ是剪切应变，σ是剪切应力即颗粒间能量密度相对于剪切应变的导数，
2.根据权利要求1所述的磁流变弹性体的链簇模型构建方法，其特征在于，包括步骤四、通过磁性颗粒粒径和粒距的联合概率密度函数修正磁致剪切模量，修正后的磁致剪切模量表示为
Gm＝∫∫G·f(d,h)dd dh
f(d,h)为磁性颗粒粒径和粒距的联合概率密度函数，d和h分别代表粒径和粒距，μd和μh分别为粒径和粒距的均值，σd和σh分别为粒径和粒距的标准差，ρ为粒径和粒距的相关系数。
3.根据权利要求1所述的磁流变弹性体的链簇模型构建方法，其特征在于，所述步骤一中对磁流变弹性体样品进行微观结构扫描具体包括步骤：首先将磁流变弹性体样品在液氮中脆断；然后对脆断的磁流变弹性体样品进行表面喷金；最后将喷金的磁流变弹性体样品置于电子扫描显微镜中进行结构扫描。
4.根据权利要求1所述的磁流变弹性体的链簇模型构建方法，其特征在于，所述步骤二建立磁流变弹性体的链簇模型时随机生成若干组(c,d)值并在该尺寸的矩形区域内随机生成粒径为e和粒距为f的粒子，按生成的(c,d)值调整矩形区域尺寸直至矩形区域内粒子体积分数达到所述设定的磁性颗粒体积分数。
5.根据权利要求4所述的磁流变弹性体的链簇模型构建方法，其特征在于，所述按生成的(c,d)值调整矩形区域尺寸时使单个所述随机生成的粒子不被矩形区域边界所切割。
6.根据权利要求1所述的磁流变弹性体的链簇模型构建方法，其特征在于，所述总能量Et表示为
7.一种磁流变弹性体的性能估测方法，其特征在于，包括步骤：根据权利要求1至6任意一项所述的磁流变弹性体的链簇模型构建方法构建磁流变弹性体的链簇模型后，由磁流变弹性体中的磁性颗粒的体积分数以及应用中的温度变化代入所述链簇模型获得所述磁流变弹性体的磁致剪切模量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110889330.0A CN113515868A (zh) | 2021-08-04 | 2021-08-04 | 一种磁流变弹性体的链簇模型构建方法及性能估测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110889330.0A CN113515868A (zh) | 2021-08-04 | 2021-08-04 | 一种磁流变弹性体的链簇模型构建方法及性能估测方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113515868A true CN113515868A (zh) | 2021-10-19 |
Family
ID=78068088
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110889330.0A Pending CN113515868A (zh) | 2021-08-04 | 2021-08-04 | 一种磁流变弹性体的链簇模型构建方法及性能估测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113515868A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114042912A (zh) * | 2021-11-12 | 2022-02-15 | 哈尔滨工业大学 | 一种通过粉末粒径精细化控制NiAl基复合材料力学性能的方法 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008147080A1 (en) * | 2007-05-25 | 2008-12-04 | Sk Energy Co., Ltd. | Magnetic composite particles and magnetorheological fluid using the same |
| US20080318045A1 (en) * | 2004-08-27 | 2008-12-25 | Fraunhofer-Gesellschaft Zur Forderung Der Angewandten Forschung E.V. | Magnetorheological Elastomers and Use Thereof |
| CN102175572A (zh) * | 2011-02-24 | 2011-09-07 | 同济大学 | 磁流变液动态屈服变异性标定的微观尺度方法 |
| EP2541230A1 (en) * | 2011-06-30 | 2013-01-02 | Koninklijke Philips Electronics N.V. | Detection of clusters of magnetic particles |
| US20130060349A1 (en) * | 2011-09-06 | 2013-03-07 | Ossur Hf | Prosthetic and orthotic devices having magnetorheological elastomer spring with controllable stiffness |
| CN103435822A (zh) * | 2013-07-05 | 2013-12-11 | 西南交通大学 | 一种用脉动磁场制备磁流变弹性体的方法 |
| CN103909468A (zh) * | 2014-04-04 | 2014-07-09 | 常熟理工学院 | 基于磁流变弹性体的并联力控制装置 |
| CN112395759A (zh) * | 2020-11-16 | 2021-02-23 | 同济大学 | 一种磁流变阻尼器建模方法 |
-
2021
- 2021-08-04 CN CN202110889330.0A patent/CN113515868A/zh active Pending
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080318045A1 (en) * | 2004-08-27 | 2008-12-25 | Fraunhofer-Gesellschaft Zur Forderung Der Angewandten Forschung E.V. | Magnetorheological Elastomers and Use Thereof |
| WO2008147080A1 (en) * | 2007-05-25 | 2008-12-04 | Sk Energy Co., Ltd. | Magnetic composite particles and magnetorheological fluid using the same |
| CN102175572A (zh) * | 2011-02-24 | 2011-09-07 | 同济大学 | 磁流变液动态屈服变异性标定的微观尺度方法 |
| EP2541230A1 (en) * | 2011-06-30 | 2013-01-02 | Koninklijke Philips Electronics N.V. | Detection of clusters of magnetic particles |
| US20130060349A1 (en) * | 2011-09-06 | 2013-03-07 | Ossur Hf | Prosthetic and orthotic devices having magnetorheological elastomer spring with controllable stiffness |
| CN103435822A (zh) * | 2013-07-05 | 2013-12-11 | 西南交通大学 | 一种用脉动磁场制备磁流变弹性体的方法 |
| CN103909468A (zh) * | 2014-04-04 | 2014-07-09 | 常熟理工学院 | 基于磁流变弹性体的并联力控制装置 |
| CN112395759A (zh) * | 2020-11-16 | 2021-02-23 | 同济大学 | 一种磁流变阻尼器建模方法 |
Non-Patent Citations (2)
| Title |
|---|
| 朱俊涛等: "磁流变弹性体主链吸附模型", 东南大学学报(自然科学版), no. 02, 20 March 2011 (2011-03-20), pages 1 * |
| 袁飞洋等: "磁流变弹性体力磁耦合本构关系的研究进展", 材料导报, no. 19, 10 October 2017 (2017-10-10) * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114042912A (zh) * | 2021-11-12 | 2022-02-15 | 哈尔滨工业大学 | 一种通过粉末粒径精细化控制NiAl基复合材料力学性能的方法 |
| CN114042912B (zh) * | 2021-11-12 | 2022-07-29 | 哈尔滨工业大学 | 一种通过粉末粒径精细化控制NiAl基复合材料力学性能的方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kazaryan et al. | Generalized phase-field model for computer simulation of grain growth in anisotropic systems | |
| Li et al. | A state-of-the-art review on magnetorheological elastomer devices | |
| Pham et al. | Cyclic deformation response of AISI 316L at room temperature: Mechanical behaviour, microstructural evolution, physically-based evolutionary constitutive modelling | |
| Koo et al. | Dynamic characterization and modeling of magneto-rheological elastomers under compressive loadings | |
| Sitharam et al. | Micromechanical modelling of monotonic drained and undrained shear behaviour of granular media using three‐dimensional DEM | |
| Radjai et al. | Fabric evolution and accessible geometrical states in granular materials | |
| Yang et al. | Performance tests and microstructure‐based sigmoid model for a three‐coil magnetorheological damper | |
| CN111666671A (zh) | 一种围岩体蠕变参数的实时反演方法 | |
| CN113515868A (zh) | 一种磁流变弹性体的链簇模型构建方法及性能估测方法 | |
| Lu et al. | A new invertible model of magnetorheological damper based on sigmoid function | |
| Delannay et al. | Predictions of inter-granular cracking and dimensional changes of irradiated polycrystalline graphite under plane strain | |
| Bao et al. | Evolution of the piezomagnetic field of ferromagnetic steels subjected to cyclic tensile stress with variable amplitudes | |
| Zang et al. | The performance of nonlinear vibration control via NiTiNOL–Steel wire ropes | |
| CN109186967B (zh) | 一种基于bp神经网络的橡胶减振器性能预测及选型方法 | |
| Erkus et al. | Dissipativity analysis of the base isolated benchmark structure with magnetorheological fluid dampers | |
| Carmichael | Remanent and transitory effects of elastic deformation of magnetite crystals | |
| Kolesnikova et al. | Magnetoactive elastomers based on ferromagnetic and ferroelectric particles: A FORC approach | |
| Dinesh et al. | Dynamic properties and liquefaction behaviour of granular materials using discrete element method | |
| Wang et al. | Investigation of a magnetorheological elastomer metamaterial sandwich beam with tunable graded stiffness for broadband vibration attenuation | |
| Singh et al. | Hysteretic response of bulk magnetostrictive material employing a novel hyperbolic vector generalized magneto-thermoelastic constitutive model | |
| Kiarie et al. | Modeling of magnetic properties of magnetorheological elastomers using JA hysteresis model | |
| Wang et al. | Modelling of hysteresis in smart actuators using the generalised Prandtl-Ishlinskii operator | |
| Deng et al. | Nonlinear dynamics of magnetorheological whole-satellite with variable parameters under small amplitude and medium-high frequency vibration | |
| Khan et al. | Investigation on the performance of MR damper with various piston configurations | |
| Li et al. | Fabrication and characterisation of patterned magnetorheological elastomers |
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 |