第五节 用整体坐标表示单元刚度矩阵.pptVIP

第五节 用整体坐标表示单元刚度矩阵 在装配结构总刚度矩阵时，必须把单元坐标表示的单元刚度矩阵变换到整体坐标下，本节将给出二者之间的关系。 1.公式推导 2.平面桁架单元整体坐标表示的单元刚度矩阵 3.平面刚架单元整体坐标表示的单元刚度矩阵 3.平面刚架单元整体坐标表示的单元刚度矩阵 平面刚架整体坐标表示的单元刚度矩阵程序设计 void change(double *ke, double *T) { int i, j,m; double k[6][6]； for(i=0; i6; i++) { for(j=0; j6; j++) { k[i][j]=0; for(m=0; m6; m++) k[i][j]+=T[m*6+i]*ke[m*i+j] ; } } for(i=0; i6; i++) { for(j=0; j6; j++) { ke[i*6+j]=0; for(m=0; m6; m++) ke[i*6+j]+= k[i][ m]* T[m*6+j] ; } } * * 单元坐标下的杆端力、杆端位移 ： 整体坐标表示的单元刚度方程： 单元刚度方程： 中的杆端位移向量 和杆端力向量 均应进行坐标变换，换成按整体坐标表示 、 。 单元坐标下的刚度方程 化为： 上式两边分别左乘 ， 且由 ： 得整体坐标表示的单元刚度矩阵： ? 由前面求得的平面桁架单元的刚度矩阵和坐标变换矩阵 令： Cx=cosα,Cy=sinα，且把上式代入 得： 上式即为平面桁架单元整体坐标表示的单元刚度矩阵 [例2-1]平面桁架如图所示，各杆截面EA均为常数。已知P1=15kN,P2=20kN，试桁架各杆轴力。 1.对结点和单元编号如图示; 2.? 列表表示各单元参数; 2.00 0.5 0.866 30° 3→2 ② 1.732 0 1 0 1→2 ① 单元长度l(m) sinα 即Cy cosα 即Cx α 单元坐标 x轴方向 单元 3.列出各单元刚度矩阵 x y 4.列出整体坐标表示的各单元刚度矩阵; 第一种方法： 直接代入公式： 单元（1）的单元坐标和整体坐标一致，所以 单元（2）的单元坐标和整体坐标不一致，必须经过以下变换 第二种方法： 利用坐标变换公式： 以上代入公式： 由前面求得的平面刚架单元的刚度矩阵和坐标变换矩阵 单元刚度矩阵和坐标变换矩阵是单元的固有属性，仅有单元的E, I, l, a有关。利用下式即可计算出整体坐标表示的单元刚度矩阵。 上式即为平面刚架整体坐标表示的单元刚度矩阵。一般不需专门记忆，只要记住单元坐标表示的单元刚度矩阵和二者之间的关系。 令： ， ，Cx=cosα,Cy=sinα，且把上式代入 得： 6m 6m 6m =20KN/m [例 2-3] 平面刚架如图所示，各杆截面相同。E=1×107kN/m2, A=0.24m2,I=0.0072m4，求各杆端力，并画出内力图。 1 2 3 1 2 3 4 [解] 1.对应结点及各单元编号如图所示; 0.12×105 4×105 0 1 0° 3→4 ③ 0.0849×105 2.8285×105 0.7071 0.7071 45° 2→3 ② 0.12×105 4×105 0 1 0° 1→3 ① Cy Cx α 单元坐标x轴 单元 2.列出单元参数表;