excite-pu&nastran&virtual.lab振动噪声仿真

一、参考资料

1、nastran动力学培训PPT

2、nastran频响

TC ，TD

3、基于LMS Virtual Lab Acoustics的发动机结构噪声预测

二 、excite-pu生成nastran频响计算文件

1、excite-pu动力学计算

注意，时域数据取一个循环，频域数据取2个循环。

2、生成nastran频响计算文件

3、nastran频响计算文件

$############################################################
$ SIMULATION RESULTS
$ Start angle (time): 720 deg (0.0393 s)
$ End angle (time): 2160 deg (0.1176 s)
$ Number of time steps: 1102
$ Mean speed: 3062 rpm
$ Sampling frequency: 25.5 Hz
$############################################################
$ Excite FEA Dynamic Response Analysis
$ Powerunit_dresp2
$ NASTRAN 2013
$############################################################
$
NASTRAN SYSTEM(1)=65537
NASTRAN SYSTEM(151)=1
NASTRAN Q4TAPER=45.
$
$
$ ------------ FILE MANAGEMENT SECTION ----------------------------------------
$
INIT DBALL LOGICAL=(DB1(100MB),DB2(100MB))
ASSIGN DB1=’.\00_Powerunit_dresp2.DBALL’
ASSIGN DB2=’.\01_Powerunit_dresp2.DBALL’
$
INIT MASTER(S)
$
INIT SCRATCH LOGICAL=(SL1(500MB)),
SCR300=(SC1(10GB))
ASSIGN SL1=’.\00_Powerunit_dresp2.SCRATCH’
ASSIGN SC1=’.\01_Powerunit_dresp2.SCR300’
$
$
$ ----------- EXECUTIVE CONTROL SECTION ---------------------------------------
$
ID NASTRAN,EXCITE
SOL 111
$
CEND
$
$
$
$ ------------ CASE CONTROL CARDS ----------------------------------------------
$
TITLE=Excite FEA Dynamic Response Analysis
SUBTITLE=Job: Powerunit_dresp2
$
METHOD=103
DLOAD=111
FREQUENCY=122
SDAMPING=144
SPC=888 $多个约束的ID
VELOCITY(SORT1,PLOT,PHASE)=ALL
$
$
$ ------------------------------------------------------------------------------
$ B U L K D A T A S T A T E M E N T S
$ ------------------------------------------------------------------------------
BEGIN BULK
PARAM,DDRMM,-1
$
$ PARAM, WTMASS, 1E-3
PARAM,POST,-1
$
$ Include FE Model File:
INCLUDE ‘…/…/…/fem/Powerunit_nas/Powerunit.bdf’
$
INCLUDE ‘…/…/…/fem/Powerunit_nas/bulk.inc’
$
$
EIGRL,103,0.,3000.,
$############################################################
$ Number of frequency steps: 117
$############################################################
FREQ,122,0.00,25.54,51.09,76.63,102.18,127.72,153.27,
。。。。。。
$ Modal Damping:
TABDMP1,144,CRIT
,10.00,0.17,20.00,0.09,30.00,0.06,40.00,0.05,
,50.00,0.04,60.00,0.03,70.00,0.03,80.00,0.03,
。。。。。。
,ENDT
$
DLOAD,111,1.0,1.0,1001,1.0,1002,1.0,1003,
,1.0,1004,1.0,1005,1.0,1006,1.0,1007,
。。。。。。
$ Number of frequency steps: 117
$
DAREA,1001,3376,3,1.0 $3376节点z向，比例系数1.0
RLOAD1,1001,1001,1,2
TABLED1* 1
*

•   0.00000000E+00  9.84588520E+02  2.55449599E+01  1.59271474E+03  
  

•   5.10899198E+01  1.34930993E+03  7.66348797E+01  1.08823536E+03  
  

• 。。。。。。
  $以上是载荷的实部表，以下是虚部表

•   2.96321535E+03  3.84938442E-01              ENDT
  

TABLED1* 2
*

•   0.00000000E+00  0.00000000E+00  2.55449599E+01  -2.39799310E+02 
  

•   5.10899198E+01  -3.28499694E+02 7.66348797E+01  -3.67624485E+02 
  

。。。。。。

•   2.96321535E+03  6.18266694E-01              ENDT
  

$

。。。。。。
。。。。。。

4、nastran频响计算文件激励力校核

取前左悬置点x向力最后一个循环做FFT，不加窗的情况和result中最后两个循环加海宁窗的结果基本一致，但是result中最后两个循环加海宁窗的结果频率是0时幅值为0，应该与实际不符合，nastran频响计算文件激励力频率为0时幅值不为0.

三、nastran及lms virtual.lab的详细过程

一、nastran结构频响仿真
1、Hypermesh建模

2、加载频域激励

3、bdf文件
$
$ NASTRAN Input Deck Generated by HyperMesh Version :14.0.0.88
$ Generated using HyperMesh-Nastran Template Version : 14.0
$ $
$ Template: NastranNX general $
$ $

$ Executive Control Cards $
SOL 111
CEND

$$ Case Control Cards $

TITLE = Nastran Dynaic Response Anaysis
SUBTITLE = Job:Test
LABEL = Test
$
$HMNAME LOADSTEP 1"SubCase1"
SUBCASE 1
LABEL= SubCase1
METHOD(STRUCTURE) = 103
FREQUENCY = 122
SDAMPING(STRUCTURE) = 144
DLOAD = 5
ANALYSIS = MFREQ
DISPLACEMENT(SORT1,PLOT) = ALL

$$ Bulk Data Cards $

BEGIN BULK
PARAM,DDRMM,-1
PARAM,POST,-1
$
$ GRID Data
$
GRID 2 50.0 -50.0 0.0
GRID 3 50.0 -40.0 0.0
GRID 4 50.0 -30.0 0.0
。。。。。。
$$

$$ Group Definitions $

$
$ CHEXA Elements: First Order
$
CHEXA 101 1 56 59 60 55 123 124+

•        125     126
  

CHEXA 102 1 55 60 51 52 126 125+

•        127     128
  

。。。。。。
PSOLID 1 1
MAT1 1212000.0 0.3 7.8-9
$$

$$ Loads and Boundary Conditions $

$
$HyperMesh name and color information for generic loadcollectors
$
$HMNAME LOADCOL 1"Darea1"
$HWCOLOR LOADCOL 1 3
$
$
$ FREQi cards
$
$HMNAME LOADCOL 122"Freqi"
$HWCOLOR LOADCOL 122 17
FREQ 122 0.0 10.0 2000.0 4000.0 5000.0
$
$
$ RLOAD1 cards
$
$HMNAME LOADCOL 4"Rload1"
$HWCOLOR LOADCOL 4 11
RLOAD1 4 1 2 3
$
$
$ DLOAD cards
$
$HMNAME LOADCOL 5"Dload"
$HWCOLOR LOADCOL 5 11
DLOAD 5 1.0 1.0 4
$
$ EIGRL cards
$
$HMNAME LOADCOL 103"Eigrl"
$HWCOLOR LOADCOL 103 6
EIGRL 103 20 MASS
$
$ DAREA Data
$
DAREA 1 131 3 1.0
$
$ TABDMP1 cards
$
$HMNAME LOADCOL 144"TabDMP1"
$HWCOLOR LOADCOL 144 7
$
TABDMP1 144 CRIT

•       10.0    0.17  1000.0    0.18  3000.0     0.2ENDT    
  

$
$
$ TABLED1 cards
$
$HMNAME LOADCOL 2"Tabeld1_R"
$HWCOLOR LOADCOL 2 11
TABLED1 2 LINEAR LINEAR

•        0.0     0.0  3999.0     0.0  4000.0100000.0  4001.0     0.0
  

•     5000.0     0.0ENDT    
  

$$
$HMNAME LOADCOL 3"Tabeld1_I"
$HWCOLOR LOADCOL 3 11
TABLED1 3 LINEAR LINEAR

•        0.0     0.0  3999.0     0.0  4000.0     0.0  4001.0     0.0
  

•     5000.0     0.0ENDT    
  

$$
ENDDATA

二、LMS Virtual.Lab声学仿真
1、 Import nastran结果op2文件，导入结构网格和load vector;
2、 生成结构网络的包络面Skin Mesh,并设置成声学网络类型；如果从文件中导入声学网络，可能需要进行网格节点和单元编号冲突检查：Tools->ID and Model Checks->Check ID Conflicts;

3、 进行声学网络前处理，定义对称面模拟半消声室的地面，定义流体材料和属性，定义场点；

4、 新建结构网络外表面的单元节点组，用于将Load Vector从结构网络转移到声学网络；

5、 数据映射转移Insert->Vibro->Acoustic Mapping Set;

6、 新建**声学网络**的网格节点组，用于制作面板；

7、 制作面板

8、 定义边界条件；

9、 新建声学响应分析CASE

10、 在面板贡献结果中查看声功率的板块贡献量；

11、 在面板贡献结果中查看场点19 在4000HZ的声压板块贡献量；

12、 在声学响应结果中查看场点19的声压级频率响应曲线。