有限元基础(期末考试题)

《有限元基础》期末测试

一、结构线性静力分析

如图所示的托架，其顶面承受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。试研究该梁板结构的瞬

态动力响应。

图

图2

建立有限元分析模型并附加动力

节点146的位移时间历程结果

三、非线性屈曲分析

如图，一根长200L in =，截面高度0.5h in =，截面面积20.25

A in =，惯性矩24/120.0052083J Ah in ==的细长杆受轴向载荷的作用，若沿X 方向取10个主

自由度，求其屈曲模态。

变形图

进行建模，施加边界条件和载荷，进行分析得：

从输出窗口中，可查到屈曲荷载系数为 38.552928。

四、热分析

MCM的结构如图所示：其中大芯片热流密度60.1×106W/m3，小芯片为61.54×106W/m3，

周围空气温度25，对流换热系数10W/(m2.K)。

MCM截面图MCM俯视图

1、建模网格划分

MCM网格划分图

分析处理(温度云图)结果

有限元基础(期末考试题)

《有限元基础》期末测试一、结构线性静力分析如图所示的托架，其顶面承受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

北京科技大学有限元考试试题

一．是非题（认为该题正确，在括号中打√；该题错误，在括号中打×。）(每小题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、简述有限单元法的收敛性准则。

有限元法课后习题答案

1、有限元是近似求解一般连续场问题的数值方法 2、有限元法将连续的求解域离散为若干个子域,得到有限个单元,单元和单元之间用节点连接 3、直梁在外力的作用下,横截面的内力有剪力和弯矩两个. 4、平面刚架结构在外力的作用下,横截面上的内力有轴力、剪力、弯矩. 5、进行直梁有限元分析,平面刚架单元上每个节点的节点位移为挠度和转角 6、平面刚架有限元分析，节点位移有轴向位移、横向位移、转角。 7、在弹性和小变形下，节点力和节点位移关系是线性关系。 8、弹性力学问题的方程个数有15个，未知量个数有15个。 9、弹性力学平面问题方程个数有8，未知数8个。 10、几何方程是研究应变和位移之间关系的方程 11、物理方程是描述应力和应变关系的方程 12、平衡方程反映了应力和体力之间关系的 13、把经过物体内任意一点各个截面上的应力状况叫做一点的应力状态 14、9形函数在单元上节点上的值,具有本点为_1_.它点为零的性质,并且在三角形单元的任一节点上,三个行函数之和为_1_ 15、形函数是_三角形_单元内部坐标的_线性_函数,他反映了单元的_位移_状态 16、在进行节点编号时,同一单元的相邻节点的号码差尽量小. 17、三角形单元的位移模式为_线性位移模式_- 18、矩形单元的位移模式为__双线性位移模式_

19、在选择多项式位移模式的阶次时,要求_所选的位移模式应该与局部坐标系的方位无关的性质为几何_各向同性 20、单元刚度矩阵描述了_节点力_和_节点位移之间的关系 21、矩形单元边界上位移是连续变化的 1.诉述有限元法的定义答：有限元法是近似求解一般连续场问题的数值方法 2.有限元法的基本思想是什么答：首先，将表示结构的连续离散为若干个子域，单元之间通过其边界上的节点连接成组合体。其次，用每个单元内所假设的近似函数分片地表示求解域内待求的未知厂变量。 3.有限元法的分类和基本步骤有哪些答：分类：位移法、力法、混合法；步骤：结构的离散化，单元分析，单元集成，引入约束条件，求解线性方程组，得出节点位移。 4.有限元法有哪些优缺点答：优点：

有限元法试题

《汽车有限元基础》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. 请简述针对动力问题的有限元分析的基本步骤。

重庆大学研究生有限元复习题及答案(2013)

1.结点的位置依赖于形态，而并不依赖于载荷的位置（×） 2.对于高压电线的铁塔那样的框架结构的模型化处理使用梁单元。√ 3.平面应变单元也好，平面应力单元也好，如果以单位厚来作模型化处理的话会得到一样的答案（×） 4.用有限元法不可以对运动的物体的结构进行静力分析（×） 5.一般应力变化大的地方单元尺寸要划的小才好（√） 6.四结点四边形等参单元的位移插值函数是坐标x、y的一次函数√ 7.在三角形单元中其面积坐标的值与三结点三角形单元的结点形函数值相等。√ 8.等参单元中Jacobi行列式的值不能等于零。√ 9.四边形单元的Jacobi行列式是常数。× 10.等参元是指单元坐标变换和函数插值采用相同的结点和相同的插值函数。√ 11.有限元位移模式中，广义坐标的个数应与单元结点自由度数相等√ 12.为了保证有限单元法解答的收敛性，位移函数应具备的条件是位移函数必须能反映单元的刚体位移和常量应变以及尽可能反映单元间的位移连续性。√ 13.在平面三结点三角形单元中，位移、应变和应力具有位移呈线形变化，应力和应变为常量特征。√ 1.梁单元和杆单元的区别？(自己分析:自由度不同)杆单元只能承受拉压荷载，梁单元则可以承受拉压弯扭荷载。具体的说，杆单元其实就是理论力学常说的二力杆，它只能在结点受载荷，且只有结点上的荷载合力通过其轴线时，杆件才有可能平衡，像均布荷载、中部集中荷载等是无法承担的，通常用于网架、桁架的分析；而梁单元则基本上适用于各种情况（除了楼板之类），且经过适当的处理（如释放自由度、耦合等），梁单元也可以当作杆单元使用。 2.有限单元法结构刚度矩阵的特点？对称性，奇异性，主对角元恒正，稀疏性，非零元素呈带状分布。 3.有限单元法的收敛性准则？完备性要求，协调性要求。位移模式要满足以下三个条件包含单元的刚体位移。当结点位移由体位移引起时，弹性体内不会产生应变。包含单元的常应变。与位置坐标无关的应变。位移模式在单元内要连续，在相邻单元之间的位移必须协调。当选择多项式来构成位移模式时，单元的连续性总得到满足，单元的协调性就是要求单元之间既不会出现开裂也不会出现重叠的现象。。 4.任何一个有限元分析问题都是空间问题，什么情况下可以简化为平面问题？轴对称问题？空

《有限元基础教程》_【MATLAB算例】3.3.7(2)__三梁平面框架结构的有限元分析(Beam2D2Node)

【MA TLAB 算例】3.3.7(2) 三梁平面框架结构的有限元分析 (Beam2D2Node) 如图3-19所示的框架结构，其顶端受均布力作用，结构中各个截面的参数都为：113.010Pa E =?，746.510I m -=?，426.810A m -=?。试基于MA TLAB 平台求解该结构的节点位移以及支反力。图3-19 框架结构受一均布力作用解答：对该问题进行有限元分析的过程如下。（1）结构的离散化与编号将该结构离散为3个单元，节点位移及单元编号如图3-20所示，有关节点和单元的信息见表3-5。（a ）节点位移及单元编号

（b）等效在节点上的外力图3-20 单元划分、节点位移及节点上的外载（2）各个单元的描述首先在MA TLAB环境下，输入弹性模量E、横截面积A、惯性矩I和长度L，然后针对单元1，单元2和单元3，分别二次调用函数Beam2D2Node_ElementStiffness，就可以得到单元的刚度矩阵k1(6×6)和k2(6×6)，且单元2和单元3的刚度矩阵相同。 >> E=3E11; >> I=6.5E-7; >> A=6.8E-4; >> L1=1.44; >> L2=0.96; >> k1=Beam2D2Node_Stiffness(E,I,A,L1); >> k2=Beam2D2Node_Stiffness(E,I,A,L2); （3）建立整体刚度方程将单元2和单元3的刚度矩阵转换成整体坐标下的形式。由于该结构共有4个节点，则总共的自由度数为12，因此，结构总的刚度矩阵为KK(12×12)，对KK清零，然后两次调用函数Beam2D2Node_Assemble进行刚度矩阵的组装。 >> T=[0,1,0,0,0,0;-1,0,0,0,0,0;0,0,1,0,0,0;0,0,0,0,1,0;0,0,0,-1,0,0;0,0,0,0,0,1] ; >> k3=T'*k2*T; >> KK=zeros(12,12); >> KK=Beam2D2Node_Assemble(KK,k1,1,2);

(完整版)福州大学有限元考试题

一判断题（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.1 有限单元法中“离散”的含义是什么？有限单元法是如何将具有无限自由度的连续介质问题转变成有限自由度问题的？位移有限元法的标准化程式是怎样的？（1）离散的含义即将结构离散化，即用假想的线或面将连续体分割成数目有限的单元，并在其上设定有限个节点；用这些单元组成的单元集合体代替原来的连续体，而场函数的节点值将成为问题的基本未知量。（2）给每个单元选择合适的位移函数或称位移模式来近似地表示单元内位移分布规律，即通过插值以单元节点位移表示单元内任意点的位移。因节点位移个数是有限的，故无限自由度问题被转变成了有限自由度问题。（3）有限元法的标准化程式：结构或区域离散，单元分析，整体分析，数值求解。 1.3 单元刚度矩阵和整体刚度矩阵各有哪些性质？各自的物理意义是什么?两者有何区别?单元刚度矩阵的性质：对称性、奇异性（单元刚度矩阵的行列式为零）。整体刚度矩阵的性质：对称性、奇异性、稀疏性。单元 Kij 物理意义 Kij 即单元节点位移向量中第 j 个自由度发生单位位移而其他位移分量为零时，在第 j 个自由度方向引起的节点力。整体刚度矩阵 K 中每一列元素的物理意义是：要迫使结构的某节点位移自由度发生单位位移，而其他节点位移都保持为零的变形状态，在所有个节点上需要施加的节点荷载。 2.2 什么叫应变能？什么叫外力势能？试叙述势能变分原理和最小势能原理，并回答下述问题：势能变分原理代表什么控制方程和边界条件？其中附加了哪些条件？ (1)在外力作用下，物体内部将产生应力σ和应变ε，外力所做的功将以变形能的形式储存起来，这种能量称为应变能。 (2)外力势能就是外力功的负值。 (3)势能变分原理可叙述如下：在所有满足边界条件的协调位移中，那些满足静力平衡条件的位移使物体势能泛函取驻值，即势能的变分为零 δ∏p=δ Uε+δV=0 此即变分方程。对于线性弹性体，势能取最小值，即 δ2∏P=δ2Uε+δ2V≥0 此时的势能变分原理就是著名的最小势能原理。势能变分原理代表平衡方程、本构方程和应力边界条件，其中附加了几何方程和位移边界条件。 2.3 什么是强形式？什么是弱形式？两者有何区别？建立弱形式的关键步骤是什么？等效积分形式通过分部积分，称式 ∫ΩCT(v)D(u)dΩ+∫ΓET(v)F(u)dΓ 为微分方程的弱形式，相对而言，定解问题的微分方程称为强形式。区别：弱形式得不到解析解。建立弱形式的关键步骤：对场函数要求较低阶的连续性。2.4 为了使计算结果能够收敛于精确解，位移函数需要满足哪些条件？为什么？只要位移函数满足两个基本要求，即完备性和协调性，计算结果便收敛于精确解。 2.6 为什么采用变分法求解通常只能得到近似解？变分法的应用常遇到什么困难？Ritz 法收敛的条件是什么？（1）在 Ritz 法中，N 决定了试探函数的基本形态，待定参数使得场函数具有一定的任意性。如果真实场函数包含在试探函数之内，则变分法得到的解答是精确的；如果试探函数取自完全的函数序列，则当项数不断增加时，近似解将趋近于精确解。然而，通常情况下试探函数不会将真实场函数完全包含在内，实际计算时也不可能取无穷多项。因此，试探函数只能是真实场函数的近似。可见，变分法就是在某个假定的范围内找出最佳解答，近似性就源于此。（2）采用变分法近似求解，要求在整个求解区域内预先给出满足边界条件的场函数。通常情况下这是不可能的，因而变分法的应用受到了限制。（3）Ritz 法的收敛条件是要求试探函数具有完备性和连续性，也就是说，如果试探函数满足完备性和连续性的要求，当试探函数的项数趋近于无穷时，则 Ritz 法的近似解将趋近于数学微分方程的精确解。 3.1 构造单元形函数有哪些基本原则？形函数是定义于单元内坐标的连续函数。单元位移函数通常采用多项式，其中的待定常数应该与单元节点自由度数相等。为满足完备性要求，位移函数中必须包括常函数和一次式，即完全一次多项式。多项式的选取应由低阶到高阶，尽量选择完全多项式以提高单元的精度。若由于项数限制而不能选取完全多项式时，也应使完全多项式具有坐标的对称性，并且一

有限元法基础试题

有限元法基础试题（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分）

Matlab有限元分析操作基础共11页

Matlab有限元分析20140226 为了用Matlab进行有限元分析，首先要学会Matlab基本操作，还要学会使用Matlab进行有限元分析的基本操作。 1. 复习：上节课分析了弹簧系统 x 推导了系统刚度矩阵

2. Matlab有限元分析的基本操作（1）单元划分（选择何种单元，分成多少个单元，标号）（2）构造单元刚度矩阵（列出…）（3）组装系统刚度矩阵（集成整体刚度矩阵）（4）引入边界条件（消除冗余方程）（5）解方程（6）后处理（扩展计算）

3. Matlab有限元分析实战【实例1】

分析：步骤一：单元划分

>>k1=SpringElementStiffness(100)

a) 分析SpringAssemble库函数 function y = SpringAssemble(K,k,i,j) % This function assembles the element stiffness % matrix k of the spring with nodes i and j into the % global stiffness matrix K. % function returns the global stiffness matrix K % after the element stiffness matrix k is assembled. K(i,i) = K(i,i) + k(1,1); K(i,j) = K(i,j) + k(1,2); K(j,i) = K(j,i) + k(2,1); K(j,j) = K(j,j) + k(2,2); y = K; b) K是多大矩阵？今天的系统刚度矩阵是什么？因为 11 22 1212 k k k k k k k k - ?? ?? - ????--+ ?? 所以 1000100 0200200 100200300 - ?? ?? - ????-- ??？

(完整版)有限元第二章课后题答案

2 弹性力学问题的有限单元法思考题 2.1 有限元法离散结构时为什么要在应力变化复杂的地方采用较密网格，而在其他地方采用较稀疏网格？答：在应力变化复杂的地方每一结点与相邻结点的应力都变化较大，若网格划分较稀疏，则在应力突变处没有设置结点，而使得所求解的误差很大，若网格划分较密时，则应力变化复杂的地方可以设置更多的结点，从而使得所求解的精度更高一些。 2.2 因为应力边界条件就是边界上的平衡方程，所以引用虚功原理必然满足应力边界条件，对吗？答：对。 2.3 为什么有限元只能求解位移边值问题和混合边值问题？弹性力学中受内压和外压作用的圆环能用有限元方法求解吗？为什么？答：有限元法是一种位移解法，故只能求解位移边值问题和混合边值问题。而应力边值问题没有确定的位移约束，不能用位移法求解，所以也不能用有限元法求解。 2.4 矩形单元旋转一个角度后还能够保持在单元边界上的位移协调吗？答：能。矩形单元的插值函数满足单元内部和单元边界上的连续性要求，是一个协调元。矩形的插值函数只与坐标差有关，旋转一个角度后各个结点的坐标差保持不变，所以插值函数保持不变。因此矩形单

元旋转一个角度后还能够保持在单元边界上的位移协调。 2.5 总体刚度矩阵呈带状分布，与哪些因素有关？如何计算半带宽？答：因素：总体刚度矩阵呈带状分布与单元内最大结点号与最小结点号的差有关。计算：设半带宽为B ，每个结点的自由度为n ，各单元中结点整体码的最大差值为D ，则B=n(D+1)，在平面问题中n=2。 2.6 为什么单元尺寸不要相差太大，如果这样，会导致什么结果？答：由于实际工程是一个二维或三维的连续体，将其分为具有简单而规则的几何单元，这样便于网格计算，还可以通过增加结点数提高单元精度。在几何形状上等于或近似与原来形状，减小由于形状差异过大带来的误差。若形状相差过大，使结构应力分析困难加大，误差同时也加大。 2.7 剖分网格时，在边界出现突变和有集中力作用的地方要设置结点或单元边界，试说明理由。答：有限元处于弹性力学问题的方法是离散法。它将一个受外力作用的连续弹性体离散成一定数量的有限小的单元集合体，单元之间只在结点上相互联系，即只有结点才能传递力。所以在边界出现突变和有集中力作用的地方要设置结点和单元边界。 2.8 为什么说三角形三结点单元是常应变单元，如果在每边中点增加一个结点，那么单元内应力如何分布？答：(1)应变矩阵[B]中的参数m j i m j i c c c b b b 、、、、、由坐标变量x 、y 之差确定。当单元的坐标差确定之后，这些参数与坐标变量x 、y 无关，

《有限元基础教程》_【ANSYS算例】4.7.1(3) 基于3节点三角形单元的矩形薄板分析(GUI)及命令流

【ANSYS 算例】4.7.1(3) 基于3节点三角形单元的矩形薄板分析如图4-20所示为一矩形薄平板，在右端部受集中力100 000N F =作用，材料常数为：弹性模量7110Pa E =?、泊松比1/3μ=，板的厚度为0.1m t =，在ANSYS 平台上，按平面应力问题完成相应的力学分析。 (a) 问题描述 (a) 有限元分析模型图4–20 右端部受集中力作用的平面问题(高深梁) 解答在ANSYS 平台上，完成的分析如下。 1. 基于图形界面的交互式操作(step by step) (1) 进入ANSYS(设定工作目录和工作文件) 程序 → ANSYS Interactive →Working directory (设置工作目录) →Initial jobname (设置工作文件名): 2D3Node →Run → OK (2) 设置计算类型 ANSYS Main Menu : Preferences… → Structural → OK (3) 选择单元类型 ANSYS Main Menu : Preprocessor →Element Type →Add/Edit/Delete… →Add… →Solid ：Quad 4node 42 →OK (返回到Element Types 窗口) → Options… →K3: Plane Strs w/thk(带厚度的平面应力问题) →OK →Close (4) 定义材料参数 ANSYS Main Menu : Preprocessor →Material Props →Material Models →Structural →Linear →Elastic → Isotropic: EX:1.0e7 (弹性模量)，PRXY: 0.33333333 (泊松比) → OK → 鼠标点击该窗口右上角的“ ”来关闭该窗口 (5) 定义实常数以确定平面问题的厚度 ANSYS Main Menu: Preprocessor →Real Constant s… →Add/Edit/Delete →Add →Type 1→ OK →Real Constant Set No: 1 (第1号实常数), THK: 0.1 (平面问题的厚度) →OK →Close (6) 生成单元模型生成4个节点 ANSYS Main Menu: Preprocessor →Modeling → Create → Nodes → On Working Plane →输入节点1的x,y,z 坐标(2,1,0)，回车→输入节点2的x,y,z 坐标(2,0,0)，回车→输入节点3的x,y,z 坐标(0,1,0)，回车→输入节点4的x,y,z 坐标(0,0,0)，回车→OK 定义单元属性 ANSYS Main Menu: Preprocessor →Modeling → Create → Elements → Elem Attributes →Element type number:1 →Material number:1→Real constant set number:1 →OK 生成单元 ANSYS Main Menu: Preprocessor →Modeling → Create → Elements → User Numbered → Thru Nodes →Number to assign to element:1→Pick nodes:2,3,4→OK →Number to assign to element:2→Pick nodes:3,2,1→OK (7) 模型施加约束和外载左边两个节点施加X,Y 方向的位移约束 ANSYS Main Menu: Solution → Define Loads → Apply →Structural → Displacement → On

西工大有限元试题附答案68872

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. 设上题中的桁架的支承情况与载荷情况如下图所示,按有限元素法求出各结点的位移与各杆的内力。

结构分析及有限元分析基础知识

第一章结构分析及有限元分析基础知识注：摘自《NX知识工程应用技术——CAD/CAE篇》洪如瑾编译清华大学出版社 [目标] 本章将简述结构分析及有限元分析的基础知识，为学习与应用结构分析做好准备，包括： ※ 结构与结构分析定义 ※ 结构的线性静态分析 ※ 材料行为与故障 ※ 有限元分析的基本概念 ※ 有限元模型 1.1结构分析基础知识 1.1.1结构基本概念 1.结构定义结构可以定义为一个正承受作用的载荷处于平衡中的系统。平衡条件意味着结构是不移动的。一个自由的支架不是一个结构，它未被连接到任一物体上并无载荷作用与它。仅当它附着到外部世界，并且有作用力、压力或力矩时，支架成为一个结构。例如横跨江面的大桥就是一个普通的结构，一个支架通过它的支撑连接到地面上，桥的重量是在结构上的一种载荷（力）。当汽车通过桥时，附加的力作用于桥的不同位置。一个好的结构必须满足以下标准：（1）当预期的载荷作用时，结构必须不出现故障。这个似乎是显而易见的，并意味着结构必须是“强度足够的”。故障意味着结构破裂、分离、弯曲，以及支撑作用载荷失败。注意：考虑到意外的载荷，通常在设计中提供安全余量。余量常常利用安全因素来描述。例如，如果在结构上期待载荷是10 000磅，规定安全因素是2.0，则结构将设计成能经受住20 000磅载荷。（2）当载荷作用时，结构必须不产生过分变形。这意味着结构必须“刚度足够”。变形可接受的极限（弯曲度、挠度、拉伸等）取决于特定情况。例如，在通常住宅中的地板由足够的吊带支撑，以防止当人在地板岸上行走时有“柔软”的感觉。（3）在它的服务生命周期，结构的行为应不会恶化。这意味着结构必须“足够耐用”，必须考虑环境影响和“磨损与破裂”。如果一座桥假定维持50年，则桥的设计必须提供整个50年寿命的结构完整性与充分的安全余量。2.结构分析结构分析是用于决定一个结构是否将正确完成任务的工程分析过程。结构将在某些方式中进行模拟和求解描述它的行为的数学方程。分析可以人工方法或用计算机方法来完成。结构分析的结果（答案）用于评估性能，摘要如下：（1）“强度足够吗？”：应力必须是在一可接受的范围内。（2）“刚度足够吗？”：位移必须是在一可接受的范围内。（3）“耐用度足够？”：对一个长的疲劳周期应力必须足够低。

有限元分析基础教程(ANSYS算例)(曾攀)

有限元分析基础教程Fundamentals of Finite Element Analysis (ANSYS算例) 曾攀清华大学 2008-12

有限元分析基础教程曾攀有限元分析基础教程 Fundamentals of Finite Element Analysis 曾攀 (清华大学) 内容简介全教程包括两大部分，共分9章；第一部分为有限元分析基本原理，包括第1章至第5章，内容有：绪论、有限元分析过程的概要、杆梁结构分析的有限元方法、连续体结构分析的有限元方法、有限元分析中的若干问题讨论；第二部分为有限元分析的典型应用领域，包括第6章至第9章，内容有：静力结构的有限元分析、结构振动的有限元分析、传热过程的有限元分析、弹塑性材料的有限元分析。本书以基本变量、基本方程、求解原理、单元构建、典型例题、MATLAB程序及算例、ANSYS算例等一系列规范性方式来描述有限元分析的力学原理、程序编制以及实例应用；给出的典型实例都详细提供有完整的数学推演过程以及ANSYS实现过程。本教程的基本理论阐述简明扼要，重点突出，实例丰富，教程中的二部分内容相互衔接，也可独立使用，适合于具有大学高年级学生程度的人员作为培训教材，也适合于不同程度的读者进行自学；对于希望在MATLAB程序以及ANSYS平台进行建模分析的读者，本教程更值得参考。本基础教程的读者对象：机械、力学、土木、水利、航空航天等专业的工程技术人员、科研工作者。

目录 [[[[[[\\\\\\ 【ANSYS算例】3.3.7(3) 三梁平面框架结构的有限元分析 1 【ANSYS算例】4.3.2(4) 三角形单元与矩形单元的精细网格的计算比较 3 【ANSYS算例】5.3(8) 平面问题斜支座的处理 6 【ANSYS算例】6.2(2) 受均匀载荷方形板的有限元分析9 【ANSYS算例】6.4.2(1) 8万吨模锻液压机主牌坊的分析(GUI) 15 【ANSYS算例】6.4.2(2) 8万吨模锻液压机主牌坊的参数化建模与分析(命令流) 17 【ANSYS算例】7.2(1) 汽车悬挂系统的振动模态分析(GUI) 20 【ANSYS算例】7.2(2) 汽车悬挂系统的振动模态分析(命令流) 23 【ANSYS算例】7.3(1) 带有张拉的绳索的振动模态分析(GUI) 24 【ANSYS算例】7.3(2) 带有张拉的绳索的振动模态分析(命令流) 27 【ANSYS算例】7.4(1) 机翼模型的振动模态分析(GUI) 28 【ANSYS算例】7.4(2) 机翼模型的振动模态分析(命令流) 30 【ANSYS算例】8.2(1) 2D矩形板的稳态热对流的自适应分析(GUI) 31 【ANSYS算例】8.2(2) 2D矩形板的稳态热对流的自适应分析(命令流) 33 【ANSYS算例】8.3(1) 金属材料凝固过程的瞬态传热分析(GUI) 34 【ANSYS算例】8.3(2) 金属材料凝固过程的瞬态传热分析(命令流) 38 【ANSYS算例】8.4(1) 升温条件下杆件支撑结构的热应力分析(GUI) 39 【ANSYS算例】8.4(2) 升温条件下杆件支撑结构的热应力分析(命令流) 42 【ANSYS算例】9.2(2) 三杆结构塑性卸载后的残余应力计算(命令流) 45 【ANSYS算例】9.3(1) 悬臂梁在循环加载作用下的弹塑性计算(GUI) 46 【ANSYS算例】9.3(2) 悬臂梁在循环加载作用下的弹塑性计算(命令流) 49 附录 B ANSYS软件的基本操作52 B.1 基于图形界面(GUI)的交互式操作(step by step) 53 B.2 log命令流文件的调入操作(可由GUI环境下生成log文件) 56 B.3 完全的直接命令输入方式操作56 B.4 APDL参数化编程的初步操作57

有限元考试试题.(优选)

一．是非题（认为该题正确，在括号中打；该题错误，在括号中打×。）(每小题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）非零元素带状分布

有限元复习题答案

1、何为有限元法？其基本思想是什么？有限元法是一种基于变分法而发展起来的求解微分方程的数值计算方法，该方法以计算机为手段，采用分片近似，进而逼近整体的研究思想求解物理问题。基本思想是化整为零集零为整。 2、为什么说有限元法是近似的方法，体现在哪里？有两点：用离散单元的组合体来逼近原始结构，体现了几何上的近似；而用近似函数逼近未知变量在单元内的真实解，体现了数学上的近似。 3、单元、节点的概念？节点：表达实际结构几何对象之间相互连接方式的概念单元：网格划分中的每一个小部分称为单元，网格间相互联结点称为节点 4、有限元法分析过程可归纳为几个步骤？结构离散化、单元分析、整体分析 5、有限元方法分几种？本课程讲授的是哪一种？位移法、力法、混合法本课程讲授位移法 6、弹性力学的基本变量是什么？何为几何方程、物理方程及虚功方程？弹性矩阵的特点？弹性力学变量:外力、应力、应变和位移。描述弹性体应变分量与位移分量之间的方程称为几何方程；物理方程描述应力分量与应变分量之间的关系；弹性体上外力在虚位移发生过程中所做的虚功与储存在弹性体内的需应变能相等。弹性矩阵由材料的弹性模量和泊松比确定，与坐标位置无关。 7、何为平面应力问题和平面应变问题？平面应力问题：在结构上满足a几何条件：研究对象是等厚度薄板。b载荷条件：作用于薄板上的载荷平行于板平面且沿厚度方向均匀分布，而在两板面无外力作用。平面应变问题：满足a几何条件：长柱体，即长度方向的尺寸远远大于横截面的尺寸，且横截面沿长度方向不变。b载荷条件：作用于长柱体结构上的载荷平行于横截面且沿纵向方向均匀分布，两端面不受力两条件的弹性力学问题。 1、何为结构的离散化？离散化的目的？何为有限元模型？ ①离散化:把连续的结构看成由有限个单元组成的集合体。②目的：建立有限元计算模型③通常把由节点,单元及相应的节点载荷和节点约束构成的模型称为有限元模型2、结构离散化时，划分单元数目的多少以及疏密分布，将直接影响到什么？确定单元数量的原则？通常如何设置节点？