当前位置:文档之家› 有限元原理 结构矩阵分析(平面桁架 平面应力) 变分

有限元原理 结构矩阵分析(平面桁架 平面应力) 变分

有限元原理 结构矩阵分析(平面桁架 平面应力) 变分
有限元原理 结构矩阵分析(平面桁架 平面应力) 变分

ansys桁架和梁的有限元分析

ansys桁架和梁的有限元分析

————————————————————————————————作者:————————————————————————————————日期:

桁架和梁的有限元分析 第一节基本知识 一、桁架和粱的有限元分析概要 1.桁架杆系的有限元分析概要 桁架杆系系统的有限元分析问题是工程中晕常见的结构形式之一,常用在建筑的屋顶、机械的机架及各类空间网架结构等多种场合。 桁架结构的特点是,所有杆件仅承受轴向力,所有载荷集中作用于节点上。由于桁架结构具有自然离散的特点,因此可以将其每一根杆件视为一个单元,各杆件之间的交点视为一个节点。 2.梁的有限元分析概要 梁的有限元分析问题也是是工程中最常见的结构形式之一,常用在建筑、机械、汽车、工程机械、冶金等多种场合。 梁结构的特点是,梁的横截面均一致,可承受轴向、切向、弯矩等载荷。根据梁的特点,等截面的梁在进行有限元分析时,需要定义梁的截面形状和尺寸,用创建的直线代替梁,在划分网格结束后,可以显示其实际形状。 二、桁架和梁的常用单元 桁架和梁常用的单元类型和用途见表7-1。 通过对桁架和粱进行有限元分析,可得到其在各个方向的位移、应力并可得到应力、位移动画等结果。 第128页

第二节桁架的有限元分析实例案例1--2D桁架的有限元分析 问题 人字形屋架的几何尺寸如图7—1所示。杆件截面尺寸为0.01m^2,试进行静力分析,对人字形屋架进行静力分析,给出变形图和各点的位移及轴向力、轴力图。 条件 人字形屋架两端固定,弹性模量为2.0x10^11N/m^2,泊松比为0.3。 解题过程 制定分析方案。材料为弹性材料,结构静力分析,属21)桁架的静力分析问题,选用Link1单元。建立坐标系及各节点定义如图7-1所示,边界条件为1点和5点固定,6、7、8点各受1000N的力作用。 1.ANSYS分析开始准备工作 (1)清空数据库并开始一个新的分析选取Utility Menu>File>Clear&Start New,弹出Clears database and Start New对话框,单击OK按钮,弹出Verify对话框,单击OK按钮完成清空数据库。 (2)指定新的工作文件名指定工作文件名。选取Utility Menu>File>Change Jobname,弹出Change Jobname对话框,在Enter New Jobname项输入工作文件名,本例中输入的工作文件名为“2D-spar”,单击OK按钮完成工作文件名的定义。 (3)指定新的标题指定分析标题。选取Ufility Menu>File>Change Title,弹出ChangeTitle对话框,在Enter New Tifie项输入标题名,本例中输入“2D-spar problem'’为标题名,然后单击OK按钮完成分析标题的定义。 (4)重新刷新图形窗9 选取Utility Menu>Plot>Replot,定义的信息显示在图形窗口中。 (5)定义结构分析运行主菜单Main Menu>Preferences,出现偏好设置对话框,赋值分析模块为Structure结构分析,单击OK按钮完成分析类型的定义。 2.定义单元类型 运行主菜单Main Menu>Preprocessor>Element Type>Add/Edit/Delete命令,弹出Element Types对话框,单击Add按钮新建单元类型,弹出Library of Element Types对话框,先选择

c语言计算平面桁架内力计算程序

#include #include #include #define M 5 int n,nc,nn,m,e,f;//节点总数,固定节点数,自由度数,杆件数int io,jo;//单根杆对号指示数 int ihl[M],ihr[M];//杆件左右节点号 double a[M];//各杆截面积 double mm[M];//杆件质量 double ea[M];//杆件EA的值 double x[M],y[M];//节点坐标 double dp[M];//总体系下的节点载荷 double t[2];//0,1分别为坐标转换矩阵的cos(),sin() double c[2][2];//总体系下的单刚 double clxy[3];//0,1,2分别为杆长,正弦,余弦 double h[M];//杆件轴力 double r[M][M];//总刚度阵 double rd;//桁架轴力杆局部系单刚 double u[M];//桁架节点位移 double v[2];//存放节点位移差 double d[M];//LDLT分解时的D矩阵的对角线元素 double l[M][M];////LDLT分解时的D矩阵的对角线元素double fdp[M];//总体系下支座反力 void iojo(int k)//计算对号指示数io,jo { int i,j; i=ihl[k-1];//k号杆左节点号进入i j=ihr[k-1];//k号杆节点右号进入i io=2*(i-nc-1);//uxi前未知位移的个数 jo=2*(j-nc-1);//uyi前未知位移的个数 } void ch(int k)//计算杆长与方向余弦函数 { int i,j; i=ihl[k-1];//k号杆左节点进入i j=ihr[k-1];//k号杆右节点进入j clxy[1]=x[j-1]-x[i-1];//k号杆x坐标差 clxy[2]=y[j-1]-y[i-1];//k号杆y坐标差 clxy[0]=sqrt(clxy[1]*clxy[1]+clxy[2]*clxy[2]);//k号杆长 clxy[1]=clxy[1]/clxy[0];//k号杆件x轴余弦 clxy[2]=clxy[2]/clxy[0];//k号杆件y轴余弦

有限元分析 均布荷载作用下深梁的变形和应力

有 限 元 分 析 上 级 报 告 学院: 专业: 姓名: 班级: 学号:

均布荷载作用下深梁的变形和应力 两端简支,长度l=5m,高度h=1m的深梁,在均布荷载q =5000N/m作用下发生平面弯曲(如图4.1所示)。已知弹性模量为30Gpa,泊松比为0.3,试利用平面应力单元PLANE82,确定跨中的最大挠度,和上下边缘的最大拉压应力。 4.1 均布荷载作用下深梁计算模型 1.理论解 具有两个简支支座支承的简支梁,它的变形和应力分布在理论上是没有解析表达式。 在一般的弹性力学教科书中,只有将两边支座简化为等效力的条件,即在两个支座的侧表面上作用有均匀分布的剪力情况,才可以得到理论解答。 (1) 设定应力函数。 获得这种情况下的解答的主要思路是:按照应力解法,考虑到应力分量关于该梁中心 位置(x=2.5,y=0.5)有对称和反对称关系。可以首先假定一个应力函数为: Φ = A(y - 0.5)5+ B(x - 2.5)2 (y -0.5)3 +C(y -0.5)3+ D(x- 2.5)2+ E(x -2.5)2 (y - 0.5) (4.1) 依据这个应力函数,可以获得各个应力分量,按照上表面受均布压力作用简支梁的上 下表面和左右侧表面的应力边界条件,确定出应力函数(4.1)中的各个待定系数A,B,C,D和E。 按照应力求解平面应力问题方法,应力函数应该满足双调和函数: ?2?2Φ = 0 (4.2) 将(4.1)应力函数代入上式后,得到: 24 B( y - 0.5) +120A(y - 0.5) = 0 (4.3) 即: B = -5A (4.4) (2)确定应力分量。 应力函数与应力分量之间的关系为: (3) 利用梁的上下表面边界条件确定积分常数。 上表面受均布压力作用简支梁的上表面(y=h=1m)的应力边界条件:

第9章 桁架和梁的有限元分析

第9章桁架和梁的有限元分析 第1节基本知识 一、桁架和梁的有限元分析概要 1.桁架杆系的有限元分析概要 桁架杆系系统的有限元分析问题是工程中最常见的结构形式之一,常用在建筑的屋顶、机械的机架及各类空间网架结构等多种场合。 桁架结构的特点是,所有杆件仅承受轴向力,所有载荷集中作用于节点上。由于桁架结构具有自然离散的特点,因此可以将其每一根杆件视为一个单元,各杆件之间的交点视为一个节点。 2.梁的有限元分析概要 梁的有限元分析问题也是是工程中最常见的结构形式之一,常用在建筑、机械、汽车、工程机械、冶金等多种场合。 梁结构的特点是,梁的横截面均一致,可承受轴向、切向、弯矩等载荷。根据梁的特点,等截面的梁在进行有限元分析时,需要定义梁的截面形状和尺寸,用创建的直线代替梁,在划分网格结束后,可以显示其实际形状。 二、桁架和梁的常用单元 桁架和梁常用的单元类型和用途见表9-1。 通过对桁架和梁进行有限元分析,可得到其在各个方向的位移、应力并可得到应力、位

移动画等结果。 第2节桁架的有限元分析实例 一、案例1——2D桁架的有限元分析 图9-1 人字形屋架的示意图 问题 人字形屋架的几何尺寸如图9-1所示。杆件截面尺寸为0.01m2,试进行静力分析,对人字形屋架进行静力分析,给出变形图和各点的位移及轴向力、轴力图。 条件 人字形屋架两端固定,弹性模量为2.0×1011 N/m2,泊松比为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按钮完成清空数据库。 (2)指定新的工作文件名指定工作文件名。选取Utility>Menu> File>Change Jobname,弹出Change Jobname对话框,在Enter New Jobname项输入工作文件名,本例中输入的工作文件名为“2D-spar”,单击OK按钮完成工作文件名的定义。 (3)指定新的标题指定分析标题。选取Utility>Menu>File>Change Title,弹出Change Title对话框,在Enter New Title项输入标题名,本例中输入“2D-spar problem”为标题名,然

基于有限元ANSYS的压力容器应力分析报告

压力容器分析报告

目录 1 设计分析依据 (1) 1.1 设计参数 (1) 1.2 计算及评定条件 (1) 1.3 材料性能参数 (1) 2 结构有限元分析 (2) 2.1 理论基础 (2) 2.2 有限元模型 (2) 2.3 划分网格 (3) 2.4 边界条件 (5) 3 应力分析及评定 (5) 3.1 应力分析 (5) 3.2 应力强度校核 (6) 4 分析结论 (8) 4.1 上封头接头外侧 (9) 4.2 上封头接头内侧 (11) 4.3 上封头壁厚 (13) 4.4 筒体上 (15) 4.5 筒体左 (17) 4.6 下封头接着外侧 (19) 4.7 下封头壁厚 (21)

1 设计分析依据 (1)压力容器安全技术监察规程 (2)JB4732-1995 《钢制压力容器-分析设计标准》-2005确认版 1.1 设计参数 表1 设备基本设计参数 正常设计压力MPa 7.2 最高工作压力MPa 6.3 设计温度℃0~55 工作温度℃5~55 工作介质压缩空气46#汽轮机油 焊接系数φ 1.0 腐蚀裕度mm 2.0 容积㎡ 4.0 容积类别第二类 计算厚度mm 筒体29.36 封头29.03 1.2 计算及评定条件 (1)静强度计算条件 表2 设备载荷参数 设计载荷工况工作载荷工况 设计压力7.2MPa 工作压力6.3MPa 设计温度55℃工作温度5~55℃ 注:在计算包括二次应力强度的组合应力强度时,应选用工作载荷进行计算,本报告中分别选用设计载荷进行计算,故采用设计载荷进行强度分析结果是偏安全的。 1.3 材料性能参数 材料性能参数见表3,其中弹性模型取自JB4732-95表G-5,泊松比根据JB4732-95的公式(5-1)计算得到,设计应力强度分别根据JB4732-95的表6-2、表6-4、表6-6确定。 表3 材料性能参数性能

四杆桁架结构的有限元分析命令流

四杆桁架结构的有限元分析 在ANSYS 平台上,完成相应的力学分析。即如图1所示的结构,各杆的弹性模量和横截面积都为4229.510N/mm E =?, 2 100mm A =,基于ANSYS 平台,求解该结构的节点位移、单元应力以及支反力。 图1四杆桁架结构 完整的命令流 !直接生成有限元模型 / PREP7 !进入前处理 /PLOPTS,DA TE,0 !设置不显示日期和时间 !设置单元、材料,生成节点及单元 ET,1,LINK1 !选择单元类型 UIMP,1,EX, , ,2.95e11, !给出材料的弹性模量 R,1,1e-4, !给出实常数(横截面积) N,1,0,0,0, !生成1号节点,坐标(0,0,0) N,2,0.4,0,0, !生成2号节点,坐标(0.4,0,0) N,3,0.4,0.3,0, !生成3号节点,坐标(0.4,0.3,0) N,4,0,0.3,0, !生成4号节点,坐标(0,0.3,0) E,1,2 !生成1号单元(连接1号节点和2号节点) E,2,3 !生成2号单元(连接2号节点和3号节点) E,1,3 !生成3号单元(连接1号节点和3号节点) E,4,3 !生成4号单元(连接4号节点和3号节点) FINISH !前处理结束 !在求解模块中,施加位移约束、外力,进行求解 /SOLU !进入求解状态(在该状态可以施加约束及外力) ANTYPE,0 !定义分析类型为静力分析 D,1,ALL !将1号节点的位移全部固定 D,2,UY , !将2号节点的Y 方向位移固定 D,4,ALL !将4号节点的位移全部固定 F,2,FX,20000, !在2号节点处施加X 方向的力(20000)

基于有限元的钢板弹簧应力分析

基于有限元的钢板弹簧应力分析 蒋阳 西华大学交通与汽车工程学院 摘 要:本文讨论了利用ANSYS 软件对钢板弹簧进行映射网格划分,并在两簧片的接触区域生成ANSYS 软件所提供的接触单元,建立起多片钢板弹簧的有限元模型。分析了施加预负荷和工作负荷时,板簧应力值显著增长的部位,从而预测板簧产生断裂的部位,可为改进设计提供指导作用。关键词:板簧、仿真、模态 1引言 钢板弹簧具有结构简单,制造、维修方便,除了作为弹性元件外,还可兼起导向和传递侧向、纵向力和力矩的作用,其片间的接触、摩擦在弹簧振动时还将起到阻尼的作用,是重要的高负荷安全部件,目前在商用车上仍被广泛采用[1]。 传统的钢板弹簧设计方法分为:三角形板计算法,板端接触法,共同曲率法[2]。上述三种计算方法对实际工作中的钢板弹簧进行简化,并不能反应实际工作中存在的复杂的非线性状态以及接触情况。 本文利用有限元分析软件ANSYS,对十片钢板弹簧的装配过程和工作过程进行计算分析的基础上, 求得在预负荷和工作负荷作用下的应力与位移等响应情况,为实际钢板弹簧的设计中确定参数提供了依据。 2计算模型的建立 某车型的板簧总成的三维实体结构见图1。板簧建模时,考虑板簧总成对称性,同时为了方便建模,取其1/4为研究对象。通过单元solid45划分网格之后,得到12390个单元和12348个节点。 图1 板簧的三维模型图 3钢板弹簧的材料属性和网格划分 钢板弹簧的材料为60CrMnBa,弹性模量为206GPa,泊松比为0.26。 solid45单元用于构造三维实体结构。单元通过8个节点来定义,每个节点有3个沿着xyz方向平移的自由度。该六面体单元有8个节点,每个节点具有X、Y、Z 三个方向的平动自由度,可以进行塑性分析、蠕变分析、膨胀分析、应力硬化分析、大变形分析和大应变能力[1]。 用单元solid45划分网格之后,得到12390个单元和12348个节点。 4钢板弹簧的片间接触单元的建立 钢板弹簧总成的片与片之间,接触与否事先未知,而且接触后存在着滑移,所以在片与片的节点间建立接触单元,模拟片间的作用力。在ANSYS 中,接触单元是覆盖在分析模型接触面之上的一层单元,主要通过识别接触对与生成接触单元,设置单元关键字和实常数来创建。利用接触单元可以跟踪接触位置、保证接触协调性防止接触表面相互穿透;并在接触表面之间传递接触应力(正压力和摩擦力)。本文选用的接触单元分别是TARGE170 ( 三维目标单元) 与CONTA173(三维8节点面与面接触单元)[2]。本文单元的实常数采用默认值。摩擦采用库仑模型,钢板弹簧之间的摩 101

平面桁架结构的有限元分析

运用ANSYS进行平面刚架模拟建模及误差分析 摘要 有限单元法(或称有限元法)是在当今工程分析中获得最广泛应用的数值计算方法。由于它的通用性和有效性,受到工程技术界的高度重视。伴随着计算机科学和技术的快速发展,现已成为计算机辅助设计和计算机辅助制造的重要组成部分。ANSYS软件是目前世界范围内增长最快的计算机辅助工程(CAE)软件,能与多数计算机辅助设计软件接口,实现数据的共享和交换。本文主要分析平面刚架在均布荷载作用下模拟的有限元模型计算与手工计算之间的误差。 关键字:ANSYS软件有限元平面刚架 PIANE STEEL FRAME WITH ANSYS SIMULATION MODELING AND ERROR ANALYSIS ABSTRACT Finite element method (or finite element method) is the most widely used in modern engineering analysis of numerical calculation method. Because of its versatility and effectiveness, attaches great importance to by the engineering and technology. Along with the rapid development of computer science and technology, has now become a computer aided design and computer aided manufacturing is an important part .At present,the software of ANSY is the fastest growing computer aided engineering (CAE) software on the world, interfacing with the majority of computer aided design software, realizing the sharing and exchange of data. This paper mainly analyzes the model of planar frame software of ANSYS. KEYWARDS:software of ANSYS,finite element,planar frame

基于ANSYS的齿轮应力有限元分析

本科毕业设计 论文题目:基于ansys的齿轮应力有限元分析 学生姓名: 所在院系:机电学院 所学专业:机电技术教育 导师姓名: 完成时间:

摘要 本文主要分析了在ansys中齿轮参数化建模的过程。通过修改参数文件中的齿轮相关参数,利用APDL语言在ANSYS软件中自动建立齿轮的渐开线。再利用图形界面操作模式,通过一系列的镜像、旋转等命令,生成两个相互啮合的大小齿轮。运用有限元分析软件ANSYS对齿轮齿根应力和齿轮接触应力进行分析计算,得出两个大小齿轮的接触应力分布云图。通过与理论分析结果的比较,验证了ANSYS在齿轮计算中的有效性和准确性。 关键词 :ANSYS,APDL,有限元分析,渐开线,接触应力。

Modeling and Finite Element Analysis of Involute Spur Gear Based on ANSYS Abstract We have mainly analyzed spur gear parametrization modelling process in the ansys software. using the APDL language through revises the gear related parameter in the parameter document,we establishes gear's involute automatically in the ANSYS software.Then, using the graphical interface operator schema, through a series of orders ,mirror images, revolving and so on, we produce the big and small gear which two mesh mutually. Carring on the stress analysis of the gear by using the finite element analysis software-- ANSYS, we obtain two big and small gear's contact stress distribution cloud charts. through with the theoretical analysis result's comparison,we explain ANSYS in the gear computation validity and the accuracy. Keywords: ANSYS; APDL;finite element analysis;involute line;contact stress

空间桁架结构程序设计(Fortran)

空间桁架静力分析程序及算例1、变量及数组说明

2、空间桁架结构有限元分析程序源代码 !主程序(读入文件,调用总计算程序,输出结果) CHARACTER IDFUT*20,OUTFUT*20 WRITE(*,*) 'Input Data File name:' READ (*,*)IDFUT OPEN (11,FILE=IDFUT,STATUS='OLD') WRITE(*,*) 'Output File name:' READ (*,*)OUTFUT OPEN(12,FILE=OUTFUT,STATUS='UNKNOWN') WRITE(12,*)'*****************************************' WRITE(12,*)'* Program for Analysis of Space Trusses *' WRITE(12,*)'* School of Civil Engineering CSU *' WRITE(12,*)'* 2012.6.25 Designed By MuZhaoxiang *' WRITE(12,*)'*****************************************' WRITE(12,*)' ' WRITE(12,*)'*****************************************' WRITE(12,*)'*************The Input Data****************' WRITE(12,*)'*****************************************' WRITE(12,100) READ(11,*)NF,NP,NE,NM,NR,NCF,ND WRITE(12,110)NF,NP,NE,NM,NR,NCF,ND 100 FORMAT(6X,'The General Information'/2X,'NF',5X,'NP',5X,'NE',5X,'NM',5X,'NR',& 5X,'NCF',5X,'ND') 110 FORMAT(2X,I2,6I7) NPF=NF*NP NDF=ND*NF CALL ANALYSE(NF,NP,NE,NM,NR,NCF,ND,NPF,NDF) END !******************************************************************** !总计算程序 SUBROUTINE ANALYSE(NF,NP,NE,NM,NR,NCF,ND,NPF,NDF)

圆孔应力有限元分析

圆孔应力有限元分析 陈春山 (安徽工业大学工商学院机械工程系) 摘要:ANSYS软件的应用领域非常广泛,可应用在以下领域:建筑、勘查、地质、水利、交通、电力、测绘、国土、环境、林业、冶金等方面,应用ANSYS软件,对平板中心圆孔的应力集中进行了有限元分析,对圆孔平板在单向和双向应力条件下的应力状况进行了计算和分析,并将有限元结果与解析解进行了比较。 关键词: 平板开小圆孔; 应力集中; 有限元分析 Round hole stress finite element analysis CHEN Chunshan (Industrial & commercial college , anhui university of technology department of mechanical engineering) Abst ract : ANSYS soft ware has a very wide range of applicat ions, can be used in t he following areas: construct ion, exp lorat ion, geology, survey ing an d mapp ing, land, wat er conservancy, t ransport at ion, elect ric p ower, environment, forestry, met allurgy, et c., t he app licat ion of ANSYS software, t he flat round hole at t he centre of the finit e element analysis of st ress concent rat ion of circle hole p lat e under t he condit ion of unidirect ional and bidirect ional st ress calculat ion and analysis, t he stress condit ion and t he finit e element result s are comp ared wit h those of t he analyt ical solut ion Key words: flat open small round hole; Stress concentration; The f inite element analysis l 前言

有限元分析(桁架结构)

有限元上机分析报告 学院:机械工程 专业及班级:机械设计及其自动化08级7班姓名:王浩煜 学号:20082798 题目编号: 2

1.题目概况 1.1 结构组成和基本数据 结构:该结构为一个六根杆组成的桁架结构,其中四根杆组成了直径为800cm的正方形,其他两根杆的两节点为四边形的四个角。 材料:该六根杆截面面积均为100cm2,材料均为Q235,弹性模量为200GPa,对于直径或厚度大于100mm的截面其强度设计值为190Mpa。 载荷:结构的左上和左下角被铰接固定,限制了其在平面内x和y方向的位移,右上角受到大小为2000KN的集中载荷。 结构的整体状况如下图所示: 1.2 分析任务 该分析的任务是对该结构的静强度进行校核分析以验算该结构否满足强度要求。 2.模型建立 2.1 物理模型简化及其分析 由于该结构为桁架结构,故认为每根杆件只会沿着轴线进行拉压,而不会发

生弯曲和扭转等变形。结构中每根杆为铰接连接,有集中载荷作用于最上方的杆和最右方杆的铰接点。 2.2单元选择及其分析 由于该结构的杆可以认为是只受拉压的杆件,故可以使用LINK180单元,该单元是有着广泛工程应用的杆单元,它可以用来模拟桁架、缆索、连杆、弹簧等等。这种三维杆单元是杆轴方向的拉压单元,每个节点具有三个自由度:沿节点坐标系X、Y、Z方向的平动。就像铰接结构一样,不承受弯矩。输入的数据有:两个节点、横截面面积(AREA)、单位长度的质量(ADDMAS)及材料属性。输出有:单元节点位移、节点的应力应变等等。由此可见,LINK180单元适用于该结构的分析。 3.3 模型建立及网格划分 (1)启动Ansys软件,选择Preferences→Structural,即将其他非结构菜单过滤掉。 (2)选择单元类型:选择Preprocessor→Element Type→Add/Edit/Delete→Add,在出现的对话框中选择Link→3d finit stn 180,即LINK180,点击“OK”

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:2 for n2=1:2 KK(DOF(n1), DOF(n2))= KK(DOF(n1), DOF(n2))+k(n1, n2); end

计算结构力学自编平面桁架

平面桁架 程序: #include "stdio.h" #include "math.h" #define unitmax 30 #define pointmax 30 #define matermax 10 #define Pmaxnum 20 #define bindmax 20 main() { int unit[unitmax][4],bind[bindmax][3],number[5]; float point[pointmax][2],material[matermax][2],P[Pmaxnum][3]; float allarray[pointmax*2][pointmax*2]; void readdata(int unit[unitmax][4],float point[pointmax][2],int number[5], int bind[bindmax][3],float material[matermax][2],float P[Pmaxnum][3]); void all(float point[pointmax][2],int unit[unitmax][4],float material[matermax][2], int number[5],float allarray[][pointmax*2]); void equa(float allarray[][pointmax*2],float P[Pmaxnum][3],int bind[][3], int number[]); void result(int unit[unitmax][4],float point[pointmax][2],int number[5], float material[matermax][2],float allarray[][pointmax*2]); readdata(unit,point,number,bind,material,P); all(point,unit,material,number,allarray); equa(allarray,P,bind,number); result(unit,point,number,material,allarray); } /******************************************************************/ void readdata(int unit[unitmax][4],float point[pointmax][2],int number[5], int bind[bindmax][3],float material[matermax][2],float P[Pmaxnum][3])

桁架结构的有限元分析MATLAB

力09创新实践 桁架结构有限元分析 学号 20092715 班级力0901-2 姓名魏强 指导教师房学谦 完成日期 2012/6/26 桁架结构有限元分析 摘要

从系统物理概念和力学原理推导有限元计算格式的方法叫做直接刚度法。本文利用推导出得有限元计算格式,通过MATLAB软件进行矩阵运算,对5杆桁架结构进行了内力分析。利用对比的方法,对照多组荷载,分析其受力的情况,为实际问题提供参考。 关键词:有限元法、MATLAB、桁架结构、内力分析 一、引言 1.工程背景及重要性 桁架结构(Truss structure)中的桁架指的是桁架梁,是格构化的一种梁式结构。桁架结构常用于大跨度的厂房、展览馆、体育馆和桥梁等公共建筑中。由于大多用于建筑的屋盖结构,桁架通常也被称作屋架。 各杆件受力均以单向拉、压为主,通过对上下弦杆和腹杆的合理布置,可适应结构内部的弯矩和剪力分布。由于水平方向的拉、压内力实现了自身平衡,整个结构不对支座产生水平推力。结构布置灵活,应用范围非常广。桁架梁和实腹梁(即我们一般所见的梁)相比,在抗弯方面,由于将受拉与受压的截面集中布置在上下两端,增大了内力臂,使得以同样的材料用量,实现了更大的抗弯强度。在抗剪方面,通过合理布置腹杆,能够将剪力逐步传递给支座。这样无论是抗弯还是抗剪,桁架结构都能够使材料强度得到充分发挥,从而适用于各种跨度的建筑屋盖结构。更重要的意义还在于,它将横弯作用下的实腹梁内部复杂的应力状态转化为桁架杆件内简单的拉压应力状态,使我们能够直观地了解力的分布和传递,便于结构的变化和组合。 在建筑结构中,桁架结构是一种应用比较普遍的结构形式,在桥梁工程、大型建筑、船舶工程、港口机械等工程领域均有广泛应用。在我国桁架结构发展迅速且应用最为广泛,如屋架、网架结构等。为了增加建筑的表现力,近些年来管桁架结构得到了许多业主的青睐,在大量的屋面结构中采用。 2.目前问题的研究现状 目前在普遍刚桁架的结构设计中,工程中普遍采用的发放时按理想铰接模型进行计算,并很据计算出的杆件界面应力选择合适的杆件型号。计算桁架结构内力时,一般采用如下基本假定:(1)接单均为铰接;(2)杆件轴线平直相交于节点中心;(3)荷载作用线通过桁架的节点。对于平面桁架还要求所有杆件轴线及荷载作用线在同一平面内。 对于桁架结构的应力分析,在方法上,结构力学中有结点法和截面法,另外

ansys桁架和梁的有限元分析

桁架和梁的有限元分析 第一节基本知识 一、桁架和粱的有限元分析概要 1.桁架杆系的有限元分析概要 桁架杆系系统的有限元分析问题是工程中晕常见的结构形式之一,常用在建筑的屋顶、机械的机架及各类空间网架结构等多种场合。 桁架结构的特点是,所有杆件仅承受轴向力,所有载荷集中作用于节点上。由于桁架结构具有自然离散的特点,因此可以将其每一根杆件视为一个单元,各杆件之间的交点视为一个节点。 2.梁的有限元分析概要 梁的有限元分析问题也是是工程中最常见的结构形式之一,常用在建筑、机械、汽车、工程机械、冶金等多种场合。 梁结构的特点是,梁的横截面均一致,可承受轴向、切向、弯矩等载荷。根据梁的特点,等截面的梁在进行有限元分析时,需要定义梁的截面形状和尺寸,用创建的直线代替梁,在划分网格结束后,可以显示其实际形状。 二、桁架和梁的常用单元 桁架和梁常用的单元类型和用途见表7-1。 通过对桁架和粱进行有限元分析,可得到其在各个方向的位移、应力并可得到应力、位移动画等结果。 第128页

第二节桁架的有限元分析实例案例1--2D桁架的有限元分析 问题 人字形屋架的几何尺寸如图7—1所示。杆件截面尺寸为0.01m^2,试进行静力分析,对人字形屋架进行静力分析,给出变形图和各点的位移及轴向力、轴力图。 条件 人字形屋架两端固定,弹性模量为2.0x10^11N/m^2,泊松比为0.3。 解题过程 制定分析方案。材料为弹性材料,结构静力分析,属21)桁架的静力分析问题,选用Link1单元。建立坐标系及各节点定义如图7-1所示,边界条件为1点和5点固定,6、7、8点各受1000N的力作用。 1.ANSYS分析开始准备工作 (1)清空数据库并开始一个新的分析选取Utility Menu>File>Clear&Start New,弹出Clears database and Start New对话框,单击OK按钮,弹出Verify对话框,单击OK按钮完成清空数据库。 (2)指定新的工作文件名指定工作文件名。选取Utility Menu>File>Change Jobname,弹出Change Jobname对话框,在Enter New Jobname项输入工作文件名,本例中输入的工作文件名为“2D-spar”,单击OK按钮完成工作文件名的定义。 (3)指定新的标题指定分析标题。选取Ufility Menu>File>Change Title,弹出ChangeTitle对话框,在Enter New Tifie项输入标题名,本例中输入“2D-spar problem'’为标题名,然后单击OK按钮完成分析标题的定义。 (4)重新刷新图形窗9 选取Utility Menu>Plot>Replot,定义的信息显示在图形窗口中。 (5)定义结构分析运行主菜单Main Menu>Preferences,出现偏好设置对话框,赋值分析模块为Structure结构分析,单击OK按钮完成分析类型的定义。 2.定义单元类型 运行主菜单Main Menu>Preprocessor>Element Type>Add/Edit/Delete命令,弹出Element Types对话框,单击Add按钮新建单元类型,弹出Library of Element Types对话框,先选择

平面桁架

常见问题解答七 第七章平面桁架 1.屋盖结构体系分为哪两大类?各有何适用? 答:屋盖结构体系分为无檩屋盖与有檩屋盖。 (1)无檩屋盖:无檩屋盖一般用于预应力混凝土大型屋面板等重型屋面,将屋面板直接放在屋架或天窗架上。 预应力混凝土大型屋面板的跨度通常采用6m,有条件时也可采用12m。当柱距大于所采用的屋面板跨度时,可采用托架支承中间屋架。 (2)有檩屋盖:有檩屋盖常用于轻型屋面材料的情况。如压型钢板、压型铝合金板、石棉瓦、瓦楞铁皮等。石棉瓦和瓦楞铁皮屋面,屋架间距通常为6m;当柱距大于或等于12m时,则用托架支承中间屋架。钢板和压型铝合金板屋面,屋架间距常大于或等于12m。 2.常用的屋架形式有哪几种?各有何特点和适用? 答:屋架外形常用的有三角形、梯形、平行弦和人字形等。 三角形桁架适用于陡坡屋面(i>1/3)的有檩屋盖体系,通常与柱子铰接,房屋的整体横向刚度较低,简支屋架弦杆支座处内力大,跨中内力小,截面不能充分发挥作用。支座处上下弦杆交角过小内力又较大,使支座节点构造复杂。 三角形屋架的腹杆布置有芬克式和人字式。芬克式的腹杆较多,但压杆短、拉杆长,受力相对合理,且可分为两个小桁架制作与运输,较为方便。人字式腹杆的节点较少,但受压腹杆较长,适用于跨度较小(L≤18m)的情况。人字式屋架的抗震性能好,所以在强地震烈度地区,,常用人字式腹杆的屋架。单斜式腹杆的屋架,腹杆和节点数目较多,只适用于下弦需要设置天棚的屋架,一般情况较少采用。

梯形屋架适用于屋面坡度较为平缓的无檩屋盖体系,与简支受弯构件的弯矩图形比较接近,弦杆受力较为均匀。屋架与柱的连接可以铰接也可以刚接。刚接可提高建筑物的横向刚度。腹杆体系可采用单斜式、人字式和再分式。屋架支承分为下承式和上承式。与柱刚接的屋架宜采用下承式。下承式使柱计算高度减小又便于在下弦设置屋盖纵向水平支撑,以往多采用之,但上承式使屋架重心降低,支座斜腹杆受拉,且安装方便,近年来逐渐推广使用。 图7-2 梯形屋架 人字形屋架的上、下弦可以是平行的,坡度为1/20~1/10,节点构造较为统一;也可以上、下弦具有不同坡度或者下弦有一部分水平段。多用于较大跨度。宜采用上承式。 人字形和梯形屋架的中部高度主要取决于经济要求,一般为(1/10~1/8)L,与柱刚接的梯形屋架,端部高度一般为(1/16~1/12)L,通常取为2.0~2.5m。与柱铰接的梯形屋架,端部高度可按跨中经济高度和上弦坡度来决定。人字形屋架因中高度一般为2.0~2.5m,跨度大于36m时可取较大高度但不宜超过3m;端部高度一般为跨度的1/18~1/12,人字形屋架可适应不同的屋面坡度,但与柱刚接时,屋架轴线坡度大于1/7,就应视为折线横梁进行框架分析;与柱铰接时,即使采用了上承式也应考虑竖向荷载作用下折线拱的推力对柱的不利影响,设计时要求在屋面板及檩条等安装完毕后再将屋架支座焊接固定。 图7-3 人字形屋架和平行弦屋架 平行弦桁架在构造方面有突出的优点,弦杆及腹杆分别等长、节点形式相同、能保证桁架的杆件重复率最大,且可使节点构造形式统一,便于制作工业化。 3.屋盖支撑有何作用? 答:屋盖支撑的作用: ①保证结构的空间整体作用。 ②避免压杆侧向失稳,防止拉杆产生过大的振动。 ③承担和传递水平荷载。 ④保证结构安装时的稳定与方便

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