做形状优化时候报错

wangyuewu

做一个freeshap的优化，check的时候可以过去，计算的过程中迭代几步就报错了，请给位朋友帮看看，非常感谢啊！！！
out文件如下：
  
************************************************************************
**                                                                    **
**                                                                    **
**                          OptiStruct 11.0                           **
**                                                                    **
**                  Design and Optimization Software                  **
**                   from Altair Engineering, Inc.                    **
**                                                                    **
**                                                                    **
**   Windows XP x64 Professional  SP1 (Build 3790)  ZYYJ-YF-YH-WYW    **
**           4 CPU:  Intel(R) Xeon(R) CPU  W3520  @ 2.67GHz           **
**                        CPU speed  2650 MHz                         **
**                     4095 MB RAM,  4095 MB swap                     **
**                                                                    **
** Build tag: 0569254_5333110_Ci32RH11M:61051-000  4000               **
************************************************************************
**  COPYRIGHT (C) 1996-2011                 Altair Engineering, Inc.  **
** All Rights Reserved.  Copyright notice does not imply publication. **
**         Contains trade secrets of Altair Engineering, Inc.         **
** Decompilation or disassembly of this software strictly prohibited. **
************************************************************************

  
*** OptiStruct defaults set from:
   install config file: C:\Program Files (x86)\Altair\11.0/hwsolvers/hwsolver.cfg.
  
1 PARAM(s) were set:
    CHECKEL ,NO
*** WARNING # 1265
PSHELL 2 has thickness 0.001 or less and bending properties
defined. This can lead to matrix singularities, causing message 153.
If this element is intended to be only a membrane element, please leave
MID2 and MID3 blank on the PSHELL data.
   Detected while reading line 480237 from file
     D:/WORK2012/A3/freeshape/00001.fem:
  "PSHELL         2       10.001          1               1        0.0"
*** WARNING 1932:    PMASS    3 is not referenced.
*** INFORMATION #  742
The dependent rotational d.o.f. of this rigid element is removed.
     RBE2 element id =    1774515
independent grid id =     306780
a dependent grid id =     222007
This is because there is no need to constrain the rotational d.o.f. of
any of the dependent grids.
*** INFORMATION #  743
The total number of rigid elements, whose rotational dependent d.o.f.
are removed because there is no need to constrain those d.o.f., is         1
NOTE #  1847
AUTOMATIC SCREENING is activated.
Use "DSCREEN,AUTO,OFF" to turn off this effect.

************************************************************************
OPTIMIZATION FILE AND PARAMETER INFORMATION :
---------------------------------------------
  FEM model file               : D:/WORK2012/A3/freeshape/00001.fem
  Output files prefix          : D:/WORK2012/A3/freeshape/00001
  Model read from files:
    D:/WORK2012/A3/freeshape/00001.fem
************************************************************************

************************************************************************
FINITE ELEMENT MODEL DATA INFORMATION :
---------------------------------------
  Total # of Grids (Structural)        :    89677
  Total # of Elements                  :   350842
  Total # of Rigid Elements            :        1
  Total # of Rigid Element Constraints :     3315
  Total # of Degrees of Freedom        :   278205
  Total # of Non-zero Stiffness Terms  :  5492586

  Element Type Information
  ------------------------
   CTETRA   Elements   :   344714
   CONM2    Elements   :        1     
   CTRIA3   Elements   :     6127     
  Load and Boundary Information
  -----------------------------
   FORCE Sets        :        4
   SPC   Sets        :        1
  Material and Property Information
  ---------------------------------
   PSOLID Cards      :        1
   PMASS  Cards      :        1
   PSHELL Cards      :        1
   MAT1   Cards      :        1
************************************************************************

************************************************************************
OPTIMIZATION PROBLEM PARAMETERS :
---------------------------------
Objective Function  : Minimize Weighted Compliance
------------------
Response Summary :
------------------
Number of compliance responses        :         4
Number of equation responses          :         1
   ------------------------
   Static Subcase Summary :
   ------------------------
   ----------   --------  --------  --------------
   Subcase ID    SPC ID   LOAD ID   Weight = w(i)
   ----------   --------  --------  --------------
          1            5         1           1.000
          2            5         2           1.000
          3            5         3           1.000
          4            5         4           1.000
   ----------   --------  --------  --------------
   ---------------------------
   Design Parameters Summary :
   ---------------------------
   Total # of free-shape cards                     :        1

   ---------------------------------
   Optimization Parameters Summary :
   ---------------------------------
   Maximum Number of Iterations    :       30

   Convergence Tolerance           :   5.0000E-03
   Step Size (Shape)               :   0.2000
   Run Type                        :   Free-Shape Optimization
   Restart from previous solution  :   No
   Scratch file directory          : ./
                                     Free space:  31767 MB
   Number of CPU processors        :         1

************************************************************************

************************************************************************
MEMORY ESTIMATION INFORMATION :
-------------------------------
Solver Type is:  Sparse-Matrix Solver
                  Direct Method

Current Memory (RAM)                                    :     261 MB
Estimated Minimum Memory (RAM) for Out of Core Solution :     261 MB
Recommended Memory (RAM) for Out of Core Solution       :     261 MB
Recommended Memory (RAM) for In-Core Solution           :    1259 MB
Recommended Number of Nodes for OS SPMD Parallel Run    :       1
(Note: The Minimum Memory Requirement is limited by Output Module.)         

DISK SPACE ESTIMATION INFORMATION :
-----------------------------------
Estimated Disk Space for Output Data Files              :    1565 MB
Estimated Scratch Disk Space for In-Core Solution       :     225 MB
Estimated Scratch Disk Space for Out of Core Solution   :    1534 MB
************************************************************************

BEGINNING OPTIMIZATION SOLUTION ....

************************************************************************
OPTIMIZATION HISTORY INFORMATION :
----------------------------------

ITERATION   0
          Auto-SPC d. o. f.s for loadcase        1
          Number of Auto-SPC d. o. f.s  =     3315
          ------------------------------
            Grid No.       Component  
          ------------------------------
              222007             4 5 6
              222008             4 5 6
              222009             4 5 6
              222010             4 5 6
              222011             4 5 6
              222012             4 5 6
              222013             4 5 6
              222014             4 5 6
              222015             4 5 6
              222016             4 5 6
              222017             4 5 6
              222018             4 5 6
              222019             4 5 6
              222020             4 5 6
              222021             4 5 6
              222022             4 5 6
              222023             4 5 6
              222024             4 5 6
              222025             4 5 6
              222026             4 5 6
              222027             4 5 6
              222028             4 5 6
              222029             4 5 6
              222030             4 5 6
              222031             4 5 6
              222032             4 5 6
              222033             4 5 6
              222034             4 5 6
              222035             4 5 6
              222036             4 5 6
              222037             4 5 6
              222038             4 5 6
              222039             4 5 6
              222040             4 5 6
          ------------------------------
            By default, Auto-SPC printing is limited to 100 dofs.

(Scratch disk space usage for starting iteration = 1134 MB)
(Running out-of-core solution)
Objective Function (Minimize WCOMP) =  2.92171E-02
Maximum Constraint Violation %      =  0.00000E+00
Volume                              =  2.41482E-02   Mass     =  1.73868E+02

Subcase  Weight     Compliance                            Weight*Comp.
       1  1.000E+00  2.646300E-02                          2.646300E-02
       2  1.000E+00  1.265212E-03                          1.265212E-03
       3  1.000E+00  1.809356E-07                          1.809356E-07
       4  1.000E+00  1.488681E-03                          1.488681E-03
                                                           ------------
Sum of Weight*Compliance                                   2.921708E-02

                           RETAINED RESPONSES TABLE
------------------------------------------------------------------------------------
Response Type  Response    Subcase Grid/     DOF/   Response     Objective    Viol.
User-ID          Label     /RANDPS Element/  Comp     Value      Reference/     %  
                           +Frqncy MID/PID/  /Reg                Constraint
                                   Mode No.                        Bound   
------------------------------------------------------------------------------------
       1 WCOMP wcomp            --       --    --   2.922E-02 MIN
------------------------------------------------------------------------------------
ITERATION   1
Objective Function (Minimize WCOMP) =  2.88334E-02   % change =        -1.31
Maximum Constraint Violation %      =  0.00000E+00
Volume                              =  2.41740E-02   Mass     =  1.74054E+02

Subcase  Weight     Compliance                            Weight*Comp.
       1  1.000E+00  2.612493E-02                          2.612493E-02
       2  1.000E+00  1.245180E-03                          1.245180E-03
       3  1.000E+00  1.795414E-07                          1.795414E-07
       4  1.000E+00  1.463066E-03                          1.463066E-03
                                                           ------------
Sum of Weight*Compliance                                   2.883336E-02

                           RETAINED RESPONSES TABLE
------------------------------------------------------------------------------------
Response Type  Response    Subcase Grid/     DOF/   Response     Objective    Viol.
User-ID          Label     /RANDPS Element/  Comp     Value      Reference/     %  
                           +Frqncy MID/PID/  /Reg                Constraint
                                   Mode No.                        Bound   
------------------------------------------------------------------------------------
       1 WCOMP wcomp            --       --    --   2.883E-02 MIN
------------------------------------------------------------------------------------
ITERATION   2
Objective Function (Minimize WCOMP) =  2.85516E-02   % change =        -0.98
Maximum Constraint Violation %      =  0.00000E+00
Volume                              =  2.41926E-02   Mass     =  1.74188E+02

Subcase  Weight     Compliance                            Weight*Comp.
       1  1.000E+00  2.587615E-02                          2.587615E-02
       2  1.000E+00  1.231201E-03                          1.231201E-03
       3  1.000E+00  1.790000E-07                          1.790000E-07
       4  1.000E+00  1.444060E-03                          1.444060E-03
                                                           ------------
Sum of Weight*Compliance                                   2.855159E-02

                           RETAINED RESPONSES TABLE
------------------------------------------------------------------------------------
Response Type  Response    Subcase Grid/     DOF/   Response     Objective    Viol.
User-ID          Label     /RANDPS Element/  Comp     Value      Reference/     %  
                           +Frqncy MID/PID/  /Reg                Constraint
                                   Mode No.                        Bound   
------------------------------------------------------------------------------------
       1 WCOMP wcomp            --       --    --   2.855E-02 MIN
------------------------------------------------------------------------------------
ITERATION   3
*** ERROR #   15 ***
CTETRA          element    1391820 is distorted (zero or negative Jacobian).
This error occurs in module "asmbdrv".
************************************************************************
RESOURCE USAGE INFORMATION
--------------------------
MAXIMUM MEMORY USED                                       261 MB
MAXIMUM DISK SPACE USED                                  2699 MB
************************************************************************

************************************************************************
COMPUTE TIME INFORMATION
------------------------
EXECUTION STARTED                         Fri May 25 10:15:42 2012   
EXECUTION COMPLETED                       Fri May 25 10:20:43 2012   
ELAPSED TIME                                         00:05:00
CPU TIME                                             00:04:17
************************************************************************

*** PROGRAM STOPPED: FATAL ERROR(s) ENCOUNTERED.
                        ***** END OF REPORT *****
For Useful OptiStruct Tips and Tricks, go to the URL:
http://www.altairhyperworks.com/tips.aspx
************************************************************************
RADIOSS/OptiStruct error termination report printed to file
"00001.stat".
PROGRAM STOPPED DUE TO ERROR.
************************************************************************
  
A fatal error has occurred during computations:
  
  *** ERROR #   15 ***
  CTETRA          element    1391820 is distorted (zero or negative Jacobian).
  
************************************************************************

zkong

经过几个迭代变形后，网格质量太差。 可以把最终的结果重新划分网格，改进网格质量后重新优化。

zkong

方便的话，把模型贴上来，我可以帮你看看。

lbsys_3122

zkong 发表于 2012-5-25 11:12
经过几个迭代变形后，网格质量太差。 可以把最终的结果重新划分网格，改进网格质量后重新优化。 ...

您好~~~你的意思是把几何导出来，再划分网格么？

zkong

apply result 把优化的结果(新的节点坐标)施加到网格上，就可以得到新的网格。  网格重新划分不需要几何。

wangyuewu

zkong 发表于 2012-5-25 11:12
经过几个迭代变形后，网格质量太差。 可以把最终的结果重新划分网格，改进网格质量后重新优化。 ...

谢谢Zkong版主
您说的这个办法具体如何操作呢？
出错的模型我变量选择的是三个面（一起选择）上的接点，如果三个面单独作为变量计算都没有问题，怎么回事呢？
模型还是不太方便上传啊！非常感谢！

wangyuewu

zkong 发表于 2012-5-27 20:42
apply result 把优化的结果(新的节点坐标)施加到网格上，就可以得到新的网格。  网格重新划分不需要几何。 ...

具体是怎么操作呢，求版主解答  谢谢

wangyuewu

zkong 发表于 2012-5-25 11:12
方便的话，把模型贴上来，我可以帮你看看。

模型已经传上来了 麻烦有时间帮看看吧  非常感谢了

zkong

其实你的问题是，形状优化后怎么得到优化的设计。

两种方式得到:

方式1:  
HyperMesh -> post -> apply result ->

"Simulation" 选择优化最后一个迭代步的结果
"Data type" 选择 "shape change",
"Nodes" 选择所有节点 "all",

点"apply"按钮，就会用形状优化的结果更新节点坐标，然后remesh就是用HyperMesh的问题了。



方式2:
在HyperMesh打开原模型的情况下，导入.grid 或者_err.grid文件， 覆盖原模型的坐标(选中FE overwrite, 见下图)。*.grid和 *_err.grid的区别是，.grid存的是最后一个通过单元质量检查的迭代步的坐标，_err.grid存的是最后坏掉的坐标。 根据需要，你可以选择其中任意一个作为remesh的起点。

wangyuewu

zkong 发表于 2012-5-30 11:58
其实你的问题是，形状优化后怎么得到优化的设计。

两种方式得到:

谢谢版主的解答，还有一个问题，我利用3D>tetra mesh>volume tetra的办法画的四面体网格，如何进行remesh呢？谢谢！

zkong

wangyuewu 发表于 2012-5-30 14:39
谢谢版主的解答，还有一个问题，我利用3D>tetra mesh>volume tetra的办法画的四面体网格，如何进行remesh ...

抱歉，我HyperMesh用得很烂，这个问题可以到HyperMesh版去问。

wangyuewu

zkong 发表于 2012-5-30 14:47
抱歉，我HyperMesh用得很烂，这个问题可以到HyperMesh版去问。

谢谢斑竹了