当前位置:文档之家› 有限元讲义 2-弹性力学平面问题有限单元法(2.6四结点四边形等参元,2.7八结点曲线四边形等参元,2.8问题补充)

有限元讲义 2-弹性力学平面问题有限单元法(2.6四结点四边形等参元,2.7八结点曲线四边形等参元,2.8问题补充)

有限元讲义 2-弹性力学平面问题有限单元法(2.6四结点四边形等参元,2.7八结点曲线四边形等参元,2.8问题补充)
有限元讲义 2-弹性力学平面问题有限单元法(2.6四结点四边形等参元,2.7八结点曲线四边形等参元,2.8问题补充)

2.6 四结点四边形单元

(The four-node quadrilateral element)

前面介绍了四结点的矩形单元

其位移函数:

xy y x U 4321αααα+++=

xy y x V

8765αααα+++=

为双线性函数,应力,应变在单元内呈线性变化,比常应力三角形单元精度高。但它对边界要求严格。本节介绍的四结点四边形等参

元,它不但具有较高的精度,而且其网格划分也不受边界的影响。

对任意四边形单元(图见下面)若仍直接采用前面矩形单元的位移函数,在边界上它便不再是线性的(因边界不与x,y 轴一致),这样会使得相邻两单元在公共边界上的位移可能会出现不连续现象(非协调元),而使收敛性受到影响。可以验证,利用坐标变换就能解决这个问题,即可以通过坐标变换将整体坐标中的四边形(图a )变换成在局部坐标系中与四边形方向无关的边长为2的正方形。

正方形四个结点i,j,m,p 按反时钟顺序对应四边形的四个结点i j m p 。

正方形的 1-=η 和 1=η 二条边界,分别对应四边形的i ,j 边界和p,m 边界;ξ=-1和ξ=+1分别对应四边形的i ,p 边界和j ,m 边界。

如果用二组直线等分四边形的四个边界线段,使四边形绘成一个非正交网格,那么该非正交网格在正方形上对应着一个等距离的规则网格(见图a, b )。 当然, 局部坐标上的A 点与整体坐标的A 点对应。

一、四结点四边形等参单元的形函数及坐标变换

由于可以将整体坐标下的四边形单元变换成局部坐标下的正方形单元,对于这种正方形单元,自然仍取形函数为: ξηαηαξαα2321+++=U ξηαηαξαα8765+++=V

引入边界条件,即可得位移函数:

=ijmp

i i U N U

i ijmp

i V N V ∑==

写成矩阵形式:

{}

{}[]{}e

e p i p

i e

d N d N N N N V U f =??

????=??????=0

00

式中形函数: ()()()ηηξξηξi i i N ++=

114

1

, ()p m j i ,,, 按照等参元的定义,我们将坐标变换式亦取为:

p p m m j j i i i ijmp

i x N x N x N x N x N x +++==

p p m m j j i i i

ijmp

i

y N y N y N y N y

N y +++==

∑ ()162--

式中形函数N 与位移函数中的完全一致。

可以验证,利用坐标变换式(2-6-1),可以把整体坐标系中的任意四边形单元(图a )变换成在局部坐标系中与四边形对应的边长为2的正方形。因此可以将上述位移函数和形函数用于任意四边形单元,并将形函数中的ξ,η理解为任意四边形单元的局部坐标。

这样由位移函数可以得到单元各点的位移。在四条边界上分别有ξ=±和η=±1,故边界上的位移呈线性变化,位移的连续性可得到保证。

于是,我们可以理解为:任意四边形单元是从基本的正方形单元变换过来的实际单元。因此又称正方形单元为母体单元,或基本单元。 例题:

为了加深理解,现考察实际单元为矩形单元的坐标变换,在2.4节中,我们定义局部坐标与整体坐标的关系是:

()01

x x a

-=

ξ ()01y y b -=η

式中(x 0 , y 0 )为局部坐标原点。 由上第一式()01

x x a

-=

ξ得: ()()j i i j x x x x x a x ++-=+=2

1

210ξξ

将其重新组合:

()()j i x x x ξξ++-=12

1121 ()()()()()()()()p m j i x x x x ηξηξηξηξ+-++++-++--=

114

1114111411141

对照2.4中的形函数表达式,便知:

p p m m j j i i x N x N x N x N x +++= 自然同理可得: p p m m j j i i y N y N y N y N y +++=

由此知,矩形单元可以看作是四结点四边形单元的特例,自然,它也是等参元。

《有限元法概论》(第二版)P 172 中,是这样解释等参元的基本概念和推导方法的: 图形变换

四结点正方形(母元) 图形变换 四结点四边形(等参元) (ξ-η平面内) ───→ (x,y 平面内)

进行图形变换的关键是进行图形结点坐标之间的变换:

正方形结点坐标 坐标变换 四边形结点坐标 (ξi ,ηi ) ────→ (x,y) i=i,j,m,p i=i,j,m,p

为了实现上述结点坐标之间的变换,可利用母元的形函数,得出(ξ,η)和(x,y)之间的坐标变换式。

图形变换具有如下性质:

1. 母元中的坐标线对应于等参元的直线;

2. 四结点正方形母元对应于四个结点可以任意布置的直边四边形等参元;

3. 变换式(2-6-1)能保证相邻等参元的边界位移彼此协调。

二、几何矩阵[B]

已知单元的应变与结点位移之间的关系是:

{}{}d N N N N x y y x

p i p

i ??

?

?

??????

??????????????=0

00

00

ε ()262-- 形函数矩阵[N]只是局部坐标ξ,η的显函数,为求形函数对整体坐标x,y 的偏导数,必须

用复合函数求导公式:

ξ

ξξ????+

????=??y

y N x x N N i i i ηηη????+

????=??y

y N x x N N i i i ()362-- 或写成: []??

??

???

???????????=???

???

????????????y N x N J N N i i i i ηξ ()a 362--

式中: []???

?

?

?????????????=ηη

ξξx y

x y

x J ()b 362-- 称为雅可比矩阵,而把它的行列式称为雅可比行列式。把式(2-6-1)代入[J]得:

[]???

????

???????????????????????????????????=???????????

?????????=∑∑∑∑p p

m m

j j

i i

p m j i

p m j i

ijmp i i ijmp i

i ijmp i i ijmp i

i

y x y x y x y x N N N N N N N N y N x

N

y N x N J ηη

η

η

ξξξξηηξξ 将形函数 ()()ηηξξi i i N ++=1141 ()p m j i ,,, 代入,分别对ξ,η求偏导,即可得到四结点四边形等参元的雅可比矩阵:

[]??

?

?

??

++++=ξβξβηβηβB A B A J 432141 ()562--

式中常数记为: ∑=

ijmp

i

i i

x A ηξ ∑=ijmp

i

i i y

B ηξ

∑=

ijmp

i

i x

ξβ1 ∑=

ijmp

i

i y

ξβ2

∑=

ijmp

i

i x ηβ3 ∑=ijmp

i i y ηβ

4

该雅可比矩阵的逆:

[]()()??

????++-+-+=

-ηβξβηβξ

βA A B B J

J 1324141 雅可比行列式:

()()()()[]ξβηββηβηβA A A J ++-++=

324116

1

()()()[]ηββξββββββ3421324116

1

B A A B -+-+-=

可以证明,如果四结点四边形的四个内角都小于180°的话,雅可比行列式|J|大于零,

其逆阵[J]-1

是存在的。换句话说,为了使上述等参元能保持较好的精度,整体坐标系下所划分的任意四边形单元必须是凸四边形,即任意内角都不能大于180°。四边形也不能太歪斜,否则会影响其精度。

利用雅可比的逆矩阵,即可求出整体坐标系下形函数的偏导数:

[][]()()),,,(11 1141p m j i i i i i i i i i i J N N J y N x N ??????++=??

?

????????

???????=?????

??

???????????--ξξηηηξηξ ()662-- 求出全部偏导,即代回(2-6-2)右侧,即可得到几何矩阵[B], [B]是ηξ,的函数,即:

[][]??

????

??

?????

?????????????????????????????=??????????????????????

=x N y N x N y N x N y N y N y N y N x N x N x N N x y

y x B p p j j i

i p j i p

j i 0

000000

0 将(2-6-6)代入即可获得[B], [B]是ηξ,的函数。

三、单元刚度矩阵

获得[B]后,便可由单刚的一般表达式:

[][][][]??=dxdy B D B t K T

求出四结点四边形的单元刚度矩阵。 在按上述公式作积分运算时,必须把面积元dxdy 变换成ηξd d ,图a 上的面积元abdc 的面积等于矢量→

ab 与矢量→

ac 的矢量积的模,

即微元

?=ac ab dA

沿ξ轴对应于d ξ的矢量增量是: →

→→

??+??=j d y i d x ab ξξ

ξξ 沿η轴对应于η的矢量增量是: →→→

??+??=j d y

i d x ac ηηηη

式中

j

i , 是坐标x,y 的单位矢,注意到:

0=?=?→

→j j i i 1=?→

→j i 则有:

??? ????+?????? ????+??=→→→→j d y

i d x j d y i d x dA ηηηηξξξξ

ηξηξηξηξd d J j i d d x y y x =????

?

??????-????=→→

因此刚度矩阵的积分式:

[][][]??--=111

1 ηξd td J D B K T ()762--

在计算单元刚度矩阵[K]中元素时,由于被积函数中出现了雅可比行列式,使得它用解

析法很难求其积分,故常采用高斯数值积分法.

四、数值积分

1.一维数值积分

()?b

a d F ξξ

基本思想:构造一个多项目式()ξψ,使在i ξ(i=1,2……n)上有 ()()i i F ξξψ= , 然后用近似函数()ξψ的积分

()?b

a

d ξξψ来近似原被积函数()ξF 的积分()?b

a d F ξξ。

i

ξ称为积分点或取样点,积分点i ξ的数值和位置决定了()ξψ近似()ξF 的程度,亦即决定数值积分的精度。

对于n 个积分点,按照积分点位置的不同选择,通常采用两种不同的数值积分方法,Newton-Cotes 积分和高斯积分方案。

二者方法基本相同,只是前者的积分点ξi 是等间距分布,而后者不是等间距分布。高斯积分的积分点位置由下述方法确定:

① 定义n 次多项式()()()()n P ξξξξξξξ---= 21 ② 由下列条件确定n 个积分点位置

()?

=b

a

i d P 0ξξξ 11,0-=n i

由上二式可见,()ξP

有以下性质:

① 在积分点上()0=i P ξ;

② 多项式()ξP 与1

210 ... ,,,-n ξξξξ在(a, b )域内正交。由此可见n 个积分点的位置 ξi 是在求积域(a, b )内与

1

210 ... ,,,-n ξξξξ正交的n 次多项式()ξP 构成方程

()?

=b

a

i d P 0ξξξ的解。

于是,被积函数()ξF 可由2n-1次多项式()ξψ来近似。 ()()()()∑∑=-=-+=n

i n i i i n i

P F l

11

1

ξξβξξξψ

()ξi l 是n-1阶Lagrange 插值函数 ()()()()()()()n n i i i l ξξξξξξξξξξξξξ------=

121

()?b a

d ξξψ近似()?b

a d F ξξ, 并考虑上面确定n 个积分点位置的约定条件得:

()()?∑=+=b

a

n

i i

i

R F H d F 1

ξξξ

式中 ()?-=b

a n i i d l H ξξ1

H i 称为积分的权系数(加权系数, 权),可见加权系数H i 与被积函数()ξF 无关,只与积分点个数和位置有关。

为便于计算积分点位置ξi 和加权系数H i , 常将上式中的积分限规格化,即令: 1-=a 1=b

由此计算出的i ξ , H i 对应于原积分域(a, b)的关系为:

i b a b a ξ22--+ i H a

b 2

- 对于多重积分,按重积分规则,计算内层时,保持外层为常数,逐层计算即可。 有关数

值积分更详尽的资料可参阅其教科书。

2.等参元计算中数值积分阶次的选择

数值积分中,如何选择积分阶次将直接影响计算精度、工作量,选择不当则有可能使计算失败。选择积分阶次的原则首先是要保证积分的精度能满足所求问题的要求。 如一维问题:

设:插值函数中的多项式阶数为p , 微分算子中导数的最高阶为m ,则有限元得到的近似能量是2(p-m)次多项式。若被积函数是2(p-m)次多项式,应选高斯积分的阶次为n=p-m+1。 三结点三角形单元(线性单元)刚度矩阵中被积函数是常数,故只需一个积分点。

双线性单元 p-m+1=1, 但插值函数将包含有二次项(ξ2 ,ξη,η2

)中的ξη项,所以要达到精确积分应采用2×2阶的高斯积分。有的有限元书中列出了不同插值函数的高斯积分点数n 及相应积分点坐标i ξ和权系数H i 的值,可供编写有限元程序时参考。

五、四边形等参元的荷载等效变换

四边形等参元等效结点荷载的计算,仍然利用局部坐标体系。

1.集中力

设在单元上任意一点C 作用有一集中力{}?

??

???=y x P P P ,根据荷载等效变换的一般公式,其等效结点荷载计算公式是:

{}[]{}P N R T C =

式中,

[]T

C N 是形函数[N]在集中力作用点C 上的值。

2.体积力

设在单元上作用有单位体积力

{}?

??

???=y x W W W ,其等效结点荷载的计算公式是:

1

1

11

{}[]{}T

R N W t J d --=?

?

ξdη

3.分布力

设单元的某一边界上承受的分布表面力是

{}?

??

???=y x q q q ,其等效结点荷载的计算公式是,

当分布力作用在单元ij, mp 边界上时:

ξξ

ξd y x t q N R T ?-??+??=1

1

2

2)()(

}{][}{ 当分布力作用在单元ij, mp 边界上时:

ηη

ηd y

x t q N R T ?-??+??=1

1

22)()(

}{][}{

2.7 八结点曲线四边形等参元

在常规有限元程序的单元库内,目前工程上对平面问题最适用的是八结点等参元。下面作一简要介绍。

单元及结点编号如图a 所示,位移和力列向量仍采用与前面类似的排列方式。为了构造其位移函数,按照前面等参元的概念,可先构造八结点正方形母体单元的位移函数(图b ),并取为双二次插值位移函数:

282726524321ξηαηξαηαξηαξαηαξαα+++++++=U

2162152141321211109ξηαηξαηαξηαξαηαξαα+++++++=V

代入边界条件,得到用结点位移{d}表示的位移场:

??

?

??

??

==∑∑==818

1

i i i i i i V N V U N U ()172-- 式中的形函数:

()()()()

()?

?

?

??

=++--==-+++=

8,7,6,5112

14,3,2,11114

1

2

222j N i N i i j j j i i i i i ηηξξηξξηηηξξηηξξ ()272--

按照等参元的定义,实际单元映射成正方的母体单元的坐标变换式与式2-7-2等同,得:

??

?

??

??

==∑∑==8

18

1

i i i i i i y N y x N x ()372-- 前面我们在推导单元刚度矩阵时,都是假定单元结点编号从左下角开始,逆时针编号。

从有限元理论上说,其编号应该是可以任意的。下面给出一种不同结点编号的例子。

实际单元基本单元

2.8 几个问题的补充

一、变厚问题

各单元取不同t 。

二、不同材料问题

各单元取不同E,μ

三、平面应力与平面应变问题

在前面三角形单元的推导中,我们假定其为平面应力问题。工程中,还有另一类情形─平面应变状态。例如,在对坝体或遂洞等长柱体进行分析时,如果取xoy 坐标平面与其横截平面行,而

Z 轴与其长度方向一致(如图)

那么,由于所考察物体在Z 方向的尺寸很大,且又受到平行于xoy 平面,且不沿长度方向变化的荷载作用,就可认为各个横截面应处于同样的状况,即近似认为Z 方向的位移分量W=0, (位移与Z 无关)

于是,由弹性力学知,在六个应力分量中也仅有三

个独立分量σx ,σy 和τxy, 而()y x z

σσμσ+= 不

独立。

并可得到平面应变问题的物理方程。

???? ??---=y x x E σμμσμε112 , ???

? ??---=x y y

E σμμ

σμε112

2

1112μμμτγ-???

?

?

?

-+

=E xy xy

比较平面应力问题:

()y x x E

μσσε-=1

, ()x y y

E

μσσε-=1

, )1(2μτγ-=E xy xy 得知只需把应力问题中的E 换成μ-1E ,μ换成μμ

-1即得应变问题。所以在这类问题的

程序设计中,通常可以同时求解应力和应变问题,只需设置一开关变量便可以实现。

四、各向异性材料

在弹性矩阵[D]中反应。如正交各向异性时的弹性矩阵[D]为:

有限元2-弹性力学平面问题有限单元法(2.1三角形单元,2.2几个问题的讨论)综述

第2章 弹性力学平面问题有限单元法 2.1 三角形单元(triangular Element) 三角形单元是有限元分析中的常见单元形式之一,它的优点是: ①对边界形状的适应性较好,②单刚形式及其推导比较简单,故首先介绍之。 一、结点位移和结点力列阵 设右图为从某一结构中取出的一典型三角形单元。 在平面应力问题中,单元的每个结点上有沿x 、y 两个方向的力和位移,单元的结点位移列阵规定为: 相应结点力列阵为: (式2-1-1) 二、单元位移函数和形状函数 前已述及,有限单元法是一种近似方法,在单元分析中,首先要求假定(构 造)一组在单元内有定义的位移函数作为近似计算的基础。即以结点位移为已知量,假定一个能表示单元内部(包括边界)任意点位移变化规律的函数。 构造位移函数的方法是:以结点(i,j,m)为定点。以位移(u i ,v i ,…u m v m )为定点上的函数值,利用普通的函数插值法构造出一个单元位移函数。 在平面应力问题中,有u,v 两个方向的位移,若假定单元位移函数是线性的,则可表示成: (,)123 u u x y x y ααα==++ 546(,)v v x y x y ααα==++ (2-1-2)a 式中的6个待定常数α1 ,…, α6 可由已知的6个结点位移分量(3个结点的坐标) {}??? ?? ?????=????? ???? ?????????????=m j i m e d d d d m j j i v u v u v u i {} i i j j m X Y X (2-1-1)Y X Y i e j m m F F F F ?? ?? ???? ???? ??==??????????????????

基于弹性力学理论和有限元法分析应力集中问题的讨论

基于弹性力学理论和有限元法分析应力集中问题的讨论 材料在外形急剧变化的部位,局部应力可以超出名义应力的数倍,对于脆性材料局部过早开始破坏,从而,削弱了构件的强度,降低了构件的承载能力。因此在工程實际中,为了确保构件的安全使用,必须科学合理的分析计算应力集中现象,以便找寻到更好的避免措施。本文首先基于弹性力学理论分析带孔无限宽板的应力分布情况,将对象的受力转化成数学表达,结论应证了应力集中的几个特性。 标签:应力集中系数;有限元分析;无限宽板;弹性力学;Inventor运用;ANSYS 1、应力集中 1.1弹性力学中概念,指物体形状、材料性质不均匀导致的局部应力急剧增高的现象。 1.2应力集中系数 最大局部应力与名义应力的比值称为理论应力集中系数ɑ。可以明确地反应应力集中的程度。 最大局部应力σmax可根据弹性力学理论、有限元法计算得到,也可由实验方法测得;名义应力σn是假设构件的应力集中因素(如孔、缺口、沟槽等)不存在,构件截面上的应力。 2、孔周应力在理想状态下的弹性力学理论分析 2.1定义受单向均匀拉伸荷载的无限宽平板,孔径2α圆孔,建立如图一理想模型。 由于结构的对称性,仅分析图一上半段1/4部分x轴正向的状态: 1)圆孔右顶点单元,即当θ=0,r=α时,代入式(2)解算得σy=3σ; 2)距孔0.2倍孔半径外,即当θ=0,r=1.2α时,代入式(2)解算得σy=2.071σ; 3)距孔1倍孔半径外,即当θ=0,r=2α时,代入式(2)解算得σy=1.221σ; 4)距孔1.5倍孔半径外,即当θ=0,r=2.5α时,代入式(2)解算得σy=1.122σ; 5)距孔2倍孔半径外,即当θ=0,r=3α时,代入式(2)解算得σy=1.074σ;

弹性力学与有限元理论部分考试题2123

弹性力学与有限元(理论部分)考试题 姓名: (90分钟) 一、填空题(25分) 1.弹性力学的基本任务是:。(1分)2.弹性力学的基本假设是:,, ,,。(2.5分)3.应力分量包括:。(2分)应变分量包括:。(2分)位移分量包括:。(2分)4.平衡微分方程反映的是分量和分量的关系,有个方程。(1.5分)几何方程反映的是分量和分量的关系,有个方程。(1.5分)物理方程反映的是分量和分量的关系,有个方程。(1.5分)5.在对受力体进行有限元分析划分网格时,网格划分较密时,优点是:,缺点是:,反之亦然,因此划分适合的网格密度十分重要。(2分) 6.对于平面问题,三角形单元是最简单、最常用的单元,在平面应力问题中单元形状为,在平面应变问题中单元形状为。(2分) 7.在有限元分析中,有六个常用矩阵:D,B,S,k,K,N,它们分别叫做矩阵,矩阵,矩阵,矩阵,矩阵和函数。(6分) 8.刚度矩阵的半带宽B与有关,B=2(d+1)。(1分) 二、简答题(33分) 1.弹性力学与材料力学有何异同?(5分) 2.什么叫单元的位移模式?分别写出三节点三角形单元,四节点矩形单元,六节点三角形单元,八节点矩形单元的位移模式。(10分)

y 3.平面应力问题和平面应变问题的特点(应力、应变)各是什么?(6分) 4.圣维南原理?并举例说明如何应用之(5分) 5.轴对称问题的特点是什么?(3分) 6.如何理解等参元法,简述用等参元法进行空间问题有限元分析的过程(4分) 三.如下图为一个受力体,其划分的单元和节点编号如图1所示,求出半带宽,写出其整体刚度矩阵。(10分) 图1 受力体单元的划分和编号图 四.在单元e 中,三角形单元三个节点分别为i 、j 、m ,,请把图2和图3的力简化到各节点上。(12分) 图2 单元受力图 图3单元受力图 y

弹性力学及有限元法学习总结

弹性力学及有限元法学习总结 摘要:本文就弹性力学的研究对象与方法,弹性力学的基本假设,研究方法,有限元法的基本思想,数学基础,有限元分析的基本步骤进行阐述。 正文:弹性力学是固体力学的一个分支学科,是研究固体材料在外部作用下(外 部作用一般包括:荷载、温度变化以及固体边界约束改变),弹性变形及应力状态的一门学科。 弹性力学的研究对象: 材料力学--研究杆件(如梁、柱和轴)材料力学的拉压、弯曲、剪切、扭转和组合变形等问题。 结构力学--在材料力学基础上研究杆系结构结构力学(如桁架、刚架等)。弹性力学--研究各种形状的弹性体,如杆弹性力学件、平面体、空间体、板壳、薄壁结构等问题。 弹性力学研究方法: 在研究方法上,弹力和材力也有区别:弹力研究方法:在区域V内严格考虑静力学、几何学和物理学三方面条件,建立三套方程; 三套方程在边界s 上考虑受力或约束条件,建立边界条件并在边界条件下求解上边界条件; 边界条件述方程,得出较精确的解答。 弹性力学的基本假设: 1)连续性,假定物体是连续的。连续性因此,各物理量可用连续函数表示。 2)均匀性与各向同性假设假定固体材料是均匀的,并且在各个方向上物理特性相同,也即材料的物理性质在空间分布上是均匀的(或不变的)3)小变形假设假定固体材料在受到外部作用(荷载、温度等)后的位移(或变形)与物体的尺寸相比是很微小的,在研究物体受力后的平衡状态时,物体尺寸及位置的改变可忽略不计,物体位移及形变的二次项可略去不 计,由此得到的弹性力学微分方程将是线性的。 4)完全弹性假设假设固体材料是完全弹性的。 5)无初始应力假设假定外部作用(荷载、温度等)之前,物体处于无应力状态,由弹性力学所求得的应力仅仅是由外部作用(荷载、温度等)所 引起的。 有限元法的基本思想: 有限元是一种结构分析的方法,先把所有系统分解为他们的元件或单元,这些元件的行为已经被充分的了解,再把元件重新组装成原来的系统。及将连续的求解区域离散为一组由有限个单元组成并按一定方式相互连接在一起的单元组

弹性力学及有限元基础复习权威版(最新)

《弹性力学及有限元基础》复习思考题 ★1.对弹性体所做的基本假设? 答:连续性假设;均匀性假设;各向同性假设;弹性假设;小变形假设; ★2.用D'Alember 原理由平衡方程推导运动微分方程? 答:微元体的平衡微分方程的表达式为: 31 112111 2332 122221 23 132333 31 23000f x x x f x x x f x x x σσσσσσσσσ????+++=?????????+++=? ????????+++=? ???? 根据D'Alember 原理,将运动物体看成是静止的,将惯性力22()u t ρ?-?当作体力加到微元体上,由上式 可以直接写出弹性动力学问题的运动微分方程: 23111211 12123232 12222221 2321323333321 23()()() u f x x x t u f x x x t u f x x x t σσσρσσσρσσσρ?????+++=????????????+++=? ???????????+++=?????? ☆3.什么是应力张量? 我们说一点的应力状态是什么涵义? 答:应力张量是一点应力状态的完整描述,它有面元方向和分解方向两个方向性,共有九个分量,由于存在对称性,其独立分量只有六个。应力张量是与坐标选择无关的不变量,但其分量与坐标有关,当已知某坐标系中的九个分量时,其他坐标系中的分量均可由应力转换公式确定。 一点的应力状态是一个具有双重方向性的物理量,其中第一个是面元的方向,用其法矢量ν表示,第二个是作用在该面元上的应力矢量方向,一般用其三个分量来表示。 4.在引出 Cauchy 应力公式时, 我们假设四面体处于平衡状态, 如不处在平衡状态则如何? 答:如果不处在平衡状态,Cauchy 应力公式仍然满足,关系式的成立与是否平衡无关。 5.在什么情况下剪应力互等定律不成立? 答:无论在变形体的内部或者表面上,若存在体力偶时,剪应力互等定律不成立。 6.任意斜截面上的正应变和剪应变的意义是什么? 答:应变张量的三个对角分量x ε、y ε、z ε称为正应变,分别等于坐标轴方向三个线元的单位伸长率,伸长为正,缩短为负。应变张量的三个非对角分量xy ε、yz ε、zx ε称为剪应变,分别等于变形前沿该分量下标所示两坐标方向的、相互正交的线元在变形后的夹角减小量之半。 7.刚性位移,刚性转动,刚体位移,刚体转动有何区别? 答:(1)刚性位移:物体内任意两点间无相对位移;(2)刚性转动:应变张量为0,转动张量不为0;(3)刚体位移:运动分为变形运动和刚体运动,每点都发生相同的位移就叫作刚体位移;(4)刚体转动:用刚性

弹性力学与有限元法分析及实例讲解

弹性力学与有限元法分析 弹性力学是固体力学的一个重要分支,是研究弹性固体在受外力作用、温度改变、边界约束或其他外界因素作用下而发生的应力、形变和位移状态的科学。有限单元法是力学、数学、物理学、计算方法、计算机技术等多种学科综合发展和结合的产物,是随着计算机技术的广泛应用而迅速发展起来的一种数值分析方法。有限元法的基本思想就是化整为零,分散分析,再集零为整。即用结构力学方法求解弹性力学问题,实质是将复杂的连续体划分为有限多个简单的单元体,单元体之间仅仅通过结点相连,实现化无限自由度问题为有限稀有度问题,将连续场函数的(偏)微分方程的求解问题转化为有限个参数的代数方程组的求解问题。 有限元方法经过近半个世纪的发展,目前已经成为各种工程问题特别是结构分析问题的标准分析方法,而有限元软件也已成为现代结构设计中不可缺少的工具。有限元软件是有限元理论通向实际工程应用的桥梁,它的应用极大地提高了力学学科解决自然科学和工程实际问题的能力,进一步促进了有限元方法的发展。ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件,广泛用于机械制造、石油化工、航空航天、汽车交通、土木工程、造船、水利等一般工业及科学研究。 ANSYS软件的组成: (一)前处理模块 该模块为用户提供了一个强大的实体建模及网格划分工具,可以方便的构造有限元模型,软件提高了100种以上的单元类型,用来模拟工程中的各种结构和材料。包括: 1.实体建模:参数化建模,布尔运算及体素库,拖拉、旋转、拷贝、蒙皮、倒角等。 2.自动网格划分,自动进行单元形态、求解精度检查及修正。 3.在集合模型上加载:点加载、分布载荷、体载荷、函数载荷。 4.可扩展的标准梁截面形状库。 (二)分析计算模块 该模块包括结构分析(可进行线性分析、非线性分析和高度非线性分析)、流体动力学分析、电磁场分析、声场分析、压电分析以及多物理场的耦合分析,可模拟多种物理介质的相互作用,具有灵敏度分析及优化分析能力。 (三)后处理模块 将计算结果以彩色等值线、梯度、矢量、粒子流、立体切片、透明及半透明等图形方式显示出来,也可以用图表、曲线形式显示或输出。 由于现在只是对ANSYS工程软件有初步的了解和掌握,所以本次作业仅以(1)结构静力学分析为例,运用ANSYS软件对汽车连杆进行受力分析;(2)

试题及其答案--弹性力学与有限元分析(DOC)

如下图所示三角形薄板,按三结点三角形单元划分后,对于与局部编码ijm 对应的整体编码,以下叙述正确的是( D )。 ① I 单元的整体编码为162 ② II 单元的整体编码为426 ③ II 单元的整体编码为246 ④ III 单元的整体编码为243 ⑤ IV 单元的整体编码为564 A. ①③ B. ②④ C. ①④ D. ③⑤ 一、填空题 1、弹性力学研究弹性体由于受外力作用、边界约束或温度改变等原因而发生的应力、 形变和位移。 2、在弹性力学中规定,线应变以伸长时为正,缩短时为负,与正应力的正负号规定相 适应。 3、在弹性力学中规定,切应变以直角变小时为正,变大时为负,与切应力的正负号规 定相适应。 4、物体受外力以后,其内部将发生内力,它的集度称为应力。与物体的形变和材料强度直接有关的,是应力在其作用截面的法线方向和切线方向的分量,也就是正应力和切应力。应力及其分量的量纲是L -1MT -2。 5、弹性力学的基本假定为连续性、完全弹性、均匀性、各向同性。 6、平面问题分为平面应力问题和平面应变问题。 7、已知一点处的应力分量100=x σMPa ,50=y σMPa ,5010=xy τ MPa ,则主应力 =1σ150MPa ,=2σ0MPa ,=1α6135' 。 8、已知一点处的应力分量, 200=x σMPa ,0=y σMPa ,400-=xy τ MPa ,则主应力=1σ512 MPa ,=2σ-312 MPa ,=1α-37°57′。 9、已知一点处的应力分量,2000-=x σMPa ,1000=y σMPa ,400-=xy τ MPa ,则主应力 =1σ1052 MPa ,=2σ-2052 MPa ,=1α-82°32′。 10、在弹性力学里分析问题,要考虑静力学、几何学和物理学三方面条件,分别建立三 套方程。 11、表示应力分量与体力分量之间关系的方程为平衡微分方程。 12、边界条件表示边界上位移与约束,或应力与面力之间的关系式。分为位移边界条件、 应力边界条件和混合边界条件。 13、按应力求解平面问题时常采用逆解法和半逆解法。 14、有限单元法首先将连续体变换成为离散化结构,然后再用结构力学位移法进行求解。 其具体步骤分为单元分析和整体分析两部分。 15、每个单元的位移一般总是包含着两部分:一部分是由本单元的形变引起的,另一部

弹性力学及有限元试题

弹性力学及有限元试题 (一) 问答题(20分) 1、什么是圣维南原理?举例说明怎样把它应用于工程问题 的简化中。 2、什么叫做一点的应力状态?如何表示一点的应力状态(要 求具体说明或表达)。 3、何谓逆解法和半逆解法?它们的理论依据是什么? 4、什么是平面应力问题?什么是平面应变问题?分别写出弹性力学平面应力问题和平面应变问题的物理方程。 5、要保证有限元方法解答的收敛性,位移模式必须满足那些条 件? (二) (10分) 1.利用坐标变换从直角坐标的平衡方程推导极坐标下平衡方程(无体力)。 2.利用坐标变换从直角坐标下几何方程推导极坐标下几何方程。 (三)已知,其他应力分量为零,求位移场。(10分) (四)设有矩形截面的悬臂粱,在 自由端受有集中荷载F;体力可以不

计。试根据材料力学公式,写出弯应力σx和切应力τxy的表达式,并取挤压应力σy=0,然后证明,这些表达式满足平衡微分方程和相容方程,再说明,这些表达式是否就表示正确的解答(10分)。 (五)设半平面体在直边界上受有集中力偶,单位宽度上力偶矩为M,试求应力分量(10分)。 提示:单位厚度上的力偶矩M的量纲是LMT-2,应力只能是M/ρ2的形式,所以可假设应力函数由:Φ=Φ(φ). (六) 铅直平面内的正方形薄板,边长为2a,四边固定,图5—18,只受重力的作用。设μ=0,试取位移分量的表达式为 用瑞利—里茨法求解(15分)。

(七)试按图示网格求解结点位移,取t =1m,μ= 0(15分)。 (八)用刚度集成法求下图所示结构的整体刚度矩阵K。(10分) 要求:单元刚度矩阵元素用e k形式表示;单元刚度矩阵用e K形式表 ij 示,其中e为单元号。

弹性力学与有限元分析试题答案

最新弹性力学与有限元分析复习题及其答案 一、 填空题 1、弹性力学研究弹性体由于受外力作用、边界约束或温度改变等原因而发生的应力、 形变和位移。 2、在弹性力学中规定,线应变以伸长时为正,缩短时为负,与正应力的正负号规定相 适应。 3、在弹性力学中规定,切应变以直角变小时为正,变大时为负,与切应力的正负号规 定相适应。 4、物体受外力以后,其部将发生力,它的集度称为应力。与物体的形变和材料强度直接有关的,是应力在其作用截面的法线方向和切线方向的分量,也就是正应力和切 应力。应力及其分量的量纲是L -1MT -2。 5、弹性力学的基本假定为连续性、完全弹性、均匀性、各向同性。 6、平面问题分为平面应力问题和平面应变问题。 7、已知一点处的应力分量100=x σMPa ,50=y σMPa ,5010=xy τ MPa ,则主应力 =1σ150MPa ,=2σ0MPa ,=1α6135' 。 8、已知一点处的应力分量, 200=x σMPa ,0=y σMPa ,400-=xy τ MPa ,则主应力 =1σ512 MPa ,=2σ-312 MPa ,=1α-37°57′。 9、已知一点处的应力分量,2000-=x σMPa ,1000=y σMPa ,400-=xy τ MPa ,则主应 力=1σ1052 MPa ,=2σ-2052 MPa ,=1α-82°32′。 10、在弹性力学里分析问题,要考虑静力学、几何学和物理学三方面条件,分别建立三 套方程。 11、表示应力分量与体力分量之间关系的方程为平衡微分方程。 12、边界条件表示边界上位移与约束,或应力与面力之间的关系式。分为位移边界条件、 应力边界条件和混合边界条件。

(绝密试题)弹性力学与有限元分析试题及其答案

2012年度弹性力学与有限元分析复习题及其答案 (绝密试题) 一、填空题 1、弹性力学研究弹性体由于受外力作用、边界约束或温度改变等原因而发生的应力、形变和位移。 2、在弹性力学中规定,线应变以伸长时为正,缩短时为负,与正应力的正负号规定相适应。 3、在弹性力学中规定,切应变以直角变小时为正,变大时为负,与切应力的正负号规定相适应。 4、物体受外力以后,其内部将发生内力,它的集度称为应力。与物体的形变和材料强度直接有关的,是应力在其作用截面的法线方向和切线方向的分量,也就是正应力和切应力。应力及其分量的量纲是L -1MT -2。 5、弹性力学的基本假定为连续性、完全弹性、均匀性、各向同性。 6、平面问题分为平面应力问题和平面应变问题。 7、已知一点处的应力分量100=x σMPa ,50=y σMPa ,5010=xy τ MPa ,则主应力=1σ150MPa ,=2σ0MPa ,=1α6135' 。 8、已知一点处的应力分量, 200=x σMPa ,0=y σMPa ,400-=xy τ MPa , 则主应力=1σ512 MPa ,=2σ-312 MPa ,=1α-37°57′。 9、已知一点处的应力分量,2000-=x σMPa ,1000=y σMPa ,400-=xy τ MPa ,则主应力=1σ1052 MPa ,=2σ-2052 MPa ,=1α-82°32′。 10、在弹性力学里分析问题,要考虑静力学、几何学和物理学三方面条件,分别建立三套方程。 11、表示应力分量与体力分量之间关系的方程为平衡微分方程。 12、边界条件表示边界上位移与约束,或应力与面力之间的关系式。分为位移边界条件、应力边界条件和混合边界条件。 13、按应力求解平面问题时常采用逆解法和半逆解法。 14、有限单元法首先将连续体变换成为离散化结构,然后再用结构力学位移法进行求解。其具体步骤分为单元分析和整体分析两部分。 15、每个单元的位移一般总是包含着两部分:一部分是由本单元的形变引起的,另一部分是由于其他单元发生了形变而连带引起的。 16、每个单元的应变一般总是包含着两部分:一部分是与该单元中各点的位置坐标有关的,是各点不相同的,即所谓变量应变;另一部分是与位置坐标无关的,是各点相同的,即所谓常量应变。 17、为了能从有限单元法得出正确的解答,位移模式必须能反映单元的刚体位移和常量应变,还应当尽可能反映相邻单元的位移连续性。 18、为了使得单元内部的位移保持连续,必须把位移模式取为坐标的单值连续函数,为

弹性力学与有限元分析试题及参考答案

弹性力学与有限元分析试题及参考答案 四、分析计算题 1、试写出无体力情况下平面问题的应力分量存在的必要条件,并考虑下列平面问题的应力分量是否可能在弹性体中存在。 (1)By Ax x +=σ,Dy Cx y +=σ,Fy Ex xy +=τ; (2))(22y x A x +=σ,)(22y x B y +=σ,Cxy xy =τ; 其中,A ,B ,C ,D ,E ,F 为常数。 解:应力分量存在的必要条件是必须满足下列条件:(1)在区域内的平衡微分方程 ????? ??=??+??=??+??0 0x y y x xy y yx x τστσ;(2)在区域内的相容方程()02222=+??? ? ????+??y x y x σσ;(3)在边界上的应力边界条件()()()() ???? ?=+=+s f l m s f m l y s xy y x s yx x τστσ;(4)对于多连体的位移单值条件。 (1)此组应力分量满足相容方程。为了满足平衡微分方程,必须A =-F ,D =-E 。此外还应满足应力边界条件。 (2)为了满足相容方程,其系数必须满足A +B =0;为了满足平衡微分方程,其系数必须满足A =B =-C /2。上两式是矛盾的,因此,此组应力分量不可能存在。 2、已知应力分量312x C Qxy x +-=σ,22 23xy C y -=σ,y x C y C xy 2 332--=τ,体力不计,Q 为常数。试利用平衡微分方程求系数C 1,C 2,C 3。 解:将所给应力分量代入平衡微分方程 ???? ?? ?=??+??=??+??00x y y x xy y yx x τστσ 得 ?? ?=--=--+-0 230 33322322212xy C xy C x C y C x C Qy 即 ()()()?? ?=+=+--0 230 333222231xy C C y C Q x C C 由x ,y 的任意性,得

弹性力学与有限元分析试题及其答案

一、填空题 1、弹性力学研究弹性体由于受外力作用、边界约束或温度改变等原因而发生的应力、形变和位移。 2、在弹性力学中规定,线应变以伸长时为正,缩短时为负,与正应力的正负号规定相适应。 3、在弹性力学中规定,切应变以直角变小时为正,变大时为负,与切应力的正负号规定相适应。 4、物体受外力以后,其内部将发生内力,它的集度称为应力。与物体的形变和材料强度直接有关的,是应力在其作用截面的法线方向和切线方向的分量,也就是正应力和切应力。应力及其分量的量纲是L -1MT -2。 5、弹性力学的基本假定为连续性、完全弹性、均匀性、各向同性。 6、平面问题分为平面应力问题和平面应变问题。 7、已知一点处的应力分量100=x σMPa , 50=y σMPa ,5010=xy τ MPa ,则主应 力=1σ150MPa ,=2σ0MPa , =1α6135' 。 8、已知一点处的应力分量, 200=x σMPa , 0=y σMPa ,400-=xy τ MPa ,则主应力=1σ512 MPa ,=2σ-312 MPa , =1α-37°57′。 9、已知一点处的应力分量, 2000-=x σMPa ,1000=y σMPa , 400-=xy τ MPa ,则主应力=1σ1052 MPa ,=2σ-2052 MPa ,=1α-82°32′。 10、在弹性力学里分析问题,要考虑静力学、几何学和物理学三方面条件,分别 建立三套方程。 11、表示应力分量与体力分量之间关系的方程为平衡微分方程。 12、边界条件表示边界上位移与约束,或应力与面力之间的关系式。分为位移边界条件、应力边界条件和混合边界条件。 13、按应力求解平面问题时常采用逆解法和半逆解法。 14、有限单元法首先将连续体变换成为离散化结构,然后再用结构力学位移法进行求解。其具体步骤分为单元分析和整体分析两部分。 15、每个单元的位移一般总是包含着两部分:一部分是由本单元的形变引起的,另一部分是由于其他单元发生了形变而连带引起的。 16、每个单元的应变一般总是包含着两部分:一部分是与该单元中各点的位置坐标有关的,是各点不相同的,即所谓变量应变;另一部分是与位置坐标无关的,是各点相同的,即所谓常量应变。 17、为了能从有限单元法得出正确的解答,位移模式必须能反映单元的刚体位移和常量应变,还应当尽可能反映相邻单元的位移连续性。 18、为了使得单元内部的位移保持连续,必须把位移模式取为坐标的单值连续函数,为了使得相邻单元的位移保持连续,就不仅要使它们在公共结点处具有相同的位移时,也能在整个公共边界上具有相同的位移。 19、在有限单元法中,单元的形函数N i 在i 结点N i =1;在其他结点N i =0及∑N i =1。 20、为了提高有限单元法分析的精度,一般可以采用两种方法:一是将单元的尺寸减小,以便较好地反映位移和应力变化情况;二是采用包含更高次项的位移模式,使位移和应力的精度提高。 二、判断题(请在正确命题后的括号内打“√”,在错误命题后的括号内打“×”)

弹性力学基础及有限单元法

第一章 1、弹性力学的任务是什么 弹性力学的任务是分析各种结构物或其构件在弹性阶段的应力和位移,校核它们是否具有所需的强度和刚度,并寻求或改进它们的计算方法。 2、弹性力学的基本假设是什么?为什么要采用这些假设? (1)假设物体是连续的——物体内部由连续介质组成,物体中没有空隙,因此物体中的应力、应变、位移等量是连续的.可以用坐标的连续函数表示。实际上,所有的物体均由分子构成,但分子的大小及分子间的距离与物体的尺寸相比是很微小的,故可以不考虑物体内的分个构造。根据这个假设所得的结果与实验结果是符合的。 (2)假设物体是匀质的和各向同性的——物体内部各点与各方向上的介质相同,因此,物体各部分的物理性质是相同的。这样,物体的弹性常数(弹性模量、泊松比)不随位置坐标和方向而变化。钢材由微小结晶体组成,晶体本身是各向异性的、但由于晶体很微小而排列又不规则,按其材料的平均性质,可以认为钢材是各向同性的。木材不是各向同性的。 (3)假设物体是完全弹性的—一物体在外加因家(裁荷、温度变化等)的作用下发生变形,在外加固素去除后,物体完全恢复其原来形状而没有任何剩余变形。同时还假定材料服从胡克定律,即应力与形变成正比。 (4)假设物体的变形是很小的——在载荷或温度变化等的作用下,物体变形而产生的位移,与物体的尺寸相比,是很微小的。在研究物体受力后的平衡状态时,可以不考虑物体尺寸的改变。在研究物体的应变时,可以赂去应变的乘积,因此,在微小形变的情况下弹性理论中的微分方程将是线性的。 (5)假设物体内无初应力一一认为物体是处于自然状态,即在载荷或温度变化等作用之前,物体内部没合应力。也就是说,出弹性理论所求得的应力仅仅是由于载荷或温度变化等所产生的。物体中初应力的性质及数值与物体形成的历史有关。若物体中有韧应力存在,则由弹性理论所求得的应力加上初应力才是物体中的实际应力。 上面基本假设中.假设(4)是属于几何假设,其他假设是属于物理假设。 3、举例说明各向同性的物体和各向异性的物体。 钢材由微小结晶体组成,晶体本身是各向异性的、但由于晶体很微小而排列又不规则,按其材料的平均性质,可以认为钢材是各向同性的。木材是各异性的。 4、弹性力学和材料力学相比,其研究方法和对象有什么区别? P3 弹性力学具体的研究对象主要为梁、校、坝体、无限弹性体等实体结构以及板、壳等受 力体。 在材料力学课程中,基本上只研究所谓杆状构件,也就是长度远大干高度和觅度的构 件。这种构件在拉压、剪切、弯曲、扭转作用下的应力和位移,是材料力学的主要研究内 容。

弹性力学与有限元分析复习题(含答案)

分析计算题 1、试写出无体力情况下平面问题的应力分量存在的必要条件,并考虑下列平面问题的 应力分量是否可能在弹性体中存在。 (1)By Ax x +=σ,Dy Cx y +=σ,Fy Ex xy +=τ; (2))(22y x A x +=σ,)(22y x B y +=σ,Cxy xy =τ; 其中,A ,B ,C ,D ,E ,F 为常数。 解:应力分量存在的必要条件是必须满足下列条件:(1)在区域内的平衡微分方程 ?? ? ? ???=??+??=??+??00x y y x xy y yx x τστσ;(2)在区域内的相容方程()02 222=+??? ? ????+??y x y x σσ;(3)在边界上的应力 边界条件()()()() ?? ?? ?=+=+s f l m s f m l y s xy y x s yx x τστσ;(4)对于多连体的位移单值条件。 (1)此组应力分量满足相容方程。为了满足平衡微分方程,必须A =-F ,D =-E 。此 外还应满足应力边界条件。 (2)为了满足相容方程,其系数必须满足A +B =0;为了满足平衡微分方程,其系数必须满足A =B =-C /2。上两式是矛盾的,因此,此组应力分量不可能存在。 2、已知应力分量312x C Qxy x +-=σ,222 3xy C y -=σ,y x C y C xy 2 332--=τ,体力不计,Q 为常数。试利用平衡微分方程求系数C 1,C 2,C 3。 解:将所给应力分量代入平衡微分方程 ?? ? ? ?? ?=??+??=??+??00x y y x xy y yx x τστσ 得 ? ? ?=--=--+-0230 33322322212xy C xy C x C y C x C Qy 即 ()()()?? ?=+=+--0 230333222231xy C C y C Q x C C 由x ,y 的任意性,得

材料力学弹性力学有限元法的异同--tl

材料力学、弹性力学、有限元法的异同 力学是研究力对物体的效应的一门学科。力对物体的效应有两种:一种是引起物体运动状态的变化,称为外效应;另一种是引起物体的变形,称为内效应。材料力学研究力的内效应,即物体的变形和破坏的规律。材料力学主要研究物体受力后发生的变形、由于变形而产生的内力以及物体由此而产生的失效和控制失效的准则。工程中各种结构或机械都是由许多杆件或零部件组成。这些杆件或零部件统称为构件。工程上构件的几何形状是各种各样的,可分为杆件、板(或壳)、实体。材料力学主要的研究对象是杆状构件。材料力学的任务,就是在分析构件内力和变形的基础上,给出合理的构件计算准则,满足既安全又经济的工程设计要求,并为后续课程如机械设计、结构力学、弹性力学和复合材料力学等提供必要的理论基础。 弹性力学又称弹性理论,是固体力学的一个分支学科。它是研究可变形固体在外部因素(力、温度变化、约束变动等)作用下所产生的应力、应变和位移的经典科学。确定弹性体的各质点应力、应变、和位移的目的就是确定构件设计中的强度和刚度指标,以此来解决实际工程结构中的强度、刚度和稳定性问题。弹性力学具体的研究对象主要为梁、柱、坝体、无限弹性体等实体结构以及板、壳等受力体。 弹性力学的研究内容和目的的任务原则上与材料力学相同,但其学科所研究的对象不同,研究方法也不完全相同。 (1)在材料力学课程中,基本上只研究杆状构件(直杆、小曲率杆),也就是长度远大于高度和宽度的构件。这种构件在拉压、剪切、弯曲、扭转作用下的应力和位移,是材料力学的主要研究内容。弹性力学解决问题的范围比材料力学要大得多。如孔边应力集中、深梁的应力分析等问题用材料力学的理论是无法求解的,而弹性力学则可以解决这类问题。如板和壳以及挡土墙、堤坝、地基等实体结构,则必须以弹性力学为基础,才能进行研究。如果要对于杆状构件进行深入的、较精确的分析,也必须用到弹性力学的知识。同时弹性力学又为进一步研究板、壳等空间结构的强度、振动、稳定性等力学问题提供理论依据,它还是进一步学习塑性力学、断裂力学等其他力学课程的基础

《弹性力学及有限元》教学大纲

《弹性力学及有限元》教学大纲 大纲说明 课程代码:5125004 总学时:40学时(讲课32学时,上机8学时) 总学分:2.5学分 课程类别:必修 适用专业:土木工程专业(本科) 预修要求:高等数学、理论力学、材料力学 课程的性质、目的、任务: 本课程是土木工程专业限选修的一门专业基础课。本课程的教学目的,是使学生在理论力学和材料力学等课程的基础上进一步掌握弹性力学的基本概念、原理和方法,了解弹性力学问题的求解思路、方法和解答,为学习相关专业课程打下初步的弹性力学基础。在此基础上,使学生掌握有限单元法的基本概念、理论、方法,了解和应用ANSYS大型结构分析程序求解简单的弹性力学问题。 课程教学的基本要求: 本课程教学环节主要包括:课堂讲授、习题课、作业、答疑、上机计算、考试。采用课堂授课方式,重点章节安排习题课。课后布置一定量的习题,以便掌握弹性力学与有限单元法的基本概念、原理和方法,用弹性力学的求解方法及大型结构分析有限单元程序求解简单的弹性力学问题。考试采用开卷方式。 大纲的使用说明: 本大纲适用于土木工程本科专业40课时的《弹性力学及有限元》课程. 大纲正文 第一章绪论学时:6学时(讲课6学时) 本章讲授要点:了解弹性力学的研究内容,理解体力、面力、应力、应变和位移等基本概念,熟悉体力、面力、应力、应变、位移等力学量的记号和符号的有关规定,理解弹性力学的基本假定;了解有限单元法的发展,掌握泛函、变分和泛函极值等基本概念;了解加权残值、里兹与伽辽金等方法。 重点:弹性力学中的应力、应变和位移等基本概念;泛函、变分、驻值等基本概念;加权残值、里兹与伽辽金等方法。 难点:应力、应变;泛函、变分、驻值;加权残值法、里兹法与伽辽金法。 第一节弹性力学的内容 第二节弹性力学中的几个基本概念 第三节弹性力学中的基本假定 第四节有限单元法的发展简介 第五节变分原理.泛函.变分.驻值 第六节加权残值法、里兹法与伽辽金法

弹性力学及有限元课程大纲

《弹性力学及有限元》课程大纲课程代码EM316 课程名称中文名:弹性力学及有限元 英文名:Elasticity and Finite Element Method 课程类别专业基础课修读类别必修 学分 2 学时32 开课学期第5学期 开课单位船舶海洋与建筑工程学院土木工程系 适用专业土木工程专业 先修课程《高等数学》、《理论力学》、《材料力学》、《结构力学》 教材及主要参考书教材: 徐芝纶. 弹性力学简明教程(第四版),北京:高等教育出版社,2013年6月。ISBN: 9787040373875 参考书: 1. 王润富.弹性力学简明教程学习指导. 北京:高等教育出版社, 2004. ISBN: 7040130815 2. 吴家龙. 弹性力学(新一版). 北京:高等教育出版社,2001. ISBN: 7560812457. 3. S.Timoshenko &J. N. Goodier. Theory of Elasticity.(Third edition) McGraw-hill Book Co.,1970. ISBN-13: 978-0070647206 4. 丁科,陈月顺. 有限单元法. 北京大学出版社,2006. ISBN: 9787301104354 一课程简介 弹性力学及有限元是土木工程专业必修的一门专业基础课。课程主要研究弹性体受外力作用或温度改变等原因而产生的应力、位移和变形。本课程的教学目的,是使学生在理论力学和材料力学等课程的基础上,进一步掌握弹性力学与有限元的基本概念、基本原理和基本方法,提高分析与计算的能力。使学生掌握有限单元法及其工程适用性,为学生从事与土木工程相关的专业技术工作、科学研究工作等打下坚实的基础。 二本课程所支撑的毕业要求 本课程支撑的毕业要求及比重如下: 序号毕业要求指标点毕业要求指标点具体内容支撑比重 1 毕业要求1.3 具有必备的土木工程专业基础知识及 在复杂土木工程问题中应用能力 65% 2 毕业要求5.2 具有至少应用一种土木工程方面的大 型分析软件能力,并了解工程适用性。 35%

弹性力学同有限元分析的关系

弹性力学同有限元分析的关系 弹性力学:弹性力学也称弹性理论,主要研究弹性体在外力作用或温度变化等外界因素下所产生的应力、应变和位移,从而解决结构或机械设计中所提出的强度和刚度问题。 是固体力学的重要分支,它研究弹性物体在外力和其它外界因素作用下产生的变形和内力。 研究对象:包括杆状构件在内的各种形状的弹性体。 弹性力学基本规律:变形连续规律、应力-应变关系和运动(或平衡)规律,它们有时被称为弹性力学三大基本规律。弹性力学中许多定理、公式和结论等,都可以从三大基本规律推导出来。 弹性力学同材料力学的比较 1、研究内容:基本上没有什么区别。弹性力学也是研究弹性体在外力作用下的平衡和运动,以及由此产生的应力和变形。 2、研究的对象:材料力学基本上只研究杆、梁、柱、轴等杆状构件,即长度远大于宽度和厚度的构件;弹性力学虽然也研究杆状构件,但还研究材料力学无法研究的板与壳及其它实体结构,即两个尺寸远大于第三个尺寸,或三个尺寸相当的构件。 3、研究的方法: 相同点:静力学、几何学与物理学三方面进行研究; 不同点:材料力学:

对构件的整个截面建立分析方程,引用一些截面的变形状况或应力情况的假设,因而得出的结果往往是近似的,不精确。 弹性力学: 对构件采用无限小单元体来建立分析方程的,因而无须引用那些假设,分析的方法比较严密,得出的结论也比较精确。所以,可以用弹性力学的解答来估计材料力学解答的精确程度,并确定它们的适用范围。 从几何形状复杂程度来考虑可以分为: 1)简单形状变形体—材料力学 2)任意形状变形体—弹性力学 任意变形体是有限元方法处理的对象,因而,弹性力学中有关变量和方程的描述是有限元方法的重要基础。 弹性力学的弱点:由于研究对象的变形状态较复杂,处理的方法又较严谨,因而解算问题时,往往需要冗长的数学运算。但为了简化计算,便于数学处理,它仍然保留了材料力学中关于材料性质的假定。

弹性力学及有限元考试复习简答题

1、简述有限单元法常分析的问题。 答:有限单元法是一种用于连续场分析的数值模拟技术,他不仅可以对机械、建筑结构的位移场和应力场进行分析,还可以对电磁学中的电磁场、传热学中的温度场、流体力学中的流体场进行分析。 2、在有限单元法中,位移模式应满足哪些基本条件。 答:1位移函数在单元节点的值应等于节点位移(即单元部是连续的) 2所选位移函数必须保证有限元的解收敛于真实解 3、简述有限单元法结构刚度矩阵的特点。 答:对称矩阵奇异矩阵稀疏矩阵具有相对独立性4、简述有限单元法中单元刚度矩阵的性质。 答:1.单元刚度矩阵是对阵矩阵 2.单元刚度矩阵的主对角线元素恒为正值 3.单元刚度矩阵是奇异矩阵 4.单元刚度矩阵仅与本身有关 5、简述有限元法中选取单元位移函数(多项式)的一般原则。 答:必须假定一个函数,所假定的位移函数必须满足两个条件:其一,它在单元节点上的值应等于节点位移;其二,由该函数出发得到的有限元解收敛于真实解。 6、要保证有限单元法计算结果的收敛性,位移函数必须满足那些条件?

答:1、完备性条件:要求单元的位移函数必须能够满足刚性位移和常量应变状态 2、协调性条件:要求单元的位移函数在单元部必须是连续函数,且必须保证相邻单元间位移协调 9、用有限元法分析实际工程问题有哪些基本步骤?需要注意什么问题? 1)建立实际工程问题的计算模型 2)选择适当的分析工具侧重考虑以下几个方面 1)前处理(Preprocessing) 2)求解(Solution) 3)后处理(Postprocessing 10、在弹性力学中根据什么分别推导出平衡微分方程、几何方程、物理方程,这三个方程分别表示什么关系? 答:根据静力学、几何学和物理学三方面条件,分别推导出平衡方程、几何方程和物理方程;三组方程分别表示:应力与载荷关系、应变与位移关系、应力与应变关系。 11、什么是平面应力问题?什么是平面应变问题?分别写出平面应力问题和平面应变问题的物理方程。 答:平面应力问题是指很薄的等厚度薄板,只在板边上受有平行于板面并且不沿厚度变化的面力,同时,体力也平行于板面并且不沿厚度变化。对应的应力分量只有σx,σy,τxy。而平面应变问题是指很长的柱形体,在柱面上受有平行于横截面并且

弹性力学与有限元分析考试试题及其答案

2012年某高校度弹性力学与有限元分析复习题及其答案 (内部资料) 一、填空题 1、弹性力学研究弹性体由于受外力作用、边界约束或温度改变等原因而发生的应力、形变和位移。 2、在弹性力学中规定,线应变以伸长时为正,缩短时为负,与正应力的正负号规定相适应。 3、在弹性力学中规定,切应变以直角变小时为正,变大时为负,与切应力的正负号规定相适应。 4、物体受外力以后,其内部将发生内力,它的集度称为应力。与物体的形变和材料强度直接有关的,是应力在其作用截面的法线方向和切线方向的分量,也就是正应力和切应力。应力及其分量的量纲是L -1MT -2。 5、弹性力学的基本假定为连续性、完全弹性、均匀性、各向同性。 6、平面问题分为平面应力问题和平面应变问题。 7、已知一点处的应力分量100=x σMPa ,50=y σMPa ,5010=xy τ MPa ,则主应力=1σ150MPa ,=2σ0MPa ,=1α6135' 。 8、已知一点处的应力分量, 200=x σMPa ,0=y σMPa ,400-=xy τ MPa , 则主应力=1σ512 MPa ,=2σ-312 MPa ,=1α-37°57′。 9、已知一点处的应力分量,2000-=x σMPa ,1000=y σMPa ,400-=xy τ MPa ,则主应力=1σ1052 MPa ,=2σ-2052 MPa ,=1α-82°32′。 10、在弹性力学里分析问题,要考虑静力学、几何学和物理学三方面条件,分别建立三套方程。 11、表示应力分量与体力分量之间关系的方程为平衡微分方程。 12、边界条件表示边界上位移与约束,或应力与面力之间的关系式。分为位移边界条件、应力边界条件和混合边界条件。 13、按应力求解平面问题时常采用逆解法和半逆解法。 14、有限单元法首先将连续体变换成为离散化结构,然后再用结构力学位移法进行求解。其具体步骤分为单元分析和整体分析两部分。 15、每个单元的位移一般总是包含着两部分:一部分是由本单元的形变引起的,另一部分是由于其他单元发生了形变而连带引起的。 16、每个单元的应变一般总是包含着两部分:一部分是与该单元中各点的位置坐标有关的,是各点不相同的,即所谓变量应变;另一部分是与位置坐标无关的,是各点相同的,即所谓常量应变。 17、为了能从有限单元法得出正确的解答,位移模式必须能反映单元的刚体位移和常量应变,还应当尽可能反映相邻单元的位移连续性。 18、为了使得单元内部的位移保持连续,必须把位移模式取为坐标的单值连续函数,为

相关主题
文本预览
相关文档 最新文档