The Structural Mechanics Module User's Guide
The Structural Mechanics Module User's Guide
Users Guide
VERSION 4.3
Structural Mechanics Module Users Guide
19982012 COMSOL
Protected by U.S. Patents 7,519,518; 7,596,474; and 7,623,991. Patents pending.
This Documentation and the Programs described herein are furnished under the COMSOL Software License
Agreement (www.comsol.com/sla) and may be used or copied only under the terms of the license agree-
ment.
COMSOL, COMSOL Desktop, COMSOL Multiphysics, and LiveLink are registered trademarks or trade-
marks of COMSOL AB. Other product or brand names are trademarks or registered trademarks of their
respective holders.
Version: May 2012 COMSOL 4.3
Contact Information
Visit www.comsol.com/contact for a searchable list of all COMSOL offices and local
representatives. From this web page, search the contacts and find a local sales
representative, go to other COMSOL websites, request information and pricing, submit
technical support queries, subscribe to the monthly eNews email newsletter, and much
more.
Chapter 1: Introduction
Applying Loads 39
Units, Orientation, and Visualization . . . . . . . . . . . . . . . 39
Load Cases . . . . . . . . . . . . . . . . . . . . . . . . . 40
Singular Loads . . . . . . . . . . . . . . . . . . . . . . . . 41
Moments in the Solid Mechanics Interface . . . . . . . . . . . . . 42
Pressure . . . . . . . . . . . . . . . . . . . . . . . . . . 42
Acceleration Loads . . . . . . . . . . . . . . . . . . . . . . 43
Temperature LoadsThermal Expansion. . . . . . . . . . . . . . 43
Total Loads . . . . . . . . . . . . . . . . . . . . . . . . . 43
Defining Constraints 44
Orientation . . . . . . . . . . . . . . . . . . . . . . . . . 44
Symmetry Constraints . . . . . . . . . . . . . . . . . . . . . 44
Kinematic Constraints . . . . . . . . . . . . . . . . . . . . . 46
Rotational Joints . . . . . . . . . . . . . . . . . . . . . . . 46
CONTENTS |3
Calculating Reaction Forces 47
Using Predefined Variables to Evaluate Reaction Forces . . . . . . . . 47
Using Weak Constraints to Evaluate Reaction Forces . . . . . . . . . 48
Using Surface Traction to Evaluate Reaction Forces . . . . . . . . . . 49
Eigenfrequency Analysis 69
4 | CONTENTS
Linear Viscoelastic Materials . . . . . . . . . . . . . . . . . . 77
Rayleigh Damping. . . . . . . . . . . . . . . . . . . . . . . 77
Equivalent Viscous Damping. . . . . . . . . . . . . . . . . . . 78
Loss Factor Damping . . . . . . . . . . . . . . . . . . . . . 79
Explicit Damping . . . . . . . . . . . . . . . . . . . . . . . 80
Piezoelectric Losses 81
About Piezoelectric Materials . . . . . . . . . . . . . . . . . . 81
Piezoelectric Material Orientation . . . . . . . . . . . . . . . . 82
Piezoelectric Losses. . . . . . . . . . . . . . . . . . . . . . 88
No Damping . . . . . . . . . . . . . . . . . . . . . . . . 91
References for Piezoelectric Damping . . . . . . . . . . . . . . . 91
CONTENTS |5
The Solid Mechanics Interface 111
Domain, Boundary, Edge, Point, and Pair Features for the Solid Mechanics
Interface . . . . . . . . . . . . . . . . . . . . . . . . . 113
Linear Elastic Material . . . . . . . . . . . . . . . . . . . . 115
Change Thickness . . . . . . . . . . . . . . . . . . . . . 119
Damping . . . . . . . . . . . . . . . . . . . . . . . . . 120
Initial Values. . . . . . . . . . . . . . . . . . . . . . . . 123
About the Body, Boundary, Edge, and Point Loads . . . . . . . . . 123
Body Load . . . . . . . . . . . . . . . . . . . . . . . . 123
Boundary Load . . . . . . . . . . . . . . . . . . . . . . 124
Edge Load . . . . . . . . . . . . . . . . . . . . . . . . 126
Point Load . . . . . . . . . . . . . . . . . . . . . . . . 126
Free. . . . . . . . . . . . . . . . . . . . . . . . . . . 127
Fixed Constraint . . . . . . . . . . . . . . . . . . . . . . 127
Prescribed Displacement . . . . . . . . . . . . . . . . . . . 128
Symmetry . . . . . . . . . . . . . . . . . . . . . . . . 130
Antisymmetry . . . . . . . . . . . . . . . . . . . . . . . 131
Roller . . . . . . . . . . . . . . . . . . . . . . . . . . 131
Periodic Condition . . . . . . . . . . . . . . . . . . . . . 132
Perfectly Matched Layers . . . . . . . . . . . . . . . . . . . 133
Linear Viscoelastic Material . . . . . . . . . . . . . . . . . . 133
Thermal Effects . . . . . . . . . . . . . . . . . . . . . . 135
Rigid Connector . . . . . . . . . . . . . . . . . . . . . . 137
Harmonic Perturbation . . . . . . . . . . . . . . . . . . . 139
Applied Force . . . . . . . . . . . . . . . . . . . . . . . 139
Applied Moment . . . . . . . . . . . . . . . . . . . . . . 139
Mass and Moment of Inertia. . . . . . . . . . . . . . . . . . 140
Rigid Domain . . . . . . . . . . . . . . . . . . . . . . . 140
Pairs for the Solid Mechanics Interface. . . . . . . . . . . . . . 141
Contact . . . . . . . . . . . . . . . . . . . . . . . . . 141
Friction . . . . . . . . . . . . . . . . . . . . . . . . . 143
Thermal Expansion . . . . . . . . . . . . . . . . . . . . . 144
Initial Stress and Strain. . . . . . . . . . . . . . . . . . . . 146
Perfectly Matched Layers . . . . . . . . . . . . . . . . . . 147
Phase . . . . . . . . . . . . . . . . . . . . . . . . . . 147
Prescribed Velocity . . . . . . . . . . . . . . . . . . . . . 148
Prescribed Acceleration . . . . . . . . . . . . . . . . . . . 149
Spring Foundation . . . . . . . . . . . . . . . . . . . . . 150
6 | CONTENTS
Pre-Deformation . . . . . . . . . . . . . . . . . . . . . . 154
Thin Elastic Layer. . . . . . . . . . . . . . . . . . . . . . 155
Added Mass . . . . . . . . . . . . . . . . . . . . . . . . 157
Low-Reflecting Boundary . . . . . . . . . . . . . . . . . . . 158
CONTENTS |7
Thermal Expansion . . . . . . . . . . . . . . . . . . . . . 221
Initial Stress and Strain. . . . . . . . . . . . . . . . . . . . 222
Damping . . . . . . . . . . . . . . . . . . . . . . . . . 223
Change Thickness . . . . . . . . . . . . . . . . . . . . . 224
Initial Values. . . . . . . . . . . . . . . . . . . . . . . . 225
About the Body Load, Face Load, Edge Load, and Point Load Features . 226
Body Load . . . . . . . . . . . . . . . . . . . . . . . . 226
Face Load . . . . . . . . . . . . . . . . . . . . . . . . 227
Edge Load . . . . . . . . . . . . . . . . . . . . . . . . 228
Point Load . . . . . . . . . . . . . . . . . . . . . . . . 229
Phase . . . . . . . . . . . . . . . . . . . . . . . . . . 229
Pinned . . . . . . . . . . . . . . . . . . . . . . . . . . 230
No Rotation . . . . . . . . . . . . . . . . . . . . . . . 231
Prescribed Displacement/Rotation . . . . . . . . . . . . . . . 231
Prescribed Velocity . . . . . . . . . . . . . . . . . . . . . 233
Prescribed Acceleration . . . . . . . . . . . . . . . . . . . 235
Symmetry . . . . . . . . . . . . . . . . . . . . . . . . 236
Antisymmetry . . . . . . . . . . . . . . . . . . . . . . . 237
Rigid Connector . . . . . . . . . . . . . . . . . . . . . . 238
Results Evaluation . . . . . . . . . . . . . . . . . . . . . 240
Chapter 5: Beams
8 | CONTENTS
Section Orientation . . . . . . . . . . . . . . . . . . . . . 266
About the Edge Load and Point Load Features. . . . . . . . . . . 268
Edge Load . . . . . . . . . . . . . . . . . . . . . . . . 268
Point Load . . . . . . . . . . . . . . . . . . . . . . . . 269
Phase . . . . . . . . . . . . . . . . . . . . . . . . . . 270
Prescribed Displacement/Rotation . . . . . . . . . . . . . . . 271
Prescribed Velocity . . . . . . . . . . . . . . . . . . . . . 273
Prescribed Acceleration . . . . . . . . . . . . . . . . . . . 274
Pinned . . . . . . . . . . . . . . . . . . . . . . . . . . 276
No Rotation . . . . . . . . . . . . . . . . . . . . . . . 277
Symmetry . . . . . . . . . . . . . . . . . . . . . . . . 277
Antisymmetry . . . . . . . . . . . . . . . . . . . . . . . 279
Point Mass . . . . . . . . . . . . . . . . . . . . . . . . 280
Point Mass Damping . . . . . . . . . . . . . . . . . . . . 281
C h a p t e r 6 : Tr u s s e s
CONTENTS |9
Cross Section Data . . . . . . . . . . . . . . . . . . . . . 309
Initial Values. . . . . . . . . . . . . . . . . . . . . . . . 309
About the Edge Load and Point Load . . . . . . . . . . . . . . 310
Edge Load . . . . . . . . . . . . . . . . . . . . . . . . 310
Phase . . . . . . . . . . . . . . . . . . . . . . . . . . 311
Straight Edge Constraint . . . . . . . . . . . . . . . . . . . 311
Pinned . . . . . . . . . . . . . . . . . . . . . . . . . . 311
Prescribed Displacement . . . . . . . . . . . . . . . . . . . 312
Prescribed Velocity . . . . . . . . . . . . . . . . . . . . . 313
Prescribed Acceleration . . . . . . . . . . . . . . . . . . . 314
Symmetry . . . . . . . . . . . . . . . . . . . . . . . . 315
Antisymmetry . . . . . . . . . . . . . . . . . . . . . . . 316
Point Mass . . . . . . . . . . . . . . . . . . . . . . . . 316
Point Mass Damping . . . . . . . . . . . . . . . . . . . . 317
Chapter 7: Membranes
10 | C O N T E N T S
Chapter 8: Multiphysics Interfaces
CONTENTS | 11
Theory for the Piezoelectric Devices Interface 382
The Piezoelectric Effect . . . . . . . . . . . . . . . . . . . 382
Piezoelectric Constitutive Relations . . . . . . . . . . . . . . . 383
Piezoelectric Material . . . . . . . . . . . . . . . . . . . . 385
Piezoelectric Dissipation . . . . . . . . . . . . . . . . . . . 385
Initial Stress, Strain, and Electric Displacement. . . . . . . . . . . 386
Geometric Nonlinearity for the Piezoelectric Devices Interface . . . . 386
Damping and Losses Theory . . . . . . . . . . . . . . . . . 389
References for the Piezoelectric Devices Interface . . . . . . . . . 392
Chapter 9: Materials
12 | C O N T E N T S
1
Introduction
This chapter introduces you to the capabilities of this module and includes a
summary of the physics interfaces as well as information about where you can find
additional documentation and model examples. The last section is a brief overview
with links to each chapter in this guide.
13
About the Structural Mechanics
Module
In this section:
The physics interfaces in this module are fully multiphysics enabled, making it possible
to couple them to any other physics interfaces in COMSOL Multiphysics or the other
modules. Available physics interfaces include:
Solid mechanics for 2D plane stress and plane strain, axial symmetry, and 3D solids
Piezoelectric modeling
Beams in 2D and 3D, Euler theory
Truss and cable elements
Shells and plates
14 | CHAPTER 1: INTRODUCTION
The modules study capabilities include static, eigenfrequency, time dependent
(transient), frequency response, and parametric studies, as well as contact and friction.
There are also predefined interfaces for linear elastic and viscoelastic materials.
Materials can be isotropic, orthotropic, or fully anisotropic, and you can use local
coordinate systems to specify material properties. Large deformations can also be
included in a study.
Coupling structural analysis with thermal analysis is one example of multiphysics easily
implemented with the module, which provides predefined multiphysics interfaces for
thermal stress and other types of multiphysics. Piezoelectric materials, coupling the
electric field and strain in both directions are fully supported inside the module
through special multiphysics interfaces solving for both the electric potential and
displacements. Piezoelectric materials can also be analyzed with the constitutive
relations on either stress-charge or strain-charge form. Structural mechanics couplings
are common in simulations done with COMSOL Multiphysics and occur in interaction
with, for example, fluid flow (fluid-structure interaction, FSI), chemical reactions,
acoustics, electric fields, magnetic fields, and optical wave propagation.
STATIC ANALYSIS
In a static analysis the load and constraints are fixed in time.
EIGENFREQUENCY ANALYSIS
An eigenfrequency analysis finds the undamped eigenfrequencies and mode shapes of
a model. Sometimes referred to as the free vibration of a structure.
PARAMETRIC ANALYSIS
A parametric analysis finds the solution dependence due to the variation of a specific
parameter, which could be, for instance, a material property or the position of a load.
THERMAL STRESS
In a transient thermal stress study, the program neglects mass effects, assuming that the
time scale in the structural mechanics problem is much smaller than the time scale in
the thermal problem.
LARGE DEFORMATIONS
You can model large deformations via a number of available Hyperelastic Materials.
You can also enable the geometric nonlinearity for the Linear Elastic Material under all
solid mechanics interfaces except the Beam interface. The interface then replaces the
small strain with the Green-Lagrange strain and the stress with the second
Piola-Kirchhoff stress. Such material is suitable to study deformations accompanied by
possible large rotations but small to moderate strains in the material, and it is
sometimes referred to as Saint Venant-Kirchhoff hyperelastic material. To solve the
problem, the program uses a total Lagrangian formulation..
16 | CHAPTER 1: INTRODUCTION
Additional functionality and material models for geomechanics and soil
mechanicssoil plasticity, concrete, and rock material modelsis
Tip available with the Geomechanics Module.
ELASTOPLASTIC MATERIALS
An elastoplastic analysis involves a nonlinear material with or without hardening. Three
different hardening models are available:
Isotropic
Kinematic
Perfectly plastic
HYPERELASTIC MATERIALS
In hyperelastic materials the stresses are computed from a strain energy density
function. They are often used to model rubber, but also used in acoustic elasticity. Four
different models are available:
Neo-Hookean
Mooney-Rivlin
Murnaghan
User-defined by providing the strain energy function
CONTACT MODELING
You can model contact between parts of a structure. The Solid Mechanics interface
supports contact with or without friction. The contact algorithm is implemented based
on the augmented Lagrangian method.
At any time, a new model can be created or physics added. Right-click the Root (top)
node and select Add Model to open the Model Wizard, or right-click a Model node and
select Add Physics.
18 | CHAPTER 1: INTRODUCTION
domains, except for shells, beams, and trusses, where they are defined on the boundary
or edge level.
Fluid Flow
20 | CHAPTER 1: INTRODUCTION
PHYSICS ICON TAG SPACE PRESET STUDIES
DIMENSION
STATIONARY STUDY
A static analysis solves for stationary displacements, rotations, and temperature
(depending on the type of physics interface). All loads and constraints are constant in
time. The equations include neither mass nor mass moment of inertia.
EIGENFREQUENCY STUDY
An eigenfrequency study solves for the eigenfrequencies and the shape of the
eigenmodes. When performing an eigenfrequency analysis, specify whether to look at
the mathematically more fundamental eigenvalue, , or the eigenfrequency, f, which is
more commonly used in a structural mechanics context.
f = ---------
2i
TIME-DEPENDENT STUDY
A time-dependent (transient) study solves a time-dependent (unsteady) problem
where loads and constraints can vary in time. Time dependent studies can be
performed using either a direct or a modal method.
C = dM M + dK K
For a frequency domain study, the module divides harmonic loads into two parts:
The amplitude, F
The phase (FPh)
Together they define a harmonic load whose amplitude and phase shift can depend on
the excitation angular frequency or excitation frequency f.
F freq = F cos t + F Ph ----------
180
= 2f
COMSOL Multiphysics allows the visualization of the amplitudes and phases as well
as the solution at a specific angle (time). The Solution at angle parameter makes this
i
task easy. When plotting the solution, the program multiplies it by e , where is the
angle in radians that corresponds to the angle (specified in degrees) in the Solution at
angle field. COMSOL Multiphysics plots the real part of the evaluated expression:
22 | CHAPTER 1: INTRODUCTION
The angle is available as the variable phase (in radians) and is allowed in plotting
expressions. Both freq and omega are available variables.
Another way to calculate the critical load is to include large deformation effects and
increase the load until the load has reached its critical value. Linear buckling is available
in the Solid Mechanics, Shell, Plate, and Truss interfaces.
These studies involve two study steps for the solver (a Static study step
plus an Eigenfrequency or Frequency Domain study step). You need to
add a new study to the model to get access to such combined study types,
and they cannot be added directly as new study steps to the existing study
Note
(solver sequence). You also have them available when starting a new
model.
Eigenfrequency Analysis
See Also
V I S C O E L A S T I C TR A N S I E N T I N I T I A L I Z A T I O N
A viscoelastic transient initialization precomputes initial states for transient and
quasi-static transient analyses when the Linear Viscoelastic Material is used. It is a
regime of instantaneous deformation and/or loading.
THERMAL COUPLINGS
Solids expand with temperature, which causes thermal strains to develop in the
material. These thermal strains combine with the elastic strains from structural loads
to form the total strain:
= el + th
th = T T ref
Thermal expansion affects displacements, stresses, and strains. This effect is added
automatically in The Thermal Stress Interface and The Joule Heating and Thermal
Expansion Interface. Also add thermal expansion to the other interfaces. Only the
coefficient of thermal expansion needs to be specified and the two temperature fields,
24 | CHAPTER 1: INTRODUCTION
T and Tref. The temperature field is a model input that typically is computed by a Heat
Transfer interface. Temperature coupling can be used in any type of study.
For each physics interface, the table indicates the availability of various analysis
capabilities. Finally the table lists the geometry levels (where data such as material
properties, loads, and constraints are specified). Edges exist only in 3D geometries.
NONLINEAR MATERIALS
DEPENDENT VARIABLES
FREQUENCY RESPONSE
LARGE DEFORMATION
EIGENFREQUENCY
LINEAR BUCKLING
DEFAULT NAME
STATIONARY
BOUNDARIES
PARAMETRIC
TRANSIENT
DOMAINS
POINTS
EDGES
STRUCTURAL MECHANICS
26 | CHAPTER 1: INTRODUCTION
Multiphysics Users Guide or COMSOL Multiphysics Reference Guide where
additional information is available.
To display additional features for the physics interfaces and feature nodes, click the
Show button ( ) on the Model Builder and then select the applicable option.
After clicking the Show button ( ), some sections display on the settings window
when a node is clicked and other features are available from the context menu when a
node is right-clicked. For each, the additional sections that can be displayed include
Equation, Advanced Settings, Discretization, Consistent Stabilization, and Inconsistent
Stabilization.
You can also click the Expand Sections button ( ) in the Model Builder to always show
some sections or click the Show button ( ) and select Reset to Default to reset to
display only the Equation and Override and Contribution sections.
For most physics nodes, both the Equation and Override and Contribution sections are
always available. Click the Show button ( ) and then select Equation View to display
the Equation View node under all physics nodes in the Model Builder.
Availability of each feature, and whether it is described for a particular physics node, is
based on the individual physics selected. For example, the Discretization, Advanced
Show More Options and Showing and Expanding Advanced Users Guide
Expand Sections Physics Sections
The Model Builder Window
Discretization Show Discretization Users Guide
Element Types and Discretization
Finite Elements Reference Guide
Discretization of the Equations
Discretization - Splitting Compile Equations Reference Guide
of complex variables
Pair Selection Identity and Contact Pairs Users Guide
Specifying Boundary Conditions for
Identity Pairs
Consistent and Show Stabilization Users Guide
Inconsistent Stabilization
Stabilization Techniques Reference Guide
Numerical Stabilization
Geometry Working with Geometry Users Guide
Constraint Settings Using Weak Constraints Users Guide
28 | CHAPTER 1: INTRODUCTION
documentation, Dynamic Help, and the Model Library are all accessed through the
COMSOL Desktop.
If you are reading the documentation as a PDF file on your computer, the
blue links do not work to open a model or content referenced in a
different users guide. However, if you are using the online help in
Important COMSOL Multiphysics, these links work to other modules, model
examples, and documentation sets.
THE DOCUMENTATION
The COMSOL Multiphysics Users Guide and COMSOL Multiphysics Reference
Guide describe all interfaces and functionality included with the basic COMSOL
Multiphysics license. These guides also have instructions about how to use COMSOL
Multiphysics and how to access the documentation electronically through the
COMSOL Multiphysics help desk.
and then either enter a search term or look under a specific module in the
documentation tree.
SI units are used to describe the relevant properties, parameters, and dimensions in
most examples, but other unit systems are available.
To open the Model Library, select View>Model Library ( ) from the main menu, and
then search by model name or browse under a module folder name. Click to highlight
any model of interest, and select Open Model and PDF to open both the model and the
documentation explaining how to build the model. Alternatively, click the Dynamic
The model libraries are updated on a regular basis by COMSOL in order to add new
models and to improve existing models. Choose View>Model Library Update ( ) to
update your model library to include the latest versions of the model examples.
If you have any feedback or suggestions for additional models for the library (including
those developed by you), feel free to contact us at [email protected].
To receive technical support from COMSOL for the COMSOL products, please
contact your local COMSOL representative or send your questions to
[email protected]. An automatic notification and case number is sent to you by
email.
Typographical Conventions
All COMSOL users guides use a set of consistent typographical conventions that make
it easier to follow the discussion, understand what you can expect to see on the
graphical user interface (GUI), and know which data must be entered into various
data-entry fields.
30 | CHAPTER 1: INTRODUCTION
In particular, these conventions are used throughout the documentation:
CONVENTION EXAMPLE
Caution A Caution icon is used to indicate that the user should proceed
carefully and consider the next steps. It might mean that an
action is required, or if the instructions are not followed, that
there will be problems with the model solution.
Important An Important icon is used to indicate that the information
provided is key to the model building, design, or solution. The
information is of higher importance than a note or tip, and the
user should endeavor to follow the instructions.
Note A Note icon is used to indicate that the information may be of
use to the user. It is recommended that the user read the text.
Tip A Tip icon is used to provide information, reminders, short
cuts, suggestions of how to improve model design, and other
information that may or may not be useful to the user.
See Also The See Also icon indicates that other useful information is
located in the named section. If you are working on line, click
the hyperlink to go to the information directly. When the link is
outside of the current PDF document, the text indicates this,
for example See The Laminar Flow Interface in the
COMSOL Multiphysics Users Guide. Note that if you are in
COMSOL Multiphysics online help, the link will work.
Model The Model icon is used in the documentation as well as in
COMSOL Multiphysics from the View>Model Library menu. If
you are working online, click the link to go to the PDF version
of the step-by-step instructions. In some cases, a model is only
available if you have a license for a specific module. These
examples occur in the COMSOL Multiphysics Users Guide.
The Model Library path describes how to find the actual model
in COMSOL Multiphysics, for example
If you have the RF Module, see Radar Cross Section: Model
Library path RF_Module/Tutorial_Models/radar_cross_section
Space Dimension Another set of icons are also used in the Model Builderthe
model space dimension is indicated by 0D , 1D , 1D
axial symmetry , 2D , 2D axial symmetry , and 3D
icons. These icons are also used in the documentation to
clearly list the differences to an interface, feature node, or
theory section, which are based on space dimension.
32 | CHAPTER 1: INTRODUCTION
ABOUT THE STRUCTURAL MECHANICS MODULE | 33
Overview of the Users Guide
The Structural Mechanics Module Users Guide gets you started with modeling using
COMSOL Multiphysics. The information in this guide is specific to this module.
Instructions how to use COMSOL in general are included with the COMSOL
Multiphysics Users Guide.
TA B L E O F C O N T E N T S , G L O S S A R Y, A N D I N D E X
To help you navigate through this guide, see the Contents, Glossary, and Index.
34 | CHAPTER 1: INTRODUCTION
THE BEAM INTERFACE
The Beams chapter describes The Beam Interface, which models Euler
(Euler-Bernoulli) beams for modeling slender 3D and 2D structures. Typical examples
are frameworks and latticeworks. The underlying theory for the interface is also
included at the end of the chapter.
T H E TR U S S I N T E R F A C E
The Trusses chapter describes The Truss Interface, which models slender 3D and 2D
structures with components capable to withstand axial forces only. Typical example is
latticeworks. The underlying theory for the interface is also included at the end of the
chapter.
The Thermal Stress Interface combines a Solid Mechanics interface with a Heat
Transfer interface. The coupling appears on the domain level, where the
temperature from the Heat Transfer interface acts as a thermal load for the Solid
Mechanics interface, causing thermal expansion.
The Joule Heating and Thermal Expansion Interface combines solid mechanics
using a thermal linear elastic material with an electromagnetic Joule heating model.
This is a multiphysics combination of solid mechanics, electric currents, and heat
transfer for modeling of, for example, thermoelectromechanical (TEM)
applications.
The Piezoelectric Devices Interface include a piezoelectric material but also full
functionality for Solid Mechanics and Electrostatics. Piezoelectric materials in 3D,
2D plane strain and plane stress, and axial symmetry, optionally combined with
other solids and air, for example.
The Fluid-Structure Interaction Interface, which is found under the Fluid Flow branch
of the Model Wizard, is also described in this chapter. The interface combines fluid
flow with solid mechanics to capture the interaction between the fluid and the solid
structure. A Solid Mechanics interface and a single-phase flow interface model the
O V E R V I E W O F T H E U S E R S G U I D E | 35
solid and the fluid, respectively. The fluid-structure interactions appear on the
boundaries between the fluid and the solid.
MATERIALS
The Materials chapter describes the material databases included with the Structural
Mechanics ModuleLiquids and Gases, with temperature-dependent fluid dynamic
and thermal properties, MEMS, an extended solid materials library with metals,
semiconductors, insulators, and polymers common in MEMS devices, and a
Piezoelectric database with over 20 common piezoelectric materials. The materials
include temperature-dependent fluid dynamic and thermal properties.
36 | CHAPTER 1: INTRODUCTION
2
The goal of this chapter is to give you an insight on how to approach the modeling
of various structural mechanics problems. In this chapter:
Applying Loads
Defining Constraints
Calculating Reaction Forces
Introduction to Material Models
Defining Multiphysics Models
Modeling with Geometric Nonlinearity
Linearized Buckling Analysis
Introduction to Contact Modeling
Eigenfrequency Analysis
Using Modal Superposition
Modeling Damping and Losses
Piezoelectric Losses
Springs and Dampers
37
Tips for Selecting the Correct Solver
Using Perfectly Matched Layers
Loads can be applied in the structural mechanics interfaces on the body, face, edge, or
point levels. Add The Solid Mechanics Interface ( ) to the Model Builder, then
right-click the node to select Body Load, Face Load, Edge Load, or Point Load from the
context menu. This guide includes a detailed description of the functionality for each
physics interface.
In this section:
USING UNITS
Enter loads in any unit, independently of the base SI unit system in the model, because
COMSOL Multiphysics automatically converts any unit to the base SI unit system. To
use the feature for automatic unit conversion, enter the unit in square brackets, for
example, 100[lbf/in^2].
APPLYING LOADS | 39
Custom coordinate systems are also available and are useful, for example, to specify a
load in any direction without splitting it into components. Right-click the Definitions
node ( ) in the Model Tree, to select a Coordinate System from the context menu.
Using Units
See Also
Coordinate Systems
Load Cases
Similar to the familiar concept of load cases, but more powerful, is the parametric
solver. The Parametric feature is available as an attribute to the Stationary solver
operation feature. Under the General section of the Parametric nodes settings window,
name parameters and specify the corresponding list of values. With the parameters
defined, the magnitude, distribution, and location of loads can easily be controlled.
Moreover, the parameters are available for use in any expression of a model.
Figure 2-1 shows a plate with a hole in plane stress loaded with a distributed load and
a point load of the same magnitude. The mesh consists of triangular elements with
quadratic shape functions. The high stress around the point load is dissipated within
the length of a few elements for both mesh cases. The stresses in the middle of the plate
and around the hole are in agreement for the distributed load and the point load. The
problem is that due to the high stress around the singular load it is easy to overlook
the high stress region around the hole. When the point load is applied, the range must
be manually set for the stress plot to get the same visual feedback of the high stress
region around the hole in the two cases. This is because the default plot settings
automatically set the range based on the extreme values of the expression that is
plotted.
Despite these findings it is good modeling practice to avoid singular loads because it is
difficult to estimate the size of the singular region. In the Structural Mechanics
Module it is possible to define loads on all boundary types. However, avoid singular
loads altogether with elastoplastic materials.
APPLYING LOADS | 41
normal mesh size
Figure 2-1: A plate with a hole subject to a distributed load (left) and a point load (right).
Pressure
A pressure is a load acting toward the normal of a face of the structure. If there are
large deformations in the model and the Include geometric nonlinearity check box is
selected under the Study settings section of the current study step, the pressure acts as
a follower load. The pressure is then defined with respect to the geometry and, as the
geometry deforms locally, the orientation of the load changes.
For modeling rotating parts under static conditions, use centrifugal acceleration loads.
The body load in the radial direction is
K r = 2 r (2-1)
where is the density of the material, is the angular frequency, and r is the radial
distance from the axis of rotation. A cylindrical coordinate system is often useful here.
When a separate physics interface is used to model heat transfer in the material, the
entry for the temperature is the dependent variable for the temperature from that
physics interface, typically T. In most cases, possible temperature variables from other
interfaces can be directly selected from a list.
For more information about how to couple heat transfer analysis with
structural mechanics analysis, see Thermal-Structural Interaction. This
Note module also includes The Thermal Stress Interface.
Total Loads
You can specify a load either as a distributed load per unit length, area, or volume, or
as a total force to be uniformly distributed on a boundary.
APPLYING LOADS | 43
Defining Constraints
Defining the proper constraints for structural mechanics models is just as important as
defining the loads as together they make up the model boundary conditions. This
module has many useful predefined features to define the constraints or to create user
defined expressions that define constraints.
In this section:
Orientation
Symmetry Constraints
Kinematic Constraints
Rotational Joints
Orientation
You can specify constraints in global as well as in any previously defined local
coordinate system.
Symmetry Constraints
In many cases symmetry of the geometry and loads can be used to your advantage in
modeling. Symmetries can often greatly reduce the size of a model and hence reduce
the memory requirements and solution time. When a structure exhibits axial
symmetry, use the axisymmetric physics interfaces. A solid that is generated by rotating
a planar shape about an axis is said to have axial symmetry. In order to make use of the
axisymmetric physics interfaces, all loads and constraints must also be the same around
the circumference.
For other types of symmetry, use the predefined symmetry and antisymmetry
constraints. This means that no expressions need to be enteredinstead just add the
type of constraint to apply to the model.
Figure 2-2: If the geometry exhibits two symmetry planes, model a quarter of the geometry
by using the Symmetry feature for the two selected surfaces.
Both geometric symmetry and loads are important when selecting the
correct constraints for a model.
Important
DEFINING CONSTRAINTS | 45
of the object would be Antisymmetry in the case of antisymmetric loading and
Symmetry in the case of symmetric loading of the object.
Kinematic Constraints
Kinematic constraints are equations that control the motion of solids, faces, edges, or
points. Add a Prescribed Displacement constraint to enter expressions for constraints.
You can define the equations using predefined coordinate systems as well as custom
coordinate systems. Special constraints, for instance to keep an edge of body straight
or to make a boundary rotate, require such constraint equations.
Rotational Joints
Joints between elements in The Truss Interface are automatically rotational joints
because the truss elements have no rotational degrees of freedom. For beams, however,
the rotational degrees of freedom are by default coupled between elements. To create
a rotational joint between two beam elements, add one additional 2D or 3D Beam
interface to a geometry. Make sure that it is only active for the boundary that includes
the point where the joint is positioned and that no other physics interface is active here.
Couple the translational degrees of freedom and leave the rotational degrees of
freedom uncoupled at the joint.
Except for the special situations, use the predefined variables as this is the simplest and
most accurate method. While weak constraints give equally accurate results as the
predefined variables, the method requires adding extra degrees of freedom to the
model. Using weak constraints, however, does have the advantage that it is possible to
visualize the reactions as a continuous fields. Evaluating the surface traction on
constrained boundaries, on the other hand, is more direct but provides only
approximate values. The following sections describe these methods.
Specify the coordinates of the Reference Point for Moment Computation at the top level
of the physics interface main nodes settings window. After editing the reference point
coordinates, right-click the study node and select Update Solution for the change to
take effect on the reaction moment calculation.
Reaction forces are not available for eigenfrequency analysis or when weak
constraints are used.
Note
With weak constraints activated, COMSOL Multiphysics adds the reaction forces to
the solution components. The variables are denoted X_lm, where X is the name of the
constrained degree of freedom (as, for example, u_lm and v_lm). The extension lm
stands for Lagrange multipliers. It is only possible to evaluate reaction forces on
constrained boundaries in the constraint directions. To compute the total reaction
force on a boundary, integrate one of the variables X_lm using Volume Integration,
Surface Integration, or Line Integration under Derived Values.
Because the reaction force variables are added to the solution components, the number
of DOFs for the model increases slightly, depending on the mesh size for the
boundaries in question. Boundaries that are adjacent to each other must have the same
constraint settings. The reason for this is that adjacent boundaries share a common
node.
For 2D models, multiply the surface traction by the cross section thickness
before integrating to calculate the total reaction force. Because the surface
traction vector is based on computed stress results, this method is less
2D accurate than solving for the reaction forces using the reaction force
variables or weak constraints.
In this section:
Ex
xy = yx ------
Ey
Mixed Formulation
As described in the Theory for the Solid Mechanics Interface, the negative mean stress
becomes an additional dependent variable when the Nearly incompressible material
check box is selected in the settings window for the material. Select this setting when
the Poissons ratio of a material is close to 0.5, which means that the material is nearly
incompressible. The mixed formulation is useful not only for linear elastic materials but
also for elastoplastic materials, hyperelastic materials, and linear viscoelastic materials.
The Hyperelastic Material and Plasticity features are available with the
Nonlinear Structural Materials Module.
Tip
Not all iterative solvers work together with mixed formulation because the stiffness
matrix becomes indefinite.
Because the shape function for the pressure must be one order less than
the shape functions for the displacements, it is not possible to use linear
elements for the displacement variables on the domains where mixed
Important
formulation is turned on.
In this section:
Thermal-Structural Interaction
Acoustic-Structure Interaction
Thermal-Electric-Structural Interaction
Thermal-Structural Interaction
The Thermal Stress Interface included with this module has a predefined one-way
coupling for thermal-structure interaction (thermal stress), which combines a Solid
Mechanics interface with a Heat Transfer interface from the Heat Transfer Module or
COMSOL Multiphysics.
This approach, using a single iteration, does not produce a correct result
if there are thermal properties that depend on the displacements, for
example, a heat source caused mechanical losses (damping) in the
Note
material.
Acoustic-Structure Interaction
Model acoustic-structure interaction when the Structural Mechanics Module is used
with an acoustics interface from the Acoustics Module. Additional tools are made
available, including a transient solver and the means to simulate absorbing or radiation
boundary conditions.
Coupled acoustic-structure models have symmetric matrices, which means you can
take advantage of the SPOOLES solver to reduce memory requirements. There are
several interfaces available that are documented and described in the Acoustic Module
Users Guide.
Thermal-Electric-Structural Interaction
The Joule Heating and Thermal Expansion Interface enables
thermal-electric-structural interaction to combine with thermal expansion, which is a
one-way coupling that includes a solid mechanics interface and a heat transfer interface
from the Heat Transfer Module or COMSOL Multiphysics, with Joule heating and
temperature-dependent electrical conductivity, which is a two-way coupling that
includes a solid mechanics interface from the Structural Mechanics Module and a heat
transfer interface from the Heat Transfer Module or COMSOL Multiphysics.
By default, COMSOL Multiphysics takes advantage of the one-way coupling and solve
the problem sequentially using the segregated solver: solve for temperature and electric
potential using a coupled approach and then perform the stress-strain analysis using
the computed temperature field from the heat transfer equation.
These settings, using a single iteration, does not produce a correct result
if there are thermal properties that depend on the displacements, making
Note the thermal-structure part into a two-way coupling.
Thin structures, where the deflection is of the same order of magnitude as the
thickness.
Where the structure exhibits large rotations. A rigid body rotation of only a few
degrees causes significant strains and stresses in a material where a linear strain
representation is used.
Where the strains are larger than a few percent.
Contact problems.
Where a prestress must be taken into account for computing the dynamic response
of a structure.
Buckling problems.
Where a deformed mesh is used.
Fluid-structure interaction problems.
In this section:
If any active feature in the model requires the analysis to be geometrically nonlinear,
the Include geometric nonlinearity check box is selected automatically, and it cannot be
cleared. The features which force this behavior are:
Contact, because the deformation between the contacting boundaries must be part
of the solution
The Fluid-Structure Interaction interface.
Usually you would also want to use geometric nonlinearity when a Moving Mesh
interface is present, but this is not forced by the program.
When you select a geometrically nonlinear study step, the behavior of several features
differs from that in a geometrically linear case:
See Also Adding Study Steps in the COMSOL Multiphysics Users Guide
The effect of using geometric nonlinearity in these interfaces is limited to the stress and
strain representation as the other effects described in Geometric Nonlinearity for the
Solid Mechanics Interface are not relevant.
See Also Adding Study Steps in the COMSOL Multiphysics Users Guide
Usually, a study of a prestressed problem uses two study steps. The first step is a
Stationary step in which the static preload is applied. The effects of the preload can be
computed with or without taking geometric nonlinearity into account. In the second
study step, where the you compute the Eigenfrequency or the Frequency Response, it is
necessary to take geometric nonlinearity into account. Even if the displacements and
strains are small, this is what gives the prestress contribution to the equations.
The same principles apply also to a linear buckling analysis, except that both study steps
should be geometrically linear. The nonlinear contribution is included in the
formulation of the buckling eigenvalue itself.
There are three Preset study types which can be used to set up these two
study steps: Prestressed Analysis, Eigenfrequency; Prestressed Analysis,
Tip Frequency Domain; and Linear Buckling.
Figure 2-4: An example of the deformation of a beam showing the undeformed state (left)
and the deformed state (right) of the beam itself with an element near its tip highlighted
(top), of the element (center) and of lines parallel to the x-axis in the undeformed state
(bottom).
For example, in the Eulerian system it is easy to define forces per unit area (known as
tractions) that act within the solid, and to define a stress tensor that represents all of
these forces that act on a volume element. Such forces could be physically measured,
for example using an implanted piezoresistor. The stress tensor in the Eulerian frame
is called the Cauchy or true stress tensor (in COMSOL this is referred to as the spatial
stress tensor). To construct the stress tensor in the Lagrangian frame a tensor
transformation must be performed on the Cauchy stress. This produces the second
Piola-Kirchhoff (or material) stress, which can be used with the Lagrange or material
If the strains are small (significantly less than 10 percent), and there are no
significant rotations involved with the deformation (significantly less than
10 degrees), geometric nonlinearity can be disabled, resulting in a linear
equation system which solves more quickly (Ref. 1). This is often the case
for many practical MEMS structures.
Tip
Geometric nonlinearity can be enabled or disabled within a given model
by changing the Include geometric nonlinearity setting in the relevant
solver step.
In the case of solid mechanics, the material and spatial frames are associated directly
with the Lagrangian and Eulerian frames, respectively. In a more general case (for
example, when tracking the deformation of a fluid, such as a volume of air surrounding
a moving structure) tying the Lagrangian frame to the material frame becomes less
desirable. Fluid must be able to flow both into and out of the computational domain,
without taking the mesh with it. The Arbitrary Lagrangian-Eulerian (ALE)
method allows the material frame to be defined with a more general mapping to the
spatial or Eulerian frame. In COMSOL, a separate equation is solved to produce this
mappingdefined by the mesh smoothing method (Laplacian, Winslow, or
hyperelastic) with boundary conditions that determine the limits of deformation (these
are usually determined by the physics of the system, whilst the domain level equations
are typically being defined for numerical convenience). The ALE method offers
significant advantages since the physical equations describing the system can be solved
in a moving domain.
REFERENCE
1. A. F. Bower, Applied Mechanics of Solids, CRC Press, Boca Raton, FL, 2010
(http://www.solidmechanics.org).
See Also Solver Studies and Study Types in the COMSOL Multiphysics Users
Guide
The idea behind this type of analysis, can be described in the following way:
Ku = K L + K NL u = f
Here the total stiffness matrix, K, has been split into a linear part, KL, and a nonlinear
contribution, KNL.
In a first approximation, KNL is proportional to the stress in the structure and thus to
the external load. So if the linear problem is solved first,
KL u0 = f0
for an arbitrary initial load level f0, then the nonlinear problem can be approximated as
K L + K NL u 0 u = f 0
An instability is reached when this system of equations becomes singular so that the
displacements tend to infinity. The value of the load at which this instability occurs can
be determined by, in a second study step, solving an eigenvalue problem for the load
multiplier .
K L + K NL u 0 u = 0
Be aware that for some structures, the true buckling load may be
significantly smaller that what is computed using a linearized analysis.
This phenomenon is sometimes called imperfection sensitivity. Small
deviations from the theoretical geometrical shape can then have a large
Note
impact on the actual buckling load. This is especially important for curved
shells.
This section has some suggestions regarding important aspects of creating models
involving contact, with or without friction, between parts. Contact is implemented
based on the augmented Lagrangian method, which is described in About Contact
Modeling in the theory section for The Solid Mechanics Interface. When modeling
contact between structural parts, set up contact pairs, which define where the parts
may come into contact. A contact pair consists of two sets of boundaries, which are the
source and destination boundaries. The 2D and 3D structural interfaces use the pairs
to set up equations that prevent the destination boundaries to penetrate the source
boundaries.
In this section:
Constraints
Contact Pairs
Boundary Settings for Contact Pairs
Time-Dependent Analysis
Multiphysics Contact
Solver and Mesh Settings for Contact Modeling
Monitoring the Solution
Constraints
Make sure that the bodies are sufficiently constrained, also in the initial position. If the
bodies are not in contact in the initial configuration, and there are no constraints on
p 10
k x = k 1 p 2
and similarly in any other direction that needs to be restrained. It is not important
whether the spring is applied to domains, boundaries, or edges, but a value for the
parameter k must be chosen so that the generated extra force balances the external load
at a displacement that is of the same order of magnitude as a characteristic element size
in the contact region.
Contact Pairs
For efficiency, only include those boundaries that may actually come in contact in the
destination. For the source, it is often a bit more efficient to make it so large that every
destination point has a corresponding source point. The corresponding source point
is obtained by following the normal to the destination until it reaches the source.
Make sure that the source boundary stiffness in the normal direction is higher than
the destination boundary stiffness. This is especially important if the difference in
stiffness is quite large (for example, over ten times larger). Also keep in mind that
for elastoplastic or hyperelastic materials there can be a significant change in stiffness
during the solution process, and choose the source and destination boundaries
accordingly. For such materials the penalty factor might have to be adjusted as the
solution progresses.
The Hyperelastic Material and Plasticity features are available with the
Nonlinear Structural Materials Module.
Tip
Once the source and destination boundaries are chosen, mesh the destination finer
than the source. Do not make the destination mesh just barely finer than the source
because this often causes nonphysical oscillations in the contact pressure. Make the
destination at least two times finer than the source.
PENA LT Y FACT OR S
When running into convergence problems, check the penalty parameters. If the
iteration process fails in some of the first augmented iterations, lower the penalty
parameters. If the model seems to converge but very slowly, consider increasing the
maximum value of the penalty parameters.
Increasing the penalty factor can lead to an ill-conditioned Jacobian matrix and
convergence problems in the Newton iterations. This is often seen by noting that the
damping factor becomes less than 1 for many Newton iterations. If this occurs,
decrease the penalty factors.
The default values for the penalty factors is based on an equivalent Youngs modulus
Eequiv of the material on the destination side. For linear elastic isotropic materials
Eequiv is equal to the actual Youngs modulus. For other types of materials Eequiv is
defined an estimate of a similar stiffness at zero strain. For elastoplastic and other
nonlinear materials where the stiffness changes with deformation it may be found that
the default value works fine until there is a significant change in stiffness. This can give
rise to convergence problems for the nonlinear solver because the penalty factor
The Hyperelastic Material and Plasticity features are available with the
Nonlinear Structural Materials Module.
Tip
For materials which are user defined or in other senses nonstandard, for example
anisotropic with large differences in stiffness in different directions, Eequiv may need
to be replaced with another estimate.
INITIAL VALUE
In force-controlled contact problems where no other stiffness prohibits the
deformation except the contact, the initial contact pressure is crucial for convergence.
If it is too low the parts might pass through each other in the first iteration. If it is too
high they never come into contact.
Time-Dependent Analysis
The contact formulation is strictly valid only for stationary problems. It is still possible
to use contact modeling also in a time-dependent analysis, as long as inertial effects are
not important in the contact region. In practice, this means that you can solve both
quasi-static problems and truly dynamic problems, as long as situations with impact are
avoided. When in doubt, try to do an a posteriori check of conservation of momentum
and energy to ensure that the solution is acceptable.
Multiphysics Contact
Multiphysics contact problems are often very ill conditioned, which leads to
convergence problems for the nonlinear solver. For example, take heat transfer
through the contact area, where initially only one point is in contact. The solution for
the temperature is extremely sensitive to the size of the contact area (that is, the
problem to determine the temperature is ill conditioned).
It is important to resolve the size of the contact area accurately, that is, to
use a very fine mesh in the contact area.
Important
You can use the contact variables (gap and contact pressure) in expressions for
quantities in other physics interfaces. As an example, a thermal resistance in the contact
region can depend on the contact pressure.
Use a direct solver instead of an iterative solver if the problem size allows it. Direct
solvers are less sensitive and can provide better convergence.
As a default, the Double dogleg nonlinear solver is selected when a stationary study
is generated and Contact features are present in the model. For the majority of
contact problems this solver has more stable convergence properties than the
Newton solver, which is the general default solver. Using the same settings, the
Double dogleg solver tends to be somewhat slower than the Newton solver on
problems where both solvers converge. It is however often possible to user larger
load steps when using the Double dogleg solver. For some problems, the Newton
solver may still be the better choice, so if you experience problems using the default
settings, try to switch solver.
For some contact problems, it may be necessary to let the parametric solver use a
more defensive strategy when going to the next parameter step. This can be
controlled by setting the value of Predictor in the Parametric feature to Constant.
A coarse mesh on a curved contact surface might lead to convergence problems, so
make sure that the mesh is sufficiently fine on the contact surface.
If the model includes friction, solve the problem without friction first. When that
model seems to work without friction, friction can then be added.
Always solve contact problems with friction incrementally using a parametric or
time-dependent solver because the development of friction forces is history
dependent. For contact problems without friction an incremental strategy is not
necessary but often a good choice.
Using the Results while solving functionality in the study step is a good practice. You
can either use a stress plot, or a plot of the contact pressure. If there are more than one
contact pair, superimpose several plots. In most cases, the scale of a deformed plot
should be set to 1 when monitoring contact problems. Note that if you select Results
while solving also in the Segregated feature, the plot will be updated after each iteration,
thus allowing you to monitor the convergence in detail.
For each contact pair, two global variables which can be used in probe plots are
available. These are the maximum contact pressure (solid.Tnmax_pair) and the
minimum gap distance (solid.gapmin_pair).
f = ---------
2i
Because only the shape and not the size of the modes (eigenvectors) have physical
significance, the computed modes can be scaled arbitrarily. Select the method for
scaling in the Eigenvalue Solver node of the solver sequence. If output of mass
participation factors is required, then Scaling of eigenvectors must be set to Mass matrix.
This means that the eigenmodes U are orthogonalized with respect to the mass matrix
M so that
T
U i MU i = 1 (2-2)
This is a common choice for the scaling of eigenvectors within the structural mechanics
field.
The mass (or modal) participation factor for mode i in direction j, rij, is defined as
T
r ij = U i Md j
The mass participation factors have the important property that when their squares for
a certain direction are summed over all modes, this sum approaches the total mass of
the model:
EIGENFREQUENCY ANALYSIS | 69
n
rij =
2
m tot
i
In practice you seldom solve for all possible modes but just a limited number. Then
this property can be used for investigating how well a certain number of selected
modes represent the total mass of the system.
The mass participation factors are available as a global variables, and these can be
shown in a table using a Global Evaluation node under Derived Values in the Results
branch, for example. The participation factor variables are available as predefined
variables in the Solver submenu.
Eigenvalue Solver
See Also
Study Types
Defining Derived Values and Tables
The analysis is linear. It is however possible that the structure has been subjected to
a preceding non-linear history. The modal response can then be a linear
perturbation around that state.
There are no non-zero prescribed displacements.
The important frequency content of the load is limited to a range which is small
when compared to all the eigenfrequencies of the model, so that its response can be
approximated with a small number of eigenmodes. In practice, this excludes wave
and shock type problems.
If the modal solution is performed in the time domain, all loads must have the same
dependency on the time.
When using the Structural Mechanics Module, there are two predefined study types
for modal superposition: Time-Dependent Modal and Frequency-Domain Modal. Both
these study types consist of two study steps: One step for computing the
eigenfrequencies and one step for the modal response.
In practice, you will often have computed the eigenfrequencies already, and then want
to use them in a modal superposition. In this case, start by adding an empty study, and
then add a Time-Dependent Modal or Frequency-Domain Modal study step to it. After
having added the study step this way, you must point the modal solver to the solution
containing the eigenfrequencies and eigenmodes. You do this by first selecting Show
default solver at the study level, and then selecting the eigenfrequency solution to be
used in the Eigenpairs section of the generated modal solver.
2
d u du
m ---------- + c ------- + ku = f t (2-3)
dt dt
In this equation u is the displacement of the degree of freedom, m is its mass, c is the
damping parameter, and k is the stiffness of the system. The time (t) dependent forcing
term is f(t). This equation is often written in the form:
where c2m0 and 02km. In this case is the damping ratio (1 for critical
damping) and 0 is the resonant frequency of the system. In the literature it is more
common to give values of than c. can also be readily related to many of the various
measures of damping employed in different disciplines. These are summarized in
Table 2-1.
TABLE 2-1: RELATIONSHIPS BETWEEN MEASURES OF DAMPING
Logarithmic u t0 d 2
decrement d = ln ---------------------- 1
u t0 +
where t0 is a reference time and is the
period of vibration for a decaying, unforced
degree of freedom.
Quality factor Q = Q 1 2
where is the bandwidth of the amplitude 1
resonance measured at 1 2 of its peak.
Loss factor 1 Qh At the resonant
= ------ --------
2 W h frequency:
where Qh is the energy lost per cycle and Wh 2
is the maximum potential energy stored in the 1
cycle. The variables Qh and Wh are available
in COMSOL as: solid.Qh and solid.Wh.
In the frequency domain, the time dependence of the force and the displacement can
be represented by introducing a complex force term and assuming a similar time
dependence for the displacement. The equations
jt jt
f t = Re Fe and u t = Re Ue
are written where is the angular frequency and the amplitude terms U and F can in
general be complex (the arguments provide information on the relative phase of
signals). Usually the real part is taken as implicit and is dropped subsequently.
Equation 2-3 takes the following form in the frequency domain:
2
mU + jcU + kU = F (2-5)
2 2 F
U + j 0 U + 0 U = ----- (2-6)
m
There are two basic damping models available in COMSOL Rayleigh damping and
models based on introducing complex quantities into the equation system.
Rayleigh Damping introduces damping in a form based on Equation 2-3. This means
that the method can be applied generally in either the time or frequency domain. The
parameter c in Equation 2-3 is defined as a fraction of the mass and the stiffness using
two parameters, dM and dK, such that
c = dM m + dK k (2-7)
Although this approach seems cumbersome with a one degree of freedom system,
when there are many degrees of freedom m, k, and c become matrices and the
technique can be generalized. Substituting this relationship into Equation 2-3 and
rearranging into the form of Equation 2-4 gives:
2
d u- 2 du 2 ft
--------- + dM + dK 0 ------- + 0 u = ---------
dt dt m
1 dM
= --- ----------- + dK 0 (2-8)
2 0
Note that Equation 2-8 holds separately for each vibrational mode in the system at its
resonant frequency. In the frequency domain it is to use frequency dependent values
of dM and dK. For example setting dM0 and dK2/0 produces a Equivalent
Viscous Damping model at the resonant frequency.
While Rayleigh damping is numerically convenient, the model does not agree with
experimental results for the frequency dependence of material damping over an
extended range of frequencies. This is because the material damping forces behave
more like frictional forces (which are frequency independent) than viscous damping
forces (which increase linearly with frequency as implied by Equation 2-5). In the
frequency domain it is possible to introduce loss factor damping, which has the desired
property of frequency independence.
2
mU + jkU + kU = F (2-9)
2 2 2 F
mU + j 0 U + 0 U = -----
m
Comparing these equations with Equation 2-5 and Equation 2-6 shows that the loss
factor is related to and c by:
= 2 ------ = ---- c
0 k
Equation 2-9 shows that the loss factor has the desired property of frequency
independence. However it is clear that this type of damping cannot be applied in the
time domain. In addition to using loss factor damping the material properties can be
entered directly as complex values in COMSOL, which results in Explicit Damping
Piezoelectric Losses are more complex and include coupling and electrical
losses in addition to the material terms.
Rayleigh Damping
As discussed for a model with a single degree of freedom, the Rayleigh damping model
defines the damping parameter c in terms of the mass m and the stiffness k as
c = dM m + dK k
where dM and dK are the mass and stiffness damping parameters, respectively. At any
resonant frequency, f, this corresponds to a damping factor, given by:
1 dM
= --- ----------- + dK 2f (2-10)
2 2f
Using this relationship at two resonant frequencies f1 and f2 with different damping
factors 1 and 2 results in an equation system
1 -
----------- f
4f 1 1 dM 1
=
1 - dK 2
----------- f
4f 2 2
As a result of its non-physical nature, the Rayleigh damping model can only be tuned
to give the correct damping at two independent resonant frequencies or to give an
approximately frequency independent damping response (which is physically what is
usually observed) over a limited range of frequencies.
Specified damping
f1 f2
f
For many applications it is sufficient to leave dM as zero (the default value) and to
define damping only using the dK coefficient. Then according to Equation 2-10
linearly increasing damping is obtained. If the damping ratio f0 or loss factor f0
is known at a given frequency f0, the appropriate value for dK is:
dK = f 0 = 2f 0
This model results in a well-defined, linearly increasing damping term that has the
defined value at the given frequency.
All physics interfaces under the Structural Mechanics branch use zero
default values (that is, no damping) for dM and dK. These default values
Tip must be changed to meet the specific modeling situation.
G = G + jG = 1 + j s G
where G' is the storage modulus, G'' is the loss modulus, and their ratio sG''G' is
the loss factor. The term G' defines the amount of stored energy for the applied strain,
whereas G'' defines the amount of energy dissipated as heat; G', G'', and s can all be
frequency dependent.
= D 1 + j s th 0 + 0 .
For hyperelastic materials the loss information appears as a multiplier in strain energy
density, and thus in the second Piola-Kirchhoff stress, S:
W s
S = 1 + j s ----------
E
Loss factor damping is available for frequency response analysis and damped
eigenfrequency analysis in all interfaces, but it is not defined for elastoplastic materials
Within a piezoelectric there is a coupling between the strain and the electric field,
which is determined by the constitutive relation:
T
S = sE T + d E
(2-11)
D = dT + T E
Here, S is the strain, T is the stress, E is the electric field, and D is the electric
displacement field. The material parameters sE, d, and T, correspond to the material
compliance, the coupling properties and the permittivity. These quantities are tensors
of rank 4, 3, and 2 respectively, but, since the tensors are highly symmetric for physical
reasons, they can be represented as matrices within an abbreviated subscript notation,
which is usually more convenient. In the Piezoelectric Devices interface, the Voigt
notation is used, which is standard in the literature for piezoelectricity but which differs
from the defaults in the Solid Mechanics interface. Equation 2-11 is known as the
strain-charge form of the constitutive relations. The equation can be re-arranged into
the stress-charge form, which relates the material stresses to the electric field:
PIEZOELECTRIC LOSSES | 81
T
T = cE S e E
(2-12)
D = eS + S E
The material properties, cE, e, and S are related to sE, d, and T. Note that in
COMSOL it is possible to use either form of the constitutive relations. In addition to
Equation 2-11 or Equation 2-12, the equations of solid mechanics and electrostatics
must also be solved within the material.
The stiffness, compliance, coupling, and dielectric material property matrices are
defined in COMSOL with the crystal axes aligned with the local coordinate axes. In
the absence of a user defined coordinate system, the local system corresponds to the
global X, Y, and Z coordinate axes.
The material properties are defined in the material frame, so that if the
solid rotates during deformation the material properties rotate with the
Note solid. See Modeling with Geometric Nonlinearity.
The crystal axes used to define material properties in COMSOL correspond to the
1987 IEEE standard. All piezoelectric material properties are defined using the Voigt
For some materials, the crystal X, Y, and Z axes are defined differently
between the 1987 IEEE standard and the 1949 I.R.E. standard.
Figure 2-6 shows the case of right-handed quartz (which is included in
the COMSOL material library as quartz; see Piezoelectric Materials
Database), which has different axes defined within the two standards.
Note The different axes sets result in different material properties so, for
example, the elasticity or stiffness matrix component cE14 of quartz takes
the value 18 GPa in the 1987 standard and 18 GPa in the 1949
standard.
The crystal cuts are also defined differently within the 1949 and 1987 standards. Both
standards use a notation that defines the orientation of a virtual slice (the plate)
through the crystal. The crystal axes are denoted X, Y, and Z and the plate, which is
usually rectangular, is defined as having sides l, w, and t (length, width, and thickness).
Initially the plate is aligned with respect to the crystal axes and then up to three
rotations are defined, using a right-handed convention about axes embedded along the
l, w, and t sides of the plate. Taking AT cut quartz as an example, the 1987 standard
defines the cut as: (YXl) 35.25. The first two letters in the bracketed expression
always refer to the initial orientation of the thickness and the length of the plate.
Subsequent bracketed letters then define up to three rotational axes, which move with
the plate as it is rotated. Angles of rotation about these axes are specified after the
bracketed expression in the order of the letters, using a right-handed convention. For
AT cut quartz only one rotation, about the l axis, is required. This is illustrated in
Figure 2-7. Note that within the 1949 convention AT cut quartz is denoted as: (YXl)
35.25, since the X-axis rotated by 180 in this convention and positive angles
therefore correspond to the opposite direction of rotation (see Figure 2-6).
PIEZOELECTRIC LOSSES | 83
Figure 2-6: Crystallographic axes defined for right-handed quartz in COMSOL and the
1987 IEEE standard (color). The 1949 standard axes are shown for comparison (gray).
Figure 2-7: Definition of the AT cut of quartz within the IEEE 1987 standard. The AT
cut is defined as: (YXl) 35.25. The first two bracketed letters specify the initial
orientation of the plate, with the thickness direction, t, along the crystal Y axis and the
length direction, l, along the X axis. Then up to three rotations about axes that move with
the plate are specified by the corresponding bracketed letters and the subsequent angles. In
this case only one rotation is required about the l axis, of 35.25 (in a right-handed
sense).
PIEZOELECTRIC LOSSES | 85
coordinate system in this case. Figure 2-9 shows how to define the local system such
that the plate has its normal parallel to the global Z axis.
Usually it is most convenient to define the local coordinates with a Rotated System
node, which defines three Euler angles according to the ZXZ convention (rotation
about Z, then X, then Z again). Note that these Euler angles define the local (crystal)
axes with respect to the global axesthis is distinct from the approach of defining the
cut (global) axes with respect to the crystal (local) axes.
Figure 2-8: Defining an AT cut crystal plate within COMSOL, with normal in the global
Y-direction. Within the 1987 IEEE standard the AT cut is defined as (YXl) 35.25.
Start with the plate normal or thickness in the Ycr direction (a) and rotate the plate
Figure 2-9: Defining an AT cut crystal plate within COMSOL, with normal in the global
Z-direction. Within the 1987 IEEE standard the AT cut is defined as (YXl) 35.25.
Begin with the plate normal in the Zcr-direction, so the crystal and global systems are
coincident. Rotate the plate so that its thickness points in the Ycr-direction (the starting
point for the IEEE definition), the global system rotates with the plate (b). Rotate the plate
35.25 about the l axis (d). Finally rotate the entire system so that the global coordinate
system is orientated as it appears in COMSOL (d). The local coordinate system should be
PIEZOELECTRIC LOSSES | 87
defined with the Euler angles (ZXZ: 0, -54.75, 0). (e) shows a coordinate system for this
system in COMSOL.
Piezoelectric Losses
Losses in piezoelectric materials can be generated both mechanically and electrically.
In the time domain, material damping can be added using the Rayleigh Damping
option in the Damping and Loss node. Electrical damping is currently not available in
the time domain.
Rayleigh Damping
HYSTERETIC LOSS
In the frequency domain, the dissipative behavior of the material can be modeled using
complex-valued material properties, irrespective of the loss mechanism. Such hysteretic
losses can be applied to model both electrical and mechanical losses. For the case of
piezoelectric materials, this means that the constitutive equations are written as
follows:
T
= s E + d E
D = d + E
T
where c E , d , and are complex-valued matrices, where the imaginary part defines the
dissipative function of the material.
Both the real and complex parts of the material data must be defined so as to respect
the symmetry properties of the material being modeled and with restrictions imposed
by the laws of physics.
A key requirement is that the dissipation density is positive; that is, there
is no power gain from the passive material. This requirement sets rules for
the relative magnitudes for all material parameters. This is important
Important
when defining the coupling losses.
X = real X
X = imag X real X
the ratio of the imaginary and real part, and the complex data is then:
X = X 1 j X
where the sign depends on the material property used. The loss factors are specific to
the material property, and thus they are named according to the property they refer to,
for example, cE. For a structural material without coupling, simply use s, the
structural loss factor.
PIEZOELECTRIC LOSSES | 89
The Piezoelectric Devices interface defines the loss factors such that a positive loss
factor corresponds to a positive loss. The complex-valued data is then based on sign
rules. For piezoelectric materials, the following equations apply (m and n refer to
elements of each matrix):
m n m n m n
c E = c E 1 + j cE
m n m n m n
e =e 1 j e
m n m n m n
S = S 1 j S
(2-13)
m n m n m n
s E = s E 1 j sE
m n m n m n
d =d 1 j d
m n m n m n
T = T 1 j T
The losses for non-piezoelectric materials are easier to define. Again, using the
complex stiffness and permittivity, the following equations describe the material:
m n m n m n
D = 1 + j D
(2-14)
m n m n m n
e = 1 j e e
Jp
D = r 0 E j ------
(the actual displacement variables do not contain any conductivity effects) or in the
total current expression J = Jd Jp where Jp = eE is the conductivity current and
Jd is the electric displacement current.
No Damping
By default, there is no damping until a Damping or Damping and Loss node is added. In
The Piezoelectric Devices Interface an undamped model can be created by selecting
No damping from the Damping type list in the settings window for the Damping and Loss
node.
PIEZOELECTRIC LOSSES | 91
5. P.C.Y. Lee, N.H. Liu, and A. Ballato, Thickness Vibrations of a Piezoelectric Plate
With Dissipation, IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency
Control, vol. 51, no. 1, January 2004.
6. P.C.Y. Lee and N.H. Liu, Plane Harmonic Waves in an Infinite Piezoelectric Plate
With Dissipation, Frequency Control Symposium and PDA Exhibition, pp. 162
169, IEEE International, 2002.
10. B. A. Auld, Acoustic Fields and Waves in Solids, Krieger Publishing Company,
1990.
The features are completely analogous, with the difference that a Spring Foundation
feature connects the structural part on which it is acting to a fixed ground, while the
Thin Elastic Layer acts between two parts, either on an internal boundary or on a pair.
Spring Data
Loss Factor Damping
Viscous Damping
SPRING DATA
The elastic properties can be defined either by a spring constant or by a force as
function of displacement. The force as a function of displacement may be more
convenient for nonlinear springs. Each spring feature has three displacement variables
defined, which can be used to describe the dependency on deformation. These
variables are named uspring1_tag, uspring2_tag, and uspring3_tag for the three
directions given by the local coordinate system. In the variable names, tag represents
the tag of the feature defining the variable The tag could for example be spf1 or tel1
for a Spring Foundation or a Thin Elastic Layer respectively. These variables measure the
relative extension of the spring after subtraction of any predeformation.
f sl = 1 + i f s
Loss factor damping is only applicable in for eigenfrequency and frequency domain
analysis. In time dependent analysis the loss factor is ignored.
The Spring Foundation features is most commonly used for simulating boundary
conditions with a certain flexibility, such as the soil surrounding a construction. An
other important use is for stabilizing parts that would otherwise have a rigid body
singularity. This is a common problem in contact modeling before an assembly has
actually settled. In this case a Spring Foundation acting on the entire domain is useful
because it avoids the introduction of local forces.
A Thin Elastic Layer between used as a pair condition can be used to simulate thin
layers with material properties which differ significantly from the surrounding
domains. Common applications are gaskets and adhesives.
When a Thin Elastic Layer is applied on an internal boundary, it usually simulates one
local flexibility, such as a fracture zone in a geological model.
In this section:
Symmetric Matrices
Selecting Iterative Solvers
Specifying Tolerances and Scaling for the Solution Components
Symmetric Matrices
The Matrix symmetry list displays in the General section of the settings window for the
Advanced subnode under a solver node such as Stationary Solver. Specify if the
assembled matrices (stiffness matrix, mass matrix) resulting from the compiled
equations are symmetric or not.
Multiphysics models solving for several physics simultaneously, for example, heat
transfer and structural mechanics. Solving for several structural mechanics physics
interfaces, such as shells combined with beams, does not create unsymmetric
matrices.
Linear viscoelastic materials
Elastoplastic analysis
One of the benefits of using the symmetric solvers is that they use less memory and are
faster. The default option is Automatic, which means the solver automatically detects if
the system is symmetric or not. Some solvers do not support symmetric matrices and
always solve the full system regardless of symmetry.
For stationary and time-dependent studies, use the GMRES iterative solver with
geometric multigrid (GMG) as the preconditioner.
Check the mesh quality when using the GMG preconditioner. It does not work well
when using the option to scale the geometry before meshing. When using extruded
meshes, the mesh cases might need to be created manually.
The default scaling for the displacement components is based on the size of the
geometry in the model, and certain reasonable scales are used for the pressure and
contact force variables, if any. You are encouraged to change these scales as soon as
better values are known or can be guessed or estimated from the applied forces, yield
stress, reaction forces, maximum von Mises stress. The same suggestion applies to the
displacement scale, which can be estimated easily if the problem is displacement
controlled. This approach can significantly improve the robustness of the solution. The
scales need to be entered using the main unit system within the model.
In this section:
PML Implementation
Known Issues When Modeling Using PMLs
PML Implementation
This PML implementation uses the following coordinate transform for the general
coordinate variable t.
t n
t' = ------- 1 i F (2-15)
w
The coordinate, t, and the width of the PML region, w, are geometrical parameters
that are automatically extracted for each region. The other parameters are the PML
scaling factor F and the PML order n that can be modified in the PML feature (both
default to unity). The software automatically computes the value for w and the
orientation of the transform for PML regions that are Cartesian, cylindrical, or
spherical.
Figure 2-10: A cube surrounded by typical PML regions of the Cartesian type.
A separate Perfectly Matched Layers node must be used for each isolated PML
domain. That is, to use one and the same Perfectly Matched Layers node, all PML
domains must be in contact with each other. Otherwise the PMLs do not work
properly.
The coordinate scaling resulting from PMLs also yields an equivalent scaling of the
mesh that may effectively result in a poor element quality. (The element quality
displayed by the mesh statistics does not account for this effect.) This typically
happens when the geometrical thickness of the PML deviates much from one
wavelength (local wavelength rather than free space wavelength). The poor element
quality causes poor convergence for iterative solvers and make the problem
ill-conditioned in general.
The expressions resulting from the stretching get quite complicated for spherical
PMLs in 3D. This increases the time for the assembly stage in the solution process.
After the assembly, the computation time and memory consumption is comparable
to a problem without PMLs. The number of iterations for iterative solvers might
increase if the PML regions have a coarse mesh.
PML regions deviating significantly from the typical configurations shown in the
beginning of this section can cause the automatic calculation of the PML parameter
Also use parts of the shapes shown, but the PML scaling does probably not work for
complex shapes that deviate significantly from these shapes.
Solid Mechanics
This chapter describes the Solid Mechanics interface, which is found under the
Structural Mechanics branch ( ) in the Model Wizard. Solid mechanics in this
context means that no simplifications are available and that you solve for the
displacements without involving rotations.
In this chapter:
103
Solid Mechanics Geometry and
Structural Mechanics Physics Symbols
The Solid Mechanics interface in the Structural Mechanics Module is available for
these space dimensions, which are described in this section:
3D Solid Geometry
2D Geometry (plane stress and plane strain)
Axisymmetric Geometry
There are also physics symbols available with structural mechanics features as described
in these sections:
3D Solid Geometry
The degrees of freedom (dependent variables) in 3D are the global displacements u, v,
and w in the global x, y, and z directions, respectively, and the pressure help variable
(used only if a nearly incompressible material is selected), and the viscoelastic strains
(used only for viscoelastic materials).
Figure 3-1: Loads and constraints applied to a 3D solid using the Solid Mechanics
interface.
PLANE STRESS
The plane stress variant of the 2D interface is useful for analyzing thin in-plane loaded
plates. For a state of plane stress, the out-of-plane components of the stress tensor are
zero.
Figure 3-2: Plane stress models plates where the loads are only in the plane; it does not
include any out-of-plane stress components.
The 2D interface for plane stress allows loads in the x and y directions, and it assumes
that these are constant throughout the materials thickness, which can vary with x and
y. The plane stress condition prevails in a thin flat plate in the xy-plane loaded only in
its own plane and without any z direction restraint.
PLANE STRAIN
The plane strain variant of the 2D interface that assumes that all out-of-plane strain
components of the total strain z, yz, and xz are zero.
Loads in the x and y directions are allowed. The loads are assumed to be constant
throughout the thickness of the material, but the thickness can vary with x and y. The
Axisymmetric Geometry
The axisymmetric variant of the Solid Mechanics interface uses cylindrical coordinates r,
(phi), and z. Loads are independent of , and the axisymmetric variant of the
interface allows loads only in the r and z directions.
Antisymmetry1 Antisymmetry
No Rotation1 No Rotation
Pinned1 Pinned
Roller Roller
1
Requires the Structural Mechanics Module
In both cases, the coordinate directions describe the direction in which a positive value
acts. The coordinate direction symbols do not change with the values actually entered
for the data.
Physics symbols are displayed even if no data values have been entered in the feature.
In some cases a single feature can display more than one feature. An example is the
Point Load feature in the Beam interface, which can display either the Point Force
symbol ( ), the Point Moment symbol ( ), or both, depending on the data
actually entered.
The physics symbols also display for any multiphysics interface that
includes Structural Mechanics feature nodes.
Note
3 Add any of the feature nodes listed in Table 3-1 to the interface. Availability is based
on license and interface.
4 When adding the boundary, edge, or point (a geometric entity) to the Selection list
in the feature settings window, the symbol displays in the Graphics window. See
Figure 3-5.
Figure 3-5: Example of the Boundary Load physics symbols as displayed in the COMSOL
Multiphysics Model Library model Deformation of a Feeder Clamp.
Figure 3-6: Example of Roller and Boundary Load physics symbols as displayed in the
COMSOL Multiphysics Model Library model Tapered Membrane End Load.
When this interface is added, these default nodes are also added to the Model Builder
Linear Elastic Material, Free (a boundary condition where boundaries are free, with no
loads or constraints), and Initial Values. Right-click the Solid Mechanics node to add
nodes that implement other solid mechanics features.
Model
Eigenvalue Analysis of a Crankshaft: Model Library path
COMSOL_Multiphysics/Structural_Mechanics/crankshaft
INTERFACE IDENTIFIER
The interface identifier is a text string that can be used to reference the respective
physics interface if appropriate. Such situations could occur when coupling this
interface to another physics interface, or when trying to identify and use variables
defined by this physics interface, which is used to reach the fields and variables in
expressions, for example. It can be changed to any unique string in the Identifier field.
The default identifier (for the first interface in the model) is solid.
DOMAIN SELECTION
The default setting is to include All domains in the model to define the displacements
and the equations that describe the solid mechanics. To choose specific domains, select
Manual from the Selection list.
From the 2D approximation list select Plane stress or Plane strain (the
default). For more information see the theory section.
When modeling using plane stress, the Solid Mechanics interface solves for
2D
the out-of-plane strain, ez, in addition to the displacement field u.
THICKNESS
S T R U C T U R A L TR A N S I E N T B E H AV I O R
From the Structural transient behavior list, select Quasi-static or Include inertial terms to
treat the elastic behavior as quasi-static (with no mass effects; that is, no second-order
time derivatives) or as a mechanical wave in a time-dependent study. The default is to
include the inertial terms to model the mechanical wave.
TY P I C A L WA V E S P E E D
The typical wave speed cref is a parameter for the perfectly matched layers (PMLs) if
used in a solid wave propagation model. The default value is solid.cp, the
pressure-wave speed. If you want to use another wave speed, enter a value or
expression in the Typical wave speed for perfectly matched layers field.
DISCRETIZATION
To display this section, click the Show button ( ) and select Discretization. Select a
Displacement fieldLinear, Quadratic (the default), Cubic, Quartic, or (in 2D) Quintic.
Specify the Value type when using splitting of complex variablesReal or Complex (the
default).
Domain, Boundary, Edge, Point, and Pair Features for the Solid
Mechanics Interface
The Solid Mechanics Interface has these domain, boundary, edge, point, and pair
features listed in alphabetical order. The list also includes subfeatures.
Added Mass
Applied Force
Applied Moment
Antisymmetry
Body Load
Boundary Load
Change Thickness
Contact
Damping
Edge Load
Fixed Constraint
Also right-click the Linear Elastic Material node to add Thermal Expansion
and Initial Stress and Strain subnodes.
DOMAIN SELECTION
From the Selection list, choose the domains to define a linear elastic solid and compute
the displacements, stresses, and strains.
MODEL INPUTS
Define model inputs, for example, the temperature field of the material uses a
temperature-dependent material property. If no model inputs are required, this section
is empty.
Solid Model
Select a linear elastic Solid modelIsotropic, Orthotropic, or Anisotropic. Select:
Isotropic for a linear elastic material that has the same properties in all directions.
Orthotropic for a linear elastic material that has different material properties in
orthogonal directions, so that its stiffness depends on the properties Ei, ij, and Gij.
Anisotropic for a linear elastic material that has different material properties in
different directions, and the stiffness comes from the symmetric elasticity matrix, D.
Youngs modulus and Poissons ratio to specify Youngs modulus (elastic modulus)
E (SI unit: Pa) and Poissons ratio (dimensionless). For an isotropic material
Youngs modulus is the spring stiffness in Hookes law, which in 1D form is E
where is the stress and is the strain. Poissons ratio defines the normal strain in
the perpendicular direction, generated from a normal strain in the other direction
and follows the equation = .
Youngs modulus and shear modulus to specify Youngs modulus (elastic modulus)
E (SI unit: Pa) and the shear modulus G (SI unit: Pa). For an isotropic material
Youngs modulus is the spring stiffness in Hookes law, which in 1D form is E
This is the wave speed for a solid continuum. In plane stress, for example,
the actual speed with which a longitudinal wave travels is lower than the
Note value given.
For each pair of properties, select from the applicable list to use the value From material
or enter a User defined value or expression.
Each of these pairs define the elastic properties and it is possible to convert
from one set of properties to another (see Table 3-2).
Tip
Density
The default Density (SI unit: kg/m3) uses values From material. If User defined is
selected, enter another value or expression.
GEOMETRIC NONLINEARITY
In this section there is always one check box. Either Force linear strains or Include
geometric nonlinearity is shown.
If a study step is geometrically nonlinear, the default behavior is to use a large strain
formulation in all domains. There are however some cases when you would still want
to use a small strain formulation for a certain domain. In those cases, select the Force
linear strains check box. When selected, a small strain formulation is always used,
independently of the setting in the study step. The default value is that the check box
is cleared, except when opening a model created in a version prior to 4.2a. In this case
the state is chosen so that the properties of the model are conserved.
The Include geometric nonlinearity check box is displayed only if the model was created
in a version prior to 4.2a, and geometric nonlinearity was originally used for the
selected domains. It is then selected and forces the Include geometric nonlinearity check
box in the study step to be selected. If the check box is cleared, the check box is
ENERGY DISSIPATION
To display this section, click the Show button ( ) and select Advanced Physics Options.
Select the Calculate dissipated creep energy check box as required.
This section is available when you also have the Geomechanics Module or
the Nonlinear Structural Materials Module. The Creep node and links to
the theory about this section is described in the Geomechanics Module
Note Users Guide and the Nonlinear Structural Materials Module Users
Guide, respectively.
Change Thickness
The Change Thickness feature is available in 2D for the Solid Mechanics and
Plate interfaces. It is available in 3D for the Membrane interface on
2D boundaries instead of domains.
Use the Change Thickness feature to model domains with a thickness other than the
overall thickness defined in the physics interfaces Thickness section.
DOMAIN SELECTION
From the Selection list, choose the domains use a different thickness.
CHANGE THICKNESS
Enter a value for the Thickness d (SI unit: m). This value replaces the overall thickness
for the domains selected above.
When you have the Nonlinear Structural Materials Module, you can also
add Damping to the Hyperelastic Material feature, which adds loss-factor
Note damping.
DOMAIN SELECTION
From the Selection list, choose the domains to add damping.
By default, this node inherits the selection from its parent node, and only
a selection that is a subset of the parent nodes selection can be used.
Note
DAMPING SETTINGS
Select a Damping typeRayleigh damping (the default), Isotropic loss factor or
Anisotropic loss factor.
In the Rayleigh damping model, the damping parameter is expressed in terms of the
mass m and the stiffness k as
= dM m + dK k
That is, Rayleigh damping is proportional to a linear combination of the stiffness and
mass, There is no direct physical interpretation of the mass damping parameter dM
(SI unit: 1/s) and the stiffness damping parameter dM (SI unit: s).
The use of loss factor damping traditionally refers to an scalar-valued loss factor s. But
there is no reasonthat s must be a scalar. Because the loss factor is a value deduced
from true complex-valued material data, it can be represented by a matrix of the same
dimensions as the dimensions of the anisotropic stiffness matrix. Especially for
orthotropic material, there should set of loss factors of all normal and shear elasticity
modulus components. the following loss-elasticity combinations are available:
Isotropic Loss Factor From the Isotropic structural loss factor list, the default uses
values From material. If User defined is selected, enter another value or expression. The
default value is 0.
Orthotropic Loss Factor From the Loss factor for orthotropic Youngs modulus list,
select From material (the default) to use the material value or select User defined and
enter values for E.
From the Loss factor for orthotropic shear modulus list, select From material (the
default) to use the material value or select User defined and enter values for G or
GVo.The values for the shear modulus loss factors are ordered in two different ways,
consistent with the selection of either Standard (XX, YY, ZZ, XY, YZ, XZ) or Voigt (XX,
YY, ZZ, YZ, XZ, XY) notation in the corresponding Linear Elastic Model. The default
values are 0. If the values are taken from the material, these loss factors are found in
the Orthotropic or Orthotropic, Voigt notation property group for the material.
Anisotropic Loss Factor From the Loss factor for elasticity matrix D list, the default uses
values From material. If User defined is selected, enter other values or expressions for D
or DVo. Select Symmetric to enter the components of D in the upper-triangular part
of a symmetric 6-by-6 matrix or Isotropic to enter a single scalar loss factor. The values
for the loss factors are ordered in two different ways, consistent with the selection of
either Standard (XX, YY, ZZ, XY, YZ, XZ) or Voigt (XX, YY, ZZ, YZ, XZ, XY) notation
is used in the corresponding Linear Elastic Model. The default values are 0. If the
values are taken from the material, these loss factors are found in the Anisotropic or
Anisotropic, Voigt notation property group for the material.
DOMAIN SELECTION
From the Selection list, choose the domains to define an initial value.
IN IT IA L VA LUES
Enter values or expressions for the initial values of the Displacement field u (SI unit: m)
(the displacement components u, v, and w in 3D) and Velocity field ut (SI unit: m/
s). The defaults are 0.
For all of these loads, right-click and choose Phase to add a phase for
harmonic loads in frequency-domain computations.
Tip
Body Load
Add a Body Load to domains for gravity effects, for example. Right-click to add a Phase
for harmonic loads in frequency-domain computations.
FORCE
Select a Load type for 2D models and then enter values or expressions for
the components in the matrix:
Load defined as force per unit area FA (SI unit: N/m2). The body load
as force per unit volume is then the value of F divided by the thickness.
2D Load defined as force per unit volume FV (the default) (SI unit: N/m3)
Total force Ftot (SI unit: N). COMSOL then divides the total force by
the volume of the domains where the body load is active.
Select a Load type for 3D models and then enter values or expressions for
the components in the matrix:
Load defined as force per unit volume FV (the default) (SI unit: N/m3)
3D Total force Ftot (SI unit: N). COMSOL then divides the total force by
the volume of the domains where the load is active.
Boundary Load
Add a Boundary Load to boundaries for a pressure acting on a boundary, for example.
Right-click to add a Phase for harmonic loads in frequency-domain computations.
BOUNDARY SELECTION
From the Selection list, choose the boundaries to define a load.
PAIR SELECTION
If Boundary Load is selected from the Pairs menu, choose the pair to define. An identity
pair has to be created first. Ctrl-click to deselect.
FORCE
After selecting a Load type, the Load list normally only contains User
defined. When combining the Solid Mechanics interface with, for example,
Note film damping, it is also possible to choose a predefined load from this list.
Select a Load type for 2D models and then enter a value or expression in
the matrix:
Select a Load type for 3D models and then enter a value or expression in
the matrix:
EDGE SELECTION
From the Selection list, choose the edges to define an edge load.
FORCE
Select a Load type for 3D models and then enter values or expressions for
the components in the matrix:
Load defined as force per unit length FL (SI unit: N/m). When
combining the Solid Mechanics interface with, for example, film
3D damping, it is also possible to choose a predefined load from this list.
Total force Ftot (SI unit: N). COMSOL then divides the total force by
the volume where the load is active.
Point Load
Add a Point Load to points for concentrated forces at points. Right-click to add a Phase
for harmonic loads in frequency-domain computations.
POINT SELECTION
From the Selection list, choose the points to define a point load.
FORCE
Enter values or expressions for the components of the Point load Fp (SI unit: N).
BOUNDARY SELECTION
For the Shell and Membrane interfaces, the Free feature is applied to edges.
For the Beam and Truss interfaces, it is applied to points.
Note
From the Selection list, choose the boundaries that are free.
PAIR SELECTION
If Free is selected from the Pairs menu, choose the pair to define. An identity pair has
to be created first. Ctrl-click to deselect.
Fixed Constraint
The Fixed Constraint node adds a condition that makes the geometric entity fixed (fully
constrained); that is, the displacements are zero in all directions. For domains, this
condition is found under More on the context menu.
D O M A I N , B O U N D A R Y, E D G E , O R P O I N T S E L E C T I O N
From the Selection list, choose, the geometric entity (domains, boundaries, edges, or
points) that are fixed.
PAIR SELECTION
If Fixed Constraint is selected from the Pairs menu, choose the pair to define. An
identity pair has to be created first. Ctrl-click to deselect.
See Also Using Weak Constraints in the COMSOL Multiphysics Users Guide
Prescribed Displacement
The Prescribed Displacement feature adds a condition where the displacements are
prescribed in one or more directions to the geometric entity (domain, boundary, edge,
or point).
If a displacement is prescribed in one direction, this leaves the solid free to deform in
the other directions. Also define more general displacements as a linear combination
of the displacements in each direction.
D O M A I N , B O U N D A R Y, E D G E , O R PO I N T S E L E C T I O N
From the Selection list, choose the geometric entity (domains, boundaries, edges, or
points) to prescribe a displacement.
PAIR SELECTION
If Prescribed Displacement is selected from the Pairs menu, choose the pair to define.
An identity pair has to be created first. Ctrl-click to deselect.
Standard Notation
To define the displacements individually, click the Standard notation button (the
default).
General Notation
To specify the displacements using a General notation that includes any linear
combination of displacement components, click the General notation Hu=R button. For
example, for 2D models, use the relationship
H u = R
v
Enter values in the H matrix and R vector fields. For the H matrix, also select an
Isotropic, Diagonal, Symmetric, or Anisotropic matrix and enter values as required. For
example, to achieve the condition u = v, use the settings
H = 1 1 R = 0
0 0 0
See Also Using Weak Constraints in the COMSOL Multiphysics Users Guide
Symmetry
The Symmetry feature adds a boundary condition that represents symmetry in the
geometry and in the loads. A symmetry condition is free in the plane and fixed in the
out-of-plane direction.
BOUNDARY SELECTION
From the Selection list, choose the boundaries that are symmetry boundaries.
PAIR SELECTION
If Symmetry is selected from the Pairs menu, choose the pair to define. An identity pair
has to be created first. Ctrl-click to deselect.
CONSTRAINT SETTINGS
To display this section, select click the Show button ( ) and select Advanced Physics
Options. Select a Constraint typeBidirectional, symmetric or Unidirectional. If required,
select the Use weak constraints check box.
See Also Using Weak Constraints in the COMSOL Multiphysics Users Guide
BOUNDARY SELECTION
From the Selection list, choose the boundaries that are antisymmetry boundaries.
PAIR SELECTION
If Antisymmetry is selected from the Pairs menu, choose the pair to define. An identity
pair has to be created first. Ctrl-click to deselect.
CONSTRAINT SETTINGS
To display this section, select click the Show button ( ) and select Advanced Physics
Options. Select a Constraint typeBidirectional, symmetric or Unidirectional. If required,
select the Use weak constraints check box.
See Also Using Weak Constraints in the COMSOL Multiphysics Users Guide
Roller
The Roller node adds a roller constraint as the boundary condition; that is, the
displacement is zero in the direction perpendicular (normal) to the boundary, but the
boundary is free to move in the tangential direction.
BOUNDARY SELECTION
From the Selection list, choose the boundaries that have roller constraints.
CONSTRAINT SETTINGS
To display this section, select click the Show button ( ) and select Advanced Physics
Options. Select a Constraint typeBidirectional, symmetric or Unidirectional. If required,
select the Use weak constraints check box.
See Also Using Weak Constraints in the COMSOL Multiphysics Users Guide
Periodic Condition
The Periodic Condition feature adds a periodic boundary condition. This periodicity
makes uix0uix1 for a displacement ui. Control the direction that the periodic
condition applies to. Right-click the Periodic Condition node to add a Destination
Selection boundary condition. If the source and destination boundaries are rotated
with respect to each other, this transformation is automatically performed, so that
corresponding displacement components are connected.
BOUNDARY SELECTION
From the Selection list, choose the boundaries to define a periodic boundary condition.
The software automatically identifies the boundaries as either source boundaries or
destination boundaries.
This works fine for cases like opposing parallel boundaries. In other cases
use a Destination Selection subnode to control the destination. By default
it contains the selection that COMSOL Multiphysics has identified.
If Floquet Periodicity is selected, enter a k-vector for Floquet periodicity kF (SI unit:
rad/m) for the X, Y, and Z coordinates (3D models), or the R and Z coordinates
(2D axisymmetric models), or X and Y coordinates (2D models).
If Cyclic symmetry is selected, select a Sector angleAutomatic (the default), or User
defined. If User defined is selected, enter a value for S (SI unit: rad). For any
selection, also enter a Mode number m (unitless).
For information about this feature, see About Infinite Element Domains
and Perfectly Matched Layers in the COMSOL Multiphysics Users
Note Guide.
DOMAIN SELECTION
From the Selection list, choose the domains to define a viscoelastic solid and compute
the displacements, stresses, and strains.
MODEL INPUT
Define model inputs, for example, the temperature field if the material model uses a
temperature-dependent material property. If no model inputs are required, this section
is empty.
L O N G - TE R M E L A S T I C P R O P E R T I E S
Define the long-term elastic properties of the viscoelastic material.
From the Specify list, select the applicable long-term elastic property pair for an
isotropic materialYoungs modulus and Poissons ratio, Bulk modulus and shear
modulus, Lam constants, or Pressure-wave and shear-wave speeds. For each set of
properties, select From material or enter a User defined value or expression in each of
the following fields as required:
Youngs modulus (elastic modulus) E (SI unit: Pa) and Poissons ratio
(dimensionless).
Bulk modulus K (SI unit: Pa) and Shear modulus G (SI unit: Pa).
Density
By default, the material Density (SI unit: kg/m3) uses values From material. Select
User defined to enter a different value or expression.
GEOMETRIC NONLINEARITY
Thermal Effects
Right-click the Linear Viscoelastic Material node to add a Thermal Effects feature. This
feature defines these thermal effects:
Temperature.
Temperature Effects
See Also
DOMAIN SELECTION
From the Selection list, choose the domains to define the thermal effects.
By default, this node inherits the selection from its parent node, and only
a selection that is a subset of the parent nodes selection can be used.
Note
MODEL INPUTS
From the Temperature T (SI unit: K) list, select an existing temperature variable from
a heat transfer interface (for example, Temperature (ht/sol1)), if any temperature
variables exist, or select User defined to enter a value or expression for the temperature
(the default is 293.15 K).
THERMAL EFFECTS
Specify the thermal properties that define the thermal effects.
Thermal Expansion
Select the Thermal expansion check box to include thermal expansion for the selected
domains and to define the Coefficient of thermal expansion and the Strain reference
temperature Tref to determine the thermal expansion together with the actual
temperature.
WLF reference temperature (or glass transition temperature) Twlf (SI unit: K)
WLF constant 1, and WLF constant 2 (material-dependent constants) C1wlf and C2wlf.
If the Thermal Stress interface is used, the thermal effects are included in
the Thermal Linear Viscoelastic Material node.
Tip
Rigid Connector
The Rigid Connector is a boundary condition for modeling rigid regions and kinematic
constraints such as prescribed rigid rotations. If the study step is geometrically
nonlinear, the rigid connector takes finite rotations into account. By coupling selective
degrees of freedom of two rigid connectors it is also possible to create various types of
mechanisms such as hinges, joints, and other mechanical systems.
Right-click to add Harmonic Perturbation, Applied Force, Applied Moment, Mass and
Moment of Inertia, or Rigid Domain nodes to the rigid connector.
BOUNDARY SELECTION
From the Selection list, choose the boundaries to connect this rigid connector.
PAIR SELECTION
If Rigid Connector is selected from the Pairs menu, choose the pair to define. An identity
pair has to be created first. Ctrl-click to deselect.
CENTER OF ROTATION
Select a Center of rotationAutomatic (the default) or User defined. If Automatic is
selected, the center of rotation is at the geometrical center of the selected boundaries.
The constraints are applied at the center of rotation. Any mass is also considered to be
located there.
If User defined is selected, enter x, y, and (in 3D) z coordinates (SI unit: m) for the
center of rotation XC in the Global coordinates of center of rotation table.
If Constrained rotation is selected, select one or more of the Constrain rotation about X,
Constrain rotation about Y, and (in 3D) Constrain rotation about Z axis check boxes in
If Prescribed rotation at center of rotation is selected, enter an Axis of rotation (in 3D)
and an Angle of rotation (SI unit: rad). In 3D, the axis of rotation is given in the
selected coordinate system. In 2D, the axis of rotation is always the z-axis.
Harmonic Perturbation
Applied Force
Right-click the Rigid Connector node to add the Applied Force feature, which adds a
force to the rigid connector. The force can act at an arbitrary position in space.
LOCATION
The Center of rotation is selected by default. Select User defined to explicitly specify the
point Xp where the force is applied. The coordinates are always given in the global
coordinate system.
APPLIED FORCE
Enter values or expressions for the components of the Applied force F (SI unit: N).
Applied Moment
Right-click the Rigid Connector node to add the Applied Moment subnode, which adds
a moment at the center of rotation.
APPLIED MOMENT
Enter values or expressions for the components of the Applied moment M in 3D or the
applied moment in the z direction Mz in 2D (SI unit: Nm).
CENTER OF MASS
The Center of rotation is selected by default. Select User defined to explicitly specify the
point Xm (SI unit: m) where the mass is located. The coordinates are always given in
the global coordinate system.
Rigid Domain
Right-click the Rigid Connector node to add the Rigid Domain subnode, which calculates
the mass and moment of inertia properties of one or more domains. The purpose is to
be able to use such domains as rigid bodies.
RIGID DOMAIN
From material is selected as the default Density (SI unit: kg/m3). In this case the
material assignment for the domain supplies the mass density. Select User defined to
enter another value or expression.
Contact
The Contact node defines boundaries where the parts can come into contact but
cannot penetrate each other under deformation. Use it for modeling structural contact
and multiphysics contact.
PAIR SELECTION
When Contact is selected from the Pairs menu, choose the pair to define. An identity
pair has to be created first. Ctrl-click to deselect.
NORMAL CONTACT
Edit or use the default Contact normal penalty factor pn (SI unit: N/m3). The default
value is min(1e-3*5^segiter, 1)*solid.Eequiv/solid.hmin_dst. The penalty
factor controls how hard the interface surface is. A very low penalty factor would
allow the boundaries to penetrate each other somewhat, while a very large value may
give convergence difficulties. The default value causes the penalty factor to be
increased during the iterations and takes material stiffness and element size at the
contact surface into account. Eequiv is an equivalent Youngs modulus for the material
on the destination, and hmin_dst is the minimum element size on the destination.
Enter a value or expression for Contact surface offset from geometric destination surface
offset (SI unit: m). The default is 0. The offset is added to the in the normal direction
of the destination surface.
Enter a value or expression for Contact surface offset from geometric source surface
offset (SI unit: m). The default is 0. The offset is added to the in the normal direction
of the source surface.
Use this offset property to adjust initial clearances (negative values) or interference fits
(positive values) without having to change the geometry. The property is also useful
for studying the effects of geometrical tolerance when the structure is still modeled
using its nominal size.
INITIAL VALUES
Enter an initial value for the Contact pressure Tn (SI unit: N/m2). The default is 0.
Friction
The Friction feature is a subnode to the Contact node, which is added to model contact
with friction, for example to model the friction as static Coulomb friction or
exponential dynamic Coulomb friction.
BOUNDARY SELECTION
This section can not be edited, but can be used to highlight the surfaces that are part
of the pairs selected.
FRICTION
Select a Friction model. This can be either Static Coulomb friction or Exponential dynamic
Coulomb friction.
Edit or use the default Contact tangent penalty factor pt (SI unit: N/m3). The default
value is min(1e-3*5^segiter, 1)*solid.Eequ/3/solid.hmin_dst. The penalty
factor controls how hard the surfaces are bonded when sticking to each other. A very
low penalty factor allows the surfaces to slide somewhat relative to each other, while a
very large value might give convergence difficulties. The default value causes the
penalty factor to be increased during the iterations and takes material stiffness and
element size at the contact surface into account. Eequ is an equivalent Youngs
modulus for the material on the destination, and hmin_dst is the minimum element
size on the destination.
Edit the Static frictional coefficient stat (dimensionless) to give the coefficient of
friction.
Edit or use the default Cohesion sliding resistance Tcohe (SI unit: N/m2). Supply a
traction which must be overcome before sliding can occur.
Edit or use the default Maximum tangential traction Tt,max (SI unit: N/m2). The
default value is Inf, indicating that no limit on the tangential traction is active.
In the case of Exponential dynamic Coulomb friction, two more parameters must be
given. Enter the Dynamic friction coefficient dyn (dimensionless) and the Friction decay
coefficient dcf (SI unit: s/m).
To display this section, select click the Show button ( ) and select
Advanced Physics Options. Initial data is not often required for the friction
Note variables.
Enter values or expressions for the components of the initial force acting on the
destination surface as Friction force Tt (SI unit: N/m2).
For the Previous contact state, select either Not in contact or In contact. This determines
whether friction effects are active already when starting the solution or not. When
using In contact, also give values for Previous mapped source coordinates (SI unit: m).
The Previous mapped source coordinates serve as initial values for computing the
tangential slip. The default value is (X, Y, Z) and indicates that the contacting
boundaries are perfectly coincident in the initial state. The mapped source coordinates
are defined as the location on the source boundary where it is hit by a certain point on
the destination boundary.
ADVANCED
For the Contact tolerance, select either Automatic or Manual. The contact tolerance is
the gap distance between the surfaces at which friction starts to act. The default value
is 106 times smaller than the overall size of the model. If Manual is selected, enter an
absolute value for the Contact tolerance contact (SI unit: m).
Thermal Expansion
Thermal Expansion is an internal thermal strain caused by changes in temperature
according to the following equation for the thermal strain:
th = T T ref
where is the coefficient of thermal expansion (CTE), T is the temperature, and Tref
is the strain-free reference temperature.
When you have the Nonlinear Structural Materials Module, you can also
Note
right-click Hyperelastic Material to add the Thermal Expansion node.
DOMAIN SELECTION
From the Selection list, choose the domains to define the coefficient of thermal
expansion and the different temperatures that cause thermal stress.
By default, this node inherits the selection from its parent node, and only
a selection that is a subset of the parent nodes selection can be used.
Note
MODEL INPUTS
From the Temperature T (SI unit: K) list, select an existing temperature variable from
a heat transfer interface (for example, Temperature (ht/sol1)), if any temperature
variables exist, or select User defined to enter a value or expression for the temperature
(the default is 293.15 K).
THERMAL EXPANSION
For the Linear Elastic Material, select Isotropic, Diagonal or Symmetric to enter one or
more components for a general coefficient of thermal expansion vector vec.
The Shell interface Elastic Material Model has slightly different Initial Stress
Note
and Strain settings and this feature is discussed in that section.
When you have the Nonlinear Structural Materials Module, this model is
a good example.
DOMAIN SELECTION
From the Selection list, choose the domains to define the initial stress or strain.
By default, this node inherits the selection from its parent node, and only
a selection that is a subset of the parent nodes selection can be used.
Note
For a 3D Initial stress model, diagonal components 0x, 0y, and 0z and
off-diagonal components 0xy, 0yz, and 0xz, for example.
For a 3D Initial strain model, diagonal components 0x, 0y, and 0z and off-diagonal
components 0xy, 0yz, and 0xz, for example.
Phase
Add a Phase feature to a Body, Boundary, Edge, or Point Load. For modeling the
frequency response the physics interface splits the harmonic load into two parameters:
The amplitude, F, which is specified in the feature node for the load.
The phase (FPh)
Together these define a harmonic load, for which the amplitude and phase shift can
vary with the excitation frequency, f
F freq = F f cos 2f + F Ph f
D O M A I N , B O U N D A R Y, E D G E , O R P O I N T S E L E C T I O N
From the Selection list, choose the geometric entity (domains, boundaries, edges, or
points) to define the phase.
By default, this node inherits the selection from its parent node, and only
a selection that is a subset of the parent nodes selection can be used.
Note
Typically the load magnitude is a real scalar value. If the load specified in
the parent feature contains a phase (using a complex-valued expression),
the software adds the phase from the Phase node to the phase already
Note
included in the load.
Prescribed Velocity
The Prescribed Velocity feature adds a boundary or domain condition where the
velocity is prescribed in one or more directions. The prescribed velocity condition is
applicable for time-dependent and frequency-domain studies. With this boundary or
domain condition it is possible to prescribe a velocity in one direction, leaving the solid
free in the other directions.
P R E S C R I B E D VE L O C I T Y
Prescribed Acceleration
The Prescribed Acceleration feature adds a boundary or domain condition, where the
acceleration is prescribed in one or more directions. The prescribed acceleration
condition is applicable for time-dependent and frequency-domain studies. With this
boundary condition, it is possible to prescribe a acceleration in one direction, leaving
the solid free in the other directions.
PRESCRIBED ACCELERATION
Spring Foundation
The Spring Foundation feature has elastic and damping boundary conditions for
domains, boundaries, edges, and points. To select this feature for the domains, it is
selected from the More submenu. Also right-click to add a Pre-Deformation subnode.
The Spring Foundation and Thin Elastic Layer features are similar, with the difference
that a Spring Foundation connects the structural part on which it is acting to a fixed
ground, while a Thin Elastic Layer acts between two parts, either on an interior
boundary or on a pair boundary.
D O M A I N , B O U N D A R Y, E D G E , O R PO I N T S E L E C T I O N
From the Selection list, choose the geometric entity (domains, boundaries, edges, or
points) for the spring foundation.
SPRING
Select the spring type and its associated spring constant of force. The first option is the
default spring type for the type of geometric entity and space dimension.
then enter values or expressions into the table for each coordinate.
2D Axi
Force per volume as function of extension FV (SI unit: N/m3), or
Total force as function of extension Ftot (SI unit: N)
then enter values or expressions into the table for each coordinate.
Edges
Points
2D
Enter values or expressions in the table for each coordinate based on space
dimension for the Loss factor for spring k. The loss factors act on the
2D Axi corresponding components of the spring stiffness.
3D
VISCOUS DAMPING
Domains
Boundaries
Points
2D
Enter a value or expression for the Total viscous force dtot (SI unit: Ns/
m).
2D Axi
3D
Pre-Deformation
Right-click the Spring Foundation or Thin Elastic Layer nodes to add a
Pre-Deformation feature as a subnode and define the coordinates for the Spring
Pre-Deformation. By including a pre-deformation, you can model cases where the
unstressed state of the spring is in another configuration than the one modeled.
D O M A I N , B O U N D A R Y, E D G E , O R PO I N T S E L E C T I O N
From the Selection list, choose the geometric entity (domains, boundaries, edges, or
points) for the spring foundation.
SPRING PRE-DEFORMATION
Based on space dimension, enter the coordinates for the Spring Pre-Deformation u0
(SI unit: m)
The Thin Elastic Layer and Spring Foundation features are similar, with the difference
that a Spring Foundation connects the structural part on which it is acting to a fixed
ground.
BOUNDARY SELECTION
From the Selection list, choose the boundaries for the thin elastic layer.
PAIR SELECTION
If Thin Elastic Layer is selected from the Pairs menu, choose the pair to define. An
identity pair has to be created first. Ctrl-click to deselect.
A default Free node is added when a Thin Elastic Layer pair node is added.
Note
2D
Enter values or expressions in the table for each coordinate based on space
dimension for the Loss factor for spring k. The loss factors act on the
2D Axi corresponding components of the spring stiffness.
3D
Added Mass
The Added Mass feature is available on domains, boundaries, and edges and can be used
for supplying inertia which is not part of the material itself. Such inertia does not need
to be isotropic, in the sense that the inertial effects are not the same in all directions.
To select this feature for the domains, it is selected from the More submenu for the Solid
Mechanics interface.
D O M A I N , B O U N D A R Y, O R E D G E S E L E C T I O N
From the Selection list, choose the geometric entity (domains, boundaries, or edges)
for the added mass.
M A S S TY P E
Domains
Select a Mass typeMass per unit volume V (SI unit: kg/m3), Mass typeMass per unit
area A (SI unit: kg/m2), or Total mass m (SI unit: kg). Then enter values or
expressions into the table for each coordinate based on space dimension.
Boundaries
Select a Mass typeMass per unit area A (SI unit: kg/m2) or Total mass m (SI unit:
kg). Then enter values or expressions into the table for each coordinate based on space
dimension.
Edges
Select a Mass typeMass per unit length L (SI unit: kg/m) or Total mass m (SI unit:
kg). Then enter values or expressions into the table for each coordinate based on space
dimension.
Low-Reflecting Boundary
Use the Low-Reflecting Boundary feature to let waves pass out from the model without
reflection in time-dependent analysis. As a default, it takes material data from the
domain in an attempt to create a perfect impedance match for both pressure waves and
shear waves. It may be sensitive to the direction of the incoming wave.
BOUNDARY SELECTION
From the Selection list, choose the boundaries for the low-reflecting boundary.
DAMPING
Select a Damping typeP and S waves (the default) or User defined.
If User defined is selected, enter values or expressions for the Mechanical impedance di
(SI unit: Pas/m).
Continuum mechanics theory also makes use of a second set of coordinates, known as
material (or reference) coordinates. These are normally denoted by uppercase
variables X, Y, and Z (or R, , and Z) and are used to label material particles. Any
material particle is uniquely identified by its position in some given initial or reference
configuration. As long as the solid stays in this configuration, material and spatial
coordinates of every particle coincide and displacements are zero by definition.
When the solid objects deform due to external or internal forces and constraints, each
material particle keeps its material coordinates X (bold font is used to denote
coordinate vectors), while its spatial coordinates change with time and applied forces
such that it follows a path
x = x X t = X + u X t (3-1)
in space. Because the material coordinates are constant, the current spatial position is
uniquely determined by the displacement vector u, pointing from the reference
position to the current position. The global Cartesian components of this displacement
vector in the spatial frame, by default called u, v, and w, are the primary dependent
variables in the Solid Mechanics interface.
By default, the Solid Mechanics interface uses the calculated displacement and
Equation 3-1 to define the difference between spatial coordinates x and material
coordinates X. This means the material coordinates relate to the original geometry,
while the spatial coordinates are solution dependent.
Loads and constraints, on the other hand, are applied in spatial directions, by default
in the canonical spatial coordinate system. This system has basis vectors in the x, y, and
z directions, which are forever fixed in space. Both the material and spatial default
coordinate system are referred to as the Global coordinate system in the user interface
and are denoted by (G) in the following theory section.
Vector and tensor quantities defined in the global coordinate system on either frame
use the frames coordinate variable names as indices in the tensor component variable
names. For example, SXY is the material frame XY-plane shear stress, also known as a
second Piola-Kirchhoff stress, while sxy is the corresponding spatial frame stress, or
Cauchy stress. There are also a few mixed tensors, most notably the deformation
gradient FdxY, which has one spatial and one material index because it is used in
converting quantities between the material and spatial frames.
It is possible to define any number of user coordinate systems on the material and
spatial frames. Most types of coordinate systems are specified only as a rotation of the
basis with respect to the canonical basis in an underlying frame. This means that they
can be used both in contexts requiring a material system and in contexts requiring a
spatial one.
The coordinate system can be selected separately for each added material model, load,
and constraint. This is convenient if, for example, an anisotropic material with different
orientation in different domains is required. The currently selected coordinate system
is known as the local coordinate system and is denoted by (L) throughout this text.
Lagrangian Formulation
The formulation used for structural analysis in COMSOL Multiphysics for both small
and finite deformations is total Lagrangian. This means that the computed stress and
deformation state is always referred to the material configuration, rather than to
current position in space.
The gradient of the displacement, which occurs frequently in the following theory, is
always computed with respect to material coordinates. In 3D:
u u u
X Y Z
u = v v v
X Y Z
w w w
X Y Z
1 T
= --- u + u
2
or in components as
1 u m u n
mn = --- + (3-2)
2 x n x m
The Duhamel-Hookes law relates the stress tensor to the strain tensor and
temperature:
s = s 0 + C 0
where C is the 4th order elasticity tensor, : stands for the double-dot tensor product
(or double contraction), s0 and 0 are initial stresses and strains, TTref, and is
the thermal expansion tensor.
1
W s = --- 0 C 0 (3-3)
2
1 ijmn
0 0
Ws = --- C ij ij ij mn mn mn
2
i j m n
TE N S O R V S . M A T R I X F O R M U L A T I O N S
Because of the symmetry, the strain tensor can be written as the following matrix:
x xy xz
xy y yz
xz yz z
Similar representation applies to the stress and the thermal expansion tensors:
s x s xy s xz x xy xz
s xy s y s yz xy y yz
s xz s yz s z xz yz z
D 15 D 25 D 35 D 45 D 55 D 56 C
1123
C
2223
C
3323
C
1223
C
2323
C
2313
0 0
0 0
0 0
Different pairs of elastic moduli can be used, and as long as two moduli are defined,
the others can be computed according to Table 3-2.
TABLE 3-2: EXPRESSIONS FOR THE ELASTIC MODULI.
Bulk modulus K E - 2
----------------------- + -------
3 1 2 3
Shear modulus G E -
--------------------
21 +
Lam constant E 2G
--------------------------------------
- K --------
1 + 1 2 3
Lam constant E - G
--------------------
21 +
Pressure-wave cp K + 4G 3-
speed -------------------------
Shear-wave cs G
speed
+ 2 0 0 0
+ 2 0 0 0
D = + 2 0 0 0
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
4G 2G 2G
K + -------- K -------- K -------- 0 0 0
3 3 3
2G 4G 2G
K -------- K + -------- K -------- 0 0 0
3 3 3
D = 2G 2G 4G
K -------- K -------- K + -------- 0 0 0
3 3 3
0 0 0 G 0 0
0 0 0 0 G 0
0 0 0 0 0 G
11 1 x
22 2 y
33 3 z
12 21 4 xy
23 32 5 yz
13 31 6 xz
thus, the Hookes law is presented in the form involving the elasticity matrix D and the
following vectors:
Beside the Standard (XX, YY, ZZ, XY, YZ, XZ) Material data ordering, the elasticity
coefficients are entered following the Voigt notation. In the Voigt (XX, YY, ZZ, YZ, XZ,
XY) Material data ordering, the sorting of indices is:
11 1 x
22 2 y
33 3 z
23 32 4 yz
13 31 5 xz
12 21 6 xy
thus the last three rows and columns in the elasticity matrix D are swapped.
ORTHOTROPIC MATERIAL
The elasticity matrix for orthotropic material in the Standard (XX, YY, ZZ, XY, YZ, XZ)
Material data ordering has the following structure:
D 11 D 12 D 13 0 0 0
D 12 D 22 D 23 0 0 0
D 13 D 23 D 33 0 0 0
D =
0 0 0 D 44 0 0
0 0 0 0 D 55 0
0 0 0 0 0 D 66
2 2
E x E y E y xy yz + xz E y E z xz E x
D 13 = ---------------------------------------------------------- , D 22 = ----------------------------------------
D denom D denom
2
E y E z E y xy xz + E x yz E y E z E y xy E x
D 23 = ----------------------------------------------------------------- , D 33 = -----------------------------------------------
D denom D denom
D 44 = G xy , D 55 = G yz , and D 66 = G xz
where
2 2 2 2
D denom = E y E z xz E x E y + 2 xy yz xz E y E z + E x E z yz + E y xy
The values of Ex, Ey, Ez, xy, yz, xz, Gxy, Gyz, and Gxz are supplied in designated
fields in the user interface. COMSOL deduces the remaining componentsyx, zx,
and zyusing the fact that the matrices D and D1 are symmetric. The compliance
matrix has the following form:
1 yx zx
------ -------- -------- 0 0 0
Ex Ey Ez
xy 1 zy
-------- ------ -------- 0 0 0
Ex Ey Ez
xz yz 1
1
-------- -------- ------ 0 0 0
D = Ex Ey Ez
1
0 0 0 --------- 0 0
G xy
1
0 0 0 0 --------- 0
G yz
1
0 0 0 0 0 ---------
G xz
The elasticity matrix in the Voigt (XX, YY, ZZ, YZ, XZ, XY) Material data ordering changes
the sorting of the last three elements in the elasticity matrix:
D 44 = G yz , D 55 = G xz , and D 66 = G xy
ANISOTROPIC MATERIAL
In the general case of fully anisotropic material, provide explicitly 21 components of
the symmetric elasticity matrix D, in either Standard (XX, YY, ZZ, XY, YZ, XZ) or Voigt
(XX, YY, ZZ, YZ, XZ, XY) Material data ordering, and 6 components of the symmetric
thermal expansion matrix.
F = f 0 T + W s T
where WsT is given by Equation 3-3. Hence, the stress can be found as
F W
s = = = C 0
T T
F
= C p log T T 0 + S elast
T
S elast = s
S elast = s x + s y + s z
Q h =
where is the strain-rate tensor and the tensor represents all possible inelastic stresses
(for example, a viscous stress).
Using the tensor components, the heat balance can be rewritten as:
T
C p
t
+ Tmn t smn = k T + Q h (3-4)
m n
In many cases, the second term can be neglected in the left-hand side of Equation 3-4
because all Tmn are small. The resulting approximation is often called uncoupled
thermoelasticity.
Strain-Displacement Relationship
The strain conditions at a point are completely defined by the deformation
componentsu, v, and w in 3Dand their derivatives. The precise relation between
strain and deformation depends on the relative magnitude of the displacement.
SMALL DISPLACEMENTS
Under the assumption of small displacements, the normal strain components and the
shear strain components are related to the deformation as follows:
u xy 1 u v
x = xy = ------- = --- +
x 2 2 y x
v yz 1 v w
y = yz = ------- = --- + (3-5)
y 2 2 z y
w xz 1 u w
z = xz = ------- = --- + .
z 2 2 z x
To express the shear strain, use either the tensor form, xy, yz, xz, or the engineering
form, xy, yz, xz.
The symmetric strain tensor consists of both normal and shear strain components:
The strain-displacement relationships for the axial symmetry case for small
displacements are
u u w u w
r = , = --- , z = , and rz = +
r r z z r
LARGE DEFORMATIONS
As a start, consider a certain physical particle, initially located at the coordinate X.
During deformation, this particle follows a path
x = x X t
For simplicity, assume that undeformed and deformed positions are measured in the
same coordinate system. Using the displacement u, it is then possible to write
x = X+u
x
dx = ------- dX = F dX
X
is used.
The deformation gradient contains the complete information about the local straining
and rotation of the material. It is a positive definite matrix, as long as material cannot
be annihilated. The ratio between current and original volume (or mass density) is
dV- 0
---------- = ------ = det F = J
dV 0
T
C = F F
1 1 T
= --- C I = --- F F I
2 2
1 u i u j u k u k
ij = --- -------- + -------- + --------- --------- (3-6)
2 X j X i X i X j
The deformation state characterized by finite (or large displacements) but small to
moderate strains is sometimes referred to as geometric nonlinearity or nonlinear
geometry. This typically occurs when the main part of the deformations presents a finite
rigid body rotation
v
L kl = r t
xl k
where v k r t is the spatial velocity field. It can be shown that L can be computed in
terms of the deformation gradient as
dF 1
L = F
dt
The velocity gradient can be decomposed into symmetric and skew-symmetric parts
L = Ld + Lw
where
1 T
L d = --- L + L
2
and
1 T
L w = --- L L
2
is called the spin tensor. Both tensors are defined on the spatial frame.
It can be shown that the material time derivative of the Green-Lagrange strain tensor
can be related to the rate of strain tensor as
d T
= F Ld F
dt
The spin tensor Lw(x,t) accounts for an instantaneous local rigid-body rotation about
an axis passing through the point x.
Components of both Ld and Lw are available as results and analysis variables under the
Solid Mechanics interface.
Stress-Strain Relationship
The symmetric stress tensor describes stress in a material:
x xy xz
= yx y yz xy = yx xz = zx yz = zy
zx zy z
This tensor consists of three normal stresses (x, y, z) and six (or, if symmetry is used,
three) shear stresses (xy, yz, xz).
For large deformations and hyperelastic material models there are more than one stress
measure:
1
S = F P
1 T 1 T
= J PF = J FSF
In case of plane strain, set i3 for i1, 2, 3. The out-of-plane stress components
si3 are results and analysis variables.
Axial Symmetry
The axially symmetric geometry uses a cylindrical coordinate system. Such a coordinate
system is orthogonal but curvilinear, and one can choose between a covariant basis
e1, e2, e3 and a contravariant basis e1, e2, e3.
1 0 0
g
ij = 0 r2 0
0 0 1
The metric tensor plays the role of a unit tensor for a curvilinear coordinate system.
For any vector or tensor A, the metric tensor can be used for conversion between
covariant, contravariant, and mixed components:
Aim g
j mj
Ai =
m
ni mj
Anm g
ij
A = g
m n
In both covariant and contravariant basis, the base vector in the azimuthal direction
has a nonunit length. To cope with this issue, the so called physical basis vectors of unit
length are introduced. These are
1 1 3
e r = e 1 = e e = --- e = re 2 e z = e 3 = e
r 2
The corresponding components for any vector or tensor are called physical.
phys ij
A ij = g ii g jj A
1
Ai
i i
I 1 A = trace A i = = A 11 + A 22 ----2- + A 33
i
r
1 j
Aj Ai
2 i
I 2 A = --- I 1 A
2
i j
i
I 3 A = det A i
STRAINS
The right Cauchy-Green deformation tensor is defined as
T T
C = u + u + u u + g
1
= --- C g
2
u u 2 w 2
C 11 = 2 + + + 1
r r r
C 12 = C 23 = 0
u w u -----
u- ------
w w
C 13 = + + + - -------
z r r z r z
2 2
C 22 = 2ru + u + r
w u 2 w 2
C 33 = 2 ------- + ------ + ------- + 1
z z z
For the Linear Elastic Material, drop the nonlinear terms inside square brackets in the
above expressions.
1
r = 11 = --- C 11 1
2
r = z = 0
1
rz = 13 = --- C 13
2
22 1- 2
= -------
2
- = --------
2
C 22 r
r 2r
1
z = 33 = --- C 33 1
2
vol = I 1 = r + + z
STRESSES
The second Piola-Kirchhoff stress tensor is computed as
W s
S = 2
C
ij W s
S = 2
C ij
The corresponding physical components are defined on the global coordinate system:
12
S r = rS
13
S rz = S
2 22
S = r S
23
S z = rS
33
Sz = S
Note that
S
ij
I1 S = g ij = S r + S + S z
i j
S
ij
Stest = test ij
i j
For the linear elastic material, the stress components in coordinate system are
ij ij ijkl
s = s0 + C kl kl 0kl
where TTref.
For anisotropic and orthotropic materials, the 4th-order elasticity tensor is defined
from D matrix according to:
The user input D matrix always contains the physical components of the elasticity
tensor
phys
C ijkl
and the corresponding tensor components are computed internally according to:
phys
ijkl C ijkl
C = -----------------------------------------------
g ii g jj g kk g ll
ijkl ij kl ik jl il jk
C = g g + g g + g g
where and are the first and second Lam elastic constants.
Loads
Specify loads as
In the same way, choose between defining the load as force per unit
volume or force per unit area for body loads acting in a domain. Also
2D
define a total force (SI unit: N) as required.
Table 3-3 shows how to define distributed loads on different geometric entity levels;
the entries show the SI unit in parentheses.
TABLE 3-3: DISTRIBUTED LOADS
For general cases, if the problem is linear in all other respects, the solution
can be made more efficient by forcing the solver to treat the problem as
linear. See Stationary Solver in the COMSOL Multiphysics Reference
Note
Guide.
THEORY
The direction of an explicitly applied distributed load must be given with reference to
a local or global coordinate system in the spatial frame, but its magnitude must be with
reference to the undeformed reference or material area. That is, the relation between
the true force f acting on the current area da and the specified distributed load F
acting on the material area dA is f = FdA. When the solid is subjected to an external
pressure, p, the true force on a surface element acts with magnitude p in the current
area da in the normal direction n:
f = pnda
da
F = pn --------
dA
Plane Stress
In a plane stress condition the out-of-plane deformation causes the thickness to
change, and this area effect is included explicitly. The equation transforms to
w
F = pn -------- 1 +
da
dA z
Axial Symmetry
To account for the radial deformation changing the circumference and therefore the
area element, the distributed load is applied as
da R + u
F = pn -------- -------------------
dA R
Equation Implementation
The COMSOL Multiphysics implementation of the equations in the Solid Mechanics
interface is based on the principle of virtual work.
The principle of virtual work states that the sum of virtual work from internal strains
is equal to work from external loads.
The total stored energy, W, for a linear material from external and internal strains and
loads equals:
W = s + u FV dv
V
U
t
+ u F S ds + u F L dl + Fp
S L p
Utest Fp
t
+ u test F S ds + u test F L dl +
S L p
STATIONARY STUDIES
COMSOL Multiphysics uses an implementation based on the stress and strain
variables. The normal and shear strain variables depend on the displacement
derivatives.
Using the tensor strain, stress, and displacement variables, the principle of virtual work
is expressed as:
W = test s + utest FV dv
V
Utest Fp
t
+ u test F S ds + u test F L dl +
S L p
V I S C O E L A S T I C TR A N S I E N T I N I T I A L I Z A T I O N
The equations using the linear viscoelastic material model corresponds to a stationary
study, with added viscoelasticity.
t
+ u test F S ds + u test F L dl + U test Fp
S L p
where the terms proportional to dM and dK appear if the Rayleigh damping is used.
For more information about the equation form in case of geometric nonlinearity see
FREQUENCY-DOMAIN STUDIES
In the frequency domain the frequency response is studied when applying harmonic
loads. Harmonic loads are specified using two components:
To derive the equations for the linear response from harmonic excitation loads
F xfreq = F x f cos t + F xPh f ----------
180
F xfreq
F freq = F yfreq
F zfreq
assume a harmonic response with the same angular frequency as the excitation load
u = u amp cos t + u
u
u= v
w
j u jt jt j
u = Re u amp e e = Re u e where u = u amp e u
e jt
u = Re u
jF xPh f ----------
180 jt jt
F xfreq = Re F x e e = Re Fx e
where
jF xPh f ----------
180
Fx = F x f e
Fx
F = F
y
Fz
EIGENFREQUENCY STUDIES
The eigenfrequency equations are derived by assuming a harmonic displacement field,
similar as for the frequency response formulation. The difference is that this study type
uses a new variable j explicitly expressed in the eigenvalue jThe
eigenfrequency f is then derived from j as
f = Im j -
------------------
2
Im
Q = -------------------
2Re
= Re
RAYLEIGH DAMPING
To model damping effects within the material, COMSOL Multiphysics uses Rayleigh
damping, where two damping coefficients are specified.
The weak contribution due to the alpha-damping is always accounted for as shown in
Equation 3-7. The contribution from the beta-damping that shown in Equation 3-7
corresponds to the case of small strains. In case of geometric nonlinearity, it becomes
dM utest Pt dv
V
where P is the first Piola-Kirchhoff stress tensor, for more information see
To further clarify the use of the Rayleigh damping, consider a system with a single
degree of freedom. The equation of motion for such a system with viscous damping is
2
d u du
m ---------2- + c ------- + ku = f t
dt dt
In the Rayleigh damping model the damping coefficient c can be expressed in terms
of the mass m and the stiffness k as
c = dM m + dK k
The Rayleigh damping proportional to mass and stiffness is added to the static weak
term.
A complication with the Rayleigh damping model is to obtain good values for the
damping parameters. A more physical damping measure is the relative damping, the
ratio between actual and critical damping, often expressed as a percentage of the critical
damping. Commonly used values of relative damping can be found in the literature.
1 dM
= --- ----------- + dK 2f
2 2f
Using this relationship at two frequencies, f1 and f2, with different relative damping,
1 and 2, results in an equation system that can be solved for dM and dK:
1 -
----------- f
4f 1 1 dM 1
=
1 - dK 2
----------- f
4f 2 2
Using the same relative damping, 1 = 2, does not result in a constant damping factor
inside the interval f1 f f2. It can be shown that the damping factor is lower inside
the interval, as Figure 3-7 shows.
Relative damping
Rayleigh damping
Specified damping
f1 f2
f
s = s 0 + 1 + j s C 0
W s
S = 1 + j s ----------
E
If modeling the damping in the structural analysis via the loss factor, use the following
definition for the elastic part of the entropy:
S elast = s j s C
This is because the entropy is a function of state and thus independent of the strain
rate, while the damping represents the rate-dependent effects in the material (for
example, viscous or viscoelastic effects). The internal work of such inelastic forces
averaged over the time period 2 can be computed as:
1
Q h = --- s Real Conj C (3-8)
2
Equation 3-8 can be used as a heat source for modeling of the heat generation in
vibrating structures, when coupled with the frequency-domain analysis for the stresses
and strains.
The Linear Viscoelastic Material feature typically does not require any
additional damping, as all the damping effects are contained in the
Note material data (such as the relaxation times).
x = x X t
For simplicity, assume that undeformed and deformed positions are measured in the
same coordinate system.Using the displacement u, it is then possible to write
x = X + u X t
x
dx = ------- dX = F dX
X
The deformation gradient F contains the complete information about the local
straining and rotation of the material. It is a two-point tensor (or a double vector),
which transforms as a vector with respect to each of its indexes. It involves both the
reference and present configurations. In terms of displacement gradient, it can be
written as F u I. The ratio between current and initial volume (or mass density) is
dV- 0
---------- = ------ = det F = J
dV 0
T T T
C = F F = u + u + u u + I (3-9)
u i u j u k u k
C ij = -------- + -------- +
X j X i --------
- --------- + ij
X i X j
(3-10)
k
1
= --- C I (3-11)
2
The first one is to include the geometric nonlinearity within the Linear Elastic
Material.
The second option is to use the Hyperelastic Material feature. This approach is
suitable for the case of large strains, but the choice of a particular hyperelastic
material is usually relevant only for a particular type of materials (such as rubbers,
for example).
s = pI + s d
p = K vol 3 T T ref
where K is the bulk modulus, is the coefficient of thermal expansion, and the strain
is decomposed as
1
= --- vol I + d
3
vol = trace ij
t
d
s d = 2 t t' -------- dt'
t'
0
where the function (t) is called the relaxation shear modulus function that can be
found by measuring the stress evolution in time when the material is held at a constant
strain.
G m exp -------
t
t = G + m
m=1
Physical interpretation of this approach, which is often called the generalized Maxwell
model, is shown in Figure 3-8
s
G1 G2 Gi
G ..... .....
1 2 i
1 2 i
s
Hence, m are the relaxation time constants of the spring-dashpot pairs in the same
branch, and Gm represent the stiffness of the spring in branch m.
Following the analogy, the abstract variables qm are introduced to represent the
extension of the corresponding springs. Thus,
s d = 2G d + 2G m q m (3-12)
m=1
Each of the variables qm is a symmetric tensor (which has as many components as the
number of strain components of the problem class).
for each branch. This assumes that traceq0m0. This is equivalent to the
assumption that any possible elastic prestrain of the material (also called viscoelastic
transient initialization) is always fast enough to equally affect all the branches. Hence,
trace q m 0
which implies that both viscoelastic strain and stress tensors are deviatoric, and it allows
to express:
q zm = q xm q ym
G0 = G + Gm
m=1
A commonly used parameter, the elastoviscosity, can be expressed for each branch in
terms of the relaxation time and shear modus as
m = Gm m
TE M P E R A T U RE E F F E C T S
For many materials, the viscoelastic properties have a strong dependence on the
temperature. A common assumption is that the material is thermorheologically simple
(TRS). In a material of this class, a change in the temperature can be transformed
directly into a change in the time scale. The reduced time is defined as
t
dt'
tr = -----------------------
a T T t'
-
0
C1 T T0
log a T = -----------------------------------
C 2 + T T0
Note that aT(T0) = 1, so that T0 is the temperature at which the original material data
is given. Usually T0 is taken as the glass transition temperature of the material. If the
temperature drops below T0C2, the WLF equation is no longer valid. The constants
C1 and C2 are material dependent, but with T0 as the glass temperature the values
C1 = 17.4 and C2 = 51.6 K are reasonable approximations for many polymers.
Think of the shift function as a multiplier to the viscosity in the dashpot in the Maxwell
model. Equation 3-13 for a TRS material is modified to
1
q m + ----------------------- q m = d
a T T m
The viscoelastic strain variables qm are treated as additional degrees of freedom. The
shape functions must be of one order less than those used for the displacements
because these variables add to the strains and stresses computed from displacement
derivatives. The viscoelastic strain variables do not require continuity, so it is possible
to use discontinuous shape functions.
s d = real s d exp jt
d = real d exp jt
s d = 2 G' + jG'' d
where
N
m
G'' = G m ----------------------------2
1 + m
m=1
are often referred to as the storage modulus and the loss modulus, respectively.
The internal work of viscous forces averaged over the time period 2 can be
computed as:
Q h = G''d conj d
AXIAL SYMMETRY
The volumetric strain is
g
ij
vol = I 1 = ij
i j
1
d = --- vol g
3
1
dij = ij --- vol g ij
3
g
ij
q mij = 0
i j
q 22 m
q 33 m = q 11 m + --------------
2
-
r
and all the other components of qm are treated as extra degrees of freedom.
N
2G m q ij m g g
ni kj
ij
sq =
k n m = 1
11
s qr = s q
12
s q r = rsq
13
s q rz = s q
2 22
sq = r sq
23
s q z = rsq
33
sqz = sq
ij ij ijkl ij
s = s0 + C kl kl T T ref 0kl + s q
ijkl ij kl ik jl il jk
C = g g + g g + g g
Both the contact map operator map (E) and the gap distance variable are defined by
the contact element elcontact. For each destination point where the operator or gap
is evaluated, a corresponding source point is sought by searching in the direction
normal to the destination boundary.
Source
m(x)
g
x
Destination
Before the boundaries come in contact, the source point found is not necessarily the
point on the source boundary closest to the destination point. However, as the
boundaries approach one another, the source point converges to the closest point as
the gap distance goes to zero.
It is possible to add an offset value to either the source, the destination or both. If an
offset is used, the geometry is treated as larger (or smaller) than what is actually
modeled when the gap is computed.
Using the special gap distance variable (solid.gap), the penalized contact pressure
Tnp is defined on the destination boundary as
where dg is the gap distance between the destination and source boundary, and pn is
the user-defined normal penalty factor.
The penalized friction traction Ttp is defined on the destination boundary as:
T tcrit
T tp = min -------------------- 1 T ttrial (3-15)
T ttrial
and
x m = map x
In Equation 3-16, pt is the user-defined friction traction penalty factor, and xm,old is
the value of xm in the last time step, and
map F x m x m old
is the vector of slip since the last time step (approximated using a backward Euler step).
Ttcrit is defined as
T np g + T tp m F x m dA + w cn T n + w ct T t dA
destination destination
where wcn and wct are contact help variables defined as:
w ct = fric T tp n T tp n i T t i + 1
where i is the augmented solver iteration number and fric is a Boolean variable stating
if the parts are in contact.
Since the contact model always makes use of the geometric nonlinearity assumption,
the contact pressure and the friction force represent, respectively, the normal and
tangential components of the nominal traction that is the force with components in
the actual configuration (that is, on the spatial frame) but related to the undeformed
area of the corresponding surface element.
FRICTION
The friction model is either no friction or Coulomb friction.
+ s d exp dcf v s if dynamic friction
= d
s otherwise
The only degrees of freedom needed to represent this assembly are the ones needed to
represent the movement of a rigid body. In 2D this is simply two in-plane translations,
and the rotation around the z-axis.
2 2 2 2
a +b +c +d = 1
is added.
2 2 2 2
a +b c d 2bc 2 ad 2ac + 2bd
R = 2ad + 2bc
2 2
a b +c d
2 2
2cd 2ab
2 2 2 2
2bd 2 ac 2ab + 2cd a b c +d
Under pure rotation, a vector from the center of rotation (Xc) of the rigid connector
to a point X on the undeformed solid will be rotated into
x Xc = R X Xc
u = x X = R I X Xc
When the center of rotation of the rigid connector also has a translation uc, then the
complete expression for the displacements on the solid is
u = R I X Xc + uc
The total rotation of the rigid connector can be also presented as a rotation vector. Its
definition is
The parameter a can be considered as measuring the rotation, while b, c, and d can be
interpreted as the orientation of the rotation vector. For small rotations, this relation
simplifies to
b
= 2 c
d
The rotation vector is available as the variables thx_tag, thy_tag, and thz_tag.
Both the initial stress and strains are tensor variables defined via components on the
local coordinate system for each domain. Input these as the following matrices:
1
p 0 = --- I 1 s 0
3
AXIAL SYMMETRY
User inputs the physical components of 0 and s0:
0r 0r 0rz s 0r s 0r s 0rz
0r 0 0z s 0r s 0 s 0z
0rz 0z 0z s 0rz s 0z s 0z
SPRING FOUNDATION
A spring gives a force that depends on the displacement and acts in the opposite
direction (in the case of a force that is proportional to the displacement, this is called
Hookes law). In a suitable coordinate system, a spring condition can be represented as
fs = f u u0
In the same way, a viscous damping can be described as a force proportional to the
velocity
f v = Du
Structural (loss factor) damping is only relevant for frequency domain analysis and
is defined as
f l = iK u u 0
where is the loss factor and i is the imaginary unit (in this case, a constant or a
diagonal matrix). If the elastic part of the spring definition is given as a force versus
displacement relation, the stiffness K is taken as the stiffness at the linearization point
at which the frequency response analysis is performed. Since the loss factor force is
proportional to the elastic force, the equation can be written as
f sl = f s + f l = 1 + i f s
T
W = u f sl + f v dA
A
T H I N E L A S T I C L A Y E R B E T W E E N TW O P A R T S
A spring or damper can also act between two boundaries of an identity pair. The spring
force then depends on the difference in displacement between the surfaces.
f sD = f sS = K u D u S u 0
The uppercase indices refer to source and destination. When a force versus
displacement description is used,
f sD = f sS = f u u 0
u = uD uS
f vD = f vS = D u D u S
f lD = f lS = iK u D u S u 0
so that
f slD = f sD + f lD = 1 + i f sD
The virtual work expression is formulated on the destination side of the pair as
uD uS
T
W = f slD + f vD dA D
AD
Here the displacements from the source side are obtained using the src2dst operator
of the identity pair.
f sd = f su = K u d u u u 0
or
f sd = f su = f u u 0
u = ud uu
f vd = f vu = D u d u u
f ld = f lu = iK u d u u u 0
f sld = f sd + f ld = 1 + i f sd
The subscripts u and d denote the up and down sides of the interior boundary.
ud us
T
W = f sld + f vd dA D
AD
Other uses for added mass are when sheets or strips of a material that is heavy, but
having a comparatively low stiffness, are added to a structure. The data for the base
material can then be kept unaltered, while the added material is represented purely as
added mass.
The value of an added mass can also be negative. You can use such a negative value for
adjusting the mass when a part imported from a CAD system does not get exactly the
correct total mass due so simplifications of the geometry.
Added mass is an extra mass distribution that can be anisotropic. It can exist on
domains, boundaries, and edges. The inertial forces from added mass can be written as
2
u
fm = M 2
t
where M is a diagonal mass distribution matrix. The contribution to the virtual work is
T
W = u f m dA
A
for added mass on a boundary, and similarly for objects of other dimensions.
For other materials, it can be activated via the solver setting. Note however that even
together with the geometric nonlinearity, the validity of any linear material model is
usually limited to the situation of possibly large displacements but small to moderate
strains. A typical example of use is to model large rigid body rotations. The
implementation is similar to that for the geometrically linear elastic material, but with
the strain tensor replaced with the Green-Lagrange strain tensor, and the stress tensor
replaced with the second Piola-Kirchhoff stress tensor, defined as:
0 0 i
S = S + C I
where TTref, and i represents all possible inelastic strains (such as plastic or
creep strains). In components, it is written as:
0 0 i
S ij = S ij + C ijkl kl kl kl kl
where the elasticity tensor C ijkl is defined from the D matrix (user input). The 2nd
Piola-Kirchhoff stress is a symmetric tensor.
1 0 i
W s = --- S I
2
which is a variable defined in the physics interface. Other stress variables are defined as
follows.
The first Piola-Kirchhoff stress P is calculated from the second Piola-Kirchhoff stress
as P FS. The first Piola-Kirchhoff stress relates forces in the present configuration
with areas in the reference configuration, and it is sometimes called the nominal stress.
Using the 1st Piola-Kirchhoff stress tensor, the equation of motion can be written in
the following form:
where the density corresponds to the material density in the initial undeformed state,
the volume force vector FV has components in the actual configuration but given with
respect to the undeformed volume, and the tensor divergence operator is computed
with respect to the coordinates on the material frame. Equation 3-18 is the strong
form that corresponds to the weak form equations solved in case of geometric
nonlinearity within the Solid Mechanics interface (and many related multiphysics
interfaces) in COMSOL. Using vector and tensor components, the equation can be
written as
2
ux P xX P xY P xZ
0 = F Vx + +
2 X Y Z
t
2
uy P yX P yY P yZ
0 = F Vy + +
t
2 X Y Z
2
uz P zX P zY P zZ
0 = F Vz + +
t
2 X Y Z
The components of 1st Piola-Kirchhoff stress tensor are non symmetric in the general
case, thus
P iJ P Ij
because the component indexes correspond to different frames Such tensors are called
two-point tensors.
The boundary load vector FA in case of geometric nonlinearity can be related to the
1st Piola-Kirchhoff stress tensor via the following formula:
FA = P n0
where the normal n0 corresponds to the undeformed surface element. Such force
vector is often referred to as the nominal traction. In components, it can be written as
1 T 1 T
s = J PF = J FSF
The Cauchy stress is a true stress that relates forces in the present configuration (spatial
frame) to areas in the present configuration, and it is a symmetric tensor.
Equation 3-18 can be rewritten in terms of the Cauchy stress as
2
u
2
= fV x s
t
where the density corresponds to the density in the actual deformed state, the volume
force vector fV has components in the actual configuration (spatial frame) given with
respect to the deformed volume, and the divergence operator is computed with respect
to the spatial coordinates.
1
p = --- I 1 s
3
which corresponds to the volumetric part of the Cauchy stress. The deviatoric part is
defined as
s d = s + pI
1
J 2 s = --- s d :s d
2
s mises = 3J 2 s
In such cases a mixed formulation can be used that represents the pressure as a
dependent variable in addition to the displacement components. This formulation
removes the effect of volumetric strain from the original stress tensor and replaces it
with an interpolated pressure, pw. A separate equation constrains the interpolated
pressure to make it equal (in a finite-element sense) to the original pressure calculated
from the strains.
For an isotropic linear elastic material, the second Piola-Kirchhoff stress tensor S,
computed directly from the strains, is replaced by a modified version:
s = s + pI p w I
where I is the unit tensor and the pressure p is calculated from the strains and material
data as
1
p = --- trace s
3
pw p
------- --- = 0 (3-19)
The modified stress tensor s is used then in calculations of the energy variation.
For orthotropic and anisotropic materials, the auxiliary pressure equation is scaled to
make the stiffness matrix symmetric. Note, however, that the stiffness matrix in this
formulation is not positive definite and even contains a zero block on the diagonal in
As a default, the low-reflecting boundary condition takes the material data from the
adjacent domain in an attempt to create a perfect impedance match for both pressure
waves and shear waves, so that
u u
n = c p n n c s t t
t t
where n and t are the unit normal and tangential vectors at the boundary, respectively,
and cp and cs are the speeds of the pressure and shear waves in the material. This
approach works best when the wave direction in close to the normal at the wall.
u
n = d i c p c s
t
where the mechanical impedance di is a diagonal matrix available as the user input, and
cp + cs
d i = ---------------- I
2
The problem can be modeled using the full solution without applying the above
described multiplicative decomposition. For such a solution, the Floquet periodicity
conditions at the corresponding boundaries of the periodicity cell are expressed as
The cyclic symmetry boundary condition presents a special but important case of
Floquet periodicity, for which the unit periodicity cell is a sector of a structure that
possesses rotational symmetry. The frequency response problem can be solved then in
one sector of periodicity by applying the periodicity condition. The situation is often
referred to as dynamic cyclic symmetry.
For an eigenfrequency study, all the eigenmodes of the full problem can be found by
performing the analysis on one sector of symmetry only and imposing the cyclic
symmetry of the eigenmodes with an angle of periodicity = m , where the cyclic
The Floquet periodicity conditions at the sides of the sector of symmetry can be
expressed as
i T
u destination = e R u source
where the u represents the displacement vectors with the components given in the
default cartesian coordinates. Multiplication by the rotation matrix given by
cos sin 0
R = sin cos 0
0 0 1
The angle represents either the periodicity of the eigenmode for an eigenfrequency
analysis or the periodicity of the excitation signal in case of a frequency-response
analysis. In the latter case, the excitation is typically given as a load vector
F = F 0 exp im atan Y X
when modeled using the cartesian coordinates; parameter m is often referred to as the
azimuthal wave-number.
This chapter describes the Shell and Plate interfaces, which are found under the
Structural Mechanics branch ( ) in the Model Wizard.
In this chapter:
213
T he S he ll a nd Pl at e In t erfaces
The Shell ( ) and Plate ( ) interfaces, found under the Structural Mechanics
branch ( ) in the Model Wizard, have the equations, elements, and functionality (in
the form of material properties, constraints, loads, thermal expansion, and initial
stresses and strains) for stress analysis and general structural mechanics in shells and
plates. Both interfaces use shell elements of the MITC type.
Model
Pinched Hemispherical Shell: Model Library path
Structural_Mechanics_Module/Verification_Models/
pinched_hemispherical_shell
The Linear Elastic Material is the only available material model. It adds a linear elastic
equation for the displacements and has a settings window to define the elastic material
properties.
When this interface is added, these default nodes are also added to the Model Builder
Linear Elastic Material, Free (a boundary condition where edges are free, with no loads
or constraints), and Initial Values. Right-click the Shell or Plate node to add other
nodes.
INTERFACE IDENTIFIER
The interface identifier is a text string that can be used to reference the respective
physics interface if appropriate. Such situations could occur when coupling this
The default identifier (for the first interface in the model) is shell or plate.
THICKNESS
Define the Thickness d by entering a value or expression (SI unit: m) in the field. The
default is 0.01 m. Use the Change Thickness feature to define a different thickness in
parts of the shell or plate. The thickness can be variable if an expression is used.
Offset Definition
If the actual shell midsurface is not the boundary on which the mesh exists, it is
possible to prescribe an offset in the direction of the surface normal by using an offset
definition. The offset is defined as positive if the shell midsurface is displaced from the
meshed boundary in the direction of the positive shell normal.
Select an option from the Offset definition listNo offset (the default), Physical offset,
or Relative offset. The default No offset means that the modeled boundary coincides
with the shell midsurface.
If Physical offset is selected, enter a value or expression in the zoffset field (SI unit:
m) for the actual distance from the meshed boundary to the shell midsurface. The
default is 0.
If Relative offset is selected, enter a value or expression in the zrel_offset field (
unitless) for the offset as the ratio between the offset distance and half the shell
thickness. The default is 0. A value of +1 means that the actual shell bottom surface
This section is available for the Shell interface only. Also see The MITC
Shell Formulation.
Note
The fold-line limit angle (SI unit: radians) is the smallest angle for a fold line; that is,
if the angle where two or more shell elements meet is larger than this value, it becomes
a fold line. The normal to the shell is discontinuous along a fold-line. Enter a value or
expression in the field. The default value is 0.05 radians (approximately 3 degrees).
ADVANCED SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Normally these settings do not need to be changed.
The Use MITC interpolation check box is selected by default, and this interpolation,
which is part of the MITC shell formulation, should normally always be active.
For the Plate interface, the Use 3D formulation check box is used to select whether six
or three variables are used in the formulation. For geometrically nonlinear analyses, or
when in-plane (membrane) forces are active, six variables must be used. This check box
is selected by default.
DISCRETIZATION
Show button ( ) and select Discretization. Select Quadratic (the default) or Linear for
the Displacement field. Specify the Value type when using splitting of complex variables
Real or Complex (the default).
Antisymmetry
Body Load
Change Thickness
Damping
Edge Load
Face Load
Initial Values
Initial Stress and Strain
Linear Elastic Material
No Rotation
Phase
Pinned
Point Load
Prescribed Acceleration (for time-dependent and frequency-domain studies)
Prescribed Displacement/Rotation
Prescribed Velocity (for time-dependent and frequency-domain studies)
Rigid Connector
Added Mass
Applied Force
Applied Moment
Free
Mass and Moment of Inertia
Fixed Constraint
Pre-Deformation
Spring Foundation
If Orthotropic is selected, due to the shell assumptions, no values for Ez, yz, or xz
need to be entered.
If a User defined Anisotropic model is selected, a 6-by-6 symmetric matrix is
displayed. Due to the shell assumptions, only enter values for D11, D12, D22, D14,
D24, D55, D66, and D56.
GEOMETRIC NONLINEARITY
In this section there is always one check box. Either Force linear strains or Include
geometric nonlinearity is shown.
If a study step is geometrically nonlinear, the default behavior is to use a large strain
formulation in all domains. There are however some cases when you would still want
to use a small strain formulation for a certain domain. In those cases, select the Force
linear strains check box. When selected, a small strain formulation is always used,
independently of the setting in the study step. The default value is that the check box
is cleared, except when opening a model created in a version prior to 4.3. In this case
the state is chosen so that the properties of the model are conserved.
The Include geometric nonlinearity check box is displayed only if the model was created
in a version prior to 4.3, and geometric nonlinearity was originally used for the selected
domains. It is then selected and forces the Include geometric nonlinearity check box in
the study step to be selected. If the check box is cleared, the check box is permanently
Thermal Expansion
Right-click the Linear Elastic Material node to add a Thermal Expansion node, which is
an internal thermal strain caused by changes in temperature. The temperature is
assumed to vary linearly through the thickness of the shell.
By default, this node inherits the selection from its parent node, and only
a selection that is a subset of the parent nodes selection can be used.
Note
THERMAL EXPANSION
Specify the thermal properties that define the thermal stress.
From the Coefficient of thermal expansion (SI unit: 1/K) list, select From material to
use the coefficient of thermal expansion from the material, or User defined to enter a
value or expression for . Select Isotropic, Diagonal or Symmetric to enter one or more
components for a general coefficient of thermal expansion vector vec.
Enter a value or expression of the Strain reference temperature Tref (SI unit: K), which
is the reference temperature where the thermal strain is zero.
From the Temperature T (SI unit: K) list, select an existing temperature variable from
a heat transfer interface (for example, Temperature (ht/sol1)), if any temperature
variables exist, or select User defined to enter a value or expression for the temperature
(the default is 293.15 K). This is the mid-surface temperature of the shell.
THERMAL BENDING
Enter the Temperature difference in thickness direction T, (SI unit: K), which affects
the thermal bending. This is the difference between the temperatures at the top and
bottoms surfaces.
For more information see Stress and Strain Calculations and Initial Values
and Prescribed Values in the theory section.
See Also
By default, this node inherits the selection from its parent node, and only
a selection that is a subset of the parent nodes selection can be used.
Note
INITIAL STRESS
Specify the initial stress as the Initial in-plane force, the Initial moment, and the Initial
out-of-plane force. The default values are zero, which mean no initial stress. Enter
values or expressions in the applicable fields for the:
For definitions of the generalized strains, see Theory for the Shell and
Plate Interfaces.
See Also
Damping
Right-click the Linear Elastic Material node to add the Damping node, which adds
Rayleigh damping or an isotropic loss-factor damping. In time-dependent,
eigenfrequency, and frequency domain response studies, model undamped or damped
problems.
By default, this node inherits the selection from its parent node, and only
a selection that is a subset of the parent nodes selection can be used.
Note
DAMPING SETTINGS
Select a Damping typeRayleigh damping (the default), Isotropic loss factor, Orthotropic
loss factor, or Anisotropic loss factor.
From the Loss factor for orthotropic shear modulus list, select From material (the
default) to use the material value or select User defined and enter values for G or GVo
If Voigt notation is selected, the loss factors for the shear moduli are given in that order
too.
The default values are 0. If the values are taken from the material, these loss factors are
found in the orthotropic material model for the material.
The defaults are 0. If the values are taken from the material, these loss factors are found
in the anisotropic material model for the material.
Change Thickness
Use the Change Thickness feature to model boundaries with a thickness or offset other
than the overall values defined in the interface Thickness section.
Offset Definition
If the actual shell midsurface is not the boundary on which the mesh exists, it is
possible to prescribe an offset in the direction of the surface normal by using an offset
definition. The offset is defined as positive if the shell midsurface is displaced from the
meshed boundary in the direction of the positive shell normal.
Select an option from the Offset definition listFrom parent (the default), No offset,
Physical offset, or Relative offset. The default From parent means that the offset is as
defined in the Thickness section of the shell interface.
If No offset is selected, it means that the modeled boundary coincides with the shell
midsurface.
If Physical offset is selected, enter a value or expression in the zoffset (SI unit: m) field
for the actual distance from the meshed boundary to the shell midsurface. The
default is 0.
If Relative offset is selected, enter a value or expression in the zrel_offset (unitless)
field for the offset as the ratio between the offset distance and half the shell
thickness. The default is 0. A value of +1 means that the actual shell bottom surface
is located on the meshed boundary, and a value of -1 means that the shell top surface
is located on the meshed boundary.
Initial Values
The Initial Values feature adds the initial values for the translational displacement and
velocity field as well as the for the normal displacement and velocity field. It can serve
as an initial condition for a transient simulation or as an initial guess for a nonlinear
analysis. Right-click to add additional Initial Values nodes.
INITIAL VALUES
Enter values or expressions for the initial values. The default value is 0 for all
displacements and velocities:
About the Body Load, Face Load, Edge Load, and Point Load Features
Add force loads acting on all levels of the geometry. Add a:
Body Load (such as self-weight loads) for the shells on the boundary (plate:
domain) level.
Face Load to boundaries (plate: domains), for example describing a pressure acting
on a surface.
Edge Load as a force or moment distributed along an edge (plate: boundary).
Point Load as concentrated forces or moments at points.
For all of these loads, right-click and choose Phase to add a phase for harmonic loads
in frequency-domain computations. In this way a harmonic load can be defined where
the amplitude and phase shift can vary with the excitation frequency f:
Body Load
Add a Body Load to boundaries (for the Plate interface add it to domains) and use it for
self-weight loads, for example. Right-click and add a Phase for harmonic loads in
frequency-domain computations.
FORCE
Enter values or expressions for the components of the Body load FV (SI unit: N/m3).
MOMENT
Enter values or expressions for the components of the Moment body load ML (SI unit:
Nm/m3).
Face Load
Add a Face Load to boundaries (for the Plate interface add it to domains), to use it as a
pressure or tangential force acting on a surface. Right-click and add a Phase for
harmonic loads in frequency-domain computations.
FORCE
Select a Load typeFace load defined as force per unit area or Pressure.
If Face load defined as force per unit area is selected, enter values or expressions for
the components of the Face load FA (SI unit: N/m2).
MOMENT
Enter values or expressions for the components of the Moment face load MA (SI unit:
Nm/m2).
Edge Load
Add an Edge Load as a force or moment distributed along an edge (for the Plate
interface add it to boundaries). Right-click and add a Phase for harmonic loads in
frequency-domain computations.
MOMENT
From the Load type list, select the definition of the moment load. Select Edge load
defined as moment per unit length ML (the default) (SI unit: N) or Edge load defined as
moment per unit area MA (SI unit: Nm/m2). Enter values or expressions for the
components of the Moment edge load in the table.
Point Load
Add a Point Load to points for concentrated forces or moments at points. Right-click
and add a Phase for harmonic loads in frequency-domain computations.
PO IN T S EL EC TIO N
From the Selection list, choose the points to define a point load.
FORCE
Enter values or expressions for the components of the Point load FP (SI unit: N).
MOMENT
Enter values or expressions for the components of the Point moment MP (SI unit:
Nm).
Phase
Right-click any load node (see About the Body Load, Face Load, Edge Load, and
Point Load Features) to add a Phase node, which adds a phase for harmonic loads in
frequency-domain computations.
By default, this node inherits the selection from its parent node, and only
a selection that is a subset of the parent nodes selection can be used.
Note
PHASE
Add the phase angle Fph (SI unit: rad) for harmonic loads. Enter the phase for each
component of the load in the corresponding fields.
Pinned
The Pinned node adds an edge, boundary, domain, or point condition that fixes the
translations in all directions, that is, all displacements are zero.
E D G E , B O U N D A R Y, D O M A I N , O R PO I N T S E L E C T I O N
From the Selection list, choose the geometric entity (edges, boundaries, domains or
points) to prescribe a pinned condition.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Select a Constraint typeBidirectional, symmetric or Unidirectional. If required, select
the Use weak constraints check box.
E D G E , B O U N D A R Y, D O M A I N , O R P O I N T S E L E C T I O N
From the Selection list, choose the geometric entity (edges, boundaries, domains or
points) to prescribe no rotations.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Select a Constraint typeBidirectional, symmetric or Unidirectional. If required, select
the Use weak constraints check box.
Prescribed Displacement/Rotation
The Prescribed Displacement/Rotation node adds an edge, boundary, domain, or point
condition to a model where the displacements and rotations are prescribed in one or
more directions.
E D G E , B O U N D A R Y, D O M A I N , O R P O I N T S E L E C T I O N
From the Selection list, choose the geometric entity (edges, boundaries, domains or
points) where the displacement or rotation is to be prescribed.
Points
From the Coordinate system list select from:
PRESCRIBED DISPLACEMENT
To define a prescribed displacement for each space direction (x, y, and z), select one or
all of the Prescribed in x direction, Prescribed in y direction, and Prescribed in z direction
check boxes. Then enter a value or expression for the prescribed displacements u0, v0,
or w0 (SI unit: m).
PRESCRIBED ROTATIONS
Select a prescribed rotation from the By listFree, Rotation, or Normal vector. Select:
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Select a Constraint typeBidirectional, symmetric or Unidirectional. If required, select
the Use weak constraints check box.
Prescribed Velocity
The Prescribed Velocity node adds an edge, boundary, domain, or point condition
where the translational or rotational velocity is prescribed in one or more directions.
The prescribed velocity condition is applicable for time-dependent and
frequency-domain studies. With this condition it is possible to prescribe a velocity in
one direction, leaving the shell free in the other directions.
E D G E , B O U N D A R Y, D O M A I N , O R P O I N T S E L E C T I O N
From the Selection list, choose the geometric entity (edges, boundaries, domains or
points) to prescribe a translational or rotational velocity.
Boundaries
From the Coordinate system list select from:
Edges
From the Coordinate system list select from:
Points
From the Coordinate system list select from:
P R E S C R I B E D VE L O C I T Y
To define a prescribed velocity for each space direction (x, y, and z), select one or more
of the Prescribed in x direction, Prescribed in y direction, and Prescribed in z direction
check boxes. Then enter a value or expression for the prescribed velocity components
vx, vy, and vz (SI unit: m/s).
P R E S C R I B E D A N G U L A R VE L O C I T Y
To define a prescribed angular velocity for each space direction (x, y, and z), select one
or all of the Prescribed around x direction, Prescribed around y direction, and Prescribed
Prescribed Acceleration
The Prescribed Acceleration node adds an edge, boundary, domain, or point condition
where the translational or rotational acceleration is prescribed in one or more
directions. The prescribed acceleration condition is applicable for time-dependent and
frequency-domain studies. With this condition it is possible to prescribe an
acceleration in one direction, leaving the shell free in the other directions.
E D G E , B O U N D A R Y, D O M A I N , O R P O I N T S E L E C T I O N
From the Selection list, choose the geometric entity (edges, boundaries, domains or
points) to prescribe a translational or rotational velocity.
Boundaries
From the Coordinate system list select from:
Edges
From the Coordinate system list select from:
PRESCRIBED ACCELERATION
To define a prescribed acceleration for each space direction (x, y, and z), select one or
more of the Prescribed in x direction, Prescribed in y direction, and Prescribed in
z direction check boxes. Then enter a value or expression for the prescribed
acceleration components ax, ay, and az (SI unit: m/s2).
Symmetry
The Symmetry node adds an edge or boundary condition that defines a symmetry edge
(Plate: boundary) or boundary (Plate: domain).
E D G E , B O U N D A R Y, O R D O M A I N S E L E C T I O N
From the Selection list, choose the geometric entity (edges, boundaries or domains) to
prescribe symmetry.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Select a Constraint typeBidirectional, symmetric or Unidirectional. If required, select
the Use weak constraints check box.
Antisymmetry
The Antisymmetry node adds an edge or boundary condition that defines an
antisymmetry edge (Plate: boundary) or boundary (Plate: domain).
D O M A I N , B O U N D A R Y, O R E D G E S E L E C T I O N
From the Selection list, choose the geometric entity (edges, boundaries or domains) to
prescribe antisymmetry.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Select a Constraint typeBidirectional, symmetric or Unidirectional. If required, select
the Use weak constraints check box.
Rigid Connector
The Rigid Connector is a boundary condition for modeling between shell edges. The
feature is similar to the rigid connectors in the Solid Mechanics interface, except that no
rigid domain feature is available. Rigid connectors from Shell and Solid Mechanics
interfaces can be attached to each other.
EDGE SELECTION
From the Selection list, choose the edges to connect this rigid connector.
CENTER OF ROTATION
Select a Center of rotationAutomatic (the default) or User defined. The center of
rotation is at the geometrical centre of the selected edges. The constraints are applied
at the center of rotation. Any mass is also considered to be located there.
If User defined is selected, enter x, y, and z coordinates (SI unit: m) for the Global
coordinates of center of rotation Xc (SI unit: m) in the table.
If Constrained rotation is selected, select one or more of the Constrain rotation about
X, Constrain rotation about Y, and Constrain rotation about Z axis check boxes in order
Results Evaluation
For visualization and results evaluation, predefined variables include all nonzero stress
and strain tensor components, principal stresses and strains, in-plane and out-of-plane
forces, moments, and von Mises and Tresca effective stresses. It is possible to evaluate
the stress and strain tensor components and effective stresses at an arbitrary distance
from the midsurface. The parameter zshell (shell.z) is found in, for example, the
Parameters table of the surface plots settings window. Its value can be changed
between 1 (downside) to 1 (upside). A value of 0 means the midsurface of the shell.
The default value is given in the Height of Evaluation in Shell, [-1,1] section of the Shell
interface.
Stresses and strains are available both in the global coordinate system, and in the shell
local system as described in Local Coordinate Systems.
For a shell to give accurate results it is important that the structure can
really be described as being thin-walled. When modeling using shells it is
Important also important to define the faces as the midplane of the real geometry.
Plates are similar to shells but act in a single plane and usually with only out-of-plane
loads. The plate and shell elements in COMSOL are based on the same formulation.
The Plate interface for 2D models is a specialization of the Shell interface. In the
following, the text fully describes the Shell interface, and the Plate interface is
mentioned only where there are nontrivial differences.
To describe a shell, you provide its thickness and the elastic material properties.
The element used for the shell interface is of Mindlin-Reissner type, which means that
transverse shear deformation is accounted for. It can thus also be used for rather thick
shells. It has an MITC formulation where MITC means mixed interpolation of
tensorial components. A general description of this element family is in Ref. 1.
For a standard plate analysis only three degrees of freedom are needed: the
out-of-plane displacement w and the displacements of the shell normals ax and ay, but
it is also possible to activate all six degrees of freedom, so that any type of analysis of a
shell initially positioned in the xy-plane can be performed using the Plate interface.
Using six degrees of freedom is the default, but three degrees of freedom can be
selected instead for efficiency.
When six degrees of freedoms are used in the Plate interface, there must
be enough constraints to suppress any in-plane rigid body modes.
Note
Also here, the rotations x, y (and possibly z) are used to a large extent.
The position of the deformed midsurface is ru, and the normal after deformation is
na. To keep the normal a unit vector requires that
n+a = 1 (4-1)
na = 0
The vectors r, u, n, and a are interpolated by the nth-order Lagrange basis functions.
The basic assumption is that the position of a point within the shell after deformation
has a linear dependence of the thickness coordinate, and thus is
1 2 3 1 2 1 2 3 1 2 1 2
x = r ( , ) + u ( , ) + n ( , ) + a ( , )
STRAINS
The in-plane Green-Lagrange strain in the local covariant components can then be
written as
1 3 r u 3
= --- r + u + n + a + + n + a
2
r 3 r 3 3 3 2
+ n + n = + +
The indices and range from 1 to 2. The transverse shear strains in local covariant
components are
3 = 3 =
--- r + u + 3 n + a n + a r + 3 n n = + 3
1
2
The constitutive relation for the shell elements is a plane stress assumption, as is
customary in shell theory. The strain component in the normal direction 33 is thus
irrelevant. The different parts of the strain tensors above can be written out as
1 u r r u u u
= --- --------- -------- + --------- -------- + --------- --------
2
1 r a a r u n n u u a a u
= --- --------- -------- + --------- -------- + --------- -------- + --------- -------- + --------- -------- + --------- --------
2
1 a n n a a a
= --- --------- -------- + --------- -------- + --------- --------
2
1 n a a
= --- --------- a + --------- n + --------- a
2
OFFSET
It is possible to model a shell with a midsurface that is not located at the meshed
surface but at a certain offset from it. The offset is assumed to occur along the normal
of the shell surface. In this case,
1 2 1 2 1 2
r ( , ) = r R ( , ) + o n ( , )
where rR is the position of the meshed reference surface and o is the offset distance.
Since all geometric derivatives are computed on the mesh on the reference surface, the
following type of expressions are used when evaluating the strains:
r- r R n
--------
= --------- + o ---------
The degrees of freedom are referred to the midsurface even when there is an offset.
All loads and boundary conditions are assumed to be applied at the midsurface, so a
force acting in the plane of the shell will not cause any bending action if there is an
offset.
The numerical integration of the element is performed over the reference surface. If
the shell is curved, the area of the actual midsurface and the reference surface will
differ. This is compensated for by multiplying the weak expressions with an area scale
factor, defined as
Any expressions depending on the coordinates are evaluated on the shell reference
surface.
An offset applied to a flat shell does not have any effect, since the area
scale factor is 1, and the normal direction is constant.
Note
ROTATION REPRESENTATION
In a geometrically linear analysis, a rotation vector is defined as
= na
na
e = -----------------
na
= acos 1 + n a
This representation is unique only for rotations up to 180 degrees, but because the
rotation vector representation is only an output convenience, is has no impact on the
analysis.
Each contribution to the virtual work of the element is numerically integrated over the
midsurface while the integration in the thickness direction is performed analytically.
In regions where the discretized surface is smooth (which is always the case for plates),
the normal of the shell surface is uniquely defined. When two or more shell element
meet at an angle, each element must, however, keep its own normal direction. It is thus
not possible to have the same set of degrees of freedom for the normal deformations
in such a point. This is automatically handled by the program.The automatic search for
these fold lines compares the normals of all surfaces sharing an edge. If the angle
between them differs more than a certain angle (with a default of 3 degrees) it is
considered as a fold line. In this case the assumption is that the angle between the shell
normals is to remain constant. This gives
nj nk = nj + aj nk + ak
or
nj ak + aj nk + aj ak = 0 (4-2)
where the values or j and k range over the number of shells elements with different
normals. The third term in Equation 4-2 is relevant only in a large deformation analysis
because it is nonlinear. A special case occurs when two adjacent surfaces are parallel but
their normal vectors have opposite directions. In this case the shell element applies the
special constraint
ak = aj
I N I T I A L VA L U E S A N D P RE S C R I B E D VA L U E S
Because the normal vector displacements are quantities that are less intuitive than the
more customary nodal rotations, it is possible to specify the initial conditions and
prescribed values in terms of nodal rotations as well as in terms of the normal vector
displacement. The representation by normal vector direction is insensitive to whether
the analysis is geometrically nonlinear or not. The direction of the shell normal is
prescribed in the deformed state, n0. The prescribed values for the actual degrees of
freedom, a0, are internally computed as
n0
a 0 = --------- n
n0
a 0 = 0 n
where 0 is the vector of prescribed nodal rotations. This relation is fully defined only
when all three components of 0 are given. It is also possible to prescribe only one or
two of the components of 0, while leaving the remaining components free. Because
it has no relevance to prescribe the rotation about the normal direction of the shell, it
is only possible prescribe individual rotations in a shell local system. In this case, the
result becomes one or two constraint relations between the components of a0. The
directions are the edge local coordinate system where t1 is the tangent to the edge and
t2 is perpendicular and inward from the edge, in the plane of the shell. These
constraints are formulated as
t 2 a 0 = 01
t 1 a 0 = 02
a 0 = R 0 I n
where R is a standard rotation matrix, representing the finite rotation about the given
rotation vector.
a = n + a0
un = 0
a t2 = 0
u t1 = 0
u t2 = 0
for the antisymmetry case. Here t1 and t2 are two perpendicular directions in the plane
of the shell.
When applied to an edge, there is a local coordinate system where t1 is the tangent to
the edge, and t2 is perpendicular and in the plane of the shell. The assumption is then
that t2 is the normal to the symmetry or antisymmetry plane. The constraints are
u t2 = 0
a t2 = 0
un = 0
u t1 = 0
a t1 = 0
u e1 = 0
a e2 n = 0
a e3 n = 0
u e3 = 0
a e1 n = 0
for the antisymmetry case. Using a general coordinate system sometimes leads to
higher accuracy, since there is no element interpolation of the constraint directions
involved.
EXTERNAL LOADS
Contributions to the virtual work from the external load is of the form
u test F + a test M n + a
where the forces (F) and moments (M) can be distributed over a boundary or an edge
or concentrated in a point. The contribution from a is only included in a geometrically
nonlinear analysis. Loads are always referred to the midsurface of the element. In the
special case of a follower load, defined by its pressure p, the force intensity is
F = p n + a
For a follower load, the change in midsurface area is not taken into account, in order
to be consistent with the assumption that thickness changes are ignored.
Ni
m = D i T m T ref + ------
h
where D is the plane stress constitutive matrix, Ni are the initial membrane forces, and
i the initial membrane strains. The influence of thermal strains is included through
Dh T 6M i
b = -------- i -------- + ----------
-
2 h h
2
where i is the initial value of the bending part of the strain tensor, and Mi are the
initial bending and twisting moments. T is the temperature difference between the
top and bottom surface of the shell. The bending stress is the stress at the top surface
of the shell if no membrane stress is present.
5 Qi
s = --- 2G i + ------
6 h
where G represent the transverse shear moduli, i is the initial average shear strain, and
Qi are the initial transverse shear forces. The correction factor 56 ensures that the
stresses are averaged so that they correspond to the ratio between shear force and
thickness. The corresponding strains and i are averaged in an energy sense.
= m + z b
2
3 s 1 z
-----------------------------
2
The shell section forces (membrane forces, bending moments, and shear forces) are
computed from the stresses as
N = h m
2
h
M = ------
6 b
Q = h s
The local system for stress output coincides with the orientations defined for the
material input. Stresses are also available transformed to the global coordinate system.
1. The local z direction ezl is the positive normal of the shell surface.
2. The local x direction exl is the projection of the first direction in the material
coordinate system (ex1) on the shell surface
e x1 e x1 e zl e zl
e xl = ----------------------------------------------------
e x1 e x1 e zl e zl
This projection cannot be performed if ex1 is normal to the shell. In that case the
second axis ex2 of the material system will instead define exl using the same procedure.
Thus, if
e x1 e zl 0.99
then
e x2 e x2 e zl e zl
e xl = ----------------------------------------------------
e x2 e x2 e zl e zl
e yl = e zl e xl
When using an isotropic material, the only effect of selecting a local coordinate
system is that the definition directions of local stresses change.
When defining orthotropic and anisotropic materials, local coordinate systems do
not need to be created that exactly follow the shell surface. It is sufficient that the
local system when projected as described above gives the intended in-plane
directions.
Beams
This chapter describes the Beam interface, which you find under the Structural
Mechanics branch ( ) in the Model Wizard.
In this chapter:
253
The Beam Interface
3D
The Beam interface ( ), found under the Structural Mechanics branch ( ) in the
Model Wizard, has the equations, elements, and features for stress analysis and general
linear structural mechanics in beams and uses beam elements of the Euler (or
Euler-Bernoulli) type.
The Linear Elastic Material node is the default material model, which adds a linear elastic
equation for the displacements and has a settings window to define the elastic material
properties.
When this interface is added, these default nodes are also added to the Model Builder:
Linear Elastic Material, Cross Section Data, Free (a condition where points are free, with
no loads or constraints), and Initial Values.
Right-click the Beam node to add other features that implement, for example, loads
and constraints. The following sections provide information about all features available
in this interface.
INTERFACE IDENTIFIER
The interface identifier is a text string that can be used to reference the respective
physics interface if appropriate. Such situations could occur when coupling this
interface to another physics interface, or when trying to identify and use variables
defined by this physics interface, which is used to reach the fields and variables in
expressions, for example. It can be changed to any unique string in the Identifier field.
The displacement field u, which has two components (u, v) in 2D and three
components (u, v, and w) in 3D.
The plane rotation angle , which has one component in 2D (th) and three
components in 3D (thx, thy, and thz).
The names can be changed but the names of fields and dependent variables must be
unique within a model.
DISCRETIZATION
The beam elements are of a first-order type, so only first-order shape functions are
available in the Beam interface. Actually, the axial displacement and twist are
represented by linear shape functions, while the bending is represented by a cubic
shape function (Hermitian element).
See Also About the Edge Load and Point Load Features
Theory for the Beam Interface
These are described in this section (listed in alphabetical order). The list also includes
subfeatures:
Antisymmetry
Cross Section Data
Damping
Edge Load
Initial Stress and Strain
Initial Values
Linear Elastic Material
No Rotation
Prescribed Acceleration (for time-dependent and frequency-domain studies)
Prescribed Displacement/Rotation
Prescribed Velocity (for time-dependent and frequency-domain studies)
Pinned
Point Load
Point Mass
Point Mass Damping
Section Orientation
Symmetry
Thermal Expansion
Added Mass
Fixed Constraint
MODEL INPUTS
Define model inputs, for example, the temperature field of the material uses a
temperature-dependent material property. If no model inputs are required, this section
is empty.
Youngs modulus and Poissons ratio to specify Youngs modulus (elastic modulus) E
(SI unit: Pa) and Poissons ratio (unitless).
Bulk modulus and shear modulus to specify the bulk modulus K (SI unit: Pa) and the
shear modulus G (SI unit: Pa).
Lam constants to specify the Lam constants (SI unit: Pa) and (SI unit: Pa).
Pressure-wave and shear-wave speeds to specify the pressure-wave speed
(longitudinal wave speed) cp (SI unit: m/s) and the shear-wave speed (transverse
wave speed) cs (SI unit: m/s).
This is the wave speed for a solid continuum. In a truss or beam element,
the actual speed with which a longitudinal wave travels is lower than the
Note value given. When using this type of input, the density must also be given.
For each pair of properties, select from the applicable list to use the value From material
or enter a User defined value or expression.
Each of these pairs define the elastic properties, and it is possible to convert from one
set of properties to another.
Density
Define the Density (SI unit: kg/m3) of the material. Select From material to take the
value from the material or select User defined to enter a value for the density.
Thermal Expansion
Right-click the Linear Elastic Material node to add the Thermal Expansion node. Thermal
expansion is an internal thermal strain caused by changes in temperature. The
temperature is assumed to vary linearly across the beams cross section.
By default, this node inherits the selection from its parent node, and only
a selection that is a subset of the parent nodes selection can be used.
Note
THERMAL EXPANSION
Specify the thermal properties that define the thermal stress.
From the Coefficient of thermal expansion (SI unit: 1/K) list, select From material to
use the coefficient of thermal expansion from the material, or User defined to enter a
value or expression for .
Enter a value or expression of the Strain reference temperature Tref (SI unit: K), which
is the reference temperature where the thermal strain is zero.
From the Temperature T (SI unit: K) list, select an existing temperature variable from
a heat transfer interface (for example, Temperature (ht/sol1)), if any temperature
variables exist, or select User defined to enter a value or expression for the temperature
(the default is 293.15 K). This is the centerline temperature of the beam.
THERMAL BENDING
Enter the Temperature gradient in local y-direction Tgy (in 2D and 3D) and in the
Temperature gradient in local z-direction Tgz (in 3D) (SI unit: K/m), which affects the
thermal bending. If beam cross section dimensions have been defined at Bending stress
evaluation pointsFrom section heights, these could be used in an expression
containing the temperature difference.
For more information about the stress and strain quantities, see Stress and
Strain Calculations and Initial Values and Prescribed Values as described
See Also for the Shell and Plate interfaces.
By default, this node inherits the selection from its parent node, and only
a selection that is a subset of the parent nodes selection can be used.
Note
INITIAL STRESS
Specify the initial stress as the Initial axial force, the Initial bending moment, and for 3D,
the Initial torsional moment. The default values are zero, which means no initial stress.
For 2D models, enter values or expressions for the Initial axial force N (SI
unit: N) and Initial bending moment Miz (SI unit: Nm).
2D
For 2D models, enter values or expressions for Initial axial strain eni
(unitless) and Initial curvature siz (SI unit: rad/m).
2D
Damping
Right-click the Linear Elastic Material node to add the Damping node, which adds
Rayleigh damping or an isotropic loss-factor damping. In time-dependent,
eigenfrequency, and frequency domain response studies, undamped or damped
problems can be modeled.
By default, this node inherits the selection from its parent node, and only
a selection that is a subset of the parent nodes selection can be used.
Note
For Rayleigh damping, enter the Mass damping parameter dM and the Stiffness damping
parameter dK in the corresponding fields. The default values are 0, which is no
damping.
For Isotropic loss factor damping, from the Isotropic structural loss factor list, select
From material (the default) to use the material value or select User defined to enter a
value or expression. The default value is 0.
Initial Values
The Initial Values feature adds an initial value for the displacement field u (the
displacement components u, v, and, in 3D, w), the velocity field, the rotations, and
the angular velocity. It can serve as an initial condition for a transient simulation or as
an initial guess for a nonlinear analysis. Right-click to add additional Initial Values
features.
INITIAL VALUES
For 2D models, enter values or expressions for the following. The default
value is 0 for all initial values:
Rectangular Section
Box Section
Circular Section
Rectangle
Enter values or expressions for the following.
Box
Enter values or expressions for the following.
Circular
Enter a value or expression for the Diameter do (SI unit: m).
Pipe
Enter values or expressions for the following.
This section is only available if User defined is selected as the Cross Section
Definition.
Note
The following table lists the basic section properties (some apply in 3D only). Enter
values for these properties in the associated fields. The default values correspond to a
circular cross section with a diameter of 0.1 m:
This section is only available if User defined is selected as the Cross Section
Definition.
Note
Stress evaluation using only section heights is meaningful only when the cross section
is symmetric.
The max shear stress factor determines the ratio between the peak and the average
shear stress over the cross section as described by Equation 5-1 and Equation 5-2.
Section Orientation
ORIENTATION METHOD
Select the Reference point (the default) or Orientation vector. If Reference point is
selected, enter a Reference point defining local y direction P (SI unit: m).
The local x direction is in the edge direction. The positive edge direction can be
checked by vector plotting the local edge tangent direction. The coordinates of the
reference point define the local xy-plane with the positive y direction defined such that
the point lies in the positive quadrant (see Figure 5-1).
For the creation of a local coordinate system to be possible, the point cannot coincide
with the edge or the edge extension. If this is attempted, an error message is generated.
The default settings for the global coordinates of the point are
[1000,1000,1000]. This is useful only for symmetric cross sections.
Note
Often a number of edges in a plane have the same orientation. It is then easy to select
all edges and specify a point anywhere in the same plane, not coinciding with an edge
or an edge extension.
The Rotation of vector around beam axis has the effect of rotating the given vector
around the beam axis before it is used to define the local xy-plane. This simplifies the
input for some cross sections, such as L-shaped profiles, where the principal axes have
a direction that is skewed relative to a more natural modeling position.
For these loads, right-click and choose Phase to add a phase for harmonic loads in
frequency-domain computations. In this way a harmonic load can be defined where
the amplitude and phase shift can vary with the excitation frequency f:
Edge Load
Add an Edge Load as a force distributed along an edge. Right-click and add a Phase
node for harmonic loads in frequency-domain computations.
3D
Local edge system. This is the system defined by the beam cross section
orientation.
Any additional user-defined coordinate system
FORCE
From the Load type list, select the definition of the force load. Select Load defined as
force per unit length FL (the default) (SI unit: N/m) or Edge load defined as force per
unit volume F (SI unit: N/m3). In the latter case the given load is multiplied by the
cross section area.
Enter values or expressions for the components of the Edge load in the table.
MOMENT
Enter values or expressions for the components of the Moment edge load ML (3D) or
Mlz (2D) (SI unit: Nm/m).
Point Load
Add a Point Load to points for concentrated forces or moments at points. Right-click
to add a Phase node for harmonic loads in frequency-domain computations.
PO IN T S EL EC TIO N
From the Selection list, choose the points to define a point load.
FORCE
Enter values or expressions for the components of the Point load FP (SI unit: N).
MOMENT
Enter values or expressions for the components of the Point moment MP (3D) or Mlz
(2D) (SI unit: Nm).
Phase
Right-click an Edge Load or Point Load node to add a Phase node, which adds a phase
for harmonic loads in frequency-domain computations.
By default, this node inherits the selection from its parent node, and only
a selection that is a subset of the parent nodes selection can be used.
Note
PHASE
Add the phase load Fph (unit: rad) for harmonic loads. Enter the phase for each
component of the load in the corresponding fields.
E D G E , B O U N D A R Y, O R P O I N T S E L E C T I O N
From the Selection list, choose the geometric entity (boundaries, edges, or points) to
prescribe a displacement or rotation.
PAIR SELECTION
If Prescribed Displacement/Rotation is selected from the Pairs menu, choose the pair to
define. An identity pair has to be created first. Ctrl-click to deselect.
For 2D models:
For 3D models:
PRESCRIBED DISPLACEMENT
To define a prescribed displacement for each space direction (x, y, and z),
select one or all of the Prescribed in x direction, Prescribed in y direction,
and Prescribed in z direction check boxes. Then enter a value or expression
3D
for the prescribed displacements u0, v0, or w0 (SI unit: m).
PRESCRIBED ROTATION
To define a prescribed rotation for each space direction (x, y, and z), select
one or all of the Prescribed around x direction, Prescribed around
y direction, and Prescribed around z direction check boxes and enter a value
3D
or expression for in each thx0, thy0, or thz0 (SI unit: rad) field.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Select a Constraint typeBidirectional, symmetric or Unidirectional. If required, select
the Use weak constraints check box.
E D G E , B O U N D A R Y, O R P O I N T S E L E C T I O N
From the Selection list, choose the geometric entity (boundaries, edges, or points) to
prescribe a velocity.
Coordinate systems with directions that change with time should not be
used.
Note
Specify the coordinate system to use for specifying the prescribed translational/
rotational velocity. The coordinate system selection is based on the geometric entity
level. Select from the following based on space dimension:
For 2D models:
2D
Boundary System (a predefined normal-tangential coordinate system)
Any additional user-defined coordinate system (boundaries and points)
For 3D models:
To define a prescribed velocity for each space direction (x and y), select
one or both of the Prescribed in x direction and Prescribed in y direction
check boxes. Then enter a value or expression for the prescribed velocity
2D
components vx and vy (SI unit: m/s).
To define a prescribed velocity for each space direction (x, y, and z), select
one or more of the Prescribed in x direction, Prescribed in y direction, and
Prescribed in z direction check boxes. Then enter a value or expression for
3D
the prescribed velocity components vx, vy, and vz (SI unit: m/s).
P R E S C R I B E D A N G U L A R VE L O C I T Y
To define a prescribed angular velocity for each space direction (x, y, and
z), select one or all of the Prescribed around x direction, Prescribed around
y direction, and Prescribed around z direction check boxes and enter a value
3D
or expression for in each vthx0, vthy0, or vthz0 (SI unit: rad/s) field.
Prescribed Acceleration
The Prescribed Acceleration node adds a boundary or domain condition where the
acceleration is prescribed in one or more directions. The prescribed acceleration
condition is applicable for time-dependent and frequency-domain studies. With this
boundary condition it is possible to prescribe a acceleration in one direction, leaving
the beam free in the other directions.
E D G E , B O U N D A R Y, O R PO I N T S E L E C T I O N
From the Selection list, choose the geometric entity (boundaries, edges, or points) to
prescribe an acceleration.
Coordinate systems with directions which change with time should not
be used.
Note
Specify the coordinate system to use for specifying the prescribed translational/
rotational acceleration. The coordinate system selection is based on the geometric
entity level. Select from the following based on space dimension:
For 2D models:
For 3D models:
PRESCRIBED ACCELERATION
To define a prescribed velocity for each space direction (x and y), select
one or both of the Prescribed in x direction and Prescribed in y direction
check boxes. Then enter a value or expression for the prescribed
2D
acceleration components ax and ay (SI unit: m/s2).
Pinned
The Pinned feature adds an edge condition that makes the edge fixed; that is, all
displacements are zero. This node is also available on boundaries and at points.
E D G E , B O U N D A R Y, O R PO I N T S E L E C T I O N
From the Selection list, choose the geometric entity (boundaries, edges, or points) to
prescribe a pinned condition.
PAIR SELECTION
If Pinned is selected from the Pairs menu, choose the pair to define. An identity pair
has to be created first. Ctrl-click to deselect.
No Rotation
The No Rotation node adds an edge condition that prevents all rotation at the edge.
This node is also available on boundaries and at points.
E D G E , B O U N D A R Y, O R P O I N T S E L E C T I O N
From the Selection list, choose the geometric entity (boundaries, edges, or points) to
prescribe no rotations.
PAIR SELECTION
If No Rotation is selected from the Pairs menu, choose the pair to define. An identity
pair has to be created first. Ctrl-click to deselect.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Select a Constraint typeBidirectional, symmetric or Unidirectional. If required, select
the Use weak constraints check box.
Symmetry
The Symmetry node adds an edge condition that defines a symmetry edge. This node
is also available on boundaries and at points.
PAIR SELECTION
If Symmetry is selected from the Pairs menu, choose the pair to define. An identity pair
has to be created first. Ctrl-click to deselect.
For 2D models:
For 3D models:
Local edge system (the default). This is the system defined by the beam
cross section orientation.
3D Global coordinate system (edges and points)
Any additional user-defined coordinate system (edges and points)
SYMMETRY
Select an Axis to use as normal direction. This specifies the direction of the normal to
the symmetry plane in the selected coordinate system.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Select a Constraint typeBidirectional, symmetric or Unidirectional. If required, select
the Use weak constraints check box.
E D G E , B O U N D A R Y, O R P O I N T S E L E C T I O N
From the Selection list, choose the geometric entity (boundaries, edges, or points) to
prescribe antisymmetry. Select from the following based on space dimension:
For 2D models:
Local edge system (the default). This is the system defined by the beam
cross section orientation.
2D Global coordinate system (boundaries and points)
Any additional user-defined coordinate system (boundaries and points)
For 3D models:
Local edge system (the default). This is the system defined by the beam
cross section orientation.
3D Global coordinate system (edges and points)
Any additional user-defined coordinate system (edges and points)
PAIR SELECTION
If Antisymmetry is selected from the Pairs menu, choose the pair to define. An identity
pair has to be created first. Ctrl-click to deselect.
ANTISYMMETRY
Select an Axis to use as normal direction. This specifies the direction of the normal to
the symmetry plane in the selected coordinate system.
Point Mass
Add a Point Mass node to model a discrete mass or mass moment of inertia that is
concentrated at a point in contrast to distributed mass modeled through the density
and area of the beam. Right-click to add a Point Mass Damping node.
POINT SELECTION
From the Selection list, choose the points to define the point mass.
PAIR SELECTION
If Point Mass is selected from the Pairs menu, choose the pair to define. An identity pair
has to be created first. Ctrl-click to deselect.
For 3D models only, from the Coordinate system list select from:
For 3D models, enter a Mass moment of inertia J (m2kg) in the table for
each axis x, y, and z.
3D
PO IN T S EL EC TIO N
From the Selection list, choose the points to define the point mass damping.
By default, this node inherits the selection from its parent node, and only
a selection that is a subset of the parent nodes selection can be used.
Note
PO IN T MASS DAMPIN G
Enter a Mass damping parameter dM (SI unit: 1/s).
About Beams
In-Plane Euler Beams
3D Euler Beam
Strain-Displacement/Rotation Relation
Stress-Strain Relation
Thermal Strain
Initial Load and Strain
Implementation
Stress Evaluation
Thermal Coupling
Coefficient of Thermal Expansion
Common Cross Sections
About Beams
A beam is a slender structure that can be fully described by the propertiesarea,
moments of inertia, and densityof the cross section. Beams are the choice for
modeling reinforcements in 3D solids and shell structures, as well as in 2D solids under
the plane stress assumption. Naturally, they can also model lattice works, both planar
and three-dimensional.
Beams can sustain loads and moments in any direction, both distributed and on
individual nodes. The beams ends and interconnections can be free, simply supported,
or clamped. In fact, the simplified boundary conditions are usually responsible for
most of the difference that can be found between a beam solution and a full 3D solid
simulation of the same structure. Point constraints on beams are well-behaved, in
contrast to the solid case, and it is possible to use discrete point masses and mass
moments of inertia.
The Beam interface is based on the principle of virtual work. The resulting equation
can equivalently be viewed as a weak formulation of an underlying PDE. The Beam
In-plane Euler beams are defined on edges in 2D. They can be used separately or as
stiffeners to 2D solid elements.
VA R I A BL E S A N D S P A C E D I M E N S I O N S
The degrees of freedom (dependent variables) are the global displacements u and v in
the global x and y directions and the rotation about the global z-axis.
3D Euler beams are defined on edges in 3D. They can be used separately or as
stiffeners to 2D solid or shell elements.
Strain-Displacement/Rotation Relation
The axial strain depends on the rotation derivative and axial displacement derivative
defined by the shape function and the transversal coordinate in the beam. For the 3D
case it becomes
ly lz u axi
= z l ---------- y l ---------- + -------------
s s s
The coordinates from the beam center line in the local transversal directions are
denoted zl and yl respectively. In the 2D case, the second term is omitted, and the local
z direction is always directed out of the plane.
The total strain consists of thermal (th), initial (i), and elastic strains(el)
Stress-Strain Relation
The stress-strain relation in the beam is described by
= E
The stress-strain relation for linear conditions including initial stress and strain and
thermal effects reads:
= E el + i = E th i + i
Thermal Strain
The temperature is assumed to vary linearly across the beams cross section. For the
3D beam it becomes
T = T m + T gz z l + T gy y l
Tm is the temperature at the beam center line while Tgz and Tgy are the temperature
gradients in the two local transversal directions. The thermal strain is thus
th = T m + T gz z l + T gy y l T ref
The initial stress distribution is given as initial moments and initial normal force
N i = i dA
A
M xi = ixz y l iyz z l dA
A
The initial strain is the strain before any loads, displacements, and initial stresses have
been applied. The initial axial strain distribution is given as initial curvature and initial
axial strain
ly lz u axi
i = z l ---------- y l ---------- + -------------
s i s i s i
In 2D the zl dependent term disappears. As initial strain for the torsional degree of
freedom, the derivative of the twist angle,
lx
----------
s i
is used.
Implementation
The implementation is based on the principle of virtual work, which states that the sum
of virtual work from internal strains and external loads equals zero:
T
W = el + u F dV = 0
V
The beam elements are formulated in terms of the stress resultants (normal force,
bending moments and twisting moment).
T m + T gz z l + T gy y l T ref + i dA =
u axi u axi
E -------------
s
------------- T m T ref dA + i dA =
s i
A A
u axi u axi
EA ------------- ------------- T m T ref + N i
s s i
Because the local coordinates are defined with their origin at the centroid of the cross
section, any surface integral of an odd power of a local coordinate evaluates to zero.
M ly = zl dA = zl Eel + i dA =
A A
ly lz u axi ly lz u axi
z l E z l ---------- y l ---------- + ------------- z l ---------- y l ---------- + -------------
s s s s i s i s i
A
T m + T gz z l + T gy y l T ref + i dA =
2 ly ly
zl E ---------- ---------- T gz dA + i z l dA =
s s i
A A
ly ly
EI yy ---------- ---------- T gz + M iy
s s i
T m + T gz z l + T gy y l T ref + i dA =
lz lz
y l E ---------- + ---------- T gy dA i y l dA =
2
s s i
A A
lz lz
EI zz ---------- ---------- + T gy + M iz
s s i
Mly is present only in 3D, and so is the torsional moment Mlx described below. The
torsional stiffness of the beam is defined using the torsional constant J given by
M
J = -------- l
G
In a similar way as for the bending part a torsional moment is then defined as
lx lx
M lx = GJ ---------- ---------- + M ix
s s i
Using the beam moment and normal force the expression for the virtual work becomes
very compact:
ly lz u axi lx
W = Mly ----------
s test
+ M lz ----------
s test
+ N -------------
s test
+ M lx ---------- dx
s test
L
A special feature of some unsymmetrical cross sections is that they twist under a
transversal load that is applied to beam centerline. As an example, this would be the
case for a U-profile under self-weight, loaded in the stiff direction. It is only a load
applied at the shear center, which causes a pure deflection without twist. This effect
can be incorporated by supplying the coordinates of the shear center in the local
coordinate system (ey, ez). A given transversal load (flx, fly, flz), which is defined as
acting along the centerline, is then augmented by a twisting moment given by
m lx = f ly e z f lz e y
The normal stress from axial force is constant over the section, and computed as
N
n = ----
A
The normal stress from bending is computed in four user-selected points (ylk, zlk) in
the cross section as
M ly z lk M lz y lk
bk = ----------------- -----------------
I yy I zz
In 2D, only two points, specified by their local y-coordinates are used.
k = bk + n
max = max k
The shear stress from twist in general has a complex distribution over the cross section.
The maximum shear stress due to torsion is defined as
M lx
t, max = -------------
Wt
where Wt is the torsional section modulus. This result is available only in 3D.
2
ly
T lz = EI yy ------------
2
-
s
2
lz
T ly = E I zz ------------
2
-
s
T lz
sz, ave = --------
A
(5-1)
T ly
sy, ave = --------
A
Since the shear stresses are not constant over the cross section, the maximum shear
stresses are also available, using section dependent correction factors:
As the directions and positions of maximum shear stresses from shear and twist are not
known in a general case, upper bounds to the shear stress components are defined as
The maximum von Mises effective stress for the cross section is the defined as
2 2 2
mises = max + 3 xy, max + 3 xz, max
Since the maximum values for the different stress components in general occur at
different positions over the cross section, the effective stress thus computed is a
conservative approximation.
Thermal Coupling
Material expands with temperature, which causes thermal strains to develop in the
material. The beams can handle any temperature variation along the beam, and linear
variation across the beam. The thermal strains together with the initial strains and
elastic strains from structural loads form the total strain.
= el + th + i
where
where th is the thermal strain, Ty and Tz are the temperature difference over the
cross section of the beam in the y and z directions, and is the coefficient of thermal
expansion. Tm is the temperature in the middle and Tref is the stress free reference
temperature.
RECTANGULAR SECTION
Figure 5-2: Rectangular section diagram for common cross section. The diagram also
displays in COMSOL Multiphysics when this option is selected.
A hy hz
Izz 3
hy hz
------------
12
ez 0
y 1.5
Iyy 3
hz hy
------------
-
12
ey 0
z 1.5
J 2 q=
h y h z q
2 2
192q ------------------------
2n 1 -
----------------- 1
3 5
------------------------------
2n 1
5
tanh
2q
min(hy/hz,hz/hy)
n=1
Wt J q=
--------------------------------------------------------------------------------------------------------------------------
-
2 min(hy/hz,hz/hy)
8
hy hz q 1
2 2 2n 1
-
--------------------------------------------------------------------------
n = 1 2n 1 cosh -------------------------
2q
p1 h y h z
--------
2 2
- ---------
p2
h h
-----y- --------z-
2 2
p3
h h
-----y- -----z-
2 2
p4 h y h z
--------
2 2
- ------
Figure 5-3: Box section diagram for common cross section. The diagram also displays in
COMSOL Multiphysics when this option is selected.
A 2 hy tz + hz ty 4 ty tz
Izz 3 3 2
t z h y + t y h z 2t z t y h z 2t z h y t y
------------------------------------------------ + -------------------------------------------------------
6 2
ez 0
y 2 2
h y h z h y 2t y h z 2t z A
-------------------------------------------------------------------------------
-
16t z I zz
Iyy 3 3 2
t y h z + t z h y 2t y t z h y 2t y h z t z
------------------------------------------------- + -------------------------------------------------------
6 2
ey 0
z 2 2
h z h y h z 2t z h y 2t y A
-------------------------------------------------------------------------------
-
16t y I yy
J 2 2 Thin-walled
2 hy ty hz tz
----------------------------------------------------- approximation
hy ty hz tz
----------------- + ----------------
tz ty
Wt 2 h y t y h z t z min t y t z Thin-walled
approximation
p1 h y h z
--------
- ---------
2 2
p2
h h
-----y- --------z-
2 2
p3
h h
-----y- -----z-
2 2
p4 h y h z
--------
- ------
2 2
CIRCULAR SECTION
Figure 5-4: Circular section diagram for common cross section. The diagram also displays
in COMSOL Multiphysics when this option is selected.
A 2
d o
----------
4
Izz 4
d
---------o-
64
ez 0
y 4
---
3
Iyy I zz
ey 0
z y
J 4
d
---------o-
32
Wt 3
d
---------o-
16
p1 do
--------
- 0
2
p2 d o
0 --------
-
2
p3
d
-----o-
2 0
p4
0 d
-----o-
2
PIPE SECTION
Figure 5-5: Pipe section diagram for common cross section. The diagram also displays in
COMSOL Multiphysics when this option is selected.
A 2 2
do di
---------------------------
4
Izz 4 4
do di
---------------------------
64
ez 0
y 3 3
d o d i A
------------------------------------
-
12 d o d i I zz
Iyy I zz
ey 0
z y
J
4 4
do di
---------------------------
32
Wt 4 4
do di
---------------------------
16d o
p1 do
--------
2 0
-
p2 d o
0 --------
-
2
p3 do
-----
- 0
2
p4
0 d
-----o-
2
Figure 5-6: H-profile section diagram for common cross section. The diagram also displays
in COMSOL Multiphysics when this option is selected.
A 2h z t y + t z h y 2t y
Izz 3 3 2
2h z t y + t z h y 2t y ty hz hy ty
- + -----------------------------------
------------------------------------------------------ -
12 2
ez 0
y 2
4h z t y h y t y + t z h y 2t y A
--------------------------------------------------------------------------------------
8t z I zz
Iyy 3 3
2t y h z + t z h y 2t y
----------------------------------------------------
12
ey 0
z 2 2
h z t z A
--------------------------
8I yy
J 3 3 Thin-walled
2t y h z + t z h y 2t y
---------------------------------------------------- approximation
3
Wt J Thin-walled
-----------------------------
-
max t y t z approximation
p1 h y h z
--------
- ---------
2 2
p2
h h
-----y- --------z-
2 2
p3
h h
-----y- -----z-
2 2
p4 h y h z
--------
- ------
2 2
U-PROFILE SECTION
Figure 5-7: U-profile section diagram for common cross section. Also displays in COMSOL
Multiphysics when this option is selected.
A h y t z + 2 h z t z t y
Izz 3 3 2
t z h y + 2t y h z t z t y h z t z h y t y
- + ----------------------------------------------------
-----------------------------------------------
12 2
zCG 2
h y t z 2h z t z + 2t y h z t z
---------------------------------------------------------------------------
-
2A
ez 2 2
t z t y 2h z t z h y t y
h z ---- + ---------------------------------------------------------- z CG
2 16I zz
y 2 2
h y h z h y 2t y h z t z A
---------------------------------------------------------------------------
-
8t z I zz
Iyy 3 3 2
8t y h z + t z h y 2t y + 3t z h y 2t y 2h z t z
-------------------------------------------------------------------------------------------------------------------------
-
12
2
z CG A
ey 0
z 2
z CG A
-----------------
2I yy
J 3 3 Thin-walled
2t y h z + t z h y 2t y
---------------------------------------------------- approximation
3
Wt J Thin-walled
-----------------------------
-
max t y t z approximation
p1 hy
--------
- z CG
2
p2
h
-----y- z
2 CG
p3
h
-----y- h
2 z z CG
p4 hy
--------
- h z z CG
2
Figure 5-8: T-profile section diagram for common cross section. Also displays in COMSOL
Multiphysics when this option is selected.
A h z t y + h y t y t z
yCG 2
t z h y t y + t y h z 2h y t y
-------------------------------------------------------------------------
2A
Izz 3 3 2
4t z h y t y + h z t y + 3t y h z 2h y t y 2
----------------------------------------------------------------------------------------------------- y CG A
12
ez 0
y 2
y CG A
------------------
2I zz
Iyy 3 3
ty hz + tz hy ty
---------------------------------------------
12
ey ty
h y ---- y CG
2
z 2 2
h z t z A
--------------------------
8I yy
J 3 3 Thin-walled
ty hz + tz hy ty
--------------------------------------------- approximation
3
Wt J Thin-walled
---------------------------
-
max t y t z approximation
p1 t z
y -------
CG 2
p2 h z
h y --------
-
y CG 2
p3
h y h
-----z-
y CG 2
p3
y t---z-
CG 2
Trusses
This chapter describes the Truss interface, which you find under the Structural
Mechanics branch ( ) in the Model Wizard.
In this chapter:
303
The Truss Interface
3D
The Truss interface ( ) is found under the Structural Mechanics branch ( ) in the
Model Wizard. Trusses are elements that can only sustain axial forces; therefore use
trusses to model truss works where the edges are straight but also to model sagging
cables like the deformation of a wire exposed to gravity.
The Linear Elastic Material node is the default material model, which adds a linear elastic
equation for the displacements and has a settings window to define the elastic material
properties.
When this interface is added, these default nodes are also added to the Model Builder:
Linear Elastic Material, Cross Section Data, Free (a condition where points are free, with
no loads or constraints), Straight Edge Constraint (to ensure that the points lie on a
straight line between the end points of the edge or boundary), and Initial Values.
Right-click the Truss node to add other features that implement, for example, loads and
constraints. The following sections provide information about all features available in
this interface.
INTERFACE IDENTIFIER
The interface identifier is a text string that can be used to reference the respective
physics interface if appropriate. Such situations could occur when coupling this
interface to another physics interface, or when trying to identify and use variables
304 | C H A P T E R 6 : TR U S S E S
defined by this physics interface, which is used to reach the fields and variables in
expressions, for example. It can be changed to any unique string in the Identifier field.
The default identifier (for the first interface in the model) is truss.
DISCRETIZATION
To display this section, click the Show button ( ) and select Discretization. Select
Quadratic (the default), Linear, Cubic, or Quartic for the Displacement field. Specify the
Value type when using splitting of complex variablesReal or Complex (the default).
T H E TR U S S I N T E R F A C E | 305
Boundary, Edge, Point, and Pair Conditions for the Truss Interface
The Truss Interface has these boundary, edge, point, and pair conditions available as
indicated.
These are described in this section (listed in alphabetical order). The list also includes
subfeatures:
Antisymmetry
Cross Section Data
Edge Load
Initial Stress and Strain
Initial Values
Linear Elastic Material
Phase
Pinned
Prescribed Acceleration (for time-dependent and frequency-domain studies)
Prescribed Displacement
Prescribed Velocity (for time-dependent and frequency-domain studies)
Point Mass
Point Mass Damping
Straight Edge Constraint
Symmetry
Thermal Expansion
Added Mass
Damping (described for the Beam interface)
Fixed Constraint
Free
306 | C H A P T E R 6 : TR U S S E S
Point Load
Pre-Deformation
Spring Foundation
MODEL INPUTS
Define model inputs, for example, the temperature field of the material uses a
temperature-dependent material property. If no model inputs are required, this section
is empty.
These settings are the same as described under Linear Elastic Material for
The Beam Interface.
Note
T H E TR U S S I N T E R F A C E | 307
GEOMETRIC NONLINEARITY
In this section there is always one check box. Either Force linear strains or Include
geometric nonlinearity is shown.
If a study step is geometrically nonlinear, the default behavior is to use a large strain
formulation in all domains. There are however some cases when you would still want
to use a small strain formulation for a certain domain. In those cases, select the Force
linear strains check box. When selected, a small strain formulation is always used,
independently of the setting in the study step. The default value is that the check box
is cleared, except when opening a model created in a version prior to 4.3. In this case
the state is chosen so that the properties of the model are conserved.
The Include geometric nonlinearity check box is displayed only if the model was created
in a version prior to 4.3, and geometric nonlinearity was originally used for the selected
domains. It is then selected and forces the Include geometric nonlinearity check box in
the study step to be selected. If the check box is cleared, the check box is permanently
removed and the study step assumes control over the selection of geometric
nonlinearity.
Thermal Expansion
Right-click the Linear Elastic Material node to add the Thermal Expansion node. Thermal
expansion is an internal thermal strain caused by changes in temperature. The
temperature is assumed to be constant over the cross section of the truss element.
The settings for the Truss interface are the same as described for the Beam
interface (excluding the thermal bending options). See Thermal
Note Expansion.
308 | C H A P T E R 6 : TR U S S E S
BOUNDARY OR EDGE SELECTION
From the Selection list, choose the geometric entity (boundaries or edges) to define the
initial stress or strain.
By default, this node inherits the selection from its parent node, and only
a selection that is a subset of the parent nodes selection can be used.
Note
Use the Cross Section Data node to enter the cross section area for the truss elements.
Initial Values
The Initial Values feature adds an initial value for the displacement field u (the
displacement components u, v, and, in 3D, w), and the velocity field. It can serve as
an initial condition for a transient simulation or as an initial guess for a nonlinear
analysis. Right-click to add additional Initial Values features.
T H E TR U S S I N T E R F A C E | 309
BOUNDARY OR EDGE SELECTION
From the Selection list, choose the geometric entity (boundaries or edges) to define
initial values.
INITIAL VALUES
Enter values or expressions in the applicable fields based on the space dimension. The
default value is 0 for all initial values of the Displacement field u (SI unit: m) and the
du
Velocity field (SI unit: m/s).
dt
For these loads, right-click and choose Phase to add a phase for harmonic loads in
frequency-domain computations. In this way, define a harmonic load where the
amplitude and phase shift can vary with the excitation frequency f:
Edge Load
Add an Edge Load as a force distributed along an edge (3D models) or boundary (2D
models). Also right-click and add a Phase for harmonic loads in frequency-domain
computations.
FORCE
From the Load type list, select the definition of the force load. Select Load defined as
force per unit length FL (the default) (SI unit: N/m) or Edge load defined as force per
unit volume F (SI unit: N/m3). In the latter case the given load intensity is multiplied
310 | C H A P T E R 6 : TR U S S E S
by the cross section area. Enter values or expressions for the components of the Edge
load in the table.
Phase
Right-click an Edge Load or Point Load node to add a Phase node, which adds a phase
for harmonic loads in frequency-domain computations.
By default, this node inherits the selection from its parent node, and only
a selection that is a subset of the parent nodes selection can be used.
Note
PHASE
Add the phase load Fph (unit: rad) for harmonic loads. Enter the phase for each
component of the load in the corresponding fields.
Pinned
The Pinned node adds an edge (3D), boundary (2D), or point (2D and 3D) condition
that makes the edge, boundary, or point fixed; that is, the displacements are zero in all
directions.
T H E TR U S S I N T E R F A C E | 311
B O U N D A R Y, E D G E , O R PO I N T S E L E C T I O N
From the Selection list, choose the geometric entity (boundaries, edges, or points) that
are fixed.
PAIR SELECTION
If Pinned is selected from the Pairs menu, choose the pair to define. An identity pair
has to be created first. Ctrl-click to deselect.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Select a Constraint typeBidirectional, symmetric or Unidirectional. If required, select
the Use weak constraints check box.
Prescribed Displacement
The Prescribed Displacement node adds an edge (3D), boundary (2D), or point (2D
and 3D) condition where the displacements are prescribed in one or more directions.
B O U N D A R Y, E D G E , O R PO I N T S E L E C T I O N
From the Selection list, choose the geometric entity (boundaries, edges, or points) to
prescribe a displacement or rotation.
PAIR SELECTION
If Prescribed Displacement is selected from the Pairs menu, choose the pair to define.
An identity pair has to be created first. Ctrl-click to deselect.
312 | C H A P T E R 6 : TR U S S E S
PRESCRIBED DISPLACEMENT
To define a prescribed displacement for each space direction (x, y, and z),
select one or all of the Prescribed in x, y, and z direction check boxes. Then
enter a value or expression for the prescribed displacements u0, v0, or w0
3D
(SI unit: m).
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Select a Constraint typeBidirectional, symmetric or Unidirectional. If required, select
the Use weak constraints check box.
Prescribed Velocity
The Prescribed Velocity node adds an edge, boundary, or point condition that
prescribes the velocity in one or more directions. The prescribed velocity condition is
applicable for time-dependent and frequency-domain studies. With this condition it is
possible to prescribe a velocity in one direction, leaving the truss free in the other
directions.
B O U N D A R Y, E D G E , O R P O I N T S E L E C T I O N
From the Selection list, choose the geometric entity (boundaries, edges, or points) to
prescribe a velocity.
T H E TR U S S I N T E R F A C E | 313
selection is based on the geometric entity level. Coordinate systems with directions
which change with time should not be used.
P R E S C R I B E D VE L O C I T Y
To define a prescribed velocity for each space direction (x and y), select
one or both of the Prescribed in x direction and Prescribed in y direction
check boxes. Then enter a value or expression for the prescribed velocity
2D
components vX and vY (SI unit: m/s).
To define a prescribed velocity for each space direction (x, y, and z), select
one or more of the Prescribed in x, y, and z direction check boxes. Then
enter a value or expression for the prescribed velocity components vX, vY,
3D
and vZ (SI unit: m/s).
Prescribed Acceleration
The Prescribed Acceleration node adds a boundary, edge, or point condition where the
acceleration is prescribed in one or more directions. The prescribed acceleration
condition is applicable for time-dependent and frequency-domain studies. With this
boundary condition it is possible to prescribe a acceleration in one direction, leaving
the truss free in the other directions.
B O U N D A R Y, E D G E , O R PO I N T S E L E C T I O N
From the Selection list, choose the geometric entity (boundaries, edges, or points) to
prescribe an acceleration.
PRESCRIBED ACCELERATION
314 | C H A P T E R 6 : TR U S S E S
Select one or more of the Prescribed in x, y, and z direction check boxes
then enter a value or expression for the prescribed acceleration
3D components aX, aY, and aZ (SI unit: m/s2).
Symmetry
The Symmetry node adds an edge (3D), boundary (2D), or point (2D and 3D)
condition that defines a symmetry edge, boundary, or point.
B O U N D A R Y, E D G E , O R P O I N T S E L E C T I O N
From the Selection list, choose the geometric entity (boundaries, edges, or points) to
prescribe symmetry.
PAIR SELECTION
If Symmetry is selected from the Pairs menu, choose the pair to define. An identity pair
has to be created first. Ctrl-click to deselect.
SYMMETRY
Select an Axis to use as normal direction. This specifies the direction of the normal to
the symmetry plane in the selected coordinate system.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Select a Constraint typeBidirectional, symmetric or Unidirectional. If required, select
the Use weak constraints check box.
T H E TR U S S I N T E R F A C E | 315
Antisymmetry
The Antisymmetry node adds an edge (3D), boundary (2D), or point (2D and 3D)
condition that defines an antisymmetry edge, boundary, or point.
B O U N D A R Y, E D G E , O R PO I N T S E L E C T I O N
From the Selection list, choose the geometric entity (boundaries, edges, or points) to
prescribe antisymmetry.
PAIR SELECTION
If Antisymmetry is selected from the Pairs menu, choose the pair to define. An identity
pair has to be created first. Ctrl-click to deselect.
ANTISYMMETRY
Select an Axis to use as normal direction. This specifies the direction of the normal to
the symmetry plane in the selected coordinate system.
CONSTRAINT SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Select a Constraint typeBidirectional, symmetric or Unidirectional. If required, select
the Use weak constraints check box.
Point Mass
The Point Mass node adds a discrete mass that is concentrated at a point in contrast to
distributed mass modeled through the density and area of the truss element.
Right-click to add a Point Mass Damping subnode.
POINT SELECTION
From the Selection list, choose the points to define the point mass.
316 | C H A P T E R 6 : TR U S S E S
PAIR SELECTION
If Point Mass is selected from the Pairs menu, choose the pair to define. An identity pair
has to be created first. Ctrl-click to deselect.
PO IN T MA SS
Enter a Point mass m (SI unit: kg).
PO IN T S EL EC TIO N
From the Selection list, choose the points to define the point mass damping.
By default, this node inherits the selection from its parent node, and only
a selection that is a subset of the parent nodes selection can be used.
Note
PO IN T MASS DAMPIN G
Enter a Mass damping parameter dM (SI unit: 1/s).
T H E TR U S S I N T E R F A C E | 317
Theory for the Truss Interface
The Truss Interface theory is described in this section:
About Trusses
Theory Background for the Truss Interface
About Trusses
Truss elements are elements that can only sustain axial forces. They have displacements
as degrees of freedom. Truss elements are sometimes referred to as bars or spars. They
live on boundaries in 2D and edges in 3D. The Truss interface supports the same study
types as the Solid Mechanics interface. Use trusses to model truss works where the
edges are straight but also to model sagging cables like the deformation of a wire
exposed to gravity. In such applications trusses are often referred to as cable elements.
I N - P L A N E TR U S S
Use the Truss interface in 2D to analyze planar lattice trusses or sagging cable-like
structures. The Truss interface is defined on edges in 2D.
318 | C H A P T E R 6 : TR U S S E S
TR U S S ( 3 D )
Use the Truss interface to model three-dimensional trusses or sagging cable-like
structures. The Truss interface is defined on edges in 3D.
STRAIN-DISPLACEMENT RELATION
The axial strain n is calculated by expressing the global strains in tangential derivatives
and projecting the global strains on the edge.
t
n = t gT t
xT xyT xzT
gT = xyT yT yzT
xzT yzT zT
T H E O R Y F O R T H E TR U S S I N T E R F A C E | 319
The strains can be expressed as either engineering strains for small displacements or
Green-Lagrange strains for large displacements. The Green-Lagrange strain tensor
used for large displacements is defined as
1 u i u j u k u k
ijT = --- + +
2 xj xi xi xj
T T T T
1 u i u j
ijT = --- + (6-1)
2 xj xi
T T
n = t x xT t x + xyT t y + xzT t z +
t y xyT t x + yT t y + yzT t z +
t z xzT t x + yzT t y + zT t z
STRESS-STRAIN RELATION
The constitutive relation for the axial stress including thermal strain and initial stress
and strain is
n = E n T T ref ni + ni (6-2)
S n = E n T Tref ni + S ni
For output, the First Piola-Kirchhoff stress Pn is then computed from the Second
Piola-Kirchhoff stress using
P n = S n s'
where s is the ratio between current and initial length. The axial force in the element
is then computed as
N = Pn A0
320 | C H A P T E R 6 : TR U S S E S
A0
n = P n -------
A
where A is the deformed area of the element. For a linearly elastic material, the area
change is
A0 2
------- = 1 n
A
This is the only occasion where the Truss interface uses the Poissons ratio .
In a geometrically linear analysis all the stress representations have the same value, as
defined by Equation 6-2.
THERMAL COUPLING
Material expands with temperature, which causes thermal strains to develop in the
material. The trusses can handle any temperature variation along the truss. The
thermal strains together with the initial strains and elastic strains from structural loads
form the total strain.
= el + th + i
where
th = T Tref
Thermal coupling means that the thermal expansion is included in the analysis.
IMPLEMENTATION
Using the principle of virtual work results in the following weak formulation
u FPi
t t
W = d n n + u F V dV +
V i
where the summation stands for summation over all points in the geometry. Replacing
the integration over the cross section with the cross-sectional area (A) and the volume
forces with line forces, the equation becomes
T H E O R Y F O R T H E TR U S S I N T E R F A C E | 321
STRAIGHT EDGE OPTION
The optional constraint to enforce the nodes to lie on the straight line between the end
points of the edge are formulated as follows:
Starting with the large displacement case, let xd1 and xd2 be the deformed position of
the two end points of the edge
x di = u i + x i (6-3)
x + u = x d1 + ta (6-4)
where t is a parameter along the line, and a is the direction vector for the line. a is
calculated from the deformed position of the end points as
a = x d2 x d1
The constraints for the edge is derived by substituting the parameter t from one of the
scalar equations in Equation 6-4 into the remaining ones. In 2D the constraint
equations become
x + u x d1 a y y + v y d1 a x
x + u x d1 a z z + w z d1 a x
y + v y d1 a z z + w z d1 a y
To avoid problems when the edge is directed in one of the coordinate axes directions,
a third constraint is added. This constraint is a linear combination of the two earlier
constraints:
y + v y d1 a x x + u x d1 a y
A linear constraint is needed in order for the solution of the small displacement
problem to become independent of the solver. The linear relation for the displacement
is
u 1 x n2 x n + u 2 x n x n1
u = ------------------------------------------------------------------------ + u ax x 2 x 1 (6-5)
x n2 x n1
322 | C H A P T E R 6 : TR U S S E S
where uax is the axial displacement along the edge, and xn are a linear parameter along
the edge
x x 2 x1 + y y2 y 1 + z z2 z1
x n = --------------------------------------------------------------------------------------------
2 2 2
x2 x1 + y2 y1 + z2 z1
Eliminating uax from Equation 6-5 results in the following linear constraint in 2D
u 1 x n2 x n + u 2 x n x n1
------------------------------------------------------------------------ u y 2 y 1
x n2 x n1
v 1 x n2 x n + v 2 x n x n1
---------------------------------------------------------------------- v x 2 x 1 = 0
x n2 x n1
u 1 x n2 x n + u 2 x n x n1
------------------------------------------------------------------------ u z 2 z 1
x n2 x n1
w 1 x n2 x n + w 2 p x n1
----------------------------------------------------------------------- w x 2 x 1 = 0
x n2 x n1
v 1 x n2 x n + v 2 x n x n1
---------------------------------------------------------------------- v z 2 z 1
x n2 x n1
(6-6)
w 1 x n2 x n + w 2 x n x n1
-------------------------------------------------------------------------- w y 2 y 1 = 0
x n2 x n1
v 1 x n2 x n + v 2 x n x n1
---------------------------------------------------------------------- v x 2 x 1
x n2 x n1
u 1 x n2 x n + u 2 x n x n1
------------------------------------------------------------------------ u y 2 y 1 = 0
x n2 x n1
T H E O R Y F O R T H E TR U S S I N T E R F A C E | 323
324 | C H A P T E R 6 : TR U S S E S
7
Membranes
This chapter describes the Membrane interface, which you find under the
Structural Mechanics branch ( ) in the Model Wizard.
In this chapter:
325
The Membrane Interface
The Membrane interface ( ) is found under the Structural Mechanics branch ( ) in
the Model Wizard. Membranes can be considered as plane stress elements in 3D with a
possibility to deform both in the in-plane and out-of-plane directions. The difference
between a shell and a membrane is that the membrane does not have any bending
stiffness. When a membrane is used by itself, a tensile prestress is necessary in order to
avoid singularity, since a membrane with no stress or compressive stress has no
transverse stiffness.
This interface can be used to model prestressed membranes, but can also be used to
model a thin cladding on solids.
Since prestress is used in almost all cases, the interface is formulated so that it always
includes geometric nonlinearity. This property does not interact with the settings for
geometric nonlinearity within the study settings, which apply to other interfaces under
the Structural Mechanics branch.
The Linear Elastic Material node is the only available material model. It adds a linear
elastic equation for the displacements and has a settings window to define the elastic
material properties.
When this interface is added, these default nodes are also added to the Model Builder:
Linear Elastic Material, Free (a condition where edges are free, with no loads or
constraints), and Initial Values.
Right-click the Membrane node to add other features that implement, for example,
loads and constraints. The following sections provide information about all features
available in this interface.
INTERFACE IDENTIFIER
The interface identifier is a text string that can be used to reference the respective
physics interface if appropriate. Such situations could occur when coupling this
interface to another physics interface, or when trying to identify and use variables
The default identifier (for the first interface in the model) is mem.
BOUNDARY SELECTION
The default setting is to include All boundaries in the model to define the membrane.
To choose specific boundaries, select Manual from the Selection list.
S T R U C T U R A L TR A N S I E N T B E H AV I O R
Select an option from the Structural transient behavior listInclude inertial terms (the
default) or Quasi-static.
THICKNESS
Define the Thickness h by entering a value or expression (SI unit: m) in the field. The
default is 0.0001 m. Use the Change Thickness feature to define a different thickness
in parts of the membrane. The thickness can be variable if an expression is used.
DISCRETIZATION
To display this section, click the Show button ( ) and select Discretization. Select
Quadratic (the default), Linear, Cubic, or Quartic for the Displacement field. Specify the
Value type when using splitting of complex variablesReal or Complex (the default).
Edge Load
Face Load
Initial Stress and Strain
Initial Values
Linear Elastic Material
Prescribed Displacement
These features are described for the Solid Mechanics interface. The list also includes
subfeatures:
Added Mass
Body Load
Change Thickness
Damping
Fixed Constraint
Free
Phase
Point Load
Pre-Deformation
Spring Foundation
Thermal Expansion
BOUNDARY SELECTION
From the Selection list, choose the boundaries to define a linear elastic membrane and
compute the displacements, stresses, and strains.
Youngs modulus and Poissons ratio to specify Youngs modulus (elastic modulus)
E (SI unit: Pa) and Poissons ratio (dimensionless). For an isotropic material
Youngs modulus is the spring stiffness in Hookes law, which in 1D form is E
where is the stress and is the strain. Poissons ratio defines the normal strain in
the perpendicular direction, generated from a normal strain in the other direction
and follows the equation = .
Youngs modulus and shear modulus to specify Youngs modulus (elastic modulus)
E (SI unit: Pa) and the shear modulus G (SI unit: Pa). For an isotropic material
Youngs modulus is the spring stiffness in Hookes law, which in 1D form is E
where is the stress and is the strain. The shear modulus is a measure of the solids
resistance to shear deformations.
Bulk modulus and shear modulus to specify the bulk modulus K (SI unit: Pa) and the
shear modulus G (SI unit: Pa). The bulk modulus is a measure of the solids
resistance to volume changes. The shear modulus is a measure of the solids
resistance to shear deformations.
Lam constants to specify the first and second Lam constants (SI unit: Pa) and
(SI unit: Pa).
Pressure-wave and shear-wave speeds to specify the pressure-wave speed
(longitudinal wave speed) cp (SI unit: m/s) and the shear-wave speed (transverse
wave speed) cs (SI unit: m/s).
This is the wave speed for a solid continuum. In plane stress, for example,
the actual speed with which a longitudinal wave travels is lower than the
Note value given.
For each pair of properties, select from the applicable list to use the value From material
or enter a User defined value or expression.
By default, this node inherits the selection from its parent node, and only
a selection that is a subset of the parent nodes selection can be used.
Note
Initial Values
The Initial Values feature adds an initial value for the displacement field and the
structural velocity field. It can serve as an initial condition for a transient simulation or
as an initial guess for a nonlinear analysis. Right-click to add additional Initial Values
features.
BOUNDARY SELECTION
From the Selection list, choose the boundaries to define initial values.
IN IT IA L VA LUES
Enter values or expressions in the applicable fields based on the space dimension. The
default value is 0 for all initial values of the Displacement field u (SI unit: m) and
du
Structural velocity field (SI unit: m/s).
dt
Face Load
Add a Face Load to boundaries to use it as a pressure or tangential force acting on a
surface. Right-click and add a Phase for harmonic loads in frequency-domain
computations.
FORCE
Select a Load typeLoad defined as force per unit area, Total force, or Pressure.
If Load defined as force per unit area is selected, enter values or expressions for the
components of the Load FA (SI unit: N/m2).
If Pressure is selected, enter a value or expression for the Pressure p (SI unit: Pa). A
positive pressure is directed in the negative shell normal direction.
If Total force is selected, enter values or expressions in the components of the Total
force Ftot (SI unit: N).
Edge Load
Add an Edge Load as a force distributed along an edge. Also right-click and add a Phase
for harmonic loads in frequency-domain computations.
EDGE SELECTION
From the Selection list, choose the edges to define an edge load.
Prescribed Displacement
The Prescribed Displacement feature adds a condition where the displacements are
prescribed in one or more directions to the geometric entity (boundary, edge, or
point).
B O U N D A R Y, E D G E , O R P O I N T S E L E C T I O N
From the Selection list, choose the geometric entity (boundaries, edges, or points) to
prescribe a displacement.
PAIR SELECTION
If Prescribed Displacement is selected from the Pairs menu, choose the pair to define.
An identity pair has to be created first. Ctrl-click to deselect.
CONSTRAINT SETTINGS
To display this section, select click the Show button ( ) and select Advanced Physics
Options. Select a Constraint typeBidirectional, symmetric or Unidirectional. If required,
select the Use weak constraints check box.
About Membranes
Theory Background for the Membrane Interface
About Membranes
Membranes are can be considered as plane stress elements in 3D with a possibility to
deform both in the in-plane and out-of-plane directions. The difference between a
shell and a membrane is that the membrane does not have any bending stiffness. The
Membrane Interface supports the same study types as the Solid Mechanics interface
except it does not include the Linear Buckling study type.
To describe a membrane, provide its thickness and the elastic material properties. All
properties may be variable over the element. All elemental quantities are integrated
only at the midsurface and this is a good approximation since by definition a membrane
is thin.
MEMBRANES - 3D MODELS
The Membrane interface in 3D can be active on internal and external boundaries of a
domain, as well as on boundaries not adjacent to any domain.
The dependent variables are the displacements u, v, and w in the global x, y, and z
directions, and the displacement derivative unn in the direction normal to the
membrane.
The dependent variables are the displacements u and w in the global r and z directions,
and the displacement derivative unn in the direction normal to the membrane in r-z
plane.
The thickness of the membrane is h, which can vary over the element. The
displacements are interpolated by the n:th order Lagrange basis functions.
STRAIN-DISPLACEMENT RELATION
The local tangential strains in the membrane are calculated by expressing the global
strains in tangential derivatives and projecting the global strains on the membrane.
t
= t gT t (7-1)
xT xyT xzT
gT = xyT yT yzT
xzT yzT zT
rT 0 rzT
gT = u-
0 --- 0
R
rzT 0 zT
The strains can are expressed as Green strains for large displacements. The Green strain
tensor is defined as
1 u i u j u k u k
ijT = --- + +
2 xj xi xi xj
T T T T
The local strain tensor (3x3) is defined using local in-plane strain tensor () and given
by
xT xyT 0
l =
xyT yT 0
0 0 en
The thermal strains and initial stresses-strains (only in in-plane directions of the
membrane) are added in the constitutive relation in a similar manner as it is done in
linear elastic continuum mechanics.
The weak equations in the Membrane interface are similar to that of linear elastic
continuum mechanics.
EXTERNAL LOADS
Contributions to the virtual work from the external load are of the form
u test F
F = pn
For a follower load, the change in midsurface area is taken into account, and
integration of the load is done in spatial frame.
STRESS CALCULATIONS
The strains are calculated in the element as described above The stresses are computed
by applying the constitutive law to the computed strains.
The membrane does not support transverse and bending forces and the only section
forces it support is the membrane force and defined as:
N = hs
where s is the local stress tensor and contains only in-plane stress components.
This local membrane surface coordinate system is defined by the boundary coordinate
system (t1,t2,n).
The quantities like stresses and strains are also available in global coordinate systems
after doing a transformation from a local (boundary) to a global coordinate system.
Multiphysics Interfaces
If you also have the Subsurface Flow Module or Geomechanics Module, the
Poroelasticity interface is also available. See the Subsurface Flow Module Users
Guide or Geomechanics Module Users Guide for details.
339
If you also have the Acoustics Module, the Acoustic-Solid Interaction, Acoustic-Shell
Interaction, and Acoustic-Piezoelectric Interaction interfaces are also available. See the
Acoustics Module Users Guide for details.
When this interface is added, these default nodes are also added to the Model Builder
Thermal Linear Elastic Material (the default material model, which adds a linear elastic
equation for the displacements and the heat equation for the temperature), Free (a
boundary condition where boundaries are free, with no loads or constraints), Thermal
Insulation, and Initial Values. Right-click the Thermal Stress node to add other features
that implement, for example, loads, constraints, heat sources, and nonlinear material
models.
INTERFACE IDENTIFIER
The interface identifier is a text string that can be used to reference the respective
physics interface if appropriate. Such situations could occur when coupling this
interface to another physics interface, or when trying to identify and use variables
defined by this physics interface, which is used to reach the fields and variables in
expressions, for example. It can be changed to any unique string in the Identifier field.
The default identifier (for the first interface in the model) is ts.
2D APPROXIMATION
From the 2D approximation list select Plane stress or Plane strain. When
modeling using plane stress, the interface solves for the out-of-plane strain
2D components in addition to the displacement field u.
THICKNESS
Enter a value or expression for the Thickness d (SI unit: m). The default
value of 1 m is suitable for plane strain models, where it represents a a
unit-depth slice, for example. For plane stress models, enter the actual
thickness, which should be small compared to the size of the plate for the
plane stress assumption to be valid. In rare cases, when thickness is
2D
changed in parts of the geometry; then use the Change Thickness feature.
This thickness also controls the thickness dz, active in the separate Heat
Transfer interface for 2D out-of-plane heat transfer.
S T R U C T U R A L TR A N S I E N T B E H AV I O R
From the Structural transient behavior list, select Quasi-static or Include inertial terms to
treat the elastic behavior as quasi-static (with no mass effects; that is, no second-order
time derivatives) or as a mechanical wave in a time-dependent study. The default is to
use the quasi-static behavior for a time-dependent thermal stress study. Select Include
inertial terms to model the structural behavior in a time-dependent study as a
mechanical wave.
TY P I C A L WA V E S P E E D F O R P E R F E C T L Y M A T C H E D L A Y E R S
The typical wave speed cref is a parameter for the perfectly matched layers (PMLs) if
used in a solid wave propagation model. The default value is ts.cp, the pressure-wave
PHYSICAL MODEL
If you also have the Heat Transfer Module, the out-of-plane heat transfer,
surface-to-surface radiation, and radiation in participating media options
are available in this section and described in the Heat Transfer Module
Tip
Users Guide.
RADIATION SETTINGS
If you have the Heat Transfer Module, then this section can be made
available. To display this section select the Surface-to-surface radiation
Note check box under Physical Model.
For either method, also select the Use radiation groups check box to enable the ability
to define radiation groups, which can, in many cases, speed up the radiation
calculations.
ADVANCED SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Normally these settings do not need to be changed.
DISCRETIZATION
To display this section, click the Show button ( ) and select Discretization.
Select an Element typeMixed order (the default) or Equal order. Mixed order means
that the Thermal Stress interface uses shape functions that are one order higher for the
displacements than for the temperature. Select U1+T1, U2+T1 (the default), U3+T2, or
U4+T3 for the Thermal stress fields for mixed-order elements or the corresponding
element-order combinations for equal-order elements. U2+T1, for example, means
second-order elements for the displacements and first-order elements for the
temperature.
For the thermal radiation, if applicable, select Linear (the default), Quadratic, Cubic,
Quartic, or (in 2D) Quintic for the Surface radiosity. Specify the Value type when using
splitting of complex variablesReal (the default) or Complex.
Initial Values
Thermal Hyperelastic Material
Thermal Linear Elastic Material
Thermal Linear Viscoelastic Material
These features are described for the Heat Transfer interface in the COMSOL
Multiphysics Users Guide (listed in alphabetical order):
These features are described for Solid Mechanics interface (listed in alphabetical order):
Added Mass
Antisymmetry
Body Load
Boundary Load
Contact
Edge Load
Fixed Constraint
Free
Linear Elastic Material
Linear Viscoelastic Material
Periodic Condition
Point Load
Pre-Deformation
Prescribed Acceleration
Prescribed Displacement
Prescribed Velocity
Rigid Connector
Roller
Spring Foundation
Symmetry
Thin Elastic Layer
Initial Values
The Initial Values feature adds initial values for the displacement field (the displacement
components u, v, and w in 3D), the temperature, and the surface radiosity (applicable
for surface-to-surface radiation only) that can serve as an initial condition for a
DOMAIN SELECTION
From the Selection list, choose the domains to define initial values.
IN IT IA L VA LUES
Enter values or expressions for the initial values of the Displacement field u (default
value 0 m) and of the Temperature T (default value 293.15 K).
DOMAIN SELECTION
From the Selection list, choose the domains to define initial values.
See the Linear Elastic Material node for details about this section and as
described for the Solid Mechanics interface.
Note
THERMAL EXPANSION
By default, the Coefficient of thermal expansion (SI unit: 1/K) is taken From material.
Select User defined to enter a value or expression for , then select Isotropic, Diagonal,
Symmetric, or Anisotropic and enter one or more components for a general thermal
expansion coefficient vector vec.
HEAT CONDUCTION
The default Thermal conductivity (SI unit: W/(mK)) uses values From material. If User
defined is selected, choose Isotropic to define a scalar value, or Diagonal, Symmetric, or
Anisotropic to enter other values or expressions in the field or matrix.
THERMODYNAMICS
The default Heat capacity at constant pressure Cp (SI unit: J/(kgK)) is taken From
material. If User defined is selected, enter a value.
GEOMETRIC NONLINEARITY
With the addition of the Heat Transfer Module or the CFD Module,
right-click the Thermal Linear Elastic Material node to add a Pressure Work
Note node.
With the Heat Transfer Module or the CFD Module, right-click the
Thermal Hyperelastic Material node to add a Pressure Work node.
Note
HYPERELASTIC MATERIAL
Select a hyperelastic Material modelNeo-Hookean, Mooney-Rivlin, Murnaghan, or User
defined. Then see below for details about each selection.
For any selection, the material Density (SI unit: kg/m3) uses values From material by
default. Select User defined to enter a different value or expression.
Neo-Hookean
If Neo-Hookean is selected as the Material model, the default values for both Lam
constant (SI unit: Pa) and Lam constant (SI unit: Pa) use values From material.
Select User defined to enter different values or variables.
Mooney-Rivlin
If Mooney-Rivlin is selected as the Material model, the Model parameters C10 (SI unit: Pa)
and C01 (SI unit: Pa) use values From material. Select User defined to enter different
values or variables. Enter the Initial bulk modulus (SI unit: Pa).
Murnaghan
If Murnaghan is selected as the Material model, the Murnaghan third-order elastic moduli
constants l (SI unit: Pa), m (SI unit: Pa), and n (SI unit: Pa) and the Lam constants
(SI unit: Pa) and (SI unit: Pa) use values From material. Select User defined to enter
different values or variables for the constants as required.
User defined
If User defined is selected as the Material model, enter an expression for the Strain energy
density Ws (SI unit: J/m3).
For the rest of the settings, see Thermal Linear Elastic Material for details
about the Thermal Expansion, Heat Conduction, and Thermodynamics
Note sections.
Right-click the Thermal Linear Viscoelastic Material node to add an Initial Stress and
Strain node. The settings window for this node is the same as for the Solid Mechanics
interface.
L O N G - TE R M E L A S T I C P R O P E R T I E S A N D G E N E R A L I Z E D M A X W E L L M O D E L
THERMAL EFFECTS
Thermal Effects
See Also
For the rest of the settings, see Thermal Linear Elastic Material for details
about the Heat Conduction, Thermodynamics, and Geometric Nonlinearity
Note sections.
When this interface is added, these default nodes are also added to the Model Builder
Fluid Properties, Linear Elastic Material, and Free Deformation (for the mesh movement
and default boundary conditions) in the domains; Wall (for the fluid), Prescribed Mesh
Displacement (for the mesh movement), and Free (for the solid mechanics, which
initially is not applicable to any boundary because the default settings assume a fluid
domain) as default boundary conditions; and Initial Values. In addition, for the
fluid-solid boundary, a Fluid-Solid Interface Boundary node adds the fluid-structure
interaction. This node is only applicable to interior fluid-solid boundaries.
Right-click Fluid-Structure Interaction to add other nodes that implement, for example,
loads, constraints, and nonlinear materials for the solid domain. The Fluid-Structure
Interaction page contains the following sections plus additional sections that are similar
INTERFACE IDENTIFIER
The interface identifier is a text string that can be used to reference the respective
physics interface if appropriate. Such situations could occur when coupling this
interface to another physics interface, or when trying to identify and use variables
defined by this physics interface, which is used to reach the fields and variables in
expressions, for example. It can be changed to any unique string in the Identifier field.
The default identifier (for the first interface in the model) is fsi.
DOMAIN SELECTION
The default setting is to include All domains in the model to define fluid-structure
interaction (domains representing the fluid and the solid). To choose specific domains,
select Manual from the Selection list.
Depending or your license, the settings in the Physical Model section can also include
selections of a turbulence model, Stokes flow, and channel flow approximation.
TY P I C A L W AV E S P E E D F O R PE R F E C T L Y M A T C H E D L A Y E R S
The typical wave speed cref is a parameter for the perfectly matched layers (PMLs) if
used in a solid wave propagation model. The default value is fsi.cp, the
pressure-wave speed. If you want to use another wave speed, enter a value or
expression in the Typical wave speed for perfectly matched layers field.
For modeling using a memory-efficient form, there are also dependent variables for the
corrected velocity field and its components and for the corrected pressure.
DISCRETIZATION
To display this section, click the Show button ( ) and select Discretization. Select a
Discretization of fluidsP1+P1, P2+P1 (the default), or P3+P2. Select a Displacement
fieldLinear, Quadratic (the default), Cubic, or Quartic. Specify the Value type when
using splitting of complex variablesReal or Complex (the default).
The consistent stabilization methods are available for the Navier-Stokes equations.
The Isotropic diffusion inconsistent stabilization method is available for the
Navier-Stokes equations.
These features are unique to the interface and described in this section:
Initial Values
Fluid-Solid Interface Boundary
These features are described for the Solid Mechanics interface (listed in alphabetical
order):
Added Mass
Antisymmetry
Body Load
Boundary Load
Contact
Edge Load
Fixed Constraint
Free
Fixed Mesh
Free Deformation
Prescribed Deformation
Prescribed Mesh Displacement
These features are described for the Laminar Flow interface in the COMSOL
Multiphysics Users Guide (listed in alphabetical order):
Boundary Stress
Flow Continuity
Initial Values
The Initial Values node adds initial values for pressure, turbulent dissipation rate,
turbulent kinetic energy, reciprocal wall distance, displacement field and velocity field.
These variables can serve as an initial condition for a transient simulation or as an initial
guess for a nonlinear analysis. If more than one set of initial values is required,
right-click to add additional Initial Values nodes.
DOMAIN SELECTION
From the Selection list, choose the domains to define an initial value.
INITIAL VALUES
Enter the initial values as values or expressions. The variables for turbulence are only
valid for fluid flow using a turbulence model.
BOUNDARY SELECTION
BOUNDARY SELECTION
From the Selection list, choose the boundaries to define a prescribed mesh
displacement. The default setting is to include all boundaries in the model.
Interior Wall
If you have a CFD Module or Heat Transfer Module, the Interior Wall
boundary condition is also available and is documented in the CFD
Module Users Guide or Heat Transfer Module Users Guide,
respectively. This boundary condition is useful for avoiding meshing of
Tip
thin wall structures by using no-slip conditions on interior curves and
surfaces.
1 In the Model Wizard, select the Fluid-Structure Interaction interface ( ) on the Add
Physics page; then click Next ( ).
2 Select Stationary ( ) or Time Dependent ( ) from the Preset Studies branch on
the Select Study Type list; then click Finish ( ).
3 Create the geometry, which should contain a fluid domain and a solid domain.
4 Add the materials, typically a fluid and a solid. Then assign the material added last
to the domains that represents the solid (or the fluid, if the solid material was added
first).
5 By default, the Fluid-Structure Interaction interface adds a Fluid Properties node for
the fluid domain as well as a Free Deformation node for the mesh displacements in
the moving mesh to all domains in the geometry. For the solid domain, the default
is the Linear Elastic Material node with the setting to include geometric nonlinearity
and an initially empty selection. To use another material model for the solid,
right-click the Fluid-Structure Interaction node and from the Solid Mechanics
submenu, in the top section of the context menu, select a material node for the
solid. Add the domains that represent the solid to its selection. That selection
automatically overrides the Fluid Properties and Free Deformation nodes in the solid.
6 Verify that the default boundary conditions are correctly assigned for the three types
of boundaries in the model: the Wall node for all fluid boundaries (and a Prescribed
Mesh Displacement node for zero mesh displacements on the same boundaries), the
Free node for all solid boundaries, and the Fluid-Solid Interface Boundary node on the
interior boundaries between the fluid and the solid. The Fluid-Solid Interface
Boundary node implements the coupling from the force exerted on the solid
boundary by the fluid as well as the as the structural velocities acting on the fluid as
a moving wall.
The fluid flow is described by the Navier-Stokes equations, which provide a solution
for the velocity field ufluid. The total force exerted on the solid boundary by the fluid
is the negative of the reaction force on the fluid,
f = n pI + u fluid + u fluid T --- u fluid I
2
(8-1)
3
where p denotes pressure, the dynamic viscosity for the fluid, n the outward normal
to the boundary, and I the identity matrix. Because the Navier-Stokes equations are
solved in the spatial (deformed) frame while the solid mechanics interfaces are defined
in the material (undeformed) frame, a transformation of the force is necessary. This is
done according to
dv
F = f --------
dV
where dv and dV are the mesh element scale factors for the spatial frame and the
material (reference) frame, respectively.
u solid
t
(the rate of change for the displacement of the solid), which act as a moving wall for
the fluid domain. The predefined Fluid-Solid Interface Boundary condition includes
these couplings for bidirectionally coupled FSI simulations.
You can use this multiphysics interface for coupled thermal, electrical, and structural
analysis of, for example, the movement of some actuator, where an electric current
causes a temperature increase, which in turn leads to a displacement through thermal
expansion.
In addition, all other features from the individual physics interfaces are accessible.
INTERFACE IDENTIFIER
The interface identifier is a text string that can be used to reference the respective
physics interface if appropriate. Such situations could occur when coupling this
interface to another physics interface, or when trying to identify and use variables
defined by this physics interface, which is used to reach the fields and variables in
expressions, for example. It can be changed to any unique string in the Identifier edit
field.
The default identifier (for the first interface in the model) is tem.
DOMAIN SELECTION
The default setting is to include All domains in the model to define Joule heating and
thermal expansion. To choose specific domains, select Manual from the Selection list.
From the 2D approximation list select Plane stress or Plane strain. When
modeling using plane stress, the interface solves for the out-of-plane strain
2D components in addition to the displacement field u.
THICKNESS
S T R U C T U R A L TR A N S I E N T B E H AV I O R
From the Structural transient behavior list, select Quasi-static or Include inertial terms to
treat the elastic behavior as quasi-static (with no mass effects; that is, no second-order
time derivatives) or as a mechanical wave in a time-dependent study. The default is to
use the quasi-static behavior for a time-dependent study. Select Include inertial terms
to model the structural transient behavior as a mechanical wave.
SWEEP SETTINGS
Select the Activate terminal sweep check box to switch on the sweep and invoke a
parametric sweep over the terminals. Enter a Sweep parameter name to assign a specific
name to the variable that controls the terminal number solved for during the sweep.
The generated lumped parameters are in the form of capacitance matrix elements. The
terminal settings must consistently be of either fixed voltage or fixed charge type.
Select an Output format for the Touchstone exportMagnitude angle, Magnitude (dB)
angle, or Real imaginary. Enter a Reference impedance Zref (SI unit: ). The default is
50 .
PHYSICAL MODEL
If you have the Heat Transfer Module, the out-of-plane heat transfer,
surface-to-surface radiation, and radiation in participating media options
are available in this section and described in the Heat Transfer Module
Tip
Users Guide.
RADIATION SETTINGS
ADVANCED SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Normally these settings do not need to be changed.
DISCRETIZATION
To display this section, click the Show button ( ) and select Discretization. Select
Linear, Quadratic, Cubic, Quartic, or (in 2D) Quintic for the order of the elements for
each of these variablesElectric potential, Displacement field, Temperature, Surface
radiosity, and Radiative intensity. The default is to use quadratic elements for the electric
potential, displacements, and temperature, and to use linear elements for the surface
radiosity. Specify the Value type when using splitting of complex variablesReal or
Complex (the default).
See Also Domain, Boundary, Edge, Point, and Pair Features for the Joule
Heating and Thermal Expansion Interface
Initial Values
The Initial Values node adds initial values for electric potential, displacement field,
temperature, and surface radiosity. These variables can serve as an initial condition for
DOMAIN SELECTION
From the Selection list, choose the domains to define an initial value.
INITIAL VALUES
Enter the initial values as values or expressions for Electric potential V (SI unit: V),
Displacement field u (SI unit: m), and Temperature T (SI unit: K). The default values
are 0 except for the temperature, which has a default initial value of 293.15 K.
Domain, Boundary, Edge, Point, and Pair Features for the Joule
Heating and Thermal Expansion Interface
Because The Joule Heating and Thermal Expansion Interface is a multiphysics
interface, every feature (except Initial Values) is shared with, and described for, other
interfaces. Below are links to the domain, boundary, edge, point, and pair features as
indicated.
These features are described for the Solid Mechanics interface (listed in alphabetical
order):
Added Mass
Antisymmetry
Body Load
Boundary Load
Contact
Edge Load
Fixed Constraint
Free
Linear Elastic Material
These features are described for the Heat Transfer and Joule Heating interfaces in the
COMSOL Multiphysics Users Guide:
Using Symmetries
Continuity on Interior Boundaries
See Also
Identity and Contact Pairs
Specifying Boundary Conditions for Identity Pairs
When this interface is added, these default nodes are also added to the Model Builder
Piezoelectric Material, Free (for the solid mechanics and default boundary conditions),
Zero Charge (for the electric potential), and Initial Values. Right-click the Piezoelectric
Devices node to add other features that implement, for example, loads, constraints, and
solid mechanics and electric materials. In 2D and 2D axial symmetry, adding a
Piezoelectric Devices interface also adds predefined base-vector coordinate systems for
the materials (in the plane 2D case) XY-, YZ-, ZX-, YX-, XZ-, and XY-planes. These
additional coordinate systems are useful for simplifying the material orientation for the
piezoelectric material.
All functionality from the Solid Mechanics and Electric Current interfaces is
accessible for modeling the solid and electric properties and
non-piezoelectric domains. Only the features unique to this interface are
described in this section. For details about the shared features see:
The default identifier (for the first interface in the model) is pzd.
DOMAIN SELECTION
The default setting is to include All domains in the model to define the dependent
variables and the equations. To choose specific domains, select Manual from the
Selection list.
2D APPROXIMATION
From the 2D approximation list select Plane stress or Plane strain (the
default). When modeling using plane stress, the Piezoelectric Devices
interface solves for the out-of-plane strain components in addition to the
2D
displacement field u.
THICKNESS
Enter a value or expression for the Thickness d (SI unit: m). The default
value of 1 m is suitable for plane strain models, where it represents a a
unit-depth slice, for example. For plane stress models, enter the actual
thickness, which should be small compared to the size of the plate for the
2D
plane stress assumption to be valid. In rare cases, use a Change Thickness
node to change thickness in parts of the geometry.
S T R U C T U R A L TR A N S I E N T B E H AV I O R
From the Structural transient behavior list, select Quasi-static or Include inertial terms to
treat the elastic behavior as quasi-static (with no mass effects; that is, no second-order
time derivatives) or as a mechanical wave in a time-dependent study. The default is to
use the quasi-static behavior for a time-dependent study. Select Include inertial terms
to model the structural transient behavior as a mechanical wave.
ADVANCED SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Normally these settings do not need to be changed.
DISCRETIZATION
To display this section, click the Show button ( ) and select Discretization. Select
Linear, Quadratic (the default), Cubic, Quartic, or (in 2D) Quintic for the Displacement
field and Electric potential. Specify the Value type when using splitting of complex
variablesReal or Complex (the default).
These features are described for the Solid Mechanics interface (listed in alphabetical
order):
Added Mass
Antisymmetry
Body Load
Boundary Load
Contact
Edge Load
Fixed Constraint
Free
Linear Elastic Material
Linear Viscoelastic Material
Point Load
Pre-Deformation
Prescribed Acceleration
Prescribed Displacement
Prescribed Velocity
Rigid Connector
Roller
Spring Foundation
Symmetry
Thin Elastic Layer
These features are described for the Electrostatics interface in the COMSOL
Multiphysics Users Guide:
Piezoelectric Material
Use the Piezoelectric Material to define the piezoelectric material properties on
stress-charge form using the elasticity matrix and the coupling matrix or on
strain-charge form using the compliance matrix and the coupling matrix. The default
settings is to use material data defined for the material in the domain. Right-click
Piezoelectric Material to add Electrical Conductivity (Time-Harmonic), Initial Stress and
Strain, and Damping and Loss nodes as required.
For entering these matrices, the ordering is different from the standard
ordering used in COMSOL Multiphysics. Instead, use the following
order (Voigt notation), which is the common convention for piezoelectric
Important
materials: xx, yy, zz, yz, xz, zy.
MODEL INPUTS
This section has field variables that appear as model inputs, if the current settings
include such model inputs. By default, this section is empty.
For Stress-charge form, select an Elasticity matrix (ordering: xx, yy, zz, yz, xz, xy) (cE)
(SI unit: 1/Pa).
For a Strain-charge form, select a Compliance matrix (ordering: xx, yy, zz, yz, xz, xy)
(sE) (SI unit: 1/Pa).
Select a Coupling matrix (ordering: xx, yy, zz, yz, xz, xy) (d) (SI unit: C/m2 or C/N).
Select a Relative permittivity (erS or erT) (unitless).
Select a Density (p) (SI unit: kg/m3).
GEOMETRIC NONLINEARITY
If a study step is geometrically nonlinear, the default behavior is to use a large strain
formulation in all domains. There are however some cases when you would still want
to use a small strain formulation for a certain domain. In those cases, select the Force
linear strains check box. When selected, a small strain formulation is always used,
independently of the setting in the study step.
See Also See The Solid Mechanics Interface for details about this section.
DOMAIN SELECTION
From the Selection list, choose the domains to define.
MODEL INPUTS
This section contains field variables that appear as model inputs, if the current settings
include such model inputs. By default, this section is empty.
ELECTRIC FIELD
Select a Constitutive relationRelative permittivity, Polarization, or Remanent
displacement.
GEOMETRIC NONLINEARITY
If a study step is geometrically nonlinear, the default behavior is to use a large strain
formulation in all domains. There are however some cases when you would still want
to use a small strain formulation for a certain domain. In those cases, select the Force
See Also See The Solid Mechanics Interface for details about this section.
DOMAIN SELECTION
From the Selection list, choose the domains to define.
By default, this node inherits the selection from its parent node, and only
a selection that is a subset of the parent nodes selection can be used.
Note
CONDUCTION CURRENT
Select an Electrical conductivity (SI unit: S/m). Select:
From material to use the conductivity value from the domain material.
DOMAIN SELECTION
From the Selection list, choose the domains to add damping.
By default, this feature inherits the selection from its parent node, and use
a selection that is a subset of the parent nodes selection can be used.
Note
DAMPING SETTINGS
Select a Damping typeRayleigh damping, Loss factor for cE, Loss factor for sE, No
damping, or Isotropic loss factor:
No damping
For Rayleigh damping, enter the Mass damping parameter dM and the Stiffness
damping parameter in the dM corresponding fields. The default values are 0, which
means no damping.
For Loss factor for cE, select From material (the default) from the Loss factor for
elasticity matrix cE list to use the value from the material or select User defined to
enter values or expressions for the loss factor in the associated fields. Select
Symmetric to enter the components of cE in the upper-triangular part of a
For Loss factor for e and Loss factor for d, select a Loss factor for coupling matrix e or d
from the list. Select User defined to enter values or expressions for the loss factor in the
associated fields. Select Symmetric to enter the components of e or d in the
upper-triangular part of a symmetric 6-by-6 matrix or select Isotropic to enter a single
scalar loss factor. The default values are 0.
For Loss factor for S and Loss factor for T, select a Loss factor for permittivity. Select
From material (the default) to use the value from the material or select User defined to
enter values or expressions for the loss factor in the associated fields. Select Symmetric
to enter the components of eS or eT in the upper-triangular part of a symmetric
6-by-6 matrix, select Isotropic to enter a single scalar loss factor, or select Diagonal. The
default values are 0.
By default, this node inherits the selection from its parent node, and only
a selection that is a subset of the parent nodes selection can be used.
Note
Dielectric Loss
Right-click the Electrical Material Model node to add a Dielectric Loss subnode to
include a dielectric loss using a dielectric loss factor.
DOMAIN SELECTION
From the Selection list, choose the domains to define.
By default, this node inherits the selection from its parent node, and only
a selection that is a subset of the parent nodes selection can be used.
Note
Initial Values
The Initial Values feature adds an initial value for the displacement field and the electric
potential. Right-click to add additional Initial Values nodes.
DOMAIN SELECTION
From the Selection list, choose the domains to define. The default setting is to include
all domains in the model.
Periodic Condition
The Periodic Condition feature adds a periodic boundary condition. This periodicity
make uix0uix1 for a displacement component ui or similarly for the electric
potential. Control the direction that the periodic condition applies to and if it applies
to the electric potential. Right-click the Periodic Condition node to add a Destination
Selection boundary condition. If the source and destination boundaries are rotated
with respect to each other, this transformation is automatically performed, so that
corresponding displacement components are connected.
This feature works well for cases like opposing parallel boundaries. In
other cases use a Destination Selection subnode to control the destination.
Note By default it contains the selection that COMSOL Multiphysics identifies.
PERIODICITY SETTINGS
Select a Type of periodicityContinuity (the default), Antiperiodicity, Floquet periodicity,
Cyclic symmetry, or User defined.
If Floquet periodicity is selected, enter a k-vector for Floquet periodicity kF (SI unit:
rad/m) for the X, Y, and Z coordinates (3D models), or the R and Z coordinates
(2D axisymmetric models), or X and Y coordinates (2D models).
If Cyclic symmetry is selected, select a Sector angleAutomatic (the default), or User
defined. If User defined is selected, enter a value for S (SI unit: rad). For any
selection, also enter a Mode number m (unitless).
If User defined is selected, select the Periodic in u, Periodic in v (for 3D and 2D
models), and Periodic in w (for 3D and 2D axisymmetric models) check boxes as
required. For all dimensions the Periodic in V check box is also available. Then for
Periodic Condition
The inverse piezoelectric effect, on the other hand, constitutes the opposite of the
direct effect. This means that an applied potential difference induces a deformation of
the crystal.
PIEZOELECTRICITY CONVENTIONS
The documentation and the Piezoelectric Devices interface use piezoelectricity
conventions as much as possible. These conventions differ from those used in other
structural mechanics interfaces. For instance, the numbering of the shear components
in the stress-strain relation differs, as the following section describes. However, the
names of the stress and strain components remain the same as in the other structural
mechanics interfaces.
STRESS-CHARGE
T
T = cE S e E
D = eS + S E
STRAIN-CHARGE
T
S = sE T + d E
D = dT + T E
The numbering of the strain and stress components is also different in piezoelectricity
theory and structural mechanics theory, and it is quite important to keep track of this
aspect in order to provide material data in the correct order. In structural mechanics
the following is the most common numbering convention, and it is also the one used
in the structural mechanics interfaces:
xx xx xx
yy yy yy
zx zz zz
= = =
xy xy 2 xy
yz yz 2 yz
xz xz 2 xz
The Piezoelectric Devices interface uses the immediately preceding piezo numbering
convention (Voigt notation) to make it easier to work with material data and to avoid
mistakes.
The constitutive relation using COMSOL Multiphysics symbols for the different
constitutive forms are thus:
STRESS-CHARGE
T
= cE e E
D = e + 0 rS E
STRAIN-CHARGE
T
= sE + d E
D = d + 0 rT E
Most material data appears in the strain-charge form, and it can be easily transformed
into the stress-charge form. In COMSOL Multiphysics both constitutive forms can be
used; simply select one, and the software makes any necessary transformations. The
following equations transform strain-charge material data to stress-charge data:
1
cE = sE
1
e = d sE
1 T
S = 0 rT d s E d
The piezoelectric material operates as described above, whereas using the two other
materials, structural and electrical problems can be modeled, together or either of
them independently.
Piezoelectric Dissipation
In order to define dissipation in the piezoelectric material for a time-harmonic analysis,
all material properties in the constitutive relations can be complex-valued matrices
where the imaginary part defines the dissipative function of the material.
X = X 1 j X
See also the same references for an explanation of the sign convention. It is also
possible to define the electrical conductivity of the piezoelectric material: S or T
depending on the constitutive relation. Electrical conductivity does not appear directly
in the constitutive equation, but it appears as an additional term in the variational
formulation (weak equation).
T
= cE 0 e E + 0
D = e 0 + 0 rS E + D 0
When solving the model, these program does not interpret these fields as a constant
initial state, but they operate as additional fields that are continuously evaluated. Thus
use these initial field to add, for example, thermal expansion or pyroelectric effects to
models.
1 u i u j u k u k
ij = --- -------- + -------- + --------- --------- (8-2)
2 X j X i X i X j
The first two items above result in another set of constitutive equations for large
deformation piezoelectricity:
T
S = cE e Em
P m = e + 0 rS 0 I E m
= 0 rS 0 I
Electric displacement field in the material orientation results from the following
relation
1
D m = P m + 0 JC E m
T
C = F F
Fields in the global orientation result from the following transformation rules:
T
E = F Em
1
P = J FP m
(8-3)
1
D = J FD m
1
v = V J
where F is the deformation gradient; J is the determinant of F; and v and V are the
volume charge density in present and material coordinates. The deformation gradient
is defined as the gradient of the present position of a material point xX + u:
The electrical part separates into two different cases: For solid domains the electric
energy consists of polarization energy within the solid and the electric energy of free
space occupied by the deformed solidthe same as for the piezoelectric materials. For
nonsolid domains this separation does not occur, and the electric displacement in these
domains directly results from the electric fieldthe electric displacement relation:
Dm = 0 r Em
The use of ALE has impacts on the formulation of the electrical large deformation
equations. The first impact is that with ALE, the gradient of electric potential directly
results in the electric field in the global orientation, and the material electric field
results after transformation.
The most visible impact is on the boundary conditions. With ALE any surface charge
density or electric displacement is defined per the present deformed boundary area,
whereas for the case without ALE they are defined per the undeformed reference area.
HYSTERETIC LOSS
The equivalent viscous and loss factor damping are special cases of a more general way
of defining damping: hysteretic loss. Generally, and independently of the microscopic
origin of the loss, the dissipative behavior of the material can be modeled using
complex-valued material properties. For the case of piezoelectric materials, this means
that the constitutive equations are written as follows:
T
= c E e E
D = e + E
S
T
= s E + d E
D = d + ET
where c E , d , and are complex-valued matrices, where the imaginary part defines
the dissipative function of the material.
Similarly to the real-valued material data, it is not possible to freely define the
complex-valued data. Instead the data must fulfill certain requirement to represent
physically proper materials. A key requirement is that the dissipation density is positive;
that is, there is no power gain from the passive material. This requirement sets rules for
the relative magnitudes for all material parameters, this is important to know, especially
when defining the coupling losses.
X = real X
X = imag X real X
the ratio of the imaginary and real part, and the complex data is then
where the sign depends on the material property used. The loss factors are specific to
the material property, and thus they are named according to the property they refer to,
for example, cE. For a structural material without coupling, simply use s, the
structural loss factor.
Depending on the field, different terminology is in use. For example, the loss tangent
tan might be referred to when working with electrical applications. The loss tangent
has the same meaning as the loss factor. Often the quality factor Qm is defined for a
material. The quality factor Qm and the loss factor i are inversely related: i1 Qm,
where i is the loss factor for cE, sE, or the structural loss factor depending on the
material.
The Piezoelectric Devices interface uses a formulation that assumes that a positive loss
factor corresponds to a positive loss. The complex-valued data is then based on sign
rules. For piezoelectric materials, the following equations apply (m and n refer to
elements of each matrix):
m n m n m n
c E = c E 1 + j cE
m n m n m n
e =e 1 j e
m n m n m n
S = S 1 j S
(8-4)
m n m n m n
s E = s E 1 j sE
m n m n m n
d =d 1 j d
m n m n m n
T = T 1 j T
The losses for other than piezoelectric materials are more straightforward to define.
Again, using the complex stiffness and permittivity, the following equations describe
the lossy material:
m n m n m n
D = 1 + j D
(8-5)
m n m n m n
e = 1 j e e
Often fully defined complex-valued data is not accessible. In the Piezoelectric Devices
interface the loss factors can be defined as full matrices or as scalar isotropic loss factors
independently of the material and the other coefficients.
T H E L O S S F A C T O R U S I N G D I F F E R E N T D A M P I N G TY P E S
The following damping types use an isotropic loss factor s:
In each case the meaning of the loss factor is the same: the fractional loss of energy per
cycle.
The difference between these damping types is how the loss enters the equation
system. Using the isotropic loss, s is used to build complex-valued material
properties, whereas when using the loss factor damping, s appears in a
complex-valued multiplier in the stress-strain relation. In the equivalent viscous
damping, s appears in a complex-valued and frequency-dependent expression for dK
of the Rayleigh damping model.
Jp
D = r 0 E j ------
(the actual displacement variables do not contain any conductivity effects) or in the
total current expression
J = Jd + Jp
Both a dielectric loss factor (Equation 8-4 and Equation 8-5) and the electrical
conductivity can be defined at the same time. In this case, ensure that the loss factor
refers to the alternating current loss tangent, which dominates at high frequency,
where the effect of ohmic conductivity vanishes (Ref. 7).
5. P.C.Y. Lee, N.H. Liu, and A. Ballato, Thickness Vibrations of a Piezoelectric Plate
with Dissipation, IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency
Control, vol. 51, no. 1, 2004.
6. P.C.Y. Lee and N.H. Liu, Plane Harmonic Waves in an Infinite Piezoelectric Plate
with Dissipation, Frequency Control Symposium and PDA Exhibition, IEEE
International, pp. 162169, 2002.
Materials
This chapter describes the materials databases included with the Structural
Mechanics Module.
In this chapter:
393
Material Library and Databases
The Structural Mechanics Module includes these materials databases: Liquids and
Gases ( ), with temperature-dependent fluid dynamic and thermal properties, MEMS
( ), an extended solid materials library with metals, semiconductors, insulators, and
polymers common in MEMS devices, and a Piezoelectric ( ) database with over 20
common piezoelectric materials.
For detailed information about all the other materials databases and the
separately purchased Material Library, see the section Materials in the
Note COMSOL Multiphysics Users Guide.
In this section:
Material LibraryPurchased
separately. Select from over 2500
predefined materials.
All the material databases (including the Material Library) are accessed from the
Material Browser. These databases are briefly described below.
RECENT MATERIALS
From the Recent Materials folder ( ), select from a list of recently used materials, with
the most recent at the top. This folder is available after the first time a material is added
to a model.
MATERIAL LIBRARY
An optional add-on database, the Material Library ( ), contains data for over 2500
materials and 20,000 property functions.
BUILT-IN
Included with COMSOL Multiphysics, the Built-In database ( ) contains common
solid materials with electrical, structural, and thermal properties.
AC/DC
Included in the AC/DC Module, the AC/DC database ( ) has electric properties for
some magnetic and conductive materials.
MEMS
Included in the MEMS Module and Structural Mechanics Module, the MEMS
database ( ) has properties for MEMS materialsmetals, semiconductors,
insulators, and polymers.
PIEZOELECTRIC
Included in the Acoustics Module, MEMS Module, and Structural Mechanics
Module, the Piezoelectric database ( ) has properties for piezoelectric materials.
PIEZORESISTIVITY
Included in the MEMS Module, the Piezoresistivity database ( ) has properties for
piezoresistive materials, including p-Silicon and n-Silicon materials.
USER-DEFINED LIBRARY
The User-Defined Library folder ( ) is where user-defined materials databases
(libraries) are created. When any new database is created, this also displays in the
Material Browser.
Figure 9-1: Assigning materials to a heat sink model. Air is assigned as the material to
the box surrounding the heat sink, and aluminum to the heat sink itself.
If a geometry consists of a heat sink in a container, Air can be assigned as the material
in the container surrounding the heat sink and Aluminum as the heat sink material itself
(see Figure 9-1). The Conjugate Heat Transfer interface, selected during model set-up,
has a Fluid flow model, defined in the box surrounding the heat sink, and a Heat
Transfer model, defined in both the aluminum heat sink and in the air box. The Heat
Transfer in Solids 1 settings use the material properties associated to the Aluminum
3003-H18 materials node, and the Fluid 1 settings define the flow using the Air material
properties. The other nodes under Conjugate Heat Transfer define the initial and
boundary conditions.
There are also some physics interface properties that by default define a material as the
Domain material (that is, the materials defined on the same domains as the physics
interface). For such material properties, select any other material that is present in the
model, regardless of its selection.
To access a material property throughout the model (across several materials) and
not just in a specific material, use the special material container root.material. For
example, root.material.rho is the density as defined by the materials in each
domain in the geometry. For plotting, you can type the expression material.rho
to create a plot that shows the density of all materials.
To access a material property from a specific material, you need to know the tags for
the material and the property group. Typically, for the first material (Material 1) the
tag is mat1 and most properties reside in the default Basic property group with the
tag def. The variable names appear in the Variable column in the table under Output
properties in the settings window for the property group; for example, Cp for the
When using the Material Browser, the words window and page are
interchangeable. For simplicity, the instructions refer only to the Material
Note Browser.
To clear the search field and browse, delete the search term and click
Search to reload all the databases.
Tip
- Click to open each database and browse for a specific material by class (for
example, in the Material Library) or physics module (for example, MEMS
Materials).
Always review the material properties to confirm they are applicable for
the model. For example, Air provides temperature-dependent properties
Important that are valid at pressures around 1 atm.
GROUP MATERIAL
The MEMS Materials database has data for these materials and a default
appearance for 3D models is applied to each material as indicated.
3D
SPECULAR EXPONENT
SURFACE ROUGHNESS
NORMAL INCIDENCE
NOISE FREQUENCY
REFLECTANCE AT
NORMAL VECTOR
NORMAL VECTOR
COLOR OPACITY
NOISE SCALE
All gases Air Simple 0.08 3 0.1 - - -
All liquids Water Cook-Torrance 0.2 0.2 0.2 - 0.7 0.05
(except Engine
oil, Mercury, and
Transformer oil)
Engine oil Custom Cook-Torrance 0.2 0.2 0.2 - 0.7 0.05
Mercury Custom Cook-Torrance 0 1 1 - 0.9 0.1
Transformer oil Plastic Blinn-Phong 0 1 1 64 - -
6. C.F. Spencer and B.A. Adler, A Critical Review of Equations for Predicting
Saturated Liquid Density, J. Chem. Eng. Data, vol. 23, no. 1, pp. 8288, 1978.
8. R.C.Weast (editor), CRC Handbook of Chemistry and Physics, 69th ed., CRC
Press, 1988.
10. M. Zabransky, V. Ruzicka, Jr., and E.S. Domalski, Heat Capacity of Liquids:
Critical Review and Recommended Values. Supplement I, J. Phys. Chem. Ref. Data,
vol. 30, no. 5, pp. 11991397, 2002.
11. W. Wagner, and H-J Kretzschmar, International Steam Tables, 2nd ed., Springer,
2008.
12. F.P. Incropera and D.P. DeWitt, Fundamentals of Heat and Mass Transfer, Fifth
ed. John Wiley & Sons, 2002.
The basic structure of this library comes from the book Microsensors, MEMS, and
Smart Devices (Ref. 13). The material properties come from two primary sources: the
CRC Handbook of Chemistry and Physics (Ref. 14) and MacMillans Chemical and
Physical Data (Ref. 15). Some of the mechanical properties in the library are instead
more MEMS-specific values from The MEMS Handbook (Ref. 16), and most of the
semiconductor properties are values from Ref. 17. Ref. 18 provides a valuable resource
for cross-checking the insulation material properties.
The table below lists the materials and the corresponding groups:
TABLE 9-3: MEMS MATERIALS
MATERIAL GROUP
MATERIAL GROUP
C [100] Semiconductors
GaAs Semiconductors
Ge Semiconductors
InSb Semiconductors
Si(c) Semiconductors
Poly-Si Semiconductors
Silicon (single-crystal) Semiconductors
Al2O3 Insulators
SiC (6H) Insulators
Si3N4 Insulators
SiO2 Insulators
ZnO Insulators
Borosilicate Insulators
Nylon Polymers
PMMA Polymers
Polymide Polymers
Polyethylene Polymers
PTFE Polymers
PVC Polymers
The MEMS Materials database has data for these materials and a default
appearance for 3D models is applied to each material as indicated.
3D
SPECULAR EXPONENT
SURFACE ROUGHNESS
NORMAL INCIDENCE
NOISE FREQUENCY
REFLECTANCE AT
NORMAL VECTOR
NORMAL VECTOR
COLOR OPACITY
NOISE SCALE
Aluminium (Al) Aluminum Cook-Torrance 0 1 1 - 0.9 0.1
Silver (Ag) Custom Cook-Torrance 0 1 1 - 0.9 0.1
Gold (Au) Gold Cook-Torrance 0 1 1 - 0.9 0.1
Chrome (Cr) Custom Cook-Torrance 0 1 1 - 0.9 0.1
Copper (Cu) Copper Cook-Torrance 0 1 1 - 0.9 0.17
Indium (In) Custom Cook-Torrance 0 1 1 - 0.9 0.1
Titanium (Ti) Titanium Cook-Torrance 0 1 1 - 0.9 0.1
Iron (Fe) Iron Cook-Torrance 0 1 1 - 0.99 0.14
Nickel (Ni) Custom Cook-Torrance 0 1 1 - 0.9 0.1
Lead (Pb) Lead Cook-Torrance 0 1 1 - 0.3 0.1
Palladium (Pd) Custom Cook-Torrance 0 1 1 - 0.9 0.1
Platine (Pt) Custom Cook-Torrance 0 1 1 - 0.9 0.1
Antimon (Sb) Custom Cook-Torrance 0 1 1 - 0.9 0.1
Tungsten (W) Custom Cook-Torrance 0 1 1 - 0.9 0.1
C [100] Custom Blinn-Phong 0 1 1 64 - -
GaAs Custom Cook-Torrance 0 1 1 - 0.9 0.1
Ge Custom Cook-Torrance 0 1 1 - 0.9 0.1
InSb Custom Cook-Torrance 0 1 1 - 0.9 0.1
Si(c) Custom Cook-Torrance 0 1 1 - 0.7 0.5
Poly-Si Custom Cook-Torrance 0 1 1 - 0.7 0.5
Silicon Custom Cook-Torrance 0 1 1 - 0.7 0.5
(single-crystal)
Al2O3 Plastic Blinn-Phong 0 1 1 64 - -
SiC (6H) Plastic Blinn-Phong 0 1 1 64 - -
Si3N4 Plastic Blinn-Phong 0 1 1 64 - -
SiO2 Plastic Blinn-Phong 0 1 1 64 - -
SPECULAR EXPONENT
SURFACE ROUGHNESS
NORMAL INCIDENCE
NOISE FREQUENCY
REFLECTANCE AT
NORMAL VECTOR
NORMAL VECTOR
COLOR OPACITY
NOISE SCALE
ZnO Custom Cook-Torrance 0 1 1 - 0.9 0.1
Borosilicate Plastic Blinn-Phong 0 1 1 64 - -
Nylon Custom Blinn-Phong 0 1 1 500 - -
PMMA Plastic Blinn-Phong 0 1 1 64 - -
Polymide Plastic Blinn-Phong 0 1 1 64 - -
Polyethylene Plastic Blinn-Phong 0 1 1 64 - -
PTFE Plastic Blinn-Phong 0 1 1 64 - -
PVC Plastic Blinn-Phong 0 1 1 64 - -
14. D.R. Lide (editor), CRC Handbook of Chemistry and Physics, 84th edition, CRC
Press, 2003.
15. A.M. James and M.P. Lord, MacMillans Chemical and Physical Data,
MacMillans Press, 1992.
The Piezoelectric materials database included with this module contains the following
materials:
MATERIAL
cE Elasticity matrix
e Coupling matrix, stress-charge
rS Relative permittivity, stress-charge
sE Compliance matrix
d Coupling matrix, strain-charge
rT Relative permittivity, strain-charge
Density
The Piezoelectric Materials database has data for these materials and a
default appearance for 3D models is applied to each material as indicated.
3D
SURFACE ROUGHNESS
SPECULAR EXPONENT
NORMAL INCIDENCE
NOISE FREQUENCY
REFLECTANCE AT
NORMAL VECTOR
NORMAL VECTOR
COLOR OPACITY
NOISE SCALE
Glossary
413
Glossary of Terms
anisotropy Variation of material properties with direction. Both global and local user
defined coordinate systems can be used to define anisotropic material properties.
axial symmetry Symmetry in both load and geometry, solves for the radial (r) and
axial (z) displacement.
bar A line element that only has translational degrees of freedom, capable of
sustaining axial forces, with no bending moments, torsional moments, or shear forces.
Can be used on lines in 2D and 3D.
beam A line element having both translational and rotational degrees of freedom.
Capable of sustaining axial forces, bending moments, torsional moments, and shear
forces. Can be used on curves in 2D and 3D.
cable A tension-only truss member used to model large deformation including sag.
Cauchy stress The most fundamental stress measure defined as force/deformed area
in fixed directions not following the body.
compliance matrix The inverse of the elasticity matrix. See elasticity matrix.
contact model The mathematical method to model bodies that come into contact
with each other.
eigenfrequency study Solving for the undamped natural frequencies and vibration
modes of a structure.
= D
equilibrium equation The equation expressing the equilibrium formulated in the stress
components.
fatigue A term describing the phenomena where a component fails after repeated
loadings and unloadings.
flexibility matrix The inverse of the elasticity matrix. See elasticity matrix.
free vibration The undamped vibration of a structure after it is displaced from the
equilibrium position and released. See also eigenfrequency analysis.
G L O S S A R Y O F TE R M S | 415
frequency response A harmonic analysis solving for the steady-state response from a
harmonic excitation. Typically a frequency sweep is performed, solving for many
excitation frequencies at one time.
hyperelastic material Material where the stresses are computed from a strain energy
density function. Often used to model rubber, but also used in acousto-elasticity.
isotropic hardening A hardening model for an elasto-plastic material where the yield
surface increases in size but maintains its original shape.
large deformation The deformations are so large so the nonlinear effect of the change
in geometry or stress stiffening need to be accounted for.
linear buckling analysis Solves for the linear buckling load using the eigenvalue solver.
mass damping parameter Rayleigh damping parameter, the coefficient in front of the
mass matrix.
mixed formulation A formulation where the pressure have been added as a dependent
variable, used for nearly incompressible materials to avoid numerical problems.
parametric study A study that finds the solution dependence due to the variation of a
specific parameter.
pinned A constraint condition where the displacement degrees of freedom are fixed
but the rotational degrees of freedom are free, typically used for frames modeled using
beams and truss elements.
plane strain An assumption on the strain field where all out-of-plane strain
components are assumed to be zero.
plane stress An assumption on the stress field, all out-of-plane stress components are
assumed to be zero.
primary creep The initial creep stage where the strain rate is decreasing with time.
principle of virtual work States that the variation in internal strain energy is equal to
the work done by external forces.
quasi-static transient study The loads vary slowly so inertia terms can be neglected. A
transient thermal study coupled with a structural analysis can often be treated as
quasi-static.
Rayleigh damping A viscous damping model where the damping is proportional to the
mass and stiffness through the mass and stiffness damping parameters.
secondary creep A creep regime where the strain rate is almost constant.
G L O S S A R Y O F TE R M S | 417
second Piola-Kirchhoff stress Conjugate stress to Green-Lagrange strain used in
large deformation analysis.
shell elements A thin element where both bending and membrane effects are
included.
source boundary One side of a contact pair; the destination boundary is prohibited
to penetrate the source boundary.
stationary study A study where the loads and constraints are constant in time. Also
called static.
stress Internal forces in the material, normal stresses are defined as forces/area normal
to a plane, and shear stresses are defined as forces/area in the plane. A fundamental
concept in structural mechanics.
tertiary creep The creep stage where the strain rate increases very rapidly, followed by
eventual failure.
Tresca stress An effective stress measure that is equal to the maximum shear stress.
INDEX| 419
boundary loads theory 180 converse piezoelectric effect 81
box sections, beams 263 coordinate systems
buckling 16 constraints and 44
built-in materials database 395 in physics symbols 108
bulk modulus loads and 39
elastic moduli 165 local edge system 234
PMLs and 99
C cable elements 318
solid mechanics theory 162
cables 319
coordinate systems, membranes 338
calculating stress and strain 249
Coulomb friction 143
canonical systems 162
critical load factor 62
Cartesian coordinate system and PMLs
cross section (node) 263
99
cross section data (node) 309
Cauchy stress 173
crystal cut standards 82
Cauchy stress tensor 59
curvature 286
Cauchy-Green tensor 171
cyclic symmetry, settings 133
Cauchy-Green tensors 189
cyclic symmetry, theory 210
centrifugal acceleration loads 43
cylindrical coordinate systems 99, 106
change thickness (node)
shell and plate interfaces 224 D damped eigenfrequency study 79
solid mechanics interface 119 damping
circular sections, beams 263 equation of motion and 73
coefficient of thermal expansion linear viscoelastic material and 77
beams 259, 291 loss factors 121
shells and plates 221 losses and 73
thermal couplings and 24 matrix 22
cohesion sliding resistance 197 piezoelectric devices 377, 389
common sections, beam interface 263 point mass 281
complex modulus 79 solid mechanics 120
consistent stabilization settings 28 types 391
constitutive relation, membranes 337 viscoelastic materials and 193
constraint settings 28 damping (node)
constraints 44 beam interface 261
contact (node) 141 shell and plate interfaces 223
contact formulation 66 solid mechanics interface 120
contact help variables 197 damping and loss (node) 377
contact modeling 195 damping models 186187
contact modeling, friction 198 defining
contact pairs 6364 anisotropic materials 169
420 | I N D E X
constraints 44 membrane interface 328
contact models 63 piezoelectric devices interface 371
isotropic materials 164 solid mechanics interface 113
multiphysics models 52 thermal stress interface 345
orthotropic materials 167 truss interface 306
thermoelastic materials 169 edge load (node)
deformation gradient 171 beam interface 268
deformations, modeling 189 membrane interface 332
destinations and sources 64 shell and plate interfaces 228
dielectric loss 391 solid mechanics interface 126
dielectric loss (node) 379 truss interface 310
dielectric loss factor 91 eigenfrequency study 21, 69
direct piezoelectric effect 81, 382 solid mechanics 185
discretization settings 27 eigenvalue solvers 69
displacement field, defining 58 elastic material properties 116118, 330
displacement gradients 163 elastic moduli 164
dissipation, piezoelectric materials 385 elasticity matrix 164
distributed loads, theory 179 elcontact variable 196
documentation, finding 29 electrical conductivity (time-harmonic)
domain features (node) 376
fluid-structure interaction interface electrical material model (node) 375
354 emailing COMSOL 30
joule heating and thermal expansion equation of motion, damping and 73
interface 366 equation view 27
piezoelectric devices interface 371 equivalent viscous damping 78
solid mechanics interface 113 equivalent Youngs modulus 143
thermal stress interface 345 Eulerian frame 59
domain material 398 excitation frequency 22, 78
Duhamel-Hookes law 163 expanding sections 27
Dulong-Petit law 169 explicit damping 80
dynamic cyclic symmetry 210 external loads, shell and plate interfaces
dynamic frictional coefficients 198 249
INDEX| 421
fluid-structure interaction interface 351 hysteretic loss 88, 389
setting up a model 358
I I.R.E. standard, for material orientation
theory 360
82
fold lines 246
IEEE standard, for material orientation
fold-line limit angle 216
82
follower loads 181, 338
IEEE standard, piezoelectric materials
free (node) 127
383
frequency domain study 22
imperfection sensitivity 62
solid mechanics 184
implementation
frequency response study 7879
beams 286
friction (node) 143
trusses 321
friction forces 144
implementing
friction in contact modeling 198
PMLs 98
friction models 143
incompressible deformation state 189
friction traction penalty factor 197
inconsistent stabilization settings 28
G gap distance variable 196 inertial effects, contact modeling 66
generalized Maxwell model 191 initial loads and strains, beams 285
geometric nonlinearity 55, 172 initial stress and strain 386
membranes and 335 beams 285
micromechanics and 57, 60 theory 200
piezoelectric devices 386 trusses 320
solid mechanics theory 205 initial stress and strain (node)
geometry, working with 28 beam interface 259
glass transition temperature, viscoelastic membrane interface 330
materials 193 shell and plate interfaces 222
global coordinate systems 162 solid mechanics interface 146
GMG preconditioners 96 truss interface 308
GMRES iterative solvers 96 initial values (node)
gradient displacements 163 beam interface 262
Green strains 171 fluid-structure interaction interface
Green-Lagrange strain 58 356
Green-Lagrange strain tensors 189 joule heating and thermal expansion
Green-Lagrange strains 171 interface 365
membrane interface 331
H harmonic loads 184
piezoelectric interface 379
heat dissipation 79
shell and plate interfaces 225
Hermitian matrices 96
solid mechanics interface 123
hide button 27
thermal stress interface 346
H-profile sections, beams 263
422 | I N D E X
truss interface 309 linearized buckling analysis 61
initial values, theory for shells and plates liquids and gases materials 401
246 load cases 40
interior wall (node) 357 load multiplier 61
Internet resources 28 loads
inverse piezoelectric effect 382 acceleration 43
isotropic materials applied to beams 280
defining 164 pressure 42
elastic properties 116, 330 singular 41
loss damping and 188 solid mechanics theory 179
loss factor damping and 121 total 43
iterative solvers 96 trusses 310
local coordinate systems 162
J Joule heating and thermal expansion in-
local edge system 234
terface 53, 362
logarithmic decrement 74
K kinematic constraints 46
loss factor damping
knowledge base, COMSOL 30
modeling 79
L Lagrange shape functions, trusses and solid mechanics and 121
319 solid mechanics theory 187
Lagrangian formulations 162 springs and 93
Lagrangian frame 58 loss modulus 79, 194
Lam constants 165 loss tangents 390
large deformations 16, 171 losses and damping 73
piezoelectric materials 386 low-reflecting boundary (node) 158
lattice trusses 318 low-reflecting boundary, theory 209
linear buckling study 16, 23, 57
M mass and moment of inertia (node) 140
linear elastic material (node) 115
mass damping parameter 77
beam interface 257
mass matrix 22
membrane interface 329
mass moment of inertia 280
shell and plate interfaces 219
Material Browser
truss interface 307
opening 399
linear elastic materials 50, 163
material coordinates 161
linear viscoelastic material (node) 133
material frame 58
linear viscoelastic material, damping and
Material Library 395
77
material models 50
linear viscoelastic materials
materials
theory 190
databases 395
linear viscoelasticity 190
domain, default 398
INDEX| 423
linear elastic 163 thermal stress interface 341
linear viscoelastic 190 truss interface 304
liquids and gases 401 viscoelastic material 134
MEMS 405 modeling fluid-structure interaction 358
nearly incompressible 207 modeling, large deformations 189
piezoelectric 383, 409 moment computations 47
piezoelectric devices 385 moments
properties, evaluating and plotting 398 beams 265, 285
membrane interface 326 shells and plates 249
theory 335 solid mechanics and 42
MEMS materials 405 moments of inertia 282
Mindlin-Reissner type shell 241 MPH-files 29
MITC shell formulation 241, 245 multiphysics modeling 52
mixed formulations 51, 208 MUMPS direct solvers 96
Model Builder settings 27
N Navier-Stokes equations 360
Model Library 29
nearly incompressible materials 51, 200,
Model Library examples
207
beam interface 254
no rotation (node)
cross section data 263
beam interface 277
damping 120
shell and plate interfaces 231
fluid-structure interaction interface
nominal stress 205
352
nonlinear geometry 172, 205
geometric nonlinearity 57
normal forces 285
initial stress and strain 146
normal stress 173
Joule heating and thermal expansion
O orientation, piezoelectric material 82
interface 362
orthotropic materials
linear elastic material 257
defining 167
load cases 40
elastic properties 117
membrane interface 326
loss damping and 188
MEMS materials database 405
loss factor damping and 121
model mass 70, 72
override and contribution settings 27
piezoelectric devices interface 369
piezoelectric materials database 409 P pair conditions
prestressed analysis study 24 beam interface 256
rigid connector 137 fluid-structure interaction interface
shell interface 214 354
solid mechanics interface 111 joule heating and thermal expansion
thermal expansion 145 interface 366
424 | I N D E X
membrane interface 328 truss interface 311
piezoelectric devices interface 371 pipe sections, beams 263
solid mechanics interface 113 plane stress and strain 105, 112
thermal stress interface 345 planes, symmetry and constraints 45
truss interface 306 plate interface 214
pair selection 28 theory 241
parametric analysis 16 plates
parametric solvers 40 damping 223
penalized friction traction 197 external loads 249
penalty factors initial and prescribed values 246
contact node and 142 intial stresses and strains 222
contact pairs and 65 linear elastic material 219
theory 197 MITC shell formulation 245
perfectly matched layers. see PML. prescribed acceleration 235
periodic condition (node) prescribed displacement/rotation 231
piezoelectric devices interface 380 prescribed velocity 233
solid mechanics interface 132 stress and strain calculations 249
periodic conditions, theory 210 symmetry and antisymmetry 247
phase (node) thermal expansion 221
beam interface 270 plotting, material properties 398
shell and plate interfaces 229 PML 98, 100
solid mechanics interface 147 point conditions
truss interface 311 beam interface 256
physics settings windows 27 fluid-structure interaction interface
physics symbols 354
coordinate directions for 108 joule heating and thermal expansion
showing 106 interface 366
piezoelectric 82 membrane interface 328
piezoelectric coupling 369 piezoelectric devices interface 371
piezoelectric crystal cut 82 solid mechanics interface 113
piezoelectric devices interface 369 thermal stress interface 345
theory 382 truss interface 306
piezoelectric losses 88 point load (node) 126
piezoelectric material (node) 373 beam interface 269
piezoelectric materials 409 shell and plate interfaces 229
pinned (node) point loads example 41
beam interface 276 point mass (node)
shell and plate interfaces 230 beam interface 280
INDEX| 425
truss interface 316 Rayleigh damping 75, 77, 186
point mass damping (node) reaction forces, evaluating 47
beam interface 281 rectangle sections, beams 263
truss interface 317 reference coordinates 161
Poissons ratio 51, 165 reference point for moment computa-
pre-deformation (node) 154 tion 47
prescribed acceleration (node) refpnt variable 112
beam interface 274 remanent electric displacement (node)
shell and plate interfaces 235 378
solid mechanics interface 149 resonant frequency 75
truss interface 314 results evaluation, for shells 240
prescribed displacement 46 rigid connector 198
prescribed displacement (node) rigid connector (node) 137, 238
solid mechanics interface 128, 333 rigid connectors
truss interface 312 kinematic constraints and 46
prescribed displacement/rotation (node) moments and 42
beam interface 271 rigid domain (node) 140
shell and plate interfaces 231 roller (node) 131
prescribed mesh displacement (node) rotated coordinate system 86
357 rotational degrees of freedom 283
prescribed values, shells and plates 246 rotational joints, beams 46
prescribed velocity (node)
S sagging cables 318
beam interface 273
scaling of eigenvectors 69
shell and plate interfaces 233
second Piola Kirchoff stress 59
solid mechanics interface 148
second Piola-Kirchhoff stress 174
truss interface 313
section orientation (node) 267
pressure loads 42, 181
selecting
pressure-wave speeds 165
solvers 95
prestressed analysis, eigenfrequency
shear modulus expression 165
study 57
shear strain 170
prestressed analysis, frequency domain
shear stress 173
study 57
shear-wave speeds 165
principle of virtual work 182
shell interface 214
Prony series 191
theory 241
Q quality factors and losses 390 shells
quaternion representation, of rigid con- damping 223
nector 199 external loads 249
initial and prescribed values 246
R rate of strain tensor 173
426 | I N D E X
intial stresses and strains 222 St. Venants principle 41
linear elastic material 219 stabilization settings 28
MITC shell formulation 245 static frictional coefficients 198
prescribed acceleration 235 stationary solvers 40, 95
prescribed displacement/rotation 231 stationary study 21
prescribed velocity 233 stiffness damping parameter 77
stress and strain calculations 249 stiffness matrix 22
symmetry and antisymmetry 247 storage modulus 79, 194
thermal expansion 221 straight edge constraint (node) 311
show button 27 straight edges 322
showing strain 170, 173
physics symbols 106 axial symmetry 171
singular loads 41 engineering form 170
skew-symmetric part 172 shear 170
solid mechanics tensor form 170
damping 120 strain tensor 170
edge loads 126 strain-displacement relation 170
intial stresses and strains 146 large displacement 171
prescribed acceleration 149 small displacement 170
prescribed velocity 148 strain-displacement, trusses 319
solid mechanics interface 111 strain-displacement/rotation 284
theory 160 strains, membranes 336
solver methods, augmented Lagrangian stress 173
196 Cachy 173
solver parameters 95 first Piola-Kirchhoff 174
solver settings 95 normal 173
SOR line solvers 97 second Piola-Kirchhoff 174
sources and destinations 64 shear 173
spatial coordinates 161 tensor 173
spatial frame 59 stress and strain tensors 163
spatial stress tensor 59 stress and strain, piezoelectric devices
spherical coordinate systems 99 383
spin tensor 173 stress evaluation, beams 289
SPOOLES solvers 53 stress stiffening 335
spring constant 93 stresses, membranes 336
spring foundation (node) 150 stress-strain relation 173
spring foundation, solid mechanics 93 beams 285
spring foundation, theory 201 trusses 320
INDEX| 427
structural damping 79 beams 290
Structural Mechanics Module 14 trusses 321
study steps, geometric nonlinearity and thermal effects (node) 135
55 thermal expansion
study types 21 loads and 43
eigenfrequency 185 thermal expansion (node)
frequency domain, solid mechanics in- beam interface 258
terface 184 shell and plate interfaces 221
parametric 16 solid mechanics interface 144
stationary, solid mechanics interface truss interface 308
183 thermal expansion, Joule heating and 362
viscoelastic transient initialization 183 thermal hyperelastic material (node) 348
surface traction and reaction forces 49 thermal linear elastic material (node) 347
symbols for physics 106 thermal linear viscoelastic material
symmetric matrices 95 (node) 350
symmetry (node) thermal strain
beam interface 277 beams 285
shell and plate interfaces 236 study types and 24
solid mechanics interface 130 thermal stress interface 341
truss interface 315 modeling 52
symmetry constraints 44 thermal-electric-structural interaction
53
T tangential strains 336
thermal-structural analysis 24
technical support, COMSOL 30
thermoelastic materials, defining 169
temperature loads 43
thermorheologically simple viscoelastic
tensors
materials 192
linear elastic materials 163
thin elastic layer (node) 155
the moving mesh interface, piezoelectric
thin elastic layer, solid mechanics 93
devices and 388
thin elastic layer, theory 202
theory
time-dependent study 21
beam interface 282
torsional constants and moments 288
fluid-structure interaction interface
total force loads 180
360
total Lagrangian formulation 16
membrane interface 335
total loads 43
piezoelectric devices interface 382
t-profile sections, beams 263
shell and plate interfaces 241
tractions 59
solid mechanics interface 160
transient study 21
truss interface 318
TRS material 193
thermal coupling
428 | I N D E X
true stress tensor 59 viscous damping 80, 94
truss interface 304 Voigt form 82
theory 318 Voigt notation 122, 167, 383
trusses volume ratio 171
intial stresses and strains 308, 320
W wave speeds 117, 258, 330
linear elastic material 307
weak constraint settings 28
loads 310
weak constraints, using 48
prescribed displacements 312
web sites, COMSOL 30
prescribed velocity 313
WLF shift functions 136137, 193
straight edge 322
Y Youngs modulus expression 165
strain-displacement 319
thermal coupling 321
thermal expansion 308
typographical conventions 30
U undamped models 91
units, loads and 39
u-profile sections, beams 263
user community, COMSOL 30
using
coordinate systems 162
predefined variables 47
spatial and material coordinates 161
weak constraints 48
uspring variable 93
V variables
elcontact 196
for material properties 398
material and spatial coordinates 161
predefined 47
refpnt 112
vdamper variable 94
viscoelastic materials
axial symmetry 194
definition 51
frequency domain analysis and damp-
ing 193
temperature effects 192
viscoelastic transient initialization 24, 183
INDEX| 429
430 | I N D E X