The Solid Mechanics Interface
The Solid Mechanics (solid) interface (), found under the Structural Mechanics branch () when adding a physics interface, is intended for general structural analysis of 3D, 2D, 1D, or axisymmetric bodies. In 2D, 1D, and 1D axisymmetry, plane stress, plane strain, or generalized plane strain assumptions can be used. The Solid Mechanics interface is based on solving the equations of motion together with a constitutive model for a solid material. Results such as displacements, stresses, and strains are computed.
The functionality provided by the Solid Mechanics interface depends on the products you are using. The Acoustics Module, MEMS Module, and Structural Mechanics Module add several features, for example geometric nonlinearity and advanced boundary conditions such as contact, follower loads, and nonreflecting boundaries.
The default material is a Linear Elastic Material. With either the Nonlinear Structural Materials Module or the Geomechanics Module, the physics interface is extended with more materials, for example, material models for plasticity, hyperelasticity, creep, and concrete. You can also add your own material models using an External Stress-Strain Relation
When this physics interface is added, the following 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. For axisymmetric models, an Axial Symmetry node is also added.
Then, from the Physics toolbar, you can add other nodes that implement, for example, solid mechanics material models, boundary conditions, and loads. You can also right-click Solid Mechanics to select physics features from the context menu.
Settings
The Label is the default physics interface name.
The Name is used primarily as a scope prefix for variables defined by the physics interface. Refer to such physics interface variables in expressions using the pattern <name>.<variable_name>. In order to distinguish between variables belonging to different physics interfaces, the name string must be unique. Only letters, numbers, and underscores (_) are permitted in the Name field. The first character must be a letter.
The default Name (for the first physics interface in the model) is solid.
2D Approximation
From the 2D approximation list, select Plane strain, Plane stress, or Generalized plane strain.
Select Out-of-plane mode extension (time-harmonic) to prescribe an out-of-plane wave number to be used in mode analysis, eigenfrequency, and frequency domain studies. When selected, enter the Out-of-plane wave number kz. The input value will only be taken into account in eigenfrequency and frequency domain studies. For mode analysis, the out-of-plane wave number is computed as an eigenvalue.
For more information, see Out-of-Plane Waves in the Structural Mechanics Theory chapter.
Axial Symmetry Approximation
Select Include circumferential displacement to add a dependent variable to account for the out-of-plane displacement. This can be used to model axial deformation with rotation around the axis of symmetry.
Select Circumferential mode extension to prescribe a circumferential wave number to be used in eigenfrequency or frequency domain studies. When selected, enter the Azimuthal mode number m.
For more information, see Circumferential Displacement and Out-of-Plane Waves in the Structural Mechanics Theory chapter.
Eigenfrequency Analysis of a Free Cylinder: Application Library path Structural_Mechanics_Module/Verification_Examples/free_cylinder.
1D Approximation
From the y direction and z direction lists, select Plane strain, Plane stress, or Generalized plane strain.
1D Axisymmetric Approximation
From the z direction list, select Plane strain, Plane stress, or Generalized plane strain.
Thickness
For 2D and 1D axisymmetric components, enter a value or expression for the thickness d. The default value of 1 m is suitable for plane strain models, where it represents a unit-thickness slice, for example. For plane stress models, consider entering the actual thickness.
Use a Change Thickness node to change thickness in parts of the geometry if necessary.
Cross-Sectional Area
For 1D components, enter a value or expression for the cross-section area Ac. The default value is 1 m2, which is suitable for models where plane strain assumptions are used in both transverse directions. For plane stress models, consider entering actual cross-section data.
Use a Change Cross Section node to change thickness in parts of the geometry if necessary.
Structural Transient Behavior
From the Structural transient behavior list, select Include inertial terms or Quasistatic. Use Quasistatic to treat the dynamic behavior as quasistatic (with no mass effects; that is, no second-order time derivatives). Selecting this option gives a more efficient solution for problems where the variation in time is slow when compared to the natural frequencies of the system. The default solver for the time stepping is changed from Generalized alpha to BDF when Quasistatic is selected.
For problems with creep, and sometimes viscoelasticity, the problem can be considered as quasistatic. This is also the case when the time dependence exists only in some other physics, like a transient heat transfer problem causing thermal strains.
Transient Solver Settings
In this section, you can add instructions used when generating the solver sequence for a wave-propagation transient problem. Select the Maximum frequency to resolve check box to activate this functionality.
Enter the Maximum frequency to resolve in the model, fmax,sol.
Select the Time stepping (method) as Fixed (preferred) or Free. The Free option is in general not recommended for wave problems.
The generated solver will be adequate in most situations if the computational mesh also resolves the frequency content in the model. Note that any changes made to these settings (after the model is solved the first time) will only be reflected in the solver if Show Default Solver or Reset Solver to Defaults is selected in the study.
For highly nonlinear problems with user-defined terms, manual tuning of the solver may be necessary. In nonlinear models, the maximum frequency to resolve should be selected based on the number of harmonics to be resolved.
Typical Wave Speed for Perfectly Matched Layers
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. To use another wave speed, enter a value or expression in the Typical wave speed for perfectly matched layers field.
This section is only available with COMSOL products that include PMLs (see https://www.comsol.com/products/specifications/).
Port Sweep Settings
Select to enable the Activate port sweep option. This option is used to compute the full scattering matrix when Port conditions are used. For more details see The Port Sweep Functionality subsection. The section only exists for 3D geometries.
Automated Model Setup
This section will only be displayed if a mesh on NASTRAN® format, containing RBE2 elements, has been imported in an Import node under Mesh. The purpose is to automatically create rigid connectors from RBE2 elements in the NASTRAN file.
An RBE2 element represents a rigid connection between a set of mesh nodes. This means that it can, and often does, connect elements from different physics interfaces.
In the menu in the section title, you can select Create Rigid Connectors from RBE2. The effect is that one rigid connector will be created for each RBE2 element in the imported file. This will happen for all physics interfaces in the Interfaces list. Supported interfaces are: Solid Mechanics, Shell, Beam, and Multibody Dynamics. If there are RBE2 elements spanning more than one physics interface, they will be automatically connected.
The created rigid connectors will have point, edge, and boundary selections as inferred from the nodes in the RBE2 element and the mesh connectivity. The ‘independent node’ of the RBE2 element is used as center of rotation for the rigid connector.
The Automated Model Setup section is present in the Solid Mechanics, Shell, and Beam interfaces. In a model that contains several physics interfaces, you should use the automated model setup from only one of them, and make sure that all the involved interfaces are selected in the Interfaces list.
Advanced Settings
To display this section, click the Show More Options button () and select Advanced Physics Options in the Show More Options dialog box. Normally these settings do not need to be changed.
You can choose how extra ODE variables added by some features are grouped in the Dependent Variables node of a generated solver sequence.
Select the Rigid materials check box to group variables added by Rigid Material nodes.
Select the Rigid connectors check box to group variables added by Rigid Connector nodes.
Select the Attachments check box to group variables added by Attachment nodes.
The selection made in the Advanced Settings section can be overridden by the settings in the Advanced section of the Rigid Material, Rigid Connector or Attachment features.
Discretization
In the Solid Mechanics interface, you can choose not only the order of the discretization, but also the type of shape functions: Lagrange or serendipity. For highly distorted elements, Lagrange shape functions provide better accuracy than serendipity shape functions of the same order. The serendipity shape functions will however give significant reductions of the model size for a given mesh containing hexahedral, prism, or quadrilateral elements. In 1D components, there is no difference between Lagrange and serendipity shape functions.
The default is to use Quadratic serendipity shape functions for the Displacement field. Using Linear shape functions will give what is sometimes called constant stress elements. Such a formulation will for many problems make the model overly stiff, and many elements may be needed for an accurate resolution of the stresses.
To display other settings for this section, click the Show More Options button () and select Advanced Physics Options in the Show More Options dialog box.
See Table 2-4 for links to common sections and Table 2-5 for common feature nodes. You can also search for information: press F1 to open the Help window or Ctrl+F1 to open the Documentation window.
Dependent Variables
The physics interface uses the global spatial components of the Displacement field u as dependent variables. The default names for the components are (u, v, w) in 3D. In 2D the component names are (u, v), and in 2D axisymmetry they are (u, w). In 1D and 1D axisymmetry the default component name is (u). You can however not use the ‘missing’ component names in 2D or 1D as a parameter or variable name, since they are used internally.
You can change both the field name and the individual component names. If a new field name coincides with the name of another displacement field, the two fields (and the interfaces that define them) share degrees of freedom and dependent variable component names. You can use this behavior to connect a Solid Mechanics interface to a Shell directly attached to the boundaries of the solid domain or to another Solid Mechanics interface sharing a common boundary.
A new field name must not coincide with the name of a field of another type (that is, it must contain a displacement field), or with a component name belonging to some other field. Component names must be unique within a model except when two interfaces share a common field name.
In the COMSOL Multiphysics Reference Manual see Table 2-4 for links to common sections and Table 2-5 to common feature nodes. You can also search for information: press F1 to open the Help window or Ctrl+F1 to open the Documentation window.
Stresses in a Pulley: Application Library path COMSOL_Multiphysics/Structural_Mechanics/stresses_in_pulley
Eigenvalue Analysis of a Crankshaft: Application Library path COMSOL_Multiphysics/Structural_Mechanics/crankshaft