有限元5-有限条研究.ppt

有限单元法 第五章 有限条法 5.1 引言 一、 发展概况 对于许多具有规则几何形状和简单边界条件的结构来说，完全的有限元分析常常是既浪费又不必要，有时甚至不可能（指超出计算机能力）。（张佑启《结构分析的有限条法》） 有限条法(Finite Strip Method)诞生于二十世纪60年代，一般认为主要创始人有：Y.K.Cheung(张佑启)教授和G.H.Powell(鲍威尔)、 D.W.Ogden(奥格登)两人。Y.K.Cheung在1966～1969年间首先用有限条法研究了矩形薄板弯曲问题，后两人开始于板式桥梁的研究工作。 二、有限条法的力学模型 象有限元一样，有限条法亦需将连续体离散化，所不同的是，不象有限元一样可沿任意方面离散，而只能沿某一方向。如图示矩形板，用有限元分析（矩形元）的网格划分如右图示，而有限条则是沿x方向等分成若干条带。 三、有限条法与有限元法的比较 四、有限条法的解题步骤 5.2 梁函数和基本函数 一、梁函数 二、基本函数 几种常用支承条件的基本函数 1．两端简支 边界条件： Y(0)=Y(0)=0 Y(l)=Y(l)=0 代入式(5-2-1)得: 5.3 弯曲板有限条法 一、位移函数 图示简支板，在任意荷载作用下，将产生光滑的变形曲面(位移场)， x、 y方向的变形曲线分别如图示，y方向变化规律取为基本函数，x方向变化规律取为梁函数。 现在，沿板的支承方向将其离散成有限个矩形板条。从中取出一典型条元e加以研究，当条元宽度b较窄时，其位移场可采用多项式f(x)和基本函数Y(y)的组合形式表示，可得挠度 ： 设条元结线i的挠度为wi，转角为θi；结线j挠度为wj，转角为θj，定义四个结线位移的表达式为: 将5-3-3代入5-3-1、5-3-2的左边并消去Ym(y)得: 将其代回5-3-1得用结线位移幅值表示的位移场 二、(沿)结线(i, j)的位移(函数) 实际应用时，通常只计算结线上的位移，将x=0和x=b代入得： 当x=0时, 当x=b时, 式中[I]为四阶单位矩阵 三、条元应变与几何矩阵[B] 由平板理论知： 子矩阵[B]m 的展开形式为: 四、条元内力 弹性矩阵[D]与上章相同，弯、扭矩的正方向如图示，且应理解为单位长度上的量。 五、条元刚度矩阵的一般形式 由单元刚度矩阵的一般表达式，条元的刚度矩阵: 其中子块: 式中：Dx , Dy , Dxy , D1 为各向异性板的弹性系数, 当各向同性时有: 5.4 用有限条法分析简支弯曲薄板 一、基本函数 由前已知，将符拉索夫振动梁函数微分方程： 二、条元刚度矩阵与条元刚度方程 1．条元刚度矩阵 上节已导得条刚的一般形式为： 可以证明，在简支条元中，由于基本函数 的正交性，当m≠n时，上述五个积分均为0。 因此，在简支条元中，非对角元子块均为零，即[S]mn =[0](m≠n)。此时，简支板的条刚可简化为： 其中: 由于条刚不耦联，便可一个谐波一个谐波的单独计算，然后再叠加，因此，简支板是最能体现有限条法优越性的一种情况。另外，从条刚中的元素表达式可见，对应于第m个谐波的条刚，每个元素均随着谐波号发生变化，这也是与前面单元刚度矩阵的不同之处。 3. 条元刚度方程 条刚[S]中的全部非对角元素(子块)为零，不仅简化了条刚，更重要的是，使得第m个谐波的结线力只与第m个谐波的结线位移有关，即总刚的非藕联性。因此，可对每个谐波所对应的结线力（位移）单独计算，然后再叠加，所以在下面的讨论中，我们可将条刚暂时理解为4×4阶的。对应于第m个谐波的条元刚度方程为： 三、总刚度矩阵及总刚方程的建立 总刚的形成仍然采用前述直接刚度法，即由各结线的平衡条件建立平衡方程，其过程亦与前述大致相同，现简述如下： 首先将对应第m个谐波的条刚按结线分块： 条刚中的每个子块均为2×2。子块具有两重下标，内下标ij指示子块在总刚中的结线位置，外层下标指示子块对应的谐波号。按结线将其展开得: 例：图示简支板，将其划分成五个条元。根据结线的平衡条件有：{P2}m={F2}m①+{F2}m② 由式可见，只有相交于结线i的两个板条才可能对总刚矩阵i行子块有“贡献”。因此，形成总刚的某行时，只需考虑对应结线左、右条刚的子块。而与其它条元的条刚无关，由于这一特性，使得简支板总刚的带宽很小（半带宽仅为4），与板单元比较， 若上例采用5×5网格，则未知数n=6×6×3=108, 半带宽D=3×(7-1)=18，规模要大得多。对每一谐波求解，可得第m个谐波所对应的结线位移，重复上述过程，将1～