桥梁博士系统(DR.Bridge)桥梁CAD——杨剑主讲桥梁CAD教学内容有限元概述2h直线桥设计计算输入4h直线桥设计计算输出2h设计计算工具2h上机练习6h第一章有限元法概论本章主要内容工程结构的基本概念土木工程结构的计算方法有限元概述平面杆系结构的有限元法1.1工程结构的基本概念工程结构原型,广义的看均为由无限多个质点所组成的三维连续固体,因而也就具有无限多个自由度的体系。考虑其具体的几何形状与应力-应变的特殊性,工程结构通常可划分:工程结构分类(1)杆系结构(2)二维结构(3)三维结构工程结构分类(1)杆系结构定义:由一定数量杆件通过一定数量结点相互连接而组成的结构体系○○○○○○○○○梁框架○○○○○桁架特点:横截面尺寸远比其杆长小,二结点的联结可能为铰接或刚接。分类:(a)平面杆系结构:全部杆件、支座及作用力均位于同一平面内;(b)空间杆系结构:全部杆件、支座及作用力不全位于同一平面内;工程结构分类(2)二维结构定义:当三维连续体z坐标向的应力或应变可以忽略时,按二维问题简化分析平面应力问题平面应变问题工程结构分类(3)三维结构定义:最一般的工程结构状态,其位移、应力、应变都是三维坐标x,y,z的函数。工程结构材料的本构关系(1)线弹性定义:结构物加载下的应变,在荷载卸除后将完全消除,从而恢复到结构未受载的原始状态,即:()Econst工程结构材料的本构关系(2)非线形弹性定义:结构物在整个受载变形过程中,应力-应变关曲线不再是直线而是曲线,即弹性模量E是变量。()E工程结构材料的本构关系(3)塑性定义:结构物加载时于卸载时的应力-应变关系不重合,有残余应变存在,且其应力-应变关系也为非线性。s工程结构材料的本构关系(4)其它本构关系弹塑性、粘弹性、粘塑性等多种本构关系模型江见鲸.《钢筋混凝土结构非线形有限元分析》.陕西科学出版社。工程结构的变形状态(1)小变形定义:结构受载后的位移与应变相当微小,与结构的原始尺寸相比可以忽略,则分析结构时可以以其未发生变形前的原始状态为依据,则其形变几何关系(位移-应变关系)是线性的。(2)大变形定义:结构受载后的位移与应变相当大,与结构的原始尺寸相比不可忽略,因此结构分析时必须计算结构几何形状的改变,即以结构变形后的现实状态作为计算依据,则其形变几何关系(位移-应变关系)是非线性的。EALEEXeq3222121混凝土斜拉桥的拉索一般为柔性索,在索的自重作用下有垂度,垂度对索的受拉性能有影响,同时索力大小对垂度也有影响。在实际计算中索一般采用一直杆表示,以索的弦长作为杆长。1.2土木工程结构的计算方法土木工程材料如木材、石料、混凝土、钢材等多为弹塑性材料,但从结构安全度的需要考虑,其工作状态通常都处于以弹性为主的阶段,这就使结构力学和弹性力学在土木工程结构分析中占有特别重要的地位。(1)结构力学:力法、位移法和混合法(2)弹性力学:基于二、三为连续体结构的几何条件、静力平衡与本构关系,按满足既定的边界条件来解析结构力学计算方法(1)力法-取结点力作为基本未知量111122133121122213323113223333000PPPXXXXXXXXXδδδ△δδδ△δδδ△结构力学计算方法(2)位移法-取结点位移作为基本未知量111122133121122213323113223333000PPPZZZZZZZZZrrrRrrrRrrrR结构力学计算方法(3)混合法-取一部分结点力和一部分位移作为基本未知量(4)三种方法的比较a.最终的数学表述均为多元线性代数方程组。b.超静定结构解算的繁简取决与超静定次数的多少;c.位移法将原结构最终简化为有限的几种基本杆件的集合,因而具有较强的通用性,便于实现程序标准化。(有限元法应用最广)d.力法的基本结构是与原结构形状相同的静定结构,因此,不同类型的原结构具有不同的形、载常数。弹性力学的解析方法1.研究弹性物体在外力和其它外界因素作用下产生的变形和内力2.主要应用于二、三维连续体结构问题3.其求解体系基于结构的几何条件、静力平衡与本构关系,最终演化的数学表述为偏微分方程,按满足既定的边界条件来解析徐芝纶:《弹性力学》,高等教育出版社。能量变分原理-有限元法分析的重要理论基础龙驭球:《结构力学》,清华大学出版社。1.3有限元法概述由于结构几何形状与边界条件的复杂多样性,依靠结构力学或弹性力学的直接解析来求取结果往往非常困难。通常最常用的是有限元法。有限元法进行的依据不再是结构弹性体原型,而是将其进行离散化处理的,由有限个单元在有限个结点相联结的替代结构,即“有限元模型”。1.3有限元法概述1.3.1有限元法的基本思想有限元法在20世纪50年代起源于飞机结构的矩阵分析,其基本思想是用有限个离散单元的集合体代替原连续体,采用能量原理研究单元及其离散集合体的平衡,以计算机为工具进行结构数值分析。它避免了经典弹性力学获得连续解的困难(建立和求解偏微分方程),使大型、复杂结构的计算容易地在计算机上完成,应用十分广泛。ANSYS,SAP,Marc,等。把整体结构离散为有限个单元,研究单元的平衡和变形协调;再把这有限个离散单元集合还原成结构,研究离散结构的平衡和变形协调。划分的单元大小和数目根据计算精度和计算机能力来确定。1.3有限元法概述1.3.2有限元法主要优点:(1)概念浅显,容易掌握。(离散、插值、能量原理、数学分析)(2)适用性强,应用范围广,几乎适用于所有连续体和场问题的分析。(结构、热、流体、电磁场和声学等问题)(3)计算规格化(采用矩阵表示),便于计算机编程。1.3有限元法概述1.3.3有限元法的基本构想:1.3有限元法概述(1)结构离散化:将原型结构划分为有限个“单元”(如直杆元、曲杆元、矩形元等),并通过有限个“结点”相互连接,从而形成称为“有限元模型”的替代结构。将作为今后分析的物理依据。1.3有限元法概述(2)单元特性计算:建立各单元结点的广义位移(轴向位移、切向位移、挠曲转角、扭转转角)与相应广义位移方向的结点内力(轴力、剪力、弯矩、扭矩)之间的关系。(即建立单元刚度矩阵)eeeFK1.3有限元法概述(3)有限元模型解析:把有限元模型所有单元的上述关系总集起来,将形成线性代数方程组形式的结构总刚度方程。FK(a)边界条件:在结构总刚度方程里引入原结构边界结点的相关约束条件。(b)结构连续性要求:汇交于i点各单元在结构受载变形后仍汇交于i点。即汇交于i点的各单元的任一广义位移方向上应有相同的变位。一维拉杆•图示阶梯形直杆,各段长度均为,横截面积分别为3A,2A,A,材料重度为γ,弹性模量E。•离散化:将单元划分为3个单元,4个结点。•单元刚度矩阵:2111113][)1(lAEk123211112][)2(lAEk23431111][)3(lAEk341111][lEAke•等效结点荷载:按静力等效原则,有:1123][)1(lAF1122][)2(lAF112][)3(lAF•对号入座,组成总刚,形成整体结构平衡方程:}{}]{[FK•整体结构平衡方程lAlAlAlAFuuuulEA21)2122()2223(23110011220022330033143211.3有限元法概述解综合方程[K]{⊿}={F}求结构节点位移{⊿}计算结构内力和应力系统分析(把单元刚度矩阵集合成结构刚度矩阵[K]形成等价节点荷载{F})离散结构为若干单元单元分析(建立单元刚度矩阵[k]e形成单元等价节点力)龙驭球:《有限单元法》1.3有限元法概述1.3.4有限元解的精确分析有限元法的解只是一种逼近原结构真实解的近似解答,解的精确性主要取决于以下方面:(1)有限元模型能否正确反映原结构的真实状态。(a)有限元模型是多单元组成的离散结构,其各单元之间仅通过结点连接,因此,内力从一单元向另一单元的传递也自然通过结点来进行。这种结点连接、结点传力的假设对于杆系结构来说,是复合杆件连接与传力的实际情况。故有限元法应用于分析杆系结构,其解为精确解。1.3有限元法概述(b)二、三维连续体:上述结点连接于传力的假定将歪曲连续体的连接与传力实际情况。连续体内,相邻单元内力的传递是通过其共有边界来进行,但在有限元模型里系假定为结点传递。因此,有限元法用以分析二、三维连续体,从本质上讲只能是近似解。有限元法的解能否真正收敛于原结构的精确解?1.3有限元法概述1.3.5单元的网格划分与类型选择有限元法模型里,单元网格划分的疏密将决定其单元尺寸与数目的大小。网格划分越密,有限元模型与原结构的形状、特性愈加接近,从而提高其近似解的精度。根据应力梯度使网格的布局合理化。即在梯度大的区域网格密些,梯度小的区域应稀些。密、稀网格之间应逐步过渡。1.3有限元法概述通常采用的单元类型可分为:(1)非协调单元:不能完全保证相邻单元间位移连续性的单元(2)协调单元:能完全保证相邻单元间位移连续性的单元(3)高阶协调单元:不但能完全保证相邻单元间位移连续性,而且还能保证单元间应变的连续性但有限元法的解能否真正收敛于原结构问题的精确解,单纯依靠增加单元数目、加密单元网格并不一定能实现,这就需要考虑单元的类型选择问题。1.3有限元法概述1.3.5桥梁结构的有限元程序1.桥梁博士,同济大学2.公路桥梁结构设计系统(GQJS),交通部公路科学研究所3.MIDAS/Civil(空间有限元程序),韩国1.4杆系结构有限元概述平面杆系结构有限元法在土木工程设计中运用非常广泛,连续梁、连续刚构、桁架拱、斜拉桥等,均可概化为平面杆系结构来分析。实际结构离散为有限元模型时,其结点的选取应遵循下列原则:(1)原结构的杆件自然交汇点作为结点;1.4杆系结构有限元概述(2)根据计算精度与验算截面需要,可适当增收结点,(a)变截面杆可用一系列等截面杆替代1.4杆系结构有限元概述(b)斜拉桥主梁的索间跨中截面增设结点1.4杆系结构有限元概述(3)实际结构的曲杆可使用一系列直杆元离散处理THEENDTHANKYOU!桥梁博士简介博士系统一个集可视化数据处理、数据库管理、结构分析、打印与帮助为一体的综合性桥梁结构设计与施工计算系统。其基本功能:(1)直线桥:能够计算钢筋混凝土、预应力混凝土、组合梁以及钢结构的各种结构体系的恒载与活载的各种线性与非线性结构响应(2)斜、弯和异型桥梁:采用平面梁格系分析各种平面斜、弯和异型结构桥梁的恒载与活载的结构响应桥梁博士简介(3)其它直线桥设计计算输入利用桥梁博士系统进行设计计算一般需要经过:直线桥设计计算输入(1)离散结构划分单元(2)施工分析(3)荷载分析(4)建立工程项目(5)输入总体信息、单元信息、钢束信息、施工阶段信息、使用阶段信息、输入优化信息(索结构)(6)进行项目计算(7)输出计算结果读懂图纸或明确设计要求将设计要求变成有限元模型查看结果单位约定坐标系平面杆系总体坐标系:系统默认的坐标系,节点坐标、节点位移以及反力均按总体坐标系输出。X:水平向右为正Y:垂直X轴向上为正单元局部坐标系:单元内力和应力均按单元局部坐标系输出。X:沿构件的纵轴线方向,以左节点到右节点方向为正Y:垂直X轴向上为正荷载方向水平力:沿整体坐标的x方向向右为正;竖直力:沿整体坐标的y方向向上为正;弯矩:依右手螺旋法则,垂直于整体坐标系向外(向用户方向)为正效应方向轴力:使单元受压为正,受拉为负剪力:由单元底缘向顶缘方向为正,反之为负弯矩:使单元底缘受拉为正,上缘受拉为负(平面)位移:与总体坐标系一致为正,反之为负正应力(法向应力):压应力为正,拉应力为负;剪应力:由截面底缘