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

You are not logged in.

#1 Re: Salome-Meca usage » Semi-coherent meshed part » 2022-03-22 15:42:29

Hi,

Assuming you are using Salome, just partitioning your geom,
mesh it as one whole (this makes fully conformal mesh on all internal faces of the partitioned geometry),
then use the Duplicate Nodes or/and Elements from the menu Modification->Transformation.
Either the 1st or the 4th option could work, but require different preparation (grouping) and post treatment (re-merging of some nodes).

This process will result fully identical element shapes and node location on the both sides of your interface, but with duplicated nodes! The rest of your mesh is fully merged.

BR.
dezsit

#2 Salome-Meca usage » Presentations of SMUD2021 » 2021-04-13 10:30:25

dezsit
Replies: 0

Dear SMUD Organizers,

Is it planned to share the presentations from the last user meeting somewhere?

Thank you in advance!
BR,
dezsit

#3 Re: Code_Aster usage » Different stiffness in compression and tension for K_T_D_L » 2021-01-27 14:15:40

Dear Anirudh,

Your question of course absolutely valid, and can be defined such behavior in aster.
A similar one is shown by Johannes Ackva (code-aster/dot/de), using the DIS_CONTACT material "behavior".

In the attached small example I used grounded springs (POI1). We need two of them connected parallel, one for the tension, one for the compression. The key is that the local csys (REPERE='LOCAL") should be opposite (ORIENTATION in AFFE_CARA_ELEM).
In this example we have stiffness equal to 1(force/displacement) in tension and 10 in compression. You can check it with xmgrace (res.func) or numerically (res.tab).

The line spring (SEG2) works similarly.

I would not say it will work in all cases, but in this simple one it provides the expected result. Tweeking with the parameters (eg. JEU in DEFI_MATERIAU) and combining other springs, very different characteristics can be achieved.

Regards,
dezsit

#4 Re: Code_Aster usage » ERROR in beam AFFE_CARA_MECA » 2020-09-03 08:10:11

Hi,

I checked only your comm file (and not the whole study). Two comments:

1) You have to identify in the solution sequence which element characteristics should be used:

eg.:

PP = MECA_STATIQUE(identifier=u'14:1',
                   CARA_ELEM=Horm_Th #!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
                   CHAM_MATER=MATS,
                   EXCIT=(_F(CHARGE=BC),
                          _F(CHARGE=Grav),
                          _F(CHARGE=BC_Horm)),
                   MODELE=model)

2) Since in one solution sequence you can use only one CARA_ELEM, you have to describe all your element characteristic in one AFFE_CARA_ELEM,

As I noticed, the element in GrPuntal_Edges group used in two different AFFE_CARA_ELEM (Horm_TH and Punt), with the same characteristic, as I mentioned both cannot be used in the same solution (and in this case it does not really have any sense, because the properties are the same).

Hope this helps,
Best regards,
dezsit

#5 Re: Code_Aster usage » [RESOLU] Define "GROUP_NO" for POST_RELEVE_T » 2020-09-03 07:37:28

Hi,

at this moment, if you create group by DEFI_GROUP, asterStudy does not recognize it, since it will be available during run time only. (I guess it could be handled, maybe in the future releases).
So you should put manually the new names at the bottom of the group selection widget (Manual selection) to use it.

Best regards,
dezsit

#6 Re: Salome-Meca installation » salome-meca source code? » 2020-09-03 07:28:40

Hi,

try this:
sourceforge.net/projects/codeaster/

Best regards,
dezsit

#7 Re: Code_Aster usage » [SOLVED] Postprocess: stresses in model consisting of 3D & 2D elements » 2020-06-24 14:01:57

Hi,

you can manage it with one file, just put all your resulting fields (DKT,3D,...) into one IMPR_RESU.

Then in paraVis, open the resulting med twice (or more, depending on the study), and select your 3D fields on one of the result file, than the 2D fields on the other result file.
See the attached picture,
here you can see mech.rmed file opened twice (it is one file from one IMPR_RESU),
on the first file I selected the 3D results: reslin_SIEQ_NOEU from ComSup2 in the selection window on the left bottom,
then on the second file I selected the DKT results, here the TOP____SIEQ_NOEU 8which contains the equivalent stresses on the top side of my shell),
Then visualized both together (VMIS filed on both) with adjusting the scales.

You can do that for any kind of combination (eg with beams or what ever), just sometimes a good visualization requires more tweeking in paraVis.

Best regards,
dezsit

#8 Re: Salome-Meca usage » In windows, How to run asterstudy from terminal? » 2020-06-21 16:11:41

Hi,,
with this you jump in, then jump out from the as_run directory:

kartikpandya wrote:

hey dezsit;

C:\Users\Lenovo>cd C:\SM-2019-w64\SM-2019-w64-0.3\PRODUCTS\code_aster\bin

C:\SM-2019-w64\SM-2019-w64-0.3\PRODUCTS\code_aster\bin>cd C:\Users\Lenovo\Deskto
p\FEA\study1_Files\RunCase_9\Result-Stage_1

But you would not have to, just simply use

C:\SM-2019-w64\SM-2019-w64-0.3\PRODUCTS\code_aster\bin\as_run export

from the directory, where your export file is placed.

So in your case you should do (if the export file is in the directory ..\Result-Stage_1):

cd C:\Users\Lenovo\Desktop\FEA\study1_Files\RunCase_9\Result-Stage_1

then

C:\SM-2019-w64\SM-2019-w64-0.3\PRODUCTS\code_aster\bin\as_run export


This is not a salome/code_aster problem, it is simply how to use the command line.

Best regards,
dezsit

#9 Re: Code_Aster usage » How to do these options in ParaVis » 2020-06-19 17:41:30

Hi,

Warp by Vector can do that, but you need a real vector with 3 elements!
In case of beams you have 6 in the DEPL filed: 3 translation 3 rotation, so you should produce a real vector from that.
You have two options to do that
1) the easy way: when you open your med file and select your filed (e.g. depl of course) , you have to tick the checkbox "Generate vectors"  (just below the field selection window). This will generate a new DEPL field (DEPL_Vector or something) what will contain the 3 necessary directional displacements, and can be used by Warp filter

2) the manual way: use the calculator filter, then apply the following expression: ihat*DX + jhat*DY + khat*DZ (DX,DY,DZ is the scalars of your DEPL filed),
this will produce a new vector field with 3 components under the name you provides in the Calculator filter (by default it is Result).

Note: you should use one of the option in any case where the DEPL field contains more than 3 directional displacement, so e.g. in case of shells, or if you use contacts.

Best regards,
dezsit

#10 Re: Salome-Meca usage » In windows, How to run asterstudy from terminal? » 2020-06-19 17:30:39

Hi,

Sorry, I made a mistake in the path, the bin directory is missing.
So you should try this on your system:

C:\SM-2019-w64\SM-2019-w64-0.3\PRODUCTS\bin\code_aster export

(Assuming the SalomeMeca installed here in your case : C:\SM-2019-w64\SM-2019-w64-0.3)

BR,
dezsit

#11 Re: Salome-Meca usage » OpenTurns Usage » 2020-06-18 18:06:18

Hello,

I think the problem is, that you extracted the "extrema" of the displacement DZ with this:

table2 = POST_RELEVE_T(ACTION=_F(OPERATION=('EXTREMA', ),
                                   INTITULE='DZ',
                                   RESULTAT=linel,
                                   NOM_CHAM='DEPL',
                                   PRECISION=1.E-06,
                                   CRITERE='RELATIF',
                                   FORMAT_C='MODULE',),
                         INFO=1,)

This produce not only one DZ value (as necessary for OpenTurns) but four, and actually not simply 'DZ':
#
#--------------------------------------------------------------------------------
##ASTER 14.04.00 CONCEPT table2 CALCULE LE 18/06/2020 A 16:57:37 DE TYPE         
#TABLE_SDASTER                                                                   
INTITULE         RESU     NOM_CHAM         NUME_ORDRE   INST         EXTREMA  NOEUD    CMP      VALE       
DZ               result   DEPL                        1  0.00000E+00 MAX      N98      DZ        6.39285E-02
DZ               result   DEPL                        1  0.00000E+00 MIN      N134     DY       -1.00000E+00
DZ               result   DEPL                        1  0.00000E+00 MAXI_ABS N134     DY        1.00000E+00
DZ               result   DEPL                        1  0.00000E+00 MINI_ABS N52      DX        0.00000E+00

So, that is way a warning reported in the message file:
   ! Parametre absent de la table : DZ.   

So you should try this:

#all dz value extracted from the node group 'Group_1' or whatever into the table
table2 = POST_RELEVE_T(
  ACTION=_F(
    GROUP_NO=('Group_1', ),
    INTITULE='DZ',
    NOM_CHAM='DEPL',
    NOM_CMP=('DZ', ),
    OPERATION=('EXTRACTION', ),
    RESULTAT=result
  )
)

#then it is filtered (only maximum absolute -> MAXI_ABS) during the numpy table creation, so you have only one value at the end.
IMPR_TABLE(
  FILTRE=_F(
    CRIT_COMP='MAXI_ABS',
    NOM_PARA='DZ'
  ),
  FORMAT='NUMPY',
  NOM_PARA=('DZ', ),
  TABLE=table2,
  UNITE=8
)

I did not tried with OpenTurns, but it eliminates the warning, so should work.

Best regards,
dezsit

#12 Re: Salome-Meca usage » In windows, How to run asterstudy from terminal? » 2020-06-18 17:03:02

Hello,

Try this:
SALOME_INSTALL_DIR\PRODUCTS\code_aster\as_run export_file

BR
dezsit

#13 Re: Salome-Meca usage » [SOLVED] Open paravis from terminal like code aster or salome » 2020-04-21 19:37:20

Just to note:

It does not worked to me either, if I try with SalomeMeca2019 (same error as yours)
but the salome shell method works with the native Salome platform (eg. 9.4) and also with older version of salome meca (eg. 2018). Interesting.

The problem is, if you start the paraview in your way, the MEDloader plugin won't be loaded.
In case of salome shell way, you get a standalone paraview with the capability to handle MED files.

Regards,

#15 Re: Salome-Meca installation » Salome Meca 2019 new version (2019.0.1) » 2019-11-04 22:14:48

Up!

Hi, could somebody from the developers comment that we can expect an updated version with OpenTurns?

Thank you and
Best regards,
dezsit

#16 Re: Code_Aster usage » Principal stresses in dynamic analysis » 2019-11-04 22:10:56

Hi,

You should check the CALC_CHAMP command and CRITERES -> SIEQ_NOEU for example. It provides all invariants, including principals, independently from the type of the analysis. After that,  you should use POST_RELEVE_T to produce table, and IMPR_TABLE to write it out to a txt file (you have a lot of option to format it to your needs),
also IMPR_RESU, FORMAT='RESULTAT' will provide you a text file with your specified results.

I would suggest to check Jean Pierre Aubry's book: framabook.org/beginning-with-code_aster/
you will find a lot of example in it, also in the Training section of this site, and in the case files of the aster. Just check the keywords I mentioned above.

BR,
dezsit

#17 Re: Salome-Meca usage » Error generating mesh » 2019-11-04 22:02:26

Hello,

Generally, hex outperforms tetra in mechanics (better mesh convergence with less elements), especially if bending involved. But in some cases hex suffers from various problems (volumetric locking, hourglassing, etc, depending on the numerics behind a specific hex formalism (reduced integration, full integration but linear element, etc.), and how a specific fea package handles this kind of problems (eg. assumed strain as in case of aster) ), you should check the element reference manual of a software package, and do verifications! But sometimes the geometry does not allow to make good hex mesh (or simply the job not worth to spend too much time with mesh preparation) so tetra can be used, but only quadratic for mechanics, and again do verification and mesh convergence study.

Code_aster fatigue case: as I mentioned the fea part of fatigue is "simple" mechanics (statics, dynamics), not aster specific. We do the fatigue evaluation separately, independently from the fea package, so I cant provide an aster specific fatigue evaluation (knowing that aster has in-built fatigue evaluation tool, but we have not used it yet).

BR
dezsit

#18 Re: Code_Aster usage » Mass point » 2019-11-04 21:49:02

Hi,

you defined the DIS_T in separate modelisation (model0) which is not used at all.
You have to define everything in the same model, what you want to use in an analysis, eg:

  model = AFFE_MODELE(MAILLAGE=mesh,
                      AFFE=(_F(PHENOMENE='MECANIQUE',
                              TOUT='OUI',
                              MODELISATION=('3D', ),),
                                 _F(PHENOMENE='MECANIQUE',
                               GROUP_MA='0D',
                               MODELISATION=('DIS_T', ),),),
                      INFO=1,
                      VERI_JACOBIEN='OUI',
                      DISTRIBUTION=_F(PARTITIONNEUR='METIS',
                                      METHODE='SOUS_DOMAINE',),
                      )

BR
dezsit

#19 Re: Code_Aster usage » Command to export MED file » 2019-11-04 21:43:14

Hi,

I am not sure python can handle environment variables (from shell) in this way (I suppose the dump.py thrown the exception of your try/except condition), it is more complicated I think, you cannot pass variable from shell to python script in this way especially because the dump.py is executed by salome session and not by the shell. But to be honest I never tried it. If you follow my example I gave you in an other post, the param dict file can handle this simply (as it is already done in the 144 line of salomeScript.py).

BR
dezsit

#20 Re: Code_Aster usage » [SOLVED] Stress analysis on a symmetrical shaft » 2019-11-04 21:34:07

Hi,

I did not check your case deeply, only the message file, but as JeanPierre said, it is clear: 3D solid does not have rotational DOF (it is not aster specific, 3D solid does not have rotational dof in any common mechanical fea package!) so you should remove the '-DRY-DRZ' from yout DDL_ESCL definition. A second thing related to the orphan node: yes aster does not able to handle orphan nodes. But as it was discussed many times, you can define a 0D element on your orphan node, and use a DIS_T or DIS_TR modelisation, depending on your needs, and attach a zero stiffness spring to this 0D element (AFFE_CARA_ELEM -> K_T_D_N or K_TR_D_N or whatever, again depending on the problem you modeling). So you have to provide stiffness for every node, even if it is zero stiffness.

BR,
dezsit.

#21 Re: Code_Aster usage » Mass point » 2019-10-31 09:19:08

Hi,

you have to create 0D elements (POI in aster) on the nodes, where you want to apply mass element, then use DISCRETE element modelisation on these 0D elements. (more keywords: AFFE_MODELE->DIS_T, AFFE_CARA_ELEM->DISCRETE->M_T_D_N)

Check Jean Pierre Aubry's excellent book: framabook.org/beginning-with-code_aster/, eg. you will find an example on the 31.page.

BR
dezsit

#22 Re: Code_Aster usage » [SOLVED] Stress analysis on a symmetrical shaft » 2019-10-31 09:09:44

Hello,

I did not check your model, but if you want to achieve something like Remote displacement and Remote load, you should check LIAISON_SOLIDE and LIAISON_RBE3.

BR.
dezsit

#23 Re: Salome-Meca usage » Error generating mesh » 2019-10-31 09:06:48

Hello,

hex mesh: yes, for example. You can produce with this a fully structured hex mesh, but will be hard to make the connecting regions properly partitioned and to provide equal mesh seeding on the opposite edges of a block. Really a lot of work.
The solution, what I proposed on my sketch, is giving you more freedom, you can use sweeping techniques (Extrude3D in salome) which not necessarily produce structured mesh, and not necessarily produces full hex mesh (it can, but sometimes some wedges will be also generated), but it definitely will produce tetra free mesh. And with less work, and giving more freedom to change the local mesh density.

FKM: yes, that is what the referenced paper about. You can check it also here for example: vdmashop.de/Research--FKM-/

BR,
dezsit

#24 Re: Code_Aster usage » Element Catalog Location » 2019-10-31 08:52:29

Hello,

You should check the source version (eg. download from this site) and not the installed one.

BR,
dezsit

#25 Re: Salome-Meca usage » Error generating mesh » 2019-10-24 21:13:54

The suggested partitioning is not part of the model I prepared to you, actually I did not make any new partition, I just used what you created. As I mentioned, I am not able to open your small study (I do not have access to SalomeMeca for a while), and hard to explain the strategy better by words, than I did in the sketch. The key for hex meshing is to produce six sided regions (actually boxes), so you have to cut out 3 six sided regions from the sphere at first, then to make conforming mesh, you have to follow this partitioning lines in all of the connecting regions, and ensure the same number of nodes on every connecting edges.

The fatigue evaluation is a post-processing step (aster has also a dedicated post-processing command for this), so you can produce input to it with any kind of FEA code, so depending on the request we used various FEA codes to get the stresses, then made the post processing evaluation by hand calcs., spreadsheet, python scripts, or dedicated fatigue evaluation tools. But yes, nowadays we working a lot with Aster.
Prof. Fricke is a big name in the field of weld fatigue, I also recommend to carefully read out every IIW recommendations, or the FKM Guideline, ff you have access to it, or the related chapter of pressure vessel codes, (eg. EN13445), or check the available resources of FemFat (or other fatigue tools), you can learn a lot from them.

BR.
dezsit