1、名词解释：
(大
剪应力互等定律：作用在两个互相垂直的面上并且垂直于该两面交线的剪应力是互等的。
小相等，正负号也相同)。
位移法：位移法是解决超静定结构最基本的计算方法。
虚功原理：弹性力学中的虚功原理可表达为：在外力作用下处于平衡状态的弹性体，如果发生了虚位移，那么所有的外力在虚位移上的虚功(外力功 )等于整个弹性体内应力在虚应
变上的虚功 (内力功 )
圣维南原理：对于作用于物体局部边界上的面力用另一组与之静力等效（主矢和主矩相等）并且作用于同一小块表面上的力系来代替，则在力系作用区域的附近，应力分布将有
显著的改变，但在远处所受的影响可不计。
最小势能原理：在满足位移边界条件的所有可能位移中，实际发生的位移使弹性体的势能最小。即对于稳定平衡状态，实际发生的位移使弹性体总势能取极小值。显然，最小势
能原理与虚功原理完全等价
叠加原理：在线弹性(物理线性 )和小变形（几何线性）情况下，作用于物体上几组荷载产生的应力和变形的总效应，等于每组荷载单独作用效应的总和。
位移模式 : 按弹性力学位移法求近似解的思路，位移作为基本未知量时，需要对单元上位移
的分布作出假设，即构造含待定参量的简单位移函数——位移模式。
形函数 :
等参元 : 等参数单元（简称等参元）就是对单元几何形状和单元内的参变量函数采用相同数目的节点参数和相同的形函数进行变换而设计出的一种新型单元。
节点的自由度：节点所具有的位移分量的数目。
单元自由度 : 一个单元所有节点的自由度总和称为单元自由度。
刚度集成法结构中的结点力是相关单元结点力的叠加，整体刚度矩阵的系数是相关单元的单元刚度矩阵系数的集成。
虚功等效：原单元载荷与等效节点载荷在单元任意虚位移上的虚功相等。
2、填空题
1．有限元法的实质：将复杂的连续体划分为有限个简单的单元体，化无限自由度问题为有限自
由度问题。
将连续场函数的（偏）微分方程的求解问题转化为有限个参数的代数方程组的求解问题。
2．弹性力学的基本假设：（1）(连续性)（2）均匀性（3）(各向同性)（4）完全弹性符合
（1）-（ 4）假定的称为 (理想弹性体 )。（ 5）小变形假定满足以上五个基本假设的弹性力学称
为 (弹性力学 )。
3．有限元方法的两大应用：科学计算或(数值模拟) 、数字设计或 (虚拟仿真 )。
4．弹性体受外力以后，其内部将产生(应力 )。还将发生变形。物体的变形状态，一般有两
种方式来描述： (1、给出各点的位移；（运动）三个位移分量 )，(2、给出各微元体的变形（应变）
六个应变分量 )。
5．在弹性体内部，三类基本方程：根据微分体上(平衡条件 )，建立平衡方程。根据微分线 段上应变－位移的几何条件，建立 (几何方程 )；根据 (应力－应变间的物理条件
)，建立物理
方程。
6．圣维南原理主要作用是
(静力等效 )
1
。
（i
l ， m ，n ） N
N i (a i i i
k
7．平面三角形的单元自由度是 （ 6），形函数是
2
，
矩形单元的单元自由度是 （ 8），形函数是
，
N l
空间四面体的单元自由度是
（ 8），形函数是，
N
1 ( 1
)( 1
i )
(i=1,2,3,4)i , i 为 i 节点的局部坐标。 i
i
4
N m
板壳单元的单元自由度是
，形函数是
。
8．单元刚度矩阵具有如下性质
N n
（对称性）、（奇异性）、（主对角元素恒等） 。
9．离散化的内容包括
（结构离散）、（载荷离散）。
10．单元参数只能通过 （节点） 传递到相邻单元。
11 位移函数中待定常数个数应等于
（单元节点自由度总数） 。
12．形函数决定了 （单元上位移分布的形态） 。事实上， 单元位移模式就是所有形函数的 （线
性组合），一个单元的位移模式决定了 （单元描述局部位移场）
的能力，决定求解的精度、
收敛性等，而形函数是最重要的因素。
13．单元刚度矩阵通式 k e
V e B
T
D
B dV
。
14．单元刚度矩阵中元素的物理意义为 （单元受节点力作用后抵抗变形的能力） ，它决定于
该单元的形状、大小、方位和 （弹性常数） ，而与单元的 （位置） 无关，即不随单元或坐标
轴的平行移动而改变。
15 载荷移置的原则 （虚功等效）。
16 引入约束的方法常有： （对角元素置一法） 主要用于节点固定的场合，对于给定节点位移的场合，主要用 （对角元素乘大数法） 引入约束。
17 等参单元刚度矩阵的积分式中被积函数很难导出解析表达式，因此等参单元的计算都采 用（数值积分） 求积分的近似值，考虑到减少计算点数，多采用 （高斯数值积分） 。
18 轴对称问题只需在 （子午面） 内描述， 板的弯曲问题主要是在 内描述。
19 计算稳态温度场实际上是求解 （偏微分方程的边值问题） ，采用 （加权余量法） 建立稳
态温度场分析的有限元列式。采用
（伽辽金法） 对权函数进行选择。
20 动态分析中包括 （固有特性分析） 和响应分析。固有特性分析主要是求解各级模态，包括求解 （各级固有频率） 和振型。 简答题：
1、有限元法的实质？
将复杂的连续体划分为有限个简单的单元体，化无限自由度问题为有限自由度问题。 将连续场函数的 （偏）微分方程的求解问题转化为有限个参数的代数方程组的求解问题。
2、理想弹性体的五点假设？
（ 1）连续性（ 2）均匀性 （ 3）各向同性 （ 4）完全弹性 符合（ 1）-（ 4）假定的称为 理想弹性体。 （ 5）小变形假定
满足以上五个基本假设的弹性力学称为线弹性力学。
3、有限元分析的基本步骤？ 1）建立研究对象的近似模型。
2）将研究对象分割成有限数量的单元
(结构离散化 )
3）用标准方法对每个单元提出一个近似解
（单元分析）
（ 1）选择位移模式
1
x
y (1
)(1
)
4
a
b 1
(1
x
)(1 y )
4
a b
1
(1
x
)(1 y ) 4 a b
1
(1
x
)(1 y ) 4
a b
（ 2）建立单元刚度矩阵
（ 3）计算等效节点力
4）将所有单元按标准方法组合成一个与原有系统近似的系统 （单元集成）
5）用数值方法求解这个近似系统。 （选择合适计算方法计算节点位移、应力、应变）
6）计算结果处理与结果验证
（如何显示、分析数据并找到有用的结论）
4、什么是形函数，形函数的基本性质？
形函数决定了单元上位移分布的形态。 事实上，单元位移模式就是所有形函数的线性组合。
一个单元的位移模式决定了该单元描述局部位移场的能力， 决定求解的精度、 收敛性等，
而形函数是最重要的因素。
5、简述有限元法中选取单元位移函数（多项式）的一般原则。
位移函数中待定系数个数应等于单元节点自由度总数。
位移函数的形式一般选为完全多项式，根据 Pascal 三角形由低阶到高阶按顺序、对称地选取；多项式的项数等于（或稍大于）单元节点自由度数。
6、什么是平面应力问题、平面应变问题？ 平面应力问题的基本特征： 1）几何特征
物体在一个方向（ z)的尺寸远远小于其它两个方向（ x,y) 的尺寸。几何为均匀薄板。
2）受力特征
薄板的两个侧面上无载荷作用
边缘上受到平行于板面且沿板厚均匀分布的面力作用； 体力平行于板面且不沿板厚变化（ x,y 的函数）
平面应变问题的基本特征： 1）几何特征
一个方向（ z)尺寸远远大于其它两个方向（ x,y) 的尺寸 ,呈现为无限长等截面柱体。
2）受力特征
外力（体力、面力）平行于横截面作用，且沿纵向不变化。
通常用多项式函数作位移模式， 对三节点三角形单元， 有 6 个待定节点位移分量， 所以单元上的位移函数只能是含 6 个待定系数的完全一次多项式：
u( x, y) a 1 a 2 x a 3 y
1 ~ a
v( x, y) a 4 a 6 为待定系数，称为广义坐标。
a 5 x a 6 y
7、什么是位移模式？如何构建，以三角形三节点为例分析其收敛性。
按弹性力学位移法求近似解的思路， 位移作为基本未知量时， 需要对单元上位移的分布作出假设， 即构造含待定参量的简单位移函数 —— 位移模式。 在有限元中， 当单元划分得足够小时，把位移函数设定为简单的多项式就可以获得相当好的精确度
8、单元刚度矩阵的特点是什么？整体刚度阵的特点？ 单元刚度矩阵的特点 1）物理意义 2）对称性 3）奇异性
4）主对角元素恒正
总刚矩阵具有单元刚阵的性质：对称性、奇异性、主对角元素恒正、奇偶行元素之和分别为零、稀疏性、带状性
9、为什么单元载荷需要移置？
有限元模型中单元通过节点连接形成离散结构；通过节点传递位移和力；单元和整体结构的特性主要是节点力学量之间的关系。因此边界条件必须对节点给出，所有载荷必须等效作用在节点上，这也是连续模型离散化的要求。
10等参元有限元分析三个基本步骤是什么？
1）计算用局部坐标表示的形函数N i对整体坐标x、 y 的偏导数；
2）将整体坐标系中的面积积分转换为在局部坐标系中的面积积分；
3）用数值积分计算出单元刚度矩阵中的元素。
步骤 1）和 2）是等参单元单元分析的关键步骤。
11如何得到整体刚度矩阵？总刚矩阵的叠加规律是什么？
基本方法是刚度集成法，即整体刚度矩阵是单元刚度矩阵的集成。刚体集成法即结构中的结点力
是相关单元结点力的叠加，整体刚度矩阵的系数是相关单元的单元刚度矩阵系数的集成。
总刚矩阵中各个子块是由各个单元刚度矩阵的相应子块叠加而成的。其叠加规律为：设总刚的某子块为Krs
1、当r=s ， Krs 是主对角线的子块，它们是绕节点r 的各个单元刚度矩阵相应对角线子块
的叠加，如主对角线上的子块；
2、当r≠ s ,而r、s被一个单元同时拥有（即rs 相关，为单元的边），则是拥有该单元边
rs 的所有单元刚度矩阵的相应子块的叠加
3、当r≠ s ,且 r、 s 不属于同一单元，即r、s 互不相关，则Krs 为零子块。
根据这些规律，总刚矩阵[K] 的形成可直接利用单刚矩阵的子块叠加即可；
12、什么是杆，什么是梁？平面钢架问题的单元刚度阵怎样求解？
13、薄板弯曲问题与平面应力问题的区别？
力学概念定义的板是指厚度尺寸相对长宽尺寸小很多的平板，且能承受横向或垂直于板面
的载荷。如板不是平板而为曲的 (指一个单元 )，则称为壳问题。如作用于板上的载荷仅为平行
于板面的纵向载荷，则称为平面应力问题；如作用于板上的载荷为垂直于板面的横向载荷，则称
为板的弯扭问题，常简称板的弯曲问题。
计算题：
1、在均质、等厚的三角形单元 ijm 的 ij 边上作用有沿 x 方向按三角形分布的载荷，求移置后的结
点载荷。
解：
X i N i 0
Y i 0
N i
X j N j 0 q x
tdxdy
Y j s 0
N j 0
X m N m 0
Y m
N m
取局部坐标 s ，在 i 点 s=0，在 j 点 s=l ， L 为 ij 边的长度。在 ij 边上，
以局部坐标表示的插值函数为，
N i
1
s N
j
s N
m 0
L
L
载荷为：
q
x s
q
L
L
2
3
s )q s
tds qt( s
s L
X i (1
2 )
L L 2L
3L
L
s q
s
tds
s 3
L
X j
qt
2 0
L
L
3 L
1 q tL
6 1
qtL
3
2、如图所示等边三角形单元，其厚度为 t
，弹性模量为
E ，泊松比
；单元的边长及结点
编号见图中所示。求： （ 1）构建三角形单元的位移模式
（2）求解形函数矩阵
N
3、已知如图 (a）所示的悬臂深梁，在右端面作用着均布拉力，其合力为P。采用如图 (b) 所
示简单网格，设1/ 3
,厚度为 t,试求节点位移。
4、用刚度集成法求图所示结构的整体刚度矩阵。
单元划分图
5、平面桁架由 2 根相同的杆组成（E， A ， L ）。求：
1）节点 2 位移。
2）每根杆应力。
解：
求出每个单元在总体坐标下的刚度矩阵：
单元 1
单元45o l m2 22
k
1T T k T
1
11
1100
T
0101100
1
EA22110000001100 L22001110100011 001100000011
1111
EA1111
2L1111
1111
单元 2
135o l2m2
22
k 2T2T k 2T2
1100T
0101100 1
EA22110000001100
L22001110100011 001100000011
1111
EA1111
2L1111
1111
将单元1， 2 的刚度方程扩张到系统规模（ 6 阶），相加后引入节点平衡条件：
再引入边界约束和载荷：
则上面 6 阶有限元方程凝聚为：
EA20u 2P1
2 L02v 2P2
解出未知位移得：
u v 2
2
L P 1
EA P 2
按公式计算杆应力：得：
1 2E2
11
L
L
11
2EA
E21111L
L2EA
02
( P1 P2 )
P12 A
P2
P1
P2
2
( P1 P2 )
0 2 A
1有限元是近似求解_一般连续_场问题的数值方法 2有限元法将连续的求解域离散为若干个子域_,得到有限个单元,单元和单元之间用节点相连 3从选择未知量的角度来看,有限元法分为三类位移法. 力法混合法 4以_节点位移_为基本未知量的求解方法称为位移法. 5以_节点力_为基本未知量的求解方法称为力法. 6一部分以__节点位移__,另一部分以_节点力_为基本未知量的求解方法称为混合法. 7直梁在外力的作用下,横截面的内力有剪力_和_弯矩_两个. 8平面刚架结构在外力的作用下,横截面上的内力有轴力_ 、剪力_和弯矩. 9进行直梁有限元分析,平面刚架单元上每个节点的节点位移为挠度和转角 10平面刚架结构中,已知单元e的坐标变换矩阵[T e ]和在局 部坐标系x’O’y’下的单元刚度矩阵[K’]e ，则单元在真体坐标 系xOy下的单元刚度矩阵为_ [K]e = [T e ] T [K’] e [T e ] 13弹性力学问题的方程个数有15个,未知量的个数有15个. 14弹性力学平面问题的方程个数有8_个,未知量个数有8_个15几何方程是研究__应变___和_位移之间关系的方程 16物理方程是描述_应力_和_应变_关系的方程 17平衡方程反映了_应力__和_位移_之间关系的 18把经过物体内任意一点各个_ 截面上的应力状况叫做__该点_的应力状态 19形函数在单元上节点上的值,具有本点为_1_.它点为零的性质,并且在三角形单元的任一节点上,三个行函数之和为_1_ 20 形函数是_三角形_单元内部坐标的_线性位移_函数,他反映了单元的_位移_状态 21在进行节点编号时,要尽量使用同一单元的相邻节点的狭长的带状尽可能小,以使最大限度地缩小刚度矩阵的带宽,节省存储,提高计算效率. 22三角形单元的位移模式为_线性位移模式_- 23矩形单元的位移模式为__线性位移模式_ 24在选择多项式位移模式的阶次时,要求_所选的位移模式应该与局部坐标系的方位无关的性质为几何_各向同性 25单元刚度矩阵描述了_节点力_和_节点位移之间的关系 26在选择多项式作为单元的位移模式时,多项式阶次的确定,要考虑解答的收敛性,即要满足单元的_完备性和协调性要求27三节点三角形单元内的应力和应变是_常数,四节点矩形单元内的应力和应变是线性_变化的 28在矩形单元的边界上,位移是线性_变化的 29整体刚度是一个呈_ 狭长的带状_分布的稀疏矩阵 30整体刚度[K]是一个奇异阵,在排除刚体位移_后,它正义阵1从选择未知量的角度来看,有限元法可分为三类(力法,位移法,混合法) 2下列哪有限元特点的描述中,哪种说法是错误的(D需要使用于整个结构的插值函数) 3几何方程研究的是(A应变和位移)之间关系的方程式 4物理方程是描述(D应力和应变)关系的方程 5平衡方程研究的是(C应力和位移)之间关系的方程式 6在划分单元时,下列哪种说话是错误的(A一般首选矩形单元) 7下列哪种单元的单元刚度矩阵必须通过积分才能得到(D矩形单元) 8单元的刚度矩阵不取决于下列哪种因素(B单元位置) 9可以证明,在给定载荷的作用下,有限元计算模型的变形与实际结构变形之间的关系为(B前者小于后者) 10ANSYS按功能作用可分为若干个处理器,其中(B求解器)用于施加载荷和边界条件 11下列有关有限元分析法的描述中,哪种说话是错误的(B单元之间通过其边界连接成组合体) 12下列关于等参数单元的描述中,哪些说话是错误的(C将规则单元变换为不规则单元后,易于构造位移模式) 13从选择未知量的角度来看,有限元可以分为三类,混合法的未知量是(C节点力和节点位移) 14下列对有限元特点的描述中,哪种说话是错误的(B对有限元求解域问题没有较好的处理方法) 15在划分单元时,下列哪种说话错误(D自由端不能取为节点) 16对于平面问题,选择单元一般首选(D三角形单元或等参单元) 17下列哪种说法不是形函数的性质(C三角形单元任一条边上的形函数,与三角形的三个节点坐标都有关) 18下列四种假设中,哪种分析不属于分析弹性力学的基本假设(C大变形假设) 19下列四种假设中,哪种不属于分析弹性力学的基本假设(B 有限变形假设) 20下列关于三角形单元说法中哪种是错误的(C在单元的公共边上应力和应变的值是连续的) 21下列关于矩形单元的说法哪项是错误的(D其形函数是线形的) 22应用圣维南原理简化边界条件时,静力等效是指前后的力系的(D主矢量相同,对于同一点的主矩也相同) 24描述同一点的应力状态需要的应力分量是(C6个) 25在选择多项式作为单元的位移模式时.多项式阶次的确定,要考虑解答的收敛性,哪种说法不是单元必须满足的要求(D 对称性)
一．是非题（认为该题正确，在括号中打；该题错误，在括号中打×。）(每小题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
姓名：学号：班级：
有限元分析及应用作业报告 一、问题描述 图示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析，并对以下几种计算方案进行比较： 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
西安交通大学 级研究生课程考试试题 考试（查）科目：有限元方法（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分) （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．针对下图所示的3个三角形元，写出用完整多项式描述的位移模式表达式。 2．如下图所示，求下列情况的带宽: a)4结点四边形元； b)2结点线性杆元。 3．对上题图诸结点制定一种结点编号的方法，使所得带宽更小。图左下角的四边形在两种不同编号方式下，单元的带宽分别是多大 4．下图所示，若单元是2结点线性杆单元，勾画出组装总刚后总刚空间轮廓线。系统的带宽是多大按一右一左重新编号（即6变成3等）后，重复以上运算。
5．设杆件1－2受轴向力作用，截面积为A，长度为L，弹性模量为E，试写出 杆端力F 1，F 2 与杆端位移 2 1 ,u u之间的关系式，并求出杆件的单元刚度矩阵)(] [e k 6．设阶梯形杆件由两个等截面杆件○1与○2所组成，试写出三个结点1、2、3的 结点轴向力F 1，F 2 ，F 3 与结点轴向位移 3 2 1 , ,u u u之间的整体刚度矩阵[K]。 7．在上题的阶梯形杆件中，设结点3为固定端，结点1作用轴向载荷F 1 =P，求各结点的轴向位移和各杆的轴力。
8． 下图所示为平面桁架中的任一单元，y x ,为局部坐标系，x ，y 为总体坐标系，x 轴与x 轴的夹角为θ。 （1） 求在局部坐标系中的单元刚度矩阵 )(][e k （2） 求单元的坐标转换矩阵 [T]； （3） 求在总体坐标系中的单元刚度矩阵 )(][e k 9．如图所示一个直角三角形桁架，已知27/103cm N E ?=，两个直角边长度 cm l 100=，各杆截面面积210cm A =，求整体刚度矩阵[K]。
10．设上题中的桁架的支承情况和载荷情况如下图所示，按有限元素法求出各结点的位移与各杆的力。 11．进行结点编号时，如果把所有固定端处的结点编在最后，那么在引入边界条件时是否会更简便些 12．针对下图所示的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.有限元法有哪些优缺点 答： 优点：
一 判断题（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.实现有限元分析标准化和规范化的载体就是单元
《汽车有限元基础》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．有限单元法首先求出的解是节点位移，单元应力可由它求得，其计算公式为。（用符号表示即可） 8．一个空间块体单元的节点有 3 个节点位移： u，v，w 9．变形体基本变量有位移应变应力基本方程平衡方程物理方程几何方程 10.实现有限元分析标准化和规范化的载体就是单元
1、何为有限元法？其基本思想是什么？ 有限元法是一种基于变分法而发展起来的求解微分方程的数值计算方法，该方法以计算机为手段，采用分片近似，进而逼近整体的研究思想求解物理问题。 基本思想是化整为零集零为整。 2、为什么说有限元法是近似的方法，体现在哪里？ 有两点：用离散单元的组合体来逼近原始结构，体现了几何上的近似；而用近似函数逼近未知变量在单元内的真实解，体现了数学上的近似。 3、单元、节点的概念？ 节点：表达实际结构几何对象之间相互连接方式的概念 单元：网格划分中的每一个小部分称为单元，网格间相互联结点称为节点 4、有限元法分析过程可归纳为几个步骤？ 结构离散化、单元分析、整体分析 5、有限元方法分几种？本课程讲授的是哪一种？ 位移法、力法、混合法本课程讲授位移法 6、弹性力学的基本变量是什么？何为几何方程、物理方程及虚功方程？弹性矩阵的特点？ 弹性力学变量:外力、应力、应变和位移。 描述弹性体应变分量与位移分量之间的方程称为几何方程；物理方程描述应力分量与应变分量之间的关系；弹性体上外力在虚位移发生过程中所做的虚功与储存在弹性体内的需应变能相等。 弹性矩阵由材料的弹性模量和泊松比确定，与坐标位置无关。 7、何为平面应力问题和平面应变问题？ 平面应力问题：在结构上满足a几何条件：研究对象是等厚度薄板。b载荷条件：作用于薄板上的载荷平行于板平面且沿厚度方向均匀分布，而在两板面无外力作用。 平面应变问题：满足a几何条件：长柱体，即长度方向的尺寸远远大于横截面的尺寸，且横截面沿长度方向不变。b载荷条件：作用于长柱体结构上的载荷平行于横截面且沿纵向方向均匀分布，两端面不受力两条件的弹性力学问题。 1、何为结构的离散化？离散化的目的？何为有限元模型？ ①离散化:把连续的结构看成由有限个单元组成的集合体。②目的：建立有限元计算模型③通常把由节点,单元及相应的节点载荷和节点约束构成的模型称为有限元模型2、结构离散化时，划分单元数目的多少以及疏密分布，将直接影响到什么？确定单元数量的原则？通常如何设置节点？
有限元复习 一、选择题（每题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）个人按上机指南步骤至少选择习题中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)
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. 设上题中的桁架的支承情况与载荷情况如下图所示,按有限元素法求出各结点的位移与各杆的内力。
有限元分析及其应用-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、问题的微分方程提法、等效积分提法和泛函变分提法之间有何关系？何谓“强形 式”？何谓“弱形式”，两者有何区别？建立弱形式的关键步骤是什么？
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 构造单元形函数有哪些基本原则？ 形函数是定义于单元内坐标的连续函数。单元位移函数通常采用多项式，其中的待定常数应该与单元节点自由度数相等。为满足完备性要求，位移函数中必须包括常函数和一次式，即完全一次多项式。多项式的选取应由低阶到高阶，尽量选择完全多项式以提高单元的精度。若由于项数限制而不能选取完全多项式时，也应使完全多项式具有坐标的对称性，并且一