3 梁单元(新).pptVIP

(6) 节点位移与节点力的关系 例 变截面梁 有一变截面梁，一端固定，另一端铰支。梁长为2l，固支端的截面尽寸为b×1.6h，铰支端的截面尺寸为b×h。梁上作用均布载荷p0。求梁端的约束反力。 离散化 将梁划分成2个单元，3个节点。每个单元 长度为,截面取平均截面。 对号入座，组合整体刚度矩阵 荷载等效节点力向量 引入边界条件 * 2.4 梁单元 1、两端承受剪力、弯矩的平面梁单元 i j x y i j x y ?1 ?2 ?3 ?4 l F1 F2 F3 F4 l （1）局部坐标下单元位移和单元力 ① 单元位移 （2-24） 节点位移 节点载荷 其中， v——y方向位移，即挠度。 ?——角位移。 ② 单元力 （2-26） 其中， Q——剪力 M——弯矩 （2-27） （2）位移函数和形函数 （2-28） ① 位移模式 设单元坐标位移模式为 ② 形函数 由单元两端点的节点位移条件，解出式（2-28）中的a1、a2、a3、a4。再代入该式，可将位移模式写为以下形式： i j x y ?1 ?2 ?3 ?4 l 梁单元内一点有2个位移： v、? 因为， ?=dv/dx；仅一个位移是独立的，取 v 。 （2-29） 式中 （2-30） （2-31） （3）应变矩阵 ① 单元弯曲应变?b与节点位移???e的关系。 梁单元上任一点的应变和该点挠度之间关系为： 1 x y （2-32） 1 y ? ? ? 将式（2-29）代入（2-32），得单元弯曲应变和单元位移之间关系 （2-34） （2-33） （4）应力矩阵 （2-35） (5) 等效节点力 对于梁上作用的集中力或集中力矩，在划分单元时可将其作用点取为节点，按结构的节点载荷处理。 这里仅考虑把单元上的横向分布载荷转化为等价节点力问题。 x y i j l py(x) （2-36） 功互等： 注意： 将形函数矩阵[N]代入上式，积分可得分布荷载的等效节点力。表1给出了几种特殊情况的等价节点力。 - 5ql2/96 ql/4 5ql2/96 ql/4 - ql2/20 7ql/20 ql2/30 3ql/20 - ql2/12 ql/2 ql2/12 ql/2 Mj Qj Mi Qi 荷载分布 i j q q i j q i j 几种横向分布荷载等价节点力 表 1 由梁弯曲的挠曲线方程 i j a 2 = T 确定待定系数 同理：在x=l处 由梁弯曲的剪力 弯矩微分关系得： 在节点上 ： (7) 单元坐标单元刚度矩阵 得单元坐标单元刚度矩阵[k]e 用矩阵表示的梁端力与梁端位移的关系 单元刚度矩阵式(2-38)适合于连续梁分析。 （2-38） 整体坐标与局部坐标方向一致。 x y ， 单元刚度矩阵 i j 1 2 2 3 1 2 3 约束反力向量 1 2 3 总荷载向量 将整体平衡方程中对应的1、2、5行和总刚中1、2、5列删去 ，得 解方程组，得节点位移值 将节点位移值代入整体平衡方程，可得约束反力 2、两端承受轴力、剪力、弯矩的平面梁单元 （平面刚架，BEAM3) i j x y i j x y ?2 ?3 ?5 ?6 l ?1 ?4 F2 F3 F5 F6 l F1 F4 （1）单元坐标单元位移和单元力 ① 单元位移 （2-39） 其中， u——x方向（轴向）位移。 v——y方向位移，即挠度。 ?——角位移。 ② 单元力 （2-40） 其中， N——轴向力 Q——剪力 M——弯矩 对于小变形问题，可以认为轴向变形和弯曲变形互不影响，因此，位移模式和形函数可以分别按一维拉压杆单元和弯剪平面梁单元的结果（式2-3和式2-28）简单集合而成。 （2）位移函数和形函数 ① 位移模式 i j x y ?2 ?3 ?5 ?6 l ?1 ?4 （2-41） ② 形函数 式中形函数[N]为： （2-42） （2-43） 其中, （3）应变矩阵 ① 单元弯曲