Cell Periodicity
Use the Cell Periodicity node to model a periodic unit cell or RVE (representative volume element) of a larger repetitive structure. Periodic or homogeneous boundary conditions are then used on the cell outer boundaries.
If more than one Cell Periodicity node is used, they must have disjoint selections. It is thus possible to model more than one cell in the same study step.
From the Cell Periodicity node, it is possible to compute the elasticity matrix, the coefficients of thermal expansion or hygroscopic swelling for the equivalent homogenized material.
The Cell Periodicity node is only available with some COMSOL products (see https://www.comsol.com/products/specifications/). Cell Periodicity is available for 3D and 2D.
Cell Properties
Select the Boundary conditions on the unit cell — Free expansion, Average strain, Average stress, or Mixed.
For all values of Boundary conditions, select the Type Periodic or Homogeneous.
When Periodic is selected, the average stresses or strains are periodic over the cell boundaries.
The Homogeneous option can be used to prescribe boundaries when the cell geometry does not represent a periodic microstructure.
Depending on the choice of the Boundary conditions on the cell, different options are available for the Calculate average properties list.
When Free expansion is selected from the Boundary conditions list, the cell is allowed to expand freely. This option is useful to determine the thermal expansion or hygroscopic swelling coefficients in heterogeneous media.
Select Calculate average propertiesNone, Coefficient of thermal expansion, or Coefficient of hygroscopic swelling.
In 2D, use either a Plane stress or a Generalized plane strain approximation to calculate averaged properties with the Free expansion boundary condition.
When Average strain is selected from the Boundary conditions list, it possible to derive homogenized elastic properties of media, such as a perforated plates, porous media or composites structures. In this case, six load cases need to be examined in 3D.
Enter values or expressions for the components of the Average strain tensor εavg. In a geometrically nonlinear analysis, the average strains are interpreted as Green–Lagrange strains.
Select Calculate average propertiesNone, Elasticity matrix, Standard (XX, YY, ZZ, XY, YZ, XZ), or Elasticity matrix, Voigt (XX, YY, ZZ, YZ, XZ, XY).
When Average stress is selected from the Boundary conditions list, it is possible to derive the homogenized compliance matrix. Enter values or expressions for the components of the Average stress tensor σavg. In a geometrically nonlinear analysis, the stresses are interpreted as second Piola–Kirchhoff stresses.
Select Calculate average propertiesNone, Compliance matrix, Standard (XX, YY, ZZ, XY, YZ, XZ), or Compliance matrix, Voigt (XX, YY, ZZ, YZ, XZ, XY).
In 2D, use a Generalized plane strain approximation to calculate averaged properties with the Average strain or Average stress boundary conditions.
When Mixed is selected from the Boundary conditions list, enter either the Average strain or Average stress for each tensor component.
When using a Plane stress or a Plane strain 2D approximation with the Mixed boundary condition, you only need to enter three components (XX, YY, XY) of the average strain or average stress tensors. For the Generalized plane strain approximation, you also need to enter the ZZ component.
For all values of Boundary conditions, select the Cell volumeAutomatic, Void volume fraction, or User defined.
When Automatic is selected, the cell volume is computed from the domain or domains selected in the Domain Selection section.
The Void volume fraction option can be used to scale the cell volume when there are voids that are not selected as domains. Enter the void volume fraction, f. This is the fraction of the cell that is occupied by voids. The volume of the cell is computed as
The User defined option allows to enter the cell volume V explicitly.
The variable <item>.vol contains the volume of the RVE used when computing average strains and stresses.
Load Group, Material, and Study Generation
You can create the elasticity matrix for a homogenized material, and make it accessible as a material to be used in other components. In order to do this, a number of fundamental load cases must be analyzed. The necessary definitions and steps required for such an analysis can be automated as described below.
On the Cell Properties section toolbar, there is an icon Automated Model Setup (). It has a list with three entries:
Create Load Groups and Study . This option can be selected when the Boundary conditions are Average strain or Average stress. When you select it, the following changes will be made to the model:
-
A number of load groups will be created under Global Definitions. They are collected in a group named Load Groups for Cell Periodicity. The load groups correspond to unit loads along different axes.
-
The Average strain or Average stress tensor in the Cell Properties section will be populated using the load group variables.
-
A new study, named Cell Periodicity Study, will be created. This study contains a stationary study step. In the Study Extensions section of the new Cell Periodicity Study, one load case is added for each load group.
Create Material by Reference (). This option can be selected when the Boundary conditions are Average strain or Average stress, and Calculate average properties is not set to None. When you select it, a new material will be created under Global>Materials. It contains the elasticity matrix or the compliance matrix. The name of this material is Homogeneous Material.
Create Material by Value (). This option is similar to Create Material by Reference, with the difference that the material properties are numerical values, and not references to variables in the node. The benefits of this option is faster computation when the material is used by other physics interfaces. Also, if the Cell Periodicity node is deleted, the material created from it can still be used in sequential computations.
Advanced
In this section, you can specify if the special cell periodicity study should also contain a parametric sweep. Select Add parametric sweepYes or No. The section is only shown if Boundary conditions is set to Average strain or Average stress.
If Yes is selected, then specify the Sweep typeSpecified Combinations or All combinations. Then, add parameters, their ranges, and their units on the Parameters table. The interpretation of the settings are the same as for the Parametric Sweep study step
For this type of analysis to work correctly, it is important that you do not edit the generated nodes manually. By clicking the Create button again, you can reset all settings in the generated nodes to their default values.
Load Cases and Effective Properties of Periodic Structures in the Structural Mechanics Modeling chapter.
Cell Periodicity in the Structural Mechanics Theory chapter.
Materials and Solid Mechanics Material Properties in the COMSOL Multiphysics Reference Manual.
Micromechanical Model of a Fiber Composite: Application Library path Structural_Mechanics_Module/Material_Models/micromechanical_model_of_a_fiber_composite
Micromechanics and Macromechanics of a Composite Cylinder: Application Library path Composite_Materials_Module/Tutorials/composite_cylinder_micromechanics_and_stress_analysis
Advanced
To display this section, click the Show More Options button () and select Advanced Physics Options in the Show More Options dialog box. Furthermore, this section is shown only if elasticity or compliance matrices are requested through Calculate average properties.
In that case, you can modify the behavior of how the studies for computing the homogenized matrices are generated by selecting a Study tagAutomatic or Legacy. The Legacy option is mainly for maintaining compatibility with versions prior to 6.0. With the default Automatic option, it is, for example, possible to set up nested homogenizations. A unique tag for the study will be composed by combining the physics interface tag and the tag of the Cell Periodicity node, for example solid1cp1std.
Constraint Settings
To display this section, click the Show More Options button () and select Advanced Physics Options in the Show More Options dialog box.
In the COMSOL Multiphysics Reference Manual:
Location in User Interface
Context Menus
Ribbon
Physics tab with Solid Mechanics selected: