×

Notice

The forum is in read only mode.

Syntax problem: Apply rotation

More
15 years 4 months ago #3707 by kwou
Replied by kwou on topic Re:Syntax problem: Apply rotation
Hi Jakob

still I cannot open your zip file: also the site states it is 0 bytes. So I cannot compare your comm file to mine.

Indeed I ran the comm file here, and it produces results, with no (critical) errors on the gravity and rotation part. My version of C-Aster is 9.4, ASTK is 1.7.3. There may be small differences between higher and lower versions of CA. And indeed, the CENTRE keyword may be omitted without any side effects.

And since you are a bit reluctant to believe me :dry: , I attach a picture of the post pro (I think it is the rotation load being displayed).

kind regards - kees

<br /><br />Post edited by: Kees Wouters, at: 2010/01/05 00:29

Interest: structural mechanics, solar energy (picture at 'my location' shows too little pv panels)

--
kind regards - kees
Attachments:
  • Jakob Andersen
  • Topic Author
  • Offline
  • Junior Member
  • Junior Member
More
15 years 4 months ago #3711 by Jakob Andersen
Replied by Jakob Andersen on topic Re:Syntax problem: Apply rotation
Hi Kees,

If all else fails: ctrl+c, ctrl+v :)

Here is my message file:

[code:1]==========================================
==========================================
DEBUT();

MA=DEFI_MATERIAU(ELAS=_F(E=2e11,
NU=0.3,
RHO=7800,),);



MAIL=LIRE_MAILLAGE(FORMAT='MED',);

MODE=AFFE_MODELE(MAILLAGE=MAIL,
AFFE=_F(TOUT='OUI',
PHENOMENE='MECANIQUE',
MODELISATION= '3D',),);



MAIL=MODI_MAILLAGE(reuse =MAIL,
MAILLAGE=MAIL,
ORIE_PEAU_3D=_F( GROUP_MA=('profile1',),),
);

MATE=AFFE_MATERIAU(MAILLAGE=MAIL,
AFFE=_F(TOUT='OUI',
MATER=MA,),);

CHAR=AFFE_CHAR_MECA(MODELE=MODE,
DDL_IMPO=(
_F(GROUP_MA='b1',
DX=0,
DY=0,
DZ=0,),),
PRES_REP=(
_F(GROUP_MA='profile1',
PRES=0,),),);

GRAV=AFFE_CHAR_MECA(MODELE=MODE,

PESANTEUR=(9810,0,0,-1),);


ROT=AFFE_CHAR_MECA(MODELE=MODE,

ROTATION=(37.7,0,1,0,),

CENTRE = (0,0,0),);



RESU=MECA_STATIQUE(MODELE=MODE,
CHAM_MATER=MATE,
EXCIT=(
_F(CHARGE=CHAR,),
_F(CHARGE=GRAV,),
_F(CHARGE=ROT,),),);

RESU=CALC_ELEM(reuse =RESU,
MODELE=MODE,
CHAM_MATER=MATE,
RESULTAT=RESU,
OPTION=('SIGM_ELNO_DEPL','EQUI_ELNO_SIGM',),
EXCIT=_F(
CHARGE=CHAR,),);

RESU=CALC_NO(reuse =RESU,
RESULTAT=RESU,
OPTION=('SIGM_NOEU_DEPL', 'EQUI_NOEU_SIGM', ),);

IMPR_RESU(FORMAT='MED',
UNITE=80,
RESU=_F(MAILLAGE=MAIL,
RESULTAT=RESU,
NOM_CHAM=('SIGM_NOEU_DEPL','EQUI_NOEU_SIGM','DEPL',),),);

FIN();

==========================================
==========================================
JDC.py : ERREUR DE COMPILATION DANS ACCAS - INTERRUPTION
&gt;&gt; JDC.py : DEBUT RAPPORT
CR phase d'initialisation
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Compilation impossible : File &quot;fort.1&quot;, line 37 !
! GRAV=AFFE_CHAR_MECA(MODELE=MODE, !
!
^ !
! SyntaxError: invalid syntax !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
fin CR phase d'initialisation

&gt;&gt; JDC.py : FIN RAPPORT
EXECUTION_CODE_ASTER_EXIT_13706=1

&lt;F&gt;_ABNORMAL_ABORT exit code = 1[/code:1]

I originally created the .comm file using Salomé, which means it should be in Aster 9.4. My ASTK is 1.7.1 - all standard CAELinux...

Thanks again!

Jakob
More
15 years 4 months ago #3712 by kwou
Replied by kwou on topic Re:Syntax problem: Apply rotation
Hi Jakob

I used your code (copy, paste and change MODELE=MOD--&gt;Bmod, including open lines):
[code:1]
ROT=AFFE_CHAR_MECA(MODELE=MODE,

ROTATION=(37.7,0,1,0,),

CENTRE = (0,0,0),);
[/code:1]

The result in the mess file, no problems reported:

[code:1]
ROT=AFFE_CHAR_MECA(INFO=1,
CENTRE=(0.0,0.0,0.0,),
LIAISON_XFEM='NON',
VERI_NORM='OUI',
ROTATION=(37.700000000000003,0.0,1.0,0.0,),
MODELE=Bmod,
);
[/code:1]

Did you copy your code between windows and linux systems?
Then, if you use a dos file in a linux environment, this causes problems with carriage return codes. I use the dos2unix tool, because of the different EOL symbols on windows and linux. So, if you put the command on one line, does it help (and make sure you make a new &lt;carriage return&gt; at the end of line in the file)?

ROT=AFFE_CHAR_MECA(MODELE=MODE,ROTATION=(37.7,0,1,0,),CENTRE = (0,0,0),); &lt;new carriage return&gt;

When I make my command file a ''dos'' file [todos &lt;command&gt;.comm] I get the following error, looking very much like yours:
[code:1]
JDC.py : ERREUR DE COMPILATION DANS ACCAS - INTERRUPTION
&gt;&gt; JDC.py : DEBUT RAPPORT
CR phase d'initialisation
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Compilation impossible : File &quot;fort.1&quot;, line 1 !
! DEBUT(); !
! ^ !
! SyntaxError: invalid syntax !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
fin CR phase d'initialisation

&gt;&gt; JDC.py : FIN RAPPORT
EXECUTION_CODE_ASTER_EXIT_7044=1

&lt;F&gt;_ABNORMAL_ABORT exit code = 1
[/code:1]

One other remark: since you use ISO measures, you need to change g to 9.81 [m/s2] instead of 9810 [mm/s2]. And you mesh cq dimensions need to be in [m] of course. Just to make sure your results make sense.

kind regards - kees<br /><br />Post edited by: Kees Wouters, at: 2010/01/05 15:55

Interest: structural mechanics, solar energy (picture at 'my location' shows too little pv panels)

--
kind regards - kees
  • Jakob Andersen
  • Topic Author
  • Offline
  • Junior Member
  • Junior Member
More
15 years 4 months ago #3713 by Jakob Andersen
Replied by Jakob Andersen on topic Re:Syntax problem: Apply rotation
Hi again Kees,

Due to driver problems (solved as of today) I have been forced to download everything on my windows laptop and save to a shared disk, so you pretty much hit the diagnose spot on!

I was not familiar with this problem. I even tried to use gedit in windows to make sure, that compatibility was ok.

I will try to implement the changes tonight.

You've been a great help!

Thanks a lot!

Best regards,
Jakob
  • Jakob Andersen
  • Topic Author
  • Offline
  • Junior Member
  • Junior Member
More
15 years 4 months ago #3715 by Jakob Andersen
Replied by Jakob Andersen on topic Re:Syntax problem: Apply rotation
Hi Kees,

Indeed it worked!

A new problem has appeared though :(

[code:1]!
!
! &lt;EXCEPTION&gt; &lt;FACTOR_20&gt; !
! !
! !
! -&gt; Matrice non factorisable : !
! Le pivot est presque nul à la ligne 115402 pour le noeud N38468 et !
! la composante DX. !
! Pour information, le nombre de décimales perdues est de 11. !
! !
! -&gt; Conseil &amp; Risque : !
! Il s'agit peut etre d'un mouvement de corps rigide mal bloqué. !
! Vérifiez les conditions aux limites. !
! Si vous faites du contact, il ne faut pas que la !
! structure ne &quot;tienne&quot; que par le contact. !
! !
! !
! !
! !
!
![/code:1]

I have imposed a DX=DY=DZ=0 condition on the 'b1' boundary - similar to what you did in your .comm file.

To make things more specific I have done a screen dump from Salome of the model - a single helically shaped NACA0018 profile. The 'b1' boundary consists of several groups of faces along the profile as can be seen highlighted in white.

I hate to take any more of your time, but do you have any ideas on how to solve this problem?

BTW: the .comm file is attached - hopefully this will work, now that i am posting from CAELinux...

Best regards,

Jakob

Attachment WNGASTK.zip not found

Attachments:
More
15 years 4 months ago #3716 by kwou
Replied by kwou on topic Re:Syntax problem: Apply rotation
Hi Jakob

Nice to hear it runs now.

And regarding the singular matrix , you not only have to define the boundary conditions ''CHAR'' in your comm file, but also apply them in the ''MECA_STATIQUE'', see snippet below.

[code:1]
RESU=MECA_STATIQUE(MODELE=Bmod,
CHAM_MATER=MATE,
EXCIT=(_F(CHARGE=CHAR,), # add _F(CHARGE=CHAR,),
_F(CHARGE=gravity,),
_F(CHARGE=rotation,),),);
[/code:1]


[Though you applied it in the ''CALC_ELEM'' part, I am not sure what the effect of this keyword is here. I *think* you normally can omit the CHARGE keyword here.]

If you can run it successfully, could you see whether the centre of rotation has any effect of the results? I am still puzzled with that, since I do not see any difference when changing the value of this keyword. Thank you in advance.

And yes, the comm/zip file is downloadable now :dry:.

Do not worry about my time, if I am running out of it you will notice :P.

kind regards - kees<br /><br />Post edited by: Kees Wouters, at: 2010/01/06 22:38

Interest: structural mechanics, solar energy (picture at 'my location' shows too little pv panels)

--
kind regards - kees
Moderators: catux
Time to create page: 0.140 seconds
Powered by Kunena Forum