|
楼主 |
发表于 2012-1-11 18:27
|
显示全部楼层
我贴上我的程序 请高手指点。。。实在是搞不懂为什么ansys结果和理论结果差距这么大???!。。。
我用ansys psd 分析了一个单自由度结构(用弹簧质点模型combin40做单元) 在自由度处输入白噪声谱值为1
ansys的分析程序如下:
finish
/clear
/prep7
/title,vm68,psd response of one dof spring-mass system
et,1,combin40 !displacement along X axis,mass at node I
R,1,4.05e5,,1e6 ! 刚度4.05e5,集中质量1e6
n,1
n,2,1
e,2,1
d,1,ux,0.0
outpr,all,all !controls the solution printout
finish
/solu
antype,modal
modopt,subsp,2,,,,off, !提取阶数,模态范围,,归一控制
mxpand,2,,,yes !扩展模态,计算单元应力
solve
/post1
set,list !!!这个频率结果与手算的sqrt(K/M)/2/pi 是一致的!
finish
/SOLu
ANTYPE,2 !分析类型:模态分析
MODOPT,LANB,10,,,,0
MXPAND,10,,,0
SOLVE
FINISH
/POST1
SET,LIST
/SOLU
ANTYPE,8 !8--spectrum分析
SPOPT,PSD,10,yes
DMPRAT,0.02
PSDUNIT,1,forc !定义功率谱类型
psdfrq,1,,0.0001,10
psdval,1,1,1
f,2,fx,1
eplot
pfact,1,node !计算PSD参与系数
psdres,disp,abs !设置文件输出格式为位移解,解为绝对值
PSDCOM,1E-6,10 !模态合并选项,选择合并前n阶
SOLVE
SAVE
FINISH
/post26
store,psd
nsol,2,2,u,x !得到节点2的位移-频率曲线
rpsd,3,2,,,1 !得到节点2的位移功率谱-频率曲线
!!!以下是对vari3进行积分得到均方值
*dim,parm,varry,1,1
int1,4,3,1 !对vari3进行时间轴vari1的积分 保存到vari4
*get,parm(1),vari,4,extrem,vmax !获得变量4的极值
finish
输出结果:节点2的位移-频率曲线 数量级是 10^-3
节点2的得到节点2的位移功率谱-频率曲线 数量级是10^-9
均方值是2.439*10^-11
这个结果与理论解差距太大!!!?
理论上节点2的位移功率谱-频率曲线 应该是传递函数*输入力谱
将得到的输出功率谱曲线进行频域内的积分 值的数量级是10^2... 这个到底问题出在哪里啊!!! 请求高手指点啊 不胜感激!!!
|
|