帮我看看随机振动例子

c_y_dragon

题目是：一个包装产品受到加速度谱为S(w)=0.2/(4*pi)*g^2/HZ的汽车激励．产品连同隔振衬垫的固有频率为２５HZ，阻尼比取０.１，求位移解响应的均方根值．根据公式算得的结果是0.0025m.　ANSYS算得的结果是0.78566e-3m，　相差很大，检查了好多次，找不到原因，请帮我看看，下面是原程序：

/PREP7  
ET,1,MASS21
ET,2,COMBIN14   
KEYOPT,1,1,0
KEYOPT,1,2,0
KEYOPT,1,3,3
KEYOPT,2,1,0
KEYOPT,2,2,0
KEYOPT,2,3,2
R,1,.5, ,  
R,2,1.2325e4, , ,   
N, ,,,,,,,  
N, ,,1,,,,,
E,       2  
TYPE,   2   
REAL,       2   
E,1,2  
FINISH  
/SOL
ANTYPE,2
MODOPT,LANB,1   
EQSLV,SPAR  
MXPAND,1, , ,0  
LUMPM,0
PSTRES,0
MODOPT,LANB,1,0,0, ,OFF
D,1,all
d,2,ux
SOLVE   
FINISH  
/SOL
ANTYPE,8
SPOPT,PSD,1,0   
D,1,UY,1.0   
DMPRAT,0.1,
PSDUNIT,1,ACCG,9.8
PSDFRQ,1, ,1e-2,5e2, , ,   
PSDVAL,1,0.0199,0.0199, , ,
PFACT,1,BASE,   
PSDRES,DISP,REL
PSDRES,VELO,OFF
PSDRES,ACEL,ABS
PSDCOM,0,1,
! /STATUS,SOLU
SOLVE   
FINISH  
/POST1  
set,3,1
PRNSOL,U,Y

非常感谢！

trainhead

PSDCOM,0,1,
你仅考虑1个模态的合并，肯定效果很差
你多扩展几阶，如30阶
QQ：745300845。
我也在做结构方面的psd分析
你的加速度激励的单位是g^2/H
确认一下，单位是否闭合。

c_y_dragon

首先，谢谢
１，单自由度系统只有一个振型吧，扩展那么多有用吗？
２，PSDUNIT,1,ACCG,9.8，还有就是想问问大家，取国际单位制，这条命令是对的吗？

c_y_dragon

这个程序的毛病找到了．还是单位出了问题．w和f之间的转换有个2*pi关系．另外，请大家还有单边谱和双边谱的问题，ansys中输入的应该是单边谱吧？用这个结果正好吻合．

suancaiyu

我用你的程序算的结果怎么是0.024409，是你给出的理论解写错了吗？
我也做谱分析 当时弄了好久才搞懂单位问题，PSDUNIT,1,ACCG,9.8我也是这么写的

c_y_dragon

上面的程序有点问题，谱值应该是０.２．这样结果正好与理论的差不多．
不知道你对阻尼理解的怎么样，这两天用matlab和ansys分别验证了几个两自由度随机振动的例子．发觉把两个不同的粘性阻尼系数转化成
alphad和betad阻尼正好和真实解差不多．由于以前理论学的不好，所以不知道这样做对不？请指教

suancaiyu

请问施加加速度功率谱 怎么得到理论解
我只好做施加力谱后的理论解  可以指教一下吗
QQ 181904318

c_y_dragon

对两自由度系统要想得到理论解是比较困难的，在白噪声情况下，才可以用积分得到．只要是讲随即振动的书上都有那几个公式．
d1=4*w1*w2*(u*w1*w2*(c1*w2+c2*w1)^2+c1*c2*(w2^2-(1+u)*w1^2)^2);
d2=4*w1*w2*(4*c1*c2*w1*w2*(w1*w2*(c2^2+c1^2+u*c1^2)+c1*c2*(w2^2+w1^2+u*w1^2)));
d=d1+d2;
s=0.2*9.8^2/4/pi;   

disp1=sqrt(pi*s*w2^2*(2*c2*w2*(u+(w2/w1)^2)+2*c1*w1*((1+u)^2+u*(w2/w1)^2)+...
    +8*c1*c2*(c2*w1*(1+u+(w2/w1)^2)+c1*w2*(1+u+(c2/c1)^2)))/d)

disp2=sqrt(pi*s*(2*c2*w2*w1^2*(u^2+u*(1+u)^2*(w1/w2)^2)+2*c1*w1*w2^2*((1-(1+u)^2*(w1/w2)^2)^2+...
+u*(1+u)^2*(w1/w2)^2)+8*c2*c1^2*w2*w1^2*(1+u)^2*(1+(1+u)*(w1/w2)^2)+...
    +8*c1^3*w1^3*(1+u)^2*((1+u)+(c2/c1)^2))/d)

ace1=sqrt(2*pi*w1^2*w2^2*s*(c2*w2*(w2^2+u*w1^2)+c1*w1*(u*w2^2+...
+(1+u)^2*w1^2)+4*(c2^3*w2*w1^2+c2^3*c1*(w1*w2^2+(1+u)*w1^3)+...
+c1^2*c2*(w2^3+(1+u)*w1^2*w2)+c1^3*(1+u)*w1*w2^2)+...
+16*c1^2*c2^2*w1*w2*(c1*w2+c2*w1))/d)

ace2=sqrt(2*pi*w1*w2^2*s*(u*c2*w1^3*w2+c1*((w2^2-(1+u)*w1^2)^2+...
+u*w1^2*w2^2)+4*(c2^3*u*w1^3*w2+c1*c2^2*((w1^2-w2^2)^2+...
+w1^2*w2^2+u*w1^4)+c1^2*c2*(w1*w2^3+(1+u)*w1^3*w2)+...
+c1^3*(1+u)*w1^2*w2^2)+16*w1*w2*c1*c2^2*((c1^2+c2^2)*w1*w2+...
+c1*c2*(w1^2+c2^2)))/d)
这是用matlab写的程序，你可以直接假设几个数据试验一下．写这几个公式很费劲，既然写好了，大家就可以拿来用．

suancaiyu

谢谢了

xiaoling850105

一楼计算出的只是位移的响应，为什么会是位移方差的响应？请问哪个程序指的是方差的计算？

余小歪

c_y_dragon 发表于 2007-6-20 16:29
首先，谢谢
１，单自由度系统只有一个振型吧，扩展那么多有用吗？
２，PSDUNIT,1,ACCG,9.8，还有就是想问 ...

9.8是采用国际标准单位下的值