Welcome to the forums. Please post in English or French.

You are not logged in.

#1 Re: Code_Aster usage » How to add a nonlinear behaviour in STAT_NON_LINE analysis » 2022-01-19 14:48:05

AsterO'dactyle wrote:

Hello,


Elastic non-linear with increasing stiffness is ELAS_HYPER (hyperelastic material as Mooney-Rivlin, neo-hookean or Signorini model)
You cannot use TRACTION. It's only for (visco)-plastic material (generalized standard materials) with strain hardening


Dear Aster'O'dactyle i have a question about it:
is it really ELAS_HYPER non-linear in terms of stiffness?  As written in u4.43.01.

"Le matériau est élastique incompressible en petites déformations si on prend C10 et C01 tels que 6*( C01+ C10)= E , où E est le module de Young."

So E doesn't change according to epsilon, the same goes for K parameter. Am i wrong or basically the material is treated with a constant E modulus across all the strains ?

I'm asking this since i have problems with strong compression of a rubber and the material doesn't reach the stiffness it should have.

Thank you in advance

Andrea

#2 Re: Code_Aster usage » [SOLVED] a little bug in IMPR_RESU with CONCEPT=( in 10.06 version » 2021-07-27 14:24:26

The same error is now back in CA 14.6 and CONCEPT is no more written in U7.05.21.
Is the CONCEPT command gone purposely or not?  if purposely is there any workaround to achieve the same result?

ERREUR A L'INTERPRETATION DANS ACCAS - INTERRUPTION
>> JDC.py : DEBUT RAPPORT
CR phase d'initialisation
   Etape : IMPR_RESU    ligne : 1043    fichier : 'fort.1'
      !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
      ! Mots clés inconnus : CONCEPT !
      !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
      b_fmt_med
         RESU
            !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
            ! Règle(s) non respectée(s) :                                               !
            ! - Il faut au moins un mot-clé parmi : ('CHAM_GD', 'RESULTAT', 'MAILLAGE') !
            !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
         Fin RESU
      Fin b_fmt_med
   Fin Etape : IMPR_RESU
fin CR phase d'initialisation

thanks in advance.

Andrea

#3 Re: Code_Aster usage » How to assign the thermal conductivity? » 2021-05-04 17:55:04

mathieu.courtois wrote:

Hello,

I think you must define the conductivity as a function of a "dependency variable" (a.k.a. "variable de commande").

Example: ssnp160a

LAMBDA = FORMULE(NOM_PARA=('NEUT1',),VALE='NEUT1')
ACIERTH=DEFI_MATERIAU(THER_FO=_F(LAMBDA=LAMBDA, RHO_CP=RHOCP));

=> lambda = f(NEUT1) (=NEUT1)

and then, define the field of NEUT1:

NEUTG=CREA_CHAMP( ...)

CMTH=AFFE_MATERIAU(MAILLAGE=MA,
                  AFFE=_F(GROUP_MA='CUBE_G',MATER=ACIERTH,),
                     AFFE_VARC=(_F(GROUP_MA='CUBE_G',
                                  CHAM_GD=NEUTG,
                                  NOM_VARC='NEUT1'),),
                  )

You can set the field by evaluating a formula on the geometry coordinates (see CREA_CHAMP/EVAL examples).

According to changelog 11.2.13, I think the usage of dependency variables is very limited in thermics (THER_LINEAIRE and not THER_NON_LINE) while it's widely available in mechanics.

Hope this helps.

MC


Hi Mathieu,
can you confirm that also in v15 the ther_non_line solver can't use NEUT variables ?  I'm struggling with this problem in order to understand what is feasibile and what is not.
Thank you in advance.

#4 Re: Code_Aster usage » Error with conductivity using AFFE_VARC NEUT1 and custom field » 2021-05-04 16:53:37

Dear all,
i'm still investigating such error.
Does anyone knows how to access directly to the lambda field (as in a common thermal simulation) ?
maybe understanding its real values i can understand the source of the interpolating error.


Another question, i tried to define lambda(NEUT1) field since it is possible to do it with E modulus (as stated in many documents), is it possible that thermal solver doesn't accept any kind of customized field lambda?

[EDIT] I finally found this post: www  .code-aster.org/forum2/viewtopic.php?id=19423
from 2014, so it is confirmed that the main reason is that i'm using a THER_NON_LINE simulation instead of a THER_LINEAIRE?




Thank you in advance for any suggestion,

Andrea

#5 Re: Code_Aster usage » Error with conductivity using AFFE_VARC NEUT1 and custom field » 2021-04-29 13:07:18

It seems, and i don't know why, that the problem arises in what the interpolation for the identity functions reads.

TEST1 =  no more identity function, constant function = 1e-9,1e-3,1e+9,1e-3  : every value of x (my custom field) y(x) is 1e-3   =  results are in agreement with 1e-3 values.

TEST2 =  playing with limits, identity function bounded =  1e-5,1e-5,1e-3,1e-3  left and right = 'CONSTANT'   = results are in agreement with value 1e-3

TEST3 = playing again with limits, identity function bounded =  1e-5,1e-5,1e-3,1e-3  left = 'CONSTANT' and right = 'LINEAIRE' = results totally non-sense.

Up to now non-sense means =  value of my function far larger than the maximum found in my custom field CHLAM (values >1e-3 maybe 1 or 10, i don't know). And moreover they seems pretty uniform.

currently the RESU with CHLAM has values between [0.0001 - 0.00065]. In which way can code_aster read such greater values from CHLAM ?

#6 Code_Aster usage » Error with conductivity using AFFE_VARC NEUT1 and custom field » 2021-04-29 10:41:46

Pisolino
Replies: 2

Dear all,
i'm trying to perform a particular simulation (injection molding filling) using a specific thermal conductivity (lambda) defined by a formula, evaluated as a field and then given as affe_varc + identity function (as performed i.e. here /viewtopic.php?id=20385 and in CA documentation u4.43.03) .

I can't understand why the solver doesn't take in count such conductivity. No errors are given, the simulation goes till the end but the results are completely nonsense. Moreover if i change the formula effectiveVisc (i.e. i double or tenfold the value) no numerical change at all in the result. I also exported the custom field CHLAM as a resu and the values are exactly what they should be, so the EVAL step is ok i suppose.

Can you give me any suggestion on where to look? is this because i use a THER_NON_LINE ?

code extraction:

def effectiveVisc(x1):
                         efV= (x1**2)/(12*viscosity/100000);                           
                         return efV;

FMULT1= FORMULE(NOM_PARA=('X1'),VALE='effectiveVisc(X1)',effectiveVisc=effectiveVisc, viscosity=viscosity,);

CHFMU1=CREA_CHAMP(OPERATION='AFFE',TYPE_CHAM='NOEU_NEUT_F', MODELE=mod_t, #PROL_ZERO='OUI',
                  AFFE=_F(TOUT='OUI', NOM_CMP=('X1'),VALE_F=FMULT1));

CHLAM=CREA_CHAMP( OPERATION='EVAL', TYPE_CHAM='NOEU_NEUT_R', CHAM_F=CHFMU1, CHAM_PARA=(CHTH));

lam_f= DEFI_FONCTION(NOM_PARA='NEUT1',VALE=(-1e9,-1e9,1e9,1e9),PROL_DROITE='LINEAIRE',PROL_GAUCHE='LINEAIRE');

matpar=DEFI_MATERIAU(THER_NL=_F(LAMBDA=lam_f, RHO_CP=dens_f,),);

mat_t=AFFE_MATERIAU(MODELE=mod_t,
                      AFFE= affemat,
                      AFFE_VARC=(
                                 _F(NOM_VARC='TEMP', VALE_REF=1.0,),
                                 _F(NOM_VARC='NEUT1',CHAM_GD=CHLAM),  #variabile viscosità specifica                               
                                 )
                               );

for i in range(1,len(listatempo1)-1):
  tnlin=THER_NON_LINE(reuse= tnlin,
                      MODELE=mod_t,
                      CHAM_MATER=mat_t,
                      COMPORTEMENT=_F(RELATION='THER_NL',),
                      EXCIT=load1,
                      INCREMENT=_F(LIST_INST=tempo,
                                   INST_INIT=listatempo1[ i ],
                                   INST_FIN= listatempo1[ i+1 ],),
                      ETAT_INIT=_F(EVOL_THER=tnlin, INST=listatempo1[ i ]),
                      NEWTON=_F(REAC_ITER=1,),
                      PARM_THETA=theta,
                      CONVERGENCE=_F(ITER_GLOB_MAXI=200,
                                     RESI_GLOB_RELA=5e-4,),);

#7 Re: Code_Aster usage » Heat transfer with ECHANGE_PAROI » 2021-03-17 19:44:40

dezsit wrote:

Dear Community,

I would like to rise up an old problem again, because I have not seen progress in it. (Or we missed something during the years?)

So my question is: Is it planed to handle the thermal contact (more precisely thermal resistance between walls - ECHANGE_PAROI) on based on a master-slave pairing approach? (I suppose the pairing algorithm of for example LIAISON_MAIL could be used easily in this case, which is already present in the AFFE_CHAR_THER operator). We could save a lot of partitioning/meshing work with that.



Thank you
and
Best regards,
DT.

Hi everybody,
are there any news about this topic? it would be really great to have a coupling between 2 non paired surfaces to define a thermal contact resistance. I remember something very easy in syrthes to define such contact thermal resistance between 2 surfaces.

thanks in advance,

Andrea

#8 Re: Code_Aster usage » [SOLVED] PRE_SIGM use for Preload » 2020-01-02 17:32:43

Hi all,
i've tried to follow your advice to model the pre-load of a bolt in an example very similar to the one attached before. However in the very first iteration my bolt gets elongated by the value of the stress/modulus. What i wanted to do however was to say to CA that the linear element @start has already a stress "inside" at that fixed elongation.
Am i doing something wrong ?

For example i tried to set a negative PRE_SIGM value but then the simulation does not converge in the first iteration.


Thank you in advance,

Andrea

#9 Re: Code_Aster usage » Problems with DYNA_VIBRA (<EXCEPTION> <ALGORITH3_13> ) » 2019-01-17 23:58:27

SOLVED
the error is due to the new default settings of  SCHEMA_TEMPS in DYNA_VIBRA that requires MASS_MECA_DIAG elements.
The correct settings to work with MASS_MECA and usual elements (in my case both DKT, 3D and COQUE_3D) is :

SCHEMA_TEMPS=_F(SCHEMA='NEWMARK'),

hope it can be useful.

Andrea

#10 Re: Code_Aster usage » Problems with DYNA_VIBRA (<EXCEPTION> <ALGORITH3_13> ) » 2019-01-15 21:28:51

Hi,
i'm encountering the same error. I'm stuck in the same loop of MATR_MASS and MATR_MASS_DIAG. I want to simulate a seismic sinusoidal acceleration in X direction in DYNA_TRAN using a gravity load. I'm trying to follow on a youtube tutorial:

YOUTUBELINK/watch?v=W8uDIw7ONJM&t=656s

Differently from the tutorial I'm working with shell DKT elements. I'm using Salome_meca 2018.0.

Can someone help me?

Thank you in advance

Andrea

#11 Re: Code_Aster usage » Hot isostatic pressing (HIP) process simulation » 2017-11-23 19:36:04

Hi all,
after some days and some experiments i'm starting to find a way to reach my goal.
I'm now trying with more precise questions.

1- the model i need is elasto-plasto-viscoplastic. I read that is possible to chain lemaitre behaviour to vmis_isot_line. "chaining" means that i have to perform a stat_non_line "stat1" with behaviour vmis_isot_line and then run another stat_non_line with behaviour (COMPORTMENT) lemaitre with reuse - evol_noli - "stat1" ?? is it right?

2- density. i need to set E,v,alpha, and all the other properties depending on Temperature (and this is not a problem) and relative density. Should i create a new field with CREA_CHAM ? and then is possible to link this added variable into fonctions (or nappe since its T and rho_r) ?

Thanks in advance,

Andrea

#12 Code_Aster usage » Hot isostatic pressing (HIP) process simulation » 2017-11-20 00:24:59

Pisolino
Replies: 1

Hi all,
i would like to simulate the HIP process with code aster (en.wikipedia.org/wiki/Hot_isostatic_pressing). I've been reading some scientific papers about the best ways to perform this simulation. One of the most effective way since now is to use an elasto-visco-plastic model. Many examples are done using Norton creep with data obtained from the full densified material. Moreover it is important to correctly model a variable relative density from 0 to 1. since we start from values like 0.7 and we aim to obtain 1 all over the solid. Can i ask some hints on how to consider the relative density variation in my model ? If my question is too rough or unclear i will try to explain better.

Thanks in advance.

Andrea

#13 Re: Code_Aster installation » Code_aster mpi: is it worth? general FAQs » 2016-10-11 17:19:00

Hi mathieu,
thanks for your reply. However during the run of my optimization .comm file I can see only one thread active during such small operations.  Maybe there is something I miss to set properly ?

Thank you

#14 Code_Aster installation » Code_aster mpi: is it worth? general FAQs » 2016-07-26 14:34:07

Pisolino
Replies: 3

Hi all,
before starting to install the mpi version of code_aster 12.6 (or maybe 13.1) i have some questions that i was not able to answer just by googling or using my very low coding skills. I'm on ubuntu 16.04, Intel® Xeon(R) CPU E5-2620 v2 @ 2.10GHz × 24 , 125.9 GiB RAM

1) the new salome-meca 2016 has, , differently from before , an open-mp pre compiled version of Aster? since i could use MUMPS for MECA-STATIQUE and seeing all the 24 threads (12 core) finally working.

2) unfortunately my simulation is an optimization loop (topological optimization) so the meca_statique part is very small compared to the rest of the loop. From what i read mpi-aster is considered FULL parallel. Does it means that operations like CREA_CHAMP can be parallelized with an mpi?

3) I have a workstation and not a cluster. MPI-aster is worthy useful only with an HPC / cluster with many nodes or on every PC (with more than 1 core) ?

Really thanks in advance

Andrea

#15 Re: Code_Aster usage » Topological Optimization » 2016-07-18 15:04:49

Hi TMK,
Your developing idea is pretty good, i've been quite busy during the last weeks but i want to help you for sure. I can give you some support chatting via skype (easier to check for me). Send me a PM with your skype contact if you want and i'll add you. The same for everyone that want to improve this algorithm together.

Andrea

#16 Re: Code_Aster usage » thermique : imposer un flux non linéaire dépendant de 'X' et 'Y' » 2016-05-05 14:41:11

Hi all,
is possible to define a flux depending on time and temperature? (INST,TEMP) ? i'm trying to simulate a thermal cycle but probably im'm doing in the wrong way.

#-------- HEAT EXCHANGE FUNCTION DEFINITION ---------
#hT is function of the wall temperature (TEMP), at the right and left (DROITE / GAUCHE) border the value is mantained constant

hT1=DEFI_FONCTION(
                 NOM_PARA='TEMP',
                 VALE=(300.0 ,300.0 ,
                       385.0 ,425.0 ,
                       453.0 ,2574.0 ,
                       676.0 ,709.0 ,
                       900.0 ,138.0 ,),
                 PROL_DROITE='CONSTANT',
                 PROL_GAUCHE='CONSTANT',);

hT2=DEFI_FONCTION(
                 NOM_PARA='TEMP',
                 VALE=(300.0 ,300.0 ,
                       385.0 ,425.0 ,
                       453.0 ,2574.0 ,
                       676.0 ,709.0 ,
                       900.0 ,138.0 ,),
                 PROL_DROITE='CONSTANT',
                 PROL_GAUCHE='CONSTANT',);

hT3=DEFI_FONCTION(
                 NOM_PARA='TEMP',
                 VALE=(300.0 ,300.0 ,
                       385.0 ,425.0 ,
                       453.0 ,2574.0 ,
                       676.0 ,709.0 ,
                       900.0 ,138.0 ,),
                 PROL_DROITE='CONSTANT',
                 PROL_GAUCHE='CONSTANT',);

#-------- THERMAL FLUX FUNCTION DEFINITION ---------

cost=37/2;
Pi=3.14;
T_al=600;

def h_ext(TEMP, INST):
                                if sin(INST/cost*Pi)> 0.751: return hT1(TEMP)
                if sin(INST/cost*Pi)< -0.751: return hT2(TEMP)
                return hT3(TEMP)

FH_T= FORMULE(NOM_PARA= ('TEMP','INST'), VALE='h_ext(TEMP,INST)',);

#-------  INTERPOLATION OF THE HEC FORMULA --------
#first is necessary to define the interpolation domain via a DEFI_LIST_REEL, next the values are interpolated with CALC_FONC_INTERP

temp1=DEFI_LIST_REEL(DEBUT=0.0,
                  INTERVALLE=_F(JUSQU_A=700.0,
                                NOMBRE=100,),);
time1=DEFI_LIST_REEL(DEBUT=0.0,
                  INTERVALLE=_F(JUSQU_A=1800.0,
                                PAS=0.5,),);

hec=CALC_FONC_INTERP(FONCTION=FH_T,
                      LIST_PARA=temp1,
                      NOM_PARA='TEMP',
                      INTERPOL='LIN',
                      PROL_DROITE='CONSTANT',
                      PROL_GAUCHE='CONSTANT',
                      NOM_PARA_FONC='INST',
                      LIST_PARA_FONC=time1,
                      PROL_DROITE_FONC='CONSTANT',
                      PROL_GAUCHE_FONC='CONSTANT',);

and the error is

CR d'execution de JDC en MIXTE
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   ! <S> Exception utilisateur levee mais pas interceptee.        !
   ! Les bases sont fermees.                                      !
   ! Type de l'exception : error                                  !
   !                                                              !
   !  les fonctions interprétées doivent être tabulées auparavant !
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
fin CR d'execution de JDC en MIXTE

Thanks in advance,

#17 Re: Code_Aster usage » Structural Optimization with CA - HOMARD » 2016-04-19 12:11:36

Hi Johannes,
now i'm having another problem with HOMARD still related to double faces however this time the output mess is more clear to understand:

   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   ! <S> Exception utilisateur levee mais pas interceptee.                        !
   ! Les bases sont fermees.                                                      !
   ! Type de l'exception : error                                                  !
   !                                                                              !
   !  La face numéro 1 de la maille M13000 possède 2 voisins.                     !
   !  Or, ceci n'est pas prévu : un seul voisin par face est autorisé.            !
   !  Cela signifie que la maillage comporte des mailles doubles ou qu'une maille !
   !  de bord est intercalée entre deux mailles volumiques.                       !
   !                                                                              !
   !  Conseils :                                                                  !
   !    Vérifiez votre maillage.                                                  !
   !    Ces mailles doubles ou intercalées ne sont peut-être pas nécessaire au    !
   ! calcul.                                                                      !
   !    Dans ce cas ne les affectez pas avec le modèle !                          !
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 

I have a mesh composed by 2 volume group_ma, based on what i understood there is a problem at the interface between the 2 volumes, how can i solve this issue? if i remove the internal faces i can't have 2 different volume groups.

 NEWMAIL=CO('M%d'%k)
		                		NEWCH=CO('C%d'%k)
		                		MACR_ADAP_MAIL(MAILLAGE_N=MAIL,
		                		     MAILLAGE_NP1=NEWMAIL,
		                		     ADAPTATION='RAFF_DERA',
		                		     #RESULTAT_N=res,
		                		     #NOM_CHAM='ERME_ELEM',
		                		     CHAM_GD=CHERRE,
		                		     NOM_CMP='ERREST',
		                		     MAJ_CHAM=_F(CHAM_MAJ=NEWCH,
			        		                  CHAM_GD=CHRATIO,
			        		                  TYPE_CHAM='NOEU_NEUT_R',
			        		                  #TYPE_CHAM='NOEU_ENEL_R',
			        		                  #RESULTAT=res,
			        		                  #NOM_CHAM='ENEL_NOEU',
 			        		                  #NUME_ORDRE=1,
			        		                  ),
		                		     CRIT_RAFF_PE=0.06,
 		                		     CRIT_DERA_PE=0.02,
		                		     NIVE_MAX=NIVE_MAX,
		                		     DIAM_MIN=DIAM_MIN,
		                		     INFO=1,
		                		     NOMBRE='OUI'
		                		     #QUALITE='OUI',
		                		     #CONNEXITE='OUI',
		                		     #TAILLE='OUI',
						                     ); 

 

Thanks in advance.

Andrea

#19 Re: Code_Aster usage » Topological Optimization » 2016-04-13 17:53:26

For anyone that is interested in topological optimization I have just finish to update the old-but-good comm file developed by Frederic regarding a SIMP 3d optimization, performed with energy or stress optimization. I tried to put enough comments to allow future improvements.

SIMP method useful links =
- http://www.cmap.polytechnique.fr/~allai … dSIMP.html
- https://books.google.it/books?id=NGmtmM … on&f=false

commfile is attached

#20 Re: Code_Aster usage » Probleme copie de fichier fort.xx en cours de calcul » 2016-04-11 11:02:22

HI all,
is this approach still valid in salome 2015.2 ? I mean i'm using an old script containing such instructions :

 CHCOM=CREA_CHAMP( OPERATION='EVAL', TYPE_CHAM='NOEU_NEUT_R', CHAM_F=CHFMU1, CHAM_PARA=(CHENEL,CHRATIO));
        IMPR_RESU(FORMAT='MED',UNITE=81,RESU=_F(MAILLAGE=MAIL,CHAM_GD=CHCOM,),);
        IOFile.SaveFileLinux(UNITE=81,FILENAME="CHCOM%d.med"%k) 

Here is the function:

 sortiepathlinux="/media/pisa/interno/ASTER/OPTI_TOPO/SIMP/tmp/";

def SaveFileLinux(UNITE=80,PATH=None,FILENAME="file.med"):
       print "debug savefilelinux"
       mpath="";
       if (PATH!=None):
              mpath=PATH
       else:
              mpath=sortiepathlinux       
       print 'cp fort.'+str(UNITE)+' '+mpath+FILENAME;
       EXEC_LOGICIEL(LOGICIEL='cp fort.'+str(UNITE)+' '+mpath+FILENAME,INFO=1);
       print 'cp fort.'+str(UNITE)+' '+mpath+FILENAME+'dinuovo';
       DEFI_FICHIER(ACTION='LIBERER',UNITE=UNITE,)
       EXEC_LOGICIEL(LOGICIEL='rm fort.'+str(UNITE),INFO=1)

Here is the error:

 Commande :
   ['cp fort.81 /media/pisa/interno/ASTER/OPTI_TOPO/SIMP/tmp/CHCOM1.med']
   
   !---------------------------------------------------------------------------------------------------------------------------------!
   ! <EXCEPTION> <SUPERVIS2_5>                                                                                                       !
   !                                                                                                                                 !
   ! Erreur inattendue lors de l'exécution de la commande 'EXEC_LOGICIEL'.                                                           !
   ! Merci de signaler cette anomalie.                                                                                               !
   !                                                                                                                                 !
   ! Erreur :                                                                                                                        !
   !                                                                                                                                 !
   ! File "/home/pisa/salome_meca/V2015_2/tools/Code_aster_stable-1240/lib/aster/Macro/exec_logiciel_ops.py", line 336, in           !
   ! exec_logiciel_ops                                                                                                               !
   !     action.execute()                                                                                                            !
   !   File "/home/pisa/salome_meca/V2015_2/tools/Code_aster_stable-1240/lib/aster/Macro/exec_logiciel_ops.py", line 100, in execute !
   !     self.executeCommand()                                                                                                       !
   !   File "/home/pisa/salome_meca/V2015_2/tools/Code_aster_stable-1240/lib/aster/Macro/exec_logiciel_ops.py", line 125, in         !
   ! executeCommand                                                                                                                  !
   !     output, error, exitCode = self.executeCmdLine( cmd, capture, silent )                                                       !
   !   File "/home/pisa/salome_meca/V2015_2/tools/Code_aster_stable-1240/lib/aster/Macro/exec_logiciel_ops.py", line 117, in         !
   ! executeCmdLine                                                                                                                  !
   !     process = subprocess.Popen(cmd, shell=self.shell, **options)                                                                !
   !   File "/home/pisa/salome_meca/V2015_2/prerequisites/Python-273/lib/python2.7/subprocess.py", line 679, in __init__             !
   !     errread, errwrite)                                                                                                          !
   !   File "/home/pisa/salome_meca/V2015_2/prerequisites/Python-273/lib/python2.7/subprocess.py", line 1249, in _execute_child      !
   !     raise child_exception                                                                                                       !
   ! [Errno 2] No such file or directory                                                                                             !
   !---------------------------------------------------------------------------------------------------------------------------------!
   

Any suggestion ? thx

Andrea

#21 Code_Aster usage » Structural Optimization with CA - HOMARD » 2016-04-04 14:20:43

Pisolino
Replies: 3

Hi all,
i've just found an old but excellent .comm file developed by user Bracchesimo that perform a kind of topological optimization. It was built for 2010 version of salome meca. I succeeded in updating the comm file to the 2015.2 SALOME version. Now i would like to increase the speed of such algorithm. My idea was to refine the mesh every iteration step with MACR_ADAP_MAIL, maybe with a refinement linked to the stress density.
I'm having an error with HOMARD :

Les deux elements suivants sont doubles :
Numero dans le calcul   :     115920
Famille MED :   -4 ; type :  102
  Ses noeuds (numerotation du calcul) :
      2492       112
  Ses noeuds (numerotation HOMARD) :
      2492       112
Numero dans le calcul   :     115688
Famille MED :    0 ; type :  102
  Ses noeuds (numerotation HOMARD) :
      2492       112
  Ses noeuds (numerotation du calcul) :
      2492       112
===== HOMARD ===== STOP =====
----- fin stderr -----------------------------------------------------------------
   
   !-------------------------------------!
   ! <EXCEPTION> <EXECLOGICIEL0_3>       !
   !                                     !
   ! Code retour incorrect (MAXI 0) : 10 !
   !-------------------------------------!

Can someone help me to solve such error?? thank you so much.

Moreover, the original algorithm is not very scalable, the longest steps are : POST_RELEVE_T and CREA_MAILLAGE that are using 1 single core. Is there any new command / particular way to make all the simulation go faster??

Thanks in advance

#22 Re: Code_Aster usage » Structure Optimization with Code_Aster » 2016-03-31 17:24:28

Solved ! it was necessary to use an incremental name (counter int converted to string) for NOM and PREFI_MAILLE :

#create the new mesh copying the current one             
        msh[i+1]=CREA_MAILLAGE(MAILLAGE=msh[ i ],
                    CREA_MAILLE=_F(NOM=str(i),
                    TOUT='OUI',
                                        PREF_MAILLE=str(i),
                    PREF_NUME=i+1),
                            );

#23 Re: Code_Aster usage » Structure Optimization with Code_Aster » 2016-03-31 15:19:34

Hi,
starting from the commfile written by Bracchesimo i tried to update the algorithm to the salome-meca 2015.2 version. I already come through some errors but now i'm stuck with an error (that i believe it is quite stupid but i can't solve it). Until now i had to switch from Numeric to numarray, then updating CALC_ELEM in CREA_CHAMP. The actual error is driven by CREA_MAILLAGE, for what i understood something can't update during the for loop, indeed the error appears at the second loop.

CR phase d'initialisation
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   ! <S> Exception utilisateur levee mais pas interceptee. !
   ! Les bases sont fermees.                               !
   ! Type de l'exception : error                           !
   !                                                       !
   !  Erreur données : maille déjà existante :  A2         !
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
fin CR phase d'initialisation

this is the code part:
#create the new mesh copying the current one             
        msh[i+1]=CREA_MAILLAGE(MAILLAGE=msh[ i ],
                    CREA_MAILLE=_F(NOM='msh',
                    TOUT='OUI',
                                     PREF_MAILLE='A',
                    PREF_NUME=i+1,),);

Any suggestion??
Thanks

Andrea

#24 Code_Aster usage » Topological Optimization » 2016-03-30 16:37:59

Pisolino
Replies: 11

Hi all,
i'm wondering if any progress has been made in Code Aster to manage the topological optimization. I found some interesting examples but still no particular hints about where to start.
OGACA is the most interesting http://www.code-aster.org/V2/spip.php?article314 and it is done with CA, there is also a SALOME-MECA Examples by Bracchesimo : Structure Optimization with Code_Aster but the link is now not valid so i cannot find any comm file. Any idea, help or hint is very appreciated. I know it will not be easy for me to obtain something similar but i want to try with all my efforts.

Thanks in advance.

Andrea

#25 Re: Code_Aster usage » Projection of results from a 2d axys to the 3d solid , PROJ_CHAMP » 2016-03-21 19:08:43

I just used that filter (transform) to rotate a quarter of a gear but i was wondering if there is something easier (imagine if you have to elaborate very small slices and multiply them). Moreover i was not able to find (maybe i'm not using the correct search-keys) such a function that project the 2d results around the axys to obtain the 3 solid.