平面桁架有限元分析及程序设计
桁架的有限元分析w
桁架的有限元分析问题:已知模型由轨座,桁架以及垫块三个部分组成。
载荷为186吨,加载与距中间400mm 处的轨座上。
材料为Q345钢,密度为33/108.7m kg ⨯,弹性模量为5102⨯MPa ,泊松比为0.3,摩擦系数为0.05。
要求对模型使用ANSYS 进行有限元分析,分析其安全性。
问题分析:模型是对称的,可以简化模型,用模型的四分之一进行ANSYS 有限元分析,可以减少计算量。
需要在两个剖面处施加对称约束,从而保证中间部分不发生位移。
垫块的底部需要添加一个全约束以防止刚体位移。
图1为简化的模型。
图1 模型建模以及运算1将模型导入将所给的模型文件dggl.sat 导入到ANSYS 中,补充完整模型,即添加垫块。
再将模型分割,得到简化后的模型,即图1。
2 划分网格将轨座和桁架分割成规整的方体,使用映射网格来将轨座和桁架划分成六面体网格,可以得到比较规整的网格。
使用扫掠划分将垫块划分成六面体网格,网格大小设定为15mm。
划分结果如图2。
图2 网格化分3 设置单元类型对于实体模型分析,我们可选用8节点SOLID185单元。
整个分析过程有关于非线性接触的问题,所以要设置接触对单元类型。
选择TARGET170和CONTACT174单元。
4接触对创建使用设置接触对向导Contact Manager来设置。
设置轨座下表面和垫块上表面作为接触面,桁架为两接触面所对应的目标面。
其中轨座与桁架的接触对需设置成绑定接触,以防止发生滑移。
创建的接触对如图3所示。
图3 接触对的创建图4 边界条件5 添加约束施加约束,要在整体模型的中间部分施加对称约束以及对垫块施加全约束,从而保证无刚体位移。
如图4所示。
6 添加载荷选择距YOZ平面400毫米处的线,加集中载荷力为1860KN,方向为竖直向下,即Y的负向。
加载结果如图5所示。
图5 添加集中力载荷7 设置材料参数和载荷步在Material Models中,设置弹性模量EX为2e5(单位为兆帕),泊松比PRXY 为0.3;材料的密度Density为7.8e-9(单位为千克每立方毫米);摩擦系数Friction coefficient为0.05。
结构力学5平面桁架讲解课件
桁架在动力荷载作用下的响应
瞬态响应
当桁架受到突然施加的动荷载 时,它会表现出瞬态响应。这 种响应通常包括一个短暂的过 渡过程,随后达到一个稳定的 振动状态。
频域响应
在周期性动荷载作用下,桁架 会表现出频域响应。通过频域 分析,可以研究桁架在不同频 率下的振动行为,并确定其振 幅和相位响应。
阻尼效应
高效的经济性
平面桁架能以较少的材料 用量承受较大的荷载,具 有较高的经济性。
平面桁架的应用场景
桥梁工程
在桥梁工程中,平面桁架常被用 作桥面板的支撑结构,能提供稳
定的支撑和承载能力。
建筑工程
在建筑工程中,平面桁架常被用于 楼层和屋盖的承重结构,以及建筑 物的支撑体系。
机械工程
平面桁架也被广泛应用于机械工程 领域,如起重机的梁架、设备的支 架等,其优良的受力性能使其在这 些场景中发挥重要作用。
桁架内力计算:轴力、剪力与弯矩
轴力计算
轴力是杆件沿轴线方向的拉力或压力。通过截面法可以得到杆件的轴力分布情况。根据杆 件的轴力和截面积,可以进一步计算杆件的应力状态,以评估其承载能力。
剪力计算
剪力是杆件横截面上的切向力。通过截面法可以得到杆件的剪力分布情况。剪力的大小和 方向决定了杆件的剪切变形和剪切应力,对于桁架的剪切稳定性分析至关重要。
05 平面桁架的数值模拟与实验验证
基于有限元的数值模拟方法
有限元法基本原理
有限元法将连续体离散为一系列小单元,通过节点连接,利用变分 原理建立节点力与位移的关系,进而求解整个结构的响应。
线性弹性有限元法
对于线弹性材料,采用线性弹性有限元法,通过刚度矩阵和载荷向 量的组装,求解节点位移。
非线性有限元法
02 平面桁架的静力学分析
第9章 桁架和梁的有限元分析
第9章桁架和梁的有限元分析第1节基本知识一、桁架和梁的有限元分析概要1.桁架杆系的有限元分析概要桁架杆系系统的有限元分析问题是工程中最常见的结构形式之一,常用在建筑的屋顶、机械的机架及各类空间网架结构等多种场合。
桁架结构的特点是,所有杆件仅承受轴向力,所有载荷集中作用于节点上。
由于桁架结构具有自然离散的特点,因此可以将其每一根杆件视为一个单元,各杆件之间的交点视为一个节点。
2.梁的有限元分析概要梁的有限元分析问题也是是工程中最常见的结构形式之一,常用在建筑、机械、汽车、工程机械、冶金等多种场合。
梁结构的特点是,梁的横截面均一致,可承受轴向、切向、弯矩等载荷。
根据梁的特点,等截面的梁在进行有限元分析时,需要定义梁的截面形状和尺寸,用创建的直线代替梁,在划分网格结束后,可以显示其实际形状。
二、桁架和梁的常用单元桁架和梁常用的单元类型和用途见表9-1。
通过对桁架和梁进行有限元分析,可得到其在各个方向的位移、应力并可得到应力、位移动画等结果。
第2节 桁架的有限元分析实例一、案例1——2D 桁架的有限元分析图9-1 人字形屋架的示意图 问题人字形屋架的几何尺寸如图9-1所示。
杆件截面尺寸为0.01m 2,试进行静力分析,对人字形屋架进行静力分析,给出变形图和各点的位移及轴向力、轴力图。
条件人字形屋架两端固定,弹性模量为2.0×1011 N/m 2,泊松比为0.3。
解题过程制定分析方案。
材料弹性材料,结构静力分析,属2D 桁架的静力分析问题,选用Link1单元。
建立坐标系及各节点定义如图9-1所示,边界条件为1点和5点固定,6、7、8点各受1000 N 的力作用。
1.ANSYS 分析开始准备工作(1)清空数据库并开始一个新的分析 选取Utility>Menu>File>Clear & Start New ,弹出Clears database and Start New 对话框,单击OK 按钮,弹出Verify 对话框,单击OK 按钮完成清空数据库。
弹性力学与有限元分析第二章-平面桁架有限元分析及程序设计
x
由单元①的刚度方程:
Fj
①
k
① ji
i
①
k
① jj
j
①
k
① ji
2
k
① jj
1
由单元③的刚度方程:
Fj
③
k
③ ji
i
③
k
③ jj
j
③
k
③ ji
3
k
③ jj
1
§2.3 结点平衡与整体刚度矩阵的集成
代入结点1的平衡条件:
k
l
xi
)
(dx j
dxi
)
(
yj
l
yi )
(dy j
dyi )
(dx j dxi ) (dy j dyi )
cos sin
由于杆件的变形产生位移:
ui dxi vi dyi
u j dxj v j dy j
因此,杆件应变为:
dl l
l
(ui
uj)
l
(vi
vj)
杆件轴力为:
(2k1 k2 )v4 P
结构的整体刚度系数
v4
P 2k1
k2
12 3
l2 l1 l1
4 P
N1
N1y
cos
k1v4
cos
k1P
(2k1 k2 ) cos
N2
k2v4
k2P 2k1 k2
位移法求解超静定结构。
§2.1 平面桁架单元的离散
结构的离散化:尽量将结构离散成数量最少的等截面直 杆单元
kki③ ③jii
ki③j
k
③ jj
3 3 3 3
§2.3 结点平衡与整体刚度矩阵的集成
基于ANSYS的平面桁架有限元分析.
PREP7 !* ET,1,LINK180 !* R,1,10, ,0 !* !* MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,2.0e6 MPDATA,PRXY,1,,0.3 WPSTYLE,,,,,,,,0 WPSTYLE,,,,,,,,1 WPSTYLE,,,,,,,,0 WPSTYLE,,,,,,,,1 FLST,3,1,8 FITEM,3,0,0,0 N, ,P51X FLST,3,1,8 FITEM,3,30,0,0 N, ,P51X FLST,3,1,8 FITEM,3,0,30,0 N, ,P51X FLST,3,1,8 FITEM,3,30,30,0 N, ,P51X FLST,3,1,8 FITEM,3,60,30,0
5
数值解与解析解的比较与分析
求出了平面桁架的数值解与解析解,现将两 者的结果进行列表对比
数值解与解析解的比较与分析
表2 整体坐标系下各节点的位移(in)
节点 解析解
U1x 0 0
U1y 0 0
U2x -0.0029 -0.002925
U2y -0.0085 -0.0084404
U3x 0 0
U3y 0 0
基于AN限元分析
平面桁架是工程中常见的结构,本文基于ANSYS平台对平面桁架进行有 限元分析。 首先通过有限元法的理论知识求得平面桁架在一定工况下的理论值,然 后利用ANSYS进行分析得到数值解,最后通过比较理论解与数值解得出结论。 利用ANSYS对平面桁架进行有限元分析,可以提取其他分析结果,对深 入研究平面桁架问题提供了强有力手段,也对其他结构问题的有限元分析具 有指导性意义与价值。
数值解与解析解的比较与分析
表4 单元①的内力与正应力(lb)
有限单元法电子课件(桁架)-PPT精选文档
0 1 0 0 0 0 P 0 0 0 0 0
0 1 0 0 0 0
0 0 0 0 K 33 K 34 K 43 K 44 0 0 K 63 K 64
图3 单元形函数(线性)示意图
平面桁架(Trusses)有限元分析(2)
2、单元应变
u ( u u ) / l [] B {} d x i j e e, x
--- 几何矩阵 [] B [1 / l / l ] e 1 e
3、单元应力
EE [ B ] { d } [ C ] { d }其中E为弹性模量, [C]=E[B] --- 应力矩阵
单元的结点位移
x
i
j
ux ( ) [ N ] { d } e
图2 局部坐标系中的杆单元
N1 N2
u i [ N ] [ N , { d } i N j] e u j
x x N 1 , N i j le le
形函数 (shape function)
1 1
le
1 x 2
▲ 有限元法的要点
●
将连续体(结构)离散为若干子区域
子区域由结点连接为等效的组合体
●
杆系结构
每个单元内假设场变量为多项式(系数不同) 用分区域连续场函数近似全区域的连续场函数
无穷自由度问题转化为有限自由度问题
●
利用变分原理得到离散场变量的大代数方程组 将微分方程边值问题转化为代数方程来求解
连续体
绪 论
x x e e
4、单元刚度矩阵
le Fjle EAe Fl i e EAe
单元的结点力
{ F } [ k] { d } e
平面桁架有限元C#编程
1题目结构如图所示: 杆的弹性模量E 为200000Mpa ,横截面面积A 为3250mm 2。
图 1 桁架示意图2实验材料PC 机一台,Microsoft Visual Studio 软件,Ansys 软件。
3实验原理(1)桁架结构特点桁架结构中的桁架指的是桁架梁,是格一种梁式结构。
桁架结构常用于大跨度的厂房、展览馆、体育馆和桥梁等公共建筑中。
由于大多用于建筑的屋盖结构,桁架通常也被称作屋架。
结构上由光滑铰链连接,载荷只作用于节点处,只约束线位移,各杆只有轴向拉伸和压缩。
(2)平面桁架有限元分析1、单元分析局部坐标系中的干单元如图所示:图 2 局部坐标系中的杆单元以下公式描述了整体位移和局部位移之间的关系:U=Tu 其中U=[ U ix U iy U jx U jy ],T=[cos θ−sin θ00sin θcos θ0000cos θ−sin θ00sin θcos θ],u=[u ix u iy u jx u jy ]U 和u 分别代表整体坐标系和局部坐标系XY 系和局部坐标系xy 下节点i 和节点j 的位移。
T 是变形从局部坐标转换到整体坐标系下的变换阵,类似的局部力和整体力也有以下关系:F=Tf其中F=[ F ixF iy F jx F jy ] ,是整体坐标系下施加在节点i 和j 上的力的分量而且f=[ f ix f iy f jx f jy ],代表局部坐标系下施加在节点i和j上的分量。
在假设的二力杆条件下,杆只能沿着局部坐标系的x方向变形,内力也总是沿着局部坐标系x的方向,因此将y方向的位移设置为0,局部坐标系下内力和位移通过刚度矩阵有如下关系:[f ixf iyf jxf jy]=|k0−k00000−k0k00000|=[U ixU iyU jxU jy]这里k=k eq=AE/L,写成矩阵形式有:f=Ku将f和u替换成F和U有:T-1F=KT-1U将方程两边乘以T得到:F=TKT-1U其中T-1是变换矩阵T的逆矩阵,替换方程中的TKT-1和U矩阵的值,相乘后得到:[F ixF iy F jx F jy]= k[cos2θsinθcosθ−cos2θ−sinθcosθsinθcosθsin2θ−sinθcosθ−sin2θ−cos2θ−sinθcosθcos2θsinθcosθ−sinθcosθ−sin2θsinθcosθsin2θ][U ixU iyU jxU jy]上述方程代表了施加外力、单元刚度矩阵和任意单元节点的整体位移之间的关系。
平面桁架结构的有限元分析
平面桁架结构的有限元分析平面桁架结构是一种经常在建筑和工程领域中使用的结构形式。
它由直杆组成,连接在节点上,形成一个稳定的平面结构。
平面桁架结构的设计和分析需要使用有限元分析方法来确定结构的受力状态和稳定性。
本文将介绍平面桁架结构的有限元分析方法,包括模型建立、加载条件、应力和变形分析等。
首先,建立平面桁架结构的有限元模型。
模型应包括杆件和节点两个基本元素。
杆件是结构的主要受力元素,节点是杆件的连接点。
通过连接节点和杆件,可以构建起整个桁架结构。
在有限元模型中,每个节点被赋予一个坐标,每个杆件的长度和截面积也需要定义。
通过这些信息,可以建立结构的有限元模型。
加载条件是进行有限元分析的第二个关键步骤。
加载条件包括结构所承受的外部力和约束条件。
外部力是指作用于结构上的力,包括重力、风力、地震力等。
约束条件是指限制结构自由运动的条件,例如固定节点或滑动支座等。
在有限元分析中,将这些加载条件应用到有限元模型中,以模拟真实结构的受力情况。
然后进行应力和变形分析。
在有限元分析中,结构的应力分布和变形情况可以通过求解有限元方程来得到。
有限元方程是由结构的力平衡和材料的应力-应变关系所组成的方程组。
通过求解有限元方程,可以计算出结构中每个节点的应力和变形情况。
这些结果可以用来评估结构的安全性和稳定性。
在进行有限元分析时,需要注意一些细节。
首先,选择合适的材料模型和参数。
不同的材料具有不同的力学特性,例如弹性模量、屈服强度等。
选择适当的材料模型和参数,以获得准确的分析结果。
其次,进行网格划分和单元类型选择。
将结构划分为小单元,并选择适当的单元类型,以确保每个单元的形状和大小适合结构的几何形状。
最后,进行后处理和结果分析。
得到应力和变形结果后,可以进行结果的可视化和分析,以评估结构的性能。
总之,平面桁架结构的有限元分析是一种有效的工具,可以用于评估结构的受力状态和稳定性。
通过合适的模型建立、加载条件选择以及应力和变形分析等步骤,可以得到准确的分析结果,为结构的设计和优化提供有力支持。
平面桁架的有限元法
Kz=Table[0, {i, 2nj}, {j,2nj}]; “开总刚度矩阵, nj 总节点数” ;
For[e=1, e<=ne, e++, For[i=1, i<=2, i++,
ke T ke T rans“p生os成e[T单] ;刚,变坐标系” ;
For[ii=1, ii<=2, ii++, r =2(i-1)+ii; rr=2(jm[[e, i+1]]-1)+ii;
b
ui i vi
o
x
xi 0, xj b
ui 1 0 0 0a1
vi
u j
v j
0 1 0
0 b 0
1 0 1
0 0 b
aa32 a4
ui 1 0 0 0a1
vi
u j
v j
0 1 0
0 b 0
1 0 1
0 0 b
aa32 a4
{ e} [ Ab ]{a}
解线性代数方程组,得
代入 {a} [ Ab ]1{ e}
{ f } [Hs ]{a}
{ f得}21 [Hs ]24[ Ab ]414{ e}41
a1
u 1
v
0
x 0
0 1
0 x
aa32
a4
{f
}21
[N
f
]24{
}e 41
节点位移与单元内位移的关
系
{ f } [N f ]{ e}
{ e} [T ]{ e}
[T
]
t 0
0
t
[t]1 [t]T [T ]1 [T ]T
[T ]{Re} [k ]{ e}
平面桁架程序计算原理及程序编制
求解方法
力法:以力未知数,必须预先满足平衡条件,然 后通过连续条件求解未知力;超静定基的选取。 位移法:以位移为未知数,各杆件的变形由相连 接的节点位移确定(变形协调条件),通过各个 节点的平衡方程求出未知位移,再由位移计算出 各杆件的内力;各节点的平衡方程也可由最小总 势能原理推导。
以平面桁架结构分析的程序设计为例,介绍结构分 析和程序设计的方法。
最小总势能原理与位移法都是以位移为未知 数使变形状态预先满足连续条件。现对上述 例题采用最小总势能原理进行求解。
• 总势能由两部分组成
➢结构的弹性应变能 U ➢外力由于结构变位所产生的势能 V
由于桁架结构的杆件只能承受拉压力,所以 单根杆件的应变能为
U EF(L)2 / L 2
对于整个桁架应变能是所有杆件应变能的叠 加,即
v u
1 2
w 4 v 2
现对各位移变量分别取偏微商后,得
E a Fg1212u1212v1u20P E a Fg 212u 1 1212 v100 P
E a Fg u10 1212 u2212v2 P
E a Fg00212u21212v2P
注:值得指出的是刚度矩阵中的系数只与结构本身的几何形态和约 束条件有关,而与外载无关。
4
3
平面桁架结构
将 U 和 V 相加,得到总势能。
由于V 是位移的一次函数,总势能就成为位 移的二次非齐次函数。
根据最小总势能原理,在所有可能的位移状态中, 真正发生的位移状态使总势能为最小。即函数对 自变量的偏微商为零,即
式中
(UV)0 (i1,2,3,4) wi
w1 u1
W
w2 w3
5-1 矩阵位移法
• 桁架是由离散杆件组成的构架结构,杆件的端点 借助于无摩擦的铰连接起来。
简单平面桁架的稳定性及极限承载基于有限元的几何非线性分析
 ̄I i I I II I , 杆件 处 于平衡 状态 ,  ̄ I i 1 8兀 =0, 由此得 到 l l =0, 即
=
f ( ” ) 班一 P 喀 ) 班 。
( 1 1 )
( 1 2 )
状态时, 总势能的一阶变分为零 , 或此体系的总势能为驻值 。 这就是势能驻值
原理 。
曰髓四圈
建筑理论与设计
简单平面桁架 的稳定性及 极限承 载基于有限元的几何 非线性分析
王国权 张 宝勤
北 京 市 住 宅 建 筑设 计 研 究院 有 限 公 司 中 国 电子 工 程 设 计 院
摘要 : 本 文 首先 介 绍 了分 析压 杆 稳定 的线 性有 限 单元 , 然 后 引入 对 几何 非 线性 问题 的 一般 讨 论 , 提 出 了结 构 的切 线 刚度 矩 阵。 最 后介 绍 了求 解几 何 非线 性方 程 的牛 顿—— 拉 斐逊 法 , 应 用 此方 法通 过 编程 分 析得 出简 单桁 架 荷载一 位 移 曲线 。
一
横 向 荷 载 所 作 的 功 = f E q ( ) 主 q ( ) ] = q f g ( ) ( x ) = q ( 9 )
虚 位移 是满 足 体系 支承 约束 条 件下 的一 个 微小 位 移变 化 , 是 实际 位移 的 阶变 分 , 因此 虚应 变 能 8 u 就是 实 际应 变 能 的一 阶变 分 , 一8 We 就 是 实 际外
( U — w) = O ( 3 1
( 1 0 )
力势能 ( 一 We ) 的 一 阶变 分 , 简 写 为 一8 We 。故 式 ( h ) 可以写成 : 8( u + v ) :8 即 8 l I : 0 , 式 中 Ⅱ= u + V = u — W为 体系 具有 的总势 能 。当体 系处 于 在平 衡
有限元分析(桁架结构)
有限元上机分析报告~学院:机械工程专业及班级:机械设计及其自动化08级7班姓名:***学号:题目编号: 2》1.题目概况结构组成和基本数据结构:该结构为一个六根杆组成的桁架结构,其中四根杆组成了直径为800cm的正方形,其他两根杆的两节点为四边形的四个角。
材料:该六根杆截面面积均为100cm2,材料均为Q235,弹性模量为200GPa,对于直径或厚度大于100mm的截面其强度设计值为190Mpa。
载荷:结构的左上和左下角被铰接固定,限制了其在平面内x和y方向的位移,右上角受到大小为2000KN的集中载荷。
结构的整体状况如下图所示:分析任务】该分析的任务是对该结构的静强度进行校核分析以验算该结构否满足强度要求。
2.模型建立物理模型简化及其分析由于该结构为桁架结构,故认为每根杆件只会沿着轴线进行拉压,而不会发生弯曲和扭转等变形。
结构中每根杆为铰接连接,有集中载荷作用于最上方的杆和最右方杆的铰接点。
单元选择及其分析由于该结构的杆可以认为是只受拉压的杆件,故可以使用LINK180单元,该单元是有着广泛工程应用的杆单元,它可以用来模拟桁架、缆索、连杆、弹簧等等。
这种三维杆单元是杆轴方向的拉压单元,每个节点具有三个自由度:沿节点坐标系X、Y、Z方向的平动。
就像铰接结构一样,不承受弯矩。
输入的数据有:两个节点、横截面面积(AREA)、单位长度的质量(ADDMAS)及材料属性。
输出有:单元节点位移、节点的应力应变等等。
由此可见,LINK180单元适用于该结构的分析。
模型建立及网格划分((1)启动Ansys软件,选择Preferences→Structural,即将其他非结构菜单过滤掉。
(2)选择单元类型:选择Preprocessor→Element Type→Add/Edit/Delete→Add,在出现的对话框中选择Link→3d finit stn 180,即LINK180,点击“OK”(3)选择实常数:选择Preprocessor→Real Constants→Add/Edit/Delete→Add,在出现的对话框中的Cross-sectional area中输入100,点击“OK”。
有限元原理 结构矩阵分析(平面桁架 平面应力) 变分
设平面桁架单元在总体坐标系中刚度矩阵的一般形式为
由(2-1-8),当单元结点位移为{1 0 0 0 }T时,在单元各结点上施加的力刚好为单元刚度矩阵中的第一列:{k11k21k31k41}T。对[k]的其他各列也可做出类似的解释。即单元刚度矩阵的每一列相当于一组特定位移下的结点力,如表2-1所示。由图2-4可以获得更为直观的理解。
它们将作为程序的输入数据(几何参数)。
每个结点有两个自由度,对结点1、2、3分别为
若暂时不考虑支承约束条件,整个结构的结点自由度为
3、单元分析(建立结点力与结点位移之间的关系)
取一个一般性的单元,设它的两个结点在结构中的编号为i, j(单元内部的结点序号)。由材料力学可知,杆的轴向刚度为EA/L。其中L为杆的长度:
单元结点自由度{u}={uiviujvj}T
结点给单元的力{r}={piqipjqj}T
在图2-3中,x’轴与x轴的夹角为α
结点的位移分量的坐标变换为
单元的位移分量的坐标变换为
或缩写为
类似,{r’}与{r}之间的转换关系为
由于
是正交矩阵,因此
也是正交矩阵。所以有
将(2-1-4)、(2-1-5)代入(2-1-2)有
为了描述结构的平衡需要建立一个坐标系,称为总体坐标系,以区别于以后出现的“局部坐标系”。总体坐标系的选择原则上不受限制,但希望使用方便。本节所选的总体坐标系示于图2-2,坐标原点与结点1重合。以u, v分别表示沿x, y方向的位移分量,p, q分别表示力沿x, y轴的力分量(投影)。
在总体坐标系中各结点的坐标为:
对结点1:
对结点2:
对结点3;
可以合并成
式(2-1-14)的右边为外载荷和支反力。左边则为单元给结点的力,它们是未知的,但可以借助单元刚度矩阵以结点位移来表示。
基于ANSYS分析的平面桁架结构优化设计
I v si a i n o t t ・ f t ・ r f fr ・ e it nc e i n f r c n r t t u t r s n e tg to n sa e o -he a to ie r s s a e d s g o o c研究 的努力方 向。
火灾作为一种多发 的灾害 , 对人们的生命及财产造成惨重的
3 3 混 凝 土 结 构 抗 火 设 计 方 法 的研 究 .
设 想混凝 土结构的抗火设计可从两个途径进行研究 :) 1把火 损失。建筑火灾对混凝土结构造成一定 的损伤甚至整体 的破坏。 建立混 凝土结 构的抗 火设计方 法 , 灾的高温作 用等效 为一 种荷 载 , 结构 上 的其他 荷载 ( 与 恒载 、 活 研究混凝 土结构 的抗火性 能 , 理应成为混凝 载、 风载 、 地震作用等 ) 一起 参与 荷载效 应组合 , 概率极 限状态 建立抗火混凝 土结构 的损伤评估及修 复加 固方法 , 按 土结构研究的一项重要 任务 。建 立我 国的混凝 土结构抗 火设计 设计方法进行设计 , 即建立考虑火灾 高温作用 的统 一的结构设计
NS S优化 分析 的基本概 念 计相比 , 优化设计 可以使土建工 程降低造价 5 %~3 %。2 0 0世纪 2 有 关 A Y
6 0年代以来 , 随着计算 机计算能 力的不断 提高 , 们把有 限元分 人 和成熟 的方 法 , 使得结构优化 的技术得到 了更快的发展。 A Y NS S优化分析中包括 的基 本概念有设 计变量 、 状态变量 、 1设计变 量是作为 自变量 , 过改变设 计变量的数值来实现 ) 通 析的方 法和各种数学规划方法相结合 , 逐步发展成为 一种系统 目标 函数 、 并 分析文件等。 文中以六杆平面桁架 为例 , 利用 AN Y S S的优化 分析功 能对 结 果的优化 , 设计变量的上下限决定了设计 变量 的变化范围。 坏可能引起结构 的连续倒塌和整体破坏 。研 究火灾高温下 , 同 加 固。 目前 , 不 对现 有建筑结构 加 固方 法的研究非 常活跃 , 充分研 借助 已有 的加 固方法 和手段 , 应 结构 的性能变化规律 ; 研究 火灾 高温下 , 结构 连续倒 塌 和整 体破 究混凝土结构 的火灾 损伤 特点 , 坏的机理 , 是结构抗火研究 的主要 内容 。
2D四杆桁架结构的有限元分析实例学习资料
2D四杆桁架结构的有限元分析实例实例:2D四杆桁架结构的有限元分析学习有限元方法的一个最佳途径,就是在充分掌握基本概念的基础上亲自编写有限元分析程序,这就需要一个良好的编程环境或平台。
MATLAB软件就是这样一个平台,它以功能强大、编程逻辑直观、使用方便见长。
将提供有限元分析中主要单元完整的MATLAB程序,并给出详细的说明。
1D杆单元的有限元分析MATLAB程序(Bar1D2Node)最简单的线性杆单元的程序应该包括单元刚度矩阵、单元组装、单元应力等几个基本计算程序。
下面给出编写的线性杆单元的四个MATLAB函数。
Bar1D2Node _Stiffness(E,A,L)该函数计算单元的刚度矩阵,输入弹性模量E,横截面积A和长度L,输出单元刚度矩阵k(2×2)。
Bar1D2Node _Assembly(KK,k,i,j)该函数进行单元刚度矩阵的组装,输入单元刚度矩阵k,单元的节点编号i、j,输出整体刚度矩阵KK。
Bar1D2Node _Stress(k,u,A)该函数计算单元的应力,输入单元刚度矩阵k、单元的位移列阵u(2×1)以及横截面积A计算单元应力矢量,输出单元应力stress。
Bar1D2Node_Force(k,u)收集于网络,如有侵权请联系管理员删除该函数计算单元节点力矢量,输入单元刚度矩阵k和单元的位移列阵u(2×1),输出2×1的单元节点力矢量forces。
基于1D杆单元的有限元分析的基本公式,写出具体实现以上每个函数的MATLAB程序如下。
%%%%%%%%%%% Bar1D2Node %% begin %%%%%%%%%function k=Bar1D2Node_Stiffness(E, A, L)%该函数计算单元的刚度矩阵%输入弹性模量E,横截面积A和长度L%输出单元刚度矩阵k(2×2)%---------------------------------------k=[E*A/L -E*A/L; -E*A/L E*A/L];%%%%%%%%%%%%%%%%%%%%%%%%%%function z=Bar1D2Node_Assembly(KK,k,i,j)%该函数进行单元刚度矩阵的组装%输入单元刚度矩阵k,单元的节点编号i、j%输出整体刚度矩阵KK%-----------------------------------DOF(1)=i;DOF(2)=j;for n1=1:2for n2=1:2收集于网络,如有侵权请联系管理员删除KK(DOF(n1), DOF(n2))= KK(DOF(n1), DOF(n2))+k(n1, n2);endendz=KK;%------------------------------------------------------------function stress=Bar1D2Node_Stress(k, u, A)%该函数计算单元的应力%输入单元刚度矩阵k, 单元的位移列阵u(2×1)%输入横截面积A计算单元应力矢量%输出单元应力stress%-----------------------------------stress=k*u/A;%-----------------------------------------------------------%%%%%%%%%%%%%%%%%%%%%%%%%function forces=Bar1D2Node_Force(k, u)%该函数计算单元节点力矢量%输入单元刚度矩阵k和单元的位移列阵u(2×1)%输出2×1的单元节点力分量forces%-----------------------------------------forces=k*u;%%%%%%%%%%% Bar1D2Node %% end %%%%%%%%%收集于网络,如有侵权请联系管理员删除【四杆桁架结构的有限元分析—数学推导】如图所示的结构,各杆的弹性模量和横截面积都为E=29.54×10N/mm2,A=100mm 2,试求解该结构的节点位移、单元应力以及支反力。
平面桁架ANSYS有限元法分析实例
2. 前处理 (1)定义单位
从第二章可知,ANSYS中单位可以不定义,但建模时一定要 保证单位的一致。
已知:各杆的弹性模量E=2.0×105MPa,各杆截面均为A=0.5cm2,杆13长 为100cm,载荷P=2KN,试求平面桁架的内力和位移。
本题采用单位m-kg-s-N较简便,建模过程中 的所有参数都选用m-kg-s-N,相应计算结果 应力为Pa。
改为国际单位制:各杆的弹性模量E=2.0×1011Pa, 各杆截面均为A=0.5e-4m2,杆13长为1m,载荷 P=2000N。
(2)定义单元类型
单元类型
特点
结点数 结点自由度
适用
LINK1 LINK8 LINK10
二维杆单元,只承受 轴向的拉压力,不考 虑弯矩
三维杆单元,具有塑 性、蠕变、膨胀、应 力刚化、大变形、大 应变等功能。
平面桁架ANSYS有限元法分析实例
例3-1 设平面三角结构的桁架123如 图3-4所示。已知:各杆的弹性模量 E=2.0×105MPa,各杆截面均为 A=0.5cm2,杆13长为100cm,载荷P=2KN, 试求平面桁架的内力和位移。
解:传统分析方法
设杆12、杆23和杆13的内力分别为N1、N2和N3。在总体坐标系 x-y(或U-V)中,由力的平衡方程可以得到结点的内力值。
3.求解 (1)施加约束
• 本例中,点1为固定支座,点3为活动支座。 • 在节点1上,约束UX、UY; • 在节点3上,约束UY。
• 在节点1上,约束UX、UY,如图; • 在节点3上,约束UY。
(2)施加载荷
选节点2,按图示完成;
•apply-,选FY,输入-2000,OK。 施加载荷后,结果如图
仅受拉或受压的三维 杆单元,具有应力刚 化和大变形功能。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
0 1 0 ui ui v e vi 0 0 0 i u k u 0 1 0 j j vj 0 0 0 vj
单元轴力:
N
AE 1 0 l
ui ui v v i 1 0 S i u u j j v j v j
1杆和3杆位移:
N1l1 14 34 E1 A1
N 2l2 24 E2 A2
P
第二章 平面桁架有限元分析及程序设计
超静定桁架
1杆轴力竖向分量:
E1 A1 E1 A1 cos2 N1 y N1 cos 14 cos v4 k1v4 l1 l1
2杆轴力:
第二章 平面桁架有限元分析及程序设计
超静定桁架
代入平衡方程:
2N1y N2 y P
(2k1 k2 )v4 P
P v4 2k1 k2
结构的整体刚度系数
1
l2
2
l1 l1
3
4 P
k1v4 k1P N1 cos cos (2k1 k2 ) cos
k2 P N 2 k2v4 2k1 k2
因此:
e
0 0
0 0
e
0 U i V 0 i U j V j
cos
sin
F T F
0 0
其中,[T]为转换矩阵:
转换矩阵的性质
T 0 0 0 0
式中:
S AE
l
cos
sin
§2.2 平面桁架的单元分析
§2.2.3 单元坐标转换矩阵
取任意杆件,建立如图所示的局部坐标系:
杆端位移:
杆端力:
ui Ui
vi u j Vi U j
vj Vj
x
Vj
y
Ui
F
e
vj uj
Uj
U i V i U j V j 1 AE 0 l 1 0 0 1 0 ui v 0 0 0 i 0 1 0 u j v 0 0 0 j
ui ui AE 1 1 S N l u j u j
式中 S 为单元应力(广义)矩阵;
AE 1 1 S l
式中 k
为单元刚度矩阵(局部坐标系)
e
§2.2 平面桁架的单元分析
杆端位移:
杆端力:
ui
vi v j 0
有限单元法及程序设计
第二章 平面桁架有限元分析及程序设计
§2.1 平面桁架单元的离散 §2.2 平面桁架单元分析 §2.3 结点平衡与整体刚度矩阵的集成 §2.4 边界条件的处理 §2.5 单元内力与支座反力的计算 §2.6 平面桁架有限元程序设计
第二章 平面桁架有限元分析及程序设计
回顾
静定桁架 B a A
式中 k 是整体坐标系下的单元刚度矩阵;
§2.2 平面桁架的单元分析
写成分块矩阵形式:
U i 2 V i AE 2 U l j V j
F e
2 2
2
i
e
j
Uj
符号:与坐标系的方向一致为正,反之为负。
§2.2 平面桁架的单元分析
杆的受力分为两种情况:
右结点固定 结点位移: 单元应变: 单元应力: 左结点固定
ui ui
ui l
uj 0
ui 0
uj uj
uj l
E ui l
E
U i A
E
§2.2 平面桁架的单元分析
§2.2.2 整体坐标系下的单元刚度矩阵
若局部坐标系与整体坐标系重合,则整体坐标系下的单元刚度矩 阵与局部坐标下的单元刚度矩阵相同。 若局部坐标系与整体坐标系不重合,如下图所示:
Vj
y
vj
N
杆端位移和杆端力 i 结点:
杆端位移: 杆端力:
uj U j
Ui ui
N
vi Vi
Eu j l
AE uj l
单元左端杆端力:
单元右端杆端力:
AE ui l AE ui l
U i A U j A
U j A
AE uj l
§2.2 平面桁架的单元分析
任意情况(左右结点均有变形)即为以上两种状态的叠加: 杆端力为:
AE AE U u uj i i l l AE AE U j ui uj l l
u j j v j
U i 式中: Fi Vi
U j F j V j
kii AE l
2
2
k
jj
AE 2 l 2
§2.2 平面桁架的单元分析
Ⅰ 1
D
FP
a
Ⅰa
a
C FP
解题方法
方法1:节点法 方法2:截面法
第二章 平面桁架有限元分析及程序设计
超静定桁架
解题方法:力法和位移法
1
l2
2
l1 l1
3
如图所示桁架,求各杆轴力。
力的平衡条件:
2 N1 cos N 2 P
位移的协调方程:
4 P
N1 N 2 2杆位移: N1
24 14 cos
dl ( x j xi ) l l (dxj dxi ) (dyj dyi ) (dx j dxi ) ( y j yi ) (dy j dyi )
cos sin
由于杆件的变形产生位移:
ui dxi
dl (ui u j ) (vi v j ) 因此,杆件应变为: l l l AE N EA [ (ui u j ) (vi v j )] 杆件轴力为: l
i
e
j
uj
Ui
Vi V j 0
i
e
j
Uj
单元杆端力方程:
U i V 0 i U j V j 0
AE AE ui uj l l AE AE ui uj l l
U i V i U j V j
1 AE 0 l 1 0
k k
ij ji
AE 2 l
2
单元刚度矩阵的性质
(1)单元刚度系数kij的意义 j自由度(结点)产生的单位杆端位移引起的i自由度(结点)的杆端力
(2)单元刚度矩阵是对称矩阵 (3)单元刚度矩阵一般是不可逆的
反力互等定理
杆件单元的应力矩阵为:
N2 y N2
E2 A2 E A 24 2 2 v4 k2v4 l2 l2
V4为第4节点竖向位移
k1 和 k2 为杆件的刚度系数; 式中:
E1 A1 cos2 k1 l1 E2 A2 k2 l2
物理意义: 4点产生单位位移,杆端产生的竖向杆端力; 由杆件的物理性质和几何性质决定;
9个单元,6个结点
16个单元,8个结点
§2.2 平面桁架的单元分析
§2.2.1 局部坐标系下的单元刚度矩阵 局部坐标系的建立
j
x
y
i
y
E,A,l e i
e
杆端位移:
j
x
ui
■ 原点:以第一个结点为坐标原点;
■ x 轴:沿单元的杆轴方向; ■ y 轴:从 x 轴逆时针旋转90°。
i
e
j
uj
杆端力:
Ui
0 0
T T T I T T T 1
转换矩阵是正交矩阵;
§2.2 平面桁架的单元分析
同理,位移也存在转换关系:
代入局部坐标系下 的刚度方程:
e
e
T F T F k k T
位移法求解超静定结构。
N1 y
§2.1 平面桁架单元的离散
结构的离散化:尽量将结构离散成数量最少的等截面直 杆单元 离散原则:每个结点离散后还是一个结点,每个杆件离 散后变成一个单元 4 ③
⑤ ⑥
5 ④
⑦ ⑧ ②
6
⑨
5 ④
⑨ ⑦ ⑧
6
⑩
11
⑤
12
7 ⑥
13 14 15
8
16
1①
2
3
1 ①
2
②
3
③
4
2
ui 2 v i u j 2 v j
F
e
Fi kii Fj k ji
kij i k jj j ui i vi
V j N sin N
因此,杆件结点力向量为:
U i 2 V i AE 2 U l j V j
e
F e
2 2
2
2
ui 2 vi e e k u j 2 v j
ui Ui
vi
x
j 结点:
杆端位移:
杆端力:
Vi
uj Uj
vj Vj
符号:与坐标系的方向一致为正,反之为负。