机械动力学试验报告Word下载.docx
《机械动力学试验报告Word下载.docx》由会员分享,可在线阅读,更多相关《机械动力学试验报告Word下载.docx(21页珍藏版)》请在冰豆网上搜索。
(三)实验框图
锤击实验框图:
(四)绘制实频、虚频、幅频、相频及Nyqiust圆图
1.绘制实频、虚频、幅频、相频及Nyqiust圆源程序清单:
%原始数据读入------------------------------------------------->
clc;
forj=1:
4
fid=0;
whilefid<
1
switchj
case1
filename='
wcry01.txt'
;
case2
wcry02.txt'
case3
wcry03.txt'
case4
wcry04.txt'
end
[fid,j]=fopen(filename,'
r'
);
iffid==-1
disp(message)
end
data_all=fscanf(fid,'
%g'
N=length(data_all);
%数据长度
fs=2560;
%采样频率
f=1000;
%分析频率
data_Re=data_all(1:
N/2);
%前半段数据为实频
data_Im=data_all(N/2+1:
N);
%后半段数据为虚频
fori=1:
N/2
Aw(i)=sqrt(data_Re(i)*data_Re(i)+data_Im(i)*data_Im(i));
%幅频
Qw(i)=atan(data_Im(i)/data_Re(i));
%相频
figure
subplot(221);
plot(data_Re(1:
400),'
b'
title('
—实频特性—'
)
subplot(222);
plot(data_Im(1:
g'
—虚频特性—'
subplot(223);
plot(Aw(1:
c'
—幅频特性—'
subplot(224);
plot(Qw(1:
—相频特性—'
400),data_Im(1:
k'
—Nyquist特性—'
第一组数据绘图:
第二组数据绘图:
第三组数据绘图
第四组数据绘图:
2.用五点求拟合圆,并计算各阶固有频率,阻尼比:
采用第一组数据计算:
源程序清单:
fprintf(‘----------采用第1组数据计算----------\n'
)
fid=0;
whilefid<
%原始数据文件名称
[fid,message]=fopen(filename,'
%读入数据
disp(message);
end
%fs=2560;
%f=1000;
peak_value=zeros(1,3);
subpoint=[14091130131170];
[max_data1,peak_value
(1)]=max(Aw(subpoint
(1):
subpoint
(2)));
[max_data2,peak_value
(2)]=max(Aw(subpoint(3):
subpoint(4)));
[max_data3,peak_value(3)]=max(Aw(subpoint(5):
subpoint(6)));
peak_value
(2)=peak_value
(2)+subpoint(3)-1;
peak_value(3)=peak_value(3)+subpoint(5)-1;
fprintf('
\n三阶峰值位置分别是:
'
fprintf('
%d%d%d\n'
peak_value
(1),peak_value
(2),peak_value(3));
m=5;
k=(m+1)/2;
dataX_curve=zeros(1,m);
dataY_curve=zeros(1,m);
3
peakva=peak_value(j);
[center_X,center_Y,curve_R,A_mainarray]=calc_curve(data_Re,data_Im,peakva);
ring_X=zeros(1,628);
ring_Y=zeros(1,628);
628
ring_X(i)=curve_R*cos(i/100)+center_X;
ring_Y(i)=curve_R*sin(i/100)+center_Y;
plot(ring_X,ring_Y,'
m'
holdon
m
dataX_curve(i)=data_Re(peak_value(j)-k+i);
dataY_curve(i)=data_Im(peak_value(j)-k+i);
plot(dataX_curve,dataY_curve,'
b-*'
%绘制拟合圆所取各点
一阶五点拟合园法'
二阶五点拟合园法'
三阶五点拟合园法'
x1=ring_X
(1)+curve_R/10;
%x轴绘制
y1=0;
x2=ring_X(314)-curve_R/10;
y2=0;
point_X=[x1,x2];
point_Y=[y1,y2];
plot(point_X,point_Y,'
'
LineWidth'
1);
y1=ring_Y(157)+curve_R/10;
%y轴绘制
x1=0;
y2=ring_Y(471)-curve_R/10;
x2=0;
plot(center_X,center_Y,'
mo'
%圆心坐标
%圆对称轴x
y1=center_Y;
y2=y1;
r--'
0.8);
%圆对称轴y
x1=center_X;
x2=x1;
holdoff;
axis('
image'
http:
chord_angle=zeros(1,6);
chord_anglem=0;
fori=2:
m+1
chord_angle1=atan(data_Im(peak_value(j)-k+i-1)/data_Re(peak_value(j)-k+i-1));
chord_angle2=atan(data_Im(peak_value(j)-k+i)/data_Re(peak_value(j)-k+i));
chord_angle(i)=abs(chord_angle2-chord_angle1);
[chord_anglem,peak_num]=max(chord_angle(1:
5));
%取得最大角度即为最大弦长位置
chord_num=peak_value(j)-k+peak_num;
ifpeak_num==5%如果固有频率为第5点则绘制下一点
x1=data_Re(peak_value(j)+k-1);
y1=data_Im(peak_value(j)+k-1);
x2=data_Re(peak_value(j)+k);
y2=data_Im(peak_value(j)+k);
3%绘制固有频率与圆心连线
x2=data_Re(peak_value(j)-k+peak_num-i+2);
y2=data_Im(peak_value(j)-k+peak_num-i+2);
end江西机械http:
natural_FREQ(j)=(peak_value(j)-k+peak_num)*2.56;
%得到各阶固有频率
pkvl_plus=(peak_value(j)+1)+(peak_value(j)-1);
pkvl_divid=(peak_value(j)+1)-(peak_value(j)-1);
%得到各阶阻尼比
damping_ratio(j)=(pkvl_plus-2*k+2*peak_num)*pkvl_divid/2/...
(peak_value(j)-k+peak_num)^2*(tan(chord_angle(peak_num)/2)+tan(chord_angle(peak_num+1)/2));
\n各阶固有频率分别是:
%输出各阶固有频率
disp(natural_FREQ);
\n各阶阻尼比分别是:
%输出各阶阻尼比
disp(damping_ratio);
程序清单:
function[center_X,center_Y,curve_R,A_mainarray]=calc_curve(data_Re,data_Im,peakva)
m=5;
dataX_curve(i)=data_Re(peakva-k+i);
dataY_curve(i)=data_Im(peakva-k+i);
curve_A=zeros(3,3);
curve_B=zeros(3,1);
curve_C=zeros(3,1);
curve_A(1,1)=dataX_curve*dataX_curve'
curve_A(1,2)=dataX_curve*dataY_curve'
curve_A(2,2)=dataY_curve*dataY_curve'
curve_A(1,3)=curve_A(1,3)+dataX_curve(i);
curve_A(2,3)=curve_A(2,3)+dataY_curve(i);
curve_B
(1)=curve_B
(1)-dataX_curve(i)^3-dataX_curve(i)*dataY_curve(i)^2;
curve_B
(2)=curve_B
(2)-dataX_curve(i)^2*dataY_curve(i)-dataY_curve(i)^3;
curve_B(3)=curve_B(3)-dataX_curve(i)^2-dataY_curve(i)^2;
curve_A(2,1)=curve_A(1,2);
curve_A(3,1)=curve_A(1,3);
curve_A(3,2)=curve_A(2,3);
curve_A(3,3)=m;
A_mainarray=curve_A;
curve_C=curve_A\curve_B;
center_X=-curve_C
(1)/2;
center_Y=-curve_C
(2)/2;
curve_R=sqrt(center_X^2+center_Y^2-curve_C(3));
%子程序结束
结果输出
----------采用第1组数据计算----------
三阶峰值位置分别是:
14326361
各阶固有频率分别是:
35.8400837.1200921.6000
各阶阻尼比分别是:
1.22610.00250.0017
三阶拟合圆图形:
2.求各阶主振型,仍然调用上述求圆心子程序
file_num=4;
main_array=zeros(file_num,3);
forks=1:
file_num
switchks
[fid,w]=fopen(filename,'
forj=1:
peak_flag=1;
ifdata_Im(peakva)<
0;
peak_flag=-1;
main_array(ks,j)=curve_R*peak_flag;
array_max=zeros(1,3);
array_max(i)=main_array(4*i-2);
main_array(j,i)=main_array(j,i)/array_max(i);
plot(main_array((4*i-3):
(4*i)),'
holdoff
switchi
一阶主阵型'
二阶主阵型'
三阶主阵型'
%程序结束
绘图及计算结果:
(五)简述模态分析的理论和方法
模态分析技术是用于对机械系统、土建结构、桥梁等工程结构系统进行动力学分析的现代化方法和手段。
它最早应用于航空、航天领域。
随着现代科学技术的发展,人们对工程产品的设计提出了愈来愈高的要求——如车辆、船舶的乘坐舒适性和噪声控制,产品轻量化设计的疲劳强度问题等,而产品结构的振动特性对此有着至关重要的影响,因此模态分析的应用领域日益扩大。
模态分析可以定义为,对结构动态特性的解析分析和试验分析,其结构动态特性用模态参数来表征。
表征模态的特征参数是振动系统的各阶固有频率、固有振型(主振型)、模态质量、模态刚度和模态阻尼等。
在数学上,模态参数是力学系统运行微分方程的特征值和特征向量;
而在试验方面则是试验测得的系统之极点(固有频率和阻尼)和振型(模态向量)。
模态分析技术可分为试验模态分析和计算模态分析两种方法。
试验模态分析采用理论分析与实验相结合的方法来识别系统的模态参数(模态频率、模态阻尼、振型),对系统的动态分析与优化设计具有实用价值,常用方法:
频带激振——锤击法:
锤击法是一种瞬态激振的方法,为非周期确定信号,并具有宽带连续频谱,可同时激出结构的各阶模态,是一种快速试验手段。
其缺点是输入能量有限,易导致响应数据的信噪比较低。
锤击激振是由带力传感器的敲击锤敲击结构实现的,这相当于给结构一个初速度。
锤头把宽频脉冲加给被测结构,同时激出所有各阶模态。
此法的设备简单,使用灵活方便,便于现场或在线检测;
由于在采样周期内响应信号基本衰减,所以一般没有因功率泄漏引起的偏度误差,因而锤击法在单点激振中得到广泛应用,特别适合轻型机构。