p8_12

来自「vtoolbox3.5是振动数值模拟工具箱」· 代码 · 共 38 行

node=[0 0;.1 0;.2 0;.3 0;.4 0;.5 0;.6 0;.7 0;.8 0;.9 0;1 0];
ncon=[1 2 7e10 2.988e-2 1 0 2700;
      2 3 7e10 2.690e-2 1 0 2700;
      3 4 7e10 2.407e-2 1 0 2700;
      4 5 7e10 2.140e-2 1 0 2700;
      5 6 7e10 1.889e-2 1 0 2700;
      6 7 7e10 1.653e-2 1 0 2700;
      7 8 7e10 1.433e-2 1 0 2700;
      8 9 7e10 1.229e-2 1 0 2700;
      9 10 7e10 1.041e-2 1 0 2700;
      10 11 7e10 8.679e-3 1 0 2700];
ncon(:,7)=ncon(:,7).*ncon(:,4);
zero=[1 1;
      1 2;
      1 3;
      2 2;
      2 3;
      3 2;
      3 3;
      4 2;
      4 3;
      5 2;
      5 3;
      6 2;
      6 3;
      7 2;
      7 3;
      8 2;
      8 3;
      9 2;
      9 3;
      10 2;
      10 3;
      11 2;
      11 3];
force=[];
conm=[];