首页/文章/ 详情

发电机防振橡胶支座的振动分析

19天前浏览844
摘要:本研究旨在测量电机工作周期内由振动引起的防振橡胶支座的振动和应力。通过使用加速度计(压电传感器)测量发动机壳体的加速度来实现这一目标。应力分析则采用ANSYS 15.0软件进行有限元数值计算。
关键词:防振橡胶支座、应力分析、发电机
  1. 1.引言

很多国家电力资源有限,家庭使用小型发电机来满足家庭用电需求。一般而言,所有类型的发电机都由两部分组成:负责将机械功转换为电能的部分,以及负责燃烧燃料以产生足够能量驱动轴转动进而带动电气部分的机械部分(即内燃机)。

防震橡胶支座示意图

发动机在工作过程中通常伴随着大量振动,制造商倾向于通过将整个发动机置于防振橡胶支座上来消除这些振动,如图1所示。然而,由于持续振动产生的动态负荷,这些支座在一段时间后会发生失效。因此,有必要研究这些橡胶部件中产生的应力。本研究采用可靠高效的有限元分析(FEA)方法进行分析。研究所用的橡胶支座为V形发电机支座,如图2所示,这种类型在当地市场较为常见。

V型防震橡胶支座

  1. 2.实验工作

为了进行问题分析,需要测量V形发电机支座中橡胶的性质,以便代入有限元分析。橡胶支座承载质量为m的发电机并吸收传递到底座框架的发动机振动,因此可以将其表示为质量-弹簧-阻尼系统,如图3所示。该系统包括自由振动,橡胶件的刚度系数K和阻尼系数C通过实验计算得出[1]。

质量-弹簧-阻尼系统

  1. 3.刚度系数测量

从V形发电机支座中取出橡胶材料,通过压缩试验来估算刚度,公式如下:
k = F / X ......(1)
其中k为刚度系数,F为施加的载荷(单位:N),X为变形量(单位:m)。图4展示了计算机记录数据的压缩装置。

橡胶压缩测试仪

当施加的力为0.06 kN时,在伸长量等于0.1 mm处计算k值。计算得出刚度等于600000 N/m。
  1. 4.压电加速度计

本研究使用的加速度传感器为B&K 4371型加速度计,如图5所示。它可测量一个方向的加速度。该传感器测得的数据单位为毫伏,要转换为m/s²,需将数据乘以10.2。

加速度传感器

  1. 5.阻尼系数测量

阻尼系数C通过以下公式计算:
C/M = 2ζωn ......(2)
其中M为质量,ζ为阻尼比,ωn为自然频率。图6所示的实验用于估算公式(2)中的参数。该实验模拟了图3所示的系统,其中V形发电机支座中的橡胶材料被释放并与质量为1 kg的质量块连接。加速度传感器安装在质量块上。完整的测试装置如图6和图7所示。

首先需要计算阻尼比。该实验测得的加速度数据存在噪声,因此使用Sigview软件包进行处理以消除噪声,然后乘以压电加速度计的校正因子10.2,最终信号如图8所示。
随后对加速度曲线进行两次积分,得到速度和位移曲线。

加速度-时间曲线

通过Sigview软件包发现,图10中第一个峰值X1 = 1.5587×10⁻⁵,第二个峰值X2 = 9.2479×10⁻⁶,因此可以通过以下公式计算阻尼比ζ:

其中Δ为对数衰减率。将X1和X2的值代入公式(3)并求解ζ,得出阻尼比等于ζ = 0.462。公式(2)中还需要第二个参数自然频率,可通过以下公式计算:

得出ωn = 774.59 Hz
将ζ和ωn的值代入公式(2),可计算得出阻尼系数C = 715.73 N·s/m。
  1. 6.初始载荷测量

在发电机运行过程中,测量机械和电气部件的加速度,以确定作用于V形发电机支座的力。使用天平称量机械和电气部件的质量,假设它们为一个整体。测得质量m = 32 kg。

因此,作用于V形支座的力可通过以下公式计算:
F = m × a ......(6)
其中F为作用于发电机橡胶支座的激励力。
如图11所示,使用加速度计测量机械和电气部件在x、y、z三个轴向的加速度。
压电传感器获得的数据单位为毫伏,首先乘以10.2转换为m/s²,然后与发电机机械和电气部件的质量值一起代入公式(6),以获得特定方向的激励力。图12、13和14分别显示了时域内x、y、z方向的加速度。

X方向加速度时域曲线

Y方向加速度时域曲线

Z方向加速度时域曲线

将橡胶性质和作用于支座的动态载荷作为输入数据用于ANSYS 15.0的有限元分析。ANSYS中的瞬态响应分析需要频域而非时域的力数据输入。因此,需要将测得的加速度从时域转换到频域。为此,编写了MATLAB程序并使用快速傅里叶变换(FFT)[2]进行转换。

X方向加速度频域曲线

Y方向加速度频域曲线

Z方向加速度频域曲线

转换到频域的结果如图15、16和17所示,使用高达2000 Hz的频率数据来估算x、y、z方向的远程力。
  1. 7.有限元建模

使用Design Modeler为ANSYS求解器建立了一个质量块放置在四个V形防振橡胶支座上的发电机模型,支持2D和3D绘图。除了假设为与实际发电机机械和电气部件整体质量相等的质量块外,所有尺寸均取自真实样品。图18展示了本研究使用的Astra型发电机。图19展示了发电机的3D模型。

  1. 8.防振橡胶支座的粘弹性行为

发电机支座中的橡胶材料表现出粘弹性材料的特性,这一术语通常用于振动阻尼问题。它们是同时具有弹性和粘性行为的非线性材料[3]。模拟粘弹性最简单的形式是Kelvin-Voigt模型[4],该模型通过并联或串联连接阻尼器和弹簧获得。本文选择了并联模式,如图20所示。

  1. 9.防振橡胶支座材料选择

橡胶被视为粘弹性材料,其在ANSYS中设置的属性是通过实验确定的,在模型中为每个橡胶支座定义了弹簧单元,如图21所示。属性列于下表:

表1. 橡胶属性
泊松比    | 密度(kg/m3) |  K (N/m) | C (N·s/m)
        0.3      |      1267        |   600000  | 715.73         
结构钢的材料属性设置如表2所示:
表2. 标准结构钢属性
杨氏模量(GPa) | 泊松比 | 密度(kg/m³)
   200        |     0.3   | 7800     
  1. 10.结果与讨论

使用ANSYS 15.0中的谐响应求解器分析防振橡胶支座。网格尺寸设为0.002 m,节点数为78295,单元数为28376。如图22所示,在整个四个支座上施加力载荷,频率范围为10-2000 Hz。图22中的红色 区域表示力的作用区域。

图23显示了运行过程中产生的von Mises应力,图24展示了最大主应力。图25和图26分别显示了von Mises应变和最大主应变。研究发现, 防振橡胶支座中产生的应力小于支座材料的极限应力。




来源:ABAQUS仿真世界
Design Modeler振动非线性燃烧电力MATLAB电机材料试验ANSYS电气
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-11-03
最近编辑:19天前
yunduan082
硕士 | 仿真主任工程... Abaqus仿真世界
获赞 153粉丝 218文章 309课程 0
点赞
收藏
作者推荐

论文推送|复合材料层压板三维渐进失效的落锤冲击建模

摘要复合材料层压板由于缺乏厚度方向的增强,容易受到面外冲击载荷的影响。与高速冲击导致的局部损伤不同,低速冲击损伤涉及多种失效机制,如基体开裂、纤维断裂和广泛的界面分层。损伤程度的不同可能导致结构承载能力的显著降低。使用可靠的渐进损伤分析工具预测复合材料的抗冲击损伤性能,对减少结构级别的大量昂贵认证试验至关重要。本研究通过ABAQUS的VUMAT实现了一种增强型显式三维损伤模型,用于模拟[454/04/-454/904]sHexplyAS4/8552复合材料层压板的落锤冲击。该模型采用三维连续损伤模型(CDM)耦合能量驱动失效机制来捕捉冲击引起的损伤及其范围。该模块在显式求解器中提供了统一的解决方案,可预测冲击响应及随后的压缩残余强度。选择了两个案例来展示动态和静态载荷下渐进失效分析的能力:1)落锤试验;2)开孔拉伸试验。将开发的VUMAT数值预测结果与试验数据和NASA开源CompDam代码的预测结果进行了比较。1、引言层压复合材料易受面外冲击载荷的影响。复合材料在冲击下可能表现出多种失效机制,包括各种层内和层间损伤。这些失效模式的起始和发展取决于诸如冲击体特性、冲击能量、复合材料铺层、结构配置或环境条件等各种冲击参数。低速冲击通常会导致低于勉强可见冲击损伤(BDIV)的损伤。基体开裂、分层和纤维断裂等不可察觉的冲击损伤可能导致其压缩残余强度显著降低。因此,了解复合材料的冲击行为及相关冲击引起的损伤至关重要。高保真数值工具可以有效地作为冲击分析的虚拟试验,并有助于减少试验矩阵。不同的方法已被用于在细观尺度上模拟层压复合材料的低速冲击,考虑均质化的复合材料层和界面分层。González等人进行了复合材料层压板的落锤冲击试验和冲击后压缩试验,并使用连续损伤力学(CDM)方法对试验进行了模拟。该CDM模型基于LaRC04失效准则,忽略了面外应力。Donadon等人提出了一种基于LS-DYNA的显式渐进模型,采用最大应变准则来表征纤维拉伸和压缩、基体拉伸和面内/面外剪切模式下复合材料层压板的冲击载荷。Bouvet等人使用零厚度界面单元模拟了基体开裂,并通过考虑开裂基体的塑性来解释冲击后的永久压痕。Aymerich等人使用界面内聚单元预测了交叉铺层层压板中冲击引起的分层。内聚单元作为Abaqus的用户材料子程序(VUMAT)实现,并采用损伤准则来包括面外压缩对分层起始预测的影响。Feng等人扩展了前人的工作,在VUMAT中结合了层内损伤和层间损伤的本构模型,并获得了与落锤冲击试验结果的良好相关性。在Zhang等人早期的工作中,潜在的分层和基体开裂区域通过与接触界面相关的损伤演化内聚区来表征。该损伤模型作为用户子程序VUINTER实现,并对复合材料结构的低速冲击损伤进行了良好的预测。尽管在表征冲击损伤和失效预测性能方面已经开展了大量工作,但探索一种高效准确的渐进损伤分析(PDA)工具来模拟冲击响应和预测冲击事件后的残余强度性能一直是持续的研究兴趣。实现收敛解决方案并减少计算资源需求仍然存在建模挑战。在本研究中,增强了一种显式三维损伤模型并将其作为VUMAT在Abaqus中实现,以模拟复合材料构件的冲击响应。采用最大应力失效准则预测层内损伤起始,并应用能量驱动模式进行损伤演化。使用显式求解器可以增强解的收敛性,增强的VUMAT为静态失效分析提供了高效的解决方案。除了基于CDM的层内损伤描述外,还在每个层间应用了基于内聚单元的层间损伤表征。选择了[454/04/-454/904]sHexplyAS4/8552复合材料层压板的落锤冲击来展示其性能。将模拟结果与González等人进行的冲击试验和开源CompDam代码的结果进行了比较。为进一步展示所开发VUMAT在静态失效分析中的适用性,重新分析了开孔拉伸案例,并将预测结果与TechScout1试验数据和使用CompDam的模拟结果进行了比较。2、方法概述通过Abaqus的VUMAT实现了正交各向异性本构材料模型以描述损伤起始和演化。本研究采用最大应力失效准则来表征层内损伤。最大应力准则充分考虑了三维应力状态,允许分别模拟不同的纤维和基体失效模式。相应的基体和纤维在拉伸、压缩和剪切下的失效方程如下,其中1和2分别为纵向和横向,3为厚度方向:面内和面外剪切失效:纤维拉伸和压缩失效:基体拉伸和压缩失效:面外拉伸和压缩失效:其中,τij为计算的剪应力分量,σii为计算的正应力分量,S、X、Y、Z分别为相应的强度值。对于特征长度为Lc的单元,由于每种基体和纤维失效模式而耗散的能量由相应的纤维SERRGcf和基体SERRGcm决定:假设分层起始遵循二次失效准则:其中σn、σs、σt分别为正应力和两个剪应力,N、S、T为相应的界面强度。对于混合模式分层,引入等效牵引力Teq和位移δeqint,分层扩展由Benzeggagh和Kenane混合模式能量准则表示:一旦材料点发生基体或纤维损伤,损伤演化由给定断裂能量的裂纹带方法表征。在裂纹带方法中,假设分布式裂纹在有限元内某一宽度范围内被弥散,使得渐进开裂的效应在连续体方案中表现为宏观应变软化。为恢复网格客观性,引入特征长度,使得连续体单元在失效过程中耗散的总能量等于相同尺寸内聚单元定义的断裂韧性(或临界应变能释放率,SERR)。对于复合材料,裂纹带方法假设每种失效模式下的裂纹演化由相应模式的SERR控制。与基于二维应力状态的现有方法不同,本方法考虑了完整的三维应力状态来确定失效起始。对于特征长度为Lc的单元,由于每种基体和纤维失效模式而耗散的能量由相应的纤维SERRGcf和基体SERRGcm决定。与上述层内损伤预测不同,层间分层使用Abaqus现有的基于混合模式牵引-分离定律的内聚单元来表征。假设分层起始遵循二次失效准则。对于混合模式分层,引入等效牵引力Teq和位移δeqint,分层扩展由Benzeggagh和Kenane混合模式能量准则表示。本文开发的CDM模型的损伤预测与CompDam的预测进行了比较。CompDam是NASA开发的Abaqus/Explicit用户材料(VUMAT)子程序。CompDam的一个关键特征是其能够根据变形梯度分解方法准确表示基体裂纹运动学。3、模拟与结果HexplyAS4/8552层压板落锤冲击试验模拟为展示所开发显式CDM模块的能力和精度,模拟了HexplyAS4/8552层压复合材料板的落锤冲击。复合材料板的铺层为[454/04/-454/904]s,尺寸为150x100x5.8mm3。复合材料板的每一层均使用三维实体单元(C3D8R)显式建模,而各层之间的界面使用Abaqus内聚单元(COH3D8)表征。冲击体描述为直径16mm的刚体(R3D4)。使用四个带橡胶头的刚性夹具(R3D4)来约束试样与基座支撑。在板与冲击体、刚性夹具和固定基座之间强制执行接触。复合材料与钢冲击体之间的接触摩擦系数定为0.3。图1显示了[454/04/-454/904]s复合材料板及冲击体的有限元模型。HexplyAS4/8552的材料属性和界面属性分别列于表I和表II。使用Abaqus显式求解器进行了包括初始夹紧和冲击模拟的一步分析。将本文3DCDM的结果与González等人的试验数据进行了比较,并与开源代码CompDam进行了基准对比。图1有限元模型表1复合材料属性表2界面层属性图2展示了[454/04/-454/904]sAS4/85520层压板在38.6J冲击能量下的厚度方向剖面图及其损伤预测示例。从图2可以看出,模型预测了由面外应力分量引起的初始基体压缩和剪切失效,随后是由面内应力引起的额外基体拉伸和剪切失效。图3a和3b展示了我们的3DCDM与CompDam在相同冲击能量下的基体损伤预测比较。可以观察到,CompDam预测的厚度方向基体损伤比我们的3DCDM更多,但在与冲击体接触区域周围的损伤较少。图3c绘制了两种代码预测的冲击力-时间历程曲线与落锤冲击试验结果的比较。结果表明,我们的CDM模型正确预测了复合材料板的全局力-时间响应。尽管CompDam代码很好地捕捉了冲击事件早期阶段板的行为,但它高估了试验峰值载荷,而我们的CDM预测与测得的峰值载荷更为一致。开孔拉伸试样模拟为进一步验证所开发3DCDM模型在面内载荷情况下的性能,使用开孔拉伸(OHT)试样进行了验证。该试验旨在检验所开发VUMAT使用Abaqus显式求解器进行静态失效预测的能力。复合材料板的铺层为[60/0/-60]3s,由通用碳纤维增强环氧树脂带材制成。材料属性总结于表III。层压板尺寸为138x38.1x2.3mm3,一端夹紧并施加拉力,另一端固定。图4显示了OHT试样的有限元模型。板的每一层均使用实体单元(C3D8R)建模,而在每层界面引入内聚单元以模拟分层。应用我们开发的3DCDM和CompDam来预测[60/0/-60]3sOHT板中的损伤起始和扩展。图5和图6分别展示了我们的3DCDM和CompDam预测的OHT板各层基体损伤与扫描试验数据的比较。可以观察到,预测的基体损伤主要出现在±60°层,从缺口处开始。两种代码的损伤预测结果一致,并与扫描试验数据吻合良好。此外,图7绘制了3DCDM、CompDam预测的应力-应变曲线与试验数据的比较。结果表明,3DCDM预测的OHT失效应力略高于CompDam。我们的3DCDM预测OHT试样的峰值应力为478MPa,约为试验获得失效载荷的88%。预测与试验之间的差异可能归因于我们CDM模型中使用的横向拉伸强度(Yt=44.4MPa)值较低,低于材料的真实属性。不同层损伤对比应力应变曲线对比4、总结本文展示了一种实现为AbaqusVUMAT的显式三维渐进损伤模型,用于动态和静态载荷下的两个试验案例。动态模拟案例涉及[454/04/-454/904]sHexplyAS4/3501-6复合材料层压板的低速冲击,而静态案例涉及[60/0/-60]3s层压板的开孔拉伸试验。我们的3DCDM模型Abaqus显式求解器捕捉了复合材料层压板各层的基体拉伸和压缩损伤、剪切失效、纤维拉伸和压缩损伤以及界面分层。将其数值预测结果与开源CompDam代码的预测进行了基准对比,并与试验数据进行了验证。模拟与试验结果达成了良好的一致性,表明我们的显式3DCDM能够有效预测复合材料板在面外冲击和面内载荷情况下的多种失效模式。除了通过使用显式求解器增强数值稳定性外,当前的解决方案框架还展示了其在模拟冲击响应后预测压缩残余强度的可行性,无需用户干预。来源:ABAQUS仿真世界

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈