(完整word版)有限元法课后习题答案

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. 有限元法有哪些优缺点

答：优点：有限元法可以模拟各种几何形状复杂的结构，得出其近似解；通过计算机程序，可以广泛地应用于各种场合；可以从其他CAD软件中导入建好的模型；数学处理比较方便，对复杂形状的结构也能适用；有限元法和优化设计方法相结合，以便发挥各自的优点。

缺点：有限元计算，尤其是复杂问题的分析计算，所耗费的计算时间、内存和磁盘空间等计算资源是相当惊人的。对无限求解域问题没有较好的处理办法。尽管现有的有限元软件多数使用了网络自适应技术，但在具体应用时，采用什么类型的单元、多大的网络密度等都要完全依赖适用者的经验。

5. 梁单元和平面钢架结构单元的自由度由什么确定

答：由每个节点位移分量的总和确定

6. 简述单元刚度矩阵的性质和矩阵元素的物理意义

答：单元刚度矩阵是描述单元节点力和节点位移之间关系的矩阵

单元刚度矩阵中元素aml的物理意义为单元第L个节点位移分量等于1，其他节点位移分量

等于0时，对应的第m个节点力分量。

7. 有限元法基本方程中的每一项的意义是什么P14

答：Q——整个结构的节点载荷列阵（外载荷、约束力）；整个结构的节点位移列阵；结构的整体刚度矩阵，又称总刚度矩阵。

8. 位移边界条件和载荷边界条件的意义是什么

答：由于刚度矩阵的线性相关性不能得到解，引入边界条件，使整体刚度矩阵求的唯一解。

9. 简述整体刚度矩阵的性质和特点P14

答：对称性；奇异性；稀疏性；对角线上的元素恒为正。

10 简述整体坐标的概念P25

答：在整体结构上建立的坐标系叫做整体坐标，又叫做统一坐标系。

11. 简述平面钢架问题有限元法的基本过程

答：1）力学模型的确定，2）结构的离散化，3）计算载荷的等效节点力，4）计算各单元的刚度矩阵，5）组集整体刚度矩阵，6）施加边界约束条件，7）求解降价的有限元基本方程，8）求解单元应力，9）计算结果的输出。

12. 弹性力学的基本假设是什么。

答：连续性假定，弹性假定，均匀性和各向同性假定，小变形假定，无初应力假定。

13.弹性力学和材料力学相比，其研究方法和对象有什么不同。

答：研究对象：材料力学主要研究杆件，如柱体、梁和轴，在拉压、剪切、弯曲和扭转等作用下的应力、形变和位移。弹性力学研究各种形状的弹性体，除杆件外，还研究平面体、空间体，板和壳等。因此，弹性力学的研究对象要广泛得多。研究方法：弹性力学和材料力学既有相似之外，又有一定区别。弹性力学研究问题，在弹性体区域内必须严格考虑静力学、几何学和物理学三方面条件，在边界上严格考虑受力条件或约束条件，由此建立微分方程和边界条件进行求解，得出较精确的解答。而材料力学虽然也考

虑这几方面的条件，但不是十分严格的，材料力学只研究和适用于杆件问题。

14. 简述圣维南原理。

答；把物体一小部分上的面力变换为分布不同但静力等效的面力，但影响近处的应力分量，而不影响远处的应力。“局部影响原理”

15.平面应力问题和平面应变问题的特点和区别各是什么？试各举出一个典型平面应力和平面应变的问题的实例。

答：平面应力问题的特点：长、宽尺寸远大于厚度，沿板面受有平行板的面力，且沿厚度均匀分布，体力平行于板面且不沿厚度变化，在平板的前后表面上无外力作用平面应变问题的特点：Z向尺寸远大于x、y向尺寸，且与z轴垂直的各个横截面尺寸都相同，受有平行于横截面且不沿z向变化的外载荷，约束条件沿z向也不变，即所有内在因素的外来作用都不沿长度变化。区别：平面应力问题中z方向上应力为零，平面应变问题中z方向上应变为零、应力不为零。举例：平面应力问题等厚度薄板状弹性体，受力方向沿板面方向，荷载不沿板的厚度方向变化，且板的表面无荷载作用。平面应变问题——水坝用于很长的等截面四柱体，其上作用的载荷均平行于横截面，且沿柱长方向不变法。

16. 三角形常应变单元的特点是什么？矩形单元的特点是什么？写出它们的位移模式。

答：三角形单元具有适应性强的优点，较容易进行网络划分和逼近边界形状，应用比较灵活。其缺点是它的位移模式是线性函数，单元应力和应变都是常数，精度不够理想。

矩形单元的位移模式是双线性函数，单元的应力、应变式线性变化的，具有精度较高，形状规整，便于实现计算机自动划分等优点，缺点是单元不能适应曲线边界和斜边界，也不能随意改变大小，适用性非常有限。

17. 写出单元刚度矩阵表达式、并说明单元刚度与哪些因素有关。

答：单元刚度矩阵与节点力坐标变换矩阵，局部坐标系下的单元刚度矩阵，节点位移有关的坐标变换矩阵。

18. 如何由单元刚度矩阵组建整体刚度矩阵（叠加法）？

答：（1）把单元刚度矩阵扩展成单元贡献矩阵,把单元刚度矩阵中的子块按其在整体刚度矩阵中的位置排列，空白处用零子块填充。（2）把单元的贡献矩阵的对应列的子块相叠加，即可得出整体刚度矩阵。

19. 整体刚度矩阵的性质。

答：（1）整体刚度矩阵中每一列元素的物理意义为：欲使弹性体的某一节点沿坐标方形发生单位为移，而其他节点都保持为零的变形状态，在各节点上所需要施加的节点力；（2）整体刚度矩阵中的主对角元素总是正的；（3）整体刚度矩阵是一个对称阵；（4）整体刚度矩阵式一个呈带状分布的稀疏性矩阵。（5）整体刚度矩阵式一个奇异阵，在排除刚体位移后，他是正定阵。

20. 简述形函数的概念和性质。

答：形函数的性质有：(1)形函数单元节点上的值，具有“本点为一、他点为零”的性质；（2）在单元的任一节点上，三角函数之和等于1；（3）三角形单元任一一条边上的形函数，仅与该端点节点坐标有关，而与另外一个节点坐标无关；（4）型函数的值在0~1之间变换。21. 结构的网格划分应注意哪些问题.如何对其进行节点编号。才能使半带宽最小。P50，P8相邻节点的号码差最小

答：一般首选三角形单元或等参元。对平直边界可选用矩形单元，也可以同时选用两种或两种以上的单元。一般来说，集中力，集中力偶，分布在和强度的突变点，分布载荷与自由边界的分界点，支撑点都应该取为节点，相邻节点的号码差尽可能最小才能使半带宽最小22. 为了保证解答的收敛性，单元位数模式必须满足什么条件？

答：（1）位移模式必须包含单元刚体位移；（2）位移模式必须包含单元的常应变；（3）位移模式在单元内要连续，且唯一在相邻单元之间要协调。在有限单元法中，把能够满足条件1和条件2的单元称为完备单元，把满足条件3的单元叫做协调单元或保续单元。

23 有限元分析求得的位移解收敛于真实解得下界的条件。

答：1.位移模式必须包含单元的刚体位移，2.位移模式必须包含单元的常应变，3.位移模式在单元内要连续，且位移在相邻单元之间要协调。

24. 简述等参数单元的概念。

答：坐标变换中采用节点参数的个数等于位移模式中节点参数的个数，这种单元称为等参单元。

25. 有限元法中等参数单元的主要优点是什么？

答：1）应用范围广。在平面或空间连续体，杆系结构和板壳问题中都可应用。

2）将不规则的单元变化为规则的单元后，易于构造位移模式。

3）在原结构中可以采用不规则单元，易于适用边界的形状和改变单元的大小。

4）可以灵活的增减节点，容易构造各种过度单元。

5）推导过程具有通用性。一维，二维三维的推导过程基本相同。

26. 简述四节点四边形等参数单元的平面问题分析过程。

答：（1）通过整体坐标系和局部坐标系的映射关系得到四节点四边形等参单元的母单元，并选取单元的唯一模式；（2）通过坐标变换和等参元确定平面四节点四边形等参数单元的几何形状和位移模式；（3）将四节点四边形等参数单元的位移模式代入平面问题的几何方程，得到单元应变分量的计算式，再将单元应变代入平面问题的物理方程，得到平面四节点等参数单元的应力矩阵（4）用虚功原理球的单元刚度矩

阵，，最后用高斯积分法计算完成。

27. 为什么等参数单元要采用自然坐标来表示形函数？为什么要引入雅可比矩阵？

答：简化计算得到形函数的偏导关系。

28．ANSYS软件主要包括哪些部分？各部分的作用是什么？

答：1.前处理模块：提供了一个强大的实体建模及网络划分工具，用户可以方便地构造有限元模型。2.分析计算模块：包括结构分析、流体力学分析、磁场分析、声场分析、压电分析以及多种物理场的耦合分析，可以模拟多种物理介质的相互作用，具有灵敏度分析及优化分析能力。3.后处理模块：可将计算后果以彩色等值线显示、梯度显示、矢量显示、粒子流迹显示、立体切片显示、透明及半透明显示等图形方式显示出来，也可将计算结果以图表、曲线形式显示出来或输出。

29．ANSYS软件提供的分析类型有哪些？

答：结构静力分析、机构动力分析、结构非线性分析、动力学分析、热分析、流体力学分析、电磁场分析、声场分析、压电分析。

30．简述ANSYS软件分析静力学问题的基本流程。

答：1.前处理器：1）定义单元类型，2）定义实常数，3）定义材料属性，4）创建实体几何模型，5）划分网络；

2.求解器：1）定义分析类型，2）施加载荷和位移约束条件，3）求解；

三角形三节点单元的位移是连续的，应变和应力在单元内是常数，因而其相邻单元将具有不同的应力和应变，即在单元的公共边界上和应变的值将会有突变。

矩形单元的边界上，位移是线性变化的，显然，在两个相邻矩形单元的公共边界上，其位移是连续的。

节点的选用原则：一般说，集中力、集中力偶、分布载荷强度的突变点、分布载荷与自由边界的分界点、支承点都能赢取为节点。

单元的划分原则：（1）划分单元的数目，视要求的计算精度和计算机的性能而定。（2）单元的大小，可根据部位的不同而有所不同。

1、试述街节点力和节点载荷的区别。

节点力是单元与节点之间的作用力；如果取整个结构为研究对象，节点力为内力，节点载荷是作用在节点上的外载荷。

2、试述求整体刚度矩阵的两种方法。

分别建立各节点的平衡方程式，写成矩阵形式，可求得整体刚度矩阵；将各单元刚度矩阵按规律叠加，也可得整体刚度矩阵。

3、平面问题中划分单元的数目是否越多越好?

不是越多越好。划分单元的数目，视要求的计算精度和计算机的性能而定。随着单元数目的接连多，有限元解逐步逼近于真实解，但是，单元数目接连加，刚求解的有限元线性方程组的数目接连多，需要占用更多的计算机内存资源，求解时间接连长，所以，在计算机上进行有限元分析时，还要考虑计算机的性能。单元数过多并不经济。

4、写出单元刚度矩阵的表达式，并说明单元刚度与那些因素有关？

[B]-单元应变矩阵，[D]-弹性矩阵，t-厚度）单元刚度矩阵取决于单元的大小、方向、和弹性常数，而与单元的位置无关，即不随单元或坐标轴的平移而改变。

5、选择多项式为单元的位移模式时，除了要满足单元的完备性和协调性要求，还须考虑什么因素？

还须考虑两个因素：1、所选的位移模式应该与局部坐标系的方位无关，即几何各向同性。

2、多项式位移模式中的项数必须等于或稍大于单元边界上的外节点的自由度数，通常取多项式的项数与单元的外节点的自由度数想等。

有限元分析与应用详细例题

《有限元分析与应用》详细例题 试题1：图示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常单元和六节点三角形单元对坝体进行有限元分析，并对以下几种计算方案进行比 较： 1）分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算； 2）分别采用不同数量的三节点常应变单元计算； 3）当选常应变三角单元时，分别采用不同划分方案计算。 一．问题描述及数学建模 无限长的刚性地基上的三角形大坝受齐顶的水压作用可看作一个平面问题，简化为平面三角形受力问题，把无限长的地基看着平面三角形的底边受固定支座约束的作用，受力面的受力简化为受均布载荷的作用。 二．建模及计算过程 1. 分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算 下面简述三节点常应变单元有限元建模过程（其他类型的建模过程类似）： 1.1进入ANSYS 【开始】→【程序】→ANSYS 10.0→ANSYS Product Launcher →change the working directory →Job Name: shiti1→Run 1.2设置计算类型 ANSYS Main Menu: Preferences →select Structural →OK 1.3选择单元类型 单元是三节点常应变单元，可以用4节点退化表示。 ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 4 node 42 →OK (back to Element Types window)→Options… →select K3: Plane Strain→OK→Close (the Element Type window) 1.4定义材料参数

华科大有限元分析题及大作业题答案——船海专业(DOC)

姓名：学号：班级：

有限元分析及应用作业报告 一、问题描述 图示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析，并对以下几种计算方案进行比较： 1)分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算； 2)分别采用不同数量的三节点常应变单元计算； 3)当选常应变三角单元时，分别采用不同划分方案计算。

二、几何建模与分析 图1-2力学模型 由于大坝长度>>横截面尺寸，且横截面沿长度方向保持不变，因此可将大坝看作无限长的实体模型，满足平面应变问题的几何条件；对截面进行受力分析，作用于大坝上的载荷平行于横截面且沿纵向方向均匀分布，两端面不受力，满足平面应变问题的载荷条件。因此该问题属于平面应变问题，大坝所受的载荷为面载荷，分布情况及方向如图1-2所示，建立几何模型，进行求解。 假设大坝的材料为钢，则其材料参数：弹性模量E=2.1e11，泊松比σ=0.3 三、第1问的有限元建模 本题将分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算。 1）设置计算类型：两者因几何条件和载荷条件均满足平面应变问题，故均取Preferences为Structural 2）选择单元类型：三节点常应变单元选择的类型是PLANE42（Quad 4node42），该单元属于是四节点单元类型，在网格划分时可以对节点数目控制使其蜕化为三节点单元；六节点三角形单元选择的类型是PLANE183（Quad 8node183），该单元属于是八节点单元类型，在网格划分时可以对节点数目控制使其蜕化为六节点单元。因研究的问题为平面应变问题，故对Element behavior（K3）设置为plane strain。 3）定义材料参数 4）生成几何模 a. 生成特征点 b.生成坝体截面 5）网格化分：划分网格时，拾取所有线段设定input NDIV 为10，选择网格划分方式为Tri+Mapped，最后得到200个单元。 6)模型施加约束： 约束采用的是对底面BC全约束。 大坝所受载荷形式为Pressure，作用在AB面上，分析时施加在L AB上，方向水平向右，载荷大小沿L AB由小到大均匀分布（见图1-2）。以B为坐标原点，BA方向为纵轴y，则沿着y方向的受力大小可表示为： ρ（1） = gh P- =ρ g = - 10 {* } 98000 98000 (Y ) y

北京科技大学有限元试题及答案

一 判断题（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、平面刚架有限元分析，节点位移有轴向位移、横向位移、转角。 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.有限元法有哪些优缺点 答： 优点：

有限元复习题答案

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

一维有限元法

实习三、一维问题的有限元方法 一）实习问题： 设 ''1 4(0,1) (0)0,(1)x u u xe x u u e e -?-+=-∈??==-??， ~ 1()u x e e u -=--令 将原问题的边界条件齐次化 ''~~ 1 ~~4()(0,1) (0)0,(1)0 x xe x e e x u u u u -?-+=---∈??? ?==， 二）算法描述： 1 ()21,1 01101 ()1,011 101 (),1 011 10(),111 [ ()()()]1[()()()()]1[()()()()]1 [ ()(i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i a p h h q h N d h a p h h q h N N d h a p h h q h N N d h a p h h q h h x x x x x x x x ξξξξξξξξξξξξξξξ------------=+++=- +++=- +++=+++???1 2010 )()()]N N d ξξξξ?1 ()2 1,1 0110 1[ ()()()]i i i i i i i i i a p h h q h N d h x x ξξξξ----=+++? 1 ()1 0101 ()110 ()()()()i i i i i i i i i i b h f h N d b h f h N d x x ξξξ ξξξ ---=+=+?? 1，单元剖分 (1,2,,)i i n e =L 2，i=1 ~ ~ 00A b == 3,计算数值积分：()()()()()()1,11,,1,1,,,,,i i i i i i i i i i i i i i i i a a a a b b -----即得单元上的i i A b 4，将i i A b 迭加到总的~~ A b 中 5，若i<=n,则i=i+1并转到底三步；否则继续下一步 6，根据边界条件调整~ ~ A b （掐头去尾），即得 A 和b 7，解线性方程组Au=b,得u 从而的h u

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

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

有限元分析大作业试题

有限元分析习题及大作业试题 要求：1）个人按上机指南步骤至少选择习题中3个习题独立完成，并将计算结果上交； 2）以小组为单位完成有限元分析计算； 3）以小组为单位编写计算分析报告； 4）计算分析报告应包括以下部分： A、问题描述及数学建模； B、有限元建模（单元选择、结点布置及规模、网格划分方 案、载荷及边界条件处理、求解控制） C、计算结果及结果分析（位移分析、应力分析、正确性分 析评判） D、多方案计算比较（结点规模增减对精度的影响分析、单 元改变对精度的影响分析、不同网格划分方案对结果的 影响分析等） E、建议与体会 4）11月1日前必须完成，并递交计算分析报告（报告要求打印）。

习题及上机指南：（试题见上机指南） 例题1 坝体的有限元建模与受力分析 例题2 平板的有限元建模与变形分析 例题1：平板的有限元建模与变形分析 计算分析模型如图1-1 所示, 习题文件名: plane 0.5 m 0.5 m 0.5 m 0.5 m 板承受均布载荷：1.0e 5 P a 图1－1 受均布载荷作用的平板计算分析模型 1.1 进入ANSYS 程序 →ANSYSED 6.1 →Interactive →change the working directory into yours →input Initial jobname: plane →Run 1.2设置计算类型 ANSYS Main Menu : Preferences →select Structural → OK 1.3选择单元类型 ANSYS Main Menu : Preprocessor →Element T ype →Add/Edit/Delete →Add →select Solid Quad 4node 42 →OK (back to Element T ypes window) → Options… →select K3: Plane stress w/thk →OK →Close (the Element T ype window) 1.4定义材料参数 ANSYS Main Menu : Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:2.1e11, PRXY :0.3 → OK 1.5定义实常数 ANSYS Main Menu: Preprocessor →Real Constant s… →Add … →select T ype 1→ OK →input THK:1 →OK →Close (the Real Constants Window)

有限元试题总结

一、简答题（40分，每小题5分） 1、 分别写出板弯类单元和平面应力膜单元上一个有限元节点的位移自由度 及其相对应的节点力列阵？ (1)薄板弯曲问题单元每节点三自由度，即每个结点有三个位移分量： 挠度w ,绕x 、y 轴转角 ??? ??? ?y x y x w θ θ轴转角绕轴转角绕挠度，即结点i 的位移 {}i yi xi i i x w y w w w d ?? ??? ? ??????????-??=??????????????=θθ ()4,1K =i 同理,相应的结点力 {}）轴力偶（上节中的绕）轴力偶（上节中的 绕竖向力 x y M y M x ??? ???????????=yi xi i i M M f F (2)平面应力膜单元每个节点两自由度，{},T i i u v ,对应节点力{},T xi yi f f 2、 欲求解在ay by cx R '''++=约束下的泛函(;,)b a I F x y y dx '=?极值，新泛函应 如何构造？ 答：* {(;,)()}b a I F x y y ay by cx R dx λ''''=+++-? 3、 欲求解在()(),,R P x y dx Q x y dy =+??约束下的泛函(;,)b a I F x y y dx '=?极 值，新泛函应如何构造？ 答：()()* {(;,)[,,']}b a I F x y y P x y Q x y y R dx λ'=++-? 4、 满足()f g g f ds L ''+=??条件下的泛函(;,)b a I F x y y dx '=?极值求解应如何构造新泛函？

有限单元法部分课后题答案

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 构造单元形函数有哪些基本原则？ 形函数是定义于单元内坐标的连续函数。单元位移函数通常采用多项式，其中的待定常数应该与单元节点自由度数相等。为满足完备性要求，位移函数中必须包括常函数和一次式，即完全一次多项式。多项式的选取应由低阶到高阶，尽量选择完全多项式以提高单元的精度。若由于项数限制而不能选取完全多项式时，也应使完全多项式具有坐标的对称性，并且一

有限元分析案例

有限元分析案例 图1 钢铸件及其砂模的横截面尺寸 砂模的热物理性能如下表所示: 铸钢的热物理性能如下表所示: 一、初始条件:铸钢的温度为2875o F，砂模的温度为80o F；砂模外边界的对流边界条件：对流系数0.014Btu/hr.in2.o F，空气温度80o F；求3个小时后铸钢及砂模的温度分布。 二、菜单操作: 1.Utility Menu>File>Change Title, 输入Casting Solidification; 2.定义单元类型:Main Menu>Preprocessor>Element Type>Add/Edit/Delete, Add, Quad 4node 55； 3.定义砂模热性能:Main Menu>Preprocessor>Material Props>Isotropic,默认材料编号1, 在Density(DENS)框中输入0.054,在Thermal conductivity (KXX)框中输入0.025,在S pecific heat(C)框中输入0.28; 4.定义铸钢热性能温度表:Main Menu>Preprocessor>Material Props>-Temp Dependent->Temp Table,输入T1=0,T2=2643, T3=2750, T4=2875; 5.定义铸钢热性能:Main Menu>Preprocessor>Material Props>-Temp Dependent ->Prop Table, 选择Th Conductivity,选择KXX, 输入材料编号2,输入C1=1.44, C2=1.54, C3=1.22, C4=1.22,选择Apply,选择Enthalpy,输入C1=0, C2=128.1, C3=163.8, C4=174.2; 6.创建关键点:Main Menu>Preprocessor>-Modeling->Create>Keypoints>In Active

深入浅出的讲清楚有限元法

“有限元法基础及应用”补充讲义（一） 顾克秋 （2005年3月） 一、引子——弹簧单元与弹簧系统 目标：掌握离散结构直接刚度法分析的原理和形式。了解有限元位移法列式的形 式和基本概念。 1、典型弹簧单元分析 写成矩阵符号形式： ? =k F j i i j j ku ku u u k F f +-=-==)(? ?? ?????????--=??????j i j i u u k k k k f f 写成矩阵形式： kd f =（1-1） （1-2） （1-3） 1-2

式（1-2）、（1-3）为弹簧单元的刚度方程，反映了单元特性：节点力与节点位移之间的关系。式中： （注意：单元节点力是节点对单元的作用力） f d k ——单元节点力列阵 ——单元节点位移列阵 ——弹簧单元的刚度矩阵 弹簧单元刚度方程讨论： 1） 有何特点？ 对称、奇异、主对角元素恒正 2） 中元素代表什么含义？ 刚度系数大小等于弹簧刚度；每列元素代表一端固定、另一端产生单位 位移时加在弹簧单元上的节点力。 3）上面单元刚度方程可以求解吗？为什么？ 不可以。刚度方程仅仅表征一个典型单元的弹性特性，单元水平上无法确定单元节点位移。只有把系统中所有单元特性集成后，在系统水平上才可能求出所有未知位移和反力。单元水平上，若已知单元的节点位移，可由刚 度方程求出所有单元节点力分量。若节点力已知，单元节点位移不能确定， 单元可作刚体运动（小位移） 。这也是单元刚度矩阵奇异性的物理解释。 k k 2、弹簧系统整体分析原理

以右图的一个弹簧系统为例，研究如何由单元特性集成系统特性并建立对系统进行求解的控制方程。 由前面得到的弹簧单元的刚度方程公式（1-2），分别写出2个弹簧单元的特性方程如下： 单元1 单元2 （注：右端节点力分量的下标1，2为单元节点的局部编号，上标是单元号） 下面按两个方法完成系统特性的装配和控制方程的建立。并在特定条件下求解。 1）由节点平衡方程导出： 系统处于平衡时，考虑各节点（1，2，3节点）的平衡条件： 节点受到的外载荷与节点受到与其连接的所有单元对其作用力（单元节点力 的反作用力）之和等于零。因此有下列（节点）平衡方程（组）： 把单元特性（1-4），（1-5）代入（1-6）得到： 写成矩阵形式： 2 2 32 11 2211 1f F f f F f F =+==（1-4） （1-5） （1-6） 3 2223322211122 1111)(u k u k F u k u k k u k F u k u k F +-=-++-=-=（1-7） （1-8） 图 1-3

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

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 无关，

有限元分析课后习题

3.1“强”形式相关的场变量要求强的连续性。定义这些场变量的所有函数必须可微，而可微的次数必须等于存在于强形式的系统方程中的偏微分方程的次数。“弱”形式通常是积分形式，且对场变量要求较弱的连续性，弱形式通常能得到更精确的解。 3.2 (a) 协调性方程 （b ）本质边界条件或运动边界条件 （c ）在初始刻和末时刻的条件 3.3 （a ）域的离散 （b ）位移插值 （c ）构造形函数 （d ）坐标变换 （e ）整体有限元方程的组装 （f ）位移约束的施加 （g ）求解整体有限元方程 3.4 理论上不用必须离散所求解问题的区域。把问题划分成单元的目的是更容易地假设位移场的模式。 3.5证明： （1）方程的左边为 []2 0120020120 23 012()d ()d [()()()]d 11 ()()()23l l l f x x a a x a x x a a x a x x a l a l a l δδδδδδδδ=++=++=++??? 方程的右边为 2012002301223012()d ()d 11 [] 23 11 [()()()] 23l l f x x a a x a x x a l a l a l a l a l a l δδδδδδ??=++???? =++=++?? 很显然方程的左右两边相等。 （2）方程的左边为 1212d () (2)d ()()2f x a a x x a a x δ δδδ=+=+ 方程的右边为 []201212d d ()()d d ()()2f x a a x a x x x a a x δδδδδδ=++=+ 很显然方程的左右两边相等。 3.6再生性和连续性

ansys有限元分析作业经典案例

有 限 元 分 析 作 业 作业名称 输气管道有限元建模分析 姓 名 陈腾飞 学 号 3070611062 班 级 07机制（2）班 宁波理工学院

题目描述： 输气管道的有限元建模与分析 计算分析模型如图1所示 承受内压：1.0e8 Pa R1=0.3 R2=0.5 管道材料参数：弹性模量E=200Gpa；泊松比v=0.26。 图1受均匀内压的输气管道计算分析模型（截面图） 题目分析： 由于管道沿长度方向的尺寸远远大于管道的直径，在计算过程中忽略管道的断面效应，认为在其方向上无应变产生。然后根据结构的对称性，只要分析其中1/4即可。此外，需注意分析过程中的单位统一。 操作步骤 1.定义工作文件名和工作标题 1.定义工作文件名。执行Utility Menu-File→Chang Jobname-3070611062，单击OK按钮。 2.定义工作标题。执行Utility Menu-File→Change Tile-chentengfei3070611062，单击OK按钮。 3.更改目录。执行Utility Menu-File→change the working directory –D/chen 2.定义单元类型和材料属性 1.设置计算类型 ANSYS Main Menu: Preferences →select Structural →OK

2.选择单元类型。执行ANSYS Main Menu→Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 8node 82 →apply Add/Edit/Delete →Add →select Solid Brick 8node 185 →OK Options…→select K3: Plane strain →OK→Close如图2所示，选择OK接受单元类型并关闭对话框。 图2 3.设置材料属性。执行Main Menu→Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic，在EX框中输入2e11，在PRXY框中输入0.26，如图3所示，选择OK并关闭对话框。 图3 3.创建几何模型 1. 选择ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS →依次输入四个点的坐标：input:1(0.3,0),2(0.5,0),3(0,0.5),4(0,0.3) →OK

有限元分析及其应用思考题附答案2012

有限元分析及其应用-2010 思考题： 1、有限元法的基本思想是什么？有限元法的基本步骤有那些？其中“离散”的含义是什 么？是如何将无限自由度问题转化为有限自由度问题的？ 答：基本思想：几何离散和分片插值。 基本步骤：结构离散、单元分析和整体分析。 离散的含义：用假想的线或面将连续物体分割成由有限个单元组成的集合，且单元之间仅在节点处连接，单元之间的作用仅由节点传递。当单元趋近无限小，节点无限多，则这种离散结构将趋近于实际的连续结构。 2、有限元法与经典的差分法、里兹法有何区别？ 区别：差分法：均匀离散求解域，差分代替微分，要求规则边界，几何形状复杂精度较低； 里兹法：根据描述问题的微分方程和相应的定解构造等价的泛函表达式，求得近似解； 有限元：基于变分法，采用分片近似进而逼近总体的求解微分方程的数值计算方法。 3、一根单位长度重量为q的悬挂直杆，上端固定，下端受垂直向下的外力P，试 1）建立其受拉伸的微分方程及边界条件； 2）构造其泛函形式； 3）基于有限元基本思想和泛函求极值构造其有限元的计算格式（即最小势能原理）。4、以简单实例为对象，分别按虚功原理和变分原理导出有限元法的基本格式（单元刚度矩 阵）。 5、什么是节点力和节点载荷？两者有何区别？ 答：节点力：单元与单元之间通过节点相互作用 节点载荷：作用于节点上的外载 6、单元刚度矩阵和整体刚度矩阵各有何特点？其中每个矩阵元素的物理意义是什么（按自 由度和节点解释）？ 答：单元刚度矩阵：对称性、奇异性、主对角线恒为正 整体刚度矩阵：对称性、奇异性、主对角线恒为正、稀疏性、带状性。 Kij，表示j节点产生单位位移、其他节点位移为零时作用i节点的力，节点力等于节点位移与单元刚度元素乘积之和。 7、单元的形函数具有什么特点？有哪些性质？ 答：形函数的特点：Ni为x，y的坐标函数，与位移函数有相同的阶次。 形函数Ni在i节点的值为1，而在其他节点上的值为0； 单元内任一点的形函数之和恒等于1； 形函数的值在0~1间变化。 8、描述弹性体的基本变量是什么？基本方程有哪些组成？ 答：基本变量：外力、应力、应变、位移 基本方程：平衡方程、几何方程、物理方程、几何条件 9、何谓应力、应变、位移的概念？应力与强度是什么关系？ 答：应力：lim△Q/△A=S △A→0 应变：物体形状的改变 位移：弹性体内质点位置的变化 10、问题的微分方程提法、等效积分提法和泛函变分提法之间有何关系？何谓“强形 式”？何谓“弱形式”，两者有何区别？建立弱形式的关键步骤是什么？

有限元考试试题及答案

一、 简答题（共40分，每题10分） 1. 论述单元划分应遵循的原则。 2. 说明形函数应满足的条件。 3. 说明四边形等参数单元中“等参数”的含义，即为什么要引入等参数单元。 4. 阐述边界元法的主要优缺点。 二、 计算题（共60分，每题20分） 1. 一杆件如图3所示，杆件上方固定后，在下方受垂直向下的集中力作用，已 知：杆件材料的杨氏模量2 721/100.3in lbf E E ?==，截面积2125.5in A =， 2275.3in A =，长度in L L 1221==，集中力lbf P 100=，用有限元方法求解B 点 和C 点位移。备注：（1）1 lbf （磅力，libra force ） = N 。（2）杨氏模量、弹性 模量、Young 氏弹性模量具有相同含义（10分） 2. 如图2 所示，有一正方形薄板，沿对角承受压力作用，厚度t=1m ，载荷 F=20KN/m ，设泊松比μ=0，材料的弹性模量为E ，试求它的应力分布。（15分） 学院 专业 学号 姓名 y 图1

图2 3. 图示结点三角形单元的124边作用有均布侧压力q，单元厚度为t，求单元的等效结点荷载。 图3

一、简答题 1. 答： 1）合理安排单元网格的疏密分布 2）为突出重要部位的单元二次划分 3）划分单元的个数 4）单元形状的合理性 5）不同材料界面处及荷载突变点、支承点的单元划分 6）曲线边界的处理，应尽可能减小几何误差 7）充分利用结构及载荷的对称性，以减少计算量 2. 答： 形函数应满足的三个条件： a.必须能反映单元的刚体位移，就是位移模式应反映与本单元形变无关的由 其它单元形变所引起的位移。 b.能反映单元的常量应变，所谓常量应变，就是与坐标位置无关，单元内所 有点都具有相同的应变。当单元尺寸取小时，则单元中各点的应变趋于相 等，也就是单元的形变趋于均匀，因而常量应变就成为应变的主要部分。 c.尽可能反映位移连续性；尽可能反映单元之间位移的连续性，即相邻单元 位移协调。 3. 答： 含义：所谓的等参数单元，就是在确定单元形状的插值函数和确定单元位移场的插值函数中采用了完全相同的形函数。 意义：构造出一些曲边地高精度单元，以便在给定地精度下，用数目较少地单元，解决工程实际地具体问题。 4. 答： 有限单元法是基于变分原理的里兹（Ritz）法的另一种形式，从而使里兹法分析的所有理论基础都适用子有限单元法，确认了有限单元法是处理连续介质问题的一种普遍方法．利用变分原理建立有限元方程和经典里兹法的主要区别是有限单元法假设的近似函数不是在全求解域而是在单元上规定的，面且事先不要求满足任何边界条件，因此它可以用来处理很复杂的连续介质问题。有限单元法中所利用的主要是伽辽金（Galerkin）法。它可以用于已经知道问题的微分方程和

有限元考试试题及答案第一组

有限元考试试题及答案 一、简答题（5道，共计25分）。 1.有限单元位移法求解弹性力学问题的基本步骤有哪些？（5分） 答：（1）选择适当的单元类型将弹性体离散化； （2）建立单元体的位移插值函数； （3）推导单元刚度矩阵； （4）将单元刚度矩阵组装成整体刚度矩阵； （5）代入边界条件和求解。 2. 在划分网格数相同的情况下，为什么八节点四边形等参数单元精度大于四边形矩形单元？（5分） 答：在对于曲线边界的边界单元，其边界为曲边，八节点四边形等参数单元边上三个节点所确定的抛物线来代替原来的曲线，显然拟合效果比四边形矩形单元的直边好。 3.轴对称单元与平面单元有哪些区别？（5分） 答：轴对称单元是三角形或四边形截面的空间的环形单元，平面单元是三角形或四边形平面单元；轴对称单元内任意一点有四个应变分量，平面单元内任意一点非零独立应变分量有三个。 4.有限元空间问题有哪些特征？（5分） 答：（1）单元为块体形状。常用单元：四面体单元、长方体单元、直边六面体单元、曲边六面体单元、轴对称单元。（2）结点位移3个分量。（3）基本方程比平面问题多。3个平衡方程，6个几何方程，6个物理方程。 5.简述四节点四边形等参数单元的平面问题分析过程。（5）分） 答：（1）通过整体坐标系和局部坐标系的映射关系得到四节点四边形等参单元的母单元， 并选取单元的唯一模式； （2）通过坐标变换和等参元确定平面四节点四边形等参数单元的几何形状和位移模式；

（3）将四节点四边形等参数单元的位移模式代入平面问题的几何方程，得到单元应变 分量的计算式，再将单元应变代入平面问题的物理方程，得到平面四节点等参 数单元的应力矩阵； （4）用虚功原理求得单元刚度矩阵，最后用高斯积分法计算完成。 二、论述题（3道,共计30分）。 1. 简述四节点四边形等参数单元的平面问题分析过程。（10分） 答：（1）通过整体坐标系和局部坐标系的映射关系得到四节点四边形等参单元的母单元，并选取单元的唯一模式； （2） 通过坐标变换和等参元确定平面四节点四边形等参数单元的几何形状和位移模式； （3）将四节点四边形等参数单元的位移模式代入平面问题的几何方程，得到单元应变 分量的计算式，再将单元应变代入平面问题的物理方程，得到平面四节点等参 数单元的应力矩阵； （4）用虚功原理求得单元刚度矩阵，最后用高斯积分法计算完成。 2.轴对称问题的简单三角形单元是否是常应力，常应变？为什么？（10分） 答：不是常应力和常应变。 因为应变与位移分量的关系式为: ? ?????????????? ? ?????????????? ???? =????????????????????????+??????=??? ???????????=w u 010r r u r u }{rz z r r z z r r w z u z w γεεεεθ,这里除含有微分算符外，还包含了r 的倒数项1/r ,则即使位移模式为线性的，但由于该项的存在，使得应变与坐标有关， 即不会是常应变。应力应变的物理关系为{ }[]{}εσD = ，由于应变不是常应变，则所求得的应力也不会是常应力。