杆件结构有限元分析课件.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《杆件结构有限元分析课件.ppt》由用户(三亚风情)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 结构 有限元分析 课件
- 资源描述:
-
1、第三讲第三讲杆件结构有限元分析杆件结构有限元分析元计算技术部 本讲通过对一维杆件进行理论分析和ELAB1.0有限元实现两种方式来介绍有限元的思想,并给出用有限元方法求解实际问题的流程:u一维杆件有限元理论分析u一维杆件ELAB软件实现基本方程基本方程的最终弱形式建立有限元模型 确定杆单元的形函数确定单元刚度,总体刚度带入边界求解方程工程建模前处理求解计算及后处理 杆件从构造上说是长度远大于其截面尺寸的一维构件。在结构力学中常常将承受轴力或扭矩的杆件统称为杆。根据桁架或者杆件承受的载荷不同,根据它们的变形不同可以将杆单元分为二维和三维桁架单元两种。承受轴向载荷等截面直杆的基本方程如下:0)(xf
2、dxdAx0()()0lxdAf xudxdx 其弱形式为:00llxjjjd uAdxf xudxP udx基本方程基本方程u一维杆件有限元理论分析一维杆件有限元理论分析 根据虚功原理,方程两边乘以虚位移u,平衡方程可以写为:其中,右端最后一项可以看作是节点力情况,所以可以不单独列出,同时 xduEdx所以上式可以继续写为:00lldu d uAEdxf xudxdx dx其中E表示弹性模量,A表示横截面积,方程左端得到单元的刚度矩阵。基本方程的最终弱形式基本方程的最终弱形式建立有限元模型建立有限元模型 现考虑一个由5个长度相同(le=1m)横截面积不同的杆件构成的一维杆件,各杆弹性模量都为
3、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左,分析杆件内位移分布:几何模型 将其划分为五个单元六个节点,即每根杆件作为一个单元,每个单元的节点关系如下图所示:单元拓扑关系确定杆单元的形函数确定杆单元的形函数考虑其中一个杆单元,其两个端点分别为节点1,节点2,基本变量为节点位移u1,u2::设该单元的位移场为u(x),这里要通过一个线性分布函数来进行插值以得到单元内的未知位移。定义一个局部坐标系:1212()1xxxx由该式可得,在节点x1 处,=-1,在节点x2处,
4、=1,如下图所示局部坐标系:基于上坐标系来定义用于位移插值的形状函数:121()21()2NN确定了形状函数后,单元内的线性位移场就可以用节点位移u1,u2来表示:1 122()u xN uN u通常来说形函数需满足以下条件:通常来说形函数需满足以下条件:在单元内,一阶导数必须存在;在单元之间的连接处,位移必须连续;推导该问题的单元刚度矩阵表达式推导该问题的单元刚度矩阵表达式对于基本方程的最终弱形式:00lldu d uAEdxf xudxdx dx在一个单元内:1 122()u xN uN u因此:1 1221212()d N uN udNdNduuudxdxdxdx由局部坐标系,与整体坐标
5、系的关系式及位移的插值函数可得:1121212221211211212112eedNdN ddxddxxxxxldNdN ddxddxxxxxl 21elxx11212211eeudNdNduuuudxdxdxll 11120121122111121111leeeeeeeeeeluldu d uAEdxu uAEdudx dxlllluu uAElulll由上面的关系式可得到:对u取与u相同的形函数,并将上面的关系式带入基本方程的最终弱形式中可得:确定该问题每个单元的刚度矩阵确定该问题每个单元的刚度矩阵222211111111111 1eeeeeeeeeeeelllAEkAElAEllllll
展开阅读全文