SHELL229


4-Node Coupled-field Shell

Valid Products: Pro | Premium | Enterprise | PrepPost | Solver | AS add-on

SHELL229 Element Description

SHELL229 is a 3D layered shell element having in-plane thermal and electrical capabilities (membrane behavior) and it supports the thermal-electric physics combination. The element has four nodes with no limitation on the number of material layers (defined by the SECDATA command). It has two degrees of freedom at each node, temperature and voltage. SHELL229 generates temperatures that can be passed to structural shell elements in order to model thermal strain. See SHELL229 in the Mechanical APDL Theory Reference for more details about this element. For verification examples, see VM215.

If the model containing the thermal-electric shell element is to be analyzed structurally, use an equivalent structural shell element instead (SHELL181 or SHELL281).

Figure 229.1: SHELL229 Geometry

SHELL229 Geometry

xo = element x-axis if ESYS is not supplied.

x = element x-axis if ESYS is supplied.


SHELL229 Input Data

The geometry, node locations, and coordinates systems for this element are shown in Figure 229.1: SHELL229 Geometry. The element is defined by four nodes, one thickness per layer, a material angle for each layer, and the material properties.

The cross-sectional properties are input using the SECTYPE,,SHELL and SECDATA commands. These properties are the thickness, material number, and orientation of each layer. The number of integration points through the thickness is always 1, irrespective of the number entered in the SECDATA command. Real constants are not used for this element.

The default orientation for this element has the S1 (shell surface coordinate) axis aligned with the element's first parametric direction at its center.

The default first surface direction S1 can be reoriented in the element reference plane (as shown in Figure 229.1: SHELL229 Geometry) using the ESYS command. You can further rotate S1 by angle THETA (in degrees) for each layer via the SECDATA command to create layer-wise coordinate systems. See Coordinate Systems for details.

KEYOPT(1) determines the element degree-of-freedom (DOF) set and the corresponding force labels and reaction solution. KEYOPT(1) is set equal to the sum of the field keys shown in Table 229.1: SHELL229 Field Keys. Only KEYOPT(1) = 110 is supported and it activates a thermal-electric analysis with TEMP and VOLT as DOF labels and heat flow and electric current as reaction solutions.

Table 229.1: SHELL229 Field Keys

Field Field Key DOF Label Force Label Reaction Solution
Thermal10TEMPHEATHeat Flow
Electric Conduction100VOLTAMPSElectric Current

The coupled-field analysis KEYOPT(1) settings, DOF labels, force labels, reaction solutions, and analysis types are shown in the following table.

Table 229.2: SHELL229 Coupled-Field Analysis

Coupled-Field Analysis KEYOPT(1) DOF Label Force Label Reaction Solution Analysis Type
Thermal-Electric110TEMP,

VOLT

HEAT,

AMPS

Heat Flow,

Electric Current

Static,

Full Transient


As shown in the following table, material property requirements consist of those required for the individual fields (thermal and electric conduction) and those required for field coupling.

Table 229.3: SHELL229 Material Properties and Material Models

Coupled-Field Analysis KEYOPT(1) Material Properties and Material Models
Thermal-Electric110ThermalKXX, KYY, KZZ, DENS, C, ENTH, HF
ElectricRSVX, RSVY, RSVZ
Coupling–-

The density, thermal conductivity, specific heat, and electrical resistivity can all be defined either with the MP command or the TB command.

The electrical resistivity can be defined as a function of primary variables by using tabular input on the MP command. For more information, see Defining Linear Material Properties Using Tabular Input in the Material Reference. Alternatively, TB,ELEC can be used to define electrical resistivity with TBOPT = RSV or electrical conductivity with TBOPT = COND. For more information, see Anisotropic Electrical Conductivity.

Nodal loads are defined with the D and F commands.

Element loads are described in Element Loading. Loads may be input on the element faces indicated by the circled numbers in Figure 229.1: SHELL229 Geometry using the SFE command. Body loads may be input at the element's material layers using the BFE commands.

SHELL229 surface and body loads are given in the following table. Most surface and body loads can be defined as a function of primary variables by using tabular input. For more information, see Applying Loads Using Tabular Input in the Basic Analysis Guide and the individual surface or body load command description in the Command Reference.

Table 229.4: SHELL229 Surface and Body Loads

Coupled-Field Analysis KEYOPT(1) Load Type Load Command Label
Thermal-Electric110SurfaceConvection

Heat Flux

Radiation

CONV

HFLUX

RDSF

Body

Heat generation on material layers

HG(1), HG(2), . . . , HG(number of layers)

HGEN

The surface loads are input on a per-unit-area basis at all six faces including the shell edges:

  • Face 1 (I-J-K-L) (bottom, -z side)

  • Face 2 (I-J-K-L) (top, +z side)

  • Face 3 (J-I)

  • Face 4 (K-J)

  • Face 5 (L-K)

  • Face 6 (I-L)

This is in contrast to SHELL157 where surface loads on shell edges (faces 3 to 6 above) are input on a per-unit-length basis.

Radiation is not available on the shell edges. Convection and heat flux cannot be applied both on the same face. The surface loads can be defined by the SFE command only. You can also generate film coefficients and bulk temperatures using the surface effect element SURF152. SURF152 can also be used with FLUID116.

Element body loads may be input on a per layer basis. They can be defined with the BFE command only. One body load value is applied to the entire layer. If the first layer body load is input, and all others are unspecified, they default to the value specified for the first layer.

The body loads can be applied on material layers, where HG(1), HG(2), HG(3), . . . ., HG (number of material layers) can be defined with the BFE command only.

A summary of the element input is given in "SHELL229 Input Summary". A general description of element input is given in Element Input.

SHELL229 Input Summary

Nodes

I, J, K, L

Degrees of Freedom
Material Properties
Surface Loads
Body Loads
KEYOPT(1)
Element degrees of freedom. See Table 229.2: SHELL229 Coupled-Field Analysis.
KEYOPT(2)

Coupling method between the DOFs:

0 -- 

Strong (matrix) coupling. May produce an unsymmetric matrix. In a linear analysis, a coupled response is achieved after one iteration.

1 -- 

Weak (load vector) coupling. Produces a symmetric matrix and requires at least two iterations to achieve a coupled response. 

KEYOPT(7)

Evaluation of film coefficient:

0 -- 

Evaluate film coefficient (if any) at average film temperature, (TS+TB)/2.

1 -- 

Evaluate at element surface temperature, TS.

2 -- 

Evaluate at fluid bulk temperature, TB.

3 -- 

Evaluate at differential temperature, |TS-TB|.

KEYOPT(8)

Material layer data storage:

0 -- 

Store data for bottom of bottom layer and top of top layer (default).

1 -- 

Store top and bottom data for all layers (the volume of data may be considerable).

KEYOPT(13)

Film coefficient matrix:

0 -- 

The program determines whether to use a diagonal or consistent film coefficient matrix (default).

1 -- 

Use a diagonal film coefficient matrix.

2 -- 

Use a consistent film coefficient matrix.

KEYOPT(15)

Specific heat matrix:

0 -- 

The program determines whether to use a diagonal or consistent specific heat matrix (default).

1 -- 

Use a diagonal specific heat matrix.

2 -- 

Use a consistent specific heat matrix.

SHELL229 Output Data

The solution output associated with the element is in two forms:

Output temperatures may be read by structural shell elements SHELL181 and SHELL281 via the LDREAD,TEMP command.

Convection heat flux is positive out of the element; applied heat flux is positive into the element.

The element output directions are parallel to the layer coordinate system.

A general description of solution output is given in Solution Output. See the Basic Analysis Guide for ways to view results.

To see the distribution of element quantities through the thickness for this element, enter the POST1 postprocessor (/POST1), then issue /GRAPHICS,POWER and /ESHAPE,1 followed by PLESOL.

The Element Output Definitions table uses the following notation:

A colon (:) in the Name column indicates that the item can be accessed by the Component Name method (ETABLE, ESOL). The O column indicates the availability of the items in the file jobname.out. The R column indicates the availability of the items in the results file.

In either the O or R columns, “Y” indicates that the item is always available, a letter or number refers to a table footnote that describes when the item is conditionally available, and “-” indicates that the item is not available.

Table 229.5: SHELL229 Element Output Definitions

Name Definition O R
ELElement NumberYY
NODESNodes - I, J, K, LYY
MATMaterial numberYY
VOLU:VolumeYY
XC, YC, ZCLocation where results are reportedY 2
THERMAL-ELECTRIC ANALYSIS (KEYOPT(1) = 110)
TG:X, Y, Z, SUMThermal gradient components and vector magnitude-1
TF:X, Y, Z, SUMThermal flux components and vector magnitude-1
EF: X, Y, Z, SUM Electric field components and vector magnitude -1
JC: X, Y, Z, SUM Conduction current density components and vector magnitude-1
JS: X, Y, Z, SUM Current density components (in the global Cartesian coordinate system) and vector magnitude -1
JHEATJoule heat generation per unit volume -1
AREAFace area11
HFAVGAverage film coefficient of the face-1
TAVGAverage face temperature11
TBAVGAverage bulk temperature1-
HEAT RATEHeat flow rate across face by convection11
HFLXAVGHeat flow rate per unit area across face caused by input heat flux-1
THICKNESSShell thickness-1

  1. Solution values are output only if calculated (based on input values).

  2. Available only at the centroid as a *GET item.

Table 229.5: SHELL229 Element Output Definitions lists output available through the ETABLE command using the Sequence Number method. See Element Table for Variables Identified By Sequence Number in the Basic Analysis Guide and The Item and Sequence Number Table in this reference for more information. The following notation is used in Table 229.6: SHELL229 Item and Sequence Numbers:

Name

output quantity as defined in Table 229.5: SHELL229 Element Output Definitions

Item

predetermined Item label for ETABLE command

FCn

sequence number for solution items for element Face n

Table 229.6: SHELL229 Item and Sequence Numbers

Output Quantity NameETABLE and ESOL Command Input
ItemFC1FC2FC3FC4FC5FC6
AREANMISC1713192531
HFAVGNMISC2814202632
TAVGNMISC3915212733
TBAVGNMISC41016222834
HEAT RATENMISC51117232935
HFLXAVGNMISC61218243036
THICKNESSNMISC37

SHELL229 Assumptions and Restrictions

  • Zero thickness layers are not allowed.

  • The tabular thickness specified by SECFUNCTION is not supported.

  • The Mass transport is not supported.

  • Nonlinear material properties are evaluated at each integration point.

  • This element may not be compatible with other elements with the VOLT degree of freedom. To be compatible, the elements must have the same reaction solution for the VOLT DOF. For more information, see Element Compatibility in the Low-Frequency Electromagnetic Analysis Guide.

  • Source current density, JS, and Joule heat, JHEAT, are available as both element and layer average values if KEYOPT(8) = 1, and only as element average values if KEYOPT(8) = 0. The layer average values can be accessed in /POST1 by issuing the LAYER command, and in /POST26 by issuing the LAYERP26 command. If either LAYER,0 or LAYERP26,0 is issued, then the element average value is printed.

SHELL229 Product Restrictions

There are no product-specific restrictions for this element.