V-g法计算颤振临界速度.pdf

17-5-14下午4:16C:\Users\dkz\Documents\MATLAB\V-g.m第1页，共2页

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%用V-g法计算颤振频率

%第一问对应thequestion1

%第二问对应thequestion2

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

K0=[4,2,1.2,0.8,0.6,0.5,0.4,0.34,0.3,0.27,0.24,0.2,0.16,0.12,0.1,0.08,0.06];

Lh0=[0.9848-0.2519i,0.9423-0.5129i,0.8538-0.8833i,0.7088-1.3853i,...

0.5407-1.9293i,0.3972-2.3916i,0.1752-3.1250i,-0.0221-3.8053i,-0.1950-4.4333i,...

-0.3798-5.1048i,-0.5520-5.8424i,-0.8860-7.2760i,-1.3450-9.5350i,-2.0020-13.4385i,...

-2.4460-16.6400i,-3.0100-21.5100i,-3.7530-29.7333i];

Lalpha0=[0.4218-0.9423i,0.1858-0.9841i,-0.3823-1.5949i,-1.5228-2.2712i,...

-3.1749-2.8305i,-4.8860-3.1860i,-8.1375-3.5625i,-11.7140-3.7396i,-15.4730-3.7822i,...

-20.0337-3.6847i,-25.3190-3.5260i,-37.7665-2.8460i,-61.4370-1.1288i,-114.4920+3.

2420i,...

-169.3460+7.8200i,-272.4100+16.1150i,-499.8530+32.82222i];

Mh0=0.5;

Malpha0=[0.375-0.2500i,0.375-0.5000i,0.375-0.8333i,0.375-1.2500i,...

0.375-1.6667i,0.375-2.0000i,0.375-2.5000i,0.375-2.9412i,0.375-3.3333i,...

0.375-3.7500i,0.375-4.1667i,0.375-5.0000i,0.375-6.2500i,0.375-8.3333i,...

0.375-10.0000i,0.375-12.5000i,0.375-16.6667i];

%thequestion1

a=-0.4;

xalpha=0.05;

ralpha2=0.25;

omeh_omea=0.6;

miu=20;

%%thequestion2

%a=-2.0;

%xalpha=1.8;

%ralpha2=3.48;

%omeh_omea=1.0;

%miu=60;

fid=fopen(solution.dat,w);

fprintf(fid,%20s\n,VARIABLES=v,g,omega);

fprintf(fid,%20s\n,ZONEI=34,F=POINT);

n=length(K0);

fori=1:n

k=K0(i);

Lh=Lh0(i);

Lalpha=Lalpha0(i);

Mh=Mh0;

Malpha=Malpha0(i);

17-5-14下午4:16C:\Users\dkz\Documents\MATLAB\V-g.m第2页，共2页

A=[miu+Lh,