The Layered Shell Interface
The Layered Shell (lshell) interface (), found under the Structural Mechanics branch () when adding a physics interface, is used to model layered structural shells on 3D boundaries. Shells are thin flat or curved structures, having significant bending stiffness.
The Layered Shell interface is applicable for thick and moderately thin shells. The formulation resembles that of a stack of fully 3D solid mechanics models, so all stress components, including interlaminar shear stresses can be resolved. For very thin shells this formulation tends to become numerically ill-behaved, and it is then better to use the Shell interface with the Linear Elastic Material, Layered model.
The Linear Elastic Material is the default material model. With the Nonlinear Structural Materials Module, the physics interface is extended with more material models, like hyperelasticity, plasticity, creep, and viscoplasticity.
When this interface is added, these default nodes are also added to the Model BuilderLinear Elastic Material, Free (a boundary condition where edges are free, with no loads or constraints), and Initial Values. Then, from the Physics toolbar, add other nodes that implement, for example, boundary conditions. You can also right-click Layered Shell to select physics features from the context menu.
The Layered Shell interface is only available with the Composite Materials Module (see https://www.comsol.com/products/specifications/).
Settings
The Label is the default physics interface name.
The Name is used primarily as a scope prefix for variables defined by the physics interface. Refer to such physics interface variables in expressions using the pattern <name>.<variable_name>. In order to distinguish between variables belonging to different physics interfaces, the name string must be unique. Only letters, numbers, and underscores (_) are permitted in the Name field. The first character must be a letter.
The default Name (for the first physics interface in the model) is lshell.
Shell Properties
Here you select on which layers in a layered material that the physics interface should be active. By default, the Use all layers check box is selected. This means that all layers in all layered materials on the selected boundaries are used.
If you deselect the Use all layers check box, you can select individual layers within a single layered material. This is a seldom used option, since it means that the physics interface is restricted to the boundaries on which a specific layered material is defined.
Structural Transient Behavior
From the Structural transient behavior list, select Include inertial terms (the default) or Quasistatic. Use Quasistatic to treat the dynamic behavior as quasistatic (with no mass effects; that is, no second-order time derivatives). Selecting this option gives a more efficient solution for problems where the variation in time is slow when compared to the natural frequencies of the system. The default solver for the time stepping is changed from Generalized alpha to BDF when Quasistatic is selected.
For problems with creep, and sometimes viscoelasticity, the problem can be considered as quasistatic. This is also the case when the time dependence exists only in some other physics, like a transient heat transfer problem causing thermal strains.
Fold-Line Settings
A fold-line is an edge where the boundaries that meet do not have a continuous orientation of the normal vector. Along fold lines, the degrees of freedom are in general not continuous, but rather connected by a set of constraints.
Select a Constraint typeFull or Simplified. The full fold-line constraint formulation is accurate, consistent with the degrees of freedom, and does not produce any stress singularities near the fold-lines. However, it is computationally expensive, and leads to a nonlinear problem, even in an otherwise linear model. On the other hand, the simplified fold-line constraint formulation assumes that thin shell kinematics can be used. It is computationally more efficient, does not force the problem to be nonlinear, and works very well for thin or moderately thick shells.
The fold-line limit angle α is the smallest angle between the normals of two boundaries that makes their intersection to be treated as a fold line. The normal to the layered shell is discontinuous along a fold-line. Enter a value or expression in the α field. The default value is 0.05 radians (approximately 3°). The value must be larger than 0, and less than π/2, but angles larger than a few degrees are not usually meaningful.
See Fold-Line Connection for the theory of fold line constraints in layered shell.
Default Through-Thickness Result Location
Enter a number between -1 and 1 for the Local z-coordinate [-1,1] for thickness-dependent results Z. The value can be changed from 1 (downside) to +1 (upside). The default is +1. A value of 0 means the midsurface of the layered shell. This is the default position for stress and strain evaluation during the results analysis.
Advanced Settings
To display this section, click the Show More Options button () and select Advanced Physics Options in the Show More Options dialog box. Normally these settings do not need to be changed.
You can choose how extra ODE variables added by some features are grouped in the Dependent Variables node of a generated solver sequence.
Select the Rigid materials check box to group variables added by Rigid Material nodes.
Select the Rigid connectors check box to group variables added by Rigid Connector and Rigid Connector, Interface nodes.
Select the Attachments check box to group variables added by Attachment nodes.
The selection made in the Advanced Settings section can be overridden by the settings in the Advanced section of the Rigid Material, Rigid Connector, Rigid Connector, Interface, or Attachment features.
Discretization
In the Layered Shell interface you can choose not only the order of the discretization, but also the type of shape functions: Lagrange or serendipity. For highly distorted elements, Lagrange shape functions provide better accuracy than serendipity shape functions of the same order. The serendipity shape functions will however give significant reductions of the model size for a given mesh containing hexahedral, prism, or quadrilateral elements. The default is to use Quadratic Lagrange shape functions for the Displacement field.
The order of the discretization is used not only in the reference surface but also in the thickness direction. It is possible to choose different order of discretization in the reference surface and in the thickness direction.
See Discretization for the details on discretization order for a layered shell.
Dependent Variables
The physics interface uses the global spatial components of the Displacement field u as dependent variables. The default names for the components are (u, v, w).
Since the Layered Shell Interface uses a discretization also in the thickness direction, the degrees of freedom are allocated not only in the plane of the shell, but also in that direction.
You can change both the field name and the individual component names. If a new field name coincides with the name of another displacement field, the two fields (and the interfaces which define them) share degrees of freedom and dependent variable component names.
A new field name must not coincide with the name of a field of another type (that is, it must contain a displacement field), or with a component name belonging to some other field. Component names must be unique within a model except when two interfaces share a common field name.
Micromechanics and Macromechanics of a Composite Cylinder: Application Library path Composite_Materials_Module/Tutorials/composite_cylinder_micromechanics_and_stress_analysis
Forced Vibration Analysis of a Composite Laminate: Application Library path Composite_Materials_Module/Dynamics_and_Vibration/forced_vibration_of_a_composite_laminate