FEM EquationElasticity: Difference between revisions
(update to new menu location) |
mNo edit summary |
||
(12 intermediate revisions by 3 users not shown) | |||
Line 4: | Line 4: | ||
<!--T:4--> |
<!--T:4--> |
||
{{Docnav |
{{Docnav |
||
|[[FEM_SolverZ88|Solver Z88]] |
|||
|[[FEM_EquationDeformation|Deformation equation]] |
|[[FEM_EquationDeformation|Deformation equation]] |
||
|[[FEM_CompEmEquations|Electromagnetic equations]] |
|||
|[[FEM_Workbench|FEM]] |
|[[FEM_Workbench|FEM]] |
||
|IconL= |
|IconL=FEM_SolverZ88.svg |
||
|IconR= |
|IconR=FEM_EquationDeformation.svg |
||
|IconC=Workbench_FEM.svg |
|IconC=Workbench_FEM.svg |
||
}} |
}} |
||
Line 24: | Line 24: | ||
<!--T:2--> |
<!--T:2--> |
||
This equation describes the mechanical properties of |
This equation describes the mechanical properties of solid bodies. |
||
<!--T:10--> |
<!--T:10--> |
||
Line 36: | Line 36: | ||
# Change the [[#Solver_Settings|equation's solver settings]] or the [[FEM_SolverElmer_SolverSettings|general solver settings]] if necessary. |
# Change the [[#Solver_Settings|equation's solver settings]] or the [[FEM_SolverElmer_SolverSettings|general solver settings]] if necessary. |
||
<!--T:34--> |
|||
'''Note''': For analyses of nonlinear |
'''Note''': For analyses of nonlinear elasticity you must use the [[Image:FEM_EquationDeformation.svg|32px]] [[FEM_EquationDeformation|Deformation equation]] ({{Version|0.21}}). The Elasticity equation is only for linear elasticity. |
||
<!--T:25--> |
|||
'''Note''': If you use more than one CPU core for the solver ({{Version|0.21}}), you cannot use the default solver settings. However, using just one CPU and the default solver settings is in many cases faster than using several CPUs because the elasticity solver is only fast when {{PropertyData|Linear Solver Type}} is set to ''Direct'' (the default, described [[FEM_SolverElmer_SolverSettings#Linear_System|here]]). For multi-CPU solving one can only use the {{PropertyData|Linear Direct Method}} ''MUMPS''. However, MUMPS is not freely available as a direct download. |
|||
==Solver Settings== <!--T:13--> |
==Solver Settings== <!--T:13--> |
||
Line 47: | Line 51: | ||
* {{PropertyData|Calculate Pangle}}: If the principal angles should be calculated. |
* {{PropertyData|Calculate Pangle}}: If the principal angles should be calculated. |
||
* {{PropertyData|Calculate Principal}}: If all stresses should be calculated. |
* {{PropertyData|Calculate Principal}}: If all stresses should be calculated. |
||
* {{PropertyData|Calculate Strains}}: If strains will be calculated. This will also calculate the stresses, even if {{PropertyData|Calculate Principal}} or {{PropertyData|Calculate Stresses}} is |
* {{PropertyData|Calculate Strains}}: If strains will be calculated. This will also calculate the stresses, even if {{PropertyData|Calculate Principal}} or {{PropertyData|Calculate Stresses}} is {{FALSE}}. |
||
* {{PropertyData|Calculate Stresses}}: If stresses should be calculated. Compared to {{PropertyData|Calculate Principal}} the Tresca yield criterion and the principal stress will not be calculated. |
* {{PropertyData|Calculate Stresses}}: If stresses should be calculated. Compared to {{PropertyData|Calculate Principal}} the Tresca yield criterion and the principal stress will not be calculated. |
||
* {{PropertyData|Constant Bulk System}}: See the Elmer manual for more info. |
* {{PropertyData|Constant Bulk System}}: See the Elmer manual for more info. |
||
* {{PropertyData|Displace Mesh}}: If mesh can be deformed. This is by default |
* {{PropertyData|Displace Mesh}}: If mesh can be deformed. This is by default {{TRUE}} and must be set to {{FALSE}} for eigenfrequency analyses. |
||
* {{PropertyData|Fix Displacement}}: If displacements or forces are set. thereby {{PropertyData|Model Lumping}} is automatically used. |
* {{PropertyData|Fix Displacement}}: If displacements or forces are set. thereby {{PropertyData|Model Lumping}} is automatically used. |
||
* {{PropertyData|Geometric Stiffness}}: Considers the geometric stiffness of the body. |
* {{PropertyData|Geometric Stiffness}}: Considers the geometric stiffness of the body. |
||
Line 57: | Line 61: | ||
* {{PropertyData|Model Lumping}}: Uses [https://en.wikipedia.org/wiki/Lumped-element_model model lumping]. |
* {{PropertyData|Model Lumping}}: Uses [https://en.wikipedia.org/wiki/Lumped-element_model model lumping]. |
||
* {{PropertyData|Model Lumping Filename}}: File to save the results from the model lumping. |
* {{PropertyData|Model Lumping Filename}}: File to save the results from the model lumping. |
||
* {{PropertyData|Stability Analysis}}: If |
* {{PropertyData|Stability Analysis}}: If {{TRUE}} {{PropertyData|Eigen Analysis}} becomes a stability analysis (buckling analysis). Otherwise a modal analysis is performed. |
||
* {{PropertyData|Update Transient System}}: See the Elmer manual for more info. |
* {{PropertyData|Update Transient System}}: See the Elmer manual for more info. |
||
* {{PropertyData|Variable}}: The variable for the elasticity equation. Only change this if {{PropertyData|Incompressible}} is set to |
* {{PropertyData|Variable}}: The variable for the elasticity equation. Only change this if {{PropertyData|Incompressible}} is set to {{TRUE}} in accordance to the Elmer manual. |
||
<!--T:16--> |
<!--T:16--> |
||
Eigenvalues: |
Eigenvalues: |
||
* {{PropertyData|Eigen Analysis}}: If an eigen analysis should be performed (calculation of eigenmodes and eigenfrequencies). |
* {{PropertyData|Eigen Analysis}}: If an eigen analysis should be performed (calculation of eigenmodes and eigenfrequencies). |
||
* {{PropertyData|Eigen System Complex}}: Should be |
* {{PropertyData|Eigen System Complex}}: Should be {{TRUE}} if the eigen system is complex. it must be {{FALSE}} for a damped eigen value analyses. |
||
* {{PropertyData|Eigen System Compute Residuals}}: Computes residuals of the eigen value system. |
* {{PropertyData|Eigen System Compute Residuals}}: Computes residuals of the eigen value system. |
||
* {{PropertyData|Eigen System Damped}}: Set a damped eigen analysis. Can only be used if {{PropertyData|[[FEM_SolverElmer_SolverSettings#Linear_System|Linear Solver Type]]}} is ''Iterative''. |
* {{PropertyData|Eigen System Damped}}: Set a damped eigen analysis. Can only be used if {{PropertyData|[[FEM_SolverElmer_SolverSettings#Linear_System|Linear Solver Type]]}} is ''Iterative''. |
||
Line 75: | Line 79: | ||
* {{PropertyData|Plane Stress}}: Computes solution according to the plane stress situation. Applies only for 2D geometry. |
* {{PropertyData|Plane Stress}}: Computes solution according to the plane stress situation. Applies only for 2D geometry. |
||
== |
==Analysis Feature Information== <!--T:18--> |
||
<!--T:19--> |
<!--T:19--> |
||
The elasticity equation takes the following |
The elasticity equation takes the following analysis features into account if they are set: |
||
<!--T:20--> |
<!--T:20--> |
||
* [[Image:FEM_ConstraintFixed.svg|32px]] [[FEM_ConstraintFixed| |
* [[Image:FEM_ConstraintFixed.svg|32px]] [[FEM_ConstraintFixed|Fixed boundary condition]] |
||
* [[Image:FEM_ConstraintDisplacement.svg|32px]] [[FEM_ConstraintDisplacement| |
* [[Image:FEM_ConstraintDisplacement.svg|32px]] [[FEM_ConstraintDisplacement|Displacement boundary condition]] |
||
* [[Image:FEM_ConstraintForce.svg|32px]] [[FEM_ConstraintForce| |
* [[Image:FEM_ConstraintForce.svg|32px]] [[FEM_ConstraintForce|Force load]] |
||
* [[Image:FEM_ConstraintInitialTemperature.svg|32px]] [[FEM_ConstraintInitialTemperature| |
* [[Image:FEM_ConstraintInitialTemperature.svg|32px]] [[FEM_ConstraintInitialTemperature|Initial temperature condition]] |
||
* [[Image:FEM_ConstraintPressure.svg|32px]] [[FEM_ConstraintPressure| |
* [[Image:FEM_ConstraintPressure.svg|32px]] [[FEM_ConstraintPressure|Pressure load]] |
||
* [[Image:FEM_ConstraintSelfWeight.svg|32px]] [[FEM_ConstraintSelfWeight| |
* [[Image:FEM_ConstraintSelfWeight.svg|32px]] [[FEM_ConstraintSelfWeight|Gravity load]] |
||
* [[Image:FEM_ConstraintSpring.svg|32px]] [[FEM_ConstraintSpring| |
* [[Image:FEM_ConstraintSpring.svg|32px]] [[FEM_ConstraintSpring|Spring]] |
||
===Note=== <!--T:21--> |
===Note=== <!--T:21--> |
||
<!--T:22--> |
<!--T:22--> |
||
* Except for calculations in 2D, for all above |
* Except for calculations in 2D, for all the above analysis features it is important that they act on a face. Features in 3D set to lines or vertices are not recognized by the Elmer solver. |
||
==Eigenmode Analysis== <!--T:23--> |
==Eigenmode Analysis== <!--T:23--> |
||
Line 98: | Line 102: | ||
<!--T:24--> |
<!--T:24--> |
||
To perform an eigenmode analysis (calculation if the eigenmodes and eigenfrequencies), you need to |
To perform an eigenmode analysis (calculation if the eigenmodes and eigenfrequencies), you need to |
||
# Set {{PropertyData|Eigen Analysis}} |
# Set {{PropertyData|Eigen Analysis}} to {{TRUE}} |
||
# Set {{PropertyData|Displace Mesh}} |
# Set {{PropertyData|Displace Mesh}} to {{FALSE}} |
||
# Set {{PropertyData|Eigen System Values}} |
# Set {{PropertyData|Eigen System Values}} to the highest number of eigenmodes you are interested in. The smaller this number the shorter the solver runtime since higher modes can be omitted from computation. |
||
# Add a [[FEM_ConstraintFixed| |
# Add a [[FEM_ConstraintFixed|fixed boundary condition]] and set at least one face of the body as fixed. |
||
# Run the solver. |
# Run the solver. |
||
<!--T: |
<!--T:35--> |
||
It is highly recommended to use {{PropertyData|Linear Solver Type}} set to ''Direct'' (the default) because this is much faster and the results are more accurate. |
|||
'''Note''': If you use more than one CPU core for the solver ({{Version|1.0}}), you cannot use ''Umfpack'', the only direct method for parallel solving is ''MUMPS''.</br>Also note that iterative solving is not recommended for eigenmode analysis. Therefore either only use one CPU core or install the ''MUMPS'' module to Elmer. |
|||
==Buckling Analysis== <!--T:26--> |
==Buckling Analysis== <!--T:26--> |
||
Line 111: | Line 115: | ||
<!--T:27--> |
<!--T:27--> |
||
To perform a buckling analysis, you need to do the same as for an [[#Eigenmode_Analysis|Eigenmode Analysis]], and additionally: |
To perform a buckling analysis, you need to do the same as for an [[#Eigenmode_Analysis|Eigenmode Analysis]], and additionally: |
||
* Set {{PropertyData|Stability Analysis}} to |
* Set {{PropertyData|Stability Analysis}} to {{TRUE}} |
||
==Results== <!--T:28--> |
==Results== <!--T:28--> |
||
<!--T:29--> |
<!--T:29--> |
||
The available results depend on the [[#Solver_Settings|solver settings]]. If none of the {{PropertyData|Calculate *}} settings was set to |
The available results depend on the [[#Solver_Settings|solver settings]]. If none of the {{PropertyData|Calculate *}} settings was set to {{TRUE}}, only the displacement is calculated. Otherwise also the corresponding results will be available. If {{PropertyData|Eigen Analysis}} was set to {{TRUE}} all results will be available for every calculated eigenmode. |
||
<!--T:30--> |
<!--T:30--> |
||
If {{PropertyData|Eigen Analysis}} was set to |
If {{PropertyData|Eigen Analysis}} was set to {{TRUE}}, the eigenfrequencies are output at the end of the solver log in the solver dialog and also in the document '''SolverElmerOutput''' that will be created in the tree view after the solver has finished. |
||
<!--T:31--> |
<!--T:31--> |
||
Line 131: | Line 135: | ||
<!--T:33--> |
<!--T:33--> |
||
whereas <math>\vec{u}</math> is the eigenvector and <math>c</math> is a complex number. |
whereas <math>\vec{u}</math> is the eigenvector and <math>c</math> is a complex number. |
||
<!--T:5--> |
<!--T:5--> |
||
{{Docnav |
{{Docnav |
||
|[[FEM_SolverZ88|Solver Z88]] |
|||
|[[FEM_EquationDeformation|Deformation equation]] |
|[[FEM_EquationDeformation|Deformation equation]] |
||
|[[FEM_CompEmEquations|Electromagnetic equations]] |
|||
|[[FEM_Workbench|FEM]] |
|[[FEM_Workbench|FEM]] |
||
|IconL= |
|IconL=FEM_SolverZ88.svg |
||
|IconR= |
|IconR=FEM_EquationDeformation.svg |
||
|IconC=Workbench_FEM.svg |
|IconC=Workbench_FEM.svg |
||
}} |
}} |
||
Line 145: | Line 150: | ||
{{FEM Tools navi{{#translation:}}}} |
{{FEM Tools navi{{#translation:}}}} |
||
{{Userdocnavi{{#translation:}}}} |
{{Userdocnavi{{#translation:}}}} |
||
{{clear}} |
Latest revision as of 13:25, 24 November 2023
FEM EquationElasticity |
Menu location |
---|
Solve → Mechanical equations → Elasticity equation |
Workbenches |
FEM |
Default shortcut |
None |
Introduced in version |
0.17 |
See also |
Deformation equation, FEM tutorial |
Description
This equation describes the mechanical properties of solid bodies.
For info about the math of the equation, see the Elmer models manual, section Linear Elasticity.
Usage
- After adding an Elmer solver as described here, select it in the tree view.
- Now either use the toolbar button or the menu Solve → Mechanical equations → Elasticity equation.
- Change the equation's solver settings or the general solver settings if necessary.
Note: For analyses of nonlinear elasticity you must use the Deformation equation (introduced in version 0.21). The Elasticity equation is only for linear elasticity.
Note: If you use more than one CPU core for the solver (introduced in version 0.21), you cannot use the default solver settings. However, using just one CPU and the default solver settings is in many cases faster than using several CPUs because the elasticity solver is only fast when DataLinear Solver Type is set to Direct (the default, described here). For multi-CPU solving one can only use the DataLinear Direct Method MUMPS. However, MUMPS is not freely available as a direct download.
Solver Settings
For the general solver settings, see the Elmer solver settings.
The elasticity equation provides these special settings:
- DataCalculate Pangle: If the principal angles should be calculated.
- DataCalculate Principal: If all stresses should be calculated.
- DataCalculate Strains: If strains will be calculated. This will also calculate the stresses, even if DataCalculate Principal or DataCalculate Stresses is
false
. - DataCalculate Stresses: If stresses should be calculated. Compared to DataCalculate Principal the Tresca yield criterion and the principal stress will not be calculated.
- DataConstant Bulk System: See the Elmer manual for more info.
- DataDisplace Mesh: If mesh can be deformed. This is by default
true
and must be set tofalse
for eigenfrequency analyses. - DataFix Displacement: If displacements or forces are set. thereby DataModel Lumping is automatically used.
- DataGeometric Stiffness: Considers the geometric stiffness of the body.
- DataIncompressible: Computation of incompressible material in connection with viscoelastic Maxwell material and a custom DataVariable.
- DataMaxwell Material: Compute the viscoelastic material model.
- DataModel Lumping: Uses model lumping.
- DataModel Lumping Filename: File to save the results from the model lumping.
- DataStability Analysis: If
true
DataEigen Analysis becomes a stability analysis (buckling analysis). Otherwise a modal analysis is performed. - DataUpdate Transient System: See the Elmer manual for more info.
- DataVariable: The variable for the elasticity equation. Only change this if DataIncompressible is set to
true
in accordance to the Elmer manual.
Eigenvalues:
- DataEigen Analysis: If an eigen analysis should be performed (calculation of eigenmodes and eigenfrequencies).
- DataEigen System Complex: Should be
true
if the eigen system is complex. it must befalse
for a damped eigen value analyses. - DataEigen System Compute Residuals: Computes residuals of the eigen value system.
- DataEigen System Damped: Set a damped eigen analysis. Can only be used if DataLinear Solver Type is Iterative.
- DataEigen System Select: Selection of which eigenvalues are computed. Note that the selection of Largest* cause an infinite run for recent Elmer solver (as of August 2022).
- DataEigen System Tolerance: Convergence tolerance for iterative eigensystem solve. The default is 100 times the DataLinear Tolerance.
- DataEigen System Values: The number of the highest eigenmode that should be calculated.
Equation:
- DataPlane Stress: Computes solution according to the plane stress situation. Applies only for 2D geometry.
Analysis Feature Information
The elasticity equation takes the following analysis features into account if they are set:
- Fixed boundary condition
- Displacement boundary condition
- Force load
- Initial temperature condition
- Pressure load
- Gravity load
- Spring
Note
- Except for calculations in 2D, for all the above analysis features it is important that they act on a face. Features in 3D set to lines or vertices are not recognized by the Elmer solver.
Eigenmode Analysis
To perform an eigenmode analysis (calculation if the eigenmodes and eigenfrequencies), you need to
- Set DataEigen Analysis to
true
- Set DataDisplace Mesh to
false
- Set DataEigen System Values to the highest number of eigenmodes you are interested in. The smaller this number the shorter the solver runtime since higher modes can be omitted from computation.
- Add a fixed boundary condition and set at least one face of the body as fixed.
- Run the solver.
It is highly recommended to use DataLinear Solver Type set to Direct (the default) because this is much faster and the results are more accurate.
Buckling Analysis
To perform a buckling analysis, you need to do the same as for an Eigenmode Analysis, and additionally:
- Set DataStability Analysis to
true
Results
The available results depend on the solver settings. If none of the DataCalculate * settings was set to true
, only the displacement is calculated. Otherwise also the corresponding results will be available. If DataEigen Analysis was set to true
all results will be available for every calculated eigenmode.
If DataEigen Analysis was set to true
, the eigenfrequencies are output at the end of the solver log in the solver dialog and also in the document SolverElmerOutput that will be created in the tree view after the solver has finished.
Note: The eigenmode displacement vector has an arbitrary value since the result is
whereas is the eigenvector and is a complex number.
- Materials: Solid, Fluid, Nonlinear mechanical, Reinforced (concrete); Material editor
- Element geometry: Beam (1D), Beam rotation (1D), Shell (2D), Fluid flow (1D)
Constraints
- Electromagnetic: Electrostatic potential, Current density, Magnetization
- Geometrical: Plane rotation, Section print, Transform
- Mechanical: Fixed, Displacement, Contact, Tie, Spring, Force, Pressure, Centrif, Self weight
- Thermal: Initial temperature, Heat flux, Temperature, Body heat source
- Overwrite Constants: Constant vacuum permittivity
- Solve: CalculiX Standard, Elmer, Mystran, Z88; Equations: Deformation, Elasticity, Electrostatic, Electricforce, Magnetodynamic, Magnetodynamic 2D, Flow, Flux, Heat; Solver: Solver control, Solver run
- Results: Purge, Show; Postprocessing: Apply changes, Pipeline from result, Warp filter, Scalar clip filter, Function cut filter, Region clip filter, Contours filter, Line clip filter, Stress linearization plot, Data at point clip filter, Filter function plane, Filter function sphere, Filter function cylinder, Filter function box
- Additional: Preferences; FEM Install, FEM Mesh, FEM Solver, FEM CalculiX, FEM Concrete; FEM Element Types
- Getting started
- Installation: Download, Windows, Linux, Mac, Additional components, Docker, AppImage, Ubuntu Snap
- Basics: About FreeCAD, Interface, Mouse navigation, Selection methods, Object name, Preferences, Workbenches, Document structure, Properties, Help FreeCAD, Donate
- Help: Tutorials, Video tutorials
- Workbenches: Std Base, Assembly, BIM, CAM, Draft, FEM, Inspection, Mesh, OpenSCAD, Part, PartDesign, Points, Reverse Engineering, Robot, Sketcher, Spreadsheet, Surface, TechDraw, Test Framework
- Hubs: User hub, Power users hub, Developer hub