%%%%%%%%%质量阵和刚度阵分块 mii=Me(1:2,1:2);mij=Me(1:2,3:4); mji=Me(3:4,1:2);mjj=Me(3:4,3:4); kii=Ke(1:2,1:2);kij=Ke(1:2,3:4); kji=Ke(3:4,1:2);kjj=Ke(3:4,3:4); %总质量阵和刚度阵组装 o=zeros(2,2); M=[mii mij 0 0 0 0 0 0 0 0; mji mjj+mii mij 0 0 0 0 0 0 0; 0 mji mjj+mii mij 0 0 0 0 0 0; 0 0 mji mjj+mii mij 0 0 0 0 0; 0 0 0 mji mjj+mii mij 0 0 0 0; 0 0 0 0 mji mjj+mii mij 0 0 0; 0 0 0 0 0 mji mjj+mii mij 0 0; 0 0 0 0 0 0 mji mjj+mii mij 0; 0 0 0 0 0 0 0 mji mjj+mii mij; 0 0 0 0 0 0 0 0 mji mjj]; %%%%%%%%%%%%%%总刚度阵 K=[kii kij 0 0 0 0 0 0 0 0; kji kjj+kii kij 0 0 0 0 0 0 0; 0 kji kjj+kii kij 0 0 0 0 0 0; 0 0 kji kjj+kii kij 0 0 0 0 0; 0 0 0 kji kjj+kii kij 0 0 0 0; 0 0 0 0 kji kjj+kii kij 0 0 0; 0 0 0 0 0 kji kjj+kii kij 0 0; 0 0 0 0 0 0 kji kjj+kii kij 0; 0 0 0 0 0 0 0 kji kjj+kii kij ; 0 0 0 0 0 0 0 0 0 0 0 kji kjj]; Mr=M(3:20,3:20); Kr=K(3:20,3:20); 错误使用 horzcat 要串联的数组的维度不一致。 出错 yadianliang (line 69) M=[mii mij 0 0 0 0 0 0 0 0; |
1 条回复