Virtual.lad or SYSNOISE 计算声—结构灵敏度

zhongsheng

各位朋友，大家好，这段时间想计算声—结构灵敏度，看了很多前辈的思路是用nastran+sysnoise联合使用，我在摸索的时候，发现了少许资料，提供的是在nastran中计算的部分，但现在我却没法通过nastran得到中间计算结果。请教下各位：

（1）       针对该问题，nastran和sysnoise分别该如何使用；

（2）  现在Virtual.lab很发达了，有没有前辈在该软件中实现声—结构灵敏度的计算呢？

superxjw

呵呵！是的，LMS SYSNOISE早都停止开发了，升级为LMS Virtual.Lab都已经快9年了。灵敏度具体的东西我没有做过，如果论坛里有朋友做过，欢迎来讨论。
关于灵敏度计算，你可以先看ON-LINE HELP里面的这几个算例。

顺便说一点，现在VL里面，可以直接调用Nastran了，不需要你去掌握Nastran的卡片，就直接能给边界条件，提交Nastran计算，而且也不用再通过手动导入Nastran结果的方法进行分析，直接从结构树就可以选择计算结果。

至于你说的直接计算，没有考察过，现在VL11里面确实有结构求解器的！但是具体没有操作过！有知道的朋友可以来讨论一下哈！

zhongsheng

superxjw 发表于 2012-12-25 10:08

                               
登录/注册后可看大图

呵呵！是的，LMS SYSNOISE早都停止开发了，升级为LMS Virtual.Lab都已经快9年了。灵敏度具体的东西我没有做 ...

我下载的资料中，nastran中的卡片是这样的。
NASTRAN SYSTEM(2)=4
ASSIGN OUTPUT2='myoutput2.OP2',UNIT=34,STATUS=UNKNOWN,FORM=UNFORMATTED,DELETE$
TIME 100
SOL 200                  
COMPILE SUBDMAP=SEDRCVR,SOUIN=MSCSOU,NOREF,LIST
ALTER 139
VECPLOT  UG,BGPDTS,SCSTM,CSTMS,CASEDR,,MEDGE,/UGVSB1/
         /0/1/'PARAM,'/'OUGCORD,'/'BASIC'/ALTSHAPE//SEID
SDR2     CASEDR,CSTMS,,,EQEXINS,,,OL2,BGPDTS,,,
         UGVSB1,,,OINT,PELSETS,VIEWTB,,,,/
         ,,OUGV1SYS,,,/APP1
COPY OL2/OL1S
OUTPUT2  OL1S,OUGV1SYS//0/34
COMPILE SUBDMAP=EXITOPT,SOUIN=MSCSOU,NOREF,LIST
ALTER 35
OUTPUT2  DESTAB,DSCMCOL,DSCM2//0/34
OUTPUT2  ,,,,//-9/34
CEND
TITLE=DYNAMIC FREQUENCY SENSITIVITY BOX
VELO(SORT1,punch) = ALL  
DESSUB       = 1
DLOAD        = 700
FREQ         = 740
ANALYSIS     = DFREQ
BEGIN BULK
PARAM,DDRMM,-1
PARAM,OPTEXIT,-4
$
$ Direct Text Input for Bulk Data
$ Elements and Element Properties for region : p1
PSHELL   1       1       .03     1               1
$ Pset: "p1" will be imported as: "pshell.1"
CQUAD4   26      1       6       37      42      12
CQUAD4   27      1       37      38      43      42
CQUAD4   28      1       38      39      44      43
.
.
.
$ Referenced Material Records
$ Material Record : steel
$ Description of Material : Date: 03-Sep-04           Time: 15:03:20
MAT1     1       2.1+10          .3      7850.
$ Nodes of the Entire Model
GRID     1               1.      1.      0.
GRID     2               1.      .8      0.
GRID     3               1.      .6      0.
GRID     4               1.      .4      0.
GRID     5               1.      .2      0.
GRID     6               1.      0.      0.
GRID     7               1.      1.      .2
.
.
.
$ Loads for Load Case : Default
$ Referenced Coordinate Frames
$ ...DESIGN VARIABLE DEFINITION
$ p1_Thickness
DESVAR   1      p1_Thi:1 .03     .027    .033    1.
$ p2_Thickness
DESVAR   2      p2_Thi:2 .03     .027    .033    1.
$ p3_Thickness
DESVAR   3      p3_Thi:3 .03     .027    .033    1.
$ p4_Thickness
DESVAR   4      p4_Thi:4 .03     .027    .033    1.
$ p5_Thickness
DESVAR   5      p5_Thi:5 .03     .027    .033    1.
$ p6_Thickness
DESVAR   6      p6_Thi:6 .03     .027    .033    1.
$ ...DEFINITION OF DESIGN VARIABLE TO ANALYSIS MODEL PARAMETER RELATIONS
DVPREL1  1       PSHELL  1       T
         1       1.
DVPREL1  2       PSHELL  2       T
         2       1.
DVPREL1  3       PSHELL  3       T
         3       1.
DVPREL1  4       PSHELL  4       T
         4       1.
DVPREL1  5       PSHELL  5       T
         6       1.
DVPREL1  6       PSHELL  6       T
         5       1.
$ Loads for Load Case : frf
RLOAD1   5       3                       1        
RLOAD1   7       4                       1
DLOAD    700     1.      1.       5      1.      7
$ Nodal Forces of Load Set :force
FORCE    3       37      0      30.       0.      1.     0.
$ Nodal Forces of Load Set :force2
FORCE    4       41      0      15.       0.      0.     1.
$ Referenced Dynamic Load Tables
$ Dynamic Load Table :frfrd
$TABLED1 1 1. 1. 100. 1.0，ENDT
TABLED1  1
        -102.2  2.10562  0.     2.10562 2.      1.37719 4.      .739608     
FREQ1   740      10.      1.      40
DRESP1,1       ,SYSNOISE, FRVELO, , ,       1, , 1  
DRESP1,2       ,SYSNOISE, FRVELO, , ,       2, , 1  
DRESP1,3       ,SYSNOISE, FRVELO, , ,       3, , 1  
DRESP1,4       ,SYSNOISE, FRVELO, , ,       7, , 1  
.
.
.
$
$ AND SET LOWER AND UPPER BOUNDS T
$ VALUE SO IT WILL ALWAYS BE RETAI
$
$ DCONSTR,1,   RID,    SID, LALLOW,
$
DCONSTR,1, 1     ,      +1.+35, +1.+35,
DCONSTR,1, 2     ,      +1.+35, +1.+35,
DCONSTR,1, 3     ,      +1.+35, +1.+35,
DCONSTR,1, 4     ,      +1.+35, +1.+35,
DCONSTR,1, 5     ,      +1.+35, +1.+35,
.
.
.
$
ENDDATA


不知道问题出在哪儿，就是结算不出结果。

superxjw

zhongsheng 发表于 2012-12-25 15:23

                               
登录/注册后可看大图

我下载的资料中，nastran中的卡片是这样的。
NASTRAN SYSTEM(2)=4
ASSIGN OUTPUT2='myoutput2.OP2',UNI ...

用新方法吧……