The Multibody Dynamics (mbd) interface (
), found under the
Structural Mechanics branch (
) when adding a physics interface, is intended for analysis of mechanical assemblies. The parts in the assembly can be rigid or flexible, and are connected by various types of joints, gears, chains, springs, or dampers. Flexible parts can be defined using solid, shell, or beam elements. There are many types of joints, such as hinges or ball joints, which can be used based on the type of connection required between components. There are many types of gears, such as spur, helical, bevel, worm, or rack, which can be used for power transmission. Computed results are displacements, velocities, accelerations, joint forces, gear contact forces, and — in flexible parts — stresses. The joints can be given properties such as spring constants, damping, friction, and limits on movement. The gears pairs can also be given properties such as gear mesh stiffness, mesh damping, backlash, transmission error, and friction.
When the Multibody Dynamics 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 (only applicable for flexible domains). Then, from the
Physics toolbar, add features that implement other multibody dynamics properties. You can also right-click
Multibody Dynamics to select physics features from the context menu.
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
mbd.
From the Structural transient behavior list, select
Include inertial terms or
Quasi-static. Use
Quasi-static to treat the mechanical behavior as quasi-static (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.
Enter the coordinates for the Reference point for moment computation xref (variable
refpnt). The resulting moments (applied or as reactions) are then computed relative to this reference point. During the results and analysis stage, the coordinates can be changed in the
Parameters section in the result nodes.
By using the Physics Node Generation button (
), you can automatically create
Rigid Domain and
Gear nodes.
The Physics Node Generation has the following options:
The Create Rigid Domains option automatically creates
Rigid Domain nodes for geometrically disconnected objects. The
Rigid Domain Selections input controls for which geometric objects to create a rigid domain. It is by default set to
From physics interface, which creates rigid domain for all objects in the selection of the physics interface. All available domain selections are also listed in the
Rigid Domain Selections drop-down menu, which makes it possible to create rigid domains for a subset of objects. Select the
Include mass and moment of inertia node check box to add a
Mass and Moment of Inertia subnode to each auto-generated
Rigid Domain node. The
Density of the created rigid domains is automatically set to zero when this check box is selected.
The Create Gears option automatically creates a
Gear node for every gear geometry part present in the geometry sequence. This option requires that the gear geometry is built from the COMSOL
Part Library. The domain selection for the auto-generated
Gear node is same as the gear geometry part. Relevant definitions and parameters are retrieved from the gear geometry part, and entered in the settings for the corresponding gear node. The
Rigid Domain Selections and
Include mass and moment of inertia does not affect the
Create Gears option.
You can select a Body defining reference frame from the list of all rigid domains, attachments, and gears used in the model. It is used to visualize the relative motion of objects with respect to the selected body.
The variables u_ref,
v_ref, and
w_ref, being the three components of the displacement field with respect to the reference body, are available in the result menus. Also, the variables
mbd.disp_ref and
mbd.vel_ref contain the total displacement and velocity with respect to the reference body, respectively.
The last two rows of the table contain a summary. In the Total row, the number of DOFs, prescribed conditions, and constraints are summed.
The Net row contains the net number of degrees of freedom of the model — that is, the difference between all degrees of freedom and the constraints and prescribed motions. A negative net number of degrees of freedom indicates that the mechanism is overconstrained and is not shown. In that case, the net number of constraints is instead displayed.
The physics interface uses the global spatial components of the Displacement field u as dependent variables in the flexible domains. The default names for the components are (
u,
v,
w) in 3D. In 2D the component names are (
u,
v), and in 2D axial symmetry they are (
u,
w). You can, however, not use the “missing” component name in the 2D cases as a parameter or variable name because it is still used internally.
In the Multibody Dynamics 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, however, give significant reductions of the model size for a given mesh containing hexahedral, prism, or quadrilateral elements.
The discretization order applies to the flexible bodies. The default is to use Linear shape functions for the
Displacement field. If you want to compute stresses with good accuracy, increase the shape function order to
Quadratic serendipity or
Quadratic.