PCOMPG
Bulk Data Entry Defines the structure and properties of a composite laminate material, allowing for global ply identification.
Format - Implicit
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
PCOMPG | PID | Z0 | NSM | SB | FT | TREF | GE | LAM | |
GPLYID1 | MID1 | T1 | THETA1 | SOUT1 | |||||
GPLYID2 | MID2 | T2 | THETA2 | SOUT2 | |||||
etc | etc | etc | etc | etc | |||||
DS | NRPT |
Format - Explicit
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
PCOMPG | PID | Z0 | NSM | SB | LAM | ||||
GPLYID1 | MID1 | T1 | THETA1 | SOUT1 | |||||
GPLYID2 | MID2 | T2 | THETA2 | SOUT2 | |||||
etc | etc | etc | etc | etc | |||||
EXPLICIT | ISOPE | HGID | NIP |
Example - Implicit
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
PCOMP | 100 | -0.5 | 1.E5 | STRN | 100. | ||||
101 | 120 | 0.2 | 0.0 | YES | |||||
2 | 120 | 0.6 | 0.0 | NO | |||||
103 | 120 | 0.2 | 0.0 | YES | |||||
1.0 |
Example - Explicit
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
PCOMP | 100 | -0.5 | 1.E5 | STRN | 100. | ||||
EXPLICIT | BWC | 100 | 5 |
Definitions
Field | Contents | SI Unit Example |
---|---|---|
PID | Unique composite property
identification.
No default (Integer > 0 or <String>) |
|
Z0 | Real number or character input
(Top/Bottom).
(Default = -0.5 * Thick, Thick being the composite total thickness (Real or blank)). Character Input 18 |
|
NSM | Nonstructural mass per unit area. No default (Real) |
|
SB | Allowable inter-laminar shear stress
(shear stress in the bonding material). Disregarded if blank or 0.0. No default (Real ≥ 0.0) |
|
FT |
Failure theory code. If blank, no failure
calculations are performed. The following failure theory codes are supported:
Default = no failure calculations are performed |
|
TREF | Reference (stress free) temperature.
2 Default = 0.0 (Real) |
|
GE | Damping coefficient. 10
11
Default = 0.0 (Real) |
|
LAM |
Laminate option. If blank, all plies must be
specified and all stiffness terms are developed. The following options are
supported:
Default = blank, that is all plies must be specified |
|
GPLYID# | Global Ply identification number. 12 No default (Integer > 0) |
|
MID# | Material IDs of individual plies. The
plies are identified by consecutively numbering them from 1 at the bottom layer. The
MIDs must refer to MAT1,
MAT2, MAT4, MAT5,
or MAT8 Bulk Data Entries. If
MID# is not specified, default is the last defined
MID#. Default = last defined MID# (Integer > 0 or blank, except that MID1 must be specified) |
|
T# | Thicknesses of individual plies. If
T# is not specified, default is the last defined
T#. Default = last defined T# (Real ≥ 0.0 or blank, except that T1 must be specified) |
|
THETA# | Orientation angle, in degrees, of the
longitudinal direction of each ply relative to the x-axis of the material coordinate
system associated with a given element. If no material coordinate system is
specified for the element, the angle is measured relative to side 1-2 of this
element. Default = 0.0 (Real or blank) |
|
SOUT# | Stress and failure index output request
for individual plies. 3
4
|
|
DS | Design switch. If non-zero (1.0), the
elements associated with this PCOMP data are included in the
topology design volume or space. Default = blank (Real = 1.0 or blank) |
|
NRPT | Number of repeat laminates 20 Default = blank (Integer > 0 or blank) |
|
EXPLICIT | Flag indicating that parameters for Explicit Analysis are to follow. | |
ISOPE | Element formulation flag for Explicit
Analysis. 21
22
23
|
|
HGID | Identification number of the hourglass
control (HOURGLS) entry. 24
25 Default = Blank (Integer > 0 |
|
NIP | Number of Gauss points through
thickness. Default = 3 (1 ≤ Integer ≤ 10) |
Comments
- String based labels allow for easier visual identification of properties, including when being referenced by other cards. (For example, the PID field of elements). For more details, refer to String Label Based Input File in the Bulk Data Input File.
- TREF specified on the PCOMPG entry overrides reference temperatures given for individual ply materials. If TREF is not specified (blank) on the PCOMPG card, then all the ply materials must have the same reference temperature.
- For SOUTi to take effect, CSTRESS must be requested in the I/O Options section of the input deck. Individual ply results will be available in addition to shell stresses and strains based on the homogenized composite properties.
- An additional piece of information available with ply results is "failure index for the element", which is the maximum of failure indices for individual plies in this element. Only the plies with SOUTi set to YES are considered in the evaluation of this maximum.
- If all plies specify zero transverse shear
coefficients (G1Z, G2Z on MAT8 card,
isotropic G for MAT1, not available for
MAT2), the in-plane shear modulus will be used to determine
transverse shear stiffness of the composite.Note: If just one layer has a nonzero value specified for transverse shear modulus, this substitution is not being performed, and user-specified values are being used for all plies.
- The signs given to stress limits for compression and tension (ST, SC, for MAT1; Xt, Xc, and so on for MAT8) are of no relevance. Absolute values are taken and used in the appropriate context to calculate failure indices.
- For composites with offset
(Z0 ≠ 0.5 * Thickness), correct values of shell stresses for the
bottom and top surfaces of the shell are produced.Note: These shell stresses are calculated using homogenized shell properties, and should be interpreted with caution.
- Element GRID thicknesses cannot be defined for elements that reference PCOMPG data.
- Plies are listed from the bottom surface upwards, in respect to the element's normal direction.
- If GE is specified on the PCOMP entry as a Real number, it will be used for the element, and the values supplied on material entries for individual plies are ignored. With USEMAT in this field, GE coefficients from ply material data will be used to calculate damping matrices for the composite. These matrices will, in general, be different for membrane, bending and shear states.
- To obtain the damping coefficient GE, multiply the critical damping ratio by 2.0.
- The global ply identification number must be unique with respect to other plies in the entry.
- For convenience, element output for the SMEAR and SMCORE options includes both homogenized shell stresses and individual ply stresses. However, because stacking sequence is ignored in these options, individual ply stresses will only be valid in cases of pure membrane deformation.
- Hill's failure theory does not differentiate between compressive and tensile strength. While different values of respective strength limits are accepted, it is still recommended that Xt is set to be equal to Xc, and Yt is set to be equal to Yc when this criteria is used. Xt and Xc are allowable tensile and compressive stresses in the principle x direction of the material. Yt and Yc are allowable tensile and compressive stresses in the principle y direction of the material.
- Failure index calculation according to Maximum Strain Theory is based on mechanical component of strain only, not on total strain. This is because only the mechanical strain contributes to actual damage of the respective ply (pure thermal expansion produces no damaging effects).
- According to the formula, some failure criteria (for example, Tsai-Wu and Hoffman) would produce a negative ply failure, depending on the problem.
- If PARAM, SRCOMPS,YES is added to the input file, strength ratios with respect to designated failure theory are output for composite elements that have failure indices requested.
- The following two formats are permissible for
the Z0 field:
Real Number:
It represents the distance from the shell element reference plane to the bottom surface of the shell (Default = -0.5 * Thick, Thick being the composite total thickness (Real or blank)).
Surface:
Top:
The shell reference plane, the plane defined by the grid points, and the top surface of the shell are coplanar.
This makes the effective "Real" Z0 value equal to the composite total thickness (-1.0 * Thick). See Figure 2.Bottom:
The shell reference plane, the plane defined by the grid points, and the bottom surface of the shell are coplanar.
This makes the effective "Real" Z0 value equal to 0. See Figure 3.Automatic offset control is available for ply thickness (size) optimization and for free-size optimization, where the specified offset values are automatically updated based on thickness changes. For free-size optimization, such an automatic offset, is only applicable when Z0=0.0 or BOTTOM.
- The material parameters, Xt, Xc, Yt, Yc, and S on the MAT8 Bulk Data Entry should be specified for failure criteria calculation.
- The repeat laminates are added to the bottom of the current laminate. Output for repeat laminates are supported in OPTI, H3D, and PUNCH formats.
- For CTRIA3 elements in explicit analysis, triangular shell formulation is automatically applied. Therefore, the definition of ISOPE has no effect on CTRIA3 elements in explicit analysis.
- Both Belytschko-Tsay and Belytschko-Wong-Chiang shell formulations are very effective and robust. However, the performance of Belytschko-Tsay is poor, if the elements are warped. With Belytschko-Wong-Chiang formulation, the limitations in the element warpage are fixed with 20-30% additional computational cost.
- The Belytschko-Tsay and Belytschko-Wong-Chiang shells do not possess stiffness in the normal rotational degree of freedom and this would lead to a singular stiffness matrix in case of implicit analysis. In explicit analysis, the unconstrained drilling degree of freedom usually does not create any difficulties since a stiffness matrix is not involved.
- For four-noded quadrilateral elements in explicit analysis, hourglass control is required to avoid spurious zero-energy modes. Triangular elements do not require hourglass control.
- When HGID is not specified, a default hourglass control is used.
- This card is represented as a property in HyperMesh.