#1 2020-11-16 14:02:20

JanBlokes

Member

Registered: 2018-09-06

Posts: 4

Composite failure criterions

Hello,

I'm Jan from the Czech Republic and I'm interested in CA and Python programming.

The more and more people are taking interest in composites and I haven't found any evaluation for composites like Tsai-Wu, Hill or Hoffman criterions. Has anybody tried to programme it in Code_Aster? I would like to have the results in Paravis.

Maybe it's my fault and can't search properly

Thank you for any information!

Jan

mecour

Member

From: Ostrava (Czech)

Registered: 2011-04-04

Posts: 160

Re: Composite failure criterions

Hello,

there is no such criterion in code-aster, you can just use orthotropic elastic material model or you can use the mfront and prepare your own material model.

mecour

JanBlokes

Member

Registered: 2018-09-06

Posts: 4

Re: Composite failure criterions

Thank you for response. I'll try.

bpaillard

Member

Registered: 2015-10-16

Posts: 25

Re: Composite failure criterions

Hi,

I'm doing it with a bunch of FORMULE, applied to SIGM_ELGA SIXX, SIYY and SIXY, with the proper fiber angle, and then calculating it with CALC_CHAMP.

JanBlokes

Member

Registered: 2018-09-06

Posts: 4

Re: Composite failure criterions

Hi bpaillard,

is it possible to share .comm file, please?

Thank you,
Jan

bpaillard

Member

Registered: 2015-10-16

Posts: 25

Re: Composite failure criterions

Yes sorry for that, I should have done it earlier. Here are parts of the file.

We have validated a cylinder case with a steel / composites material distribution, on Von Mises and Tsai Wu. If Code_aster developers want to include it in the distribution please feel free to contact me, it should be doable.

Regards.

Xc=..
Yc=..
Xt=..
Yt=..
for ii in range(len(plis)):
.....
    FO[ii] = FORMULE(
VALE=    '(SIXX*cos(teta)**2+SIYY*sin(teta)**2 - SIXY*2*cos(teta)*sin(teta) )*(1/'+str(Xt)+'-1/'+str(Xc)+') +  \
     (SIXX*sin(teta)**2+SIYY*cos(teta)**2 + SIXY*2*cos(teta)*sin(teta) )*(1/'+str(Yt)+'-1/'+str(Yc)+') + \
     (SIXX*cos(teta)**2+SIYY*sin(teta)**2 - SIXY*2*cos(teta)*sin(teta) )**2/('+str(Xt)+'*'+str(Xc)+') + \
     (SIXX*sin(teta)**2+SIYY*cos(teta)**2 + SIXY*2*cos(teta)*sin(teta) )**2/('+str(Yt)+'*'+str(Yc)+') +\
     ( (SIXX*sin(teta)*cos(teta)-SIYY*sin(teta)*cos(teta) + SIXY*(cos(teta)**2-sin(teta)**2) ) /'+str(S)+')**2',
NOM_PARA=('SIXX','SIYY','SIXY',), teta=-PROP_2[ii][1]*pi/180.)

    PLI1P[ii]=CALC_CHAMP(RESULTAT=plii[ii],
                 GROUP_MA='PROP_2',
                 CHAM_UTIL=_F(
                    NOM_CHAM='SIGM_ELGA',
                              FORMULE=FO[ii],
                              NUME_CHAM_RESU=1,
),);
.....

Formulas can be for instance found in D. Gay book 'matériaux composites'

JanBlokes

Member

Registered: 2018-09-06

Posts: 4

Re: Composite failure criterions

Thank you very much, it is very useful!