航空航天结构中的有限元方法专题培训(pdf 110页)
有限元-第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有限元法的实现过程
有限元法基础讲义

有限元法的基本概念
结构离散化: 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页

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课件

如何克服局限性
改进模型
通过更精确地描述实际 结构,减少模型简化带
来的误差。
优化网格生成
采用先进的网格生成技 术,提高网格质量,降
低计算误差。
采用高效算法
采用并行计算、稀疏矩 阵技术等高效算法,提
高计算效率。
误差分析和验证
对有限元法的结果进行误 差分析和验证,确保结果
的准确性和可靠性。
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
航空航天结构中的有限元方法 编程大作业
for i=1:2*nn data(10,i)=F(1,i) end for i=1:2*nn data(11,i)=bc(1,i) end filename=input('请输入保存的文件名\n','s') fid=fopen(filename,'w'); fprintf(fid,'%f %f %f %f %f %f %f %f %f %f %f %f\n',data) fclose(fid)
三、 计算程序
编写的计算程序如下, matlab 对矩阵的操作十分方便, 首先读入数据并将它们存入一个 名为 data 的矩阵,再用这个矩阵对各参数进行赋值,组装总体刚度矩阵也可以很轻松地进 行,因为编程时设置的字长比较短,为确保计算精度,在引入边界条件时用了置 1 法。 filename=input('请输入生成的数据文件的文件名\n','s') fid=fopen(filename,'r'); [data,count]=fscanf(fid,'%f %f %f %f %f %f %f %f %f %f %f',[11 inf]); fprintf(1,'%f %f %f %f %f %f %f %f %f %f %f',data); fclose(fid) %读入数据 nn=data(1,1) nm=data(2,1) for i=1:nm L(i,1)=data(3,i) end E=data(4,1) A=data(5,1) for i=1:nm Q(i,1)=data(6,i) end 2
运行这个.m 文件,在指引下输入各项数据,最后能在工作目录下生成文件名可以自定 义的 txt 数据文件。 由于时间仓促, 编写时设置的输入量写成了杆件长度和杆件的倾斜角, 今后可对这些进 行改进,输入节点坐标,然后由 L2 = (������2 − ������1 )2 + (������2 − ������1 )2 ������2 − ������1 Q = arctan ������2 − ������1 可得出相应杆件的长度和倾角。 另外, 也可以不使用这个小程序生成数据文件, 只要存储数据时按照格式存成一个矩阵, 计算程序就能读取相关数据进行计算。
有限元基本原理与概念培训课件

离散化的目的
将复杂的连续系统简化为 易于分析的离散模型,以 便进行数值计算和分析。
离散化的方法
根据实际问题的需求,可 以采用不同的离散化方法, 如四面体离散化、六面体 离散化等。
单元选择与形状函数
单元选择
选择合适的单元类型以逼 近真实形状,常用的单元 类型有四面体、六面体、 板壳等。
形状函数
描述单元内节点位移与单 元位移之间关系的函数, 用于建立节点位移与整体 位移之间的关系。
形状函数的性质
满足完备性和协调性,以 保证整体位移的连续性和 一致性。
刚度矩阵与载荷向量
刚度矩阵
描述节点力与节点位移之间关系 的矩阵,由单元刚度矩阵组装而
成。
载荷向量
作用在系统上的外力向量,包括集 中载荷、分布载荷等。
平衡方程
通过建立整体刚度矩阵和载荷向量 的平衡方程,可以求解节点的位移。
位移求解与应力分析
THANKS FOR WATCHING
感谢您的观看
SolidWorks Simulation
基于SolidWorks平台的有限元分析插件, 适合中小型企业使用。
COMSOL Multiphysics
提供多物理场耦合分析的有限元软件,适用 于多物理场仿真。
软件操作界面与基本功能
操作界面
每个软件都有自己的操作界面,用户可以通过界面进行模型建立、网格划分、边界条件设置等操作。
对非线性问题的处理能力有限
对于非线性问题,有限元法的求解过程可能变得不稳定或收敛速度变慢。
未来发展方向与挑战
发展更高效的算法
为了进一步提高有限元法的求解 效率,需要研究和发展更高效的
算法和软件实现。
处理更复杂的问题
有限元法讲稿

平面问题的有限单元法有限单元法是随着计算机的出现而发展起来的一种有效数值计算方法,有限单元法出现于40年代,被应用于飞机结构分析,有限元这个术语是1956年首先使用的。
目前已广泛地用于工程结构的力学分析中。
第一节基本概念一、实质理想化连续体―――――――单元集合体(解析模拟、逼近求解区域)无限自由度有限个自由度有限单元法首先把结构划分成许多单元,在一定的简化假设前提下,研究单元的力学特性,即单元分析;然后把各单元综合起來, 把局部的力学特性扩展到整体, 即整体分析;12最后导出一组以结构结点位移为未知量的代数方程组。
通过求解方程组而得到单元的结点位移值,就可近似计算出结构任意一点的受力状态。
这种以结点位移为基本未知量的计算方法称为有限单元位移法。
二、理论基础弹性力学:变分原理能量原理基本方程:几何方程、物理方程1. 平面问题的几何方程⎭⎬⎫⎩⎨⎧⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧∂∂∂∂∂∂∂∂=⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧∂∂+∂∂∂∂∂∂=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=v u x y y x x v y u y v x u xy y x 00γεεε3这就是弹性力学平面问题的几何方程,它给出了某一点的位移与该点应变之间的关系。
反映了变形协调关系。
2. 平面问题的物理方程⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧xy y x xy y x E γεεμμμμτσσ2100010112 即: D εζ= 给出了力与变形之间的关系,称为平面问题物理方程,是针对平面应力问题推导出的。
对于平面应变问题,只需将公式中的E 换成21μ-E,把μ换成μμ-1即可。
这样,弹性矩阵D 就变为:4⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡-----+-=)1(22100011011)21)(1()1(μμμμμμμμE D (2.8)这就是适用于平面应变问题的弹性矩阵。
3. 弹性体的能量原理(1)应变能 在弹性范围内,对于平面问题,某个面积A 的应变能可以用下式表示tdxdy V xy xy y y Ax x )(21γτεσεσ++=⎰⎰ 写成矩阵的形式为5t d x d y V AT ζε⎰⎰=21 (2)外力势能 外力势能用矩阵的形式可表示为∑-=i T i P V P d 式中 {}iy ixi P P =P −−作用在弹性体i 点的外力分量; {}i i i v u =d −−i 点的位移分量(3)弹性体的总势能弹性体在外力作用的总势能定义为应变能和荷载势能之和,即∑⎰⎰-=+=i T i AT PP tdxdy V V E P ζεd 21 (4)最小势能原理6单元的众多的结点位移)(e δ中, 须满足条件:0)(=e P E δ即 0)()(=∂∂e e P E δ的一组位移才是真正的位移。
有限元法PPT课件

Motorola– Drop Test Fujitsu-Computers Intel –Chip Integrity
电子
Baxter - Equipment J&J – Stents Medtronic - Pacemakers
医疗
Principia-spain Arup-U.K. T.Y. Lin - Bridge
有限元法
左图所示,为分析齿轮上一个齿内的应力分布,可分析图中所示的一个平面截面内位移分布.作为近似解,可以先求出图中各三角形顶点的位移.这里的 三角形就是单元,其顶点就是节点。
从物理角度理解, 可把一个连续的齿形截面单元之间在节点处以铰链相链接,由单元组合而成的结构近似代替原连续结构,在一定的约束条件下,在给定的载荷作用下,就可以求出各节点的位移,进而求出应力.
一.Abaqus公司简介
公司
’00 ’01 ’02 ’03 ’04 ‘05 ’06 ‘07
18%
18%
20%
SIMULIA公司(原ABAQUS公司)成立于1978年,全球超过600名员工,100% 专注于有限元分析领域。 全球28个办事处和9个代表处 业务迅速稳定增长,是当前有限元软件行业中唯一保持两位数增长率的公司。 2005年5月ABAQUS加入DS集团,将共同成为全球PLM的领导者
Where :
Displacement interpolation functions (位移插值函数)
13.3 Approximating Functions for Two-Dimensional Linear Triangular Elements (二维线性三角形单元的近似函数)
node (节点)
element(单元)
有限元ppt课件

因此有 y(x) (x)
试探函数中所取的项数越多,逼近的精度越高。
将试探函数代入式(1-9),可以得到关于n个待定系数
的泛函表达式,简记为 I y(x) I(1,2,3, ,n)
根据多元函数有极值的必要条件,有
1
I (1,2 ,3,
2
I (1,2 ,3,
力,它反映了内力在截面上的分布密度。
z
y
o
zx
xz
z zy
yz
切应力互等定律 xy yx , xz zx , yz zy
y
应力矩阵
x xy
yx
T
x y z xy yz zx
y
x
z
微分体的应力分量
v y w z u v
0
0
yz
zx
y x y
v
w
0
y
0
x
0
z
u v
0
w
39
厚度为1的微分体,在水平方向拉
力F的作用下发生了位移 xdx
拉力表达式:
F xdy 1
x
x dy
拉力做的功:
dx
xdx
dW
1 2
F xdx
将F代入:
dW
1 2
x
x
dxdy
40
储存在微分体内的应变能:
航空航天结构中的有限元方法编程作业资料

航空航天结构中的有限元方法编程大作业姓名熊蕾班级110514学号11051121目录一、使用程序介绍 (1)二、数据文件生成程序 (1)三、计算程序 (2)四、例题求解 (4)例1: (4)例2: (5)例3: (6)五、程序不足之处与改进方法 (7)1、数据文件的可读性 (7)2、数据长度混乱 (7)3、输入量很难从图形得出 (8)4、算法的性质 (8)六、参考文献 (8)一、使用程序介绍由于有限元刚度方程求解过程是一个矩阵运算的过程,使用matlab编写程序会比C语言更加方便。
MATLAB是近年来得到迅速发展的数学软件,它将高性能的数值计算和可视化集成在一起,并提供了大量内置函数,被广泛地应用于科学计算、控制系统、信息处理等领域的分析、仿真和实际工作。
二、数据文件生成程序为了方便计算程序读取数据,我编写了以下程序来生成txt数据文件,这样所得的数据文件中数据储存的格式正是计算时需要的格式。
format longnn = input('请输入节点总数\n')%节点数nm = input('请输入单元数\n')%单元数L = input('请输入各杆件的长度,用分号隔开,单位为:m\n')%各单元长度(单位:米)E = input('请输入弹性模量,单位为:GPa\n')%弹性模量A = input('请输入杆件截面积,单位为:m^2\n')%截面积Q = input('请输入各杆倾角,用分号隔开,单位为:度\n')%各杆倾角JM = input('请输入各杆对应的节点(同一杆的节点序号用空格隔开,不同杆之间用分号隔开)\n')%各杆对应的节点(同一杆的节点序号用空格隔开,不同杆之间用分号隔开)nf=4%外载荷数量F=input('请输入外载荷,单位为:kN\n')%外载荷矩阵bc=input('请输入边界条件,被约束自由度序号的元素置1,未被约束置0')%边界条件,有位移约束的置1data=zeros(11,2*nn)data(1,1)=nndata(2,1)=nmfor i=1:nmdata(3,i)=L(i,1)enddata(4,1)=Edata(5,1)=Afor i=1:nmdata(6,i)=Q(i,1)endfor i=1:nmdata(7,i)=JM(i,1)endfor i=1:nmdata(8,i)=JM(i,2)enddata(9,1)=nffor i=1:2*nndata(10,i)=F(1,i)endfor i=1:2*nndata(11,i)=bc(1,i)endfilename=input('请输入保存的文件名\n','s')fid=fopen(filename,'w');fprintf(fid,'%f %f %f %f %f %f %f %f %f %f %f %f\n',data)fclose(fid)运行这个.m文件,在指引下输入各项数据,最后能在工作目录下生成文件名可以自定义的txt数据文件。
有限元方法概述

主要工学硕士数学课程
工程数学 计算方法(数值分析) 随机过程 矩阵论 运筹学(最优化方法) 图论 模糊数学 有限元方法 小波分析 应用泛函分析北 Nhomakorabea航空航天大学
数学课程在研究生培养中的重要性
科技发展日新月异,数学科学地位不断提
高,在自然科学和工程技术方面广泛应用。 数学的面貌发生很大变化,现代数学在理 论上更加抽象、方法上更加综合、应用上 更加广泛。 综合运用数学的能力关系到研究生的创新 能力和研究水平的提高,对研究生的论文 质量至关重要。
X
北京航空航天大学
(2)单元分析 用单元节点位移表示单元内部位移-第i个单元 中的位移用所包含的结点位移来表示。
ui 1 ui ( x xi ) u ( x ) ui Li ui 第i结点的位移 xi 第i结点的坐标
北京航空航天大学
第i个单元的应变 应力 内力
du ui 1 ui i dx Li
自重作用下等截面直杆的解
受自重作用的等截面直杆 如图所示,杆的长度为L, 截面积为A,弹性模量为 E,单位长度的重量为q, 杆的内力为N。 试求:杆的位移分布,杆 的应变和应力。
北京航空航天大学
材料力学解答
N ( x) q( L x)
x
N ( x) q ( L x) A A
d2y EI 2 P ( x L) dx
M ( x) EI d2y dx 2
x
和边界条件
y |x 0 0 dy |x 0 0 dx
M ( x) P ( x L)
北京航空航天大学
再如对于弹性力学问题,可以建立起基本方程与 边界条件,如下: 平衡方程: 几何方程: 物理方程: 边界条件:
航空航天结构中有限元方法理论与实践的教学衔接

龙源期刊网
航空航天结构中有限元方法理论与实践的教学衔接
作者:王晓军邱志平吕峥
来源:《科教导刊》2012年第20期
摘要在航空航天结构分析中,有限元方法应用极其广泛。
本文针对航空航天背景专业,提出有限元方法课程教学应紧密围绕强化理论基础、突出工程实践、重在人才素质、面向科技创新的改革思路,结合作者多年从事航空航天有限元教学的经验,给出了在航空航天结构有限元方法教学中应涉及的几个实例。
关键词航空航天结构有限元方法实践教学。
FOD培训资料

基于激光雷达的fod检测技术
基于毫米波雷达的fod检测技术
基于超声波的fod检测技术
fod检测技术的优势与局限
优势
高精度、高灵敏度、高可靠性、非接触式测量等。
局限
受环境影响较大,如光线、空气质量等;对物体特征的提取和识别存在一定 的难度;设备成本较高,难以大规模应用。
03
fod应用场景与案例
fod在安检领域的应用
FOD技术可用于快递物流行业,检测包裹中是否含有违禁品和
危险品。
冷链物流
02
FOD技术可用于冷链物流行业,检测冷藏车和冷藏物品中是否
有异常发热和损坏等情况。
供应链管理
03
FOD技术可用于供应链管理,检测原材料、零部件等物品中是否含Fra bibliotek金属杂质等违禁品。
fod在航空领域的应用
机场航班
FOD技术可用于机场航班安检,检测乘客随身携带物品和托运行李中是否含有违禁品和危 险品。
《fod培训资料》
xx年xx月xx日
目录
• fod基本概述 • fod检测原理与技术 • fod应用场景与案例 • fod发展趋势与展望 • fod常见问题与对策
01
fod基本概述
fod是什么
fod(有限元法)是一种广泛应用于工程分析的数值计算方法 ,它通过将连续体离散化为有限个单元,并对每个单元进行 数值分析,以实现对整个系统的近似求解。
fod将复杂的物理系统抽象化为数学模型,通过计算机程序实 现对模型的分析、预测和优化。
fod的历史与发展
fod起源于20世纪50年代,最早应用于航空航天领域,后 来逐渐扩展到土木工程、机械制造、能源、交通等领域。
经过几十年的发展,fod已成为工程分析的常用工具,并 不断得到优化和完善。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 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扭转问题。