有限元 5-有限条法.docVIP

第5章 有限条法 5.1 引言 一、 发展概况 有限条法(Finite Strip Method)诞生于二十世纪60年代，一般认为主要创始人有：Y.K.Cheung(张佑启)教授和G.H.Powell(鲍威尔)、 D.W.Ogden(奥格登)两人。Y.K.Cheung在1966～1969年间首先用有限条法研究了矩形薄板弯曲问题，后两人开始于板式桥梁的研究工作。 二、有限条法的力学模型 有限条法可看作是有限元的一种特殊形式或分支，是一种（有限元）半解析法，适应于一些量大面方的，常用的规则结构形式，采用有限条法可使弹性力学中的二维问题化为一维问题（三维化二维），使总刚方程降阶，从而提高效率。 象有限元一样，有限条法亦需将连续体离散化，所不同的是，不象有限元一样可沿任意方面离散，而只能沿某一方向。如图示矩形板，用有限元分析（矩形元）的网格划分如右图示，而有限条则是沿x方向等分成若干条带。 有限条：x方向采用多项式插值函数 （梁函数） y方向采用三角级数表示： 然后板的位移函数采用一总和函数表示： 5.2 梁函数和基本函数 一、梁函数 梁函数用以表示条元的横向变化规律。图示梁有两个结点(i,j), 每个结点两个位移: 线位移(挠度)、； 角位移、 任意点的位移函数: 代入边界条件可得： [L]为在第二章中推导出的平面梁单元的形函数,此处称梁函数。 二、基本函数 基本函数用来表示条元的纵向变化规律，Y.K.Cheung将基本函数取为振动梁微分方程（按振动梁微分方程的规律变化）: 根据结构力学中“梁的无限自由度振动”理论可知，其解的一般形式为: (5-2-1) 在结构力学中，上式又称符拉索夫振动梁函数，式中C1～C4 为由两端边界条件确定的待定常数，由不同的支承条件，可得出相应的基本函数。下面介绍几种常用支承条件的基本函数。 1．两端简支 边界条件： Y(0)=Y(0)=0 Y(L)=Y(L)=0 代入式(5-2)得: 或者 =π, 2π, 3π……mπ 取, m=1, , 2. 一端简支一端固定 边界条件: Y(0)=Y(0)=0； Y(L)=Y(L)=0 特解(基本函数): ，,，, （m=5,6,…）(特征方程 tgμ=thμ)。 。谐波图与简支梁相似。 3. 两端固定 边界条件: 基本函数: ，，，，，（m=5,6,…） 5.3 弯曲板有限条法 一、位移函数 图a示简支板,在任意荷载作用下,将产生光滑的变形曲面(位移场), x, y方向的变形曲线分别如图示,y方向变化规律取为基本函数,x方向变化规律取为梁函数。 现在，沿板的支承方向将其离散成有限个矩形板条。从中取出一典型条元e加以研究，当条元宽度b较窄时, 其位移场可采用多项式f(x)和基本函数Y(y)的组合形式表示,如是可得 挠度 ： (5-3-1) 式中f(x)为描述横向位移变化规律的多项式函数，即，梁函数。Y(y)为摸拟板条纵向(y方向)挠度曲线的基本函数。代入得： 其横向转角： (5-3-2) 设条元结线i的挠度为,转角；结线j挠度为,转角(如图b)。 定义四个结线位移的表达式为: (5-3-3) 式中: ，，， 为第m个谐波所对应的结线i, j的挠度位移幅值和转角位移幅值。 将5-3-3代入5-3-1、5-3-2的左边并消去得: 由此解得: 将其代回5-3-1得用结线位移幅值表示的位移场： (5-3-4) 式中: 是对应于第m个谐波的结线位移幅值列阵。 5-3-5是对应对于第m个谐波的形函数。它与梁函数[L]相比只是增加了乘子“”。 于是板条中任意点的位移用矩阵形式表示为： (5-3-4’) 二、(沿)结线(i, j)的位移(函数) 实际应用时，通常只计算结线上的位移，将x=0和x=b代入得： 当x=0时, , x=b时, , 所以沿结线i, j的位移: 式中[I]为四阶单位矩阵:[I]=[1 1 1 1] 即:结线上沿y方向任意点的位移，是通过由第m个谐波解出的两结线的位移幅值 乘以相应基本函数 得到的，其最后的位移是r个谐波的叠加。 三、条元应变与几何矩阵[B] 由平板理论知： 式中[B]即为弯曲板条的几何矩阵： 子矩阵 的展开形式为: 若取简支条元, 则式中: , , 四、条元内力 弹性矩阵[D]与上章相同，弯、扭矩的正方向如图示，且应理解为单位长度上的量。 五、条元刚度矩阵的一