西工大结构有限元习题库

有限元法基础及应用习题集一、填空1.有限元法是求解连续场力学和物理问题的一种方法。用有限元法求解连续体或结构的力学问题的三个主要步骤是：①；②；③。2.离散化就是把连续体或结构分割成若干个在处相互连接，尺寸有限的结合体来代替原来的连续结构。3.单元分析阶段导出的单元刚度方程建立了和之间的关系。单元刚度方程的核心是矩阵。该矩阵具有性和性，且主对角元素。4.建立实体单元（一维杆单元、三节点三角形平面单元等）的刚度方程时，须应用作为平衡条件。5.弹性力学几何方程反映弹性体变形时和之间的关系。6.单元位移模式eNvu中N称为矩阵。该方程的含义是。7.单元某节点i的形函数Ni在该点的值为，在其它节点的值均为。一个单元所有节点形函数之和等于。8.作用在单元上的载荷须按的原则移置到节点上，因为。9.单元刚度矩阵奇异性的力学意义是：。10.结构有限元平衡方程QK建立了有限元离散结构中节点的和之间的关系。该方程的力学意义是有限元离散结构中节点的和之间的平衡。11.整体刚度矩阵具有如下性质：①②③④。12.对一定的有限元网格，整体刚度矩阵的半带宽与有关。半带宽越小，求解时占用计算机资源。13.为保证有限元解的收敛性，单元位移模式应满足和。14.建立任意形状和方位平面四边形单元和空间六面体单元时，需要采用与单元位移模式中相同的用局部坐标表示的节点形函数对节点坐标进行插值以获得一种坐标变换，这种变换称为，采用等参变换的单元称为。15.节点数越多的单元，其位移模式多项式，单元的能力越强，所以精度。16.弹性力学几何方程反映弹性体变形时和之间的关系。17.弹性力学边界条件包括和。18.弹性体的虚位移是假想在弹性体上发生的满足条件的微小位移场。弹性体的虚功原理可以概括为等于。219.弹性力学物理方程反映弹性体变形时和之间的关系。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．结构原始平衡方程式为什么要做约束处理？325．试述平面应力问题和平面应变问题的几何、受力和变形特征。26．平面应力问题和平面应变问题有什么区别？27．举例说明，在什么样情况下可以将工程问题转化成平面应力问题？在什么情况下可以将工程问题转化为平面应变问题？28．为什么说平面三节点三角形单元为常应力单元，如何解决由于这种单元的特点所引起的计算精度不高的问题？29．用示意图画出空间结构常用的单元类型。30．简单四面体单元为什么说是一种常应变单元？31．轴对称结构有什么特点？轴对称结构如何简化处理？三、计算与分析1．如图所示，根据弹簧单元的刚度方程推导出系统的平衡方程。2．根据弹簧单元的刚度方程，导出下列系统的整体刚度平衡方程。并代入边界条件，得出节点位移求解方程，并得出节点3的位移和节点1的支反力。3.对图示弹簧系统，k1=300N/mm，k2=k3=200N/mm，k4=200N/mm，F1=600N，F2=400N。求：（1）其总刚度矩阵；（2）节点1、2、3的位移；（3）节点4、5的反力；（4）弹簧1、2、3、4中的力。44．如下图所示，5个弹簧连接在一起，各弹簧的刚度系数如图上标出。123456kk2k3kk2F求：（1）系统刚度矩阵；（2）节点3处作用F力后，各节点的位移，固定节点1、6处的反作用力。5．如图所示一维杆系由两个材料相同截面不同的直杆单元（1）与（2）组成，弹性模量E。在节点1、2、3上作用有轴向集中载荷Q1、Q2、Q3而平衡。试求解下列各问题：（1）建立结构的有限元平衡方程；（2）如果节点1被固定（u1=0），Q2=P，Q3=0，通过建立的平衡方程求各节点位移、节点1约束反力。（3）如果Q2=0，Q3=P，其他条件不变，试根据问题（2）的解答和有关力学概念直接给出节点2、3的位移。6.图示杆-弹簧系统，材料弹性模量为E。试列出其有限元平衡方程，并进行约束处理。57.如图所示一维杆系由两个材料相同截面不同的直杆单元（1）与（2）组成，弹性模量E，节点1，3固定，节点2受集中力P。试求解下列各问题：（4）建立结构的有限元平衡方程。（5）求解节点2的位移和各杆的应力。（6）如果P=0，且所有杆上受沿x方向作用的均匀线分布力q，求未知节点位移和固定端反力。8.平面桁架由2根相同的杆组成（E，A，L）。求：（1）节点2位移；（2）每根杆应力。9．如图所示三杆钢桁架，节点1、节点3处固定，节点2处受力2xF，2yF，所有杆件材料相同，弹性模量为E，截面积均为A，求各杆受力。1232xF2yFxy45123L610．如图所示2杆结构，每根杆的弹性模量均为E，横截面积均为A。建立坐标系和节点系统如图所示，在节点1处作用x方向的力F，求1u，1v。FLL12312xy11．证明杆单元变换矩阵T1TT。12．如图所示刚架由两根等截面工字型钢构成，两端固支，系统所受载荷如图所示。梁截面积A=0.006m2，截面惯性距546.8710mI，弹性模量822.110kN/mE，每根梁长5mL。求：每根梁所示内力。5m2.5m2.5mkN23QkN/m5.2qkN12QmkN05.0MkN41Q13．试推导梁单元的坐标变换矩阵0000000000100000000000000001T其中，cos，sin。14．如图所示刚架结构，所有梁材料和截面尺寸相同，截面积为A，惯性距为I，材料弹性模量为E。试7写出每根梁单元的刚度矩阵和结构的总体刚度矩阵。LL12312xy315．有一正方形板，沿对角承受压力作用，板厚1mt，载荷20kNP，如图所示。为简化计算，设泊松比0，材料弹性模量为E，求它的应力分布。kN20P2m2mkN20P16．试证明三角形单元形状函数1,1,2,32iiiiNxyabxcyi满足下列性质：0,,1,iijNxyij17．如下图（a）所示悬臂深梁，自由端有垂向均布载荷F，梁的厚度为t，设材料弹性模量为E，泊松比13，若采用（b）所示的简单网格系统，求各节点的位移。81m2mF)(a2/F)(b2/F13241218．正方形板如图所示，边长为a，厚度为t，弹性模量为E，泊松比0.15，节点1作用集中力F，节点2、3、4被固定，若采用图示坐标系统和单元节点结构，求各节点位移和应力。xyF13241219．如图所示，用近似法取,rrzz计算两个轴对称单元a、b的单元刚度矩阵eK，设材料的弹性模量为E，泊松比0.15。LLZijmaLLijmb20．试写出如图所示5节点等参元形状函数，并求出其雅克比矩阵的表达式和单元刚度矩阵。1234521．图示悬臂梁为平面应力问题，试写出边界条件。9ACBDxOyq22．如图平面问题，以单元④为例，通过实算，讨论单元点号按顺序轮换时单元刚度矩阵K④及其变化规律。12341234567xyaaa23．如图所示平面三角形桁架，终点坐标为：1（0,0），2（/2L，/2L），3（2L，0），E、A为弹性模量及截面积。用有限元法求：（1）节点位移；（2）单元内力；（3）支座反力。123yxP24．平面桁架如图所示，62210kg/cmE，21.0cmA。求节点位移和单元内力，并利用节点1的10平衡检验计算结果。123yxkg100P40cm30cm25．下图中结构分别采用（b）、（c）两种编节点号方式，分别求其刚度矩阵带宽。123456789123456123456)(a)(b)(c26．教材P20练习题1-9中，求下列2种情况下节点位移、节点1约束反力。（1）节点1位移为0，Q2=Q3=P（2）节点1位移为0，Q2=Q3=0，整个杆受到沿轴线的均匀线分布力q，方向向右。27．根据材料力学知识和单元刚度矩阵物理意义推导出简单梁单元刚度矩阵的第三列和第四列元素。4321444342413433323124232221141312114321uuuuaaaaaaaaaaaaaaaassss28．对图示有限元模型，用符号“△”标出总刚度矩阵中非零子块的分布，并计算半带宽。1129．对图示平面问题，考虑到对称性，试用图形表示出其有限元模型，要求：（1）划分单元，单元数目适当；（2）给出节点编号方案；（3）标出节点载荷和位移约束。30.对图示平面问题，考虑到对称性，试用图形表示出其有限元模型，要求：（1）划分单元，单元数目适当；（2）给出节点编号方案；（3）标出节点载荷和位移约束。31．根据单元刚度矩阵元素的物理意义求弹簧单元和杆单元的刚度矩阵。32.通过对节点位移插值建立三节点三角形单元的位移模式和形函数。33.用虚功原理推导出三节点三角形单元刚度方程。34.对三节点三角形单元证明其形函数满足：1nmlNNN1235．图示三角形单元：①按公式求形函数和形函数矩阵；②求该单元的应变矩阵。ijmxyab36．计算图示平面三角形单元的等效节点载荷列阵。设单元厚度为h。ijmxyaa0q02q37．如图所示，两个形状相似的三节点三角形平面单元，对应边长比为2:1，材料、厚度相同，方位相同。约束左边上2个节点x，y方向位移，自由节点N1，N2均受铅直向下集中力P。两个模型分别用有限元软件计算后，发现计算结果有下列关系：1）节点N1和节点N2的位移相等；2）单元①的应力是单元②应力的二分之一。试对上述现象进行解释。单单单N1Pxyz单单单N2P38.计算图示平面三角形单元的等效节点载荷列阵。设单元厚度为h。13xyijmaa0q39.将图示水坝作为平面应变问题，试用图形表示出你的有限元模型，要求：（1）用三角形单元离散，建议单元边长1m左右或小于1m（2）