Abaqus Tutorial 11b PlyFailure
Abaqus Tutorial 11b PlyFailure
Composites
Modelling ply failure
Stephanie Miot
First matrix
crack
Fibre failure
initiation
1.
Introduction
In this tutorial, you will modify a structural model of a stiffened panel to define the
material properties including the ply failure parameters. You will then perform a
static analysis of a bending test and visualize the simulation of the damage
propagation with Abaqus/Viewer.
Abaqus/Viewer
Use the visualization module to create ply stack plots and contour plots on
different plies
Preliminaries
The stiffened panel is composed of:
-
a skin
Dimensions: 1600
16 mm x 1000 m, 6 mm thick
2 stiffeners
Dimensions: 200
20 mm x 140 mm, 1600 mm
m long, 4 mm thick
2.
the properties for each material and include the Hashins failure criteria
and
nd the progressive damage models
the mesh
Finally, you will run a static analysis and use the visualization module to postpost
process the results of the simulation.
3.
1. Define the mechanical behaviour of the UD plies with the following material
properties:
T300/M18
T800/M18
Elasticity coefficients
E1
170 GPa
192
92 GPa
E2
9 GPa
15 GPa
12
0.34
0.3
0.34
G12
4.8 GPa
5.6 GPa
G13
4.8 GPa
5.6 GPa
G23
4.5 GPa
4.9 GPa
2050 MPa
2240
0 MPa
XC
1200 MPa
1230
0 MPa
YT
62 MPa
71 MPa
YC
190 MPa
210
0 MPa
SL
81 MPa
90 MPa
ST
81 MPa
90 MPa
Fracture toughnesses
Long. tensile fracture
energy: Gftc
95 kJ/m
105
5 kJ/m
kJ/m
Long. compressive
fracture energy: Gfcc
103 kJ/m
108 kJ/m
0.2 kJ/m
0.2 kJ/m
Trans. compressive
fracture energy: Gmcc
0.2 kJ/m
0.2 kJ/m
a. Go into the Property Module and click the Create Material icon
f. Click Suboptions and select Damage Stabilization. Specify the value of the
viscosity coefficient for each failure mode: 1e-4.
1e
Name
the
new
lay
lay-up
definition: Skin lay-up.. Set the Initial ply count at 6 and select the Element
Type: Continuum Shell.
Shell Click Continue...
b. In the Edit Composite Layup dialog box, define the Layup orientation. Select
Definition: Coordinate system.
system Click the Select CSYS icon
to select
c. Accept the default selection for the Normal direction: Axis 3 and the
Stacking Direction: Element direction 3.
3
i. Use the options available in the Display tab to check the orientation of the
lay-up. Then click OK to create the new lay-up.
web, 32-ply region, (02, -452, 902, 452, 452, 902, -452, 02)S
4.
Mesh
Element Type dialog box, select Family: Continuum Shell and accept the
other default selections.
. Select Technique:
Technique Sweep and
f. Check that the stacking direction is correct. Click the Assign Stack Direction
icon
and select the top face of the skin to define the reference
orientation.
d. In the menu bar, select Object: Assembly to visualise the mesh of the
different instances.
10
5.
Static analysis
1. Create a new analysis step and define the parameters of the analysis.
analysis
a. Go into the Step Module.
Module Click the Create Step icon
c. In the Edit Step dialog box, toggle on Nlgeom. In the Incrementation tab,
define the Maximum
aximum number of increments: 100 000 and the Increment
size: Initial: 0.01, Minimum: 1e-9, Maximum: 0.1. Click OK.
11
2. Create new output requests. The preselected default output does not include
the damage variables and the default selection of the layered section points
includes only the top and bottom points.
points To visualize the stresses and the
damage evolution in each ply in the Visualization module, you will write additional
field output to the output database file.
file
a. Click the Create Field Output icon
b. In the Edit Field Output Request dialog box, select Domain: Set: Skin-1.All.
Then select Frequency: Evenly spaced time intervals and set Interval: 10.
c. In the Output Variables list, expand the Stresses list and toggle on the
variable S,, then expand the Failure/Fracture list and toggle on the variables
DAMAGEFT,
DAMAGEF
DAMAGEFC,
DAMAGEMT,
HSNFCCRT,
HSNMTCRT
HSNMTCRT,
and
DAMAGEMC,,
HSNMCCRT.
Finally,
HSNFTCRT,
expand
the
12
13
6.
Boundary conditions
14
b. Select
Step:
Initial
Initial,
Category:
Mechanical
and
Type:
c. Select the region from the list of eligible Sets: BC-1.. Toggle on U1, U2 and
U3. Click OK.
15
e. Select the reference point RP-1 then specify the applied displacements:
U2 = 14 mm and U3 = -70
- mm. Click OK.
16
7.
Analysis
nalysis and post-processing
post
c. In the Edit Job dialog box, open the tab: Parallelization.. Toggle On Use
multiple processors and set the number of processors at 4.. Accept the other
default selections and click OK. Then submit the job: Panel01.
d. Click Monitor... to monitor the job while it is running.
17
a. When the job submission has been completed, in the Job Manager dialog
box, click Results or open Panel01.odb in the Visualization Module.
Module
Create...
c. In the Create Display Group dialog box, select Part instances: SKIN-1.
Click the Replace icon
f. In the Field Output dialog box, in the Primary Variable tab, select the
Output Variable: HSNFTCRT.
HSNFTCRT Click Apply then click Section Points..
g. In the Section Points dialog box, click Selection method: Plies and select
the ply Skin - ply 0 - 1.. Click Apply.
You can then visualise the fibre failure prediction in the 0 ply.
1
2
18
h. In the Field Output dialog box, select the Output Variable: DAMAGEFT to
visualise the fibre failure propagation. Click Apply. In the Section Points
dialog box, click Selection method: Plies and select the ply Skin - ply 0 - 1.
Click Apply.
i. In the Field Output dialog box, select the Output Variable: DAMAGEMT
DAMAGE
to
visualise the matrix failure propagation. In the Section Points dialog box,
select the ply Skin - ply 90-1. Then select Skin - ply -45.. You can observe
the propagation of the matrix damage in the different plies.
plies
2
1
1
2
19
l. In the Field Output dialog box, select the Output Variable: DAMAGEFT.
DAMAGE
In
the Section Points dialog box, select the ply SRO - ply 0. Click Apply.
Select various variables and section points to plot the numerical predictions
for the damage initiation and propagation in different plies.
20
m. Click Tools XY Data Create. In the Create XY Data dialog box, select
ODB field output.. Click Continue...
n. In the XY Data from ODB Field Output dialog box, in the Variables tab,
select Position: Unique Nodal.
Nodal In the variables list, expand RF: Reaction
force and toggle on Magnitude,
Magnitude then expand U: Spatial displacement and
toggle on Magnitude.. In the Elements/Nodes tab, select Method: Node
sets: ASSEMBLY_CONSTRAINT-2_REFERENCE_POINT.
ASSEMBLY_CONSTRAINT
. Click Save. A
warning message is displayed. Click OK. Then click Dismiss.
21
o. Click Tools XY Data Create. In the Create XY Data dialog box, select
Operate on XY data.. Click Continue...
r. Double-click
click the chart. In the Chart Options dialog box, toggle off Fill then
click Dismiss. Double-click
click the X axis. In the Axis Options dialog box, in the
Title tab, click the Font icon. Set the size at 24 and click OK.. In the Axes tab,
click the Font icon
on and set the size at 18.
22
23