|
摘 要:介绍电力电子设备的抗震要求,及常用模态求解方法。重点运用ANSYS进行某框架刚度优化设计的过程,最后结合实验结果对有限元分析加以证明。
9 C4 X% J' `* B6 R7 ` 关键词:抗震;模态;ANSYS;有限元
% p, Z/ X" {+ W {; b
" l/ H {, P8 g) S 电子电力设备在使用过程中要经受恶劣的环境,在振动环境下,由于振动的疲劳效应及共振现象,可能出现电性能下降、零部件失效、疲劳损伤甚至破坏的现象;特别是户外用电子通讯设备在地震发生时,如能够抵抗地震冲击带来的破坏,则能够在抗震防灾中发挥巨大作用。因此对户外用电子电力设备的可靠性提出了严格要求,要求其具备一定抗振动与冲击的能力。在新产品的设计开发过程中通常用模态分析技术进行结构动态特性的估计以及优化设计。这一过程要得到系统的固有振动特征,然后在给定的使用或振动环境中进行评估,根据结果进行优化。 " p$ x# d, o. W" q9 J6 W
" |# u: L, F. L8 F 1 常见模态分析方法
# X* J; L: q' {4 X8 j- k- k( d: a/ q) k 1 B q8 t7 ?* i) u
1.1 理论求解法 4 s3 Y) \' s* ^3 c
理论求解是一种传统的求解途径。对于一些理想的简单弹性振动,可以用模态分析理论得到精确解。特别是可简化为单自由度的振动系统可以方便的求解。对于多自由度系统,特别是系统结构复杂,无法进行理想化时,其质量矩阵和刚度矩阵无法通过简单的理论计算来获得,也就无法使用理论法来求解。 7 b- Q+ D( K! W) e. n
1.2 试验法 ) j! q" m( H9 R; g( C( T
试验法是最直接的获得系统固有特征的方法。试验装置包括振动台、信号发生器、传感器、信号采集系统、频谱分析仪等。由于其需要实物原形,至少是模型结合相应的试验设备来进行,周期长,耗费较高的人力、物力等资源,很难在普通项目中运用。
0 ?8 _% M1 H) m1 O G6 u, S 1.3 有限元法 ! S2 C+ t; e+ T! ]0 V" t0 c; J% @# f
有限元法实质上是一种在力学模型上进行近似的数值计算方法。连续体结构体系本身内部不存在自然的连接关系,需要人为地在连续体内部和边界上划分节点,以分片(单元)连续的形式来逼近原来复杂的几何形状。线性有限元法是一种利用位能变分和分割近似原理求解线性弹性力学问题的数值方法。它首先把连续弹性体分割为在节点上相连的单元组合体,然后以节点位移为基本未知量,分别在各单元内选取位移函数,并按线性弹性力学的几何方程、本构方程和虚功方程或位能变分方程,建立并求解关于位移的线性代数方程组,把无限个自由度的问题化为有限自由度问题。随着有限元法的快速发展相继出现了许多有限元法的应用程序,如:ANSYS、I-DEAS、ADAMAS等。其中ANASYS是高校及工程界普遍使用的分析软件。
6 R8 s3 W* u) C; J5 d - @# K' h9 y9 d& K' ]7 \8 q: x
2 模态分析实例
* x: [) [7 h/ A j a * j4 V3 K% T. U4 N
某户外通信机柜(图1)的骨架由图2所示截面的杆件焊接而成,沿着杆件分布着图3所示的安装定位孔。机柜内部安装通信单元,要求能够在地震多发的环境中使用,在震灾发生时以及灾后要求设备正常工作。 0 b% v( d9 A- F1 d
7 \* x! m1 M, X0 a" n! R* }( ~) i 根据IEC防灾规范,地震的振动频率主要集中在0.1~35Hz,为避开共振产生的破坏,系统的固有振动频率应当尽可能高,特别是避开地震能量集中的低频范围0.1~6Hz。下面介绍如何使用ANSYS来进行固有频率及振型的求解,并根据分析结果进一步优化。 6 P: @$ g! t x/ Q+ _
2.1 问题分析 9 ?. p9 G* y( X0 M
此框架结构虽然是由杆焊接而成,由于沿着杆壁分布着安装孔,不能简化为梁单元。杆件是由薄钢板冲压后滚压成型,根据其截面性质(图2),可使用壳单元。ANSYS的建模功能相对比较弱,本例首先在Pro/E中完成建模,然后读入ANASYS中。
1 n- o/ m5 n2 x* l7 }( U 2.2 前处理
! G: H; J2 B3 r. o% ] 在Pro/E中长度单位为mm;为使单位统一,在ANSYS中用/UNITS命令设置系统参数。长度单位为mm,质量单位为g。在前处理器中设置单元类型为63号壳单元,设置相应的实常数并设置材料属性,泊松比=0.3,密度为 。用Modeling>Create>Glue命令把各梁结合在一起。接下来是网格划分,在Modeling>Mesh Attributes中设置单元号,材料号,单元实常数号;然后使用Mesh Tool进行网格划分。模态分析对于网格划分的质量要求并不高,可以采用先粗后细的方法来验证。即先用较少的网格数进行分析,再用两倍网格数分析。若两次比较结果差异在可接受范围内即可认为网格划分足够精确。 # U1 _) b5 c0 W
2.3 加载及求解 ; ~* b& u& \# ^1 P, j8 l4 y* i
对底部的四个顶点施加约束,选择ALL DOF,位移值设置为0。Solution>Analysis Type>New Analysis中设置模态分析;在Analysis Options中设置模态提取方法为Reduced法,并设定X为主自由度方向。运行Select Solve>Current LS,保存。接着依次设定Y、Z为主自由度,运行求解并保存。 A* Y: e5 a( N5 ^/ l- @
2.4 后处理
- d' s% E% W- J" y# F( o4 b1 | 通过General Postproc>Read Results命令依次读入各自方向的第一阶数据。通过Plot Results>Deformed Shape>Def+undeformed命令绘制变形图4~图6。) \' o- A5 D9 H
! L5 r6 d! c- c
X、Z、Y方向一阶固有频率分别是8.432,5.798,17.654;低于或是很接近地震能量集中的低频范围0.1~6Hz。由振型图可知X、Z方向振动模态主要是围绕焊接角点的摆动,而梁没有发生大的弯曲;Y方向则是在角点处发生大变形。观察图7所示的角点,它是沿着杆件的接缝直接焊而成,根据分析结果可知,角点局部结构刚度太低。改进方案应当是加强焊接角点的结构刚度。
" c7 c4 {7 G- H. H3 B) j 2.5 优化局部刚度 / _& Y X5 d9 V) _
采用图8所示的角件如图9所示,在每个角点插入角件,然后沿着接缝焊接牢固可起到此作用。 6 m. E3 v2 k9 {) ^7 P/ L7 [+ g1 e- j* o
; J4 ?- y" t( l6 P+ Q! S 在ANSYS中导入角件,修改模型,按照2.2~2.4步骤重新分析。得到的X、Z、Y方向一阶固有频率分别是:12.759,8.946,32.236。改进后的模型各方向上的一阶固有频率明显提高。
" c* @* E5 m. K2 u 3 试验分析
- m6 `7 M- b2 E ; k- `4 l. G' \5 `; E: Q) t) y; z# I* O
按照改进后的结构制作出框架原形,然后进行振动台测试。用白噪声进行固有频率扫描,获得的各方向一阶固有频率结果如图10~12所示。实验结果与有限元解接近。 2 T, C3 v- e8 o3 H6 r! a
9 d6 b& Z* M1 i% b4 m" ~# C6 |
4 结语 7 u# c! J. n' F4 n& l- Y7 f* g
- ~ `: _+ W; z9 s7 }# f 通过有限元与模态实验结果的对比,可以看出ANSYS能够为结构的刚度优化设计提供很大的帮助。并且它不像理论法,不需要烦琐的计算;不像模态实验,不需要昂贵的实验设备以及实验原型。从经济性、费用、可信度等方面来讲ANSYS不失为一种好的模态分析方法。
U4 J* H) P9 J+ B: ^ ^7 Y
0 Y3 D0 g1 W$ N, e l$ T8 ~) }) q 参考文献
) g+ @2 z" V( ^/ W$ F 〔1〕王勖成.有限元法〔M〕.北京:清华大学出版社,2003. ' G9 J( d, L/ o# R% D2 T
〔2〕师汉民.机械振动系统〔M〕.武汉:华中理工大学出版社,1990.
2 f: O9 _3 H% O 〔3〕李德葆.实验模态分析及其应用〔M〕.北京:科学出版社,2001.
1 Z& k% Q9 r; w/ i 〔4〕丰定国.抗震结构设计〔M〕.武汉:武汉工业大学出版社,2001. |
|