有限元习题册

有限元法基础及应用习题集一、填空1. 有限元法是求解连续场力学和物理问题的一种方法。用有限元法求解连续体或结构的力学问题的三个主要步骤是：; O2. 离散化就是把连续体或结构分割成若干个在处相互连接，尺寸有限的结合体來代替原來的连续结构。3. 单元分析阶段导出的单元刚度方程建立了和之间的关系。单元刚度方程的核心是矩阵。该矩阵具有性和性，且主对角元素O4. 建立实体单元(一维杆单元、三节点三角形平面单元等)的刚度方程时，须应用作为平衡条件。5. 弹性力学几何方程反映弹性体变形时和之间的关系。6. 单元位移模式：=旳5中网称为矩阵。该方程的含义是O7. 单元某节点i的形函数叫在该点的值为,在其它节点的值均为。一个单元所有节点形函数之和等于O8. 作用在单元上的载荷须按的原则移置到节点上，因为。9. 单元刚度矩阵奇异性的力学意义是：。10. 结构有限元平衡方程K5=0建立了有限元离散结构中节点的和之间的关系。该方程的力学意义是有限元离散结构中节点的和之间的平衡。11. 整体刚度矩阵具有如下性质：(3) o12. 对一定的有限元网格，整体刚度矩阵的半带宽与有关。半带宽越小，求解时占用计算机资源o13. 为保证有限元解的收敛性，单元位移模式应满足和。14. 建立任意形状和方位平面四边形单元和空间六面体单元时，需要釆用与单元位移模式中相同的用局部坐标表示的节点形函数对节点坐标进行插值以获得一种坐标变换，这种变换称为，釆用等参变换的单元称为O15. 节点数越多的单元，其位移模式多项式,单元的能力越强，所以精度。16. 弹性力学几何方程反映弹性体变形时和之间的关系。17. 弹性力学边界条件包括和o18. 弹性体的虚位移是假想在弹性体上发生的满足条件的微小位移场。弹性体的虚功原理可以概括为等于o19. 弹性力学物理方程反映弹性体变形时和之间的关系。20. 平面应力问题的典型例子是、平面应变问题的典型例子是。21. 建立平面问题或空间问题的单元特性方程（单元分析）阶段，需要用到弹性力学的方程和方程。二、简答题1. 简述弹性力学平面问题有限元法中单元特性分析的过程。2. 简述建立整体有限元平衡方程的过程。3. 平面三节点三角形单元中位移、应变和应力具有什么特征？有何优缺点？4. 四节点矩形单元中位移、应变和应力具有什么特征？有何优缺点？5. 简单三角形单元刚度矩阵元素的大小与哪些因素有关？与哪些因素无关？6. 画出三节点三角形单元形函数的图形，并分析其在边界上的分布特点。7. 对一个给定的弹性力学问题，有那些途径可以提高有限元法求解精度？8. 按位移求解的有限元法中：（1）应用了哪些弹性力学的基本方程？（2）应力边界条件及位移边界条件是如何反映的？（3）力的平衡条件是如何满足的？（4）变形协调条件是如何满足的？9. 有限元的收敛条件是什么？证明三节点三角形单元满足收敛条件。10. 平面应力三角形单元和空间轴对称三角形单元分别代表物理空间中什么样的物体？11. 试述所学各类单元节点数、节点位移分量、单元自由度数目。12. 位移函数应满足哪些要求？写出梁单元的位移函数。13. 空间轴对称问题的位移分量、应变分量、应力分量有哪些？14. 简单（纯弯）梁单元的节点位移分量、单元自由度？15. 平面梁单元的节点有几个自由度？其在局部坐标系下节点位移分量有哪些？16. 弹性力学的基本假设？弹性力学有哪些基本方程和边界条件？17. 一维杆单元、三节点三角形平面单元、三节点三角形空间轴对称单元的形函数矩阵、应变矩阵、单元刚度矩阵的行数和列数分别是多少？18. 对于平面问题简单三角形单元，为什么单元刚度矩阵是常数矩阵？19. 什么是等参变换？等参变换的基本条件是什么？哪些情况使等参变换不成立？划分等参单元时应注意哪些问题？20. 应用等参单元时，为什么要釆用高斯积分？高斯积分的数目如何确定？21. 弹性力学平面问题求解时应用的三角形单元是等参单元吗？为什么？21. 什么是等参单元，等参单元的主要优点是什么？22. 写出平面四节点等参元的坐标变换的雅克比(Jacobian)矩阵。23. 非节点载荷为什么要等效变换成节点载荷，如何变换？作变换时应注意什么问题？24. 结构原始平衡方程式为什么要做约束处理？25. 试述平面应力问题和平面应变问题的儿何、受力和变形特征。26. 平面应力问题和平面应变问题有什么区别？27. 举例说明，在什么样情况下可以将工程问题转化成平面应力问题？在什么情况下可以将工程问题转化为平面应变问题？28. 为什么说平面三节点三角形单元为常应力单元，如何解决由于这种单元的特点所引起的计算精度不高的问题？29. 用示意图画出空间结构常用的单元类型。30. 简单四面体单元为什么说是一种常应变单元？31. 轴对称结构有什么特点？轴对称结构如何简化处理？三、计算与分析1. 如图所示，根据弹簧单元的刚度方程推导出系统的平衡方程。X.k/比2-123Fu2F?i(3t凡2. 根据弹赞单元的刚度方程，导出下列系统的整体刚度平衡方程。并代入边界条件，得出节点位移求解方程，并得出节点3的位移和节点1的支反力。k、=100N/mm,仏=200N/mm,k、=100N/mmP=500N,u,=uA=03. 对图示弹赞系统，kl=300N/mm,k2=k3=200N/mm,k4=200N/mm,Fl=600N,F2=400No求：（1）其总刚度矩阵；（2）节点1、2、3的位移；（3）节点4、5的反力；（4）弹赞1、2、3、4中的力。4. 如下图所示，5个弹赞连接在一起，各弹赞的刚度系数如图上标出。求：（1）系统刚度矩阵；（2）节点3处作用F力后，各节点的位移5,固定节点1、6处的反作用力。5. 如图所示一维杆系由两个材料相同截面不同的直杆单元（1）与（2）组成，弹性模量E。在节点1、2、3上作用有轴向集中载荷0、Q：、Q,而平衡。试求解下列各问题：（1）建立结构的有限元平衡方程；（2）如果节点1被固定（u产0）,Q：=P,QfO,通过建立的平衡方程求各节点位移、节点1约朿反力。（3）如果QfO,Q尸P,其他条件不变，试根据问题（2）的解答和有关力学概念直接给出节点2、3的位移。q一分211(1)6. 图示杆-弹赞系统，材料弹性模量为E。试列出其有限元平衡方程，并进行约束处理。57. 如图所示一维杆系由两个材料相同截面不同的直杆单元（1）与（2）组成，弹性模量E,节点1,3固定，节点2受集中力P。试求解下列各问题：（4）建立结构的有限元平衡方程。（5）求解节点2的位移和各杆的应力。（6）如果P=0,且所有杆上受沿x方向作用的均匀线分布力q,求未知节点位移和固定端反力。2/1,A,EIP8. 平面桁架由2根相同的杆组成（E,A,L）o求:（1）节点2位移；（2）每根杆应力。9. 如图所示三杆钢桁架，节点1、节点3处固定，节点2处受力耳“所有杆件材料相同，弹性模量为E,截面积均为A,求各杆受力。10. 如图所示2杆结构，每根杆的弹性模量均为E,横截面积均为A。建立坐标系和节点系统如图所示，在节点1处作用x方向的力F,求,v1Oii. 证明杆单元变换矩阵rT=Tfo12. 如图所示刚架由两根等截面工字型钢构成，两端固支，系统所受载荷如图所示。梁截面积=0.006m2,截面惯性距/=6.87xW5m4,弹性模量E=2.1xlOskN/nr,每根梁长L=5mo求：每根梁所示内力。Q=2kN13试推导梁单元的坐标变换矩阵20000A0000001000000A00000000001其中，A=cos(p,/=sin(po14. 如图所示刚架结构，所有梁材料和截面尺寸相同，截面积为丛惯性距为久材料弹性模量为&试写出每根梁单元的刚度矩阵和结构的总体刚度矩阵。15. 有一正方形板，沿对角承受压力作用，板厚r=lm,载荷P=20kN,如图所示。为简9化计算，设泊松比“=0,材料弹性模量为求它的应力分布。S16. 试证明三角形单元形状函数M（兀y）=土仏+bjX+c;j）i=1,2,32A满足下列性质：N（x,y）=Q，lj17. 如下图（a）所示悬臂深梁，自由端有垂向均布载荷尸，梁的厚度为&设材料弹性模量为E,泊松比“=1/3,若采用（b）所示的简单网格系统，求各节点的位移。/iEL/fJ/2m18. 正方形板如图所示，边长为厚度为,弹性模量为泊松比0.15,节点1作用集中力尸，节点2、3、4被固定，若采用图示坐标系统和单元节点结构，求各节点位移和应力。S19. 如图所示，用近似法取r=r,z=Z计算两个轴对称单元a、b的单元刚度矩阵斤7,设材料的弹性模量为E,泊松比0.15o20. 试写出如图所示5节点等参元形状函数，并求出其雅克比矩阵的表达式和单元刚度矩阵。35221. 图示悬臂梁为平面应力问题，试写出边界条件。922. 如图平面问题，以单元为例，通过实算，讨论单元点号按顺序轮换时单元刚度矩阵K及其变化规律。23. 如图所示平面三角形桁架，终点坐标为：1（0,0）,2（L&,L/0）,3（忑厶，0）,E、A为弹性模量及截面积。用有限元法求：（1）节点位移；（2）单元内力；（3）支座反力。24. 平面桁架如图所示，；=2xlO5kg/cnr,A=l.Ocnr求节点位移和单元内力，并利用节点1的平衡检验计算结果。(1) 节点1位移为0,Q2=Qs=P(2) 节点1位移为0,Q:=Q3=0,整个杆受到沿轴线的均匀线分布力q,方向向右。27. 根据材料力学知识和单元刚度矩阵物理意义推导出简单梁单元刚度矩阵的第三列和第四列元素。28. 对图示有限元模型，用符号“标出总刚度矩阵中非零子块的分布，并计算半带宽。1135681029. 对图示平面问题，考虑到对称性，试用图形表示出其有限元模型，要求:（1）划分单元，单元数目适当；（2）给出节点编号方案；（3）标出节点载荷和位移约束。30. 对图示平面问题，考虑到对称性，试用图形表示出其有限元模型，要求:（1）划分单元，单元数目适当；（2）给出节点编号方案；（3）标出节点载荷和位移约束。31根据单元刚度矩阵元素的物理意义求弹赞单元和杆单元的刚度矩阵。32. 通过对节点位移插值建立三节点三角形单元的位移模式和形函数。33. 用虚功原理推导出三节点三角形单元刚度方程。34. 对三节点三角形单元证明其形函数满足：M+N”,+N”=l35. 图示三角形单元：按公式求形函数和形函数矩阵；求该单元的应变矩阵。36. 计算图示平面三角形单元的等效节点载荷列阵。设单元厚度为h。37. 如图所示，两个形状相似的三节点三角形平面单元，对应边长比为2:1,材料、厚度相同，方位相同。约束左边上2个节点x,y方向位移，自由节点N，N,均受铅直向下集中力P。两个模型分别用有限元软件计算后，发现计算结果有下列关系：1）节点M和节点N?的位移相等；2）单元的应力是单元应力的二分之一。试对上述现象进行解释。38. 计算图示平面三角形单元的等效节点载荷列阵。设单元厚度为h。39. 将图示水坝作为平面应变问题，试用图形表示出你的有限元模型，要求:（1）用三角形单元离散，建议单元边长lm左右或小于lm（2）给出节点编号方案（3）写出节点载荷和位移边界条件2mQoh12m40. 函数“（Q如图所示，求其在和心之间有效的一维线性插值多项式。41. 如图所示正方形桁架，周边长桁架由五条杆单元组成，弹性模量为截面积为丛求尸载荷作用下2、3点的位移。15#42. 如图所示二节点杆单元ij,沿杆轴线分别作用一均布载荷绻(如图(a)和分布载荷(如图(b)o分别求两种情况下的等效节点载荷。#43. 采用杆单元的方法，求解如图所示结构的所有节点的位移、三个杆单元的应力、支座反力。相关的材料参量和尺寸为Ej=E2=E3=2x105Pa,3At=2A2=A3=0.06m2,=I：=I3=#X/1400nun45. 如图所示的等剖而梁，弯曲刚度为EI,承受分布载荷q及q(x)作用，求:(1) 各梁的等效结构载荷;(2) 节点位移；(3) 单元的节点力。46. 一悬臂梁的一端由弹费支持，弹赞的刚度系数为在载荷尸作用下，求端点2的位移及转角。1ElVL3#47. 如图所示桁架结构，各元件的A,Z均相同，1-4杆做短了小试求(1) 节点位移；(2) 1-4杆应力。48. 利用对称条件，处理以下结构(要求画出简化图以及给出其边界条件)。49-如图所示杆结构，杆剖面面积为儿材料的弹性模量都为求2点作用载荷为P时节点的位移。50. 如图所示的自由体结构，在平衡力系作用下，用有限元分析问题时边界条件如何处理?51. 如图所示杆板结构，按下列情况划分，选取单元：（1）结构由10个两节点杆单元和8个三节点三角形板单元集合而成;（2）结构由5个节点杆单元和2个六节点三角形板单元集合而成试分析：两种分单元情况下，采用相同的节点编号，（1）总刚度矩阵大小是否相同？（2）半带宽是否一样？（3）杆板单元间位是否协调？（4）单元中内力特点是否一样？52. 试求如图所示结构的节点位移。己知：各杆元E,A,L均相同。53. 试求如图所示梁结构的节点位移。1EJ2EJ3LLr54.如果对图所示的有限元网格给出如下的单元刚度矩阵，试求引入边界条件Ui=0jVi=04=H至15以后的总刚度矩阵。重新编节点号是否能减小带宽?311113111131111319#55.正方形平板，厚度为卞边长为弹性模量鸟材料泊松比载荷E按图所示分元，求b3点的位移。156. 推导三节点平板单元在局部坐标系xOy中的单元刚度矩阵。#57. 三节点三角形单元ijm的位移函数能否选为(1)w(x9y)=珂+a2x2+a5yv(x,y)=a4+a5x+a6y2(2)u(x.y)=atx2+a2xy+ay2v(x9y)=a4x2+a5xy+a6y2#58求三节点和六节点三角形单元在自重作用下的等效节点载荷。（设比重为厂厚度为t,单元面积）。59. 证明常应变三角形单元发生刚体位移时，单元中将不产生应力。60. 以平面问题常应变三角形单元为例，证明单元刚度矩阵的任何一行（或列）元素的总和为零。61. 三角形单元ijm的jm边作用有如图所示线性分布载荷，求节点载荷向量。#XO62. 试写出图示平面梁结构的单元节点位移列阵。77777II/7/77777X77777#63. 图示桁架，各杆的拉压刚度EA=105kN,斜杆在制作时比设计尺寸长了5=OOlm,求由此产生的内力。64. 图示一个平面应力状态的直角三角形单元，设E为常量，“=丄，r=i,试求:6（1）形函数矩阵庄（2）应力矩阵S；（3）单元刚度矩阵W21