变截面梁弯曲问题的样条有限元解法.pdfVIP

第九届全国建筑工程计算机应用学术会议论文集 黑龙江．哈尔滨1998 变截面梁弯曲问题的样条有限元解法 刘叶丹 赵琼梅 (广州大学) (广东省建筑工程学校) 【摘要】本文用三次B样条有限元法，讨论了变截面粱的弯曲问匙。对刚度(它是长度x的连续函数) 用一次样条插值函数适近，得出了，挠度函的基本表达武，推导出了适用一卜菩种边界条件的统一计算格式 井得出了有关的所有系数矩阵，非常便于在计算机上实现。 关键词：样条，变截面梁 一、刖舌 对于变截面粱，由于其惯性矩为x的函数．其表达式往往是复杂的，其挠度曲线所满 足的微分方程在求解上是很麻烦的，所以用一般的数学方法很难求解。本文以变分原理和 三次B样条函数为基础，并对变刚度做一次样条插值，推导}|{了适合于求各种边界条件 的梁的挠度曲线的统一计算格式。这里梁的刚度必顽是连续变化的。本方法中，近似后的 刚度函数仍是连续的，故用不着象其他方法如【2l中『自于将刚度函数近似成不连续的函数 而采用在间断点上附加连续性条件来处理。本方诘具有讣算餐少，精度高，通用性强等特 点。另外本文的结论也可作为分析变刚度板弯曲蒯题的个引导。 二、样条函数 三次样条函数： 似+2)3 x∈[一2，一1】 (x+2)3—4(x一1)3 r∈卜1，0] 仍(Y)=÷(2一x)’一4(1一x)。 z∈(O，1] 0 f2一n’ x∈【1，2】 O H≥2 苎垄壁全璺壅苎三堡茎墨垫生璺兰查鱼堡垒墨叁 墨垄兰：!鱼型!堕 三、基本方程 设粱的纵长在X轴向，弯曲发生在XY平面内。设梁变形前，轴线与X轴一一致，两 的横向荷载线密度，并设粱的抗变刚度D是x的连续函数(如变截面梁)。那么，在小变 形情况下，梁的应变能为： t，(V)=；fD(f)矿。’(r)】。出 外功势能为： f’(V)=一l『似)V@)at 梁的总能量： 矿(V)=J(V)+F(V)=：fD矿”dx—f触 两端边界条件最常见的是： l、固支边：V=0，V’=0 2、简支边：V=0，V”=0 3、自由边：V”=O(DV”)’：：0 对于一般的边界支承，只须在应变能中增加一项 I竺∥+旦。：l L2 2 J，。 或 陟吩产L 例如：若d=0，B=常数，表示在支承处承受弹性反力；若a=常数，B=o表示 在支承受弹性反弯矩。 由最小势能原理， 平衡态的位移使总势能V(v)达到最小值 Mi疗矿(v1 (1) 、／ 、’ V(H； 其中HE2是指[o，a]中二次导数平方可积、(有限能量)并满足强加边界条件的函 数集合。在以上所列举的各种边界条件中，只有固支边V=o，V‘=o及简支边中V= 56 第九届奎国建筑工程计算机应用学术套议论文集