Elemento Tipo Shell 181 de ANSYS
Elemento Tipo Shell 181 de ANSYS
Elemento Tipo Shell 181 de ANSYS
Page: 1
SHELL181
4-Node Structural Shell
MP ME ST PR PRN DS DSS <> <> <> <> PP VT EME MFS
Product Restrictions
Element Reference
Page: 2
Single-Layer Definition
To define the thickness (and other information), use section definition, as follows:
SECTYPE,,SHELL
SECDATA,THICKNESS, ...
A single-layer shell section definition provides flexible options. For example, you can specify the number of integration
points used and the material orientation.
Multilayer Definition
The shell section commands allow for layered shell definition. Options are available for specifying the thickness, material,
orientation, and number of integration points through the thickness of the layers.
You can designate the number of integration points (1, 3, 5, 7, or 9) located through the thickness of each layer when
using section input. When only one, the point is always located midway between the top and bottom surfaces. If three
or more points, two points are located on the top and bottom surfaces respectively and the remaining points are
distributed equal distance between the two points. The default number of integration points for each layer is three;
however, when a single layer is defined and plasticity is present, the number of integration points is changed to a
minimum of five during solution.
The following additional capabilities are available when defining shell layers:
SHELL181 accepts the preintegrated shell section type (SECTYPE,,GENS).
When the element is associated with the GENS section type, thickness or material definitions are not required.
You can use the function tool to define thickness as a function of global/local coordinates or node numbers
(SECFUNCTION).
You can specify offsets (SECOFFSET).
A section can be partially defined using data from a FiberSIM .xml file (SECTYPE,,SHELL,FIBERSIM).
Other Input
The default orientation for this element has the S1 (shell surface coordinate) axis aligned with the first parametric
direction of the element at the center of the element, which connects the midsides of edges LI and JK and is shown as
xo in Figure 181.1. In the most general case, the axis can be defined as:
where:
Element Reference
Page: 3
The element supports degeneration into a triangular form; however, use of the triangular form is not recommended,
except when used as mesh filler elements or with the membrane option (KEYOPT(1) = 1). The triangle form is generally
more robust when using the membrane option with large deflections.
To evaluate stresses and strains on exterior surfaces, use KEYOPT(1) = 2. When used as overlaid elements on the faces
of 3-D elements, this option is similar to the surface stress option (described in the Theory Reference for the
Mechanical APDL and Mechanical Applications), but is more general and applicable to nonlinear analysis. The element
used with this option does not provide any stiffness, mass, or load contributions. This option should only be used in
single-layer shells. Irrespective of other settings, SHELL181 provides stress and strain output at the center of the
layer.
SHELL181 uses a penalty method to relate the independent rotational degrees of freedom about the normal (to the
shell surface) with the in-plane components of displacements. The program chooses an appropriate penalty stiffness by
default. A drill stiffness factor can be specified via the SECCONTROLS command.
Element loads are described in Node and Element Loads. Pressures may be input as surface loads on the element faces
as shown by the circled numbers on Figure 181.1. Positive pressures act into the element. Because shell edge pressures
are input on a per-unit-length basis, per-unit-area quantities must be multiplied by the shell thickness.
Temperatures may be input as element body loads at the corners of the outside faces of the element and at the corners
of the interfaces between layers. The first corner temperature T1 defaults to TUNIF. If all other temperatures are
unspecified, they default to T1. If KEYOPT(1) = 0 and if exactly NL+1 temperatures are input, one temperature is used
for the four bottom corners of each layer, and the last temperature is used for the four top corner temperatures of the
top layer. If KEYOPT(1) = 1 and if exactly NL temperatures are input, one temperature is used for the four corners of
each layer. That is, T1 is used for T1, T2, T3, and T4; T2 (as input) is used for T5, T6, T7, and T8, etc. For any other
input pattern, unspecified temperatures default to TUNIF.
Using KEYOPT(3), SHELL181 supports uniform reduced integration and full integration with incompatible modes. By
default, this element uses the uniform reduced integration for performance reasons in nonlinear applications.
Using reduced integration with hourglass control creates some usage restrictions, although minimal. For example, to
capture the in-plane bending of a cantilever or a stiffener (see Figure 181.2), a number of elements through the
thickness direction is necessary. The performance gains achieved by using uniform reduced integration are significant
enough to offset the need to use more elements. In relatively well-refined meshes, hourglassing issues are largely
irrelevant.
When the reduced integration option is used, you can check the accuracy of the solution by comparing the total energy
(SENE label in ETABLE) and the artificial energy (AENE label in ETABLE) introduced by hourglass control. If the ratio
of artificial energy to total energy is less than 5%, the solution is generally acceptable. The total energy and artificial
energy can also be monitored by using OUTPR,VENG in the solution phase.
Bilinear elements, when fully integrated, are too stiff in in-plane bending.SHELL181 uses the method of incompatible
modes to enhance the accuracy in bending-dominated problems. This approach is also called "extra shapes" or "bubble"
modes approach. SHELL181 uses the formulation that ensures satisfaction of the patch test (J. C. Simo and F. Armero,
"Geometrically nonlinear enhanced strain mixed methods and the method of incompatible modes," IJNME, Vol. 33, pp.
1413-1449, 1992).
When including incompatible modes in the analysis, you must use full integration. KEYOPT(3) = 2 implies the inclusion
of incompatible modes and the use of full (2x2) quadrature.
SHELL181, with KEYOPT(3) = 2 specified, does not have any spurious energy mechanisms. This specific form of
SHELL181 is highly accurate, even with coarse meshes. We recommend that you use KEYOPT(3) = 2 if you encounter
any hourglass-related difficulties with the default options. KEYOPT(3) = 2 is also necessary if the mesh is coarse and inplane bending of the elements dominate the response. We recommend this option with all layered applications.
KEYOPT(3) = 2 imposes the fewest usage restrictions. You can always choose this option. However, you can improve
element performance by choosing the best option for your problem. Consider the problems illustrated in Figure 181.2
Element Reference
Page: 4
The cantilever beam and the beam cross-section to be modeled with shells are typical examples of in-plane bendingdominated problems. The use of KEYOPT(3) = 2 is the most effective choice in these circumstances. Reduced
integration would require refined meshes. For example, reduced integration for the cantilever beam problem requires
four elements through the thickness, whereas the full integration with incompatible modes only requires one element
through the thickness.
For the stiffened shell, the most effective choice is to use KEYOPT(3) = 0 for the shell and KEYOPT(3) = 2 for the
stiffener.
When KEYOPT(3) = 0 is specified, SHELL181 uses an hourglass control method for membrane and bending modes. By
default, SHELL181 calculates the hourglass parameters for both metal and hyperelastic applications. To specify the
hourglass stiffness scaling factors, use the SECCONTROLS command.
SHELL181 includes the linear effects of transverse shear deformation. An assumed shear strain formulation of BatheDvorkin is used to alleviate shear locking. The transverse shear stiffness of the element is a 2x2 matrix as shown
below:
In the above matrix, k = 5/6, G = shear modulus, and h = thickness of the shell.
SHELL181 can be associated with linear elastic, elastoplastic, creep, or hyperelastic material properties. Only isotropic,
anisotropic, and orthotropic linear elastic properties can be input for elasticity. The von Mises isotropic hardening
Contains proprietary and confidential information of ANSYS, Inc.
and its subsidiaries and affiliates
Element Reference
Page: 5
plasticity models can be invoked with BISO (bilinear isotropic hardening), MISO (multilinear isotropic hardening), and
NLISO (nonlinear isotropic hardening) options. The kinematic hardening plasticity models can be invoked with BKIN
(bilinear kinematic hardening), MKIN and KINH (multilinear kinematic hardening), and CHABOCHE (nonlinear kinematic
hardening). Invoking plasticity assumes that the elastic properties are isotropic (that is, if orthotropic elasticity is used
with plasticity, ANSYS assumes the isotropic elastic modulus = EX and Poisson's ratio = NUXY).
Hyperelastic material properties (2, 3, 5, or 9 parameter Mooney-Rivlin material model, Neo-Hookean model, Polynomial
form model, Arruda-Boyce model, and user-defined model) can be used with this element. Poisson's ratio is used to
specify the compressibility of the material. If less than 0, Poisson's ratio is set to 0; if greater than or equal to 0.5,
Poisson's ratio is set to 0.5 (fully incompressible).
Both isotropic and orthotropic thermal expansion coefficients can be input using MP,ALPX. When used with
hyperelasticity, isotropic expansion is assumed.
Use the BETAD command to specify the global value of damping. If MP,DAMP is defined for the material number of
the element (assigned with the MAT command), it is used for the element instead of the value from the BETAD
command. Similarly, use the TREF command to specify the global value of reference temperature. If MP,REFT is
defined for the material number of the element, it is used for the element instead of the value from the TREF
command. But if MP,REFT is defined for the material number of the layer, it is used instead of either the global or
element value.
With reduced integration and hourglass control (KEYOPT(3) = 0), low frequency spurious modes may appear if the
mass matrix employed is not consistent with the quadrature rule. SHELL181 uses a projection scheme that effectively
filters out the inertia contributions to the hourglass modes of the element. To be effective, a consistent mass matrix
must be used. We recommend setting LUMPM,OFF for a modal analysis using this element type. The lumped mass
option can, however, be used with the full integration options (KEYOPT(3) = 2).
KEYOPT(8) = 2 stores midsurface results in the results file for single or multi-layer shell elements. If you use SHELL,
MID, you will see these calculated values, rather than the average of the TOP and BOTTOM results. You should use this
option to access these correct midsurface results (membrane results) for those analyses where averaging TOP and
BOTTOM results is inappropriate; examples include midsurface stresses and strains with nonlinear material behavior,
and midsurface results after mode combinations that involve squaring operations such as in spectrum analyses.
KEYOPT(9) = 1 reads initial thickness data from a user subroutine.
You can apply an initial stress state to this element via the INISTATE command. For more information, see "Initial
State" in the Basic Analysis Guide.
The effects of pressure load stiffness are automatically included for this element. If an unsymmetric matrix is needed for
pressure load stiffness effects, use NROPT,UNSYM.
A summary of the element input is given in "SHELL181 Input Summary". A general description of element input is
given in Element Input.
Element Reference
Page: 6
Element Reference
Page: 7
) can be recovered from stress resultants with respect to the reference plane as
Element Reference
Page: 8
SHELL181 does not support extensive basic element printout. POST1 provides more comprehensive output processing
tools; therefore, ANSYS suggests using the OUTRES command to ensure that the required results are stored in the
database.
Definition
O R
EL
Y Y
NODES
Nodes - I, J, K, L
- Y
MAT
Material number
- Y
THICK
Average thickness
- Y
VOLU:
Volume
- Y
XC, YC, ZC
Element Reference
Page: 9
PRES
- Y
TEMP
T1, T2, T3, T4 at bottom of layer 1, T5, T6, T7, T8 between layers 1-2, similarly
for between next layers, ending with temperatures at top of layer NL(4*(NL+1)
maximum)
- Y
LOC
- 1
Stresses
3 1
S:1, 2, 3
Principal stresses
S:INT
Stress intensity
- 1
S:EQV
Equivalent stress
- 1
EPEL:X, Y, Z, XY
Elastic strains
3 1
EPEL:EQV
EPTH:X, Y, Z, XY
Thermal strains
3 1
EPTH:EQV
EPPL:X, Y, Z, XY
3 2
EPPL:EQV
EPCR:X, Y, Z, XY
3 2
EPCR:EQV
EPTO:X, Y, Z, XY
3 -
EPTO:EQV
NL:EPEQ
- 2
NL:CREQ
NL:SRAT
- 2
NL:PLWK
Plastic work
- 2
NL:HPRES
Hydrostatic pressure
- 2
SEND:ELASTIC,
PLASTIC, CREEP
- Y
- 8
Q13, Q23
- 8
11, 22, 12
Membrane strains
- Y
Curvatures
- 8
13, 23
- 8
LOCI:X, Y, Z
SVAR:1, 2, ... , N
State variables
ILSXZ
Element Reference
Page: 10
ILSYZ
ILSUM
ILANG
Angle of interlaminar shear stress vector (measured from the element x-axis
toward the element y-axis in degrees)
Membrane stresses
Bending stresses
Peak stresses
St: 13, 23
1. The following stress solution repeats for top, middle, and bottom surfaces.
2. Nonlinear solution output for top, middle, and bottom surfaces, if the element has a nonlinear material.
3. Stresses, total strains, plastic strains, elastic strains, creep strains, and thermal strains in the element coordinate
system are available for output (at all section points through thickness). If layers are in use, the results are in
the layer coordinate system.
4. Available only at centroid as a *GET item.
5. Available only if OUTRES,LOCI is used.
6. Available only if the USERMAT subroutine and TB,STATE are used.
7. The equivalent strains use an effective Poisson's ratio: for elastic and thermal this value is set by the user
(MP,PRXY); for plastic and creep this value is set at 0.5.
8. Not available if the membrane element option is used (KEYOPT(1) = 1).
Table 181.2: SHELL181 Item and Sequence Numbers lists output available through ETABLE using the Sequence
Number method. See Creating an Element Table in the Basic Analysis Guide and The Item and Sequence Number
Table in this manual for more information. The following notation is used in Table 181.2: SHELL181 Item and
Sequence Numbers:
Name
output quantity as defined in the Table 181.1: SHELL181 Element Output Definitions
Item
predetermined Item label for ETABLE
E
sequence number for single-valued or constant element data
I,J,K,L
sequence number for data at nodes I, J, K, L
N11
SMISC
N22
SMISC
N12
SMISC
M11
SMISC
M22
SMISC
M12
SMISC
Q13
SMISC
Q23
SMISC
Element Reference
Page: 11
11
SMISC
22
SMISC
10
12
SMISC
11
k11
SMISC
12
k22
SMISC
13
k12
SMISC
14
13
SMISC
15
23
SMISC
16
THICK
SMISC
17
P1
SMISC
18
19
20
21
P2
SMISC
22
23
24
25
P3
SMISC
27
26
P4
SMISC
P5
SMISC
P6
SMISC
Sm: 11
SMISC
34
Sm: 22
SMISC
35
Sm: 12
SMISC
36
Sb: 11
SMISC
37
Sb: 22
SMISC
38
Sb: 12
SMISC
39
SMISC
40
SMISC
41
SMISC
42
SMISC
43
SMISC
44
SMISC
45
St: 13
SMISC
46
St: 23
SMISC
47
29
28
-
32
31
30
-
33
Bottom of Layer i
Top of Layer NL
ILSXZ
SMISC
8 * (i - 1) + 51
8 * (NL - 1) + 52
ILSYZ
SMISC
8 * (i - 1) + 53
8 * (NL - 1) + 54
ILSUM
SMISC
8 * (i - 1) + 55
8 * (NL - 1) + 56
Element Reference
Page: 12
ILANG
SMISC
8 * (i - 1) + 57
8 * (NL - 1) + 58