声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2206|回复: 10

[结构分析] 帮我看看随机振动例子

[复制链接]
发表于 2007-6-20 15:39 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
题目是:一个包装产品受到加速度谱为S(w)=0.2/(4*pi)*g^2/HZ的汽车激励.产品连同隔振衬垫的固有频率为25HZ,阻尼比取0.1,求位移解响应的均方根值.根据公式算得的结果是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

非常感谢!
回复
分享到:

使用道具 举报

发表于 2007-6-20 15:47 | 显示全部楼层

原因

PSDCOM,0,1,
你仅考虑1个模态的合并,肯定效果很差
你多扩展几阶,如30阶
QQ:745300845。
我也在做结构方面的psd分析
你的加速度激励的单位是g^2/H
确认一下,单位是否闭合。
 楼主| 发表于 2007-6-20 16:29 | 显示全部楼层
首先,谢谢
1,单自由度系统只有一个振型吧,扩展那么多有用吗?
2,PSDUNIT,1,ACCG,9.8,还有就是想问问大家,取国际单位制,这条命令是对的吗?
 楼主| 发表于 2007-6-21 10:33 | 显示全部楼层
这个程序的毛病找到了.还是单位出了问题.w和f之间的转换有个2*pi关系.另外,请大家还有单边谱和双边谱的问题,ansys中输入的应该是单边谱吧?用这个结果正好吻合.
发表于 2007-6-22 15:54 | 显示全部楼层
我用你的程序算的结果怎么是0.024409,是你给出的理论解写错了吗?
我也做谱分析 当时弄了好久才搞懂单位问题,PSDUNIT,1,ACCG,9.8我也是这么写的
 楼主| 发表于 2007-6-22 17:58 | 显示全部楼层
上面的程序有点问题,谱值应该是0.2.这样结果正好与理论的差不多.
不知道你对阻尼理解的怎么样,这两天用matlab和ansys分别验证了几个两自由度随机振动的例子.发觉把两个不同的粘性阻尼系数转化成
alphad和betad阻尼正好和真实解差不多.由于以前理论学的不好,所以不知道这样做对不?请指教

评分

1

查看全部评分

发表于 2007-6-25 16:02 | 显示全部楼层
请问施加加速度功率谱 怎么得到理论解
我只好做施加力谱后的理论解  可以指教一下吗
QQ 181904318
 楼主| 发表于 2007-6-25 19:47 | 显示全部楼层
对两自由度系统要想得到理论解是比较困难的,在白噪声情况下,才可以用积分得到.只要是讲随即振动的书上都有那几个公式.
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写的程序,你可以直接假设几个数据试验一下.写这几个公式很费劲,既然写好了,大家就可以拿来用.

评分

1

查看全部评分

发表于 2007-6-27 21:51 | 显示全部楼层
谢谢了
发表于 2010-5-15 10:41 | 显示全部楼层
一楼计算出的只是位移的响应,为什么会是位移方差的响应?请问哪个程序指的是方差的计算?
发表于 2014-7-16 16:40 | 显示全部楼层
c_y_dragon 发表于 2007-6-20 16:29
首先,谢谢
1,单自由度系统只有一个振型吧,扩展那么多有用吗?
2,PSDUNIT,1,ACCG,9.8,还有就是想问 ...

9.8是采用国际标准单位下的值
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-10 23:42 , Processed in 0.063185 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表