How calculate the Principal Stresses on a surface of a solid?

Johannes_ACKVA · 2019-05-27 09:53:53

Bonjour,

does anybody know about a way how to calculate the 2 principal stresses on a surface of a solid?

Assuming that no pressure acts on the solid it is ensured that 2 (of the 3 principal stresses) are in the tangent plane of the 3D-solid. In some standards, for example the german FKM-Richtlinie those 2 tangent principal stresses are needed for stress evaluation. I looked to CALC_CHAMP, but it seems none of the options does that

thank you for your any help

Regards
Johannes_ACKVA

_________________________________________________________
CODE-ASTER-courses at Ingenieurbüro für Mechanik, Germany

*** CODE-ASTER INTRODUCTION
14 - 18 Oct 2019

*** CODE-ASTER DYNAMIC ANALYSIS
21 - 22 Nov 2019

Volker · 2019-05-27 11:47:54

Bonjour Johannes,
which part of FKM Guideline do you mean ??
Do you mean for example calculation of G_Sigma and G_Tau according equation no. 4.3.17 ??

Kind regards Volkere

jeanpierreaubry · 2019-05-27 16:10:47

hello Johannes

if you do a CALC_CHAMP on some group you are to have in the result

'VMIS',
'PRIN_1','PRIN_2','PRIN_3', the principal stress values (modulus of the vector)
'VECT_1_X','VECT_1_Y','VECT_1_Z',
'VECT_2_X','VECT_2_Y','VECT_2_Z',
'VECT_3_X','VECT_3_Y','VECT_3_Z',
the component along X, Y, Z global of each principal stress

after that anything is possible with some post processing in code_aster or Gmsh

the real difficulty is to eliminate the one perpendicular to the plane of the element
i have been struggling at that for quite a while and never found a good solution

however beware that principal directions and principal stress may be meaningless once we do averaging at nodes in SIEQ_ELNO, SIEQ_NOEU
code-aster.org/forum2/viewtopic.php?id=16737
it may be strictly true only with ELGA field
and as far as i know in 3D modeling the Gauss point are not on the face of the elements

i cannot at the time being give a more valuable information
but this task is on my list for 2D elements

the attached screen shot shows SIEQ_NOEU as a vector on a shell model

jean pierre aubry

paster · 2019-05-28 08:22:16

Hello

Maybe you can try to put skin elements around your model. Of course thin thickness, low Young modulus, same nodes .
after that this is Aster "cooking" (only 2 principal stresses by node, element…): I'm not a Code_Aster's specialist (may be Jean Pierre?). This method is quite general and simple

Good luck

jeanpierreaubry · 2019-05-28 08:59:27

Maybe you can try to put skin elements around your model. Of course thin thickness, low Young modulus, same nodes .

i do think this is not necessary maybe giving wrong results by polluting the model
as post processing on a group of surface elements defined on the skin does the job
at least this is my experience

Volker · 2019-05-28 10:04:30

Hi All,
here a sketch with my interpretation of the local stresses according the FKM Guideline:

How can we calculate the stress gradient with Code_Aster??

Kind regards Volker
PS

jeanpierreaubry wrote:

i do think this is not necessary maybe giving wrong results by polluting the model..

I total agree with you ,
but I got also the same advice from a commercial support-team to introduce skin elements ;-)

Attachments:
FKM_stresses.jpg, Size: 47.57 KiB, Downloads: 372

jeanpierreaubry · 2019-05-28 10:28:43

Volker wrote:

but I got also the same advice from a commercial support-team to introduce skin elements ;-)

probably because their code cannot handle stress on face boundary

paster · 2019-05-28 11:24:35

Re Hello

Open your eyes. Commercial codes are not the devil! This method is not exect of course but better than any kind of extrapolation.
There is no pollution if you take care to the stiffness (lower than possible by thickness) of the skin membrane elements (same dof's like volume elements). Ex: QUA8 on HEX20, TRI3 on TET4.... (for 3D analysis)

Sincerely

jeanpierreaubry · 2019-05-28 12:57:31

paster wrote:

Re Hello
Open your eyes. Commercial codes are not the devil! This method is not exect of course but better than any kind of extrapolation.

this is not the question
code_aster does what is wanted without this !

Johannes_ACKVA · 2019-05-28 19:21:39

@paster,

skin elements for stress evaluation: with some FE-codes this is the suitable method, for example Nastran. Nastran has pure membran elements which have only the translational DOFs, so those as skin elements do not augment the number of DOFs. But Code-Aster's 2D-structural elements have all the rotational DOFs in addition. When making them very thin, there can be numerical problems with the rot DOFs.

@Volker

at the moment I haven't thoroughly studied the FKM-Guideline, I m preparing only the first steps.
For the derivative of stresses something like this (very rough) is perhaps suitable:

CALC_CHAMP( SIGM_NOEU, SIEQ_NOEU)

MACR_LIGN_COUPE to get the evolution of stresses through the thickness of a solid and put it in a table. You can collect components of SIGM_* of of SIEQ_*. Which means that you can collect also the principal stresses

RECU_FONCTION to make a function: Stress(ThicknCoor) or Principal(ThicknCoor)

CALC_FONCTION ( DERIVE..) to make the derivative of the function

When treating with princ stresses I see the problem that the directions of the principal stresses change with the abszissa, do you know a solution to that?

Regards
Johannes_ACKVA

_________________________________________________________
CODE-ASTER-courses at Ingenieurbüro für Mechanik, Germany

*** CODE-ASTER INTRODUCTION
14 - 18 Oct 2019

*** CODE-ASTER DYNAMIC ANALYSIS
21 - 22 Nov 2019

paster · 2019-05-30 18:04:21

Hello

One more I'm not a code aster guru but, elements like MEMBTR3, MEMBQU4... are not good for your application (they have 3 translational dof's). I don't know what meshing tool you use but even with Gibi, it's very easy to extract the envelop of a 3D model

After that end for me

AsterO'dactyle · 2019-05-31 21:17:21

Hello,

Try SIRO_ELEM in CALC_CHAMP

Johannes_ACKVA · 2019-06-03 22:31:27

AsterO'dactyle wrote:

Try SIRO_ELEM in CALC_CHAMP

that's it ! Many thanks, AsterO'dactyle !

Regards
Johannes_ACKVA

_________________________________________________________
CODE-ASTER-courses at Ingenieurbüro für Mechanik, Germany

*** CODE-ASTER INTRODUCTION
14 - 18 Oct 2019

*** CODE-ASTER DYNAMIC ANALYSIS
21 - 22 Nov 2019

#1 2019-05-27 09:53:53

How calculate the Principal Stresses on a surface of a solid?

#2 2019-05-27 11:47:54

Re: How calculate the Principal Stresses on a surface of a solid?

#3 2019-05-27 16:10:47

Re: How calculate the Principal Stresses on a surface of a solid?

#4 2019-05-28 08:22:16

Re: How calculate the Principal Stresses on a surface of a solid?

#5 2019-05-28 08:59:27

Re: How calculate the Principal Stresses on a surface of a solid?

#6 2019-05-28 10:04:30

Re: How calculate the Principal Stresses on a surface of a solid?

#7 2019-05-28 10:28:43

Re: How calculate the Principal Stresses on a surface of a solid?

#8 2019-05-28 11:24:35

Re: How calculate the Principal Stresses on a surface of a solid?

#9 2019-05-28 12:57:31

Re: How calculate the Principal Stresses on a surface of a solid?

#10 2019-05-28 19:21:39

Re: How calculate the Principal Stresses on a surface of a solid?

#11 2019-05-30 18:04:21

Re: How calculate the Principal Stresses on a surface of a solid?

#12 2019-05-31 21:17:21

Re: How calculate the Principal Stresses on a surface of a solid?

#13 2019-06-03 22:31:27

Re: How calculate the Principal Stresses on a surface of a solid?

Board footer