文档库 最新最全的文档下载
当前位置:文档库 › 有限元考试试题及答案第一组

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

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

有限元考试试题及答案

一、简答题(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 = ,由于应变不是常应变,则所求得的应力也不会是常应力。

3.在薄板弯曲理论中做了哪些假设?薄板单元和厚板单元的基本假设有什么不同?(10分)

答:四种假设: 1)变形前的中面法线在变形后仍为弹性曲面的法线。

2)变形前后板的厚度不变。 3)板变形时,中面无伸缩。 4)板内各水平层间互不挤压。

不同点:薄板单元假设横向纤维无挤压,板的中面法线变形后仍保持为直线,该直线垂直于变形后的中面,但是厚板单元的假设考虑横向变形的影响,板的中面法线变形后仍 基本保持为直线,但该直线不再垂直于变形后的中面,法线绕坐标轴的转角不再是挠度的导数,而是独立的变量。

三、计算题(3道,共计45分)。

1.如图所示等腰直角三角形单元,其厚度为t ,弹性模量为E ,泊松比0ν=;单元的边长及结点编号见图中所示。求

(1) 形函数矩阵N (2) 应变矩阵B 和应力矩阵S

(3) 单元刚度矩阵e K (12分)

解:设图1所示的各点坐标为点1(a ,0),点2(a ,a ),点3(0,0)

于是,可得单元的面积为 1

2

A =2a ,及

(1) 形函数矩阵N 为

12122121

(0a a )a

1

(00a )a 1

(a a 0)

a N x y N x y N x y =

+-=++=-+g g ;

[][]123123 N N N ==N I I I N N N (2) 应变矩阵B 和应力矩阵S 分别为

12a 010-a a -a a ????=??????B ,220010a a a 0????=??????B ,32-a 0100a 0-a ????=??????

B ; []12

3=B B B B

12a 00-a a 11-a a 22E ?

?????=??????S ,22000a a 1a 02E ??????=??????S ,32

-a 000a 10-a 2E ??????=??????

S ;[]

[]

123123 ==S D B B B S S S (3) 单元刚度矩阵e K

2.如图所示的四结点矩形单元,求出节点3的位移。设厚度t =1m ,μ=0,E

为常量。(13分) 注:对于四节点矩形单元有:

()()()()()()()()()?

???

???

?

??

?+-=++=-+=--=

ηξηξηξηξ1141

114

1

114

1

114

1

.

14321N N N N →)

4,3,2,1()1)(1(41=++=i N i i i ηηξξ

()[][][][]e

T A

e

k k k k k k k k k k k k k k k k y x t B D B k ?????

?

?

?????==??4443

42

41343332

31242322

21141312

11d d .2

, [][][][][][][]()

()()()())

4,3,2,1,( 3111311a 212123111311218d d d d 21

11

1

=?????

???

?????

?? ??+-+??? ??+-+-+??? ??+-+??? ??+-=

==????--j i b a b b a a b Et B D B abt y x t B D B k j i j i j i j i j i j i j i j i j i j i j i j i j

T

i

j

T

A

i

ij

ηηξξμξξηηηξμξμηηξμξμηξξηημηηξξμη

ξ解:对于四节点矩形单元有:

()()()()()()()()?

???

???

?

???+-=++=-+=--=

ηξηξηξηξ1141

114

1

1141

114

1

4321N N N N →)4,3,2,1()1)(1(41=++=i N i i i

ηηξξ

[][][][]e

T A

e

k k k k k k k k k k k k k k k k y x t B D B k ?????

??

?????==??4443

42

41343332

31242322

21141312

11d d , []{}{}e e e R k =δ,代入边界条件μ1=ν1=μ2=ν2=μ4=ν4=0,将对应的行和列划掉没剩下的

方程为:

[]?

??

???=?????

?P P k --νμ3333, 又)4,3,2,1()1)(1(4

1

=++=i N i i i ηηξξ,且1133==ηξ,,a=1,b=1

所以

所以

[]??????=??????P P k --νμ3333→?

?????=????????????P P E --νμ4114833

解得???

???=?

????

?1158-νμ33E P 3.有一如图3(a)所示的剪力墙,墙顶作用竖向荷载P 。将该剪力墙划分为两个三结点三角形常应力单元,单元和结点编号如图3(b)所示,并将荷载P 分成两个P/2作用在3、4结点。已知单元厚度为t ,弹性模量为E ,泊松比μ=1/3。求结点3和结点4的位移,以及单元①的应变和应力。(20分)

解:建立直角坐标系(注Y 轴向下为正),单元①i,j,m 对应的节点编号为3,1,4,单元②对

应的节点编号为2,4,1。

对于单元①:i(0,0),j(0,4),m(2,0) bi=yi-ym=4;bj=ym-yi=0;bm=yi-yj=-4 ci=xm-xj=2;cj=xi-xm=-2;cm=xj-xi=0 三角形面积A=1/2*2*4=4

几何矩阵[B]=????

???

???----40000204024202000481

弹性矩阵[D]=?????

???

?

??

?

--210

0010112μμ

μ

μE

单元刚度矩阵

?=t B D B k T ]][[][][1=314)1(8001684)1(4)1(8)1(40816048)1(2)1(4)1(2)23(4)1(4)9(2)1(162????

?????

???????????------------+--μμμμμμμμμμμμμμEt

然后合成总刚[K]。

整体节点力矢量为]02

02

[}{2211P P

F F F F F y

x y x =

节点位移矢量为[]44

330000

}{v u v u d =

}]{[}{d K F =,采用缩减矩阵法划去位移为零的行与列,得

解得[]T Et P v u v u 88.199.850.142.8443

3-=????

?

???????

单元①的应变{}[]T Et

P

d B 56.047.0285.0}]{[1

1

-=

=ε 单元①的应力{}[]T t

P

d B D 84.068.158.0}]{][[1

1

-=

有限元填空选择题及答案

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

华科大有限元分析题及大作业题答案——船海专业(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

2011有限元试题

西安交通大学 级研究生课程考试试题 考试(查)科目:有限元方法(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)

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

1.针对下图所示的3个三角形元,写出用完整多项式描述的位移模式表达式。 2.如下图所示,求下列情况的带宽: a)4结点四边形元; b)2结点线性杆元。 3、对上题图诸结点制定一种结点编号的方法,使所得带宽更小。图左下角的四边形在两种不同编号方式下,单元的带宽分别就是多大? 4、下图所示,若单元就是2结点线性杆单元,勾画出组装总刚后总刚空间轮廓线。系统的带宽就是多大?按一右一左重新编号(即6变成3等)后,重复以上运算。

5. 设杆件1-2受轴向力作用,截面积为A,长度为L,弹性模量为E,试写出杆端力F1,F 2与杆端位移21,u u 之间的关系式,并求出杆件的单元刚度矩阵)(][e k 6、设阶梯形杆件由两个等截面杆件\o \a c(○,1)与错误!所组成,试写出三个结点1、2、3的结点轴向力F 1,F 2,F3与结点轴向位移321,,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. 设上题中的桁架的支承情况与载荷情况如下图所示,按有限元素法求出各结点的位移与各杆的内力。

有限元分析及其应用思考题附答案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、问题的微分方程提法、等效积分提法和泛函变分提法之间有何关系?何谓“强形 式”?何谓“弱形式”,两者有何区别?建立弱形式的关键步骤是什么?

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

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

相关文档
相关文档 最新文档