西北工业大学研究生入学试题
考试科目：结构有限元素法题号：05602 说明：共 3 页第1 页
共 3 页第2 页题号：05602
2011年有限元试题 一、简答 1、有限元平面问题有几个自由度？空间问题有几个自由度？弹性力学三大方程分别应用了何种假定？ 2、位移模式的完备性条件是什么？必须要满足吗？ 3、如图所示节点编号，计算半带宽，并画出最合理的节点编号，刚度矩阵中对角元素为何都为正值？ 4、等参单元变换的必要条件是什么？ 5、简述铁木辛柯梁的基本特点，解释何谓剪切锁死现象，并给出避免出现剪切闭锁的条件。 二、推导题 1、用自然坐标推导四节点矩形单元的形函数 2、计算等效荷载，如图 3、分析下图的高斯点选取是否合理 三、计算 1、如图 （1）计算三角形单元的位移模式 （2）已知三节点的位移分别为() ,;,;,j j m i m i u u u v v v （已知数据），A 点的坐标（0.5，0.2,0）求A 点的位移分量
2、已知基本方程d 0d u u x +=，边界条件01x u ==，选取试函数1u cx =+ 用伽辽金法确定系数c 3、求雅克比矩阵J （同作业题） 4、如图所示结构当节点编号如右图所示时，可求出其单元刚度矩阵为 11121621 2226616266a a a a a a K a a a ??????=?????? （考试中均为具体数据） 试确定局部刚度矩阵[]22K 及[]45K
2005年有限元试卷 一、简答题 1、加权残数法、变分法与有限元法的联系与区别是什么？有限元方法有什么优点？ 2、单元分析的重点是什么？单元刚度矩阵有什么特征？其中每一个元素的物理意义是什么？单元分析中坐标转换的作用是什么？ 3、在有限元法中，等参单元的主要有点事什么？在应用等参单元时，坐标变换的精度与位移模式的精度是否一样？等参单元计算中积分阶次选择的原则是什么？ 4、为什么位移型有限元应力解的精度要低于位移解？应力解的近似性表现在哪些方面？应力近似解的性质是什么？试分析下列平面单元中的位移和应力的误差量级？ （1）三点三角形单元（2）四节点举行单元（3）六节点三角形单元 （4）四节点直边四边形等参单元（5）把节点曲边四边形等参单元 5、在博班弯曲问题中，何谓协调单元和非协调单元？试论证四节点矩形板12个自由度弯曲单元是完备的非协调单元 6、常用的构造单元插值函数的方法有哪三种？各有哪些优点？ 7、什么是材料非线性问题和几何非线性问题？各自的难点是什么？材料非线性和几何非线性有什么联系和区别？ 8、什么是非线性有限元方程的迭代解法和增量解法？各自的适用范围是什么？常见的迭代收敛准则有哪些？当有限元求解结果不对时，可从哪些方面进行检查？ 9、简要说明学习“有限元的”的心得体会？ 二、计算推导题 1、图示等刚度悬臂梁受均布荷载作用，假定试函数为：()21cos x x c l f π??=- ??? ，梁截面抗弯刚度为EI ，试采用最小势能原理计算梁的最大挠度 X
一．是非题（认为该题正确，在括号中打；该题错误，在括号中打×。）(每小题2分) （1）用加权余量法求解微分方程，其权函数V和场函数u的选择没有任何限制。（×）（2）四结点四边形等参单元的位移插值函数是坐标x、y的一次函数。（√）（3）在三角形单元中，其面积坐标的值与三结点三角形单元的结点形函数值相等。（√）（4）二维弹性力学问题的有限元法求解，其收敛准则要求试探位移函数C1连续。（×）（5）有限元位移法求得的应力结果通常比应变结果精度低。（×）（6）等参单元中Jacobi行列式的值不能等于零。（√）（7）在位移型有限元中，单元交界面上的应力是严格满足平衡条件的。（×）（8）四边形单元的Jacobi行列式是常数。（×）（9）利用高斯点的应力进行应力精度的改善时，可以采用与位移插值函数不同结点的形函数进行应力插值。（√）（10）一维变带宽存储通常比二维等带宽存储更节省存储量。（√） 二．单项选择题（共20分，每小题2分）C B B C B C D C C C 1 在加权余量法中，若简单地利用近似解的试探函数序列作为权函数，这类方法称为 ____C__________。 （A）配点法（B）子域法（C）伽辽金法 2 等参变换是指单元坐标变换和函数插值采用__B____的结点和______的插值函数。（A）不相同，不相同（B）相同，相同（C）相同，不相同（D）不相同，相同 3 有限元位移模式中，广义坐标的个数应与_____B______相等。 （A）单元结点个数（B）单元结点自由度数（C）场变量个数 4 采用位移元计算得到应力近似解与精确解相比较，一般______C_____。 （A）近似解总小于精确解（B）近似解总大于精确解（C）近似解在精确解上下震荡（D）没有规律 5 如果出现在泛函中场函数的最高阶导数是m阶，单元的完备性是指试探函数必须至少 是__B____完全多项式。 （A）m-1次（B）m次（C）2m-1次 6 与高斯消去法相比，高斯约当消去法将系数矩阵化成了____C_____形式，因此，不用 进行回代计算。 （A）上三角矩阵（B）下三角矩阵（C）对角矩阵 7 对称荷载在对称面上引起的________D________分量为零。 （A）对称应力（B）反对称应力（C）对称位移（D）反对称位移 8 对分析物体划分好单元后，______C____会对刚度矩阵的半带宽产生影响。 （A）单元编号（B）单元组集次序（C）结点编号 9 n个积分点的高斯积分的精度可达到__C____阶。 （A）n-1 （B）n（C）2n-1 （D）2n 10 引入位移边界条件是为了消除有限元整体刚度矩阵K的____C______。 （A）对称性（B）稀疏性（C）奇异性 三．简答题（共20分，每题5分） 1、简述有限单元法结构刚度矩阵的特点。 （1）对称性；（2）奇异性；（3）主对角元恒正；（4）稀疏性；（5）非零元素带状分布
《有限元基础》期末测试 一、结构线性静力分析 如图所示的托架，其顶面承受2 lbf in的均匀分布载荷。托架通过有孔的表面 50/ ν=，托架尺固定在墙上，托架是钢制的，弹性模量6 =?，泊松比0.3 E psi 2910 寸如图，单位为英寸。试通过ANSYS求其变形图及von Mises应力分布图。 对题目分析。进行建模，网格划分 托架网格图
施加约束后，就可以对实体进行加载求解， 托架变形图 托架变形图输出的是原型托架和施加载荷后托架变形图的对比，
虚线部分即为托架的原型，托架变形图可看出，由于载荷的作用，托架上面板明显变形了，变形最严重的就是红色部分，这是因为其离托板就远，没有任何物体与其分担载荷，故其较容易变形甚至折断。这是我们在应用托架的时候应当注意的。 节点位移图
托架von Mises 应力分布图
上面两个图为托架的应力分布图，由图可看出主要在两孔处出现应力集中，也就是说这些地方所受的应力的最大的，比较容易出现裂痕。我们在应用托架的时候，应当注意采取一些设施，以便减缓其应力集中。特别是在施加载荷时，绝对不能够超过托架所能承受的极限，否则必将导致事故的发生。 二、动力分析 如图1有一梁板结构，板的四角由四根梁固定支撑，板质量集中于中央。梁板材料相关参数为弹性模量112210/E N m =?，泊松比0.3ν=，密度 337.810/kg m ρ=?。板的厚度0.02t =，板长2000L mm =，宽1000B mm =，板的质量100M kg =。梁长1000h mm =，截面面积为42210A m -=?，惯性矩为 84210J m -=?，现在板的表面施加均匀压力载荷如图2。试研究该梁板结构的瞬 态动力响应。 图 1 图2
西安交通大学 级研究生课程考试试题 考试（查）科目：有限元方法（II ）时间 年 月 日下午 一、4 ) 4,4(),()5,5(),()2,6(),()2,2(),(4 4332211====y x ， y x ，y x ， y x 母体单元为22?的正方形，如图所示。 求：（1）单元坐标变换()(ξηξ,,, y y x x == （2）变换的Jacobi 行列式detJ 的解析表达式，并分析该变换是否存在奇异性（8分）。 二、分析以下两种单元的位移场是否具备收敛到真实解所需的各项条件。（30） （1） 13结点矩形平面应力单元 结点参数取为：)13~ 1( ,=i v u i i 位移场为： 3 132 2 123 113 102 92 83726524321xy y x y x y xy y x x y xy x y x u ααααααααααααα++++++++++++= 3 262 2 253 243 232 222 2132021918217161514xy y x y x y xy y x x y xy x y x v ααααααααααααα++++++++++++=（2） 6自由度三角形薄板弯曲单元 结点参数取为： ()3~1=i w i ()6~4=??? ????i n w i 位移场为： 2 652 4321y xy x y x w αααααα+++++= 三、13结点平面应力单元如图所示， 在计算单元刚度矩阵时取图示的9个 积分点。试分析在单元一级是否存在 出现零变形能位移模式的可能性。 ,u x 7 8 10 9 11 12 1 2 3 4 5 6
一判断题节点的位置依赖于形态而并不依赖于载荷的位置√2. 对于高压电线的铁塔那样的框架结构的模型化处理使用梁单元×3. 不能把梁单元、壳单元和实体单元混合在一起作成模型√4. 四边形的平面单元尽可能作成接近正方形形状的单元×5. 平面应变单元也好平面应力单元也好如果以单位厚来作模型化处理的话会得到一样的答案×6. 用有限元法不可以对运动的物体的结构进行静力分析√7. 一般应力变化大的地方单元尺寸要划的小才好×8. 所谓全约束只要将位移自由度约束住而不必约束转动自由度√9. 同一载荷作用下的结构所给材料的弹性模量越大则变形值越小√10一维变带宽存储通常比二维等带宽存储更节省存储量。二、填空平面应力问题与薄板弯曲问题的弹性体几何形状都是薄板但前者受力特点是平行于板面且沿厚度均布载荷作用变形发生在板面内后者受力特点是垂直于板面的力的作用板将变成有弯有扭的曲面。平面应力问题与平面应变问题都具有三个独立的应力分量三个独立的应变分量但对应的弹性体几何形状前者为薄板后者为长柱体。位移模式需反映刚体位移反映常变形满足单元边界上位移连续。单元刚度矩阵的特点有对称性奇异性还可按节点分块。轴对称问题单元形状为三角形或四边形截面的空间环形单元由于轴对称的特性任意一点变形只发生在子午面上因此可以作为二维问题处理。等参数单元指的是描述位移和描述坐标采用相同的形函数形式。等参数单元优点是可以采用高阶次位移模式能够模拟复杂几何边界方便单元刚度矩阵和等效节点载荷的积分运算。有限单元法首先求出的解是节点位移单元应力可由它求得其计算公式为。8、一个空间块体单元的节点有 3 个节点位移变形体基本变量有位移应变应力基本方程平衡方程物理方程几何方程10.实现有限元分析标准化和规范化的载体就是单元 三选择题分等参变换是指单元坐标变换和函数插值采用__B___的结点和______ 的插值函数。不相同不相同相同相同相同不相同不相同 相同2 有限元位移模式中广义坐标的个数应与_______B____相等。单元结点个数 单元结点自由度数场变量个数 3 如果出现在泛函中场函数的最高阶导数是m阶单元的完备性是指试探函数必须至少是___B___完全多项式。-1次 次-1次 4 与高斯消去法相比高斯约当消去法将系数矩阵化成了____C_____形式因此不用进行回代计算。上三角矩阵下三角矩阵对角矩阵5 对分析物体划分好单元后会对刚度矩阵的半带宽产生影响。单元编号单元组集次序结点编号6 n个积分点的高斯积分的精度可达到__C____阶。--引入位移边界条件是为了消除有限元整 体刚度矩阵的_____C_____。对称性稀疏性奇异性三简答题 共20分每题5分、简述有限单元法结构刚度矩阵的特点。2、简述有限元法中选取单元位移函数多项式的一般原则。1、答答对前3个给4分对称性 奇异性主对角元恒正稀疏性非零元素带状分布2、答一般原则有(1) 广义坐标的个数应该与结点自由度数相等选取多项式时常数项和坐标的一次项必须完备多项式的选取应由低阶到高阶尽量选取完全多项式以提高单元的精度。有限元方法分析的目的对变形体中的位移、应力、应变进行定义和表达进而建立平衡方程、几何方程和物理方程。2)针对具有任意复杂几何形状的变形体完整得获取在复杂外力作用下它内部的准确力学信息。3)力学分析的基础上对设计对象进行强度(strength)、刚度评判修改、优化参数。有限单元法分析步骤1、结构的离散化2、选择位移模式3 、分析单元的力学特性4、集合所有单元平衡方程得到整体结构的平衡方程5、由平衡方程求解未知节点位移6、单元应变和应力的计算4连续体结构分析的基本假定连续性假设完全弹性假设均匀性假设
一．是非题（认为该题正确，在括号中打√；该题错误，在括号中打×。）(每小题2分) （1）用加权余量法求解微分方程，其权函数V和场函数u的选择没有任何限制。（）（2）四结点四边形等参单元的位移插值函数是坐标x、y的一次函数。（）（3）在三角形单元中，其面积坐标的值与三结点三角形单元的结点形函数值相等。（）（4）二维弹性力学问题的有限元法求解，其收敛准则要求试探位移函数C1连续。（）（5）有限元位移法求得的应力结果通常比应变结果精度低。（）（6）等参单元中Jacobi行列式的值不能等于零。（）（7）在位移型有限元中，单元交界面上的应力是严格满足平衡条件的。（）（8）四边形单元的Jacobi行列式是常数。（）（9）利用高斯点的应力进行应力精度的改善时，可以采用与位移插值函数不同结点的形函数进行应力插值。（）（10）一维变带宽存储通常比二维等带宽存储更节省存储量。（） 二．单项选择题（共20分，每小题2分） 1 在加权余量法中，若简单地利用近似解的试探函数序列作为权函数，这类方法称为 ________________。 （A）配点法（B）子域法（C）伽辽金法 2 等参变换是指单元坐标变换和函数插值采用______的结点和______的插值函数。 （A）不相同，不相同（B）相同，相同（C）相同，不相同（D）不相同，相同 3 有限元位移模式中，广义坐标的个数应与___________相等。 （A）单元结点个数（B）单元结点自由度数（C）场变量个数 4 采用位移元计算得到应力近似解与精确解相比较，一般___________。 （A）近似解总小于精确解（B）近似解总大于精确解（C）近似解在精确解上下震荡（D）没有规律 5 如果出现在泛函中场函数的最高阶导数是m阶，单元的完备性是指试探函数必须至少 是______完全多项式。 （A）m-1次（B）m次（C）2m-1次 6 与高斯消去法相比，高斯约当消去法将系数矩阵化成了_________形式，因此，不用进 行回代计算。 （A）上三角矩阵（B）下三角矩阵（C）对角矩阵 7 对称荷载在对称面上引起的________________分量为零。 （A）对称应力（B）反对称应力（C）对称位移（D）反对称位移 8 对分析物体划分好单元后，__________会对刚度矩阵的半带宽产生影响。 （A）单元编号（B）单元组集次序（C）结点编号 9 n个积分点的高斯积分的精度可达到______阶。 （A）n-1 （B）n（C）2n-1 （D）2n 10 引入位移边界条件是为了消除有限元整体刚度矩阵K的__________。 （A）对称性（B）稀疏性（C）奇异性 三．简答题（共20分，每题5分） 1、简述有限单元法结构刚度矩阵的特点。 2、简述有限元法中选取单元位移函数（多项式）的一般原则。 3、简述有限单元法的收敛性准则。
《汽车有限元基础》2009-2010二学期考试试卷
《汽车有限元基础》2009-2010第二学期考试试卷 一、填空题 1. 有限元法的基本思想是用个单元的集合来代替原来具有个自由 度的连续体。 2. 单元刚度矩阵K中元素K ij的物理意义：当单元第j个自由度产生而其它自由度固定时，在第i个自由度产生的。 3．按照各杆轴线及外力作用线在空间的位置，杆系结构可分为： 和。4．平面刚架中各单元发生轴向拉压变形及面内的弯曲变形，而且这两种变形相互独立，因此刚架单元可以看成是由单元和单元叠加而成。因此，平面刚架单元的节点位移应包含个平动分量和个转动分量。 5．工程中常用的薄板单元有：单元和单元。6．有限元分析的主要步骤先后为：(1) 网格划分, (2) , (3) 。 7. 单元特性分析的主要内容先后为：(1) 、(2) 、(3) 应力或内力、(4) 、(5) 单元节点载荷。 8．对于弹性变形体，承受的外载荷共有三种：集中载荷、和。在有限元法中，对于没有作用在节点上的这些外载荷，是按照的原则将其移置到节点上。 9．工程中任一平板，若其厚度为t，板面宽度为b，当t/b小于时可以认为是薄板。常用的薄板单元有：单元和单元。10．薄壳单元中的应力可看成平面应力问题和问题中两种应力的叠加。 11．求解结构系统的动力响应时，常用的两种求解方法为：和 12．在有限元分析中，为了描述几何模型和有限元模型，需要用到几种坐标系： (1) (2) (3) 和(4)
《汽车有限元基础》2009-2010第二学期考试试卷 二、 问答题 1．某一薄板矩形单元的节点编号按照逆时针依次为i 、j 、m 和p 。假设该单元每个节点的位移表示为{}{}T yi xi i i w θθδ=， (i, j, m, p )；该单元每个节点的载荷表示为{}{}T iy ix i i T T Z F θθ=，(i, j, m, p )。请写出该单元的单元节点位移列阵和单元 节点载荷列阵。 2．请写出使用有限元分析软件时，进行数据前处理的主要工作内容。 3．右下图为一典型三节点三角形平面单元，节点按照逆时针依次编号为i 、j 和m ，节点的坐标依次为(x i ，y i )，(x j ，y j )、（x m ，y m ）。假设单元内任意一点的两个位移分量分别表示u 和v 。请写出该单元位移模式的多项式形式，并简述待定常数个数的确定理由。 4. 请简述针对动力问题的有限元分析的基本步骤。
一 判断题（20分） （×）1. 节点的位置依赖于形态，而并不依赖于载荷的位置 （√）2. 对于高压电线的铁塔那样的框架结构的模型化处理使用梁单元 （×）3. 不能把梁单元、壳单元和实体单元混合在一起作成模型 （√）4. 四边形的平面单元尽可能作成接近正方形形状的单元 （×）5. 平面应变单元也好，平面应力单元也好，如果以单位厚来作模型化 处理的话会得到一样的答案 （×）6. 用有限元法不可以对运动的物体的结构进行静力分析 （√）7. 一般应力变化大的地方单元尺寸要划的小才好 （×）8. 所谓全约束只要将位移自由度约束住，而不必约束转动自由度 （√）9. 同一载荷作用下的结构，所给材料的弹性模量越大则变形值越小 （√）10一维变带宽存储通常比二维等带宽存储更节省存储量。 二、填空（20分） 1．平面应力问题与薄板弯曲问题的弹性体几何形状都是 薄板 ，但前者受力特点是： 平行于板面且沿厚度均布载荷作用 ，变形发生在板面内； 后者受力特点是： 垂直于板面 的力的作用，板将变成有弯有扭的曲面。 2．平面应力问题与平面应变问题都具有三个独立的应力分量： σx ，σy ，τxy ，三个独立的应变分量：εx ，εy ，γxy ，但对应的弹性体几何形状前者为 薄板 ，后者为 长柱体 。3．位移模式需反映 刚体位移 ，反映 常变形 ，满足 单元边界上位移连续 。 4．单元刚度矩阵的特点有：对称性 ， 奇异性 ，还可按节点分块。 5．轴对称问题单元形状为：三角形或四边形截面的空间环形单元 ，由于轴对称的特性，任意一点变形只发生在子午面上，因此可以作为 二 维问题处理。 6．等参数单元指的是：描述位移和描述坐标采用相同的形函数形式。等参数单元优点是：可以采用高阶次位移模式，能够模拟复杂几何边界，方便单元刚度矩阵和等效节点载荷的积分运算。 7．有限单元法首先求出的解是 节点位移 ，单元应力可由它求得，其计算公式为 {}{} [][]e D B σδ=。（用符号表示即可） 8．一个空间块体单元的节点有 3 个节点位移： u ，v ，w 9．变形体基本变量有位移应变应力 基本方程 平衡方程 物理方程 几何方程 10.实现有限元分析标准化和规范化的载体就是单元
1.结点的位置依赖于形态，而并不依赖于载荷的位置（×） 2.对于高压电线的铁塔那样的框架结构的模型化处理使用梁单元。√ 3.平面应变单元也好，平面应力单元也好，如果以单位厚来作模型化处理的话会得到一样的答案（×） 4.用有限元法不可以对运动的物体的结构进行静力分析（×） 5.一般应力变化大的地方单元尺寸要划的小才好（√） 6.四结点四边形等参单元的位移插值函数是坐标x、y的一次函数√ 7.在三角形单元中其面积坐标的值与三结点三角形单元的结点形函数值相等。√ 8.等参单元中Jacobi行列式的值不能等于零。√ 9.四边形单元的Jacobi行列式是常数。× 10.等参元是指单元坐标变换和函数插值采用相同的结点和相同的插值函数。√ 11.有限元位移模式中，广义坐标的个数应与单元结点自由度数相等√ 12.为了保证有限单元法解答的收敛性，位移函数应具备的条件是位移函数必须能反映单元的刚体位移和常量应变以及尽可能反映单元间的位移连续性。√ 13.在平面三结点三角形单元中，位移、应变和应力具有位移呈线形变化，应力和应变为常量特征。√ 1.梁单元和杆单元的区别？(自己分析:自由度不同)杆单元只能承受拉压荷载，梁单元则可以承受拉压弯扭荷载。具体的说，杆单元其实就是理论力学常说的二力杆，它只能在结点受载荷，且只有结点上的荷载合力通过其轴线时，杆件才有可能平衡，像均布荷载、中部集中荷载等是无法承担的，通常用于网架、桁架的分析；而梁单元则基本上适用于各种情况（除了楼板之类），且经过适当的处理（如释放自由度、耦合等），梁单元也可以当作杆单元使用。 2.有限单元法结构刚度矩阵的特点？对称性，奇异性，主对角元恒正，稀疏性，非零元素呈带状分布。 3.有限单元法的收敛性准则？完备性要求，协调性要求。位移模式要满足以下三个条件包含单元的刚体位移。当结点位移由体位移引起时，弹性体内不会产生应变。包含单元的常应变。与位置坐标无关的应变。位移模式在单元内要连续，在相邻单元之间的位移必须协调。当选择多项式来构成位移模式时，单元的连续性总得到满足，单元的协调性就是要求单元之间既不会出现开裂也不会出现重叠的现象。。 4.任何一个有限元分析问题都是空间问题，什么情况下可以简化为平面问题？轴对称问题？空
1、何为有限元法？其基本思想是什么？ 有限元法是一种基于变分法而发展起来的求解微分方程的数值计算方法，该方法以计算机为手段，采用分片近似，进而逼近整体的研究思想求解物理问题。 基本思想是化整为零集零为整。 2、为什么说有限元法是近似的方法，体现在哪里？ 有两点：用离散单元的组合体来逼近原始结构，体现了几何上的近似；而用近似函数逼近未知变量在单元内的真实解，体现了数学上的近似。 3、单元、节点的概念？ 节点：表达实际结构几何对象之间相互连接方式的概念 单元：网格划分中的每一个小部分称为单元，网格间相互联结点称为节点 4、有限元法分析过程可归纳为几个步骤？ 结构离散化、单元分析、整体分析 5、有限元方法分几种？本课程讲授的是哪一种？ 位移法、力法、混合法本课程讲授位移法 6、弹性力学的基本变量是什么？何为几何方程、物理方程及虚功方程？弹性矩阵的特点？ 弹性力学变量:外力、应力、应变和位移。 描述弹性体应变分量与位移分量之间的方程称为几何方程；物理方程描述应力分量与应变分量之间的关系；弹性体上外力在虚位移发生过程中所做的虚功与储存在弹性体内的需应变能相等。 弹性矩阵由材料的弹性模量和泊松比确定，与坐标位置无关。 7、何为平面应力问题和平面应变问题？ 平面应力问题：在结构上满足a几何条件：研究对象是等厚度薄板。b载荷条件：作用于薄板上的载荷平行于板平面且沿厚度方向均匀分布，而在两板面无外力作用。 平面应变问题：满足a几何条件：长柱体，即长度方向的尺寸远远大于横截面的尺寸，且横截面沿长度方向不变。b载荷条件：作用于长柱体结构上的载荷平行于横截面且沿纵向方向均匀分布，两端面不受力两条件的弹性力学问题。 1、何为结构的离散化？离散化的目的？何为有限元模型？ ①离散化:把连续的结构看成由有限个单元组成的集合体。②目的：建立有限元计算模型③通常把由节点,单元及相应的节点载荷和节点约束构成的模型称为有限元模型2、结构离散化时，划分单元数目的多少以及疏密分布，将直接影响到什么？确定单元数量的原则？通常如何设置节点？
例1某三角形3节点单元，3个节点的坐标分别为(0, 0)、(3, 0)、(0, 4)o 作用在单元内的点(1, 1)处作用有一个大小等于10N、方向沿x轴正向的集中 力P。求该集中力的等效节点荷载。 解：(1)形函数及形函数矩阵计算 根据面积比或形函数公式，可计算得到各形函数为: Ar 5 Ar 4 3 A, = —、A 9——、A’ 3 =— 12 12 °12 形函数矩阵为: (2)等效节点荷载计算 例2图示三角形3节点单元，设13边的长度为3m,在13边作用有如图所示的分布荷载，求该分布荷载的等效节点荷载。 解：(1)形函数及形函数矩阵计算 在13边上，节点2的形函数％ = 0,设t为节点1到节点3的位置参数, 在节点1处取0，在节点3处取1,则在13边上有： A* =! —/■> M = t 1 O
-20此一10N3 - 20(1 -/- Wt lot - 20 =3f 1 -1 1 一 I 10￡ - 20 dt = 3 Jo d - t) (io* - 20)' o *(10* - 20) >d t 25 T o o o 20 5 1 - t0 0 0 Z 0 LA;J = 0 1 — f 0 0 0 匕 (2)分布荷载的参数表示 在13 边上，q v = 0 ,设 / = a + bt ,由: t = 0, q x = -20 t = 1 , q x = -10 可求得： a = -20 b = 10 于是有: lOt - 20 Cly 另外，由于在边上为一次函数，也可直接根据形函数插值建立分布函数: （3）等效节点荷载计算 \dl -25 - 20
有限元复习 一、选择题（每题1分，共10分） 二、判断题（每空1分，共10分） 三、填空题（每空1分，共10分） 三、简答题（共44分）共6题 四、综述题（共26分）两题 一．基本概念 1. 平面应力/平面应变问题；空间问题/轴对称问题；杆梁问题；线 性与非线性问题 平面应力问题 （1） 均匀薄板（2）载荷平行于板面且沿厚度方向均匀分布 在六个应力分量中，只需要研究剩下的平行于XOY 平面的三个应力分量，即x y xy yx σσττ=、、 （000z zx xz zy yz σττττ=====，，）。 一般0z σ=，z ε并不一定等于零，但可由x σ及y σ求得，在分析问题时不必考虑。于是只需要考虑 x y xy εεγ、、三个应变分量即可。 平面应变问题
（1） 纵向很长，且横截面沿纵向不变。（2）载荷平行于横截面且沿纵向 均匀分布 z yz zx εγγ===只剩下三个应变分量x y xy εεγ、、。也只需要考虑x y xy σστ、、三个应力分量即可 轴对称问题 物体的几何形状、约束情况及所受外力都对称于空间的某一根轴。 轴对称单元的特点（与平面三角形单元的区别）：轴对称单元为圆环体，单元与单元间为节圆相连接；节点力与节点载荷是施加于节圆上的均布力；单元边界是一回转面；应变不是常量。 在轴对称问题中，周向应变分量θε是与r 有关。 板壳问题 一个方向的尺寸比另外两个方向尺寸小很多，且能承受弯矩的结构称为板壳结构，并把平分板壳结构上下表面的面称为中面。如果中面是平面或平面组成的折平面，则称为平板；反之，中面为曲面的称为壳。 杆梁问题 杆梁结构是指长度远大于其横断面尺寸的构件组成的系统。在结构力学中常将承受轴力或扭矩的杆件称为杆，而将承受横向力和弯矩的杆件称为梁。 平面（应力应变）问题与板壳问题的区别与联系 平面应力问题是指很薄的等厚度薄板，只在板边上受有平行于板面并且不沿厚度变化的面力，同时，体力也平行于板面并且不沿厚度变化。而平面应变问题是指很长的柱形体，在柱面上受有平行于横截面并且不沿长度变化的面力，同时
有限元法基础试题（A ） 一、填空题（5×2分） 1.1单元刚度矩阵e T k B DBd Ω = Ω? 中，矩阵B 为__________，矩阵D 为___________。 1.2边界条件通常有两类。通常发生在位置完全固定不能转动的情况为_______边界，具体指定有限的非零值位移的情况，如支撑的下沉，称为_______边界。 1.3内部微元体上外力总虚功： ()(),,,,e x x xy y bx xy x y y by d W F u F v dxdy δστδτσδ??=+++++??+(),,,,x x y y xy y x u v u u dxdy σδσδτδδ??+++??的表达式中，第一项为____________________的虚功，第二项为____________________的虚功。 1.4弹簧单元的位移函数1N +2N =_________。 1.5 ij k 数学表达式：令j d =_____，k d =_____，k j ≠，则力i ij F k =。 二、判断题（5×2分） 2.1位移函数的假设合理与否将直接影响到有限元分析的计算精度、效率和可靠性。（ ） 2.2变形体虚功原理适用于一切结构（一维杆系、二维板、三位块体）、适用于任何力学行为的材料（线性和非线性），是变形体力学的普遍原理。 （ ） 2.3变形体虚功原理要求力系平衡，要求虚位移协调，是在“平衡、协调”前提下功的恒等关系。 （ ） 2.4常应变三角单元中变形矩阵是x 或y 的函数。 （ ） 2.5 对称单元中变形矩阵是x 或y 的函数。 （ ） 三、简答题（26分） 3.1列举有限元法的优点。（8分） 3.2写出有限单元法的分析过程。（8分） 3.3列出3种普通的有限元单元类型。（6分） 3.4简要阐述变形体虚位移原理。（4分） 四、计算题（54分） 4.1对于下图所示的弹簧组合，单元①的弹簧常数为10000N/m ，单元②的弹簧常数为20000N/m ，单元③的弹簧常数为10000N/m ，确定各节点位移、反力以及单元②的单元力。（10分） 4.2对于如图所示的杆组装，弹性模量E 为10GPa ，杆单元长L 均为2m ，横截面面积A 均为2×10-4m 2，弹簧常数为2000kN/m ，所受荷载如图。采用直接刚度法确定节点位移、作用力和单元②的应力。（10分）
《结构有限元分析》复习题（闭卷） 一、绪论 1.概述有限元法分析问题的过程。 二、平面问题 2.对平面问题T3单元，推导其位移模式。 3.对平面问题T3单元，证明形函数在本节点取值为1，在其它节点取值为0。 4.对平面问题T3单元，证明形函数在任意一点上取值之和为1。 5.对平面问题T3单元，证明边界上一点的形函数，与相对顶点的坐标无关。 6.对平面问题T3单元，证明边界上的位移协调性。 7.对平面问题T3单元，说明单元边界上无限点的约束等效于对该边节点的约束。 8.对平面问题T3单元，证明Li=Ni（i=i、j、m）。 9.对平面问题T3单元，证明∑NiXi=X，∑NiYi=Y。 10.对平面问题T3单元，利用最小势能原理，推导单元刚度矩阵的矩阵表达式。 11.说明刚度矩阵的性质和物理意义。 12.对平面问题T3单元，推导单元自重的等效节点力。 13.对平面问题T3单元，推导单元边界上均布压力的等效节点力。 14.对平面问题T3单元，推导单元边界上三角形分布压力的等效节点力。 15.对平面问题R4单元，推导其位移模式。 16.对平面问题R4单元，证明边界上的位移协调性。 17.试写出处理约束的两种方法（划0置1法，乘大数法）的过程。 三、空间问题和轴对称问题 18.对轴对称问题T3单元，推导其位移模式。 19.对轴对称问题T3单元，采用简化计算，推导单元自重的等效节点力。 20.对轴对称问题T3单元，采用简化计算，推导离心力的等效节点力。 21.对轴对称问题T3单元，采用简化计算，推导边界上梯形分布压力的等效节点力。 四、等参单元 22.对平面问题Q4等参单元，构造其位移模式。 23.对平面问题Q4等参单元，推导其几何矩阵。 24.对平面问题Q4等参单元，说明雅可比行列式的意义，并加以数学证明。 25.对平面问题Q4等参单元，证明其完备性、协调性。 26.对平面问题Q4-8变节点等参单元，构造其形函数。 27.对空间问题Hex8-20变节点等参单元，构造其形函数。
有限元试题及答案
一判断题（20分） （×）1. 节点的位置依赖于形态，而并不依赖于载荷的位置 （√）2. 对于高压电线的铁塔那样的框架结构的模型化处理使用梁单元 （×）3. 不能把梁单元、壳单元和实体单元混合在一起作成模型 （√）4. 四边形的平面单元尽可能作成接近正方形形状的单元 （×）5. 平面应变单元也好，平面应力单元也好，如果以单位厚来作模型化 处理的话会得到一样的答案 （×）6. 用有限元法不可以对运动的物体的结构进行静力分析 （√）7. 一般应力变化大的地方单元尺寸要划的小才好 （×）8. 所谓全约束只要将位移自由度约束住，而不必约束转动自由度 （√）9. 同一载荷作用下的结构，所给材料的弹性模量越大则变形值越小（√）10一维变带宽存储通常比二维等带宽存储更节省存储量。 二、填空（20分） 1．平面应力问题与薄板弯曲问题的弹性体几何形状都是薄板，但前者受力特点是：平行于板面且沿厚度均布载荷作用，变形发生在板面内； 后者受力特点是：垂直于板面的力的作用，板将变成有弯有扭的曲面。 2．平面应力问题与平面应变问题都具有三个独立的应力分量：σx，σy，τxy ，三个独立的应变分量：εx，εy，γxy，但对应的弹性体几何形状前者为薄板，后者为长柱体。3．位移模式需反映刚体位移，反映常变形，满足单元边界上位移连续。 4．单元刚度矩阵的特点有：对称性，奇异性，还可按节点分块。
5．轴对称问题单元形状为：三角形或四边形截面的空间环形单元，由于轴对称的特性，任意一点变形只发生在子午面上，因此可以作为二维问题处理。6．等参数单元指的是：描述位移和描述坐标采用相同的形函数形式。等参数单元优点是：可以采用高阶次位移模式，能够模拟复杂几何边界，方便单元刚度矩阵和等效节点载荷的积分运算。 7．有限单元法首先求出的解是节点位移，单元应力可由它求得，其计算公式为。（用符号表示即可） 8．一个空间块体单元的节点有 3 个节点位移： u，v，w 9．变形体基本变量有位移应变应力基本方程平衡方程物理方程几何方程 10.实现有限元分析标准化和规范化的载体就是单元
1.针对下图所示的3个三角形元,写出用完整多项式描述的位移模式表达式。 ２.如下图所示,求下列情况的带宽: a)4结点四边形元; b)2结点线性杆元。 ３、对上题图诸结点制定一种结点编号的方法,使所得带宽更小。图左下角的四边形在两种不同编号方式下,单元的带宽分别就是多大? 4、下图所示,若单元就是2结点线性杆单元,勾画出组装总刚后总刚空间轮廓线。系统的带宽就是多大?按一右一左重新编号(即6变成3等)后,重复以上运算。
５. 设杆件1-2受轴向力作用,截面积为Ａ,长度为Ｌ,弹性模量为Ｅ,试写出杆端力Ｆ１,F 2与杆端位移21,u u 之间的关系式,并求出杆件的单元刚度矩阵)(][e k 6、设阶梯形杆件由两个等截面杆件＼o ＼a ｃ(○,1)与错误!所组成,试写出三个结点1、2、3的结点轴向力F 1,F 2,Ｆ3与结点轴向位移321,,u u u 之间的整体刚度矩阵［K]。 7. 在上题的阶梯形杆件中,设结点3为固定端,结点1作用轴向载荷F 1=Ｐ,求各结点的轴向位移与各杆的轴力。 8、 下图所示为平面桁架中的任一单元,y x ,为局部坐标系,ｘ,y 为总体坐标系,x 轴与x 轴的夹角为 。 (1) 求在局部坐标系中的单元刚度矩阵 )(][e k (2) 求单元的坐标转换矩阵 ［T]; (３) 求在总体坐标系中的单元刚度矩阵 )(][e k
９.如图所示一个直角三角形桁架,已知27/103cm N E ?=,两个直角边长度cm l 100=,各杆截面面积210cm A =,求整体刚度矩阵［K ］ 。 10. 设上题中的桁架的支承情况与载荷情况如下图所示,按有限元素法求出各结点的位移与各杆的内力。
重庆大学研究生有限 元大作业
课程研究报告 科目：有限元分析技术教师：阎春平姓名：色学号： 2 专业：机械工程类别：学术 上课时间： 2015 年 11 月至 2016 年 1 月 考生成绩： 阅卷评语： 阅卷教师 (签名)
有限元分析技术作业 姓名: 色序号: 是学号: 2 一、题目描述及要求 钢结构的主梁为高160宽100厚14的方钢管，次梁为直径60厚10的圆钢管（单位为毫米），材料均为碳素结构钢Q235；该结构固定支撑点位于左右两端主梁和最中间。主梁和次梁之间是固接。试对在垂直于玻璃平面方向的2kPa 的面载荷（包括玻璃自重、钢结构自重、活载荷（人员与演出器械载荷）、风载荷等）作用下的舞台进行有限元分析。 二、题目分析 根据序号为069，换算得钢结构框架为11列13行。由于每个格子的大小为1×1（单位米），因此框架的外边框应为11000×13000（单位毫米）。 三、具体操作及分析求解 1、准备工作 执行Utility Menu：File → Clear＆start new 清除当前数据库并开始新的分析，更改文件名和文件标题，如图1.1。选择GUI filter，执行 Main Menu: Preferences → Structural → OK，如图1.2所示
图1.1清除当前数据库并开始新的分析 图1.2 设置GUI filter 2、选择单元类型。 执行Main Menu: Preprocessor →Element Type →Add/Edit/Delete →Add→ select→ BEAM188，如图2.1。之后点击OK(回到Element Types window) →Close
1、何为有限元法？其基本思想是什么? 有限元法是一种基于变分法而发展起来的求解微分方程的数值计算方法，该方法以计算机为手段，采用分片近似，进而逼近整体的研究思想求解物理问题。 基本思想是化整为零集零为整。 2、为什么说有限元法是近似的方法，体现在哪里? 有两点：用离散单元的组合体来逼近原始结构，体现了几何上的近似；而用近似函数逼近未知变量在单元内的真实解，体现了数学上的近似。 3、单元、节点的概念? 节点：表达实际结构几何对象之间相互连接方式的概念 单元：网格划分中的每一个小部分称为单元，网格间相互联结点称为节点 4、有限元法分析过程可归纳为几个步骤? 结构离散化、单元分析、整体分析 5、有限元方法分几种？本课程讲授的是哪一种? 位移法、力法、混合法本课程讲授位移法 6、弹性力学的基本变量是什么？何为几何方程、物理方程及虚功方程？弹性矩阵的特点？ 弹性力学变量：外力、应力、应变和位移。 描述弹性体应变分量与位移分量之间的方程称为几何方程；物理方程描述应力分量与应变分量之间的关系；弹性体上外力在虚位移发生过程中所做的虚功与储存在弹性体内的需应变能相等。 弹性矩阵由材料的弹性模量和泊松比确定，与坐标位置无关。 7、何为平面应力问题和平面应变问题? 平面应力问题：在结构上满足a几何条件：研究对象是等厚度薄板。b载荷条件: 作用于薄板上的载荷平行于板平面且沿厚度方向均匀分布，而在两板面无外力作用。 平面应变问题：满足a几何条件：长柱体，即长度方向的尺寸远远大于横截面的尺寸，且横截面沿长度方向不变。b载荷条件：作用于长柱体结构上的载荷平行于横截面且沿纵向方向均匀分布，两端面不受力两条件的弹性力学问题。 1、何为结构的离散化？离散化的目的？何为有限元模型? ①离散化：把连续的结构看成由有限个单元组成的集合体。②目的：建立有限元计算模型③通常把由节点，单元及相应的节点载荷和节点约束构成的模型称为有限元模型 2、结构离散化时，划分单元数目的多少以及疏密分布，将直接影响到什么？确定单元数量的原则？通常如何设置节点？
姓名：色学号：2 专业：机械工程类别：学术 上课时间：2015年11 月至2016 年 1 月考生成绩： 阅卷评语： 阅卷教师(签名)
有限元分析技术作业 姓名: 色序号: 是学号: 2 一、题目描述及要求 钢结构的主梁为高160宽100厚14的方钢管，次梁为直径60厚10的圆钢管（单位为毫米），材料均为碳素结构钢Q235；该结构固定支撑点位于左右两端主梁和最中间。主梁和次梁之间是固接。试对在垂直于玻璃平面方向的2kPa的面载荷（包括玻璃自重、钢结构自重、活载荷（人员与演出器械载荷）、风载荷等）作用下的舞台进行有限元分析。 二、题目分析 根据序号为069，换算得钢结构框架为11列13行。由于每个格子的大小为1×1（单位米），因此框架的外边框应为11000×13000（单位毫米）。 三、具体操作及分析求解 1、准备工作 执行Utility Menu：Fi le→Clear＆start new 清除当前数据库并开始新的分析，更改文件名和文件标题，如图 1.1。选择GUI filter，执行Main Menu: Preferences→Structural→OK，如图1.2所示 图1.1清除当前数据库并开始新的分析 1
2 图1.2 设置GUI filter 2、选择单元类型。 执行Main Menu: Preprocessor →Element Type →Add/Edit/Delete →Add→ select→ BEAM188，如图2.1。之后点击OK(回到Element Types window) →Close 图2.1选择单元 3、定义材料属性 该钢结构材料为碳素结构钢Q235，其弹性模量为210GPa ，执行Main Menu→Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic ，此处协调单位制为mmkgs ，故EX 设为2.1E8, PRXY 设置