有限元方法-02.ppt

2.2.6引入边界条件 2.2.6引入边界条件 乘大数法（对角线元素扩大法） 当边界条件不一定是零位移，而是已知值时，在总刚度矩阵中，把与给定节点位移对应的主对角线上的元素程乘以相当大的一个数，如1x1015，而该行该列上的其余元素不变。 在节点载荷列向量中，把相应的项用给定位移与相应的主对角线上的元素、同一相当大的数如1x1015这三项的乘积代替，而其余元素不变。 2.2.6引入边界条件 2.2.7解方程组求节点位移 高斯消元法 三角分解法 2.2.8求单元内力 求解得到的节点位移带到单元方程中返求由节点位移引起的节点力。 计算每个单元的内力、弯矩图和变形。 2.2.9平面刚架例题 P=20kN L2=2m L1=1m L2/2 1 2 3 4 1 2 3 2 3 2 1 2 1 3 4 3 y x x y y x (a) 平面刚架 (b)单元、节点编号与整体、局部坐标 （有限元直接法例题1） UY UX 2 3 2 1 2 1 3 4 3 UY UX 节点力和节点位移在整体坐标下的编号 F9 F7 F8 F6 F4 F5 F9 F7 F8 F6 F4 F5 F3 F1 F2 F12 F10 F11 2.2.9平面刚架例题 2.2.9平面刚架例题 由于外力作用于单元2上且不在节点上，因此只有节点2和3有等效节点荷载，根据2.2.4中等效节点载荷公式可以求得： P=20kN 1 2 3 4 1 2 3 U02 U03 V02 V03 M02 M03 节点位移引起的单元节点力之和等于节点荷载，据此建立各个节点的平衡方程，如下： 写成向量形式： 其中： 2.2.9平面刚架例题 单元1的刚度矩阵 F1 K11 K12 K13 K14 K15 K16 δ1 F2 K21 K22 K23 K24 K25 K26 δ2 F3 = K31 K32 K33 K34 K35 K36 δ3 F4 K41 K42 K43 K44 K45 K46 δ4 F5 K51 K52 K53 K54 K55 K56 δ5 F6 K61 K62 K63 K64 K65 K66 δ6 两节点(6个自由度)平面梁单元的刚度矩阵一般型式 材料力学超静定03 2)将原系统上的载荷以及多余未知力加在基本系统上，称为相当系统。 p B A Fy Fx Mz 超静定梁的相当系统 材料力学超静定04 3)要使相当系统能够代替原来系统，则两者在变形情况上应完全一致。即，相当系统在多余未知力作用处的位移能够满足一定条件——符合原系统在多余约束处的变形谐调条件。 Ux=0, Uy=0, θz=0 材料力学超静定05 线性叠加原理： p Fx Mz Fy + + + p B A Fy Fx Mz = Ubx=0, Uby=0, θbz=0 材料力学超静定06 柔度法——力法 D1=D1P+ δ11F1 + δ12 F2+ δ13 F3+… +δ1nFn Dn=DnP+ δn1F1 + δn2 F2+ δn3 F3+… +δnnFn 柔度方程与矩阵 D1 δ11δ12δ13δ14δ15…δ1n F1 D1P D2 δ21 F2 D2P = + Dn δn1 δnn Fn DnP 柔度法——力法 对于n次超静定结构，解除n个约束条件； 选择对应的n个未知力F1,F2,,,Fn； 将实际载荷作用在基本系统上，得到相应于n个未知力处的n个位移D1P,D2P,…,DnP。 将n个未知力的单位值分别作用在基本系统上。由每个单位力，确定相应于所有n个未知力处的n个柔度δij。 柔度δij：定义为由于Fj处力的单位值所引起相应于Fi力处的位移。 实际位移D1,D2,…,Dn等于载荷产生的位移与未知力产生的位移之和。（位移的协调方程） 柔度法——力法 p2 p1 p4 p3 p2 p1 p4 p3 f1 f2 f3 f4 p2 p1 p4 p3 d1p f1=1 d2p d3p d4p 4次超静定系统 相当系统 f2=1在各约束点产生的位移 f1=1在各约束点产生的位移 δ11 δ21 δ31 δ41 δ12 δ22 δ32 δ42 δ13 δ23 δ3