First step for shape optimization?

I tried shape optimization.
Original article is here

COMM file is modified to do on Code Aster13.6.

Changing items are as follows:

  • Comment out Numeric
import numpy as N
#import numeric as N
import math
  • Addition of 1 dim element

After reading mesh data,  a node group ’TOUT’ is added.

msh[0]=DEFI_GROUP(reuse =msh[0],
MAILLAGE=msh[0],
CREA_GROUP_NO=_F(GROUP_MA='model',NOM='TOUT',) ,);

Then  1 dim elements are added on the nodes shown in ‘TOUT’.

nmh[i]=CREA_MAILLAGE(MAILLAGE=msh[0],
                     CREA_POI1=_F(NOM_GROUP_MA='spElmt',
                                    GROUP_NO='TOUT',),);
  • a procedure to make new mesh data

Not msh[i]⇒msh[i+1], renew a calculated mesh group in msh[0].

msh[0]=DEFI_GROUP(reuse =msh[0],
                  MAILLAGE=msh[0],
                  DETR_GROUP_MA=_F(NOM='model',),
                  DETR_GROUP_NO=_F(NOM='TOUT',),
                  CREA_GROUP_MA=_F(NOM='model',MAILLE=NewModel,),
                 CREA_GROUP_NO=_F(GROUP_MA='model',NOM='TOUT',) ,);

 

This is a snapshot of ASTK.

Result of 1st step.

10th step.

20th step

30th step

Final (34step)

A element which are low stress are removed.





Data is Optiastkv13-6.

2 Replies to “First step for shape optimization?”

コメントを残す

メールアドレスが公開されることはありません。 * が付いている欄は必須項目です