calculation of Stress in composites

[ho1mo]

Hi ,

I am Homo and beginner in fiber composites with ABAQUS, but I have a problem,

I can't calculate stress and other element variables, can you help me ?

Too,I have to use ABAQUS/CAE for model and here is my *.inp file.

Thanks for your help

Homo

[Jorgen]

I don't quite understand your question. Do don't need to calculate the stress, Abaqus will do it for you ...

-Jorgen

[ho1mo]

Hi
Nice to meet you Mr.Jorgen,
I have simulated a solid composite model in ABAQUS, but I don't know how I can request stresses , strains and any elemental variables in it.

Can you help me?
I can send you *.inp file if requested.

Thanks for your kindly reply
Homo

[Jorgen]

Ah, I see. You need to specify to Abaqus that you want the stresses, strain, ets saved. You can do this in either CAE or by manually editing the inp-file.

Here is an example of how you can do that for an Abaqus/Standard analysis:

*Output, field, variable=PRESELECT

- Jorgen

[ho1mo]

Thanks a lot Dr.Jorgen,
But I think should be do another, I don't know.
This is my *.inp file, can you check my inp file?

Best regards
Homo
*HEADING
INFINITELY LONG CYLINDRICAL TUBE SUBJECTED TO
INTERNAL PRESSURE
UNITS: N, mm
*RESTART,WRITE
*NODE
1, 60.0, 0.0
9, 0.0, 60.0
161, 140.0, 0.0
169, 0.0, 140.0
*NGEN,LINE=C,NSET=Di
1, 9, 1, , 0.0, 0.0, 0.0
*NGEN,LINE=C,NSET=Do
161, 169, 1, , 0.0, 0.0, 0.0
*NFILL,NSET=Face1
Di, Do, 16, 10
*NCOPY,CHANGE NUMBER=2000,OLD SET=Face1,SHIFT,NEW SET=Face2
0.0, 0.0, 300.0
0.0,
*NFILL,NSET=Nall
Face1, Face2, 2, 1000
*NSET,GENERATE,NSET=Xsymm
1, 161, 10
1001, 1161, 10
2001, 2161, 10
*NSET,GENERATE,NSET=Ysymm
9, 169, 10
1009, 1169, 10
2009, 2169, 10
*NSET,GENERATE,NSET=Noutp
1, 161, 10
**
*ELEMENT,TYPE=C3D20R
1, 1, 41, 43, 3, 2001, 2041, 2043, 2003, 21, 42, 23, 2, 2021, 2042,
2023, 2002, 1001, 1041, 1043, 1003
11, 41, 161, 163, 43, 2041, 2161, 2163, 2043, 101, 162, 103, 42,
2101, 2162, 2103, 2042, 1041, 1161, 1163, 1043
*ELGEN,ELSET=Ei
1, 4, 2, 1
*ELGEN,ELSET=Eo
11, 4, 2, 1
*ELSET,GENERATE,ELSET=Inside
1, 4
*ELSET,ELSET=Eioutp
1,
*ELSET,ELSET=Eooutp
11 ,
*SOLID SECTION,ELSET=Ei,COMPOSITE,ORIENTATION=Ocircum,STA CKDIR=1
1.0, 3, Mcircum
1.0, 3, Maxial
*SOLID SECTION,ELSET=Eo,COMPOSITE,ORIENTATION=Ocircum,STA CKDIR=1
1.0, 3, Mcircum
1.0, 3, Maxial
1.0, 3, Mcircum
1.0, 3, Maxial
1.0, 3, Mcircum
1.0, 3, Maxial
*MATERIAL,NAME=Mcircum
*ELASTIC,TYPE=ENGINEERING CONSTANTS
10.0E3, 250.0E3, 10.0E3, 0.01, 0.25, 0.25, 5.0E3, 2.0E3,
5.0E3 ,
*MATERIAL,NAME=Maxial
*ELASTIC,TYPE=ENGINEERING CONSTANTS
10.0E3, 10.0E3, 250.0E3, 0.25, 0.01, 0.01, 2.0E3, 5.0E3,
5.0E3 ,
*ORIENTATION,NAME=Ocircum,SYSTEM=CYLINDRICAL
0.0, 0.0, 0.0, 0.0, 0.0, 30.0
1, 0.0
*BOUNDARY
Xsymm, 2
Ysymm, 1
Face1, 3
Face2, 3
*******************
*STEP,PERTURBATION
*STATIC
*DLOAD
Inside, P6, 50.0
*EL PRINT,ELSET=Eioutp
1,2,3,4,5,6
COORD, S11, S22, S33
*EL PRINT,ELSET=Eooutp
1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16
COORD, S11, S22, S33
*EL PRINT,ELSET=Eooutp
17,18
COORD, S11, S22, S33
*EL FILE
S, E
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,E
*NODE PRINT,NSET=Noutp
U1,
*NODE FILE
U,RF
*OUTPUT,FIELD
*NODE OUTPUT
U,RF
*END STEP

[Jorgen]

Looks OK during a quick glanze.
I suggest that you give it a try and debug it using the Abaqus manuals - which are excellent.

- Jorgen