航空航天结构中的有限元方法专题培训(pdf 110页)
- 格式:pdf
- 大小:1.66 MB
- 文档页数:110
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年。
平面问题的有限单元法有限单元法是随着计算机的出现而发展起来的一种有效数值计算方法,有限单元法出现于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 δ的一组位移才是真正的位移。
航空航天结构中的有限元方法编程大作业姓名熊蕾班级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数据文件。
龙源期刊网
航空航天结构中有限元方法理论与实践的教学衔接
作者:王晓军邱志平吕峥
来源:《科教导刊》2012年第20期
摘要在航空航天结构分析中,有限元方法应用极其广泛。
本文针对航空航天背景专业,提出有限元方法课程教学应紧密围绕强化理论基础、突出工程实践、重在人才素质、面向科技创新的改革思路,结合作者多年从事航空航天有限元教学的经验,给出了在航空航天结构有限元方法教学中应涉及的几个实例。
关键词航空航天结构有限元方法实践教学。
航空航天结构中的有限元方法邱志平 王晓军 编著北京航空航天大学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扭转问题。