基于有限元的机床导轨结合面参数的特性分析_李景奎
现代制造工程2007年第5期制造工艺/工艺装备
基于有限元的机床导轨结合面参数的特性分析
李景奎,石宏,郭建烨
(沈阳航空工业学院,沈阳110034)
摘要:机床导轨的结合面是影响机床静、动态特性的关键因素。根据直线滚动导轨结合面的静、动刚度和阻尼等的实验结果,基于ANSY S 软件进行有限元分析,得到的结果接近于机床的实际情况,为机床的动态设计提供一种有效方法。关键词:机床;直线滚动导轨;结合面;有限元
中图分类号:TH 16 文献标识码:B 文章编号1671—3133(2007)05—0092—03
0 引言
机床是制造业的基础。机床的性能直接影响加
工产品的质量。而机床导轨的结合部分是影响机床动态特性的关键因素,因此世界各国都开展了对机床导轨结合面的分析研究。本文从单根直线滚动导轨入手,通过静刚度实验结果,得出导轨的法向静刚度、切向静刚度,并通过模态实验得出导轨的动刚度、阻尼;基于有限元分析软件对直线滚动导轨结合面的参数进行分析研究,得出整机的动态分析结果[1-4]
。其结果接近于机床的实际情况,从而使设计人员在图样设计阶段就能够预测整机特性。
1 导轨参数实验
直线滚动导轨的结合面参数包括结合面静刚度、动刚度和阻尼等几个方面。由于结合面的动态特性
具有非线性、多因素相关、依附于机械系统等特点,因此增加了研究的难度。目前对导轨参数的分析方法主要有实验测定和数值解析等方法。本文采用实验测定方法得到结合面的特性参数。1.1
静刚度实验
图1 直线滚动导轨刚
度实验结构布置
静刚度是滚动导轨副的一项重要使用性能。对滚动导轨副进行静刚度实验,可以验证滚动导轨副静刚度理论分析的正确程度,并为用户提供产品的性能依据和参考。静刚度实验[5]
结构
布局如图1所示,在夹具上表面均取4个测点进行静刚度实验,实验数据见表1。由表1中数据可以计算出,导轨法向静刚度k =[
∑(f
i
/y i )]/i =298N /μm ,
式中,f 为施加的力;y 为相对位移量;i 为实验次数。同理,得导轨切向静刚度为:235N /μm 。
表1静刚度实验数据
测点第一组实验千分表示值/μm 第二组实验千分表示值/μm 0N 980N 1960N 0N 980N 1960N 13.02.51.03.02.01.1236.030.526.036.030.525.5311.06.01.512
.06.52.24
1.0
-1.0
-3.0
1.0
-1.0
-3.0
1.2 动刚度实验
建立用模态参数表示的振动系统的运动方程,称
为模态分析,振动系统模态参数主要有:固有频率、振型、模态质量、模态刚度、模态阻尼。模态分析方法有时域法和频域法。本文采用的实验模态分析法属于
频域范畴[6]
,它是通过实验测定数据来确定模态参数的。实验原理如图2所示。进行法向刚度试验时,激振器垂直联接在滑块的中心,无松动;切向刚度试验时,激振器水平垂直于滑块的侧面,无松动。启动激振器,待系统稳定后,各种仪器开始记录。电平记录仪记录的数据结果如图3所示。
图2 动刚度测定实验原理
92
制造工艺/工艺装备现代制造工程2007年第5
期
图3 电平记录仪记录的数据曲线
从图3中可以看出,法向激振在690H z 时,振幅最大,系统发生共振。在这个频率下对试件、导轨等分别进行激振,得到振型如图4所示,可以发现共振是
由导轨与滑块之间的结合面引起的。同样,切向激振在534H z 时,系统发生共振
。
图4 690H z 时激振振型
1.3 导轨结合面动刚度、阻尼的计算
振动位移激励的频响函数(即位移导纳)为:
H x (ω)=1(k -m ω2)+jc ω=1︺k 1-ωω02
+2j ξω
ω0(1)………………………………………………式中:ω为角频率;k 为刚度;c 为阻尼;ω0为固有频率;m 为滑块质量,kg ;ξ为阻尼比。
速度导纳H v (ω)、加速度导纳H a (ω)、位移导纳H x (ω)之间的关系如下:
H v (ω)=j ωH a (ω
)(2)……………………………H x (ω)=-ω2
H a (ω)(3)………………………图5 阻尼比的估计
导轨滑块以及夹具的总重量为5kg ,导轨结合部法向固有频率为690H z ,相应的振幅为20dB ,0dB 对应标定值0.1mm /s 。由公式20lg (V /0.1)=20,求得:V =1mm /s 。
速度导纳:
H v (ωn )=
V /F =1/10=0.1mm /(s N )式中:V 为速度;F 为激振力,这里为给定值。根据图5
可估算系统的固有频率,在共振频率附近作一条距ω轴为0.707A (ωn )的水平线交幅频特性曲线于a 、b 两点,对应的频率为ω1、ω2
[7]
。
系统的阻尼比:
ξ=ω2-ω12ωn =740-670
2×690=0.05(4)
………………H x (ωn )=
1
2ξk
1-ξ2≈1
2ξ
k (5)…………………动刚度:k =12ξωn H v (ωn )=2.26E7N /m 阻尼:c =1063N s /m
导轨结合面切向固有频率为534H z ,相应的振幅
为27dB ,0dB 对应标定值0.1mm /s ,由公式20lg (V /0.1)=27,求得:V =2.2387mm /s 。
速度导纳:
H v (ωn )=V /F =2.2387/10=0.2239mm /(s N )同样可计算出系统的阻尼比:
ξ=ω2-ω12ωn
=600-5002×534=0.09
动刚度:k =1
2ξωn H v (ωn )=7.252E6N /m
阻尼:c =1084N s /m
2 导轨参数的有限元分析
2.1 模态分析
在考虑导轨参数进行有限元分析时,一般有两种分析方法:接触单元法和刚度阻尼弹簧法。接触单元法分为两种基本类型:刚体-柔体的接触和半柔体-柔
体的接触,它们都是一种非线性分析。ANSYS 支持三种接触方法:点-点接触单元、点-面接触单元以及面-面接触单元
[8]
。接触单元法需要设定接触刚度,接触
刚度由K =αP β
求得,其中K 为接触刚度,α、β分别为结合面特性参数,P 为接触面压强。刚度阻尼弹簧法是在两个子结构之间用一系列弹簧和阻尼器组成弹簧单元并连接起来,这一系列弹簧单元就构成了结合面的动力学模型,模型的动力学参数即为各弹簧单元的弹簧刚度和阻尼。这是一种比较成熟的结合面分析方法。通常,对考虑导轨参数的滑动导轨机床进行有限元分析时采用接触单元法,而对采用直线滚动导轨的机床在进行有限元分析时采用刚度阻尼弹簧法。本文采用的是刚度阻尼弹簧法,直线滚动导轨的动力学模型如图6所示。有限元分析模型包括床身、导轨、床鞍、小溜板部分,有限元模型如图7所示。其中模拟刚度、阻尼的单元选用COMB I N E14弹簧单元。在进
93
现代制造工程2007年第5期
制造工艺/工艺装备
行模态分析时,定义法向COMB I N E14单元的实常数为:刚度为:298N /μm ,阻尼为:1063N s /m ;定义切向COMB I N E14单元的实常数为:刚度为:235N /μm ,阻尼为:1084N s /m 。用ANSYS 程序求解后,得到机床的固有频率,一阶为131.132H z ,二阶为188.501H z ,三阶为201.217H z
。
图6
直线滚动导轨的动力学模型
图7 直线滚动导轨有限元模型
2.2 谐响应分析
在谐响应分析时,分两次进行。一次是X 向激振的谐响应分析,一次是Z 向激振的谐响应分析。在X 向谐响应分析时,选取分析频率在690H z 附近,法向
弹簧动刚度为:2.26E7N /m ,阻尼为:1063N s /m ,激振点在小溜板中心,拾振点在导轨上贴近滑块处;在Z 向谐响应分析时,选取分析频率在530H z 附近,切向弹簧动刚度为:7.252E6N /m ,阻尼为:1084N s /m ,激振点在小溜板左侧中心,拾振点在导轨上贴近滑块处。用ANSYS 程序求解后,得到谐响应分析结果如图8所示。
3 结语
利用ANSYS 对机床导轨结合面特性进行了有限元分析
[8]
,分析过程考虑了滚动导轨结合面参数,与
机床激振实验相对照,分析结果比较准确。其中结合面的静刚度可以通过静刚度实验获得,动刚度和阻尼可以从模态实验中识别。这种基于有限元的机床导
轨结合面参数的动态特性分析方法可以使设计人员图8 谐响应分析结果
在设计阶段便可获知机床的动态特性,使预测整机特性成为可能。
参考文献:
[1] 王海滨.机床结合面实验设计及有限元建模的研究
[D ].北京:北京工业大学,1999.
[2] 洪福昌.数控机床结合面特性参数的识别研究及其
应用软件的开发[D ].北京:北京工业大学,2000.[3] 韩会斌.机床结合面阻尼系数识别研究[D ].北京:
北京工业大学,2001.
[4] 陈永亮,李云梅,徐燕申,等.基于有限元的模块化振
动料斗动态分析[J ].现代制造工程,2005(3):79-81.
[5] 刘建素,李准.滚动直线导轨副静刚度试验[J ].制
造技术与机床,1997(7):12-13.
[6] 林国光.CD6245B 马鞍车床的试验模态分析[J ].组
合机床与自动化加工技术,2004(5):63-65.[7] 张耀满,王旭东,蔡光起,滕立波.高速机床有限元分
析及其动态性能试验[J ].组合机床与自动化加工技术,2004(12):15-17.
[8] 刘国庆,杨庆东.AN SYS 工程应用教程机械篇[M ].
北京:中国铁道出版社,2003.
作者通讯地址:辽宁沈阳市沈北新区37号航院新区(110034)E -m ail :m y -lj k @t om.co m 收稿日期:2006-12-25
94
有限元分析方法和材料断裂准则
一、有限元模拟方法 金属切削数值模拟常用到两种方法，欧拉方法和拉格朗日方法。欧拉方法适合在一个可以控制的体积内描述流体变形，这种方法的有限元网格描述的是空间域的，覆盖了可以控制的体积。在金属切削过程中，切屑形状的形成过程不是固定的，采用欧拉方法要不断的调整网格来修改边界条件，因此用欧拉方法进行动态的切削过程模拟比较困难。欧拉方法适用于切削过程的稳态分析（即“Euler方法的模拟是在切削达到稳定状态后进行的”[2]），仿真分析之前要通过实验的方法给定切屑的几何形状和剪切角[1]。 而拉格朗日方法是描述固体的方法，有限元网格由材料单元组成，这些网格依附在材料上并且准确的描述了分析物体的几何形状，它们随着加工过程的变化而变化。这种方法在描述材料的无约束流动时是很方便的，有限元网格精确的描述了材料的变形情况。实际金属切削加工仿真中广泛采用的拉格朗日方法，它可以模拟从初始切削一直到稳态的过程，能够预测切屑的形状和工件的残余应力等参数[2]。但是用这种方法预定义分离准则和切屑分离线来实现切屑和工件的分离，当物质发生大变形时常常使网格纠缠，轻则严重影响了单元近似精度，重则使计算中止或者引起严重的局部变形[1]。 为了克服欧拉描述和拉格朗日描述各自的缺点，Noh和Hirt在研究有限差分法时提出了ALE(Arbitrary Lagrange-Euler)描述，后来又被Hughes，liu和Belytschko等人引入到有限元中来。其基本思想是:计算网格不再固定，也不依附于流体质点，而是可以相对于坐标系做任意运动。由于这种描述既包含Lagrange的观点，可应用于带自由液面的流动，也包括了Euler观点，克服了纯Lagrange 方法常见的网格畸变不如意之处。自20世纪80年代中期以来，ALE描述己被广泛用来研究带自由液面的流体晃动问题、固体材料的大变形问题、流固祸合问题等等。金属的高速切削过程是一个大变形、高应变率的热力祸合过程，正适合采用ALE方法。 采用ALE方法进行高速切削仿真克服了拉格朗日方法和欧拉方法需要预先定义分离线、切屑和工件分离准则，假定切屑形状等缺点，避免了网格畸变以及网格再划分等问题，使切屑和工件保持良好的接触，使计算易于收敛[1][4]。 二、材料断裂准则 在金属切削成形有限元模拟中提出了多种切屑分离准则，这些准则可以分为两种类型：物理准则和几何准则。 优点： 几何分离准则需要预定义加工路径，在加工路径上判断刀尖与刀尖前单元节点的距离变化来判断分离与否。当两点的距离小于某个临界值时，刀尖前单元的节点被分成两个，其中一个节点沿前刀面向上移动形成切屑，另一个保留在加工表面上形成己加工表面[1][2]。。 物理分离准则是基于刀尖前单元节点的应力、应变及应变能等物理量定义分离条件，当单元中的该物理量的值超过给定材料的对应值时，单元节点就会分离[2]。（物理标准主要是基于制定的一些物理量的值是否达到临界值而进行判断的，主要有基于等效塑性应变准则、基于应变能密度准则、断裂应力准则等[5]）。 Carroll和Strenkowski使用了等效塑性应变作为物理分离准则的标准，在一些有限元软件中该标准的演化得到了应用，ABAQUS/Explicit中的剪切失效准则(shear failure)就是这样一种物理准则，它根据单元积分点处的等效塑性应变值是否到达预设值来判断材料是否失效[1]。 缺点：
有限元法
【第1章思考题】 1、何为有限元法？其基本思想是什么？ 1)“有限单元法”简称“有限元法”，是借助于电子计算机解决工程问题的近似方法。 2)“化整为零，集零为整”。也就是将一个原来连续的物体假想地分割成由有限个单元所组成的集合体，简称“离散化”。然后对每个单元进行力学特征分析，即建立单元节点力和节点位移之间的关系。最后，把所有单元的这种关系式集合起来，形成整个结构的力学特性关系，即得到一组以节点位移为未知量的代数方程组。处理后即可求解，求得结点的位移，进一步求出应变和应力 2、为什么说有限元法是近似的方法，体现在哪里？p3 用离散单元的组合体来逼近原始结构，体现了几何上的近似；而用近似函数逼近未知变量在单元内的真实解，体现了数学上的近似。 3、单元、节点的概念？ 网格划分中的每一个小部分称为单元。网格间相互联结点称为节点。 4、有限元法分析过程可归纳为几个步骤？p4 结构离散化、单元分析、整体分析 5、有限元方法分几种？本课程讲授的是哪一种？ 从选择基本未知量的角度来看，可分为3类： 1、位移法：以节点位移为基本未知量的求解方法称为位移法。本课程讲授的内容 2、力法：以节点力为基本未知量的求解方法称为力法; 3、混合法：一部分以节点位移，另一部分以节点力作为基本未知量的求解方法称为混合法。位移法 6、弹性力学的基本变量是什么？p8何为几何方程p11、物理方程p12及虚功方程？p14弹性矩阵的特点？ 弹性力学变量:外力、应力、应变和位移。描述弹性体应变分量与位移分量之间的方程称为几何方程；物理方程描述应力分量与应变分量之间的关系；弹性体上外力在虚位移发生过程中所做的虚功与储存在弹性体内的需应变能相等。弹性矩阵由材料的弹性模量和泊松比确定，与坐标位置无关。 7、何为平面应力问题和平面应变问题p17 平面应力问题：在结构上满足a几何条件：研究对象是等厚度薄板。b载荷条件：作用于薄板上的载荷平行于板平面且沿厚度方向均匀分布，而在两板面无外力作用。 平面应变问题：满足a几何条件：长柱体，即长度方向的尺寸远远大于横截面的尺寸，且横截面沿长度方向不变。b载荷条件：作用于长柱体结构上的载荷平行于横截面且沿纵向方向均匀分布，两端面不受力两条件的弹性力学问题 【第2章思考题】 1、何为结构的离散化？离散化的目的？何为有限元模型？ 结构的离散化:把连续的结构看成由有限个单元组成的集合体②目的：建立有限元计算模型 ③通常把由节点,单元及相应的节点载荷和节点约束构成的模型称为有限元模型 2、结构离散化时，划分单元数目的多少以及疏密分布，将直接影响到什么？确定单元数量的原则？通常如何设置节点？ ①单元的数量要根据计算精度的要求和计算机的容量来确定，因此在保证精度的前提，力求采用较少的单元。②节点的布置：a集中载荷的作用点b分布载荷强度的突变点 c分布载荷与自由边界的分界点d支承点e厚度不同或材料不同的区域等都应取为节点。 3、节点总码的编号原则？何为半带宽？半带宽与节点总码的编号有何关系？p21 ①节点编号时，应注意尽量使同一单元的相邻节点的号码差值尽可能地小些，以便缩小刚度矩阵的带宽，节约计算机存储。节点应顺短边编号为好②包括对角线在内的半个带状区域中每行具有的元素的个数，③半带宽B=（相关节点编号最大差值+1）*2
有限元分析复习内容汇总
1、有限元是近似求解一般连续场问题的数值方法 2、有限元法将连续的求解域离散为若干个子域,得到有限个单元,单元和单元之间用节点连接 3、直梁在外力的作用下,横截面的内力有剪力和弯矩两个. 4、平面刚架结构在外力的作用下,横截面上的内力有轴力、剪力、弯矩 . 5、进行直梁有限元分析,平面刚架单元上每个节点的节点位移为挠度和转角 6、平面刚架有限元分析，节点位移有轴向位移、横向位移、转角。 7、在弹性和小变形下，节点力和节点位移关系是线性关系。 8、弹性力学问题的方程个数有15个，未知量个数有15个。 9、弹性力学平面问题方程个数有8，未知数8个。 10、几何方程是研究应变和位移之间关系的方程 11、物理方程是描述应力和应变关系的方程 12、平衡方程反映了应力和体力之间关系的 13、把经过物体内任意一点各个截面上的应力状况叫做一点的应力状态 14、9形函数在单元上节点上的值,具有本点为_1_.它点为零的性质,并且在三角形单元的任一节点上,三个行函数之和为_1_ 15、形函数是_三角形_单元内部坐标的_线性_函数,他反映了单元的_位移_状态 16、在进行节点编号时,同一单元的相邻节点的号码差尽量小. 17、三角形单元的位移模式为_线性位移模式_- 18、矩形单元的位移模式为__双线性位移模式_ 19、在选择多项式位移模式的阶次时,要求_所选的位移模式应该与局部坐标系的方位无关的性质为几何_各向同性 20、单元刚度矩阵描述了_节点力_和_节点位移之间的关系 21、矩形单元边界上位移是连续变化的 1. 诉述有限元法的定义 答：有限元法是近似求解一般连续场问题的数值方法 2. 有限元法的基本思想是什么 答：首先，将表示结构的连续离散为若干个子域，单元之间通过其边界上的节点连接成组合体。其次，用每个单元内所假设的近似函数分片地表示求解域内待求的未知厂变量。 3. 有限元法的分类和基本步骤有哪些 答：分类：位移法、力法、混合法；步骤：结构的离散化，单元分析，单元集成，引入约束条件，求解线性方程组，得出节点位移。 4. 有限元法有哪些优缺点 答：优点：有限元法可以模拟各种几何形状复杂的结构，得出其近似解；通过计算机程序，可以广泛地应用于各种场合；可以从其他CAD软件中导入建好的模型；数学处理比较方便，对复杂形状的结构也能适用；有限元法和优化设计方法相结合，以便发挥各自的优点。 缺点：有限元计算，尤其是复杂问题的分析计算，所耗费的计算时间、内存和磁盘空间等计算资源是相当惊人的。对无限求解域问题没有较好的处理办法。尽管现有的有限元软件多数使用了网络自适应技术，但在具体应用时，采用什么类型的单元、多大的网络密度等都要完全依赖适用者的经验。 5. 梁单元和平面钢架结构单元的自由度由什么确定 答：由每个节点位移分量的总和确定 6. 简述单元刚度矩阵的性质和矩阵元素的物理意义 答：单元刚度矩阵是描述单元节点力和节点位移之间关系的矩阵 单元刚度矩阵中元素aml的物理意义为单元第L个节点位移分量等于1，其他节点位移分量
利用labview进行信号的时域分析
利用labview进行信号的时域分析 信号的时域分析主要是测量测试信号经滤波处理后的特征值，这些特征值以一个数值表示信号的某些时域特征，是对测试信号最简单直观的时域描述。将测试信号采集到计算机后，在测试VI 中进行信号特征值处理，并在测试VI 前面板上直观地表示出信号的特征值，可以给测试VI 的使用者提供一个了解测试信号变化的快速途径。信号的特征值分为幅值特征值、时间特征值和相位特征值。 用于信号时域分析的函数，VIs,Express VIs主要位于函数模板中的Signal Processing子模板中，其中多数对象位于Waveform Measurements子模板，如图所示 LabVIEW8.0中用于信号分析的Waveform Measurements子模板 基本平均值与均方差VI 基本平均值与均方差VI-------Basic Averaged DC—RMS.vi用于测量信号的平均以及均方差。计算方法是在信号上加窗，即将原有信号乘以一个窗函数，窗函数的类型可以选择矩形窗、Haning窗、以及Low side lob窗，然后计算加窗后信号的均值以及均方差值。 演示程序的前面板和后面板如下图所示 Basic Averaged DC—RMS演示程序的前面板
Basic Averaged DC—RMS演示程序的后面板 平均值与均方差值 平均值与均方差值VI------Averaged DC—RMS.vi同样也是用于计算信号的平均值与均方差值，只是Averaged DC—RMS.vi的输出是一个波形函数，这里我们可以看到加窗截断后，正弦信号的平均值和均方差随时间变化的波形。 编写程序演示Average DC----Averaged—RMS.vi的使用方法，程序的后面板和前面板如下图所示 Averaged DC—RMS演示程序的后面板
《有限元分析与应用课程标准》
《有限元分析及应用》课程标准 课程代码：汽车学分：3 建议课时数：64 英文名称： 适用专业：计算机辅助设计与分析 先修课程：《计算机辅助设计》 课程团队负责人及成员：陈良萍、刘宏强、王云、赵静、李蕾、黄艺、史俊玲、 毛新 1.课程定位和设计思路 1.1课程定位 本课程是为计算机辅助设计与分析专业本科生开设的一门专业核心课程，重点介绍有限元法的基本原理和方法、一些成熟的有限元软件功能和简单的分析步骤，同时结合工程实际，为他们进一步学习或实际应用及参加科研工作开辟道路。其任务是通过先修课程中所学知识的综合运用和新知识的获取，使学生初步掌握现代设计中的一种重要方法，开阔视野，提高能力，以适应科学技术发展的要求。 1.2设计思路 在教学中，首先通过力学中的矩阵位移法思想的对比教学，引出连续介质力学有限单元法的学习重点在于单元的插值函数如何构造。这因为，虽说矩阵位移法是对杆系结构而言的，但其结构的离散化和组建整体刚度方程的思想完全可以借鉴到连续介质力学，它们的不同点只是在单元刚度矩阵的建立；而不同单元类型的单元刚度矩阵的建立，又取决于对应单元插值函数的构造。这样处理，不但使学生抓住了本课程的教学重点，而且对有限单元法的整体思想有了宏观上掌握；起到主动学习而非被动接受的作用。在单元构造的教学中，理论学习的重点在于常规单元的介绍；通过常规单元介绍插值函数的完备性与收敛性等。接之，介绍高次单元、等参单元等教学内容。在理论教学中，强调数学论证的严谨性和工程应用的适应性。
结合工程实例教学，拓宽学生数值分析方面的应用能力在课内对不同的单元类 型进行介绍时，及时抓住不同单元在应用中的对比教学与其适用性，并结合工程实例介绍单元类型的合理选取和单元网格的合理划分等。为学生在实际问题的数值分析中如何选定单元和剖分单元奠定了一定的基础和经验。 2.工作任务和课程目标 2.1工作任务 由于采用有限单元法的分析计算软件大多已商业化，而熟悉应用这些中的常规软件也应是本门课程的主要教学内容。在课内学生学会使用软件建立分析模型的基本步骤，其中包括分析模型抽象、几何模型绘制、单元网格划分、材料定义、边界条件定义、方程求解方法等。因课内教学时数的不足，学生应利用课余时间学习，以提高对实际问题的数值分析能力。 2.2课程目标 从教学思想和方法上对原课程进行改革，使学生从较高层次上理解有限元方法的实质，掌握有限元分析的工具，并具备初步处理工程问题的能力;使该课程成为具有较宽口径和较大覆盖面的、面向计算机辅助设计方面的专业基础课；注意课程体 系的整体优化，强调课程的深度、广度与应用。 3.教学方针落实情况
有限元原理与步骤
2.1.1 有限元法基本原理（Basic Theory of FEM） 有限元法的基本思想是离散的概念，它是指假设把弹性连续体分割成数目有限的单元，并认为相邻单元之间仅在节点处相连。根据物体的几何形状特征、载荷特征、边界约束特征等，选择合适的单元类型。这样组成有限的单元集合体并引进等效节点力及节点约束条件，由于节点数目有限，就成为具有有限自由度的有限元计算模型，它替代了原来具有无限多自由度的连续体[24][25]。 有限元法从选择基本未知量的角度来看，可分为三类：位移法、力法和混合法。以节点位移为基本未知量的求解方法称为位移法；以节点力为基本未知量的求解方法称为力法；一部分以节点位移，另一部分以节点力作为基本未知量的求解方法称为混合法。由于位移法通用性强，计算机程序处理简单、方便，成为应用最广泛的一种方法[26]。 有限元法的求解过程简单、方法成熟、计算工作量大，特别适合于计算机计算。再加上它有成熟的大型软件系统支持，避免了人工在连续体上求分析解的数学困难，使其成为一种非常受欢迎的、应用极广泛的数值计算方法[27]。 2.1.2 有限元法基本步骤（Basic Process of FEM） 有限元法求解各种问题一般遵循以下的分析过程和步骤[28][29]： 1. 结构的离散化 结构的离散化是进行有限元法分析的第一步，它是有限元法计算的基础。将结构近似为具有不同有限大小和形状且彼此相连的有限个单元组成的计算模型，习惯上称为有限元网格划分。离散后单元与单元之间利用单元的节点相互连接起来，而单元节点的设置、性质、数目等应视问题的性质、描述变形形态的需要和计算精度而定。所以有限元法分析的结构已不是原有的物体或结构物，而是同种材料的由众多单元以一定方式连接成的离散物体。这样，用有限元分析计算所获得的结果是近似的。显然，单元越小（网格越密）则离散域的近似程度越好，计算结果也越精确，但计算量将增大，因此结构的离散化是有限元法的核心技术之一。有限元离散过程中又一重要环节是单元类型的选择，这应根据被分析结构的几何形状特点、载荷、约束等因素全面考虑。 2. 位移模式的选择 位移模式是表示单元内任意点的位移随位置变化的函数，位移模式的选择是有限元特性分析的第一步。由于多项式的数学运算比较简单、易于处理，所以通常是选用多项式作为位移函数。选择合适的位移函数是有限元分析的关键，它将决定有限元解的性质与近似程度。位移函数的选择一般遵循以下原则（有限元解的收敛条件）：
导轨的设计与选择
一、导轨的设计与选择。 1、对导轨的要求 1）导轨精度高 导轨精度是指机床的运动部件沿导轨移动时的直线和它与有关基面之间的相互位置的准确性。无论在空载或切削工件时导轨都应有足够的导轨精度，这是对导轨的基本要求。 2）耐磨性能好 导轨的耐磨性是指导轨在长期使用过程中保持一定导向精度的能力。因导轨在工作过程中难免磨损，所以应力求减少磨损量，并在磨损后能自动补偿或便于调整。 3）足够的刚度 导轨受力变形会影响部件之间的导向精度和相对位置，因此要求轨道应有足够的刚度。 4）低速运动平稳性 要使导轨的摩擦阻力小，运动轻便，低速运动时无爬行现象。5）结构简单、工艺性好 导轨的制造和维修要方便，在使用时便于调整和维护。 2、对导轨的技术要求 1）导轨的精度要求 滑动导轨，不管是V-平型还是平-平型，导轨面的平面度通常取0.01~0.015mm,长度方面的直线度通常取0.005~0.01mm；侧导向面的直线度取0.01~0.015mm,侧导向面之间的平行度取
0.01~0.015mm,侧导向面对导轨地面的垂直度取0.005~0.01mm。2）导轨的热处理 数控机床的开动率普遍都很高，这就要求导轨具有较高的耐磨性，以提高其精度保持性。为此，导轨大多需要淬火处理。导轨淬火的方式有中频淬火、超音频淬火、火焰淬火等，其中用的较多的是前两种方式。 二、导轨的种类和特点 导轨按运动轨迹可分为直线运动导轨和圆运动导轨；按工作性质可分为主运动导轨、进给运动导轨和调整导轨；按接触面的摩擦性质可分为滑动导轨、滚动导轨和静压导轨等三大类。 1）滑动导轨：是一种做滑动摩擦的普通导轨。滑动导轨的优点是结构简单，使用维护方便，缺点是未形成完全液体摩擦时低速易爬行，磨损大，寿命短，运动精度不稳定。滑动导轨一般用于普通机床和冶金设备上。 2）滚动导轨的特点是：摩擦阻力小，运动轻便灵活；磨损小，能长期保持精度；动、静摩擦系数差别小，低速时不易出现"爬行"现象，故运动均匀平稳。缺点是：导轨面和滚动体是点接触或线接触，抗振性差，接触应力大，故对导轨的表面硬度要求高；对导轨的形状精度和滚动体的尺寸精度要求高。因此，滚动导轨在要求微量移动和精确定位的设备上，获得日益广泛的运用。 3）静压导轨是利用液压力让导轨和滑块之间形成油膜，使
时域分析
机械振动故障诊断中时域参数指标的分析 一、滚动轴承的失效形式 1.疲劳剥落 在滚动轴承的滚动或滚动体表面，由于承受交变负荷的作用是接触面表层金属呈片状玻剥落，并逐步扩大而形成凹坑。如继续运转，则将形成面积剥落区域。由于安装不当或轴承座孔与轴的中心线倾斜等原因将使轴承中局部区域承受较大负荷而出现早期疲劳破坏。 2.磨损 当滚动轴承密封不好，使灰尘或微粒物质进入轴承，或是润滑不良，将引起接触表面较严重的擦伤或磨损，并使轴承的振动和噪声增大。 3.断裂和裂纹 材料缺陷和热处理不当，配合过硬两太大，组合设计不当，如支撑面有沟槽而引起应力集中等，将形成套圈裂纹和断裂。 4.压痕 外接硬颗粒物质进入轴承中，并压在滚动体与滚道之间，可是滚动表面形成压痕。此外，过大的冲击负荷也可以使接触表面产生局部塑性变形而形成凹坑。当轴承静止时，即使负荷很小，由于周围环境的振动也将在滚道上形成均匀分布的凹坑。 5.腐蚀 电机或者机械漏电或者有部分静电时产生电流，一般轴承都是需要使用，在轴承内部可以在轴承的内圈、外圈、滚动体之间产生油膜（很薄左右），电流可以击穿轴承内部的（油膜），造成轴承内圈、外圈、滚动体之间的直接接触、在接触的表面会产生电击，对轴承的沟道造成损伤，从而引起轴承早期失效。 6.胶合 指滚道和滚动体表面由于受热而局部融合在一起的现象。常发生在润滑不良、告诉、重在、高温、启动加速度过大等情况下。由于摩擦发热，轴承零件可以在极短时间内达到很高的温度，导致表面灼伤或某处表面上的金属粘附到另一表面上。 二、时域参数主要参数指标 峰值、均值、方差、歪度、峭度、均方根值，波形指标、脉冲指标、峭度指标、歪度指标和裕度指标。其中前一类是有量纲指标，后一类是无量纲指标。 1.峰值 在某个时间段内幅值的最大值。由于它是一个时不稳参数，不同的时刻变动很大，因此常用来检测冲击振动。 2.均值 指在一组数据中所有数据之和再除以这组数据的个数，这里指所有幅值的均值。反映了数据趋势的大小。信号的均值反映信号中的静态部分，一般对诊断不起作用，但对计算其它参数有很大影响，所以，一般在计算时应先从数据中去除均值，剩下对诊断有用的动态部分。 计算表达式：11N i i x x N ==∑。 3.均方根值 也称有效值，在电路中定义为一确定的交流电相当于多大数值的交流电在相同时间内所做的功一样。它用来反应信号的能量大小，特别适用于具有随机振动的性质的轴承测量。在滚动轴承的故障诊断中，均方根值可以用来反应各个滚动体在滚道上运动时，由于制造精度差以及工作表面点蚀所产生的不规则振动状况。制造精度愈低或轴承磨损程度愈大，则均方根值值愈高。对于正常轴承以及表面发生点蚀的轴承均方根值很稳定，不受偶然因素的干扰；但对于表面剥落或局部损伤产生的冲击脉冲振动波形，脉冲幅值的大小均方根值是反映不出来的。
有限元边界条件和载荷
X边界条件和载荷 10.1边界条件 施加的力和/或者约束叫做边界条件。在HyperMesh中，边界条件存放在叫做load collectors的载荷集中。Load collectors可以通过在模型浏览器中点击右键来创建(Create > Load Collector)。 经常（尤其是刚开始）需要一个load collector来存放约束（也叫做spc-单点约束），另外一个用来存放力或者压力。记住，你可以把任何约束（比如节点约束自由度1和自由度123）放在一个load collector中。这个规则同样适用于力和压力，它们可以放在同一个load collector中而不管方向和大小。 下面是将力施加到结构的一些基本规则。 1.集中载荷（作用在一个点或节点上） 将力施加到单个节点上往往会出现不如人意的结果，特别是在查看此区域的应力时。通常集中载荷（比如施加到节点的点力）容易产生高的应力梯度。即使高应力是正确的（比如力施加在无限小的区域），你应该检查下这种载荷是不是合乎常理？换句话说，模型中的载荷代表了哪种真实加载的情形？ 因此，力常常使用分布载荷施加，也就是说线载荷，面载荷更贴近于真实情况。 2.在线或边上的力 上图中，平板受到10N的力。力被平均分配到边的11个节点上。注意角上的力只作用在半个单元的边上。
上图是位移的云图。注意位于板的角上的红色“热点”。局部最大位移是由边界效应引起的（例如角上的力只作用在半个单元的边上），我们应该在板的边线上添加均匀载荷。 上述例子中，平板依然承受10N的力。但这次角上节点的受力减少为其他节点受力的一半大小。 上图显示了由plate_distributed.hm文件计算得到的平板位移的云图分布。位移分布更加均匀。 3.牵引力（或斜压力） 牵引力是作用在一块区域上任意方向而不仅仅是垂直于此区域的力。垂直于此区域的力称为压力。
有限元复习重点
●有限元起源于20世纪50年代中期航空工程中飞机结构的矩阵分析。 ●有限元基本思想：在力学模型上将一个原来连续的物体离散成为有限个具有一定大小的单元，这些单元仅在有限个节点上相连接，并在节点上引进等效力以代替实际作用于单元上的外力。对于每个单元,根据分块近似的思想，选择一种简单的函数来表示单元内位移的分布规律，并按弹性理论中的能量原理(或用变分原理)建立单元节点力和节点位移之间的关系。最后，把所有单元的这种关系式集合起来，就得到一组以节点位移为未知量的代数方程组，解这些方程组就可以求出物体上有限个离散节点上的位移。 “一分一合”，化整为零，集零为整，把复杂的结构看成由有限个单元组成的整体。 ●单元、节点、边界：采用8节点四边形等参数单元把受力体划分成网格，这些网格称为单元；网格间互相连接的点称为节点；网格与网格的交界线称为边界。节点数和单元数目是有限的。 ●有限元法的优点：（1）理论基础简明，物理概念清晰，且可在不同的水平上建立起对该法的理解。(2) 具有灵活性和适用性，应用范围极为广泛。(3) 该法在具体推导运算中，广泛采用了矩阵方法，便于实现程序设计的自动化。 ●有限单元法分为三类:位移法（以节点位移为基本未知量）、力法（以节点力为基本未知量）和混合法（一部分以节点位移，另一部分以节点力作为基本未知量）。 ●有限元法分析计算的基本步骤可归纳如以下五点。1.结构的离散化（将某个机械结构划分为由各种单元组成的计算模型）在平面问题用三角形、矩形或任意四边形单元。在空间问题用四面体、长方体或任意六面体单元2.单元分析①选择位移模式（位移模式是表示单元内任意点的位移随位置变化的函数式，由于所采用的函数是一种近似的试函数，一般不能精确地反映单元中真实的位移分布）位移模式或位移函数：i n i i a y φ∑=②建立单元刚度方程e e e F k =δ，e 为单元编号；e δ为单元的节点位移向量；e F 为单元的节 点力向量 ；e k 为单元刚度矩阵.③计算等效节点力：用等效的节点力来代替所有作用在单元上的力。3.整体分析：整体的有限元方程F K =δ。K 为整体结构的刚度矩阵；δ为整体节点位移向量；F 为整体载荷向量。4.求解方程，得出节点位移5.由节点位移计算单元的应变与应力 ●有限元中得一个基本近似性是几何近似性 ●有限元中的变量：应力、应变、变形。基本方程有：平衡方程、物理方程、几何方程。边界条件：力边界、位移边界。 ●弹性力学的任务是分析弹性体在受外力作用并处于平衡状态下产生的应 力、应变和位移状态及其相互关系等。 ●外力：体力（分布在物体体积内的力---重力、惯性力、电磁力）、面力（分布在物体表面上的力---流体压力、接触力、风力） ●应力：物体受外力的作用，或由于温度有所改变，其内部将发生内力。
传统木结构的整体有限元分析
传统木结构的整体有限元分析 1.引言 中国古建筑是中华文明的重要组成部分，是中华民族乃至世界建筑艺术的瑰宝，具有极高的文物、历史和艺术价值。而其中的木结构古建筑，不仅蕴含了丰富的历史文化信息，由于其建筑材料和建筑方式的独特性，更有其独特而优良的力学性质。 对这些古建筑的动力特性的研究，从七十年代就已经开始了，但是由于技术的限制，这些研究还远远不够。随着社会的进步，人们也开始对古建筑的维护投入了更多的关注。因此对古建筑的研究也要求进一步的深入。 本文根据2007年1月18日木结构足尺模型振动台实验结果,采用有限元计算软件对木结构动力特性进行计算模拟,并将实验数据与计算值进行对比，希望更深刻地了解木结构建筑的抗震性能和结构耗能减震的基本原理，这对木结构建筑遗产的保护修缮具有重要的意义。 2.木结构整体有限元分析方法 早在1994年Kasal[1]等就利用大型商业有限元软件ANSYS对一层木框架房屋进线性的静力分析。在此模型中，剪力墙被简华成由刚性杆和斜向弹簧组成的桁架模型线性由斜向弹簧的单元特性来实现，而屋面和楼板被简化为超级单元。 2001年，由Slovenia的研究小组提出的Slovenia模型[2][3]将木结构房屋的整体分三个阶段：钉连接模型一墙体模型一木结构房屋整体模型。其研究思路为：先根据D分析剪力墙所得的滞回曲线，将每片墙简化成一个等效支撑框架。定义斜撑单元的参模型的滞回曲线拟合而得到，并采用CANNY-E（采用Newmark 算法）程序对整体行非线性动力时程分析。 3.木结构的整体有限元分析 3.1 足尺寸实验模型概况 本文以日本防灾科学技术研究所兵库抗震工程研究中心进行的足尺寸木结构的振动台实验为原型进行有限元分析。该振动台实验主要研究带墙体覆面板结构自振以及在不同地震波程度下的动力特性。模型标准层结构平面布置层高为2.93m，柱横向间距和纵向间距均为1.92m，采用以杉木为原材料的木框架结构。柱截面和基础梁截面均为120mm×120mm，屋面外框梁截面120mm×270mm，次梁截面为120mm×210mm，其梁和柱均为榫卯连接，墙面板为干式土壁覆面板。 3.2 有限元计算模型 本工作希望从数值方法出发，用简单有效的方法，建立木结构的有限元计算模型，对其动力特性进行计算模拟，并结合实验数据评判模型。 建立的有限元计算模型主要包括以下几个方面： （1）基础模拟。地震波在地表传播时，地基是一个变形体，地震发生时结构基础处各点的运动是不同的。但是，对于一般建筑物，其长度远小于地震波的波长(它和场地介质的情况有关)，因此通常情况下将建筑物的地基近似看作刚性盘体是合理的[8]。因此在本次实验中，基础梁是固定在振动台上，计算模型中假定基础为刚性连接。 （2）木框架模拟。实验中的木结构框架可视为一种梁柱结构体系。梁柱之间上下叉接，左右卡连，如图3所示是实验中梁柱榫卯连接。榫卯连接是介于刚接与铰接之间的半刚性连接，在进行有限元分析时，通常的方法是用空间二节点虚拟弹簧单元来模拟这种半刚性连接性质。在同一空间位置的梁柱各端部节点与相应梁柱构件各自对应，并选择合适的自由度赋予弹簧刚度参数，形成半刚性连接[5]。因此，在计算模型中，柱一柱、梁一梁和梁一柱之间用弹簧单元来实现它们之间半刚性的连接。 考虑到木构架材质主要发挥其顺纹力学性质，可以将材料近似看作各向同性。参考文献[4]本文采用的木构架材料弹性模量15.5×109Pa,密度为3766kg/m3,泊松比0.25。 （3）屋面板单元。实验模型中屋面刚度很大，可以认为是刚性的，因此用Shell63单元固接在屋面梁上模拟。屋面上的配重在剪力有限元模拟过程中，利用质量单元Mass21模拟，将屋盖配重按面积等效原则
语音信号的时域特征分析
中北大学 课程设计说明书 学生姓名：蒋宝哲学号： 24 学生姓名：瓮泽勇学号： 42 学生姓名：侯战祎学号： 47 学院：信息商务学院 专业：电子信息工程 题目：信息处理实践:语音信号的时域特征分析指导教师：徐美芳职称: 讲师 2013 年 6 月 28 日
中北大学 课程设计任务书 2012-2013 学年第二学期 学院：信息商务学院 专业：电子信息工程 学生姓名：蒋宝哲学号： 24 学生姓名：瓮泽勇学号： 42 学生姓名：侯战祎学号： 47 课程设计题目：信息处理实践:语音信号的时域特征分析起迄日期： 2013年6 月7日～2013年6月 28 日 课程设计地点：学院楼201实验室、510实验室、608实验室指导教师：徐美芳 系主任：王浩全 下达任务书日期: 2013 年 6 月 7 日
语音信号的采集与分析 摘要 语音信号的采集与分析技术是一门涉及面很广的交叉科学，它的应用和发展与语音学、声音测量学、电子测量技术以及数字信号处理等学科紧密联系。其中语音采集和分析仪器的小型化、智能化、数字化以及多功能化的发展越来越快，分析速度较以往也有了大幅度的高。本文简要介绍了语音信号采集与分析的发展史以及语音信号的特征、采集与分析方法，并通过PC机录制自己的一段声音，运用Matlab进行仿真分析，最后加入噪声进行滤波处理，比较滤波前后的变化。 关键词：语音信号，采集与分析， Matlab 0 引言 通过语音传递倍息是人类最重要、最有效、最常用和最方便的交换信息的形式。语言是人类持有的功能．声音是人类常用的工具，是相互传递信息的最主要的手段。因此，语音信号是人们构成思想疏通和感情交流的最主要的途径。并且，由于语言和语音与人的智力活动密切相关，与社会文化和进步紧密相连，所以它具有最大的信息容量和最高的智能水平。现在，人类已开始进入了信息化时代，用现代手段研究语音信号，使人们能更加有效地产生、传输、存储、获取和应用语音信息，这对于促进社会的发展具有十分重要的意义。 让计算机能听懂人类的语言，是人类自计算机诞生以来梦寐以求的想法。随着计算机越来越向便携化方向发展，随着计算环境的日趋复杂化，人们越来越迫切要求摆脱键盘的束缚而代之以语音输人这样便于使用的、自然的、人性化的输人方式。作为高科鼓应用领域的研究热点，语音信号采集与分析从理论的研究到产品的开发已经走过了几十个春秋并且取得了长足的进步。它正在直接与办公、交通、金融、公安、商业、旅游等行业的语音咨询与管理．工业生产部门的语声控制，电话、电信系统的自动拨号、辅助控制与查询以及医疗卫生和福利事业的生活支援系统等各种实际应用领域相接轨，并且有望成为下一代操作系统和应用程序的用户界面。可见，语音信号采集与分析的研究将是一项极具市场价值和挑战性的工作。我们今天进行这一领域的研究与开拓就是要让语音信号处理技术走人人们的日常生活当中，并不断朝更高目标而努力。 语音信号采集与分析之所以能够那样长期地、深深地吸引广大科学工作者去不断地对其进行研究和探讨，除了它的实用性之外，另一个重要原因是，它始终与当时信息科学中最活跃的前沿学科保持密切的联系．并且一起发展。语音信号采集与分析是以语音语言学和数字
对有限元方法的认识
我对有限元方法的认识 1有限元法概念 有限元方法（The Finite Element Method, FEM）是计算机问世以后迅速发展起来的一种分析方法。每一种自然现象的背后都有相应的物理规律，对物理规律的描述可以借助相关的定理或定律表现为各种形式的方程（代数、微分、或积分）。这些方程通常称为控制方程（Governing equation）。 针对实际的工程问题推导这些方程并不十分困难，然而，要获得问题的解析的数学解却很困难。人们多采用数值方法给出近似的满足工程精度要求的解答。 有限元方法就是一种应用十分广泛的数值分析方法。 有限元方法是处理连续介质问题的一种普遍方法，离散化是有限元方法的基础。 这种思想自古有之：古代人们在计算圆的周长或面积时就采用了离散化的逼近方法：即采用内接多边形和外切多边形从两个不同的方向近似描述圆的周长或面积，当多边形的边数逐步增加时近似值将从这两个方向逼近真解。 近年来随着计算机技术的普及和计算速度的不断提高，有限元分析在工程设计和分析中得到了越来越广泛的重视，已经成为解决复杂的工程分析计算问题的有效途径，现在从汽车到航天飞机几乎所有的设计制造都已离不开有限元分析计算，其在机械制造、材料加工、航空航天、汽车、土木建筑、电子电器、国防军工、船舶、铁道、石化、能源、科学研究等各个领域的广泛使用已使设计水平发生了质的飞跃。 国际上早在 60 年代初就开始投入大量的人力和物力开发有限元分析程序。“有限单元”是由Clough R W于1960年首次提出的。但真正的有限元分析软件是诞生于 70 年代初期，随着计算机运算速度的提高，内、外存容量的扩大和图形设备的发展，以及软件技术的进步，发展成为有限元分析与设计软件，但初期其前后处理的能力还是比较弱的，特别是后处理能力更弱。
有限元分析软件及应用
3.5 ANSYS软件加载、求解、后处理技术 3.5.1 ANSYS 3.5.1 ANSYS 荷载概述荷载概述 在这一节中将讨论: 有限元分析软件及应用 8 有限元分析软件及应用 8 A. 载荷分类 3.5 ANSYS 软件加载、求解、后处理技术 3.5 ANSYS 软件加载、求解、后处理技术 B. 加载 C. 节点坐标系 D. 校验载荷 孙瑛 孙瑛 E. 删除载荷 哈哈尔尔滨滨工工业业大学空大学空间结间结构研构研究中心究中心 2010秋 2010秋 SSRC SSRC 1/ 76 S Space pace S Stru truc ctu ture re R Res esear earc ch h C Center enter, H , HI IT, T, CH CHIN INA A
理技术 A. 载荷分类 B. 加载 A. 载荷分类 B. 加载 ANSYS中的载荷可分为: 可在实体模型或 FEA 模型节点和单元上加载自由度DOF - 定义节点的自由度( DOF )值结构分析_ 沿单元边界均布的压力 沿线均布的压力 位移集中载荷 - 点载荷结构分析_力面载荷 - 作用在表面的分布载荷结构分析_压力 在关键点处 在节点处约 约束体积载荷 - 作用在体积或场域内热分析_ 体积膨胀、内生 束 成热、电磁分析_ magnetic current density等实体模型 FEA 模型惯性载荷 - 结构质量或惯性引起的载荷重力、角速度等 在关键点加集中力在节点加集中力 SSR SSRC C SSR SSRC C 2/ 76 3/ 76 S Space pace S Stru truc ctu ture re R Res esear earc ch h C Center enter, H , HI IT, T, CH CHIN INA A S Space pace S Stru truc ctu ture re R Res esear earc ch h C Center enter, H , HI IT, T, CH CHIN INA A
整车(CAE)有限元建模通用规范
CAE Version 1.0 目录 1. 前言 (2) 2. 命名，编号 (2) 2.1. 概述 (2) 2.2. 特例 (6) 3. 单位 (6) 3.1. 单位制 (6) 3.2. 常用材料基本参数 (6) 4. 坐标系 (7) 5. 分网 (8) 5.1. 1D单元 (8) 5.1.1. 焊点单元(Beam，spotweld， ACM等) (8) 5.1.2. Rigid (8) 5.1.3. Mass (9) 5.2. 2D单元 (9) 5.3. 3D单元 (11) 5.4. 局部特征处理 (13) 5.4.1. 孔 (13) 5.4.2. 加强筋 (15) 5.4.3. 圆角\倒角 (18) 5.4.4. 法兰 (21) 5.4.5. 拼焊板处理（待定） (22) 5.4.6. 其他突出边 (22) 5.4.7. 肋板 (22) 5.5. 连接（TBD） (23) 5.6. 包边 (24) 6. 模型检查 (24) 6.1. 网格质量检查 (24) 6.2. 模型一致性检查 (25) 6.3. 边界条件检查 (25) 6.4. 部件连接检查 (25) 6.5. 整车检查 (25)
1.前言 为了保证有限元模型的通用性，减少重复性工作，特制定本规范，所有零部件建模将依据本规范所规定标准。为了便于管理和维护，现阶段模型采用HyperMehsh v9.0的hm前处理模板，生成*hm格式文件。通用的整车有限元模型包含以下信息：node、element、component、property、assembly等(不包括材料信息)。由于各区域对整车模型材料信息要求不同，共享模型建好之后应用到具体区域的时候再添加材料信息，以碰撞分析为例，专门生成material.k(material.dyn)文件，用include语句进行调用。 为了得到更好的结果，在建模过程中允许不按照本规范建立模型，但是一定要在模型卡片中写出理由，以便于本规范的更新。 2.命名，编号 2.1.概述 (1) 整车模型分为BIW、closure、chassis、trim四个子系统，各子系统又包 含相应的总成，每个总成由若干零件组成。各构成关系(整车—子系统—总成—零件，注意上下级之间的assembly)及编号如表 1所示： (2) 一个零件对应一个component，一个material，一个property，三者ID 号均为一致，同种材料共用同一材料曲线； (3) 零件的命名使用简写后的零件名，并将EPL表格中的零件号注释在comment中，常用词缩写规范如表 2所示。零件名规范为：[零件的名称名词] [零件的描述] [内/外] [前/后] [上/下] [左/右]，举例： 左A柱上内加强板： reinforcement A pillar inner upper left hand，简写后为： reinf A pilr inr upr lh。 右后控制臂支座1： bracket 1 control arm rear right hand，简写后为： brkt 1 ctrl arm rr rh。 表 1：整车构成及编号 子系统 总成 ID区间 part node element (output) 1 1‐9999 1‐9999 (spotweld) 2‐10 10000‐49999 10000‐49999 BIW roof 11‐99 50000‐99999 50000‐99999 front end 100‐199 100000‐199999 100000‐199999 inner side body 200‐299 200000‐299999 200000‐299999 outer side body 300‐399 300000‐399999 300000‐399999 front floor 400‐499 400000‐499999 400000‐499999 rear floor 500‐599 500000‐599999 500000‐599999 预留 600‐999 600000‐999999 600000‐999999 closure hood 1000‐ 1099 1000000‐ 1099999 1000000‐1099999 front door 1100‐ 1199 1100000‐ 1199999 1100000‐1199999 rear side door 1200‐1200000‐1200000‐1299999
有限元分析报告大作业
基于ANSYS软件的有限元分析报告 机制1205班杜星宇U201210671 一、概述 本次大作业主要利用ANSYS软件对桌子的应力和应变进行分析,计算出桌子的最大应力和应变。然后与实际情况进行比较,证明分析的正确性,从而为桌子的优化分析提供了充分的理论依据，并且通过对ANSYS软件的实际操作深刻体会有限元分析方法的基本思想，对有限元分析方法的实际应用有一个大致的认识。 二、问题分析 已知：桌子几何尺寸如图所示，单位为mm。假设桌子的四只脚同地面完全固定，桌子上存放物品，物品产生的均匀分布压力作用在桌面，压力大小等于300Pa，其中弹性模量E=9.3GPa，泊松比μ=0.35，密度ρ=560kg/m3，分析桌子的变形和应力。
将桌脚固定在地面，然后在桌面施加均匀分布的压力，可以看作对进行平面应力分析，桌脚类似于梁单元。由于所分析的结构比较规整且为实体，所以可以将单元类型设为八节点六面体单元。 操作步骤如下： 1、定义工作文件名和工作标题 （1）定义工作文件名：执行Utility Menu/ File/Change Jobname，在弹出Change Jobname 对话框修改文件名为Table。选择New log and error files复选框。 （2）定义工作标题：Utility Menu/File/ Change Title，将弹出Change Title对话框修改工作标题名为The analysis of table。 （3）点击：Plot/Replot。 2、设置计算类型 （1）点击：Main Menu/Preferences,选择Structural,点击OK。