《有限元方法与MATLAB程序设计》_周克民 习题解答汇总 习题3.1--5.1.doc

PAGE 42 2 PAGE 29 - PAGE 42 - 附录C 习题答案 第3章 习题3.1 求图示阶梯变截面直杆单元刚度矩阵。 解答：将杆件分为2个等截面直杆单元，按照公式可以写成分块形式： 组装在一起得到结构刚度方程： 再分组展开，并假设中间结点没有荷载： 由第2式解得： 代入第1式： 最后得到单元刚度矩阵： 至此，需要把刚度子矩阵具体带入并整理。不过这些已经不需要力学概念了，都是繁琐得数学计算过程了，计算量还是挺大。可以用MATLAB代码（Ex3_1.m）推导如下： syms EI EA L real L1 = L/2; EI0 = [2*EI,EI]; for i =1:2 EI1 = EI0(i); kii = [ EA/L1,0,0;0, 12*EI1/L1^3, 6*EI1/L1^2; 0, 6*EI1/L1^2,4*EI/L1]; kij = [-EA/L1,0,0;0,-12*EI1/L1^3, 6*EI1/L1^2; 0,-6*EI1/L1^2,2*EI/L1]; kjj = [ EA/L1,0,0;0, 12*EI1/L1^3,-6*EI1/L1^2; 0,-6*EI1/L1^2,4*EI/L1]; k(:,:,i) = [kii,kij;kij,kjj]; end Z = zeros(3,3); i = 1:3; j =4:6; ke = [k(i,i,1),Z;Z,k(j,j,2)]-[k(i,j,1);k(j,i,2)]*((k(j,j,1)+k(i,i,2))\[k(j,i,1),k(i,j,2)]); disp(ke) 结果显示如下： [ EA/L, 0, 0, -EA/L, 0, 0] [ 0, -(64*EI)/(7*L^3), (16*EI)/L^2, 0, (64*EI)/(7*L^3), (16*EI)/(7*L^2)] [ 0, (16*EI)/L^2, 0, 0, -(16*EI)/L^2, (4*EI)/L] [-EA/L, 0, 0, EA/L, 0, 0] [ 0, (64*EI)/(7*L^3), -(16*EI)/L^2, 0, -(64*EI)/(7*L^3), -(16*EI)/(7*L^2)] [ 0, (16*EI)/(7*L^2), (4*EI)/L, 0, -(16*EI)/(7*L^2), (24*EI)/(7*L)] 习题3.2 图示杆件一段铰接，为简单记，忽略轴向变形。推导刚度矩阵。提示杆端力和杆端位移分别为： 解答：刚度方程为 分解为 注意，由第2式得： 带入第1式： 整理得： 得到： 整理得： 上述运算得MATLAB代码如下(Ex3_2.m)： syms EI L k = 2*EI/L^3*[6,3*L,-6,3*L; 3*L,2*L^2,-3*L,L^2;-6,-3*L,6,-3*L; 3*L,L^2,-3*L,2*L^2]; ke = k(1:3,1:3)-k(1:3,4)*k(4,1:3)/k(4,4); disp(ke) 输出为 [ (3*EI)/L^3, (3*EI)/L^2, -(3*EI)/L^3] [ (3*EI)/L^2, (3*EI)/L, -(3*EI)/L^2] [ -(3*EI)/L^3, -(3*EI)/L^2, (3*EI)/L^3] 习题3.3 求图示杆端力。 解答：形函数为 结点力公式为 （1）在分布力作用下结点力为 上述运算得MATLAB代码（Ex3_3a.m）如下： syms x L N = [(1-x/L)^2*(1+2*x/L),x*(1-x/L)^2,(x/L)^2*(3-2*x/L),x^2/L*(x/L-1)] 显示结果为 [ (13*L)/32, (11*L^2)/192, (3*L)/32, -(5*L^2)/192] （2）在集中力作用下结点力为 上述运算得MATLAB代码（Ex3_3b.m）如下： syms N(x) [1,4] syms L N(x) = [(1-x/L)^2*(1+2*x/L),x*(1-x/L)^2,(x/L)^2*(3-2*x/L),x^2/L*(x/L-1)]; dN(