p8_12.m

振动工具箱

% This is problem 8.12 in the text solve using the FEM.
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=[];