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.