[求助]直接声振耦合出错

zhynpu

PR VERSION 11-SL1
Requested task is SYSNOISE sequential computation.
Working directory on server is C:\Users\zhy\AppData\Local\Temp
Checking that the SYSNOISE launch script (D:\lms\VirtualLab\Sysnoise\5.6\bin\sysnoise_nogui.bat) exists on server...
OK
SYSNOISE log file on server is C:\Users\zhy\AppData\Local\Temp\zhy-599277-sysnoise.log.
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxx  JOB EXECUTION START
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx

                System for Computational Vibro-Acoustics
------------------------------------------------------------------------
SYSNOISE Rev 5.6  Windows_x86_64  MAR2003 DATE: 20-OCT-2016 19:55:57
PATCH LEVEL  O  (Build Date: 15-NOV-2012)
READING THE SET-UP FILE: D:\lms\VirtualLab\Sysnoise\5.6\bin\\..\res\SYSNOISE.STP
THE SYSNOISE PROFILE FILE NAME HAS NOT BEEN DEFINED.
STANDARD GRAPHIC DEVICE SET TO DEFAULT                 
PLOTTING GRAPHIC DEVICE SET TO SCREENDUMP              

MEMORY IN USE   =         800000
AVAILABLE MEMORY=     2033500000
FREE MEMORY     =     2032700000
CREATING DEFAULT SYSNOISE MODEL FILE default.sdb
SYSNOISE ANALYSIS OPTION :
      BEM INDIRECT VARIATIONAL UNBAFFLED FREQUENCY FLUID

SYSNOISE>  SCL> ENVIRONMENT SECTION SETUP USRDIR 'C:\Users\zhy\AppData\Local\Temp\' RETURN
SYSNOISE>  SCL> ENVIRONMENT SECTION SETUP TMPDIR 'C:\Users\zhy\AppData\Local\Temp\' RETURN
SYSNOISE>  SCL> Open Model 1 File zhy-599277-Structure.sdb Original Return
OPENING SYSNOISE MODEL FILE zhy-599277-Structure.sdb
LOADING MODEL Exported VLModel
WARNING:  There are     53931 3-Dimensional Elements
SYSNOISE>  SCL> Open Model 2 File zhy-599277-Acoustic.sdb Original Return
OPENING SYSNOISE MODEL FILE zhy-599277-Acoustic.sdb
LOADING MODEL Exported VLModel
SYSNOISE>  SCL> Extract Summary Return
SYSNOISE ANALYSIS OPTION :
      FEM VARIATIONAL FREQUENCY FLUID
SYSNOISE SUMMARY INFORMATION
------------------------------------------------------------------------
Title...Exported VLModel
Date....20-OCT-2016 19:55:57
Number of dimensions......................      3
Number of nodes...........................  50280
Number of elements........................ 110185
Number of material characteristics........      2
Number of load cases......................      1
------------------------------------------------------------------------
SYSNOISE>  SCL> Environment Section SETUP GEO_TOLERANCE '1E-006' Return
SYSNOISE>  SCL> Return Environment Section SETUP BELL 'on' Return
SYSNOISE>  SCL> Environment Section SETUP SOLVITER 0 Return
SYSNOISE>  SCL> Environment Section SETUP STABLEFEM TRUE Return
SYSNOISE>  SCL> Parameter Model 1
PARAMETER>  SCL>   Save Displacements None
PARAMETER>  SCL>   Return
Updating Model Parameters
SYSNOISE>  SCL> Parameter Model 2
PARAMETER>  SCL>   Physical
PARAMETER>  SCL>   Save Potentials Step 1
PARAMETER>  SCL>   Save Results none
PARAMETER>  SCL>   Store Potentials none
STORE POTENTIALS>  SCL>   Return
PARAMETER>  SCL>   Store Results none
STORE RESULTS>  SCL>   Return
PARAMETER>  SCL>   NoFlow
PARAMETER>  SCL>   Return
Updating Model Parameters
SYSNOISE>  SCL> Activate Model 2 return
SYSNOISE>  SCL> Solve
SOLVE>  SCL> Frequency 80 To 350 LinStep 1
SOLVE>  SCL> Return
========================================================================
FREQUENCY RESPONSE MODULE
------------------------------------------------------------------------
FOR  271 FREQUENCIES FROM       80.000 TO      350.000 HZ
ANALYSIS #     1  FREQUENCY =         80.000 HZ  TYPE = NORMAL
ANALYSIS #     2  FREQUENCY =         81.000 HZ  TYPE = NORMAL
ANALYSIS #     3  FREQUENCY =         82.000 HZ  TYPE = NORMAL
ANALYSIS #     4  FREQUENCY =         83.000 HZ  TYPE = NORMAL
ANALYSIS #     5  FREQUENCY =         84.000 HZ  TYPE = NORMAL
ANALYSIS #     6  FREQUENCY =         85.000 HZ  TYPE = NORMAL
ANALYSIS #     7  FREQUENCY =         86.000 HZ  TYPE = NORMAL
ANALYSIS #     8  FREQUENCY =         87.000 HZ  TYPE = NORMAL
ANALYSIS #     9  FREQUENCY =         88.000 HZ  TYPE = NORMAL
ANALYSIS #    10  FREQUENCY =         89.000 HZ  TYPE = NORMAL
ANALYSIS #    11  FREQUENCY =         90.000 HZ  TYPE = NORMAL
ANALYSIS #    12  FREQUENCY =         91.000 HZ  TYPE = NORMAL
ANALYSIS #    13  FREQUENCY =         92.000 HZ  TYPE = NORMAL
ANALYSIS #    14  FREQUENCY =         93.000 HZ  TYPE = NORMAL
ANALYSIS #    15  FREQUENCY =         94.000 HZ  TYPE = NORMAL
ANALYSIS #    16  FREQUENCY =         95.000 HZ  TYPE = NORMAL
ANALYSIS #    17  FREQUENCY =         96.000 HZ  TYPE = NORMAL
ANALYSIS #    18  FREQUENCY =         97.000 HZ  TYPE = NORMAL
ANALYSIS #    19  FREQUENCY =         98.000 HZ  TYPE = NORMAL
ANALYSIS #    20  FREQUENCY =         99.000 HZ  TYPE = NORMAL
ANALYSIS #    21  FREQUENCY =        100.000 HZ  TYPE = NORMAL
ANALYSIS #    22  FREQUENCY =        101.000 HZ  TYPE = NORMAL
ANALYSIS #    23  FREQUENCY =        102.000 HZ  TYPE = NORMAL
ANALYSIS #    24  FREQUENCY =        103.000 HZ  TYPE = NORMAL
。。。
ANALYSIS #   264  FREQUENCY =        343.000 HZ  TYPE = NORMAL
ANALYSIS #   265  FREQUENCY =        344.000 HZ  TYPE = NORMAL
ANALYSIS #   266  FREQUENCY =        345.000 HZ  TYPE = NORMAL
ANALYSIS #   267  FREQUENCY =        346.000 HZ  TYPE = NORMAL
ANALYSIS #   268  FREQUENCY =        347.000 HZ  TYPE = NORMAL
ANALYSIS #   269  FREQUENCY =        348.000 HZ  TYPE = NORMAL
ANALYSIS #   270  FREQUENCY =        349.000 HZ  TYPE = NORMAL
ANALYSIS #   271  FREQUENCY =        350.000 HZ  TYPE = NORMAL
------------------------------------------------------------------------
SOLVE PARAMETERS
MODEL:   1 - Exported VLModel
------------------------------------------------------------------------
    PHYSICAL Coordinates
    No saving of DISPLACEMENTS
    Energy Calculation Enabled
    Surface post-processing on positive side
------------------------------------------------------------------------
SOLVE PARAMETERS
MODEL:   2 - Exported VLModel
------------------------------------------------------------------------
    PHYSICAL Coordinates
    Saving POTENTIALS every        1 steps
    No saving of RESULTS
    Energy Calculation Enabled
    Surface post-processing on positive side
------------------------------------------------------------------------
LAUNCHING COUPLED FEM-FEM (Physical Coords - Freq) ANALYSIS

Number of coupled load cases =     1
NO SUPERIMPOSED NODES WITH TOLERANCE=  3.04350E-06
NO ISOLATED NODES
CHECKING ELEMENT TOPOLOGY AND DISTORTION
NO SUPERIMPOSED ELEMENTS WITH TOLERANCE=  3.04350E-06
Checking MESH Connectivity...
     17 Part(s) and    17 Sub-Domain(s) found
      1966 ENVELOPE FACE(S) REMAINING FOR PART    1
      2444 ENVELOPE FACE(S) REMAINING FOR PART    2
      2252 ENVELOPE FACE(S) REMAINING FOR PART    3
      2264 ENVELOPE FACE(S) REMAINING FOR PART    4
      1962 ENVELOPE FACE(S) REMAINING FOR PART    5
      3352 ENVELOPE FACE(S) REMAINING FOR PART    6
      1990 ENVELOPE FACE(S) REMAINING FOR PART    7
      2446 ENVELOPE FACE(S) REMAINING FOR PART    8
      4716 ENVELOPE FACE(S) REMAINING FOR PART    9
      3296 ENVELOPE FACE(S) REMAINING FOR PART   10
      1986 ENVELOPE FACE(S) REMAINING FOR PART   11
      3620 ENVELOPE FACE(S) REMAINING FOR PART   12
      2250 ENVELOPE FACE(S) REMAINING FOR PART   13
      4176 ENVELOPE FACE(S) REMAINING FOR PART   14
      4748 ENVELOPE FACE(S) REMAINING FOR PART   15
      2246 ENVELOPE FACE(S) REMAINING FOR PART   16
      4554 ENVELOPE FACE(S) REMAINING FOR PART   17
Occupied space:             62 Mb  (  0.77 %)
COMPUTING PART INFORMATION:
      Part       1   Volume:   5.18715E-02
      Part       2   Volume:   7.50870E-02
      Part       3   Volume:   6.24870E-02
      Part       4   Volume:   6.24870E-02
      Part       5   Volume:   5.18715E-02
      Part       6   Volume:   1.04256E-01
      Part       7   Volume:   5.18715E-02
      Part       8   Volume:   7.50870E-02
      Part       9   Volume:   5.47470E-02
      Part      10   Volume:   1.04256E-01
      Part      11   Volume:   5.18715E-02
      Part      12   Volume:   3.63720E-02
      Part      13   Volume:   6.24870E-02
      Part      14   Volume:   1.25487E-01
      Part      15   Volume:   5.47470E-02
      Part      16   Volume:   6.24870E-02
      Part      17   Volume:   9.64788E+00
THIS MODEL SEEMS SUITABLE FOR A FEM FLUID APPROACH
CPU Stats --- Delta:  0:00:02 (      1.96)  Total:  0:00:02 (      1.96)

Building FLUID / STRUCTURE interface

Using Model      1 as STRUCTURAL Model
Using Model      2 as FLUID Model
There are  27445 interface node(s) on the FEM FLUID side
There are 241333 interface node(s) on the FEM STRUCTURE side
Generating INTERFACE mesh
Structural INTERFACE MESH contains      81088 elements
Fluid INTERFACE MESH contains      50268 elements
Computing FLUID Element System Matrices
Computing MASS and STIFFNESS elements matrices
using HIGH PRECISION STABILIZED scheme for linear elements
CPU Stats --- Delta:  0:00:01 (      0.90)  Total:  0:00:03 (      2.86)
Computing VELOCITY matrices on elements faces
CPU Stats --- Delta:  0:00:00 (      0.31)  Total:  0:00:03 (      3.17)
CPU Stats --- Delta:  0:00:00 (      0.00)  Total:  0:00:03 (      3.17)
CPU Stats --- Delta:  0:00:00 (      0.00)  Total:  0:00:03 (      3.17)
  Running Bacon and Dynam ...
CPU Stats --- Delta:  0:00:36 (     35.65)  Total:  0:00:39 (     38.82)
Reverse Cuthill McKee Ordering :
    Old maximum bandwidth:        21990
    New maximum bandwidth:         1472
Occupied space:            219 Mb  (  2.70 %)
Computing STRUCTURAL Element System Matrices
Reading Graph File ...
CPU Stats --- Delta:  0:00:02 (      1.52)  Total:  0:00:40 (     40.34)
Occupied space:            244 Mb  (  3.00 %)
Computing structural pressure field element matrices
CPU Stats --- Delta:  0:00:00 (      0.07)  Total:  0:00:40 (     40.41)
Occupied space:            244 Mb  (  3.01 %)
Computing COUPLING matrices for FLUID-STRUCTURE interface elements
****Interpolation between non-matching meshes:****
=== Coupling relation #           1
*The structural interface contained        81088
  faces or elements (one or two-sided).       48584
  acoustic faces have been coupled to the structure.
**WARNING**:         6303 structural
elements, faces or sides of elements where excluded
  from the interface for being beyond the tolerance.
*There were     94 orphaned  nodes out of the  27445 nodes of the interpolating
surface.
*There were      0 points interpolated with a nodal distance-based algorithm and 662875 points interpolated with an element-based algorithm.
*There were            0  structural faces with TRIM coupled to            0
  acoustic faces.
=================
****End of interpolation information:****

CPU Stats --- Delta:  0:10:40 (    640.48)  Total:  0:11:21 (    680.89)
   Compute sparse matrix coordinate profile
   Keeping "L" part of graph
    Memory allocated for the sparse matrix
       Integer*4  :   62139379
       Complex*16 :   61108890
Occupied space:           1504 Mb  ( 18.49 %)
This node (zhy) has number 0
ANALYSIS NUMBER      1 /    271                FREQUENCY =     80.000 Hz
------------------------------------------------------------------------
CPU Stats --- Delta:  0:02:19 (    138.89)  Total:  0:13:40 (    819.78)
Assembling Structural Right-Hand-Side vector
Assembling structural pressure field RHS
A FATAL ERROR OCCURRED.  PROGRAM STOPPED
**** ERROR has occurred during computation ****
You may try to compute the solution in manual mode and try to attach the results afterwards.

zhynpu

是警告引起的错误吗？
会不会是内存不够引起的？

zhynpu

superxjw

我能看到的就是Warning那一段：
**WARNING**:         6303 structural
elements, faces or sides of elements where excluded
  from the interface for being beyond the tolerance.
*There were     94 orphaned  nodes out of the  27445 nodes of the interpolating
surface.
*There were      0 points interpolated with a nodal distance-based algorithm and 662875 points interpolated with an element-based algorithm.
*There were            0  structural faces with TRIM coupled to            0
  acoustic faces.

貌似你定义的耦合面有点问题

lengxuef

这个看不出原因来，仔细检查一下模型吧，节点，单元。看看有没有问题。

zhynpu

我觉得网格应该不会有问题。我在acoustic 网格中加了一个声源后，将直接声振耦合计算原来的use an existing load vector 的改为no load vector set和use existing boundry 算就可以算。但是load vector 的数据是由fluent 导出cgns 变换出的，用简化的结构网格和声学网格做直接声振耦合计算时也可以用来算算。

zhynpu

lengxuef 发表于 2016-10-21 23:31
这个看不出原因来，仔细检查一下模型吧，节点，单元。看看有没有问题。
...

我觉得网格应该不会有问题。我在acoustic 网格中加了一个声源后，将直接声振耦合计算原来的use an existing load vector 的改为no load vector set和use existing boundry 算就可以算。但是load vector 的数据是由fluent 导出cgns 变换出的，用简化的结构网格和声学网格做直接声振耦合计算时也可以用来算算

lengxuef

zhynpu 发表于 2016-10-22 20:35
我觉得网格应该不会有问题。我在acoustic 网格中加了一个声源后，将直接声振耦合计算原来的use an exist ...

问题解决了吗？

zhynpu

lengxuef 发表于 2016-10-27 10:32
问题解决了吗？

没有，不知道该怎么办？

zhynpu

lengxuef 发表于 2016-10-27 10:32
问题解决了吗？

*_samcef.log or *_samcef_dy.res file for errors !!!
重新画的网格报的错

zhynpu

zhynpu 发表于 2016-10-27 19:31
*_samcef.log or *_samcef_dy.res file for errors !!!
重新画的网格报的错

我找到网格新画的网格的问题了，但是还是和之间一样的

superxjw

zhynpu 发表于 2016-10-27 21:49
我找到网格新画的网格的问题了，但是还是和之间一样的

我看了你发的东西，你为什么要把CFD计算结果Data Transfer到结构网格上呢？我印象中做耦合，是先转移到声学网格上吧，你看看是不是这一步做错了。

zhynpu

superxjw 发表于 2016-10-30 10:35
我看了你发的东西，你为什么要把CFD计算结果Data Transfer到结构网格上呢？我印象中做耦合，是先转移到声 ...

你的书，第388页，也是这样做的，不过用的是模态声振耦合。我用模态声振耦合也可以算，但是我算机械模态时number of modes指定为100，运算时不收敛终止了，用frequency range 0~100Hz模式算的模态。不知道这样的行不行？我对机械方面实在不了解。

zhynpu

superxjw 发表于 2016-10-30 10:35
我看了你发的东西，你为什么要把CFD计算结果Data Transfer到结构网格上呢？我印象中做耦合，是先转移到声 ...

算模态的时候，没有restrain。

superxjw

zhynpu 发表于 2016-10-31 20:42
你的书，第388页，也是这样做的，不过用的是模态声振耦合。我用模态声振耦合也可以算，但是我算机械模态 ...

对结构没有约束，就表示自由模态。有一个很关键的问题，就是模态截止频率。比如你最后要计算到150Hz的响应，那么一般来讲，模态至少要计算到300Hz，才能保证一定精度。