2-杆梁结构的有限元分析原理1

合集下载

杆梁结构的有限元分析原理

杆梁结构的有限元分析原理

e
下面考察该简单问题的FEA求解过程。 (1) 离散化
两个杆单元,即:单元①和单元②
(2) 单元的特征及表达
对于二结点杆单元,设该单元的位移场为 么它的两个结点条件为
,那
设该单元的位移场具有模式(考虑两个待定系数)
利用结点条件,可以确定系数a0和a1,即
将系数a0和a1代入
,可将
表达成结点位移(u1, u2)的关系,即
其中, 为整体坐标系下的单元刚度矩阵, 为 整体坐标系下的结点力,即
由最小势能原理(针对该单元),将 对待定的 结点位移向量 取一阶极小值,有整体坐标系中 的刚度方程
对于本节给出的杆单元,具体有
4.3.3 空间问题中杆单元的坐标变换
就空间问题中杆单元,局部坐标系下的结点位移还 是 而整体坐标系中的结点位移为
这时由全部结点位移[0 u2 u3]分段所插值 出的位移场为全场许可位移场。
由最小势能原理(即针对未知位移u2和u3求 一阶导数),有
可解出
(5) 计算每个单元的应变及应力
在求得了所有的结点位移后,由几何方程
可求得各单元的应变
由方程 可求得各单元的应力
(6) 求结点1的支反力
就单元 ①的势能,对相应的结点位移求极值,可以 建立该单元的平衡方程,即
其中
由一维问题几何方程和物理方程,则该单元 的应变和应力为
其中
单元的势能
其中 叫做单元刚度矩阵。
叫做单元结点外载。
在得到“特征单元”的单元刚度矩阵和单元 结点外载后,就可以计算该单元的势能,因 此,计算各单元的矩阵 和 是一个关 键,下面就本题给出了个单元的 和 。
具体就单元①,有 单元①的结点位移向量
(5) 单元的刚度方程

2_杆系结构有限元分析1

2_杆系结构有限元分析1

( x) Nii N j j
x x N 1 , N 其中 i 为形函数。 j l l
由材料力学扭转可知
d dN e e M GI p GI p θ GI p B θ dx dx
其中 B
dN 1 1 dx l l
§1-2 扭转杆单元
e
外力势能 V u
e

e T
fe
e
1 e T e e e T 总势能 U V u K u u f e 2
e e
§1-1 拉(压)杆单元
1 e T e e e T U V u K u u f e 2
e e e
根据最小势能原理,势能泛函取驻值的必要条件
空间杆单元坐标变换矩阵
0 T 0
单元在两个坐标系中刚度矩阵转换关系同样有
K e T T K ' T
e
矩阵中仅仅包含有坐标的倾角,仅平行移动坐标轴,刚度矩阵 中元素值不变,矩阵的阶数也不改变。
§1-2 扭转杆单元
结点位移向量θe i , j
T
结点力向量
平衡关系
杆单元结点力向量
f U i
e
Uj
T
单元在外力和内力作用下处于平衡状态,反映单元平衡状态 的关系式就是刚度方程。下面利用最小势能原理推导单元的 刚度方程。 最小势能原理:在满足连续条件和边界条件的位移中,满足 平衡条件的位移其总势能最小,反之亦然。 单元总势能
e U e V e
M e Mi , M j
T
杆件发生自由扭转时,待求位移是截面的扭转角 ( x) 在局部坐标系中,每一个点将具有一个基本未知位移,最简单 的单元位移函数可以设为

杆梁结构的有限元分析原理[详细]

杆梁结构的有限元分析原理[详细]

le
EAe
le
EAe
u1 u2
P1
le
P2
u1 u2
1 qeTK eqe PeTqe 2
刚度矩阵
节点力列阵
3)离散单元的装配
在得到各个单元的势能表达式后,需要进行离散单元的装配,以
求出整个系统的总势能,对于该系统,总势能包括两个单元部分
e 1 2
1 q1T K1q1 q2T K 2q2 P1Tq1 P2Tq2 2
第4章 杆系结构的有限元分析原理
杆梁单元概述
讨论杆梁单元和由它们组成的平面和空间杆梁结构系统. 从构造上来说其长度远大于其截面尺寸的一维构件 承受轴力或扭矩的杆件成为杆 杆梁问题都有精确解 承受横向力和弯矩的杆件称为梁 平面桁架 平面刚架 连续梁 空间刚架 空间桁架等 承受轴力或扭矩的杆件称为杆 将承受横向力和弯矩的杆件称为梁 变截面杆和弯曲杆件
单元节点条件:u(0)=u1, u(l)=u2
从而得
a0 ui ,
a1
uj
le
ui
i
1,
j
2
回代得
u(x) a0 a1x
ui
u j ui le
x
1
x le
ui
x le
u
j
Niui N ju j
其中Ni,Nj是形函数。
写成矩阵形式为
q Niu Nqe
N
ju
ui u j
1 2
u1
EA1
u2
l1 EA1
l1
EA1
l1
EA1
u1 u2
R1
l1
0
u1 u2
1 2
u2
EA2

梁的有限元分析原理

梁的有限元分析原理

梁的有限元分析原理梁的有限元分析原理是一种工程结构分析方法,广泛应用于建筑、桥梁、航空航天、汽车等领域。

它通过将连续的结构离散化为有限数量的小单元,通过数学模型进行计算,得出结构的力学性能和响应情况。

梁的有限元分析原理是有限元分析的基础,下面将对其进行详细介绍。

首先,梁的有限元分析原理基于梁理论,即在横向较小、纵向较长的情况下,结构可以近似为一维梁。

梁的有限元分析原理通过将梁划分为多个单元,每个单元内部可以看作两个节点之间的一段杆件,通过建立节点之间的力学关系方程,得到整个结构的力学性能。

其次,梁的有限元分析原理利用了变分原理,即将结构的势能取极小值,建立了结构的力学方程。

通过对于梁的弯曲、剪切和轴向力等方面的力学模型进行合理的假设与简化,可以得到结构的位移与力的关系,从而解决结构的力学问题。

在梁的有限元分析中,需要进行以下几个步骤:1.几何离散化:将梁结构划分为多个单元,每个单元具有相同的形状与尺寸,通常为矩形或三角形。

2.模型建立:根据梁理论以及力学方程的简化假设,建立节点的力学关系方程,包括位移、应力、应变等参数。

3.材料性能定义:确定梁材料的力学性能参数,如弹性模量、截面惯性矩等。

这些参数对梁结构的力学性能具有重要影响。

4.边界条件施加:根据实际问题设定边界条件,包括固定支座、约束条件等。

这些条件对于解决梁结构的位移、应力等问题至关重要。

5.方程求解:通过数学方法求解得到节点之间的力学关系方程,利用数值计算技术进行迭代求解,得到梁结构的位移、应力等参数。

6.结果分析:根据求解得到的结果,进行力学性能分析,如最大应力、挠度、模态分析等。

根据分析结果评估结构的强度与稳定性。

总结起来,梁的有限元分析原理是一种基于梁理论的工程结构分析方法,通过将结构离散化为多个小单元,利用力学关系方程和数值计算技术求解得到结构的力学性能。

通过梁的有限元分析原理,工程师可以更加准确地评估结构的强度与稳定性,对结构进行优化设计。

杆梁结构有限元分析

杆梁结构有限元分析

3.1 杆梁结构的直接解法
机械分社
(1)平面压杆有限元法的直接法
由节点平衡有: 即有:
U1(1)u1 U1(1)u2 N1
U
u (1)
21
(U
(2 2
)
U
(1) 2
)u2
U
(2 2
)u3
F1
U
(2 3
)
u2
U
(2 3
)
u3
F2
EA1 l1
u1
EA1 l1
u2
N1
EA1 l1
u1
( EA1 l1
3.1 杆梁结构的直接解法
机械分社
杆梁结构是指长度远大于其横截面尺寸的构件组成的杆 件系统,例如机床中的传动轴,厂房刚架与桥梁结构中的梁 杆等,可以用杆单元或梁单元来进行离散化。
空间杆系:平面杆系是指各杆轴线和外力作用线位于一 个平面内,若各杆轴线和外力作用线不在一个平面内。 (1)平面压杆有限元法的直接法
单元刚度矩阵每一列元素表示一组平衡力系,对于平面 问题,每列元素之和为零。
3.1 杆梁结构的直接解法
机械分社
(2)平面梁单元有限元法的直接法 2)节点位移与节点力之间的关系
Ui
Vi
k11
k21
M i U j
k31
k41
V
j
M j
k51
k61
他们在轴和轴的投影之和等于零:
vi
6EI l2
i
12EI l3
vj
6EI l2
j
M
j
6EI l2
vi
2EI l
i
6EI l2
vj
4EI l

桥梁结构分析的有限元原理及其程序简介

桥梁结构分析的有限元原理及其程序简介


e FEe = K E Rδ e
其中 R 为坐标变换矩阵。 若 e 号单元内还作用有跨间荷载以及给定的温 度分布,它们在局部坐标系下的单元等效结点荷载 分别记为 Pqe 和 PTe ,则
e e FE = ΚE Rδe − Pqe − Pte
以上即杆系结构有限元法的基本计算过程。
1.2 有限元软件简介
1.2有限元软件简介
与通用有限元的区别
ANSYS MIDAS/CIVIL
前处理 单元、材料、边界、荷载
前处理 单元、材料、边界、荷载、施工过程、 预应力、收缩徐变等 求解 静力、动力、稳定等 后处理 显示、列表、时程等 设计验算 基于规范的荷载组合、 设计验算
求解 静力、动力、稳定等
后处理 显示、列表、时程等
1. 桥梁结构分析的内容
• (1)桥梁一般是分阶段逐步施工完成的,结构最终受力 状态往往与施工过程有着很大的关系,因而结构分析必须 按实际的施工过程和结构形成的过程逐阶段进行分析,并 且能够自动累加各阶段的内力和位移等。 (2)计算成桥后在二期恒载,支座不均匀沉降、混凝土 长期收缩、徐变效应、温度变化等作用下的内力和位移。 (3)计算各种活载引起的内力和位移,包括影响线或影 响面的计算以及对它们进行纵向、横向的加载等。 (4)计算各种偶然荷载(加地震)等引起的内力和位移。 (5)按规范对上述各种荷载引起的内力和位移进行组 合,得出最不利的组合情况。 (6)按规范进行强度、刚度、抗裂性、稳定性以及动力 性能验算。
2.2 桥梁结构分析的施工过程及体系转换 • 比如,同为三跨连续梁,在合拢的先后顺 序上,先合拢边跨还是中跨对结构成桥内 力是有影响的; • 有时为了获得良好的成桥线形或内力,可 以在施工中采取一些辅助措施。

有限元-梁系结构的有限元法

有限元-梁系结构的有限元法

4x l
3x 2 l2
) i
x l
(3x l
2)
j
容易验证 : x 0: u ui v vi i x l: u u j v v j j
(3-1a),(3-1b)或(3-2a),(3-2b)称为平面梁单元的位移插值 函数
二、建立节点位移与节点力关系
1、 轴向节点力
E Fx A
拉压杆问题的回顾
1、杆的基本概念:
杆--轴线为直线的细长构件,沿轴线承受 拉(压)载荷; 杆模型--平面假设将杆简化为一维问题, 可由杆轴线代表; 杆变形特点--只与轴向位移相关;
拉压杆问题的回顾
2、杆有限元的基本概念
节点位移—轴向位移,每节点1个自由度; 节点力—轴力; 结构离散:轴线划分为若干直线段; 单元分析:建立节点力与节点位移关系; 节点平衡:对每一节点,建立相关节点力与 外力的平衡关系,得到一线性方程组; 约束处理:引入已知节点位移,使方程组可解
梁系结构实例
2、平面梁系
1、节点力平衡的需求--单元节点力(在 局部坐标系中)向整体坐标系的变换; 2、单元分析的需求--节点位移(在整体 坐标系中)向局部坐标系的变换; 3、结构对称性的利用(练习,作业3)。
l2 2EI
l
0
Vi
i
u
j
(3-4)
6EI l2
4EI
V
j j
l
(3-4)式是用矩阵表示的梁节点力与节点位移的关系
式(3-4)还可写成:
F
e
K e
e
(3-5)
e
F
——称为局部坐标下的节点力列向量
e ——称为局部坐标下的节点位移列向量
e
K

梁的有限元分析原理

梁的有限元分析原理
y
j
·
x

Chapter 5 Bernoulli-Euler Beam
z
27
福州大学研究生课程-有限元程序设计
平面桁架杆单元(2D LINK1)
空间杆单元(3D
LINK8)
平面刚架,BEAM3 空间梁单元(BEAM4)
Chapter 5 Bernoulli-Euler Beam
28
福州大学研究生课程-有限元程序设计
举例说明
Chapter 5 Bernoulli-Euler Beam
18
福州大学研究生课程-有限元程序设计
这种高斯积分阶数低于被积函数所有项次精确 积分所需要阶数的积分方案称之为减缩积分。 实际计算表明:采用缩减积分往往可以取得较 完全积分更好的精度。这是由于: 精确积分常常是由插值函数中非完全项的 最高方次要求,而决定有限元精度的是完全多 项式的方次。这些非完全的最高方次项往往不 能提高精度,反而可能带来不好的影响。取较 低阶的高斯积分,使积分精度正好保证完全多 项式方次的要求,而不包括更高次的非完全多 项式的要求,其实质是相当用一种新的插值函 数替代原来的插值函数,从而一定情况下改善 19 Chapter 5 Bernoulli-Euler Beam 了单元的精度。
福州大学研究生课程-有限元程序设计
有限元程序设计
——梁单元,静力问题
谷 音 福州大学土木工程学院
2012
1
福州大学研究生课程-有限元程序设计
§1. 介绍. 框架结构,例如桁架、桥梁 轴力构件 axial elements 杆 受弯构件 flexural elements 梁 平面梁单元 plane beam element
Chapter 5 Bernoulli-Euler Beam

杆结构 分析的有限元方法(有限元)

杆结构   分析的有限元方法(有限元)
局部坐标系中的单元述
杆单元形状函数
杆单元刚度矩阵
平面问题中的坐标变换
梁结构分析的有限元方法
梁:承受横向荷载和弯矩的杆件。
梁的主要变形为挠度v
横截面变形前后都垂直于杆变形前的轴线x轴
中性层变形=0
纯弯曲没有剪力,只有弯矩
梁截面的惯性矩
杆结构分析的有限元方法
杆:承受轴向荷载的杆件
最基本的承力结构件:杆、梁
弹簧--简单的承受轴力的结构件
有限元方法中,每一个处理步骤都是标准化和规范化的,
因而可以在计算机上通过编程来自动实现。
F=kδ
k--刚性系数
位移的绝对变化量/杆件的伸长量δ=u2—u1
应力某截面上单位面积上的内力/内力的分布集度
应变相对伸长量单位长度的伸长量
杆单元的特性是节点位移及节点力的方向都是沿轴线方向。
杆结构的力学分析
铰接的杆结构----杆只受轴力-----杆件拉伸问题---可自然离散
两端为铰接的杆件只承受轴力。
各个单元研究(基于局部坐标系的表达)
各个单元研究
离散单元的集合、组装
杆单元及坐标变换
自由度:描述物体位置状态的每个独立变量。
对于杆单元,其节点位移有两个自由度。

杆梁结构的有限元分析原理

杆梁结构的有限元分析原理

杆梁结构的有限元分析原理杆梁结构是工程中常用的一种结构形式,它由多个杆件或梁组成,用于承担载荷和传递力量。

有限元分析是一种通过将结构离散为许多小单元,利用数学方法对结构进行分析的技术。

下面将详细介绍杆梁结构的有限元分析原理。

一、杆件离散化在有限元分析中,首先需要将杆梁结构离散化为一组子结构,即离散化为一组离散的杆件。

离散后的每个杆件可以看作是一个子系统,每个子系统由两个节点组成,节点之间以杆件连接。

通过节点与杆件的连接方式,能够模拟出整个杆梁结构的受力特点。

离散化的过程中,需要确定杆件的几何形状、截面以及材料特性等参数,并根据实际情况设置合适的杆件单元数目。

通常,单元数目越多,离散程度越高,结果越接近真实情况,但计算成本也会增加。

二、有限元法的基本原理有限元方法的基本原理是将结构分成许多小的单元,每个单元内的行为可以用简单的数学函数来表示。

对于杆梁结构,常用的单元有梁单元和杆单元。

梁单元适用于承受弯曲强度较大的杆件,而杆单元适用于承受轴向载荷的杆件。

通过将结构分成小单元后,可以建立一个与原结构相似的离散模型,并在每个单元上建立相应的方程。

三、应力应变关系在进行有限元分析时,需要获得每个杆件的应变和应力。

应变与杆件的变形有关,而应力与应变之间的关系则与材料的本构关系有关。

对于线弹性材料,应力与应变之间可以通过胡克定律来描述。

胡克定律表明,应力与应变之间成线性关系,材料的弹性模量E、泊松比ν以及应变关系能够决定应力。

应根据结构中不同材料的应变特性来选择相应的材料模型。

四、施加边界条件在进行有限元分析前,需要施加适当的边界条件。

边界条件用于模拟实际情况中的约束和限制。

常见的边界条件有固定边界、弹性边界和施工阶段边界。

五、求解位移和应力当离散化杆梁结构、建立了位移和应变关系、施加了边界条件之后,可以通过数值求解方法,例如有限元法中的坐标变形法,计算得到结构的位移和应力。

坐标变形法能够基于得到的位移结果,进一步计算应力。

第二章-杆和梁结构的有限元法案例

第二章-杆和梁结构的有限元法案例

第二章
杆和梁结构的有限元法
§2.1.2 弹簧系统分析
注意: 上述弹簧系统的分析求解原理和过程就是有限元 法求解连续体力学问题时对离散后系统的分析求 解原理和过程。
第二章
杆和梁结构的有限元法
§2.1.2 弹簧系统分析
例题1:弹簧系统
已知条件:
求:(a) 系统总刚度矩阵 (b) 节点2,3的位移
单元特性
系统平衡方程
第二章 杆和梁结构的有限元法
KD F
2)单元方程扩大相加法 单元特性
F1 f11
相加
F2 f 21 f12 F3 f 22
系统节点 平衡条件
引入系统节点平衡条件
KD F
系统节点平衡方程
第二章 杆和梁结构的有限元法
2.2 杆单元和平面桁架
杆单元
2.2.1 一维等截面 杆单元
fi k f j k
第二章
k ui k u j
f kd
杆和梁结构的有限元法
2、弹簧系统的集成 1)列节点平衡方程法
F1 f11 F2 f 21 f12 F3 f 22
系统节点 平衡条件
F1 k1u1 k1u2 F2 k1u1 ( k1 k2 )u2 k2u3 F3 k2u2 k2u3
第二章 杆和梁结构的有限元法
k k k
k k
fi k f j k
k ui k u j
kii k k ji
kij k jj
§2.1.2 弹簧系统分析
求解一个弹簧系统:
1)各单元的特性分别为:
第二章 杆和梁结构的有限元法

《有限元理论与数值方法》第三讲-杆、梁结构有限元分析

《有限元理论与数值方法》第三讲-杆、梁结构有限元分析
杆件结构可分为桁杆和梁两类。 由杆件组成的结构体系称为杆系。由桁杆组成的杆系称为桁架; 由梁组成的杆系称为刚架。若杆系和作用力均位于同一平面内,则称 为平面桁架或平面刚架,否则称为空间桁架或空间刚架。
Finite Element Theory and Numerical Method
一、杆、梁的物理力学模型
拉压杆单元如图3-6所示,已知等直杆件杆长为 l 横截面面积为 A 材料弹性模量为 E 所受轴向分布载荷集度为 p(x) 杆端位移分别为 u1 u2
杆端力分别记为 F1 F2
1、建立位移场
F1, u1 xa
1
a p(x)
2 F2 , u2
x
设局部坐标系下杆中任意点a的坐标为 xa
因为只有两个边界条件 u1
形函数具有如下性质: 1)本端为1,它端为0 2)单元内任意一点总和为1
N1(0) 1
N1(1) 0
N2 (0) 0 N2 (1) 1
N1() N2 () 1
2、应变分析
du dx
dN dx
ue
dN1 dx
B为应变矩阵或者几何矩阵。
dN2 dx
u
e
1 l
1 l
ue
[B1
B2 ]ue Bue
图示所示桁架 l 2m
EA 1.2106 kN
试求1-2杆和1-4杆单元的局部坐标单元 刚度矩阵
1-2杆:抗拉刚度 EA / l 6106 kN/m
F1 10N 3
1
F2 20N 4
2
ke1
EA l
1 1
1
1
6
105
1 1
1
1
kN
/
m
1-4杆:抗拉刚度 EA /( 2l) 4.24264 105 kN/m

有限单元法课件第四章 杆件系统的有限元法

有限单元法课件第四章 杆件系统的有限元法
桁杆 梁
(a)
(b)
由杆件组成的结构体系称为杆系,如起重机,桥梁等。
由桁杆组成的杆系称为桁架。
由梁组成的杆系成为刚架。
若杆系和作用力均位于同一平面内,则称为平面桁架 或平面刚架,否则称为空间桁架或空间刚架。
由于杆件结构采用一维单元进行离散,所以杆系的网 格划分容易用半自动方法实现。当采用自动网格划 分方法时,杆系的几何模型是由杆件轴线构成的线框 模型。
R
e P
RiP R jP
R
lP
R
R
e F
RiF R jF
Rlx Rly NlT l R l
lF T l
Px dx (l i, j ) Py
e T
Bj dx
kii k ji
kij k jj
其中矩阵元素为
kst D Bt dx B as 0 EA 0 at 0 0 0 bs dx 0 EI 0 bt ct 0 cs 0 0 EAas at dx 0 EIb b EIb c s t s t EIcs bt EIcs ct 0
e
du dx e x 2 B Bi q x d v dx 2
Bj q
e
其中
ai 0 0 Bi 0 b c i i a j 0 0 Bj 0 b c j j 1 12 6 ai a j bi b j 3 x 2 l l l 4 6 2 6 ci 2 x cj 2 x l l l l

桥梁结构分析的杆系有限元法及结构模型的建立2015

桥梁结构分析的杆系有限元法及结构模型的建立2015

结构的离散化
确定了结构的全部 节点,也就确定了 结构的单元划分, 然后对结构进行单 元编号和节点编号, 通常单元编号用①, ②,……表示,节 点编号用1, 2,……表示,如图 所示。
6 67
5
4
3
5
4
1
2
1
2
3
单元杆端力与杆端位移的表示方法
• 平面桁架单元的局部坐标和整体坐标:
y
y
x
3
x2
2
y
1
结构分析的杆系有限元法
• 概述 • 有限单元法的概念及应用 • 结构的离散化 • 单元杆端力与杆端位移 • 逆步变换 • 单元刚度矩阵 • 总刚度矩阵 • 边界条件的后处理法 • 线性代数方程组的数值解法
结构分析的含义
• 结构分析的含义,不仅指在一定的已知条件下对结构的变 形和内力等进行计算,而且包括分析构件刚度变化对内力 变化的影响,对结构的几何组成进行分析,以及选择合理 的结构形式等等。
结构分析的有限元法
• 美国20世纪70年代推出的至今仍然是世界销售量最大的 NASTRAN(NAsa STRuctural Analysis,美国国家航空和 宇宙航行局结构分析程序系统)程序与当时西德推出的 ASKA(Automatic System for Kinematics Analysis,运动 分析的自动程序系统)齐名,同为当时最为著名和广泛应 用的程序,但几十年后的现在,ASKA已无法与 NASTRAN相比。原因是ASKA后来没有大规模的资金投 入,使程序不断得到滚动发展(维护)和组织推广、剌激 程序在竞争中不断改进各种功能。
向量
X
e i
Yi e
F
e
Fi e Fje

有限单元法 第2章 杆系结构的有限元法分析

有限单元法 第2章 杆系结构的有限元法分析

义 & 可以进一步求得单元刚度矩阵为 )
( & # 0# ( $’ $ % 8 . ! 1 # $ ’ 0# # 同时 & 我们可以根据式 $ % 求出等 效 结 点 荷 载 矩 阵 ’ 这 里 要 指 出 的 是 ) 分 布 荷 载 ! .$
! # !! !
! # $! !
! 第 ! 章 ! 杆系结构的有限元法分析 # #! ! """""""""""""""""""""""""""""""""""""""""""
不适定的 " 第九步 # 求解方程组 " 计算结构的整体结点位移列阵 ## 并 进一步 计算各 单元 的应力 分量及主应力 $ 主向 " 第十步 # 求单元内力 # 对计算成果进行整理 $ 分析 # 用表格 $ 图线标示出所需的位移 及应力 " 大型商业软件 % 如 )* + , + 等 & 一般都具有强大的后处理功能 # 能够 由计算 机自 动绘制彩色云图 # 制作图线 $ 表格乃至动画显示 "
矩阵 ’ $ %进行应力 ( 应变分析 ’ 根据材料力学中应变的定义 & 有 ) ! # # $’ 2 + 2 $ ( ( ( ( $’ $’ $’ . 0 ! ! . " 3 3 .% ". . ! ! ! !! "# ’ ’ 2 # 2 #

第三讲 杆件结构有限元分析

第三讲 杆件结构有限元分析
du dx

l
0
AE
l du d u dx f x udx 0 dx dx
其中E表示弹性模量,A表示横截面积,方程左端得到单元的刚度矩阵。
建立有限元模型
现考虑一个由5个长度相同(le=1m)横截面积不同的杆件构成的一维杆件,各杆弹性模量都为 E=1.0e10pa,A1=0.5m2,A2=0.4 m2,A3=0.3 m2,A4=0.2 m2,A5=0.1 m2,如图1所示,右端给定位移 u右=0.1,左端固定位移u左,分析杆件内位移分布:
根据虚功原理,方程两边乘以虚位移δu,平衡方程可以写为:

其弱形式为:
l
0
[
d ( A x ) f ( x)] udx 0 dx

l
0
A x
l d u dx f x udx Pj u j 0 dx j
基本方程的最终弱形式
其中,右端最后一项可以看作是节点力情况,所以可以不单独列出,同时 x E 所以上式可以继续写为:
网格尺寸设置
网格划分信息
网格划分
选择calculate → calculate,在弹出的对话框,点击OK,保存,前处理完毕。
工程求解
点击工具栏中“求解计算”按钮,完成模型的求解计算。
后处理
点击工具栏中的“后处理”按钮进入GID,查看计算结果,如下图所示。
结果分析: 本章针对一个变截面一维杆件,通过理论分析和ELAB1.0软件实现两种方式来分析,一方面对有限元
几何模型
将其划分为五个单元六个节点,即每根杆件作为一个单元,每个单元的节点关系如下图所示:
单元拓扑关系
确定杆单元的形函数
考虑其中一个杆单元,其两个端点分别为节点1,节点2,基本变量为节点位移u1,u2::

有限元方法第三章杆系结构有限元

有限元方法第三章杆系结构有限元
稳定性以及波浪载荷的影响。
应用实例
某大型桥梁的稳定性分析
采用杆系结构有限元对某大型桥梁进行稳定性分析,评估其在不同载 荷下的变形和承载能力。
高层建筑的抗震性能研究
利用杆系结构有限元模拟高层建筑的抗震性能,分析地震作用下结构 的响应和破坏模式。
汽车悬挂系统的优化设计
通过杆系结构有限元模拟汽车悬挂系统的运动和受力情况,优化悬挂 参数以提高车辆行驶的稳定性和舒适性。
有限元方法第三章杆系结 构有限元
• 引言 • 杆系结构有限元的基本概念 • 杆系结构有限元的建模方法 • 杆系结构有限元的求解方法 • 杆系结构有限元的应用案例 • 结论与展望
01
引言
目的和背景
杆系结构是工程中常见的一种结构形式,广泛应用于桥梁、 建筑、机械等领域。由于其具有复杂的几何形状和受力特性 ,因此需要采用有限元方法进行数值分析。
THANKS
感谢观看
04
杆系结构有限元的求解方法
求解步骤
确定边界条件
根据实际情况,确定杆系结构 的边界条件,如固定、自由、 受压等。
求解线性方程组
将所有单元的平衡方程组合成 一个线性方程组,然后使用数 值方法求解该线性方程组。
建立离散模型
首先将杆系结构离散化为若干 个小的单元,每个单元具有一 定的物理属性。
应用力学平衡方程
杆系结构有限元的优缺点
优点
能够处理复杂的几何形状和边界条件, 适用于大规模问题求解,计算精度可 调,可模拟复杂的结构和场。
缺点
需要针对不同的问题建立不同的模型, 计算量大,需要较高的计算机资源, 对于非线性问题求解较为困难。
03
杆系结构有限元的建模方法
建模步骤
确定研究问题

杆梁结构的有限元分析原理

杆梁结构的有限元分析原理

杆梁结构的有限元分析原理杆梁结构是一种常见的工程结构,广泛用于建筑、桥梁、机械等领域。

为了研究杆梁结构的力学性能和设计优化,常用的方法之一是有限元分析。

有限元分析是一种数值计算方法,通过将连续结构离散化为一个个有限的单元(元素),再通过计算单元之间的相互作用来近似表示整个结构的力学性能。

下面将逐步介绍杆梁结构的有限元分析原理。

1.离散化:首先,将杆梁结构离散化为一个个的单元,通常可以选择线性单元、二次单元等。

线性单元简单且计算效率高,而二次单元更准确但计算开销较大。

根据具体工程需求和分析要求,选择合适的单元进行离散化。

每个单元由节点和单元梁组成。

2.建立本地坐标系:为了方便计算,对于每个单元,可建立本地坐标系。

本地坐标系是以单元的一个节点为原点,并建立与该节点有关的坐标轴。

通过本地坐标系可以方便地描述单元内部的各种力和力矩。

3.单元刚度矩阵计算:对于每个单元,需要计算其刚度矩阵。

刚度矩阵描述了单元内部的相互作用,包括节点间的弯曲刚度和剪切刚度等。

通过根据材料的力学特性和几何信息,可以得到单元刚度矩阵。

4.装配全局刚度矩阵:将所有单元的刚度矩阵按照它们的几何关系组装成全局刚度矩阵。

全局刚度矩阵描述了整个杆梁结构的力学行为。

5.施加边界条件和加载情况:根据具体问题的边界条件和加载情况,在全局刚度矩阵中添加与之对应的约束和加载项。

边界条件通常涉及到约束的位移和力的平衡,加载情况则涉及到外界施加在结构上的力。

6.求解杆梁结构的位移:通过求解全局刚度矩阵与位移的乘积等式,可以得到结构的位移。

位移是描述结构变形的重要参数,可以用来计算应力、应变和变形等。

7.计算应力和应变:通过已知的位移以及杆梁的几何信息,可以计算单元内部的应力和应变。

应力和应变是评估杆梁结构受力情况的重要指标,在结构设计和安全评估中具有重要作用。

8.结果后处理:最后,可以通过后处理技术对有限元分析的结果进行处理和展示。

例如,可以绘制位移云图、应力云图等,以方便工程师对结构的力学性能进行评估和优化。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
13
福州大学研究生课程-有限元程序设计
杆梁结构的有限元分析原理
§2.1 一维单元 -轴向受力杆
Step 1 离散 Discretization 将杆划分为几个单元,用标准化的分段小单元来逼近结 构
14
福州大学研究生课程-有限元程序设计
Step 2 选择单元的位移模式 寻找能够满足位移边界条件的许可位移场 对于如图杆单元,该单元位移场为u(x),那么它 的两个节点条件为:
27
福州大学研究生课程-有限元程序设计
求解支反力: 系统方程重组,分块:
=> 未知节点位移 => 求得支反力
28
福州大学研究生课程-有限元程序设计
Trivial expression 修改对角元素, α为一个非常大的数 ,1015
Advantage is 1) easy to be programmed 2) the properties of [K] are kept.
福州大学研究生课程-有限元程序设计
Step 2 选择单元的位移模式
根据节点条件,可以确定出待定系数:则
则u(x)表达称为节点位移(u1,u2)的关系,即
矩阵表达形式
N为形函数矩阵
q为节点位移列阵
17
福州大学研究生课程-有限元程序设计
Step 3 应变和应力Strain Field and Stress Field 由一维问题的几何方程和物理方程,可以得到位移与 应变的关系,则: where [B] is (应变矩阵)strain matrix which is given as
Chapter 1 About Finite Element Method
37
福州大学研究生课程-有限元程序设计
有限元分析的基本步骤
⑤引入边界条件 Impose the boundary conditions 修改系统方程 ⑥解方程 Solve the system equations 得到未知问题的节点值 ⑦后处理 Make additional computations if desired 例如节点位移 → 应变与应力
33
福州大学研究生课程-有限元Байду номын сангаас序设计
Step 8. 求解方程 Solving equations. Assume EA=const. L i = l0/ 3 = L , q(x)=q0 In the example
8
34
福州大学研究生课程-有限元程序设计
Step 9. 其他计算 Additional calculation
在工程实际中,杆单元可能处于整体坐标系中的任意 一个位置,因此需要将原来局部坐标系中所得到的单 元表达等价地变换到整体坐标系中,这样不同位置的 单元才有公共的坐标基准,以便对各个单元进行集成 的和装配。
41
福州大学研究生课程-有限元程序设计
局部坐标系中的节点 位移:
整体坐标系中 的节点位移:
42
福州大学研究生课程-有限元程序设计
建立立单元的平衡方方程
{F}e = ∫ [ B] E[ B] Adx{u}Δ[k ]{u}
0
Li
T
其中K为杆单 元的刚度矩阵: P为节点力力列阵:
二二结点一一维单元
单元势能的表达
=Ue-We
运用最小势能原理对位移求极值 建立单元的平衡方程
K e ⋅ qe = Pe
其中K为杆单 元的刚度矩阵: P为节点力列阵:
• 有限单元法的概念
• 通过材料力学求解和有限元求解进行比较 • 例:等截面直杆在自重作用下的拉伸图(a) • 单位杆长重量为q,杆长为L,截面面积为A, 弹性模数为E
4
福州大学研究生课程-有限元程序设计
5
福州大学研究生课程-有限元程序设计
• 等截面直杆在自重作用下的有限元法解答 • 1、离散化
单元的势能是一个标量,不会由于坐标系 的不同而改变,将变换后的参数带入,则
43
福州大学研究生课程-有限元程序设计
由最小势能原理,则得到整体坐标系中的刚度 方程为:
44
福州大学研究生课程-有限元程序设计
45
福州大学研究生课程-有限元程序设计
46
Chapter 1 About Finite Element Method
38
福州大学研究生课程-有限元程序设计
39
福州大学研究生课程-有限元程序设计
杆单元及其坐标变换
(1)
(2)
杆单元扩展用于桁架构件,对于平面桁架构件, 一个节点为两个自由度 如果是桁架构件,需要坐标变换
40
福州大学研究生课程-有限元程序设计
未知场函数(导数)在各结点上数值——自由度
Chapter 1 About Finite Element Method
11
福州大学研究生课程-有限元程序设计
Ø 杆系结构:梁、拱、框架、桁架等,它们常可离 散成杆元和梁元。
○ ○ ○ ○ ○ ○ ○ ○ ○
梁 拱 框架 桁架



○ ○
12
福州大学研究生课程-有限元程序设计
10
福州大学研究生课程-有限元程序设计
有限元如何实现?
分段近似 近似函数 联接简单函数 在小段区域内定义 (element) 利用每个单元内假设的近似函数来分片地表示全求解域 上待求的未知场函数
有限元法的基本思路
将连续系统分割成有限个分区或单元,对每个单元提出一 个近似解,再将所有单元按标准方法组合成一个与原有系 统近似的系统。
有限元程序设计方法简介
程序基本框图
1、输入基本数据(结构描述): (1)控制数据:如结点总数、单 元总数、约束条件总数等; (2)结点数据:如结点编号、结 点坐标、约束条件等; (3)单元数据:如单元编号、单 元结点序号、单元的材料特性、 几何特性等; (4)载荷数据:包括集中载荷、 分布载荷等。
开始 输入基本数据 计算单元刚度矩阵 形成总体刚度矩阵 形成结点荷载向量 引入约束条件 求解方程组,输出结点位移 计算单元应力,输出结果 结束
则设该单元的位移场具有线性模式 (考虑两个待定系数)
15
福州大学研究生课程-有限元程序设计
u =[ N
i
( x) i
N
( x) 2
⎧ u1 ⎫ ⎧ u1 ⎫ ] ⎨ ⎬ = [N ] ⎨ ⎬ ⎩u2 ⎭ ⎩u2 ⎭
i
i
i
= [1 − x / Li
⎧ u1 ⎫ x / L] ⎨ ⎬ ⎩u2 ⎭
在大多数情况中,外力并不直接作用于节点,例如: 沿单元轴向施加的剪切力。因此,要将外力转换为 节点等效荷载,
F : concentrated axial loads;
24
福州大学研究生课程-有限元程序设计
Step 6 单元集成或组装Assemblage of Elements
该系统的组装程序是基于单元节点的相容性条件和平 衡条件。相容性:通过节点连接的单元的未知节点变 量是相同的。平衡条件:某一方向上的节点力之和要 等于作用在该节点上的外力合力。 离散单元进行装配,整体刚度矩阵由两个单元的 刚度叠加而成 K=K(1)+K(2)+… 与整体节点位移向量相互对应 {δ }T=[u1 u2 u3 u4]
• 2 把外载荷集中到节点上 • 把第i 单元和第i+1 单元重量的一半 • 集中到第i+1 结点上
6
福州大学研究生课程-有限元程序设计
• 3、假设线单元上的位移为线性函数 • 用单元节点位移表示单元内部位移 • 第 i 个单元中的位移用所包含的结点位移来表示,
7
福州大学研究生课程-有限元程序设计
福州大学研究生课程-有限元程序设计
有限元程序设计
谷 音 福州大学土木工程学院
2015
1
福州大学研究生课程-有限元程序设计
有限元法
n n n n
直接法(一维杆单元) 虚功原理 变分法 加权余量(残值)法
2
福州大学研究生课程-有限元程序设计
杆梁结构的有限元分析原理
福州大学研究生课程-有限元程序设计
应变矩阵(几何函数矩阵)[B]是一个与x无关的常矩 阵。 E 为杨氏弹性模量
18
福州大学研究生课程-有限元程序设计
Step 4 单元的平衡 Element Equilibrium Relation 节点力Fi作用在单元节点 虚位移
列向量表示单元节点位移
19
福州大学研究生课程-有限元程序设计
由虚功原理
8
福州大学研究生课程-有限元程序设计
9
福州大学研究生课程-有限元程序设计
有限元法的一般概念
n
n
加权余量方法和变分法通过对未知场函数进行 试函数近似,能把连续问题化为离散问题,但 求解能力有限,一个主要障碍就是试函数是在 全场范围内定义的。 有限元法的思想是把整个求解区域分成一些既 不重叠也没有裂缝的子区域(单元)-瑞利里 兹法。以未知函数在节点上的值作为基本未知 参数,运用各种方法得到单元的刚度方程,这 样就可以用离散的方法求解。
35
福州大学研究生课程-有限元程序设计
结果
补救方法: 增加自由度 increase the dof
36
福州大学研究生课程-有限元程序设计
有限元分析的基本步骤
① 离散连续体Discretize the continuum (前处理 preprocessing)划分单元 ②选择插值函数 Select interpolation functions 定义节点 Assign nodes 选择插值函数Choose the interpolation function ③单元分析Find element properties 平衡方程,应力应变描述Equilibrium, expression of strains, stresses ④ 组合集成单元Assemble the element properties 获得整体系统方程to obtain the system equations The expression of the behaviour of the system
相关文档
最新文档