kindred/create
1
1
Fork 0
Code Issues 30 Pull Requests Actions Packages Projects 9 Releases 2 Wiki Activity

[FEM] Elmer: next step to fix eigenfrequency analysis

Browse Source 
- add missing parameters to perform modal analyses
- set mandatory complex statement
This commit is contained in:
Uwe
2022-08-06 18:27:51 +02:00
parent d88729d2c4
commit 29197f503d
2 changed files with 140 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

76
src/Mod/Fem/femsolver/elmer/equations/elasticity.py
View File

@@ -32,6 +32,8 @@ from femtools import femutils
from ... import equationbase
from . import linear
EIGEN_SYSTEM_SELECT = ["Smallest Magnitude", "Largest Magnitude", "Smallest Real Part",\
"Largest Real Part", "Smallest Imag Part", "Largest Imag Part"]
def create(doc, name="Elasticity"):
return femutils.createObject(
@@ -87,13 +89,58 @@ class Proxy(linear.Proxy, equationbase.ElasticityProxy):
obj.addProperty(
"App::PropertyBool",
"EigenAnalysis",
"Elasticity",
"Eigen Values",
"If true, modal analysis"
)
obj.addProperty(
"App::PropertyBool",
"EigenSystemComplex",
"Eigen Values",
(
"Should be true if eigen system is complex\n"
"Must be false for a damped eigen value analysis."
)
)
obj.addProperty(
"App::PropertyBool",
"EigenSystemComputeResiduals",
"Eigen Values",
"Computes residuals of eigen value system"
)
obj.addProperty(
"App::PropertyBool",
"EigenSystemDamped",
"Eigen Values",
(
"Set a damped eigen analysis. Can only be\n"
"used if 'Linear Solver Type' is 'Iterative'."
)
)
obj.addProperty(
"App::PropertyIntegerConstraint",
"EigenSystemMaxIterations",
"Eigen Values",
"Max iterations for iterative eigensystem solver"
)
obj.addProperty(
"App::PropertyEnumeration",
"EigenSystemSelect",
"Eigen Values",
"Which eigenvalues are computed"
)
obj.addProperty(
"App::PropertyFloat",
"EigenSystemTolerance",
"Eigen Values",
(
"Convergence tolerance for iterative eigensystem solve\n"
"Default is 100 times the 'Linear Tolerance'"
)
)
obj.addProperty(
"App::PropertyInteger",
"EigenSystemValues",
"Elasticity",
"Eigen Values",
"Number of lowest eigen modes"
)
obj.addProperty(
@@ -160,14 +207,35 @@ class Proxy(linear.Proxy, equationbase.ElasticityProxy):
)
)
obj.addProperty(
"App::PropertyBool",
"PlaneStress",
"Equation",
(
"Computes solution according to plane\nstress situation.\n"
"Applies only for 2D geometry."
)
)
obj.EigenSystemValues = 5
obj.Priority = 10
obj.CalculatePrincipal = True
obj.DisplaceMesh = True
# according to the Elmer tutorial and forum, for stresses direct solving
# is recommended -> test showed 10 times faster and even more accurate
# according to Elmer tutorial and forum, for stresses direct solving
# is recommended -> tests showed 10 times faster and even more accurate
obj.LinearSolverType = "Direct"
obj.LinearDirectMethod = "Umfpack"
# according to Elmer manual should be true except there is damping
# we set damping by default to false
obj.EigenSystemComplex = True
# according to Elmer manual the default is 300
obj.EigenSystemMaxIterations = (300, 1, int(1e8), 1)
# according to Elmer manual the default is Smallest magnitude
obj.EigenSystemSelect = EIGEN_SYSTEM_SELECT
obj.EigenSystemSelect = "Smallest Magnitude"
# according to Elmer manual default is 100 times the Linear Tolerance
# since this is a small value we must set an expression, see linear.py for background
obj.setExpression("EigenSystemTolerance", str(100 * obj.LinearTolerance))
obj.Variable = "Displacement"

70
src/Mod/Fem/femsolver/elmer/writer.py
View File

@@ -745,6 +745,14 @@ class Writer(object):
s["Displace mesh"] = equation.DisplaceMesh
s["Eigen Analysis"] = equation.EigenAnalysis
if equation.EigenAnalysis is True:
s["Eigen System Convergence Tolerance"] = \
equation.EigenSystemTolerance
s["Eigen System Complex"] = equation.EigenSystemComplex
if equation.EigenSystemComputeResiduals is True:
s["Eigen System Compute Residuals"] = equation.EigenSystemComputeResiduals
s["Eigen System Damped"] = equation.EigenSystemDamped
s["Eigen System Max Iterations"] = equation.EigenSystemMaxIterations
s["Eigen System Select"] = equation.EigenSystemSelect
s["Eigen System Values"] = equation.EigenSystemValues
if equation.FixDisplacement is True:
s["Fix Displacement"] = equation.FixDisplacement
@@ -811,19 +819,77 @@ class Writer(object):
equation.addProperty(
"App::PropertyBool",
"EigenAnalysis",
"Elasticity",
"Eigen Values",
"If true, modal analysis"
)
if hasattr(equation, "DoFrequencyAnalysis"):
equation.EigenAnalysis = equation.DoFrequencyAnalysis
equation.removeProperty("DoFrequencyAnalysis")
if not hasattr(equation, "EigenSystemComplex"):
equation.addProperty(
"App::PropertyBool",
"EigenSystemComplex",
"Eigen Values",
(
"Should be true if eigen system is complex\n"
"Must be false for a damped eigen value analysis."
)
)
equation.EigenSystemComplex = True
if not hasattr(equation, "EigenSystemComputeResiduals"):
equation.addProperty(
"App::PropertyBool",
"EigenSystemComputeResiduals",
"Eigen Values",
"Computes residuals of eigen value system"
)
if not hasattr(equation, "EigenSystemDamped"):
equation.addProperty(
"App::PropertyBool",
"EigenSystemDamped",
"Eigen Values",
(
"Set a damped eigen analysis. Can only be\n"
"used if 'Linear Solver Type' is 'Iterative'."
)
)
if not hasattr(equation, "EigenSystemMaxIterations"):
equation.addProperty(
"App::PropertyIntegerConstraint",
"EigenSystemMaxIterations",
"Eigen Values",
"Max iterations for iterative eigensystem solver"
)
equation.EigenSystemMaxIterations = (300, 1, int(1e8), 1)
EIGEN_SYSTEM_SELECT = ["Smallest Magnitude", "Largest Magnitude", "Smallest Real Part",\
"Largest Real Part", "Smallest Imag Part", "Largest Imag Part"]
if not hasattr(equation, "EigenSystemSelect"):
equation.addProperty(
"App::PropertyEnumeration",
"EigenSystemSelect",
"Eigen Values",
"Which eigenvalues are computed"
)
equation.EigenSystemSelect = EIGEN_SYSTEM_SELECT
equation.EigenSystemSelect = "Smallest Magnitude"
if not hasattr(equation, "EigenSystemTolerance"):
equation.addProperty(
"App::PropertyFloat",
"EigenSystemTolerance",
"Eigen Values",
(
"Convergence tolerance for iterative eigensystem solve\n"
"Default is 100 times the 'Linear Tolerance'"
)
)
equation.setExpression("EigenSystemTolerance", str(100 * equation.LinearTolerance))
if not hasattr(equation, "EigenSystemValues"):
# EigenmodesCount was renamed to EigenSystemValues
# to follow the Elmer manual
equation.addProperty(
"App::PropertyInteger",
"EigenSystemValues",
"Elasticity",
"Eigen Values",
"Number of lowest eigen modes"
)
if hasattr(equation, "EigenmodesCount"):