关于冲击响应谱

dreamstone

hcharlie大哥，有没有相关冲击响应谱由时域到频域的相关计算资料

hcharlie

给你一个我以前写的冲击响应谱计算公式，（见10楼）你可以自己编程。

[ 本帖最后由 hcharlie 于 2009-11-16 16:02 编辑 ]

xyxiang1

做冲击试验时,如果用加速计测量关键部件,例如:洗衣机某部件,通过测回来的波形,进行分析计算,得出这个部件的冲击响应谱.不就可以用冲击响应谱代替冲击脆值来描述了.

hcharlie

是可以用冲击响应谱综合来进行冲击试验，但一般控制仪只能接收比较简单的冲击谱，复杂的拟合不了。所以你用此法时，必须将你测得的谱加以简化。

dreamstone

谢了，你那个式子是不是再说不断的通过频率w的变化求出对应的不同Ymax，我的意思是已经得到了一个连续结构（比如说是板或者梁）在冲击下的响应的时域，我的工件要在这个环境下工作，而我是想通过冲击响应谱的频域来控制振动台实现这个环境，在进行反馈控制中，加速度传感器测量的是一个响应的时域信号，这个时域信号就是我的控制目标，我想得到这个冲击响应的频域，如何实现

[ 本帖最后由 dreamstone 于 2009-11-17 17:25 编辑 ]

hcharlie

用10楼的公式算，有什么问题吗？

dreamstone

谢谢了，我懂了，，终于弄明白了数天来笼罩在我头上的乌云是咋回事了，谢谢了

loren_xia

我最讨厌理论的东西了，其实很简单的一件事，为什么要复杂化呢？:lol

dreamstone

最近我像了解下，国内关于冲击响应谱的试验规范，找了几个地方没找到，你有相关资料么

xingzhe

近日也在找资料，见到这篇文章详细讲解如何求冲击响应谱：
http://www.doc88.com/p-5780308071.html

zqtian

　　一、冲击响应谱试验实现方式
　　1.用振动台实现，振动台控制仪软件中专门有这个模块。
　　2.用专用的冲击响应谱试验台，国内常用的是摆锤式冲击台。
　　3.利用谐振板的方式，就是把一块钢板四点固定，垂直于板面方向撞击板，产生振动。
　　二、冲击响应谱原理
　　冲击响应谱是一系列固有频率不同的二阶系统受到同一冲击激励的总结果。冲击响应谱绘制方法为，其横坐标为频率，纵坐标为具有该固有频率的二阶系统在此激励下的最大响应加速度。
　　[ 本帖最后由 zqtian 于 2010-3-16 23:57 编辑 ]

dreamstone

hcharlie大哥，我用龙格库塔算法写了个冲击响应谱的计算程序，对于半正选来说，最大值应该出现在71hz左右，而我的出现在35hz左右，怎么也弄不明白，
clear all
t1=0:0.00001:0.005;
n1=length(t1);
x1=zeros(1,n1);
t2=0.005:0.00001:0.015;
x2=sin(pi*(t2-0.005)/0.01);
t3=0.015:0.00001:0.02;
n3=length(t3);
x3=zeros(1,n3);
t=[t1,t2,t3];
x=[x1,x2,x3];
figure(1)
plot(t,x);
nn=length(x);
h=0.00001953125;
c=0.05;
u2=x;%输入加速度
fz=1;
y=zeros(1,nn);%响应位移
y1=zeros(1,nn);%响应速度z(t)
y2=zeros(1,nn);%响应加速度
d=[];
y2max=[];
gmax=0;
dmax=0;
i=1;
wn=2*pi*fz;
while (1)
       for j=2:nn
        k1=y1(j-1)*h;
        l1=(-u2(j-1)-2*c*wn*y1(j-1)-wn*wn*y(j-1))*h;
        k2=(y1(j-1)+0.5*l1)*h;
        l2=(-(u2(j-1)+0.5*(u2(j)-u2(j-1)))-2*c*wn*(y1(j-1)+0.5*l1)-wn*wn*(y(j-1)+0.5*k1))*h;
        k3=(y1(j-1)+0.5*l2)*h;
        l3=(-(u2(j-1)+0.5*(u2(j)-u2(j-1)))-2*c*wn*(y1(j-1)+0.5*l2)-wn*wn*(y(j-1)+0.5*k2))*h;
        k4=(y1(j-1)+l3)*h;
        l4=(-u2(j)-2*c*wn*(y1(j-1)+0.5*l3)-wn*wn*(y(j-1)+k3))*h;
        y(j)=y(j-1)+1.0/6.0*(k1+2*k2+2*k3+k4);
        y1(j)=y1(j-1)+(l1+2*l2+2*l3+l4)/6;
        y2(j)=-u2(j)-2*c*wn*y1(j)-wn*wn*y(j);%相对加速度
        d(j)=y2(j)+u2(j);%绝对jiasudu
        %         y222(j+1)=abs(y22(j+1));
        %         y21(j+1)=y222(j+1)/9.81;%以g为单位
        %         if (abs(d)>abs(dmax))
        %             dmax=d;
        %             gmax=dmax/9.81;
        %         end
    end
    %     if (i>wn)
    %         i=i+1;
    %         break;
    %     end
    i=i+1;
    fzr=round(fz);
    y2max(fzr)=(max(abs(d)))/9.81;
    k(fzr)=fz;
    fz=wn/(2*pi);
    if (fz>3000)
        break;
    end
    %     fz=fz*1.122;
    wn=wn*1.122;
   
end
figure(2);
loglog(y2max,'-*')%对数坐标
grid on

seahero198432

其实响应归根结底还是个波形，就是先采集好产品受到的特定冲击波形的激励后得到的响应，然后再计算响应波形的SRS。说某产品的SRS是没有错的，它是与产品的自身特性有关的。

seahero198432

冲击响应谱里面是Hz-G， 而振动控制是Hz-G2/Hz，是不一样的。

hcharlie

基础是绝对刚体，质量是无穷大。