航空航天结构中的有限元方法专题培训(pdf 110页)

合集下载

有限元-第1章

有限元-第1章

L K 1n u1 L K 2n u 2 L L M L K in u i L L M L K nn u n
为了看出该方程能给出所需的结果,取出该方程组中第 i 个方程式
10 30 u i* = K i1u1 + K i 2 u 2 + L + 10 30 u i + L + K in u n
单元节点载荷列阵。
五、 约束处理 由于结构刚度矩阵是奇异的、不能求逆。造成刚度矩阵奇异的原因是在建立刚度矩阵 时,解除了结构的外界约束而成了自由结构,使结构可产生刚体运动。因此必须排除结构 的刚体位移,使结构刚度矩阵成为非奇异的,才能求解出节点位移。一般情况下,引进边 界位移约束条件后可排除刚体运动,否则,还应适当给定某些节点的位移值。经引进边界 约束条件后,可适当地减少待求的节点位移的数目和方程的数目。 在这里介绍一种适合计算机实施的约束处理的作法,称“置大数法” 。设给定的节点 位移为 ui= u i* ,只需在 ui 所在行中, 将结构刚度矩阵的主对角元素 K ij 置入一个大数, 如 10 30 , 同时,将对应行的载荷项 Pi 用 10 30 u i* 代替,于是结构刚度方程成为
{ } { } { } { }
e V
——单元的等效节点载荷列阵, PVe , Pse , P0e 分别是
{ }{ }{ }
e s
单元的体力、面力、节点上的集中力的等效节点载荷列阵,
{P } = ∫∫∫[N ] {f }dv
T Ve
{P } = ∫∫ [N ] {P}ds
T se p
(1-10a)
方程(1-10)表示了结构的总位能是各单元的应变能之和加上各单元等效节点载荷的 位能之和。将单元刚度矩阵 K e 和单元节点载荷列阵 P e 按结构节点位移列阵 {∆} 的自由 度数和排列顺序添零升阶,式(1-10)可进一步完成 1 1 T T T T Π = {∆} ∑ K e {∆} − {∆} ∑ P e = {∆} [K ]{∆} − {∆} {P} 2 2

有限元法课件

有限元法课件

工程有限单元法
2.3 有限元法与工程求解问题的关系
机电工程学院
通常,实际工程问题可分为线性问题和非线性问题、边界 规则与不规则问题。有限元法其实是非线性问题,如图右 所示。
工程有限单元法
三、有限元法的基本步骤
机电工程学院
无论对于什么样的结构,有限元分析过程都 是类似的。其基本步骤为: (1)研究分析结构的特点,包括结构形状与 边界、载荷工况等; (2)将连续体划分成有限单元,形成计算模 型,包括确定单元类型与边界条件、材料特性 等;
工程有限单元法
工程有限单元法 课程介绍
机电工程学院
一、课程内容: 1、有限元法理论基础; 2、有限元软件ANSYS应用。 二、学习方法: 理论与实践相结合,即通过应用有限元分析 实际问题来掌握有限元理论。 三、学时数:36学时(理论学时+上机学时) 四、考核方式:平时成绩+报告成绩
工程有限单元法 第一章 概述
机电工程学院
1.1 有限元法概述 有限元法诞生于20世纪中叶,随着计算机技 术和计算方法的发展,已成为计算力学和计算 工程科学领域里最为有效的方法,它几乎适用 于求解所有连续介质和场的问题。
工程有限单元法
一、什么是有限元法?
机电工程学院
有限元法是将连续体理想化为有限个单元集 合而成,这些单元仅在有限个节点上相连接, 即用有限个单元的集合来代替原来具有无限个 自由度的连续体。
工程有限单元法
二、有限元法的基本思想 有限元法的基本思想是:“分与合”。
机电工程学院
“分”是为了划分单元,进行单元分析; “合”则是为了集合单元,对整体结构进行综合 分析。
结构离散-单元分析-整体求解
工程有限单元法
2.1有限元法的实现过程

有限元法基础讲义

有限元法基础讲义
Tel:(025)4892202-2504 Fax:(025)4895966
有限元法的基本概念
结构离散化: 1)划分网格; 2)载荷移置; 3)简化约束。
单元刚度矩阵与刚度系数: 1)单元刚度矩阵物理意义为单元抵抗变形的能力; 2)刚度系数的物理意义是产生单位位移时需要的力的大小。
南京航空航天大学能源与动力学院机械振动冲击仿真研究室(PC:210016)
南京航空航天大学能源与动力学院机械振动冲击仿真研究室(PC:210016)
Tel:(025)4892202-2504 Fax:(025)4895966
本课程的要求
1. 做好笔记,及时复习与总结 2 . 阅读参考书籍独立上机操作 3 . 独立上机操作
南京航空航天大学能源与动力学院机械振动冲击仿真研究室(PC:210016)
南京航空航天大学能源与动力学院机械振动冲击仿真研究室(PC:210016)
Tel:(025)4892202-2504 Fax:(025)4895966
弹性力学中的基本概念
六个剪应力之间有一定的互等关系。例如,以ab为矩轴,可得:
2 zx y 2 yx z 0
由于单元可以被分割各种形状和大小不同的尺寸,所以它能很 好的适应复杂的几何形状,复杂的材料特性和复杂的边界条件,再 加上它有成熟的大型软件系统支持,使它已成为一种非常受欢迎的, 应用极广的数值计算方法。
南京航空航天大学能源与动力学院机械振动冲击仿真研究室(PC:210016)
Tel:(025)4892202-2504 Fax:(025)4895966
lim Q F (N/m2)
S S0
南京航空航天大学能源与动力学院机械振动冲击仿真研究室(PC:210016)

北航有限元分析与应用第一讲

北航有限元分析与应用第一讲

A() (k ) (k )Q 0 x x y y
• 边界条件:

B() { k q 0 q上 n
0 上
数值计算方法分类
特点
差分法
优缺点
离散求解域;差分代替微分; 要求规则边界,几何 解代数方程组 形状复杂时精度低
1-5 有限单元法的形成与发展
数学家们则发展了微分方程的近似解法,包括有限差分方法,变 分原理和加权余量法。 在 1963 年 前 后 , 经 过 J.F.Besseling, R.J.Melosh, R.E.Jones, R.H.Gallaher, T.H.Pian(卞学磺)等许多人的工作,认识到有限元法 就是变分原理中Ritz近似法的一种变形,发展了用各种不同变分原理 导出的有限元计算公式。 1965 年 O.C.Zienkiewicz 和 Y.K.Cheung (张佑启)发现只要能写成 变分形式的所有场问题,都可以用与固体力学有限元法的相同步骤 求解。 1969年B.A.Szabo和G.C.Lee指出可以用加权余量法特别是 Galerkin 法,导出标准的有限元过程来求解非结构问题。
v1 2 Fy12 u

1 2
2
X2 ②

1

3
Fy22 Fx12
• 节点1沿x方向的位 1 u 移 1 1、其余节点位 F 移全为0时轴向压力v 为: EA EA cos 1
1 1
2
2

Fx22
1 y1
Fy23 u
1 1
(
l1
)l1
Fx11
l1
3
Fx23
实例1(单元分析)
• 节点1作用于单元1上的力,在x和y方向的分量分别为:

有限元方法PPT文档共32页

有限元方法PPT文档共32页
有限元方法
56、死去何所道,托体同山阿。 57、春秋多佳日,登高赋新诗。 58、种豆南山下,草盛豆苗稀。晨兴 理荒秽 ,带月 荷锄归 。道狭 草木长 ,夕露 沾我衣 。衣沾 不足惜 ,但使 愿无违 。 59、相见无杂言,但道桑麻长。 60、迢迢新秋夕,亭亭月将圆。
谢谢你的阅读
ቤተ መጻሕፍቲ ባይዱ
❖ 知识就是财富 ❖ 丰富你的人生
71、既然我已经踏上这条道路,那么,任何东西都不应妨碍我沿着这条路走下去。——康德 72、家庭成为快乐的种子在外也不致成为障碍物但在旅行之际却是夜间的伴侣。——西塞罗 73、坚持意志伟大的事业需要始终不渝的精神。——伏尔泰 74、路漫漫其修道远,吾将上下而求索。——屈原 75、内外相应,言行相称。——韩非

北航有限元第一讲

北航有限元第一讲
北京航空航天大学
1965 年 O.C.Zienkiewicz 和 Y.K.Cheung ( 张 佑 启 ) 发现只要能写成变分形式的所有场问题,都可 以用与固体力学有限元法的相同步骤求解。 1967 1967年,Zienkiewicz和Cheung出版了第一本有 Zienkiewicz Cheung 关有限元分析的专著。
EA(un+1 − un ) qLn Nn = A n = σ = Ln 2 qL2 −un + un+1 = n 2EA
再加上约束条件
u1 = 0
因此可以得到n+1个方程构成的方程组,可解 出n+1个结点的位移。
北京航空航天大学
有限元方法的基本思想和原理是“简单”而 “朴素”的,在发展初期,许多学术权威对该 方法的学术价值有所鄙视,国际著名刊物 Journal of Applied Mechanics许多年来拒绝刊登 有关有限元方法的文章,其理由是没有新的科 学实质。 现在完全不同了,由于有限元方法在科学研究 和工程分析中的地位,有关有限元方法的研究 已经成为数值计算的主流。涉及有限元方法的 杂志有几十种之多。
EA(ui+1 − ui ) EA(ui+2 − ui+1 ) q − = (Li + Li+1 ) Li Li+1 2
Li 令 λi = Li+1
q 1 2 − ui + (1+ λi )ui+1 − λiui+2 = (1+ )Li 2EA λi
北京航空航天大学
对于第n+1个结点,第n个单元的内力与第 n+1个结点上的外载荷平衡,
北京航空航天大学
工程师方面

有限元法PPT课件

有限元法PPT课件
和时间。
如何克服局限性
改进模型
通过更精确地描述实际 结构,减少模型简化带
来的误差。
优化网格生成
采用先进的网格生成技 术,提高网格质量,降
低计算误差。
采用高效算法
采用并行计算、稀疏矩 阵技术等高效算法,提
高计算效率。
误差分析和验证
对有限元法的结果进行误 差分析和验证,确保结果
的准确性和可靠性。
05 有限元法的应用实例
有限元法ppt课件
目 录
• 引言 • 有限元法的基本原理 • 有限元法的实现过程 • 有限元法的优势与局限性 • 有限元法的应用实例 • 有限元法的前沿技术与发展趋势 • 结论
01 引言
有限元法的定义
01
有限元法是一种数值分析方法, 通过将复杂的结构或系统离散化 为有限个简单元(或称为元素) 的组合,来模拟和分析其行为。
有限元法在流体动力学分析中能够处理复杂的流体流动和 压力分布。
详细描述
通过将流体域离散化为有限个小的单元,有限元法能够模 拟流体的流动、压力、速度等状态,广泛应用于航空、航 天、船舶等领域。
实例
分析飞机机翼在不同飞行状态下的气动性能,优化机翼设 计。
热传导分析
总结词
有限元法在热传导分析中能够处理复杂的热传递过程。
实例
分析复杂电磁设备的电磁干扰问题,优化设备性能。
06 有限元法的前沿技术与发 展趋势
多物理场耦合的有限元法
总结词
多物理场耦合的有限元法是当前有限元法的重要发展方向, 它能够模拟多个物理场之间的相互作用,为复杂工程问题提 供更精确的解决方案。
详细描述
多物理场耦合的有限元法涉及到流体力学、热力学、电磁学 等多个物理场的耦合,通过建立统一的数学模型,能够更准 确地模拟多物理场之间的相互作用。这种方法在航空航天、 能源、环境等领域具有广泛的应用前景。

航空整体结构件加工变形有限元数值仿真

航空整体结构件加工变形有限元数值仿真

2009年2月第35卷第2期北京航空航天大学学报Journa l o f Be iji ng U nivers it y of A eronauti cs and A stronauti cs F ebruary 2009V o.l 35 N o 2收稿日期:2008 07 20基金项目:国防基础科研资助项目作者简介:张以都(1959-),男,北京人,教授,yd z hang @buaa .edu .cn.航空整体结构件加工变形有限元数值仿真张以都 张洪伟(北京航空航天大学虚拟现实技术与系统国家重点实验室,北京100191)摘 要:加工变形是航空整体结构件数控加工领域普遍存在的问题.选取典型航空整体结构件,基于有限元法研究了毛坯初始残余应力对于加工变形的影响,仿真结果与试验结果一致;对残余应力引起的加工变形特点进行了分析,基于不同加工路径对航空结构件加工变形进行了有限元数值仿真,结果表明,对于单一残余应力引起的加工变形分析,可以忽略加工路径的影响.考虑铣削力与残余应力的耦合作用,建立了残余应力、铣削力的耦合力学模型,给出了具体仿真流程,针对不同加工路径进行有限元仿真,研究了铣削力对于残余应力分布及最终加工变形的影响.关 键 词:残余应力;加工变形;有限元;铣削力;加工路径中图分类号:TP 391.77文献标识码:A 文章编号:1001 5965(2009)02 0188 05Finite e l e ment sm i ul a ti o n ofmachi n ing def or mati o n f oraeronauti c a lmono lit hic co mponentZhang Y idu Zhang H ong w ei(S tate K ey Laboratory ofV i rt ua lReali ty T echnology and Syste m s ,B eiji ng Un i versity ofA eronau tics and A stronauti cs ,Beiji ng 100191,Ch i na)Abstr act :M achining defor m ation has been a bo ttle neck issue in aeronauticalm anufacturi n g fie l d .B ased on finite ele m ent si m u lation ,the i n fl u ence of or i g i n al residua l stresses on the defo r m ation in the process o fm a ch i n i n g fo r aeronauticalm onolithic co m ponent w as stud ied .The red istr i b ution la w of resi d ual stresses during m ach i n i n g and the trend o f defor m ation change w ere analyzed .The si m ulation conc l u si o ns ag ree w it h experi m ental resu lts .The m achining defor m ation characteristics due to resi d ual stressesw ere analyzed and t h e si m u lation of differentm achini n g path w as carried ou.t The resu lts i n d icate that the research on defor m ation caused by residua l stresses a lone can neglect the effect o fm ach i n i n g path .The coupling e ffect bet w een residual stress and m illing force w as taken i n to account and the si m ulati o n w as perfor m ed for d ifferent m achi n i n g pat h .Key w ords :residual stress ;m ach i n i n g d istortion ;finite ele m ent m ethod;m illi n g force ;m ac h i n i n g path在整体航空结构件的铣削加工过程中普遍存在有加工变形的问题,这是数控加工领域公认的难题.针对加工变形问题,学者们在各个领域进行了广泛研究,如切削加工过程模拟、残余应力数值模拟等,其在一定程度上为优化工艺提供了理论依据.这些文献主要集中在二维切削加工机理的研究或主要研究引起加工变形的某一因素上,这对于预测航空结构件的整体变形有一定局限性[1-3].研究表明,毛坯初始残余应力是引起加工变形的重要原因[4-5].在实际加工过程中,铣削力是影响加工变形的另一因素.现有文献常常不计铣削力与初始残余应力的耦合作用,而且铣削加工是一个复杂的三维动态过程[6],在有限元数值模拟中有待对铣削加工过程的关键技术进行系统深入研究.本文基于有限元仿真技术研究了初始残余应力对于整体结构件的加工变形的影响,仿真结果与试验结果一致,并对残余应力引起的加工变形特点进行了分析;建立了包含初始残余应力、铣削力的耦合力学有限元模型,基于不同加工路径进行有限元仿真,对铣削力在加工变形中的影响作了深入研究,探索耦合条件下框类航空结构件加工变形规律.1 残余应力对加工变形的影响分析1.1 有限元模型及初始残余应力场构建选取典型整体结构件按相似性原理设计比例件进行加工过程有限元仿真.工件材料为7075铝合金,经过淬火热处理.毛坯尺寸:300mm 84mm 26mm,工件壁厚为2mm.弹性模量为71GPa ,泊松比为0.3,塑性模量为250M Pa ,屈服应力为455M Pa ,选用双线性强化模型来描述材料本构关系.初始残余应力分布如图1所示[7]包括长度方向(X 向)的残余应力及宽展方向(Y 向)的残余应力.编写残余应力施加算法,对有限元模型添加初始残余应力[8].工件理想设计形状如图2所示.图1 初始残余应力分布规律图2 工件理想形状1.2 仿真过程仿真边界条件是在工件的非加工部分选择3个点,分别约束X YZ,YZ 及Z 平移自由度.3点选择原则是:完全约束工件的刚体运动,而不影响工件的任意变形.按照加工过程,将要去除的单元按加工顺序排列,一层一层的去掉,每铣掉一层单元,将其作为一个载荷步计算一次.模拟过程中,将铣削去除部分分成11层,分11次去除,共经过11个载荷步计算,完成该整体结构件的模拟加工.1.3 结果分析分别施加如图1所示的单向(X 向)残余应力及双向(X 向、Y 向)残余应力,进行仿真计算.结果表明,在残余应力作用下,工件发生明显的变形,中部鼓起,两端上翘,整个工件向着被去除材料的一侧弯曲.试件侧壁高度较小,变形量相对较小,以底面的变形最为严重,底面的翘曲变形导致侧壁的弯曲变形.从工件选取如图2所示A,B,C 3个变形关键点,其中A 点是中部框中部的点,也是底面长度方向的中点,B 点是开口端部框的中部的点,C 点是另一端框上部的一个点.在单向残余应力作用下,随材料去除过程,3点变形量的变化趋势如图3所示.图3 单向残余应力作用下剥层过程中各点变形变化趋势在双向残余应力作用下,选择同样的3点,绘制其变形量随材料的去除变化曲线,如图4所示.图4 双向残余应力作用下剥层过程中各点变形变化趋势无论是单向残余应力还是双向残余应力作用,在材料去除过程中,都会在去除量达50%左右时变形量达到最大值,然后逐步减小,并趋向稳定,但施加双向残余应力引起的加工变形值比只189第2期 张以都等:航空整体结构件加工变形有限元数值仿真施加沿长度方向的残余应力所引起的变形值要小,这是由于宽展方向的残余应力会在一定程度上牵制变形.在变形量变化较大的前几次,残余应力变化都较大,去除层均为压应力层.而去除拉应力层过程中,由于拉应力层的反向变形,使得变形减小,同时由于拉应力层变化梯度较小,所以变形变化量也较小.仿真结果的变化规律得到文献[9]中的试验结论的支持.1.4 试验结果与仿真结果对比针对一框类试件进行铣削加工试验及仿真计算.图5为加工完成后的工件,其变形为:两端向上翘起,中间向下凹陷.最大变形量发生在试件长度方向的端部,数值大约为0.45mm.图6为试件底面长度方向的一组节点试验结果与仿真结果的变形曲线对比.结果均表明试件图5 试验结果端部翘曲变形a试验结果b 仿真结果图6 底部沿长度方向的一组节点的变形比较底面变形为两端上翘,中间部分下陷.工件最大翘曲变形发生在工件的端部;向下的最大凹陷发生在长度方向左右腔型的中间,这是因为该位置是整个工件中刚度最薄弱的地方.仿真结果大于试验结果,原因在于试验所用的板材一般是经过时效处理的,初始应力有部分释放,而仿真方法无法模拟时效处理中的应力松弛现象,所以仿真时板内初始残余应力是按理想情况添加的,会高于试验所用板内初始残余应力.2 残余应力引起加工变形特点分析初始残余应力引起的加工变形主要包括两部分:一是工件整体的翘曲变形;另一部分是由于残余应力的不断释放与重分布,造成局部的 过切 或 欠切 铣削误差变形.真实的加工变形是两种变形的综合.研究表明,残余应力引起的加工变形与工件材料的加工顺序及刀具走刀路径等有关.对某单框类结构件选用4种常用的走刀加工路径进行有限元仿真研究: 由外向内 回 字型走刀; 由内向外 回 字型走刀; 沿轧制方向 之 字型走刀; 沿横向 之 字型走刀.在残余应力作用下,仿真结果表明4种不同走刀路径对工件整体宏观变形没有任何影响.这是由于残余应力是弹性应力,以弹性势能的形式储存在工件体内,随材料去除,工件能量不断变化,由弹性力学理论可知,这种弹性能量的变化只与最终的工件形状有关.也就是说,初始残余应力引起的工件整体变形与具体的铣削过程无关.另外,随残余应力的释放,造成了局部的铣削误差变形.这种铣削误差变形幅值较小,属于弹性变形.所以研究整体工件变形时,铣削误差变形影响是很小的.根据这一结论,如果单一考虑初始残余应力对于工件整体变形的影响,则不必逐步进行铣削过程仿真,而可以将要去除的材料单元一次完全 杀死 ,这样可以大大提高计算效率.3 铣削力与残余应力耦合作用分析建立包含初始残余应力、铣削力和装夹条件的框类结构件耦合力学模型,探索多因素耦合条件下框类结构件加工变形规律.铣削加工过程是铣削力不断加载和卸载的过程.通过动态给拟去除材料单元进行加载、卸载分析来获得残余应力及加工变形.加载过程是指在铣削过程中,由于铣削力的作用,而在工件局部发生塑性变形产生残190北京航空航天大学学报 2009年余应力的过程,3向铣削力以单元面载荷形式作用在指定单元上,并沿着一定的走刀路径进行.卸载过程是指切屑层被去除、刀具离开当前工件位置,工件内部残余应力重新分布这一过程,即对所加载单元施加3向零载荷进行求解.由于载荷是一个动态变化的过程,采用APDL 语言编制运行程序,实现铣削过程.其算法流程如图7所示.图7 铣削力、残余应力耦合作用仿真流程基于上述流程,进行不同加工路径有限元加工仿真.选取某典型节点长度方向残余应力随动态铣削力的变化的响应曲线,如图8所示.仿真结果表明,4种加工路径下的最终残余应力差别较小,但铣削过程中会有较大的波动.其中路径1所引起的变形波动最小,铣削过程较平稳,这个结论得到了文献[10]的支持.加工过程的波动会引起局部的铣削误差变形.这种变形是由于铣削力的作用,造成局部的 让刀 变形,表现为 多切 或 欠切 .图9是不考虑初始残余应力影响,在铣削力作用下不同加工路径典型节点的应力变化情况,相比考虑残余应力影响,这种情况下,应力幅值小得多.对于整体宏观变形而言,铣削力的影响相对图8 不同走刀路径典型节点应力变化图9 仅考虑铣削力不同走刀路径节点应力变化于毛坯体的初始残余应力是很小的.从根本上来说,高速加工铣削力难以影响到最终的工件宏观变形,是由于铣削力的作用不足以改变残余应力分布状态,或者说铣削力产生的应力场在与原始的残余应力场发生相互作用后已经被覆盖了,即铣削力要起到较大的作用,必须首先产生较大的应力场与初始残余应力场进行耦合作用后,仍然保持相当的数值,才有可能影响到最终的应力场及最终的变形.在对航空结构件的整体变形研究中,在特定加工路径下,铣削力的作用是可以不考虑的.但在薄壁件的加工过程中,由于初始残余应力的释放,铣削力的作用需要得到关注.4 结 论1)加工过程中,双向残余应力引起的加工变形量与单向残余应力引起的加工变形量有较大差异,但残余应力分布相似.加工变形过程中,存在着变形量临界值,在该临界值处,工件变形达到最大值,随后工件变形量逐渐减小并趋向稳定.2)加工路径不同对于工件的残余应力的释放与重分布有不同影响,但若只考虑残余应力单一因素,工件加工变形只与材料去除量的多少及工件最终结构形式有关,而与加工路径、加工顺序191第2期 张以都等:航空整体结构件加工变形有限元数值仿真等去除过程无关.在仿真运算过程中,为提高计算效率,可以一次全部去除掉.3)考虑残余应力与铣削力的耦合影响,仿真结果表明,铣削力的存在会影响残余应力的释放和重分布.但高速铣削力所产生的残余应力只存在于很薄的一层内,而且时间短暂,所以对初始残余应力的影响还是比较小的,从而对最终的加工变形影响也很小.在确定加工路径的基础上,可以忽略铣削力的影响,提高计算效率,而不会影响宏观整体变形研究.4)基于有限元仿真结果,通过采取某些工艺措施减小残余应力或均匀化残余应力,或在加工过程中,基于加工变形临界值预留一定加工余量,在材料去除达到变形量稳定值时,反向加工去除材料,以控制加工变形.参考文献(References)[1]E l Ax i rM H.A m ethod of modeli ng residu al s tress d i stribu tioni n t u rn i ng for d ifferentm aterial s[J].In t ernati onal Jou rnal ofM ach i ne Tool s&M anu fact u re,2002,42:1055-1063[2]F i hri Fass iH,Bou ss h i ne L,Chaaba A,et a.l Nu m eri cal si m ulation of orthogonal cu tti ng by i n cre m en tal el ast oplas tic analys i s and fi n ite el e m en t m et hod[J].Jou r n al ofM aterial s Processing Techno l ogy,2003,14:181-188[3]M arti n B aker.F i n it e el e m en t s i m u lati on of high speed cu ttingforces[J].Jou rnal of M at eri als Processi ng Technol ogy,2006,176:117-126[4]B ri nm s m eier E,H annover J,Ca mm ett T,et a.l R es i dual stressm easure m en t and cause i n m ach i n i ng p rocess[J].Ann al s of theC I RP,1982,31(2):491-509[5]H e i z A,H aszler A,KeidelC,et a.l Recen t devel opm en t i n al um inum all oy f or aeros pace app li cati on[J].M ater SciE ngA,1999, A280:102-107[6]黄志刚,柯映林.飞机整体框类结构件铣削加工的模拟研究[J].中国机械工程,2004,15(11):991-995H uang Zh igang,Ke Y i ngli n.Study on key technologies ofm illi ngp rocess si m u l ation for aerospacem onoli th i c co m pon ents[J].Ch i na M echan ical Engi n eeri ng,2004,15(11):991-995(i n Ch i n ese)[7]王秋成,柯映林,章巧芳.7075铝合金板材残余应力深度梯度的评估[J].航空学报,2003,24(4):336-338W ang Q i ucheng,Ke Y i ngli n,Zhang Q i aof ang.E val u ati on of re si dual stress dep t h profili ng i n7075al um i nu m a ll oy p l ates[J].A ct a Aeronau ti ca et A stronauti ca S i n ica,2003,24(4):336–338(in Ch i nese)[8]张洪伟,张以都,赵晓慈,等.航空结构件加工变形仿真关键技术[J].北京航空航天大学学报,2008,34(2):239-243Zhang H ong w e,i Zhang Y i du,Zhao X i aoc,i et a.l Research on k ey tec hn iqu es i n the si m u l ati on ofm ach i n i ng d i storti on f or aero nauti calm ono lit h ic co m ponen t[J].Jou rnal of B eiji ng Un i versity ofA eron auti cal and Astronautics,2008,34(2):239-243(i n Ch i nese)[9]董兆伟.铝合金整体结构件铣削加工残余应力及变形研究[D].北京:北京航空航天大学机械工程及自动化学院,2006Dong Zhao w e.i S t udy on the m illi ng res i dual stress and def or m a tion of the al um i nu m alloy i n t egrated fra m e w ork[D].Beiji ng:S chool ofM echan i cal Engi n eeri ng and Auto m ation,Be iji ng Un i vers i ty ofA eronau tics and A stron auti cs,2006(i n Ch i nese) [10]王立涛.关于航空框类结构件铣削加工残余应力和变形机理的研究[D].杭州:浙江大学机械与能源工程学院,2003W ang L it ao.S t udy on res i dual stresses and distorti on t h eory of aeronau tical fra m e struct u re i n t he m illi ng[D].H angz h ou:C ollege ofM echan i cal and E nergy Engi n eeri ng,Zh eji ang U n i vers it y,2003(i n C h i nese)192北京航空航天大学学报 2009年。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

航空航天结构中的有限元方法邱志平 王晓军 编著北京航空航天大学2011年1月前言有限元方法经过半个多世纪的发展,现已成为当今工程问题中应用最广泛的数值计算方法。

有限元方法集多学科理论知识于一身,且有着自己的理论基础和解题方法。

有限元方法首先被航空结构工程师引入并发展,并由于其在解决工程技术问题时的灵活、快速及有效性,发展非常迅速,现在其解题范围几乎渗透到了的各个研究领域,包括固体变形场、流体场、电磁场、温度场和声场等。

近年来,由于有限元分析商业化软件的普及,有限元分析不再只为少数专业人员所掌握,转而成为高校、科技工作者和工程技术人员所广泛使用的通用分析工具。

拥有了先进的和自动化的有限元分析软硬件平台,并不意味着就掌握了有限元分析方法和能够得到正确的分析结果。

对于实际的工程结构,特别是航空航天领域复杂的组合结构,工作环境复杂严峻,技术要求苛刻,要取得合乎工程标准的可信的结构分析结果,需要工程技术人员具有较高的理论素养和实际经验。

本书拟作为航空航天院校及相关专业的教学参考书。

目前国内介绍有限元法的书籍很多,与同类教材相比,本书具有以下特点:首先,为了使本书具有相对的系统性与完整性,前几章介绍了有限元方法的基本原理和理论基础,内容简洁而重点突出,为后面的实例分析奠定基础;其次,针对航空航天领域的典型结构,如机身结构、机尾翼结构和起落架结构等,详细介绍了其有限元模型的建立、网格划分、边界条件的选取以及载荷施加、求解和结果分析,并把常见的杆件结构、板和壳体问题分析融入到航空航天典型结构的分析当中;第三,介绍了航空航天领域的突出问题的工程应用:静力分析,动力分析和复合材料结构的有限元分析技术等。

本书内容完整,具有显著的航空航天特色,可作为航空航天、力学、机械等专业学生的教材,也可作为上述专业教师和工程技术及科研开发人员的参考书。

鉴于工科院校学生理论知识与工程应用的并重性,本书分为上下两篇:第一篇主要介绍了有限元法的基本理论,共分九章。

第一章绪论部分概述了有限元法的早期工作、现状和未来,并简述了有限元在航空航天领域的应用背景和分析方法。

第二章介绍了有限元法的基本概念和理论基础,第三章至第五章分别介绍了杆系结构,平面问题和空间问题的有限元法,第六章和第七章分别介绍了轴对称问题和板壳问题的有限元法,第八章介绍了高阶单元与等参单元的理论,第九章给出了结构动力有限元分析的概念和基本理论,并在第十章介绍了作者在不确定性结构分析的区间有限元方法方面的一些工作。

第一篇中的理论部分具有普适性,既为第二篇所述的航空航天有限元奠定基础,又可作为相关专业人士在有限元理论方面的入门参考。

第二篇主要介绍了有限元法在航空航天中的应用,分为航空部分和航天部分,目的是为相关专业人士在使用有限元法分析航空航天结构时提供参考。

航空部分以大型飞机结构为原型,介绍了模型简化,约束与载荷的施加等,并分别就飞机的机身、机翼、舵面、发动机短舱、气密舱、座舱盖、货舱门、机身机翼对接区和起落架部分的有限元分析方法进行了详细的介绍。

航天部分主要介绍了卫星结构分析中的有限元法,包括有限元模型的建立,静力分析和模态分析。

最后,作者将近年来在有限元方面所作的工作及在工程中的具体应用做了较为详细的阐述,以期读者能够将理论与实际融会贯通,正确并熟练应用有限元分析方法进行工程结构分析。

作者长期以来一直从事航空航天飞行器结构分析设计研究和教学工作,积累了丰富的科研教学经验。

目前国内还没有一本书能够全面系统地介绍空天飞行器结构有限元方法的基础性和专门性的教材,本书的出版适应于我国国防类大学学生的需要。

在编写本书的过程中,博士研究生祁武超、吴迪等付出了许多辛勤的劳动,仔细录入、校审了初稿,作者在此表示衷心地感谢!感谢本教材编写所参考书目的作者,书中如有漏引之处,还请作者见谅!由于作者水平所限,书中一定存在某些不妥和需要修改之处,敬请读者批评指正。

作者 2011年1月目录第1章绪论 ..................................................................................................................1.1 有限元的发展、现状和未来.................................................................................................1.1.1 有限元法的早期工作.................................................................................................1.1.2 有限元法的发展和现状.............................................................................................1.1.3 有限元法的未来.........................................................................................................1.2 有限元法在航空航天中的应用.............................................................................................1.2.1 有限单元法在航空航天领域中的应用背景 .............................................................1.2.2 航空结构分析.............................................................................................................1.2.3 航天结构分析.............................................................................................................习题................................................................................................................................................. 第2章有限单元法的基本概念和理论基础 ..................................................................2.1 有限元法的基本思想.............................................................................................................2.2 有限元法的基本概念.............................................................................................................2.2.1 结构离散(有限元建模).........................................................................................2.2.2 插值函数.....................................................................................................................2.2.3 有限元法的收敛准则.................................................................................................2.2.4 有限元法的基本步骤.................................................................................................2.2.5 场问题的一般描述.....................................................................................................2.3 弹性力学基本理论.................................................................................................................2.4 变分原理...............................................................................................................................2.5 有限元平衡方程...................................................................................................................习题............................................................................................................................................... 第3章杆系结构有限元分析 ........................................................................................3.1 拉压杆单元...........................................................................................................................3.1.1 一般规定...................................................................................................................3.1.2 位移函数...................................................................................................................3.1.3 几何关系和物理关系...............................................................................................3.1.4 平衡关系...................................................................................................................3.1.5 坐标变换...................................................................................................................3.2 扭转杆单元...........................................................................................................................3.3 平面直梁单元.......................................................................................................................3.3.1 位移函数...................................................................................................................3.3.2 梁元的刚度矩阵.......................................................................................................3.3.3 坐标变换...................................................................................................................3.3.4 等效结点载荷...........................................................................................................3.4 总体刚度矩阵和总体载荷列向量.......................................................................................3.5 刚度矩阵的物理意义和性质...............................................................................................3.6 位移边界条件.......................................................................................................................3.7 总刚度平衡方程的求解.......................................................................................................3.8 算例.......................................................................................................................................3.9 小结.......................................................................................................................................习题............................................................................................................................................... 第4章平面问题有限元分析 ........................................................................................4.1 引言.......................................................................................................................................4.2 常应变三角形单元...............................................................................................................4.2.1 离散化.......................................................................................................................4.2.2 位移模式与形函数...................................................................................................4.2.3 基于最小势能原理的单元特性分析 .......................................................................4.3 单元等效结点载荷列阵.......................................................................................................4.4 矩形双线性单元...................................................................................................................4.4.1 位移模式与形函数...................................................................................................4.4.2 单元刚度矩阵和单元等效载荷列阵 .......................................................................4.4.3 单元等效结点荷载矩阵...........................................................................................4.5 应力计算结果的整理...........................................................................................................习题............................................................................................................................................... 第5章空间问题有限元分析 ........................................................................................5.1 三维应力状态.......................................................................................................................5.2 四面体常应变单元...............................................................................................................5.3 直六面体单元.......................................................................................................................习题............................................................................................................................................... 第6章轴对称问题的有限元分析 ................................................................................6.1 单元位移函数.......................................................................................................................6.2 单元应变与应力...................................................................................................................6.2.1 单元应变...................................................................................................................6.2.2 单元应力...................................................................................................................6.3 单元刚度矩阵.......................................................................................................................6.4 整体刚度矩阵.......................................................................................................................6.5 等效结点载荷.......................................................................................................................6.5.1 体积力.......................................................................................................................6.5.2 表面力.......................................................................................................................习题............................................................................................................................................... 第7章板壳问题有限元分析 ........................................................................................7.1 薄板问题的有限元法...........................................................................................................7.1.1 矩形单元的位移函数...............................................................................................7.1.2 矩形单元的刚度矩阵...............................................................................................7.1.3 矩形单元的等效结点载荷和内力矩 .......................................................................7.2 薄壳问题的有限元法...........................................................................................................7.2.1 结构载荷列阵...........................................................................................................7.2.2 单元刚度矩阵...........................................................................................................7.2.3 结点应力计算...........................................................................................................习题............................................................................................................................................... 第8章高阶单元与等参数单元 ....................................................................................8.1 高阶单元...............................................................................................................................8.1.1 建立形函数的方法...................................................................................................8.1.2 多项式的完备性.......................................................................................................8.1.3 矩形单元—Lagrange族单元...................................................................................8.1.4 矩形单元—Serendipity族单元 ...............................................................................8.2 平面四结点等参元...............................................................................................................8.2.1 坐标变换与等参单元...............................................................................................8.2.2 单元刚度矩阵的计算...............................................................................................8.2.3 等参变换的条件和等参单元的收敛性 ...................................................................8.3 八结点曲边等参单元...........................................................................................................8.3.1 位移函数...................................................................................................................8.3.2 等参单元等效结点力...............................................................................................习题............................................................................................................................................... 第9章结构动力有限元分析 ........................................................................................9.1 动力问题有限元列式-基本概念..........................................................................................9.2 运动方程式...........................................................................................................................9.2.1 惯性力和阻尼力.......................................................................................................9.2.2 运动方程的建立.......................................................................................................9.2.3 动力方程与静力方程的区别...................................................................................9.3 质量矩阵...............................................................................................................................9.3.1 集中质量矩阵...........................................................................................................9.3.2 一致质量矩阵...........................................................................................................9.4 阻尼矩阵...............................................................................................................................9.4.1 单元阻尼矩阵...........................................................................................................9.4.2 总体阻尼矩阵...........................................................................................................9.5 无阻尼自由振动分析-特征值问题......................................................................................9.6 振型的性质...........................................................................................................................9.6.1 振型的规格化...........................................................................................................9.6.2 振型的正交性...........................................................................................................9.7 有阻尼的自由振动分析.......................................................................................................9.8 结构动力响应分析...............................................................................................................9.8.1 振型叠加法...............................................................................................................9.8.2 直接积分法...............................................................................................................习题............................................................................................................................................. 参考文献............................................................................................................................主要符号表Ω求解域Γ求解域边界,,X Y Z表面力分量,,X Y Z体积力分量σ正应力τ剪应力ε正应变γ剪应变σ应力向量ε应变向量E弹性模量G剪切模量μ泊松比I截面惯性矩l杆长A横截面积ϕ剪切影响系数sA有效抗剪面积D弹性矩阵B单元应变矩阵S单元应力矩阵λ方向余弦矩阵T坐标变换矩阵I单位矩阵∏势能泛函W总功U内力虚功V外力虚功,单元体积(),,,rN r i j k= 单元形函数N单元形函数矩阵''''o x y z局部(单元)坐标系oxyz整体坐标系,,u v w位移场分量δ位移场向量,,,,,ii i i xi y ziu v wθθθ结点位移分量iu结点位移向量eK单元刚度矩阵eu单元结点位移向量ef单元结点载荷向量K总体刚度矩阵u总体结点位移向量f总体结点载荷向量,,,i j k 结点编号eF单元结点力P集中载荷q分布载荷T扭矩M弯矩Q剪力N 轴力θ 扭转角ν 挠度V 单元体积∆ 三角形单元面积Λ 6倍四面体单元体积M 质量矩阵C 阻尼矩阵Φ 振型矩阵Λ 特征值对角矩阵12,μμ 材料阻尼系数λ 特征值φ 特征向量,t τ 时间 ω 角频率 ,,b αβ 不确定性参数 I x 区间数 c x 区间数I x 的中点 x ∆ 区间I x 的半径 x 区间I x 的上界 x 区间I x 的下界 注:局部坐标系下的量使用上角标符号(’)第1章绪论1.1有限元的发展、现状和未来1.1.1有限元法的早期工作有限元法的根本思想可追溯到Courant在1943年的工作,他首先将在一系列三角形区域上定义的分片连续函数和最小势能原理相结合,来求解St. Venant扭转问题。

相关文档
最新文档