在NASTRAN中粘弹性阻尼结构的分析方法及计算程序
- 格式：doc
- 大小：428.50 KB
- 文档页数：6
1 在NASTRAN中粘弹性阻尼结构的分析方法
吕毅宁
1. 背景
MSC.Nastran 2005 Quick Reference Guide
MSC.Nastran Version 70 Advanced Dynamic Analysis User’s Guide 2
3
2. 粘弹性材料的复刚度矩阵，推导公式，得到NASTRAN计算所需的输入格式
41)1(ddddddKKigK VddREFVddKgK14
VddREFREFVddVddVddKfTIggifTRgKfiTIfTRKigfK141)()(1)()()1()(
)()()(fGifGfG VddREFVddKGfGifGfK1)()()(
1)(1)(REFREFGfGgfTR gGfGgfTIREFREF)(1)(
4
5 3. 建模需要的Key words
4. 参数计算Matlab程序
close all; clear all; clc;
g=0.0; g=0.02; %g: Global structural damping;;
Gref=1.0; gref=1.0; %Gref: Reference shear modulus(G); gref: Reference
structural damping(GE);
f=10.0:10:500.0; %Frequency;
f=f';
%Gp=102.39e4.*f.^0.625; eta=1.6274.* f.^0.072; Gpp= Gp.*eta; 6 Gp=102.39e-2.*f.^0.625; eta=1.6274.* f.^0.072; Gpp= Gp.*eta;
TR=1.0/gref.*(Gp./Gref-1.0); %G';
TI=1.0/gref.*(Gpp./Gref-g); %G";
plot(f,TR,'-+'); grid on;
figure; plot(f,eta,'-+'); grid on;
figure; plot(f,TI,'-+'); grid on;
fid=fopen('TRTI.dat','w');
for n=1:length(f)
fprintf(fid,', %7.4f, %7.4f',f(n),TR(n));
if(mod(n,4)==0) fprintf(fid,'\n'); end;
end
fprintf(fid,'\n\n\n')
for n=1:length(f)
fprintf(fid,', %7.4f, %7.4f',f(n),TR(n));
if(mod(n,4)==0) fprintf(fid,'\n'); end;
end
上面标示为蓝色的是NASTRAN输入参数，输入方式如下：
g PARAM, G, 0.02
Gref G
gref GE
TR TABLED1 n
TI TABLED1 n+1
