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

运用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

1 有限元法简介

有限元法是一种高效能、常用的数值计算方法。科学计算领域，常常需要求解各类微分方程，而许多微分方程的解析解一般很难得到，使用有限元法将微分方程离散化后，可以编制程序，使用计算机辅助求解。有限元法在早期是以变分原理为基础发展起来的，所以它广泛地应用于以拉普拉斯方程和泊松方程所描述的各类物理场中（这类场与泛函的极值问题有着紧密的联系）。自从1969年以来，某些学者在流体力学中应用加权余数法中的迦辽金法(Galerkin)或最小二乘法等同样获得了有限元方程，因而有限元法可应用于以任何微分方程所描述的各类物理场中，而不再要求这类物理场和泛函的极值问题有所联系。

2 ANSYS 软件简介

ANSYS 软件是美国ANSYS 公司研制的大型通用有限元分析软件，是世界范围内增长最快的计算机辅助工程(CAE)软件，是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件。在核工业、铁道、石油化工、航空航天、机械制造、能源、汽车交通、国防军工、电子、土木工程、造船、生物医学、轻工、地矿、水利、日用家电等领域有着广泛的应用。ANSYS 功能强大，操作简单方便，现在已成为国际最流行的有限元分析软件。目前，中国100多所理工院校均采用ANSYS 软件进行有限元分析或者作为标准教学软件。

3 实例简介

3.1问题描述

如下图所示的平面刚架，两杆均为45a 号工字钢，A=102错误！未找到引用源。,截面二次矩I=32240错误！未找到引用源。,用有限元方法分析该结构的变形，受力节点的位移，及各杆的轴力、剪力和弯矩，并将所得内力图与手算结果进行对比。结构中各个截面的参数都为：3.0,101.211=?=u Pa E 。

4 ANSYS结构建模分析

4.1定义参数

4.1.1指定工程名和分析标题

（1）启动ANSYS软件，选择File→Change Jobname命令，弹出如图所示的对话框：

（2）在对话框中输入”dam”,同时把【New log and error files】中的复选框选为yes,单击确定。

4.1.2定义单元类型

（1）选择Main Menu→Preprocessor→Element Type→Add/Edit/Delete 命令，弹出如图所示【Element Types】对话框：

（2）单击【Element Types】对话框中的【Add】按钮，在弹出的如下所示【Library of Element Types】对话框：

（3）选择左边文本框中的【Beam】选项，右边文本框中的[2D elastic 3]选项，单击确定。

（4）返回【Element Types】对话框。如下图所示：

(5)单击【close】按钮结束，单元定义完成。

4.1.3定义单元常数

（1）在ANSYS程序主界面中选择Main Menu→Preprocessor→Real Constants→Daa/Edit/Delete命令，弹出如下所示【Real Constants】对话框：

（2）单击【Add】按钮，进行下一个【Choose Element Types】对话框，选择【Beam 3】单元，单击确定，弹出如下对话框：

（3）定义好参数，单击ok,回到[Real Constants]对话框，单击Close。

4.1.4定义材料参数

（1）在ANSYS程序主界面，选择Main Menu→Preprocessor→Material Pros →Material Models命令，弹出如下所示【Define Material Behavior】对话框：

（2）选择对话框右侧Structural→Linear→Elastic→Isotropic命令。并单击【Isotropic】选项，接着弹出如下所示【Linear Isotropic Proprties for Material Number】对话框：

（3）输入好相应的数字，单击ok，回到【Define Material Behavior】对话框，直接关闭对话框，至此材料参数定义完毕。

4.2创建几何模型

（1）在ANSYS程序主界面，选择Main Menu→Preprocessor→Modeling→

Creat→Keypoints→In Active CS命令，弹出如下所示对话框：

（2）创建相应的关键点，单击OK按钮，得到相应的几何模型。

4.3划分网格

(1)选择Main Menu→Preprocessor→Meshing→Mesh Tool命令，弹出实体选择对话框:

(2)在Mesh Tool→Size Controls→Lines→Set命令，弹出实体对话框点击Pick All→OK，在弹出实体对话框如下图所示，在NDIV中填写20。

（3）在Mesh Tool→Mesh命令，弹出实体对话框Mesh Lines→Pick All，网格已划分完毕。

4.4加载数据

（1）选择Main Menu→Preprocessor→Solution→Define Loads→Apply→Structural→Displacement→On Keypoints命令，弹出实体选择对话框，在关键点处加载相关约束如下图所示。

（2）选择Main Menu→Preprocessor→Solution→Define Loads→Apply→Structural→Pressure→On Nodes命令，弹出实体选择对话框，在关键点处加相应荷载如下图所示。

4.5求解

（1）选择Main Menu→Solution→Solve→Current LS命令，弹出如下图所示窗口，其中【/STATUS Command】窗口里包括了所要计算模型的求解信息和荷载信息。

（2）单击上图OK按钮，程序进行计算。

（3）计算完毕，出现如图所示对话框：

4.6结果分析

4.6.1显示变形图

（1）选择Main Menu→General Postproc→Read Results→First Set命令，读入最初结果文件。

(2)选择Main Menu→General Postproc→Plot Results→Deformed Shape 命令，弹出【Plot Deformed Shape】对话框，如下所示：

（3）选择【Def+undeformed】选项，并单击OK，出现如下最终变形图。

4.6.2显示内力图（轴力图、剪力图、弯矩图）

选择Element Table→Define Table→Add命令，弹出如下对话框：

然后在SMISC,后面输入1，点Apply,然后重复输入7，2，8，12，完成之后选择命令Plot Results→Contour Plot→Line Elem Res,会弹出如下对话框：

然后在LabI后面选1，LabJ后面选7，点OK,得到的是轴力图，选2，8得到的是剪力图，选6，12得到的是弯矩图，分别如下图所示：

4.7手算结果与ansys分析结果对比分析

4.7.1手算所得的内力图如下图所示：

4.7.2对比分析

由ANSYS及人工手算所得结果进行对比，可以看出ANSYS软件计算结果的轴力图和剪力图跟手算结果是一样的，弯矩图中最大弯矩值也是相等的。某些点处弯矩值可能不一样，但误差很小，造成误差的原因可能是划分单元的等分不够。

5 结语

经过近50年特别是近30年的发展，有限元法的基础理论和方法已经比较成熟，已成为当今工程技术领域中应用最为广泛，成效最为显著的数值分析方法。ANSYS软件可直接输入结构几何图形并在图形里面进行有限元网格划分，还可

直接输入约束条件及荷载数据，这样使繁琐的数据文件直接填写工作变得直观、轻松，又容易发现输入过程中的错误，及时修正，可以应用于结构分析的各个领域及进行各个场之间的耦合，该软件是目前国际上最流行、最实用、最全面的大型通用有限元分析程序，对其深入研究必将对我国的经济建设起促进作用。

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对话框，先选择

结构力学实验-平面桁架结构的设计

结构力学实验土木建筑学院 实验名称：平面桁架结构的设计 实验题号:梯形桁架D2-76 姓名: 学号： 指导老师: 实验日期:

一、实验目的 在给定桁架形式、控制尺寸和荷载条件下，对桁架进行内力计算，优选杆件截面，并进行刚度验算。 ①掌握建立桁架结构力学模型的方法，了解静定结构设计的基本过程； ②掌握通过多次内力和应力计算进行构件优化设计的方法； ③掌握结构刚度验算的方法。 梯形桁架D ；其中结点1到结点7的水平距离为15m；结点1到结点8的距离为2m；结点7到结点14的距离为3m。选用的是Q235钢，[ɑ]=215MPa。

完成结构设计后按如下步骤计算、校核、选取、设计、优化 二、强度计算 1）轴力和应力 2）建立结构计算模型后，由“求解→内力计算”得出结构各杆件的轴力N（见图3）再由6=N/A得出各杆件应力。 表1内力计算 杆端内力值 ( 乘子 = 1) -------------------------------------------------------------------------------------------- 杆端 1 杆端 2 ------------------------------------- ------------------------------------------ 单元码轴力剪力弯矩轴力剪力弯矩 -------------------------------------------------------------------------------------------- 1 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 2 51.9230769 0.00000000 0.00000000 51.9230769 0.00000000 0.00000000 3 77.1428571 0.00000000 0.00000000 77.1428571 0.00000000 0.00000000 4 67.5000000 0.00000000 0.00000000 67.5000000 0.00000000 0.00000000 5 39.7058823 0.00000000 0.00000000 39.7058823 0.00000000 0.00000000 6 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 7 -54.0000000 0.00000000 0.00000000 -54.0000000 0.00000000 0.00000000 8 -52.0383336 0.00000000 0.00000000 -52.0383336 0.00000000 0.00000000 9 -77.3140956 0.00000000 0.00000000 -77.3140956 0.00000000 0.00000000 10 -81.1798004 0.00000000 0.00000000 -81.1798004 0.00000000 0.00000000 11 -81.1798004 0.00000000 0.00000000 -81.1798004 0.00000000 0.00000000 12 -67.6498337 0.00000000 0.00000000 -67.6498337 0.00000000 0.00000000 13 -39.7940198 0.00000000 0.00000000 -39.7940198 0.00000000 0.00000000 14 -54.0000000 0.00000000 0.00000000 -54.0000000 0.00000000 0.00000000 15 66.4939824 0.00000000 0.00000000 66.4939824 0.00000000 0.00000000 16 -41.5384615 0.00000000 0.00000000 -41.5384615 0.00000000 0.00000000 17 33.3732229 0.00000000 0.00000000 33.3732229 0.00000000 0.00000000 18 -21.8571428 0.00000000 0.00000000 -21.8571428 0.00000000 0.00000000 19 5.27613031 0.00000000 0.00000000 5.27613031 0.00000000 0.00000000 20 -18.0000000 0.00000000 0.00000000 -18.0000000 0.00000000 0.00000000 21 19.7385409 0.00000000 0.00000000 19.7385409 0.00000000 0.00000000 22 -31.5000000 0.00000000 0.00000000 -31.5000000 0.00000000 0.00000000 23 42.0090820 0.00000000 0.00000000 42.0090820 0.00000000 0.00000000 24 -47.6470588 0.00000000 0.00000000 -47.6470588 0.00000000 0.00000000 25 62.0225709 0.00000000 0.00000000 62.0225709 0.00000000 0.00000000

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

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

第七专题平面桁架结构

平面桁架结构 一、平面桁架的形式 1．屋盖结构体系 屋盖分为无檩屋盖有檩屋盖。无檩屋盖一般用于预应力混凝土大型屋面板等重型屋面，将屋面板直接放在屋架上。有檩屋盖常用于轻型屋面材料的情况。 2．屋架的形式 屋架外形常用的有三角形、梯形、平行弦和人字形等。 桁架外形应尽可能与其弯矩图接近，这样弦杆受力均匀，腹杆受力较小。腹杆的布置应尽量用长杆受拉、短杆受压，腹杆的数目宜少，总长度要短，斜腹杆的倾角一般在30°～60°之间，腹杆布置时应注意使荷载都作用在桁架的节点上。 （1）三角形桁架 三角形桁架适用于陡坡屋面（i＞1/3）的有檩屋盖体系，屋架通常与柱子只能铰接。弯矩图与三角形的外形相差悬殊，弦杆受力不均，支座处内力较大，跨中内力较小，弦杆的截面不能充分发挥作用。支座处上、下弦杆交角过小内力又较大，使支座节点构造复杂。 （2）梯形桁架 梯形屋架适用于屋面坡度较为平缓的无檩屋盖体系，它与简支受弯构件的弯矩图形比较接近，弦杆受力较为均匀。梯形屋架与柱的连接可以做成铰接也可以做成刚接。梯形屋架的中部高度一般为（1/10～1/8）L，与柱刚接的梯形屋架，端部高度一般为（1/16～1/12）L，通常取为2.0～2.5m。与柱铰接的梯形屋架，端部高度可按跨中经济高度和上弦坡度决定。 （3）人字形桁架 人字形屋架的上、下弦可以是平行的，坡度为1/20～1/10，节点构造较为统一；也可以上、下弦具有不同坡度或者下弦有一部分水平段，以改善屋架受力情况。人字形屋架因中高度一般为2.0～2.5m，跨度大于36m时可取较大高度但不宜超过3m；端部高度一般为跨度的1/18～1/12。 （4）平行弦桁架 平行弦桁架在构造方面有突出的优点，弦杆及腹杆分别等长、节点形式相同、能保证桁架的杆件重复率最大，且可使节点构造形式统一，便于制作工业化。 3．托架形式 支承中间屋架的桁架称为托架，托架一般采用平行弦桁架，其腹杆采用带竖杆的人字形体系。托架高度般取跨度的1/5～1/10，托架的节间长度一般为2m或3m。 二、屋盖支撑

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

四杆桁架结构的有限元分析 在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)

基于MATLAB的平面刚架静力分析

基于MATLAB 的平面刚架静力分析 为了进一步理解有限元方法计算的过程，本文根据矩阵位移法的基本原理应用MATLAB 编制计算程序对以平面刚架结构进行了静力分析。本文还利用ANSYS 大型商用有限元分析软件对矩阵位移法的计算结果进行校核，发现两者计算结果相当吻合，验证了计算结果的可靠性。 一、 问题描述 如图1所示的平面刚架，各杆件的材料及截面均相同，E=210GPa ，截面为0.12×0.2m 的实心矩形，现要求解荷载作用下刚架的位移和内力。 5m 4m 3m 图1 二、矩阵位移法计算程序编制 为编制程序方便考虑，本文计算中采用“先处理法”。具体的计算步骤如下。

（1） 对结构进行离散化，对结点和单元进行编号，建立结构（整体）坐标系 和单元（局部）坐标系，并对结点位移进行编号； （2） 对结点位移分量进行编码，形成单元定位向量e λ； （3） 建立按结构整体编码顺序排列的结点位移列向量δ，计算固端力e F P 、等 效结点荷载E P 及综合结点荷载列向量P ； （4） 计算个单元局部坐标系的刚度矩阵，通过坐标变换矩阵T 形成整体坐标 系下的单元刚度矩阵e T e K T K T = ； （5） 利用单元定位向量形成结构刚度矩阵K ； （6） 按式1=K P δ- 求解未知结点位移； （7） 计算各单元的杆端力e F 。 根据上述步骤编制了平面刚架的分析程序。程序中单元刚度矩阵按下式计算。 32322 23 2 32 22 0000 1261260 064620 00001261260062640 EA EA l l EI EI EI EI l l l l EI EI EI EI l l l l K EA EA l l EI EI EI EI l l l l EI EI EI EI l l l l ??- ??? ???- ?? ? ???- ??? ?=??-?? ? ???---??? ???-??? ?

桁架结构

桁架结构 桁架结构（Truss structure）中的桁架指的是桁架梁，是格构化的一种梁式结构。桁架结构常用于大跨度的厂房、展览馆、体育馆和桥梁等公共建筑中。由于大多用于建筑的屋盖结构，桁架通常也被称作屋架。 主要结构特点 各杆件受力均以单向拉、压为主，通过对上下弦杆和腹杆的合理布置，可适应结构内部的弯矩和剪力分布。由于水平方向的拉、压内力实现了自身平衡，整个结构不对支座产生水平推力。结构布置灵活，应用范围非常广。桁架梁和实腹梁（即我们一般所见的梁）相比，在抗弯方面，由于将受拉与受压的截面集中布置在上下两端，增大了内力臂，使得以同样的材料用量，实现了更大的抗弯强度。在抗剪方面，通过合理布置腹杆，能够将剪力逐步传递给支座。这样无论是抗弯还是抗剪，桁架结构都能够使材料强度得到充分发挥，从而适用于各种跨度的建筑屋盖结构。更重要的意义还在于，它将横弯作用下的实腹梁内部复杂的应力状态转化为桁架杆件内简单的拉压应力状态，使我们能够直观地了解力的分布和传递，便于结构的变化和组合。 桁架的历史演变 只受结点荷载作用的等直杆的理想铰结体系称桁架结构。它是由一些杆轴交于一点的工程结构抽象简化而成的。桁架在建造木桥和屋架上最先见诸实用。古罗马人用桁架修建横跨多瑙河的特雷江桥的上部结构（发现于罗马的浮雕中，文艺复兴时期，意大利建筑师（拔拉雕 Palladio）也开始采用木桁架建桥出现朗式、汤式、豪式桁架。英国最早的金属桁架是在1845年建成的，适合汤式木桁架相似的格构桁架，第二年又采用了三角形的华伦式桁架。 桁架种类 桁架可按不同的特征进行分类。 根据桁架的外形分为：平行弦桁架（便于布置双层结构；利于标准化生产，但杆力分布不够均匀）、折弦桁架（如抛物线形桁架梁，外形同均布荷载下简支梁的弯矩图，杆力分布均匀，材料使用经济，构造较复杂）、三角形桁架（杆力分布更不均匀，构造布置困难，但斜面符合屋顶排水需要）。

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

运用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

2016基本平面刚架各种荷载MATLAB程序

% 平面刚架MATLAB程序 % 2003.9.16 2007.2.28 2008.4.1 2009.10 2011.10 2013.9 2014.09 2016.03 %************************************************* % 变量说明 % NPOIN NELEM NVFIX NFPOIN NFPRES % 总结点数,单元数, 约束个数, 受力结点数, 非结点力数 % COORD LNODS YOUNG % 结构节点坐标数组, 单元定义数组, 弹性模量 % FPOIN FPRES FORCE FIXED % 结点力数组，非结点力数组，总体荷载向量, 约束信息数组 % HK DISP % 总体刚度矩阵，结点位移向量 %************************************************** format short e %设定输出类型 clear %清除内存变量 FP1=fopen('6-6.txt','rt') %打开初始数据文件 %读入控制数据 NELEM=fscanf(FP1,'%d',1); %单元数 NPOIN=fscanf(FP1,'%d',1); %结点数 NVFIX=fscanf(FP1,'%d',1); %约束数 NFPOIN=fscanf(FP1,'%d',1); %作用荷载的结点个数 NFPRES=fscanf(FP1,'%d',1); %非结点荷载数 YOUNG=fscanf(FP1,'%f',1); %弹性模量 % 读取结构信息 LNODS=fscanf(FP1,'%f',[6,NELEM])' % 单元定义：左、右结点号，面积，惯性矩，线膨胀系数，截面高度（共计NELEM组）COORD=fscanf(FP1,'%f',[2,NPOIN])' % 坐标：x,y坐标（共计NPOIN 组） FPOIN=fscanf(FP1,'%f',[4,NFPOIN])' % 节点力（共计NFPOIN 组）：受力结点号、X方向力（向右正）， % Y方向力（向上正），M力偶（逆时针正） FPRES=fscanf(FP1,'%f',[7,NFPRES])' % 均布力（共计 % NFPRES 组）：单元号、荷载类型、荷载大小、距离左端长度，温差=（下端-上端）梯形上边。下边(改) % 荷载类型1-均布荷载2-横向集中力3-纵向集中力4-三角形荷载5-温度荷载6-梯形荷载 FIXED=fscan f(FP1,'%f',NVFIX)' % 约束信息：约束对应的位移编码（共计NVFIX 组） %--------------------------------------------------------- HK=zeros(3*NPOIN,3*NPOIN); % 张成总刚矩阵并清零 FORCE=zeros(3*NPOIN,1); % 张成总荷载向量并清零 %形成总刚 for i=1:NELEM % 对单元个数循环

平面桁架结构matlab

桁架结构计算第四章P56 ******************************************************************************* function y=plane_truss_element_stiffness(E,A,L,theta) %平面桁架单元刚度 x=theta*pi/180; C=cos(x); S=sin(x); y=E*A/L*[ C*C C*S -C*C -C*S; C*S S*S -C*S -S*S; -C*C -C*S C*C C*S; -C*S -S*S C*S S*S];%平面桁架刚度矩阵 ******************************************************************************* function y=plane_truss_assemble(K,k,i,j) %平面桁架组装 K(2*i-1,2*i-1)=K(2*i-1,2*i-1)+k(1,1); K(2*i-1,2*i)=K(2*i-1,2*i)+k(1,2); K(2*i-1,2*j-1)=K(2*i-1,2*j-1)+k(1,3); K(2*i-1,2*j)=K(2*i-1,2*j)+k(1,4); K(2*i,2*i-1)=K(2*i,2*i-1)+k(2,1); K(2*i,2*i)=K(2*i,2*i)+k(2,2); K(2*i,2*j-1)=K(2*i,2*j-1)+k(2,3); K(2*i,2*j)=K(2*i,2*j)+k(2,4); K(2*j-1,2*i-1)=K(2*j-1,2*i-1)+k(3,1); K(2*j-1,2*i)=K(2*j-1,2*i)+k(3,2); K(2*j-1,2*j-1)=K(2*j-1,2*j-1)+k(3,3); K(2*j-1,2*j)=K(2*j-1,2*j)+k(3,4); K(2*j,2*i-1)=K(2*j,2*i-1)+k(4,1); K(2*j,2*i)=K(2*j,2*i)+k(4,2); K(2*j,2*j-1)=K(2*j,2*j-1)+k(4,3); K(2*j,2*j)=K(2*j,2*j)+k(4,4); y=K; ******************************************************************************* function y=plane_truss_element_force(E,A,L,theta,u)%力的表达式 x=theta*pi/180; C=cos(x); S=sin(x); y=E*A/L*[-C -S C S]*u; ******************************************************************************* function y=plane_truss_element_stress(E,L,theta,u) %应力表达式 x=theta*pi/180; C=cos(x); S=sin(x); y=E/L*[-C -S C S]*u; ***************************************************************************************************** *****************************************************************************************************

桁架单元例子MATLAB 1

no axial forces acting on the beam. Use two elements to solve the problem. (a) Determine the deflection and slope at x = 0.5, 1 and 1.5 m; (b) Draw the bending moment and shear force diagrams for the entire beam; (c) What are the support reactions? (d) Use the beam element shape functions to plot the deflected shape of the beam. Use EI = 1,000 Nm, L = 1 m, and F = 1,000 N. Solution: Solution: (a) Given, ?(?)=?(?)=?=1?; ??=1000??; ?=1000? For any element of length L, the structural stiffness matrix is defined as, ???=????? 126? ?126? 6?4?? ?6?2?? ?12?6? 12?6?6?2?? ?6?4?? ? The element stiffness matrix for element 1 is: ?(?) =????(?)???126?12664?62?12?612?6 62?64?=1000?126?126 64?62?12?612?662?64 ? The element stiffness matrix for element 2 is: Element 1 Element 2

有限元分析(桁架结构)

有限元上机分析报告 学院：机械工程 专业及班级：机械设计及其自动化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”

基于MATLAB的桁架结构优化设计

基于MAT LAB 的桁架结构优化设计 林　琳　张云波 (华侨大学土木系福建泉州　362011) 【摘　要】　介绍了基于BP 神经网络的全局性结构近似分析方法,解决了结构优化设计问题中变量的非线性映射问题。在此基础上,利用改进的遗传算法,对桁架结构在满足应力约束条件下进行结构最轻优化设计。利用 Matlab 的神经网络工具箱,编程求解了三杆桁架优化问题。 【关键词】　改进遗传算法;BP 神经网络;结构优化设计;满应力准则 【中图分类号】　T U20114 【文献标识码】　A 【文章编号】　100126864(2003)01-0034-03 TRUSS STRUCTURA L OPTIMIZATON BASE D ON MAT LAB LI N Lin ZH ANG Y unbo (Dept.of Civil Engineering ,Huaqiao University ,Quanzhou ,362011) Abstract :Optimal structural design method based on BP neural netw ork and m odified genetic alg orithm were proposed in this paper.The high parallelism and non -linear mapping of BP neural netw ork ,an approach to the global structural approximation analysis was introduced.It can s olve the mapping of design variables in structural optimization problems.C ombining with an im proved genetic alg orithm ,the truss structure is optimized to satis fy the full stress criteria.Under the condition of MAT LAB 5.3,an exam ple of truss structure has been s olved by this method. K ey w ords :G enetic alg orithm ;BP neural netw ork ;Structural optimization design ;Full stress principle 结构优化设计,就是在满足结构的使用和安全要求的基础上,降低工程造价,更好地发挥投资效益。传统的优化方法有工程法和数学规划法,其难以解决离散变量问题,对多峰问题容易陷入局部最优,且对目标函数要求有较好的连续性或可微性。而近年来提出的基于生物自然选择与遗传机理的随机搜索遗传算法对所解的优化问题没有太多的数学要求,可以处理任意形式的目标函数和约束,对离散设计变量的优化问题尤为有效。进化算子的各态历经性使得遗传算法能够非常有效地进行概率意义下的全局搜索,能高效地寻找到全局最优点。但采用遗传算法时,进化的每一代种群成员必须要进行结构分析,因此所需的结构分析次数较多。 1　桁架结构优化设计问题的表述 在满足应力约束条件下的桁架重量最轻优化问题为: min w (A )=Σn i =1ρA i L i s.t 1 σi ≤[σi ] (i =1,2……n ) A min ≤A i ≤A max w (A )为结构总重量,ρ为材料密度,L i 为第i 杆的长度,A i 为第i 杆件面积,σi 为第i 杆的应力,[σi ]为第i 杆的许用 应力,A min 、A max 分别为杆件面积的下界与上界;n 为杆件总数。 2　神经网络结构近似分析方法 人工神经网络是由大量模拟生物神经元功能的简单处理单元相互连接而成的巨型复杂网络,它是一个具有高度非线 性的超大规模连续时间自适应信息处理系统,易处理复杂的非线性建模问题。文献[1]在K olm og orov 多层神经网络映射存在定理的基础上,针对近似结构分析问题提出的多层神经网络映射存在定理,确定了近似结构分析的神经网络的基本模型。从理论上证明一个三层神经网络可用来描述任一弹性结构的应力、位移等变量和结构设计变量之间的映射关系,为利用人工神经网络来进行结构近似分析提供理论基础。 211　BP 神经网络及其算法改进 BP 神经网络,即误差反向传播神经网络。其最主要的 特性就是具有非线性映射功能。1989年R obert Hecht -Niel 2 s on 证明了对于任何闭区间内的一个连续函数,都可用一个 隐含层的BP 网络来逼近。因而一个三层BP 网络可完成任意的n 维到m 维的映照,它由输入层、隐层和输出层构成。 传统的BP 网络存在着局部极小问题和收敛速度较慢的问题,因此本文采用了动量法和学习率自适应调整的策略,提高了学习速度并增加了算法的可靠性。 动量法考虑了以前时刻的梯度方向,降低了网络对误差曲面局部细节的敏感性,有效地抑制了网络陷于局部极小。 w (k +1)=w (k )+α[(1-η)D (k )+ηD (k -1)] α(k )=2λα(k -1)λ=stg n[D (k )D (k -1)] w (A )为权值向量,D (k )=- 5E 5w (k ) 为k 时刻的负梯度,D (k -1)为k -1时刻的负梯度,η为动量因子,α为学习率。 4 3 低　温　建　筑　技　术 2003年第1期(总第91期)

桁架结构静力测试

桁架结构静力测试 邬雨萱1450502 金永学1550873 1.工程背景： 钢桁架桥在现实中应用广泛，工程实例中有各 种各样的钢桁架桥。钢桁架桥一般为超静定结 构，以使桥更为安全。桁架杆件主要受轴向拉 应力或压应力而不受弯矩。因此可以最大限度 发挥材料的性能，让承受更大的力，因此其十 分适合于大跨度结构。如图所示就是一座钢桁 架桥。但是实际应用中的桁架桥的结点往往并 非全铰接，其中或多或少带有刚接特性，因此实际使用时桁架的受力与理论计算并不完全相同。桁架结构是现代工程结构中最常用的结构之一。在荷载作用下，桁架杆件主要承受轴向拉力或压力，从而能充分利用材料的强度，节省材料，减轻自重和增大刚度，同时，桁架结构还具有造型优美，坚固耐用，具有艺术性等特点，在现代工程实践当中得到广泛的应用。因此，桁架的设计和测试显得尤为重要。 1.实验目的： （1）设计并组装桁架结构；图1

（2）理论分析选定杆件轴力大小和方向； （3）了解应变片测量原理及使用方式； （4）测定桁架各杆件轴力大小，并与理论值比较； 2.实验内容： （1）桁架搭建：该桁架由24根265mm×10mm×5mm和90根190mm×10mm×5mm的钢杆通过螺钉连结起来。成型后效果如下图。 图一桁架实物图 (a) (b) (c) 图二节点构造图 （2）实验方案设计：杆件选择：在实验中，为了测得杆的轴力，我们选择了三种不同的杆件粘贴应变片。杆件位置及编号如下图所示： 杆件2

每个测点在杆件的正，反两面分别粘贴应变片，编号后，再引出导线，接入DH-3818静态应变测试仪上。将应变片粘贴在杆件两侧，目的是排除由于受力不在桁架所在平面内而造成的杆件弯曲对测试的影响。在实验处理数据时，应取两个读数的平均值作为杆件的应变值。 加载设计：因简支梁的挠度在力集中在梁中点时达到最大，所以我们将荷载加在桁架的中间位置。为了加载方便，我们把加载点设计在桁架的上弦点A 处。如上图所示。 （3）受力分析：该桁架结构有一定的对称性，在作受力分析图时我们只画结构的一半受力图： （4）操作步骤：a.在需要测量的杆件上贴好应变片，将各应变片导线接入DH-3818静态应变测试仪并用电烙铁焊接牢固； b.将DH-3818静态应变测试仪各通道清零并平衡； c.加载，记录下各通道的读数，计算轴力，与理论值进行比较。 （5）实验数据处理： 测得每个杆件的横截面都是10.25×3.30mm(取横截面积为34mm 2)的矩形，取弹性模量E=210G ，重力加速度g=9.8m/s 。 数据表如下： 其它杆件受力 外载荷杆件1受力杆件2受力杆件3受力

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