变幅杆的ANSYS分析设计

1、研究生课程论文课程名称现代设计方法题目变幅杆的ANSY鲂析设计学院机电工程学院专业机械工程班级姓 名 张佳欢向义河胡鹏飞吕建军熊海亮鲁中原娄晓刚指导教师李刚炎2014年 12月 15 日变幅杆的ANSYS分析设计摘要由于很多材料的脆硬特性，传统加工方法已不能满足对这些材料零件的精密加工要求，尤其是在航空航天的需求更加强烈，超声波振动技术是日本首先提出来，然后超声振动精密切削技术便逐渐成为各国的研究热点。在设计超声波振动系统中， 变幅杆的设计是核心部分，因此变幅杆是整个系统的核心部件。本文利用ANSYS有限元分析软件对变幅杆进行模态分析 及谐响应分析，最终根据分析结果对变幅杆进行修正。关键词：变

2、幅杆，超声波，有限元分析AbstractBecause a lot of hard brittle characteristics, traditional processing method has been can not meet the requireme nt of the precisi on mach ining of these materials parts, especially in aerospace n eeds more inten se, ultras onic vibrati on tech no logy was put forward by Japa n&#

3、39;s first, and the n ultras onic vibrati on precisi on cutt ing tech no logy has gradually become the focus in the coun tries. In the desig n of ultras onic vibrati on system, the desig n is the core part of the amplitude, so the amplitude is the core component of the whole system. This paper, by u

4、sing ANSYS finite element an alysis software for the amplitude modal an alysis and harm onic resp onse an alysis, fin ally to modify the amplitude accord ing to the results of the an alysis .Keywords : Amplitude， Ultrasonic wave， Finite Element Ansys19.研究内容简介碳纤维在高温高压环境中仍然能保持较高的强度，尺寸稳定性好，并且热膨胀系数小， 质量

5、轻，耐侵蚀。由于它本身的结构性能是由碳纤维的质地柔软而且粘性大的性能与强度较 高、硬度大的纤维组织的性能共同决定的，两者呈现出二相或多相结构。其机械加工性能比较差，是一种典型的难加工的材料。主要呈现特点如下：材料产生分层破坏，刀具磨损严重，寿命短，残余应力超声波振动切削实际上我们可以理解为超声波产生的振动，然后我们将这产生的振动施加到刀具上面，而且由于超声波是正弦波且具有周期性，所以刀具能够进行周期性的振动， 正是由于刀具周期性的振动，使得刀具与工件间能够产生有规律的分离作用，而这一作用则是超声波振动切削的最大优势。在对超声波振动切削系统进行设计，并用matlab对理论计算值进行验证后。本文主

6、要利用有限元软件对变幅杆进行动力学分析，优化其结构。本课题中超声波振动的频率为 40千赫兹。图1超声波振动车削系统1 超声波发生器2 换能器3 变幅杆4 刀具在超声波振动车削过程中，与传统车削最大的不同是，施加了一定振幅的超声波振动，虽然这一振动的振幅值很小， 但是超声波振动具有 20kHz以上的高频率特性， 所以刀具在单位时 间内的振动次数多，刀具产生的速度快。在这种激励作用下，工件产生微观的裂纹，材料发 生塑性变形而被去除，工件表面粗糙度低，加工精度好。1) 切削力：切削力仅仅是普通车削的10%£右，切削力大大降低，切削时间缩短短，工件与刀具摩擦下降，切削温度降低，减少了热损伤和

7、表面残余应力，以及热变形等。2) 表面粗糙度：超声车削加工不会产生积屑瘤，表面粗糙度降低，从而提高了加工 精度。3) 刀具寿命：由于是间歇切削，切削温度低，刀具发热减少，磨损降低，寿命得到明显提高。4）切屑：切屑颗粒小，切屑温度低，切削温度低，排屑比较顺畅。5）切削液：超声振动车削，当工件与刀具分离时，切削液容易进入切削区，冷却润 滑充分。超声波振动车削加工与普通车削加工中，两者刀具均存在着振动。不同的是普通切削加工中刀具的振动无规律性，无法控制刀具的运动轨迹，使得工件的加工质量下降和刀具的磨 损严重。而超声波振动切削的刀具的振动具有一定的振幅，且具有高频特性，增大剪切角， 可以减小刀具与工件

8、的摩擦因数，使得应力更集中，工件的刚性得到加强。二动力学分析原理相比于静力学分析而言，对于系统的动力学分析，它主要是确立结构的动力学特性， 譬如振型或者振动频率等。在本次课题的超声波振动切削系统中，对于变幅杆，我们主要是对其进行动力学分析。在本课题中，超声波振动系统主是实现刀具的纵向周期性振动，而实现这一目标的核心部件则是变幅杆，变幅杆的振动性能必须良好。我们仅仅通过理论计算得到的各项参数并不能完全满足振动性能的要求，只要当激振源也就是超声波发生器与变幅杆的固有频率相匹配 时，整个系统的振动系统才能达到要求。为了达到这一目的， 我们必须要对变幅杆进行相关的动力学分析，以此来优化变幅杆的结构尺寸

9、。在这次的分析中，我们主要使用的是动力学分析中的模态分析和谐响应分析。下面将详细介绍一下两种动力学分析。1）模态分析模态分析的主要作用是用来确定系统的振动特性参数即系统结构的固有频率和振动形式。模态分析的利用系统的物理结构方程求解，得到其特征值和特征向量。M ex m x: - 心（i）上式是无阻尼结构系统的振动方程：式中：M质量矩阵；K刚度矩阵；X加速度矩阵；x位移矩阵；若令：xUTsi n（ 7）（ 2）则有：X - 2 ' ?sinc't）（3）代入式1可得： KI-：2 IM DO(4)上式就是系统的振动特征方程，我们要求解特征方程的特征解，j及其对应的特征向量/，也就

10、是模态分析的实质。 其中打就是系统结构的固有频率， 而/也就是结构对应的 振动形式。在ANSYS件中主要有以下几种提取模态的方法：分块法、子空间法、 缩减法、动态提取法、不对称法、阻尼法、QR阻尼法等，使用何种模态提取方法主要取决于系统结构的大小和具体的场合。但是大多数模态分析主要采用缩减法完成。2)谐响应分析我们运用谐响应分析，主要是用来求得系统结构的动态性能特性。系统的结构必须为线性系统结构，目前暂时不考虑非线性系统。它是在外界施加线性系统某一固定频率的正弦载荷，禾U用谐响应分析，求得系统的稳态相应。在正弦载荷的作用下，系统结构以某一频率振动，其运动方程：M / U C 卜皿 I.K 卜U

11、 - F / sin vt(5)式中：C阻尼矩阵；F 简谐载荷幅值向量二激振力频率。位移方程为：xFAlsin 讹(6)在上式中：位移幅值向量；:位移响应滞后激励载荷相位角。谐响应分析在动力学分析中，占有比较重要的位置，它可以通过在线性结构系统上施加 正弦载荷来求解系统的动态性能响应。 在有限元分析软件中， 给我们提供了三种最常见的方 法，分别是：完全法、缩减法和模态叠加法。三种方法中，功能最强大的则是完全法，它是采用完整的系统矩阵进行计算结构的响应， 适用于结构复杂的系统进行谐响应分析。它也是系统默认的方法， 对计算机的硬件要求也是比较高的。但是该方法也有其优点， 操作简便，多载荷没有明显的

12、限制和要求，因此使用比 较广泛。但是也存在一定的缺陷，就是预应力选项是不可采用的。第二种是缩减法。缩减法也是相对于完全法来说的，它主要采取自由度和矩阵的缩减和压缩解决问题。通过求解主自由度的位移后，分析软件在将其扩展到整个系统结构上。相比于完全法，其运算速度较快而且占用系统的空间比完全法小了许多。第三种则是模态叠加法。它则是需要利用模态分析所得到的系统固有频率以及振型来求解。它则是三种方法中计算最快的方法。缺点是，需要进行模态分析后求解。操作相对复杂。三. 圆锥型变幅杆动力学分析 3.1圆锥型变幅杆模态分析3.1.1变幅杆前处理我们根据变幅杆的相关结构尺寸，进行三维建模。这次我们运用的三维建模

13、软件是Solidworks。建模成功后，我们将变幅杆的三维结构直接导入到有限元分析软件中。然后我们首先对设置单元类型为 SOLID45体单元类型；然后，定义材料属性：分别是材料的弹性模 量（E）、泊松比（）、密度等，相关参数如下表1所示；最后变幅杆三维模型进行自由映射网格单元网格划分的前处理操作，划分后的模型如图2所示。0.0070)0 (mm)35.0017.5052,50表1变幅杆的材料属性材料弹性模量E/（GPa）密度P( kg/m3)泊松比声速c(m/ s)45#21078800.285170图2变幅杆网格划分3.1.2施加约束及求解根据实际使用情况， 法兰盘处我们固定在刀架上面， 属

14、于完全约束。所以在法兰盘的地 方施加完全约束。然后定义分析类型为模态分析，采用方法为分块法，频率区间为30kHz到50kHz,设置提取10阶模态。最后求解。3.1.3后处理及查看结果后处理阶段，结构的固有频率、振型、应力及位移等已经写入结果文件。如图3所示是模态分析的结果列表，可以看出，在频率为30-50kHz内共得到7阶模态，如图4所示。从图4中我们可以看出，频率为41378Hz、41406Hz及44189Hz各阶分析中，变幅杆发生的是弯曲变形；频率为38929Hz，变幅杆发生扭转变形，这些振动均不是纵向变形。再看频率为34708Hz、34058Hz，变幅杆主体振动脱离法兰盘固定，脱离实际要

15、求，这也不是我 们设计中所需要的纵向振动；最终只有频率为40193Hz，变幅杆只发生纵向变形，而且是纵向振动，是我们设计所需要的纵向振动。鹿 MfwUIT删Tn L TM :帧Wit mm3DI4JIIOI II ir图3圆锥型变幅杆分析结果列表AN ；YS14.0"OU i 0 r tnTiflhrn JFiiiq吐iey Urtrnniln mi3 OS客 ujri”43DKIIOEH ITJ 1031-3 -J St tn uh(a)阶振型(b )二阶振型Jt MhJmJgfi*. T曲MfWfTrig Frqjim-qiri-b血mm3DI4JI Uli II IIIM

16、77;/Nmk llEJI111 lain W*3?&3 Ml «l ME 315?"(c)三阶振型Jt NnJflTi加 Cxtt>Wi 5 Tm 伽i FrwijmfHjUrfi miri(e)五阶振型Tzl.slT'rti-nitah? *Fhm Frt网册m冲e FsfitfJimcY W3I7 i腌 UdE nrm 24Unir»1i.i» 222Ed L- 1弘掘 M伯氐蛇 柑亂伯 in n55IJWIIflJflNwc 1-54.11 iis-at IIS-JB 弼的 TF 1-93 知蔗id(d )四阶振型业

17、87;34iiiCielCfmfltirnCiFiBEfiinncY陶Urrm-TiANSYS14.03DIVI IDE II IIIHUVNtan 1知打131 LT = 112K EjLJ JS JTT -W G-eJ戸邛9J|>J(f )六阶振型SQ14/0(g )七阶振型图4分析得到的各阶模态振型(圆锥型)通过以上对结果得到， 虽然其中有我们所需要的振型，但是其振型的固有频率与我们设计的理论值偏差较大，我们必须对变幅杆的结构尺寸进行修改完善，已达到我们理论要求值。 多次取值计算分析找到最佳值。分别取 L = 74,75,76,77,77,78mm和 x0=27,28,29,3O,

18、31mm进行多次模态分析，查看纵向振动的固有频率以及变形情况，从 分析的结果可以看出，当L=77mm x0=30mm寸，得到的分析结果和我们的理论设计值最近。而这里主要修改的结构尺寸主要是对变幅杆的共振长度L以及节点位移 x0进行修改，通过我们对L=77mm x0=30mn进行变幅杆的模态分析，搜索频率的区间还是30kHz到50kHz。分析方法和前面一致，得到七阶模态，如图4-4所示。Ta bu I a r 匚>曰菱日图5圆锥型变幅杆各阶模态频率如图6和7分别是频率为39784Hz的总变形图和轴向位移图。 从图中可以得出在法兰盘处是位移零点，也就是该点处于固定状态。TWP 叶加皿:如険3

19、但辰曲LWtrrmMabai Caiiikniita SfiAirnHI ilrhHt 14:37他碎吐9 I'M*fl MSI1 Ei?1胸12W!血如MilANS>S14.0J5oa伽心阳WriOgg FW 1询 胆巾irralcn FiiefjAiet bfi3DI4IIOHII3<W31 9M.iUES54TLi-I u.?nJ i9 抽。H U.1H25JM图6轴向形变图图7总变形图3.2圆锥型变幅杆谐响应分析3.2.1变幅杆的三维建模模型就采用尺寸修改后的模型，用solidworks 进行建模。3.2.2 加载及求解定义分析类型为谐响应分析，其次分析方法选择为完

20、全法，然后加载，在法兰盘处施加完全约束，在输入端施加一个大小为2um的位移载荷，进行求解。323后处理1）轴向位移曲线在通用后处理器 P0ST1中查看分析结果：在输出轴向位移曲线时，选择变幅杆中心线 上均匀分布的 若干个节点作为绘图特殊点，得到输入为2 u m时轴向位移曲线（如图 8）PI37图8轴向位移曲线从上图可以看出，在小端面输出位移为6.021um，从而得知放大系数大约为3左右，位移为零处大约在x=30mm处。2）幅频响应曲线在0-50kHz范围内，得到幅频响应曲线。 有曲线可以看出，频率在40kHz附近振幅最大， 到达共振，这个与设计值相符，符合变幅杆的设计要求。图9幅频响应曲线3.

21、3优化后的结果分析经过优化后，共振频率为39784kHz,与设计出入不大， 共振长度为L=77mm优化后的图 如9所示。图10优化后圆锥型变幅杆四. 阶梯型变幅杆动力学分析4.1阶梯型变幅杆模态分析4.1.1三维模型建立0.007D.00 （耐a）17 5Q5250图11阶梯型变幅杆网格划分与圆锥型变幅杆一样，首先，设置单元类型采用SOLID45体单元类型；其次，定义材料属性与圆锥型设置相同，在此不做详细叙述；最后将建立变幅杆三维模型导入ANSYS中并对其进行网格划分，划分后的模型如图4-10所示。4.1.2施加载荷以及求解模态分析时内容圆锥型变幅杆设置情况相同，在此不做过多叙述。4.1.3后

22、处理如图12和13所示，在频率区间内共得到七阶模态。 其中二阶振型是纵向变形， 其他都 是弯曲变形。虽然二阶振型纵向变形， 但是其共振频率比我们设计的值远远要小，所以同样要对变幅杆尺寸进行优化以达到我们理论设计的工作频率。图12阶梯型变幅杆频率A： MultiTVw TW 帥e 种 ijiifriOi irM 14 41ANSS讯EH 叨白m +i Jia zr 13.9734U LJU !>m- ；TAM 口册E Pt EF Hqumty卫卯曲 UdiE 加rri111.14 BTjaia+i aii ?r的i(b )二阶振型(a) 阶振型ElTW! n R F-MquimEY WQI

23、lMg WiE HWii 21UHir»i4 44iniUMK321113胡山诃?j51S4H0l± 刘弓国 now 1I.5D34U LJd !>nr ；tltZZMH73.113 酣 JOBS 生i隔 4-5XI? 1 対电仙 2711429DHI押PPA:Mcil*lT 询 WoF?ifl(jn* hwTnWOf<t*TTirtiji FMqumEYUriiE HWii ±H4>iir»i45(c)三阶振型(d)四阶振型* HbljHt油| 毗iFiiflan. TTW 1W 口 fftym e rtjuimEY iifriiM；

24、 uwr nwri2114'iirM 14-5.WBI7M4H国1何ANSS14.0ikHoliHT：hi r*rti-naHant iw Trrfeeopwirtrn F *甲imtY 41 IflEi 吐 UriiE M14.nirMi45n.34151JEI1S1J4U加2闻21 mWBSMjr 闵UEti T5JM5111M 049 刘拙T71炉 iajw(e)五阶振型(f)六阶振型JUrifliDFV<m.!non J tHM 口理细洛川乂閒 FlREflMIEV ill liH 吃 ijrr-TtTi Di4iiOBi?arEEE11kl.J|4MdK 163K 14

25、345 1ZZM IHiT MIN t! 4M 仍珂 事HIM Ml21(g)七阶振型图13阶梯型变幅杆各阶振型经过多次修改分析，最后得到L=60mmt匕较合适。4.1.4修改后变幅杆的模态分析当L =60mm x =30mm寸，对阶梯型变幅杆的模态分析，频率搜索区间还是20kHz到50kHz之间。得到如下结果：Tabular DataModeI" Frecquency Hn1:L3393922.3394733,380944442861554437366*4439077.47982图14修正后尺寸的模态分析列表4.2阶梯型变幅杆谐响应分析4.2.1三维建模采模态分析时进行优化后的变幅

26、杆三维模型。4.2.2施加载荷及求解加载过程与圆锥型变幅杆设置相同，在此不做过多叙述，最后进行求解。通过谐响应分析后得出在频率为 40kHz时轴向位移分布图和应力分布图，如图15所示。QnfClImflMrHnuiilfnTjipe-: Drcrlkinri DrtmiitanlX.Ms；* Fraqua<K¥：4D：OD HzPium空护d BVhR. nimGleft,I CwnJ niteMHHUM 1623ErjffwWr* SJrTjipe: EquwiienlStrewFraquaiK¥：4D：0D HzPruzpXQ比 d *MF?:.侦ooQian? unwEsw OOMfllldfiiQ0WKZ1I7 -QJ90HM DJdDigiiM -0001451： 申 Q(XN|0 -DWOJKSBn2D.MiD.K imn3D.5HHi 44iiar1D2M5.1249 ILII1：inRI Mu2D.MZQIVHJ