曲线箱梁弯桥的有限单元法
【专利摘要】本发明提供一种曲线箱梁弯桥的有限单元法，这种基本方法的思想的基础即“化整为零”，就是离散化结构。首先利用单元这一基本元素，系统的来实现用有限个未知量的数量去无线逼近无限个未知量。有限元模型由一些简单的单元组成并且单元之间通过荷载传递作用是通过节点相连接的，它就是一种数学抽象方法。其中构件空间中的坐标位置由节点表示，它有自由度，并且之间还有物理作用。一组节点自由度就形成了一个单元，用来表示相互的数值，可以用矩阵的形式表示。单元有如线、面、实体和二维或三维的多种样式。线性方程可以很好的体现每个单元的特性。整个构件的数学模型就是将单元作为一个整体。
【专利说明】曲线箱梁弯桥的有限单元法
[0001]【技术领域】:
本发明涉及有限单元法，尤其涉及曲线箱梁弯桥的有限单元法，其属于数据分析【技术领域】。
[0002]【背景技术】:
在连续体上直接进行近似计算就是有限单元法，它本质是一种数值方法，它就是利用数学办法对真实的物理系统进行模拟逼近。目前分析结构最常用的就是有限单元法，因为它能够处理任意的荷载类型和任意形状的复杂结构，所以对于曲线梁桥来说，无论什么结构形式均可以使用有限单元法进行分析。
[0003]
【发明内容】
:
本发明提供一种曲线箱梁弯桥的有限单元法，其通过将有限单元法应用于曲线箱梁弯桥中计算出曲线梁的内力。
[0004]本发明采用如下技术方案:一种曲线箱梁弯桥的有限单元法，其包括如下步骤: 步骤一、结构的离散化
所述离散化就是首先将结构用一定的单元形式划分成有限个单元体，然后把单元的指定点设为相邻连接单元的结点，形成单元的集合体，以集合体来代表原结构；
步骤二、确定位移模式
离散化工作完成后，合理假设单元中的位移分布，具体做法就是假设单元中任意一点的位移可用一个简单、合理的坐标函数来表示，这个坐标函数就叫做位移函数或位移模式，将多项式作为位移函数，多项式的微积分运算较为简单并且从泰勒级数展开的意义上来说，任意局部光滑函数都可以用多项式来进行逼近，工作的结果就是建立起下面所示的矩阵方程
,[明⑷
式中一单元中任一点的位移列阵；
[况]一形函数矩阵，其元素是坐标的函数；
⑷ 一单元的结点位移列阵；
步骤三、单元特性分析
确定了单元位移函数之后，便可以对单元做下面三个方面的工作:
1、运用几何方程也就是应变-位移关系，将单元中任一点的应变用结点位移表示，即建立了下面的矩阵方程
(神]⑷
式中(2*1—单元中任一点的应变列阵；
卜形变矩阵，一般其元素是坐标的函数；
2、运用物理方程也就是应力-应变关系，推导出单元应力矩阵方程用单元结点表示的
【权利要求】
1.一种曲线箱梁弯桥的有限单元法，其特征在于:包括如下步骤: 步骤一、结构的离散化 所述离散化就是首先将结构用一定的单元形式划分成有限个单元体，然后把单元的指定点设为相邻连接单元的结点，形成单元的集合体，以集合体来代表原结构； 步骤二、确定位移模式 离散化工作完成后，合理假设单元中的位移分布，具体做法就是假设单元中任意一点的位移可用一个简单、合理的坐标函数来表示，这个坐标函数就叫做位移函数或位移模式，将多项式作为位移函数，多项式的微积分运算较为简单并且从泰勒级数展开的意义上来说，任意局部光滑函数都可以用多项式来进行逼近，工作的结果就是建立起下面所示的矩阵方程
{ 8 } = [N] { 8 } 式中(8丨一单元中任一点的位移列阵； [刚一形函数矩阵，其元素是坐标的函数； .|86| 一单元的结点位移列阵； 步骤三、单元特性分析 确定了单元位移函数之后，便可以对单元做下面三个方面的工作: ^运用几何方程也就是应变-位移关系，将单元中任一点的应变用结点位移表示，即建立了下面的矩阵方程I ^ I = [8] I 8 6| 式中(8丨一单元中任一点的应变列阵； [8]—形变矩阵，一般其元素是坐标的函数； I运用物理方程也就是应力-应变关系，推导出单元应力矩阵方程用单元结点表示的 (0 | = [0] [8] ( 8 6| = [8] ( 8 6| 式中{0} —单元中任一点的应力列阵； [0]一与单元材料有关的弹性矩阵； [3]—应力矩阵，一般其元素是坐标的函数； 匕运用最小势能原理和虚位移建立如下的刚度方程 IV』+ ^} = 0^(6』 式中IV』一单元结点力列阵；
一单元等效荷载列阵，与作用在单元上的外荷载有关； [^] 一单元刚度矩阵，按下式计算；[[.1=广.1[事 步骤四、研究集合所有单元的特性推导整个结构的平衡方程 利用对号入座的直线刚度法集成整个结构的综合等效结点荷载列阵和整体刚度矩阵，进而建立构件的整体刚度方程[幻⑷= 式中[幻一结构整体刚度矩阵； (八丨一结构整体位移列阵；一结构综合等效结点荷载列阵； 步骤五、求解方程组和计算输出的结果 求出位移后，可以进行下一步的计算即计算应力及内力，随后用图形或数表的方式输出整理后的结果，在上述这些基础上再结合具体条件及问题来进行结构的设计。
【文档编号】G06F17/50GK103838913SQ201310725636
【公开日】2014年6月4日 申请日期:2013年12月25日 优先权日:2013年12月25日
【发明者】王赞芝, 辛立凤 申请人:广西科技大学