×

Notice

The forum is in read only mode.

Why do I get 3 principal stress when solving a Plane stress (C_PLAN) problem??

  • Danilo Estay
  • Topic Author
  • Offline
  • New Member
  • New Member
More
16 years 2 days ago #4079 by Danilo Estay
This is really extrange for me, but as far as I know if the problem is plane stress I should have only 2 principal stress (I'm using C_PLAN if this is the problem please let me know), but when I extract the results to a table I get this

[code:1]
#Reac/Force at nodes
#Principal Stress
,INTITULE ,NOEUD ,RESU ,NOM_CHAM ,NUME_ORDRE ,INST ,ABSC_CURV ,COOR_X ,COOR_Y ,COOR_Z ,DX ,DY ,VMIS ,TRESCA ,PRIN_1 ,PRIN_2 ,PRIN_3 ,VMIS_SG ,TRSIG
,ReacXYZ ,N2 ,RESU ,REAC_NODA , 1 , 0.00000E+00 , 0.00000E+00 , 5.91974E-17 , 1.25767E-17 , 0.00000E+00 , 4.32859E+07 , 2.40147E+07 ,- ,- ,- ,- ,- ,- ,-
,ReacXYZ ,N1 ,RESU ,REAC_NODA , 1 , 0.00000E+00 , 1.00000E+01 , 5.00000E+00 , 8.66025E+00 , 0.00000E+00 , 1.78576E+01 ,-4.80295E+07 ,- ,- ,- ,- ,- ,- ,-
,ReacXYZ ,N3 ,RESU ,REAC_NODA , 1 , 0.00000E+00 , 2.00000E+01 , 1.00000E+01 , 0.00000E+00 , 0.00000E+00 ,-4.32859E+07 , 2.40147E+07 ,- ,- ,- ,- ,- ,- ,-
,Esfuerzos ,N137 ,RESU ,EQUI_NOEU_SIGM , 1 , 0.00000E+00 , 0.00000E+00 , 7.99715E+00 , 1.31366E+00 , 0.00000E+00 ,- ,- , 2.14783E+07 , 2.15065E+07 ,-2.14720E+07 ,-2.20640E+04 , 3.45016E+04 ,-2.14783E+07 , 0.00000E+00
,Esfuerzos ,N246 ,RESU ,EQUI_NOEU_SIGM , 1 , 0.00000E+00 , 5.50244E+00 , 2.50229E+00 , 1.60226E+00 , 0.00000E+00 ,- ,- , 1.72916E+07 , 1.74341E+07 ,-1.74341E+07 ,-2.89568E+05 , 0.00000E+00 ,-1.72916E+07 , 0.00000E+00
,Esfuerzos ,N465 ,RESU ,EQUI_NOEU_SIGM , 1 , 0.00000E+00 , 8.88197E+00 , 3.67078E+00 , 4.77336E+00 , 0.00000E+00 ,- ,- , 1.16345E+07 , 1.18679E+07 ,-1.18679E+07 ,-4.86831E+05 , 0.00000E+00 ,-1.16345E+07 , 0.00000E+00
,Esfuerzos ,N498 ,RESU ,EQUI_NOEU_SIGM , 1 , 0.00000E+00 , 1.12077E+01 , 5.49592E+00 , 6.21482E+00 , 0.00000E+00 ,- ,- , 1.98208E+07 , 1.98574E+07 ,-1.98006E+07 ,-1.65611E+04 , 5.67891E+04 ,-1.98208E+07 , 0.00000E+00
,Esfuerzos ,N508 ,RESU ,EQUI_NOEU_SIGM , 1 , 0.00000E+00 , 1.45725E+01 , 5.00000E+00 , 2.88675E+00 , 0.00000E+00 ,- ,- , 7.02033E+06 , 7.06523E+06 ,-7.06523E+06 ,-6.97455E+06 , 0.00000E+00 ,-7.02033E+06 , 0.00000E+00

[/code:1]

If I see the model with Salome postprocessor, there are also 3 principal stress, so the problem is not when I extract the values.

The problem is a triangle in the xy plane (I try other geometry but the result is the same) with a displacement imposed at each corner (to simulate a compressive force)

[code:1]
CHAR=AFFE_CHAR_MECA(MODELE=MODE,
DDL_IMPO=(_F(GROUP_NO='p1',
DX=0.001,
DY=5e-4,),
_F(GROUP_NO='p3',
DX=-0.001,
DY=5e-4,),
_F(GROUP_NO='p2',
DX=0.0,
DY=-0.001,),),);
[/code:1]

[code:1]
MODE=AFFE_MODELE(MAILLAGE=MAIL,
AFFE=_F(TOUT='OUI',
PHENOMENE='MECANIQUE',
MODELISATION='C_PLAN',),);
[/code:1]

If anybody can tell why this happend I'll appreciate it, I'm trying to aproximate this values using another method but I get large errors in the points with 3 principal stresses
  • Danilo Estay
  • Topic Author
  • Offline
  • New Member
  • New Member
More
15 years 11 months ago #4135 by Danilo Estay
Well it seems like nobody knows why this happens, anybody can tell me how to write the the stress according to xyz directions?
More
15 years 11 months ago #4141 by kwou
Hoi Danilo

I am not really sure what you mean by stresses according to xyz directions. If you mean Sxx, Syy, Xzz, ... than you can use this code snippet after MECA_STATIQUE/CALC_ELEM:

[code:1]
res=CALC_NO(reuse =resuSUP,
RESULTAT=resuSUP,
OPTION=('SIGM_NOEU_DEPL','EQUI_NOEU_SIGM',),);


TB_nodf=POST_RELEVE_T(ACTION=(_F(OPERATION='EXTRACTION',
INTITULE='local stresses',
RESULTAT=res,
NOM_CHAM='SIGM_NOEU_DEPL',
TOUT_ORDRE='OUI',
GROUP_NO='Nsupport',
NOM_CMP=('SIXX','SIYY','SIXY',),),),
TITRE='local stresses',);

IMPR_TABLE(TABLE=TB_nodf,FORMAT='TABLEAU',UNITE=26,
SEPARATEUR=' * ',TITRE='local stresses');
[/code:1]

kind regards - kees

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

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