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

Merge branch 'FreeCAD:master' into master

Browse Source
This commit is contained in:
Zolko-123
2021-08-02 22:41:11 +02:00
committed by GitHub
parent da0c6cc545 6dbd8fb54d
commit 0085f276b2
138 changed files with 9441 additions and 485 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
src/Gui/Macro.h
View File

@@ -77,7 +77,7 @@ public:
void commit(void);
/// cancels the recording session
void cancel(void);
/// indicates if a macro recording in in progress
/// indicates if a macro recording is in progress
bool isOpen(void) const {return openMacro;}
/// insert a new line in the macro
void addLine(LineType Type,const char* sLine,bool pending=false);

2
src/Mod/Arch/ArchComponent.py
View File

@@ -1537,7 +1537,7 @@ class ViewProviderComponent:
return False
def colorize(self,obj,force=False):
"""If an object is a clone, set it it to copy the color of its parent.
"""If an object is a clone, set it to copy the color of its parent.
Only change the color of the clone if the clone and its parent have
colors that are distinguishably different from each other.

2
src/Mod/Draft/DraftTools.py
View File

@@ -105,7 +105,7 @@ FreeCADGui.updateLocale()
# sets the default working plane
plane = WorkingPlane.plane()
FreeCAD.DraftWorkingPlane = plane
defaultWP = Draft.getParam("defaultWP",1)
defaultWP = Draft.getParam("defaultWP",0)
if defaultWP == 1: plane.alignToPointAndAxis(Vector(0,0,0), Vector(0,0,1), 0)
elif defaultWP == 2: plane.alignToPointAndAxis(Vector(0,0,0), Vector(0,1,0), 0)
elif defaultWP == 3: plane.alignToPointAndAxis(Vector(0,0,0), Vector(1,0,0), 0)

6
src/Mod/Draft/draftfunctions/move.py
View File

@@ -83,9 +83,9 @@ def move(objectslist, vector, copy=False):
newobj = make_copy.make_copy(obj)
else:
newobj = obj
newobj.X = obj.X + real_vector.x
newobj.Y = obj.Y + real_vector.y
newobj.Z = obj.Z + real_vector.z
newobj.X = obj.X.Value + real_vector.x
newobj.Y = obj.Y.Value + real_vector.y
newobj.Z = obj.Z.Value + real_vector.z
elif obj.isDerivedFrom("App::DocumentObjectGroup"):
pass

10
src/Mod/Draft/draftgeoutils/geometry.py
View File

@@ -165,11 +165,17 @@ def findDistance(point, edge, strict=False):
return None
if strict and ve2:
# Note 1: DraftVecUtils.angle(App.Vector(1, 1, 0)) => -0.7854
# Note 2: Angles are in the +pi to -pi range.
ang1 = DraftVecUtils.angle(ve1.sub(center))
ang2 = DraftVecUtils.angle(ve2.sub(center))
angpt = DraftVecUtils.angle(newpoint.sub(center))
if ((angpt <= ang2 and angpt >= ang1)
or (angpt <= ang1 and angpt >= ang2)):
if ang1 >= ang2: # Arc does not cross the 9 o'clock point.
if ang1 >= angpt and angpt >= ang2:
return dist
else:
return None
elif ang1 >= angpt or angpt >= ang2:
return dist
else:
return None

6
src/Mod/Draft/draftguitools/gui_base_original.py
View File

@@ -180,6 +180,7 @@ class DraftTool:
self.ui.sourceCmd = None
if self.planetrack:
self.planetrack.finalize()
App.DraftWorkingPlane.restore()
if hasattr(Gui, "Snapper"):
Gui.Snapper.off()
if self.call:
@@ -302,4 +303,9 @@ class Modifier(DraftTool):
super(Modifier, self).__init__()
self.copymode = False
def Activated(self, name="None", noplanesetup=False, is_subtool=False):
super(Modifier, self).Activated(name, noplanesetup, is_subtool)
# call DraftWorkingPlane.save to sync with
# DraftWorkingPlane.restore called in finish method
App.DraftWorkingPlane.save()
## @}

2
src/Mod/Draft/draftguitools/gui_edit_draft_objects.py
View File

@@ -344,7 +344,7 @@ class DraftCircleGuiTools(GuiTools):
]
elif node_idx == 1: # user is over arc start point
return [
("mset first angle", lambda: self.handle_set_first_angle(edit_command, obj, node_idx)),
("set first angle", lambda: self.handle_set_first_angle(edit_command, obj, node_idx)),
]
elif node_idx == 2: # user is over arc end point
return [

46
src/Mod/Draft/draftguitools/gui_groups.py
View File

@@ -254,14 +254,13 @@ class SetAutoGroup(gui_base.GuiCommandSimplest):
# and globally initialized in the `Gui` namespace to run
# some actions.
# If there is only a group selected, it runs the `AutoGroup` method.
params = App.ParamGet("User parameter:BaseApp/Preferences/Mod/Draft")
self.ui = Gui.draftToolBar
s = Gui.Selection.getSelection()
if len(s) == 1:
if (utils.get_type(s[0]) == "Layer") or \
(App.ParamGet("User parameter:BaseApp/Preferences/Mod/Draft").GetBool("AutogroupAddGroups", False)
and (s[0].isDerivedFrom("App::DocumentObjectGroup")
or utils.get_type(s[0]) in ["Site", "Building",
"Floor", "BuildingPart"])):
if (utils.get_type(s[0]) == "Layer"
or (params.GetBool("AutogroupAddGroups", False)
and groups.is_group(s[0]))):
self.ui.setAutoGroup(s[0].Name)
return
@@ -269,27 +268,26 @@ class SetAutoGroup(gui_base.GuiCommandSimplest):
# including the options "None" and "Add new layer".
self.groups = ["None"]
gn = [o.Name for o in self.doc.Objects if utils.get_type(o) == "Layer"]
if App.ParamGet("User parameter:BaseApp/Preferences/Mod/Draft").GetBool("AutogroupAddGroups", False):
if params.GetBool("AutogroupAddGroups", False):
gn.extend(groups.get_group_names())
if gn:
self.groups.extend(gn)
self.labels = [translate("draft", "None")]
self.icons = [self.ui.getIcon(":/icons/button_invalid.svg")]
for g in gn:
o = self.doc.getObject(g)
if o:
self.labels.append(o.Label)
self.icons.append(o.ViewObject.Icon)
self.labels.append(translate("draft", "Add new Layer"))
self.icons.append(self.ui.getIcon(":/icons/document-new.svg"))
self.groups.extend(gn)
self.labels = [translate("draft", "None")]
self.icons = [self.ui.getIcon(":/icons/button_invalid.svg")]
for g in gn:
o = self.doc.getObject(g)
if o:
self.labels.append(o.Label)
self.icons.append(o.ViewObject.Icon)
self.labels.append(translate("draft", "Add new Layer"))
self.icons.append(self.ui.getIcon(":/icons/document-new.svg"))
# With the lists created is uses the interface
# to pop up a menu with layer options.
# Once the desired option is chosen
# it launches the `proceed` method.
self.ui.sourceCmd = self
pos = self.ui.autoGroupButton.mapToGlobal(QtCore.QPoint(0, self.ui.autoGroupButton.geometry().height()))
self.ui.popupMenu(self.labels, self.icons, pos)
# With the lists created is uses the interface
# to pop up a menu with layer options.
# Once the desired option is chosen
# it launches the `proceed` method.
self.ui.sourceCmd = self
pos = self.ui.autoGroupButton.mapToGlobal(QtCore.QPoint(0, self.ui.autoGroupButton.geometry().height()))
self.ui.popupMenu(self.labels, self.icons, pos)
def proceed(self, labelname):
"""Set the defined autogroup, or create a new layer.

5
src/Mod/Draft/draftguitools/gui_offset.py
View File

@@ -110,7 +110,10 @@ class Offset(gui_base_original.Modifier):
self.mode = "Circle"
self.center = self.shape.Edges[0].Curve.Center
self.ghost.setCenter(self.center)
self.ghost.setStartAngle(math.radians(self.sel.FirstAngle))
if self.sel.FirstAngle <= self.sel.LastAngle:
self.ghost.setStartAngle(math.radians(self.sel.FirstAngle))
else:
self.ghost.setStartAngle(math.radians(self.sel.FirstAngle) - 2 * math.pi)
self.ghost.setEndAngle(math.radians(self.sel.LastAngle))
elif utils.getType(self.sel) == "BSpline":
self.ghost = trackers.bsplineTracker(points=self.sel.Points)

4
src/Mod/Draft/draftguitools/gui_planeproxy.py
View File

@@ -47,8 +47,8 @@ class Draft_WorkingPlaneProxy:
"""Set icon, menu and tooltip."""
d = {'Pixmap': 'Draft_PlaneProxy',
'MenuText': QT_TRANSLATE_NOOP("Draft_SetWorkingPlaneProxy","Create working plane proxy"),
'ToolTip': QT_TRANSLATE_NOOP("Draft_SetWorkingPlaneProxy","Creates a proxy object from the current working plane.\nOnce the object is created double click it in the tree view to restore the camera position and objects' visibilities.\nThen you can use it to save a different camera position and objects' states any time you need.")}
'MenuText': QT_TRANSLATE_NOOP("Draft_WorkingPlaneProxy","Create working plane proxy"),
'ToolTip': QT_TRANSLATE_NOOP("Draft_WorkingPlaneProxy","Creates a proxy object from the current working plane.\nOnce the object is created double click it in the tree view to restore the camera position and objects' visibilities.\nThen you can use it to save a different camera position and objects' states any time you need.")}
return d
def IsActive(self):

4
src/Mod/Draft/draftguitools/gui_selectplane.py
View File

@@ -127,6 +127,9 @@ class Draft_SelectPlane:
self.taskd.form.fieldGridExtension.valueChanged.connect(self.onSetExtension)
self.taskd.form.fieldSnapRadius.valueChanged.connect(self.onSetSnapRadius)
# save previous WP to ensure back to the last used when restored
FreeCAD.DraftWorkingPlane.save()
# Try to find a WP from the current selection
if FreeCADGui.Selection.getSelectionEx(FreeCAD.ActiveDocument.Name):
if self.handle():
@@ -161,7 +164,6 @@ class Draft_SelectPlane:
# Reset everything else
FreeCADGui.Control.closeDialog()
FreeCAD.DraftWorkingPlane.restore()
FreeCADGui.ActiveDocument.resetEdit()
return True

10
src/Mod/Draft/draftguitools/gui_setstyle.py
View File

@@ -192,7 +192,7 @@ class Draft_SetStyle_TaskPanel:
vobj = obj.ViewObject
if vobj:
if "LineColor" in vobj.PropertiesList:
vobj.LineColor = self.form.LineColor.property("color").rgb()<<8
vobj.LineColor = self.form.LineColor.property("color").getRgbF()
if "LineWidth" in vobj.PropertiesList:
vobj.LineWidth = self.form.LineWidth.value()
if "DrawStyle" in vobj.PropertiesList:
@@ -208,7 +208,7 @@ class Draft_SetStyle_TaskPanel:
except Exception:
pass
if "ShapeColor" in vobj.PropertiesList:
vobj.ShapeColor = self.form.ShapeColor.property("color").rgb()<<8
vobj.ShapeColor = self.form.ShapeColor.property("color").getRgbF()
if "Transparency" in vobj.PropertiesList:
vobj.Transparency = self.form.Transparency.value()
if "FontName" in vobj.PropertiesList:
@@ -218,7 +218,7 @@ class Draft_SetStyle_TaskPanel:
if "FontSize" in vobj.PropertiesList:
vobj.FontSize = FreeCAD.Units.Quantity(self.form.TextSize.text()).Value
if "TextColor" in vobj.PropertiesList:
vobj.TextColor = self.form.TextColor.property("color").rgb()<<8
vobj.TextColor = self.form.TextColor.property("color").getRgbF()
if "ArrowType" in vobj.PropertiesList:
vobj.ArrowType = ["Dot", "Circle", "Arrow", "Tick", "Tick-2"][self.form.ArrowStyle.currentIndex()]
if "ArrowSize" in vobj.PropertiesList:
@@ -243,7 +243,7 @@ class Draft_SetStyle_TaskPanel:
vobj = obj.ViewObject
vobj.FontName = self.form.TextFont.currentFont().family()
vobj.FontSize = FreeCAD.Units.Quantity(self.form.TextSize.text()).Value
vobj.LineColor = self.form.TextColor.property("color").rgb()<<8
vobj.LineColor = self.form.TextColor.property("color").getRgbF()
vobj.ArrowType = ["Dot", "Circle", "Arrow", "Tick", "Tick-2"][self.form.ArrowStyle.currentIndex()]
vobj.ArrowSize = FreeCAD.Units.Quantity(self.form.ArrowSize.text()).Value
vobj.ShowUnit = self.form.ShowUnit.isChecked()
@@ -255,7 +255,7 @@ class Draft_SetStyle_TaskPanel:
vobj = obj.ViewObject
vobj.FontName = self.form.TextFont.currentFont().family()
vobj.FontSize = FreeCAD.Units.Quantity(self.form.TextSize.text()).Value
vobj.TextColor = self.form.TextColor.property("color").rgb()<<8
vobj.TextColor = self.form.TextColor.property("color").getRgbF()
vobj.LineSpacing = self.form.LineSpacing.value()
def onLoadStyle(self,index):

19
src/Mod/Draft/draftguitools/gui_tool_utils.py
View File

@@ -72,10 +72,9 @@ def select_object(arg):
If it is of type Keyboard and the `ESCAPE` key, it runs the `finish`
method of the active command.
If it is of type Mouse button and `BUTTON1` press,
it captures the position of the cursor (x, y)
and the object below that cursor to add it to the active selection;
then it runs the `proceed` method of the active command
If Ctrl key is pressed, multiple selection is enabled until the
button is released.
Then it runs the `proceed` method of the active command
to continue with the command's logic.
"""
if arg["Type"] == "SoKeyboardEvent":
@@ -83,16 +82,8 @@ def select_object(arg):
App.activeDraftCommand.finish()
# TODO: this part raises a coin3D warning about scene traversal.
# It needs to be fixed.
elif arg["Type"] == "SoMouseButtonEvent":
if arg["State"] == "DOWN" and arg["Button"] == "BUTTON1":
cursor = arg["Position"]
snapped = gui_utils.get_3d_view().getObjectInfo((cursor[0],
cursor[1]))
if snapped:
obj = App.ActiveDocument.getObject(snapped['Object'])
Gui.Selection.addSelection(obj)
App.activeDraftCommand.component = snapped['Component']
App.activeDraftCommand.proceed()
elif not arg["CtrlDown"] and Gui.Selection.hasSelection():
App.activeDraftCommand.proceed()
selectObject = select_object

2
src/Mod/Draft/draftobjects/patharray.py
View File

@@ -328,7 +328,7 @@ class PathArray(DraftLink):
for n in edgeNames:
e = sub[0].Shape.getElement(n)
sl.append(e)
return Part.Wire(sl)
return Part.Wire(Part.__sortEdges__(sl))
def onChanged(self, obj, prop):
"""Execute when a property is changed."""

39
src/Mod/Draft/draftutils/groups.py
View File

@@ -41,6 +41,34 @@ from draftutils.translate import _tr
from draftutils.messages import _msg, _err
def is_group(obj):
"""Return True if the given object is considered a group.
Parameters
----------
obj : App::DocumentObject
The object to check.
Returns
-------
bool
Returns `True` if `obj` is considered a group:
The object is derived from `App::DocumentObjectGroup` but not
a `'LayerContainer'`.
Or the object is of the type `'Project'`, `'Site'`, `'Building'`,
`'Floor'` or `'BuildingPart'` from the Arch Workbench.
Otherwise returns `False`.
"""
return ((obj.isDerivedFrom("App::DocumentObjectGroup")
and obj.Name != "LayerContainer")
or utils.get_type(obj) in ("Project", "Site", "Building",
"Floor", "BuildingPart"))
def get_group_names(doc=None):
"""Return a list of names of existing groups in the document.
@@ -54,12 +82,10 @@ def get_group_names(doc=None):
Returns
-------
list of str
A list of names of objects that are "groups".
These are objects derived from `App::DocumentObjectGroup`
or which are of types `'Floor'`, `'Building'`, or `'Site'`
from the Arch Workbench.
A list of names of objects that are considered groups.
See the is_group function.
Otherwise, return an empty list.
Otherwise returns an empty list.
"""
if not doc:
found, doc = utils.find_doc(App.activeDocument())
@@ -71,8 +97,7 @@ def get_group_names(doc=None):
glist = []
for obj in doc.Objects:
if (obj.isDerivedFrom("App::DocumentObjectGroup")
or utils.get_type(obj) in ("Floor", "Building", "Site")):
if is_group(obj):
glist.append(obj.Name)
return glist

4
src/Mod/Draft/draftutils/gui_utils.py
View File

@@ -64,6 +64,10 @@ def get_3d_view():
Return `None` if the graphical interface is not available.
"""
if App.GuiUp:
# FIXME The following two imports were added as part of PR4926
# Also see discussion https://forum.freecadweb.org/viewtopic.php?f=3&t=60251
import FreeCADGui as Gui
from pivy import coin
if Gui.ActiveDocument:
v = Gui.ActiveDocument.ActiveView
if "View3DInventor" in str(type(v)):

5
src/Mod/Draft/importDWG.py
View File

@@ -203,14 +203,15 @@ def convertToDxf(dwgfilename):
"""
import os, tempfile, subprocess, sys
import shutil
if shutil.which("dwg2dxf"):
try:
outdir = tempfile.mkdtemp()
basename = os.path.basename(dwgfilename)
result = outdir + os.sep + os.path.splitext(basename)[0] + ".dxf"
proc = subprocess.Popen(("dwg2dxf", dwgfilename, "-o", result))
proc.communicate()
return result
except:
pass
teigha = getTeighaConverter()
if teigha:

8
src/Mod/Fem/App/FemPostPipeline.cpp
View File

@@ -60,11 +60,11 @@ const char* FemPostPipeline::ModeEnums[]= {"Serial","Parallel",NULL};
FemPostPipeline::FemPostPipeline()
{
ADD_PROPERTY_TYPE(Filter, (0), "Pipeline", App::Prop_None, "The filter used in in this pipeline");
ADD_PROPERTY_TYPE(Filter, (0), "Pipeline", App::Prop_None, "The filter used in this pipeline");
ADD_PROPERTY_TYPE(Functions, (0), "Pipeline", App::Prop_Hidden, "The function provider which groups all pipeline functions");
ADD_PROPERTY_TYPE(Mode,(long(0)), "Pipeline", App::Prop_None, "Selects the pipeline data transition mode. In serial every filter"
"gets the output of the previous one as input, in parallel every"
"filter gets the pipelien source as input.");
ADD_PROPERTY_TYPE(Mode,(long(0)), "Pipeline", App::Prop_None, "Selects the pipeline data transition mode. In serial, every filter"
"gets the output of the previous one as input. In parallel, every"
"filter gets the pipeline source as input.");
Mode.setEnums(ModeEnums);
}

86
src/Mod/Fem/CMakeLists.txt
View File

@@ -52,6 +52,8 @@ SET(FemExamples_SRCS
femexamples/ccx_cantilever_base_face.py
femexamples/ccx_cantilever_base_solid.py
femexamples/ccx_cantilever_ele_hexa20.py
femexamples/ccx_cantilever_ele_quad4.py
femexamples/ccx_cantilever_ele_quad8.py
femexamples/ccx_cantilever_ele_tetra4.py
femexamples/ccx_cantilever_ele_tria3.py
femexamples/ccx_cantilever_ele_tria6.py
@@ -78,6 +80,7 @@ SET(FemExamples_SRCS
femexamples/material_multiple_bendingbeam_fivefaces.py
femexamples/material_multiple_tensionrod_twoboxes.py
femexamples/material_nl_platewithhole.py
femexamples/mystran_plate.py
femexamples/rc_wall_2d.py
femexamples/square_pipe_end_twisted_edgeforces.py
femexamples/square_pipe_end_twisted_nodeforces.py
@@ -94,6 +97,8 @@ SET(FemExampleMeshes_SRCS
femexamples/meshes/mesh_buckling_ibeam_tria6.py
femexamples/meshes/mesh_buckling_plate_tria6.py
femexamples/meshes/mesh_canticcx_hexa20.py
femexamples/meshes/mesh_canticcx_quad4.py
femexamples/meshes/mesh_canticcx_quad8.py
femexamples/meshes/mesh_canticcx_seg2.py
femexamples/meshes/mesh_canticcx_seg3.py
femexamples/meshes/mesh_canticcx_tetra10.py
@@ -109,6 +114,7 @@ SET(FemExampleMeshes_SRCS
femexamples/meshes/mesh_flexural_buckling.py
femexamples/meshes/mesh_multibodybeam_tetra10.py
femexamples/meshes/mesh_multibodybeam_tria6.py
femexamples/meshes/mesh_plate_mystran_quad4.py
femexamples/meshes/mesh_platewithhole_tetra10.py
femexamples/meshes/mesh_rc_wall_2d_tria6.py
femexamples/meshes/mesh_section_print_tetra10.py
@@ -124,6 +130,7 @@ SET(FemExampleMeshes_SRCS
SET(FemInOut_SRCS
feminout/__init__.py
feminout/convert2TetGen.py
feminout/exportNastranMesh.py
feminout/importCcxDatResults.py
feminout/importCcxFrdResults.py
feminout/importFenicsMesh.py
@@ -250,6 +257,19 @@ SET(FemSolverFenics_SRCS
femsolver/fenics/fenics_tools.py
)
SET(FemSolverMystran_SRCS
femsolver/mystran/__init__.py
femsolver/mystran/add_con_fixed.py
femsolver/mystran/add_con_force.py
femsolver/mystran/add_femelement_geometry.py
femsolver/mystran/add_femelement_material.py
femsolver/mystran/add_mesh.py
femsolver/mystran/add_solver_control.py
femsolver/mystran/solver.py
femsolver/mystran/tasks.py
femsolver/mystran/writer.py
)
SET(FemSolverZ88_SRCS
femsolver/z88/__init__.py
femsolver/z88/solver.py
@@ -276,6 +296,7 @@ SET(FemTestsApp_SRCS
femtest/app/test_result.py
femtest/app/test_solver_calculix.py
femtest/app/test_solver_elmer.py
femtest/app/test_solver_mystran.py
femtest/app/test_solver_z88.py
)
@@ -300,6 +321,12 @@ SET(FemTestsCcx_SRCS
femtest/data/calculix/ccx_buckling_flexuralbuckling.inp
femtest/data/calculix/ccx_cantilever_faceload.inp
femtest/data/calculix/ccx_cantilever_ele_hexa20.inp
femtest/data/calculix/ccx_cantilever_ele_quad4.inp
femtest/data/calculix/ccx_cantilever_ele_quad8.inp
femtest/data/calculix/ccx_cantilever_ele_seg2.inp
femtest/data/calculix/ccx_cantilever_ele_seg3.inp
femtest/data/calculix/ccx_cantilever_ele_tria3.inp
femtest/data/calculix/ccx_cantilever_ele_tria6.inp
femtest/data/calculix/ccx_cantilever_nodeload.inp
femtest/data/calculix/ccx_cantilever_prescribeddisplacement.inp
femtest/data/calculix/constraint_centrif.inp
@@ -353,6 +380,16 @@ SET(FemTestsMesh_SRCS
femtest/data/mesh/tetra10_mesh.z88
)
SET(FemTestsMystran_SRCS
femtest/data/mystran/__init__.py
femtest/data/mystran/ccx_cantilever_ele_quad4.bdf
femtest/data/mystran/ccx_cantilever_ele_seg2.bdf
femtest/data/mystran/ccx_cantilever_ele_tria3.bdf
femtest/data/mystran/ccx_cantilever_faceload.bdf
femtest/data/mystran/ccx_cantilever_nodeload.bdf
femtest/data/mystran/mystran_plate.bdf
)
SET(FemTestsOpen_SRCS
femtest/data/open/__init__.py
femtest/data/open/all_objects_de9b3fb438.FCStd
@@ -362,6 +399,32 @@ SET(FemTestsZ88Main_SRCS
femtest/data/z88/__init__.py
)
SET(FemTestsZ88CcxcantiEleHex20_SRCS
femtest/data/z88/__init__.py
femtest/data/z88/ccx_cantilever_ele_hexa20/51.txt
femtest/data/z88/ccx_cantilever_ele_hexa20/z88.dyn
femtest/data/z88/ccx_cantilever_ele_hexa20/z88elp.txt
femtest/data/z88/ccx_cantilever_ele_hexa20/z88i1.txt
femtest/data/z88/ccx_cantilever_ele_hexa20/z88i2.txt
femtest/data/z88/ccx_cantilever_ele_hexa20/z88i5.txt
femtest/data/z88/ccx_cantilever_ele_hexa20/z88int.txt
femtest/data/z88/ccx_cantilever_ele_hexa20/z88man.txt
femtest/data/z88/ccx_cantilever_ele_hexa20/z88mat.txt
)
SET(FemTestsZ88CcxcantiEleTria6_SRCS
femtest/data/z88/__init__.py
femtest/data/z88/ccx_cantilever_ele_tria6/51.txt
femtest/data/z88/ccx_cantilever_ele_tria6/z88.dyn
femtest/data/z88/ccx_cantilever_ele_tria6/z88elp.txt
femtest/data/z88/ccx_cantilever_ele_tria6/z88i1.txt
femtest/data/z88/ccx_cantilever_ele_tria6/z88i2.txt
femtest/data/z88/ccx_cantilever_ele_tria6/z88i5.txt
femtest/data/z88/ccx_cantilever_ele_tria6/z88int.txt
femtest/data/z88/ccx_cantilever_ele_tria6/z88man.txt
femtest/data/z88/ccx_cantilever_ele_tria6/z88mat.txt
)
SET(FemTestsZ88Ccxcantifl_SRCS
femtest/data/z88/__init__.py
femtest/data/z88/ccx_cantilever_faceload/51.txt
@@ -375,19 +438,6 @@ SET(FemTestsZ88Ccxcantifl_SRCS
femtest/data/z88/ccx_cantilever_faceload/z88mat.txt
)
SET(FemTestsZ88Ccxcantihex_SRCS
femtest/data/z88/__init__.py
femtest/data/z88/ccx_cantilever_ele_hexa20/51.txt
femtest/data/z88/ccx_cantilever_ele_hexa20/z88.dyn
femtest/data/z88/ccx_cantilever_ele_hexa20/z88elp.txt
femtest/data/z88/ccx_cantilever_ele_hexa20/z88i1.txt
femtest/data/z88/ccx_cantilever_ele_hexa20/z88i2.txt
femtest/data/z88/ccx_cantilever_ele_hexa20/z88i5.txt
femtest/data/z88/ccx_cantilever_ele_hexa20/z88int.txt
femtest/data/z88/ccx_cantilever_ele_hexa20/z88man.txt
femtest/data/z88/ccx_cantilever_ele_hexa20/z88mat.txt
)
SET(FemTestsZ88Ccxcantinl_SRCS
femtest/data/z88/__init__.py
femtest/data/z88/ccx_cantilever_nodeload/51.txt
@@ -428,6 +478,7 @@ SET(FemAllScripts
${FemSolverElmer_SRCS}
${FemSolverElmerEquations_SRCS}
${FemSolverFenics_SRCS}
${FemSolverMystran_SRCS}
${FemSolverZ88_SRCS}
${FemTests_SRCS}
${FemTestsApp_SRCS}
@@ -435,10 +486,12 @@ SET(FemAllScripts
${FemTestsCcx_SRCS}
${FemTestsElmer_SRCS}
${FemTestsMesh_SRCS}
${FemTestsMystran_SRCS}
${FemTestsOpen_SRCS}
${FemTestsZ88Main_SRCS}
${FemTestsZ88Ccxcantifl_SRCS}
${FemTestsZ88Ccxcantihex_SRCS}
${FemTestsZ88CcxcantiEleHex20_SRCS}
${FemTestsZ88CcxcantiEleTria6_SRCS}
${FemTestsZ88Ccxcantinl_SRCS}
${FemTools_SRCS}
)
@@ -464,6 +517,7 @@ INSTALL(FILES ${FemSolverCalculix_SRCS} DESTINATION Mod/Fem/femsolver/calculix)
INSTALL(FILES ${FemSolverElmer_SRCS} DESTINATION Mod/Fem/femsolver/elmer)
INSTALL(FILES ${FemSolverElmerEquations_SRCS} DESTINATION Mod/Fem/femsolver/elmer/equations)
INSTALL(FILES ${FemSolverFenics_SRCS} DESTINATION Mod/Fem/femsolver/fenics)
INSTALL(FILES ${FemSolverMystran_SRCS} DESTINATION Mod/Fem/femsolver/mystran)
INSTALL(FILES ${FemSolverZ88_SRCS} DESTINATION Mod/Fem/femsolver/z88)
INSTALL(FILES ${FemTests_SRCS} DESTINATION Mod/Fem/femtest)
INSTALL(FILES ${FemTestsApp_SRCS} DESTINATION Mod/Fem/femtest/app)
@@ -471,10 +525,12 @@ INSTALL(FILES ${FemTestsFiles_SRCS} DESTINATION Mod/Fem/femtest/data)
INSTALL(FILES ${FemTestsCcx_SRCS} DESTINATION Mod/Fem/femtest/data/calculix)
INSTALL(FILES ${FemTestsElmer_SRCS} DESTINATION Mod/Fem/femtest/data/elmer)
INSTALL(FILES ${FemTestsMesh_SRCS} DESTINATION Mod/Fem/femtest/data/mesh)
INSTALL(FILES ${FemTestsMystran_SRCS} DESTINATION Mod/Fem/femtest/data/mystran)
INSTALL(FILES ${FemTestsOpen_SRCS} DESTINATION Mod/Fem/femtest/data/open)
INSTALL(FILES ${FemTestsZ88Main_SRCS} DESTINATION Mod/Fem/femtest/data/z88)
INSTALL(FILES ${FemTestsZ88CcxcantiEleHex20_SRCS} DESTINATION Mod/Fem/femtest/data/z88/ccx_cantilever_ele_hexa20)
INSTALL(FILES ${FemTestsZ88CcxcantiEleTria6_SRCS} DESTINATION Mod/Fem/femtest/data/z88/ccx_cantilever_ele_tria6)
INSTALL(FILES ${FemTestsZ88Ccxcantifl_SRCS} DESTINATION Mod/Fem/femtest/data/z88/ccx_cantilever_faceload)
INSTALL(FILES ${FemTestsZ88Ccxcantihex_SRCS} DESTINATION Mod/Fem/femtest/data/z88/ccx_cantilever_ele_hexa20)
INSTALL(FILES ${FemTestsZ88Ccxcantinl_SRCS} DESTINATION Mod/Fem/femtest/data/z88/ccx_cantilever_nodeload)
INSTALL(FILES ${FemTools_SRCS} DESTINATION Mod/Fem/femtools)

2
src/Mod/Fem/Gui/AppFemGui.cpp
View File

@@ -40,6 +40,7 @@
#include "DlgSettingsFemGmshImp.h"
#include "DlgSettingsFemInOutVtkImp.h"
#include "DlgSettingsFemMaterialImp.h"
#include "DlgSettingsFemMystranImp.h"
#include "DlgSettingsFemZ88Imp.h"
#include "ViewProviderFemMesh.h"
#include "ViewProviderFemMeshShape.h"
@@ -169,6 +170,7 @@ PyMOD_INIT_FUNC(FemGui)
new Gui::PrefPageProducer<FemGui::DlgSettingsFemGmshImp> (QT_TRANSLATE_NOOP("QObject","FEM"));
new Gui::PrefPageProducer<FemGui::DlgSettingsFemCcxImp> (QT_TRANSLATE_NOOP("QObject","FEM"));
new Gui::PrefPageProducer<FemGui::DlgSettingsFemElmerImp> (QT_TRANSLATE_NOOP("QObject","FEM"));
new Gui::PrefPageProducer<FemGui::DlgSettingsFemMystranImp> (QT_TRANSLATE_NOOP("QObject","FEM"));
new Gui::PrefPageProducer<FemGui::DlgSettingsFemZ88Imp> (QT_TRANSLATE_NOOP("QObject","FEM"));
new Gui::PrefPageProducer<FemGui::DlgSettingsFemMaterialImp> (QT_TRANSLATE_NOOP("QObject","FEM"));

4
src/Mod/Fem/Gui/CMakeLists.txt
View File

@@ -63,6 +63,7 @@ set(FemGui_UIC_SRCS
DlgSettingsFemGmsh.ui
DlgSettingsFemInOutVtk.ui
DlgSettingsFemMaterial.ui
DlgSettingsFemMystran.ui
DlgSettingsFemZ88.ui
TaskCreateNodeSet.ui
TaskObjectName.ui
@@ -129,6 +130,9 @@ SET(FemGui_DLG_SRCS
DlgSettingsFemMaterial.ui
DlgSettingsFemMaterialImp.cpp
DlgSettingsFemMaterialImp.h
DlgSettingsFemMystran.ui
DlgSettingsFemMystranImp.cpp
DlgSettingsFemMystranImp.h
DlgSettingsFemZ88.ui
DlgSettingsFemZ88Imp.cpp
DlgSettingsFemZ88Imp.h

253
src/Mod/Fem/Gui/DlgSettingsFemMystran.ui Normal file
View File

@@ -0,0 +1,253 @@
<?xml version="1.0" encoding="UTF-8"?>
<ui version="4.0">
<class>FemGui::DlgSettingsFemMystranImp</class>
<widget class="QWidget" name="FemGui::DlgSettingsFemMystranImp">
<property name="geometry">
<rect>
<x>0</x>
<y>0</y>
<width>418</width>
<height>144</height>
</rect>
</property>
<property name="windowTitle">
<string>Mystran</string>
</property>
<layout class="QVBoxLayout" name="verticalLayout">
<item>
<layout class="QVBoxLayout" name="verticalLayout_2">
<item>
<widget class="QGroupBox" name="gb_mystran_param">
<property name="sizePolicy">
<sizepolicy hsizetype="Preferred" vsizetype="Preferred">
<horstretch>0</horstretch>
<verstretch>0</verstretch>
</sizepolicy>
</property>
<property name="layoutDirection">
<enum>Qt::LeftToRight</enum>
</property>
<property name="title">
<string>Mystran binary</string>
</property>
<property name="alignment">
<set>Qt::AlignLeading|Qt::AlignLeft|Qt::AlignTop</set>
</property>
<layout class="QHBoxLayout" name="horizontalLayout">
<property name="sizeConstraint">
<enum>QLayout::SetNoConstraint</enum>
</property>
<item>
<layout class="QGridLayout" name="gl_mystran">
<item row="0" column="2">
<widget class="Gui::PrefCheckBox" name="cb_mystran_binary_std">
<property name="text">
<string>Search in known binary directories</string>
</property>
<property name="checked">
<bool>true</bool>
</property>
<property name="prefEntry" stdset="0">
<cstring>UseStandardMystranLocation</cstring>
</property>
<property name="prefPath" stdset="0">
<cstring>Mod/Fem/mystran</cstring>
</property>
</widget>
</item>
<item row="0" column="0">
<widget class="QLabel" name="l_mystran_binary_std">
<property name="text">
<string>Mystran</string>
</property>
</widget>
</item>
<item row="2" column="0">
<widget class="QLabel" name="l_mystran_binary_path">
<property name="enabled">
<bool>false</bool>
</property>
<property name="minimumSize">
<size>
<width>100</width>
<height>0</height>
</size>
</property>
<property name="text">
<string>Mystran binary path</string>
</property>
</widget>
</item>
<item row="2" column="2">
<widget class="Gui::PrefFileChooser" name="fc_mystran_binary_path" native="true">
<property name="enabled">
<bool>false</bool>
</property>
<property name="sizePolicy">
<sizepolicy hsizetype="Preferred" vsizetype="Preferred">
<horstretch>0</horstretch>
<verstretch>0</verstretch>
</sizepolicy>
</property>
<property name="minimumSize">
<size>
<width>0</width>
<height>0</height>
</size>
</property>
<property name="sizeIncrement">
<size>
<width>0</width>
<height>0</height>
</size>
</property>
<property name="baseSize">
<size>
<width>0</width>
<height>0</height>
</size>
</property>
<property name="toolTip">
<string>Leave blank to use default Mystran binary file location</string>
</property>
<property name="prefEntry" stdset="0">
<cstring>MystranBinaryPath</cstring>
</property>
<property name="prefPath" stdset="0">
<cstring>Mod/Fem/mystran</cstring>
</property>
</widget>
</item>
<item row="4" column="0">
<widget class="QLabel" name="l_mystran_write_comments">
<property name="text">
<string>Comments</string>
</property>
</widget>
</item>
<item row="4" column="2">
<widget class="Gui::PrefCheckBox" name="cb_mystran_write_comments">
<property name="text">
<string>write comments to input file</string>
</property>
<property name="checked">
<bool>true</bool>
</property>
<property name="prefEntry" stdset="0">
<cstring>writeCommentsToInputFile</cstring>
</property>
<property name="prefPath" stdset="0">
<cstring>Mod/Fem/mystran</cstring>
</property>
</widget>
</item>
</layout>
</item>
</layout>
</widget>
</item>
</layout>
</item>
<item>
<spacer name="verticalSpacer">
<property name="orientation">
<enum>Qt::Vertical</enum>
</property>
<property name="sizeHint" stdset="0">
<size>
<width>20</width>
<height>40</height>
</size>
</property>
</spacer>
</item>
</layout>
</widget>
<layoutdefault spacing="6" margin="11"/>
<customwidgets>
<customwidget>
<class>Gui::PrefCheckBox</class>
<extends>QCheckBox</extends>
<header>Gui/PrefWidgets.h</header>
</customwidget>
<customwidget>
<class>Gui::FileChooser</class>
<extends>QWidget</extends>
<header>Gui/FileDialog.h</header>
</customwidget>
<customwidget>
<class>Gui::PrefFileChooser</class>
<extends>Gui::FileChooser</extends>
<header>Gui/PrefWidgets.h</header>
</customwidget>
</customwidgets>
<resources>
<include location="Resources/Fem.qrc"/>
</resources>
<connections>
<connection>
<sender>cb_mystran_binary_std</sender>
<signal>toggled(bool)</signal>
<receiver>l_mystran_binary_path</receiver>
<slot>setEnabled(bool)</slot>
<hints>
<hint type="sourcelabel">
<x>406</x>
<y>45</y>
</hint>
<hint type="destinationlabel">
<x>148</x>
<y>68</y>
</hint>
</hints>
</connection>
<connection>
<sender>cb_mystran_binary_std</sender>
<signal>toggled(bool)</signal>
<receiver>fc_mystran_binary_path</receiver>
<slot>setEnabled(bool)</slot>
<hints>
<hint type="sourcelabel">
<x>406</x>
<y>45</y>
</hint>
<hint type="destinationlabel">
<x>406</x>
<y>68</y>
</hint>
</hints>
</connection>
<connection>
<sender>cb_mystran_binary_std</sender>
<signal>toggled(bool)</signal>
<receiver>l_mystran_binary_path</receiver>
<slot>setDisabled(bool)</slot>
<hints>
<hint type="sourcelabel">
<x>406</x>
<y>45</y>
</hint>
<hint type="destinationlabel">
<x>148</x>
<y>68</y>
</hint>
</hints>
</connection>
<connection>
<sender>cb_mystran_binary_std</sender>
<signal>toggled(bool)</signal>
<receiver>fc_mystran_binary_path</receiver>
<slot>setDisabled(bool)</slot>
<hints>
<hint type="sourcelabel">
<x>406</x>
<y>45</y>
</hint>
<hint type="destinationlabel">
<x>406</x>
<y>68</y>
</hint>
</hints>
</connection>
</connections>
</ui>

74
src/Mod/Fem/Gui/DlgSettingsFemMystranImp.cpp Normal file
View File

@@ -0,0 +1,74 @@
/***************************************************************************
* Copyright (c) 2021 FreeCAD Developers *
* Author: Bernd Hahnebach <bernd@bimstatik.ch> *
* Based on src/Mod/Fem/Gui/DlgSettingsFemGeneralImp.cpp *
* *
* This file is part of the FreeCAD CAx development system. *
* *
* This library is free software; you can redistribute it and/or *
* modify it under the terms of the GNU Library General Public *
* License as published by the Free Software Foundation; either *
* version 2 of the License, or (at your option) any later version. *
* *
* This library is distributed in the hope that it will be useful, *
* but WITHOUT ANY WARRANTY; without even the implied warranty of *
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
* GNU Library General Public License for more details. *
* *
* You should have received a copy of the GNU Library General Public *
* License along with this library; see the file COPYING.LIB. If not, *
* write to the Free Software Foundation, Inc., 59 Temple Place, *
* Suite 330, Boston, MA 02111-1307, USA *
* *
***************************************************************************/
#include "PreCompiled.h"
#include "DlgSettingsFemMystranImp.h"
#include "ui_DlgSettingsFemMystran.h"
#include <Gui/Application.h>
#include <Gui/PrefWidgets.h>
using namespace FemGui;
DlgSettingsFemMystranImp::DlgSettingsFemMystranImp( QWidget* parent )
: PreferencePage( parent )
, ui(new Ui_DlgSettingsFemMystranImp)
{
ui->setupUi(this);
}
DlgSettingsFemMystranImp::~DlgSettingsFemMystranImp()
{
// no need to delete child widgets, Qt does it all for us
}
void DlgSettingsFemMystranImp::saveSettings()
{
ui->cb_mystran_binary_std->onSave();
ui->fc_mystran_binary_path->onSave();
ui->cb_mystran_write_comments->onSave();
}
void DlgSettingsFemMystranImp::loadSettings()
{
ui->cb_mystran_binary_std->onRestore();
ui->fc_mystran_binary_path->onRestore();
ui->cb_mystran_write_comments->onRestore();
}
/**
* Sets the strings of the subwidgets using the current language.
*/
void DlgSettingsFemMystranImp::changeEvent(QEvent *e)
{
if (e->type() == QEvent::LanguageChange) {
ui->retranslateUi(this);
}
else {
QWidget::changeEvent(e);
}
}
#include "moc_DlgSettingsFemMystranImp.cpp"

54
src/Mod/Fem/Gui/DlgSettingsFemMystranImp.h Normal file
View File

@@ -0,0 +1,54 @@
/**************************************************************************
* Copyright (c) 2021 FreeCAD Developers *
* Author: Bernd Hahnebach <bernd@bimstatik.ch> *
* Based on src/Mod/Fem/Gui/DlgSettingsFemGeneralImp.h *
* *
* This file is part of the FreeCAD CAx development system. *
* *
* This library is free software; you can redistribute it and/or *
* modify it under the terms of the GNU Library General Public *
* License as published by the Free Software Foundation; either *
* version 2 of the License, or (at your option) any later version. *
* *
* This library is distributed in the hope that it will be useful, *
* but WITHOUT ANY WARRANTY; without even the implied warranty of *
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
* GNU Library General Public License for more details. *
* *
* You should have received a copy of the GNU Library General Public *
* License along with this library; see the file COPYING.LIB. If not, *
* write to the Free Software Foundation, Inc., 59 Temple Place, *
* Suite 330, Boston, MA 02111-1307, USA *
* *
***************************************************************************/
#ifndef FEMGUI_DLGSETTINGSFEMMYSTRANIMP_H
#define FEMGUI_DLGSETTINGSFEMMYSTRANIMP_H
#include <Gui/PropertyPage.h>
#include <memory>
namespace FemGui {
class Ui_DlgSettingsFemMystranImp;
class DlgSettingsFemMystranImp : public Gui::Dialog::PreferencePage
{
Q_OBJECT
public:
DlgSettingsFemMystranImp( QWidget* parent = 0 );
~DlgSettingsFemMystranImp();
protected:
void saveSettings();
void loadSettings();
void changeEvent(QEvent *e);
private:
std::unique_ptr<Ui_DlgSettingsFemMystranImp> ui;
};
} // namespace FemGui
#endif // FEMGUI_DLGSETTINGSFEMMYSTRANIMP_H

1
src/Mod/Fem/Gui/Workbench.cpp
View File

@@ -337,6 +337,7 @@ Gui::MenuItem* Workbench::setupMenuBar() const
<< "FEM_SolverCalculixCxxtools"
<< "FEM_SolverCalculiX"
<< "FEM_SolverElmer"
<< "FEM_SolverMystran"
<< "FEM_SolverZ88"
<< "Separator"
<< "FEM_EquationElasticity"

2
src/Mod/Fem/Init.py
View File

@@ -66,6 +66,8 @@ FreeCAD.addExportType("FEM mesh TetGen (*.poly)", "feminout.convert2TetGen")
FreeCAD.addImportType("FEM mesh formats (*.bdf *.dat *.inp *.med *.unv *.vtk *.vtu *.z88)", "Fem")
FreeCAD.addExportType("FEM mesh formats (*.dat *.inp *.med *.stl *.unv *.vtk *.vtu *.z88)", "Fem")
FreeCAD.addExportType("FEM mesh Nastran (*.bdf)", "feminout.exportNastranMesh")
FreeCAD.addImportType("FEM result CalculiX (*.frd)", "feminout.importCcxFrdResults")
FreeCAD.addImportType("FEM mesh Fenics (*.xml *.xdmf)", "feminout.importFenicsMesh")

11
src/Mod/Fem/ObjectsFem.py
View File

@@ -799,6 +799,17 @@ def makeSolverElmer(
return obj
def makeSolverMystran(
doc,
name="SolverMystran"
):
"""makeSolverMystran(document, [name]):
makes a Mystran solver object"""
import femsolver.mystran.solver
obj = femsolver.mystran.solver.create(doc, name)
return obj
def makeSolverZ88(
doc,
name="SolverZ88"

22
src/Mod/Fem/femcommands/commands.py
View File

@@ -741,9 +741,9 @@ class _SolverCalculix(CommandManager):
def __init__(self):
super(_SolverCalculix, self).__init__()
self.pixmap = "FEM_SolverStandard"
self.menuetext = "Solver CalculiX (experimental)"
self.menuetext = "Solver CalculiX (new framework)"
self.accel = "S, C"
self.tooltip = "Creates a FEM solver CalculiX (experimental)"
self.tooltip = "Creates a FEM solver CalculiX new framework (less result error handling)"
self.is_active = "with_analysis"
self.is_active = "with_analysis"
self.do_activated = "add_obj_on_gui_noset_edit"
@@ -775,6 +775,19 @@ class _SolverElmer(CommandManager):
self.do_activated = "add_obj_on_gui_noset_edit"
class _SolverMystran(CommandManager):
"The FEM_SolverMystran command definition"
def __init__(self):
super(_SolverMystran, self).__init__()
self.pixmap = "FEM_SolverStandard"
self.menuetext = "Solver Mystran"
self.accel = "S, M"
self.tooltip = "Creates a FEM solver Mystran"
self.is_active = "with_analysis"
self.do_activated = "add_obj_on_gui_noset_edit"
class _SolverRun(CommandManager):
"The FEM_SolverRun command definition"
@@ -797,7 +810,6 @@ class _SolverZ88(CommandManager):
def __init__(self):
super(_SolverZ88, self).__init__()
self.pixmap = "FEM_SolverZ88.svg"
self.menuetext = "Solver Z88"
self.accel = "S, Z"
self.tooltip = "Creates a FEM solver Z88"
@@ -974,6 +986,10 @@ FreeCADGui.addCommand(
"FEM_SolverElmer",
_SolverElmer()
)
FreeCADGui.addCommand(
"FEM_SolverMystran",
_SolverMystran()
)
FreeCADGui.addCommand(
"FEM_SolverRun",
_SolverRun()

2
src/Mod/Fem/femexamples/ccx_cantilever_base_edge.py
View File

@@ -71,6 +71,8 @@ def setup_cantilever_base_edge(doc=None, solvertype="ccxtools"):
elif solvertype == "ccxtools":
solver_obj = ObjectsFem.makeSolverCalculixCcxTools(doc, "CalculiXccxTools")
solver_obj.WorkingDir = u""
elif solvertype == "mystran":
solver_obj = ObjectsFem.makeSolverMystran(doc, "SolverMystran")
else:
FreeCAD.Console.PrintWarning(
"Not known or not supported solver type: {}. "

4
src/Mod/Fem/femexamples/ccx_cantilever_base_face.py
View File

@@ -61,6 +61,10 @@ def setup_cantilever_base_face(doc=None, solvertype="ccxtools"):
elif solvertype == "ccxtools":
solver_obj = ObjectsFem.makeSolverCalculixCcxTools(doc, "CalculiXccxTools")
solver_obj.WorkingDir = u""
elif solvertype == "mystran":
solver_obj = ObjectsFem.makeSolverMystran(doc, "SolverMystran")
elif solvertype == "z88":
solver_obj = ObjectsFem.makeSolverZ88(doc, "SolverZ88")
else:
FreeCAD.Console.PrintWarning(
"Not known or not supported solver type: {}. "

2
src/Mod/Fem/femexamples/ccx_cantilever_base_solid.py
View File

@@ -60,6 +60,8 @@ def setup_cantilever_base_solid(doc=None, solvertype="ccxtools"):
elif solvertype == "elmer":
solver_obj = ObjectsFem.makeSolverElmer(doc, "SolverElmer")
ObjectsFem.makeEquationElasticity(doc, solver_obj)
elif solvertype == "mystran":
solver_obj = ObjectsFem.makeSolverMystran(doc, "SolverMystran")
elif solvertype == "z88":
solver_obj = ObjectsFem.makeSolverZ88(doc, "SolverZ88")
else:

2
src/Mod/Fem/femexamples/ccx_cantilever_ele_hexa20.py
View File

@@ -38,7 +38,7 @@ def get_information():
"meshtype": "solid",
"meshelement": "Hexa20",
"constraints": ["fixed", "force"],
"solvers": ["calculix", "z88", "elmer"],
"solvers": ["calculix", "elmer", "z88"],
"material": "solid",
"equation": "mechanical"
}

101
src/Mod/Fem/femexamples/ccx_cantilever_ele_quad4.py Normal file
View File

@@ -0,0 +1,101 @@
# ***************************************************************************
# * Copyright (c) 2021 Bernd Hahnebach <bernd@bimstatik.org> *
# * *
# * This file is part of the FreeCAD CAx development system. *
# * *
# * This program is free software; you can redistribute it and/or modify *
# * it under the terms of the GNU Lesser General Public License (LGPL) *
# * as published by the Free Software Foundation; either version 2 of *
# * the License, or (at your option) any later version. *
# * for detail see the LICENCE text file. *
# * *
# * This program is distributed in the hope that it will be useful, *
# * but WITHOUT ANY WARRANTY; without even the implied warranty of *
# * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
# * GNU Library General Public License for more details. *
# * *
# * You should have received a copy of the GNU Library General Public *
# * License along with this program; if not, write to the Free Software *
# * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 *
# * USA *
# * *
# ***************************************************************************
import FreeCAD
import Fem
from . import manager
from .ccx_cantilever_base_face import setup_cantilever_base_face
from .manager import get_meshname
from .manager import init_doc
def get_information():
return {
"name": "CCX cantilever quad4 face elements",
"meshtype": "face",
"meshelement": "Quad4",
"constraints": ["fixed", "force"],
"solvers": ["calculix", "mystran"],
"material": "solid",
"equation": "mechanical"
}
def get_explanation(header=""):
return header + """
To run the example from Python console use:
from femexamples.ccx_cantilever_ele_quad4 import setup
setup()
See forum topic post:
CalculiX cantilever modeled with quad4 face elements
"""
def setup(doc=None, solvertype="ccxtools"):
# init FreeCAD document
if doc is None:
doc = init_doc()
# explanation object
# just keep the following line and change text string in get_explanation method
manager.add_explanation_obj(doc, get_explanation(manager.get_header(get_information())))
# setup CalculiX cantilever
doc = setup_cantilever_base_face(doc, solvertype)
femmesh_obj = doc.getObject(get_meshname())
geom_obj = doc.getObject("CanileverPlate")
# load the quad4 mesh
from .meshes.mesh_canticcx_quad4 import create_nodes, create_elements
new_fem_mesh = Fem.FemMesh()
control = create_nodes(new_fem_mesh)
if not control:
FreeCAD.Console.PrintError("Error on creating nodes.\n")
control = create_elements(new_fem_mesh)
if not control:
FreeCAD.Console.PrintError("Error on creating elements.\n")
# overwrite mesh with the quad4 mesh
femmesh_obj.FemMesh = new_fem_mesh
# set mesh obj parameter
femmesh_obj.Part = geom_obj
femmesh_obj.SecondOrderLinear = False
femmesh_obj.ElementDimension = "2D"
femmesh_obj.ElementOrder = "1st"
femmesh_obj.CharacteristicLengthMax = "2000.0 mm"
femmesh_obj.CharacteristicLengthMin = "2000.0 mm"
femmesh_obj.RecombinationAlgorithm = "Blossom full-quad"
femmesh_obj.RecombineAll = True
doc.recompute()
return doc

101
src/Mod/Fem/femexamples/ccx_cantilever_ele_quad8.py Normal file
View File

@@ -0,0 +1,101 @@
# ***************************************************************************
# * Copyright (c) 2021 Bernd Hahnebach <bernd@bimstatik.org> *
# * *
# * This file is part of the FreeCAD CAx development system. *
# * *
# * This program is free software; you can redistribute it and/or modify *
# * it under the terms of the GNU Lesser General Public License (LGPL) *
# * as published by the Free Software Foundation; either version 2 of *
# * the License, or (at your option) any later version. *
# * for detail see the LICENCE text file. *
# * *
# * This program is distributed in the hope that it will be useful, *
# * but WITHOUT ANY WARRANTY; without even the implied warranty of *
# * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
# * GNU Library General Public License for more details. *
# * *
# * You should have received a copy of the GNU Library General Public *
# * License along with this program; if not, write to the Free Software *
# * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 *
# * USA *
# * *
# ***************************************************************************
import FreeCAD
import Fem
from . import manager
from .ccx_cantilever_base_face import setup_cantilever_base_face
from .manager import get_meshname
from .manager import init_doc
def get_information():
return {
"name": "CCX cantilever quad8 face elements",
"meshtype": "face",
"meshelement": "Quad8",
"constraints": ["fixed", "force"],
"solvers": ["calculix"],
"material": "solid",
"equation": "mechanical"
}
def get_explanation(header=""):
return header + """
To run the example from Python console use:
from femexamples.ccx_cantilever_ele_quad8 import setup
setup()
See forum topic post:
CalculiX cantilever modeled with quad8 face elements
"""
def setup(doc=None, solvertype="ccxtools"):
# init FreeCAD document
if doc is None:
doc = init_doc()
# explanation object
# just keep the following line and change text string in get_explanation method
manager.add_explanation_obj(doc, get_explanation(manager.get_header(get_information())))
# setup CalculiX cantilever
doc = setup_cantilever_base_face(doc, solvertype)
femmesh_obj = doc.getObject(get_meshname())
geom_obj = doc.getObject("CanileverPlate")
# load the quad8 mesh
from .meshes.mesh_canticcx_quad8 import create_nodes, create_elements
new_fem_mesh = Fem.FemMesh()
control = create_nodes(new_fem_mesh)
if not control:
FreeCAD.Console.PrintError("Error on creating nodes.\n")
control = create_elements(new_fem_mesh)
if not control:
FreeCAD.Console.PrintError("Error on creating elements.\n")
# overwrite mesh with the quad8 mesh
femmesh_obj.FemMesh = new_fem_mesh
# set mesh obj parameter
femmesh_obj.Part = geom_obj
femmesh_obj.SecondOrderLinear = False
femmesh_obj.ElementDimension = "2D"
femmesh_obj.ElementOrder = "1st"
femmesh_obj.CharacteristicLengthMax = "4000.0 mm"
femmesh_obj.CharacteristicLengthMin = "4000.0 mm"
femmesh_obj.RecombinationAlgorithm = "Blossom full-quad"
femmesh_obj.RecombineAll = True
doc.recompute()
return doc

4
src/Mod/Fem/femexamples/ccx_cantilever_ele_seg2.py
View File

@@ -37,7 +37,7 @@ def get_information():
"meshtype": "edge",
"meshelement": "Seg2",
"constraints": ["fixed", "force"],
"solvers": ["calculix"],
"solvers": ["calculix", "mystran"],
"material": "solid",
"equation": "mechanical"
}
@@ -72,6 +72,7 @@ def setup(doc=None, solvertype="ccxtools"):
# setup CalculiX cantilever
doc = setup_cantilever_base_edge(doc, solvertype)
femmesh_obj = doc.getObject(get_meshname())
geom_obj = doc.getObject("CantileverLine")
# load the seg2 mesh
from .meshes.mesh_canticcx_seg2 import create_nodes, create_elements
@@ -87,6 +88,7 @@ def setup(doc=None, solvertype="ccxtools"):
femmesh_obj.FemMesh = new_fem_mesh
# set mesh obj parameter
femmesh_obj.Part = geom_obj
femmesh_obj.SecondOrderLinear = False
femmesh_obj.ElementDimension = "1D"
femmesh_obj.ElementOrder = "1st"

9
src/Mod/Fem/femexamples/ccx_cantilever_ele_tetra4.py
View File

@@ -36,7 +36,7 @@ def get_information():
"meshtype": "solid",
"meshelement": "Tetra4",
"constraints": ["fixed", "force"],
"solvers": ["calculix", "elmer", "z88"],
"solvers": ["calculix", "elmer", "mystran", "z88"],
"material": "solid",
"equation": "mechanical"
}
@@ -67,12 +67,19 @@ def setup(doc=None, solvertype="ccxtools"):
# just keep the following line and change text string in get_explanation method
manager.add_explanation_obj(doc, get_explanation(manager.get_header(get_information())))
# delete explanation object wrongly added with setup faceload
if hasattr(doc, "Explanation_Report001"):
doc.removeObject("Explanation_Report001")
doc.recompute()
# setup cantilever faceload and exchange the mesh
doc = setup_with_faceload(doc, solvertype)
femmesh_obj = doc.getObject(get_meshname())
geom_obj = doc.getObject("Box")
# clear mesh and set meshing parameter
femmesh_obj.FemMesh = Fem.FemMesh()
femmesh_obj.Part = geom_obj
femmesh_obj.SecondOrderLinear = False
femmesh_obj.ElementDimension = "3D"
femmesh_obj.ElementOrder = "1st"

4
src/Mod/Fem/femexamples/ccx_cantilever_ele_tria3.py
View File

@@ -37,7 +37,7 @@ def get_information():
"meshtype": "face",
"meshelement": "Tria3",
"constraints": ["fixed", "force"],
"solvers": ["calculix"],
"solvers": ["calculix", "mystran"],
"material": "solid",
"equation": "mechanical"
}
@@ -72,6 +72,7 @@ def setup(doc=None, solvertype="ccxtools"):
# setup CalculiX cantilever
doc = setup_cantilever_base_face(doc, solvertype)
femmesh_obj = doc.getObject(get_meshname())
geom_obj = doc.getObject("CanileverPlate")
# load the tria3 mesh
from .meshes.mesh_canticcx_tria3 import create_nodes, create_elements
@@ -87,6 +88,7 @@ def setup(doc=None, solvertype="ccxtools"):
femmesh_obj.FemMesh = new_fem_mesh
# set mesh obj parameter
femmesh_obj.Part = geom_obj
femmesh_obj.SecondOrderLinear = False
femmesh_obj.ElementDimension = "2D"
femmesh_obj.ElementOrder = "1st"

2
src/Mod/Fem/femexamples/ccx_cantilever_ele_tria6.py
View File

@@ -32,7 +32,7 @@ def get_information():
"meshtype": "face",
"meshelement": "Tria6",
"constraints": ["fixed", "force"],
"solvers": ["calculix"],
"solvers": ["calculix", "z88"],
"material": "solid",
"equation": "mechanical"
}

2
src/Mod/Fem/femexamples/ccx_cantilever_faceload.py
View File

@@ -35,7 +35,7 @@ def get_information():
"meshtype": "solid",
"meshelement": "Tet10",
"constraints": ["fixed", "force"],
"solvers": ["calculix", "z88", "elmer"],
"solvers": ["calculix", "elmer", "mystran", "z88"],
"material": "solid",
"equation": "mechanical"
}

2
src/Mod/Fem/femexamples/ccx_cantilever_nodeload.py
View File

@@ -35,7 +35,7 @@ def get_information():
"meshtype": "solid",
"meshelement": "Tet10",
"constraints": ["fixed", "force"],
"solvers": ["calculix", "z88", "elmer"],
"solvers": ["calculix", "elmer", "mystran", "z88"],
"material": "solid",
"equation": "mechanical"
}

51
src/Mod/Fem/femexamples/meshes/mesh_canticcx_quad4.py Normal file
View File

@@ -0,0 +1,51 @@
def create_nodes(femmesh):
# nodes
femmesh.addNode(0.0, 500.0, 0.0, 1)
femmesh.addNode(0.0, 500.0000000000002, 1000.0, 2)
femmesh.addNode(8000.0, 500.0, 0.0, 3)
femmesh.addNode(8000.0, 500.0000000000002, 1000.0, 4)
femmesh.addNode(0.0, 500.0000000000001, 500.0, 5)
femmesh.addNode(2000.0, 500.0, 0.0, 6)
femmesh.addNode(4000.0, 500.0, 0.0, 7)
femmesh.addNode(6000.0, 500.0, 0.0, 8)
femmesh.addNode(1000.0, 500.0, 0.0, 9)
femmesh.addNode(7000.0, 500.0, 0.0, 10)
femmesh.addNode(3000.0, 500.0, 0.0, 11)
femmesh.addNode(5000.0, 500.0, 0.0, 12)
femmesh.addNode(8000.0, 500.0000000000001, 500.0, 13)
femmesh.addNode(2000.0, 500.0000000000002, 1000.0, 14)
femmesh.addNode(4000.0, 500.0000000000002, 1000.0, 15)
femmesh.addNode(6000.0, 500.0000000000002, 1000.0, 16)
femmesh.addNode(1000.0, 500.0000000000002, 1000.0, 17)
femmesh.addNode(7000.0, 500.0000000000002, 1000.0, 18)
femmesh.addNode(3000.0, 500.0000000000002, 1000.0, 19)
femmesh.addNode(5000.0, 500.0000000000002, 1000.0, 20)
femmesh.addNode(2000.0, 500.0000000000001, 500.0, 21)
femmesh.addNode(1000.0, 500.0000000000001, 500.0, 22)
femmesh.addNode(6000.0, 500.0000000000001, 500.0, 23)
femmesh.addNode(7000.0, 500.0000000000001, 500.0, 24)
femmesh.addNode(4000.0, 500.0000000000001, 500.0, 25)
femmesh.addNode(3000.0, 500.0000000000001, 500.0, 26)
femmesh.addNode(5000.0, 500.0000000000001, 500.0, 27)
return True
def create_elements(femmesh):
# elements
femmesh.addFace([14, 17, 22, 21], 21)
femmesh.addFace([6, 21, 22, 9], 22)
femmesh.addFace([1, 9, 22, 5], 23)
femmesh.addFace([2, 5, 22, 17], 24)
femmesh.addFace([4, 18, 24, 13], 25)
femmesh.addFace([3, 13, 24, 10], 26)
femmesh.addFace([8, 10, 24, 23], 27)
femmesh.addFace([16, 23, 24, 18], 28)
femmesh.addFace([7, 25, 26, 11], 29)
femmesh.addFace([6, 11, 26, 21], 30)
femmesh.addFace([14, 21, 26, 19], 31)
femmesh.addFace([15, 19, 26, 25], 32)
femmesh.addFace([16, 20, 27, 23], 33)
femmesh.addFace([8, 23, 27, 12], 34)
femmesh.addFace([7, 12, 27, 25], 35)
femmesh.addFace([15, 25, 27, 20], 36)
return True

37
src/Mod/Fem/femexamples/meshes/mesh_canticcx_quad8.py Normal file
View File

@@ -0,0 +1,37 @@
def create_nodes(femmesh):
# nodes
femmesh.addNode(0.0, 500.0, 0.0, 1)
femmesh.addNode(0.0, 500.00000000000324, 1000.0, 2)
femmesh.addNode(8000.0, 500.0, 0.0, 3)
femmesh.addNode(8000.0, 500.00000000000324, 1000.0, 4)
femmesh.addNode(0.0, 500.0000000000001, 500.0, 5)
femmesh.addNode(0.0, 500.00000000000006, 250.0, 6)
femmesh.addNode(0.0, 500.00000000000017, 750.0, 7)
femmesh.addNode(4000.0, 500.0, 0.0, 8)
femmesh.addNode(2000.0, 500.0, 0.0, 9)
femmesh.addNode(6000.0, 500.0, 0.0, 10)
femmesh.addNode(8000.0, 500.0000000000001, 500.0, 11)
femmesh.addNode(8000.0, 500.00000000000006, 250.0, 12)
femmesh.addNode(8000.0, 500.00000000000017, 750.0, 13)
femmesh.addNode(4000.0, 500.0000000000002, 1000.0, 14)
femmesh.addNode(2000.0, 500.0000000000002, 1000.0, 15)
femmesh.addNode(6000.0, 500.0000000000002, 1000.0, 16)
femmesh.addNode(4000.0, 500.0000000000001, 500.0, 17)
femmesh.addNode(2000.0, 500.0000000000001, 500.0, 18)
femmesh.addNode(4000.0, 500.00000000000017, 750.0, 19)
femmesh.addNode(2000.0, 500.00000000000017, 750.0, 20)
femmesh.addNode(6000.0, 500.0000000000001, 500.0, 21)
femmesh.addNode(6000.0, 500.00000000000017, 750.0, 22)
femmesh.addNode(4000.0, 500.00000000000006, 250.0, 23)
femmesh.addNode(6000.0, 500.00000000000006, 250.0, 24)
femmesh.addNode(2000.0, 500.00000000000006, 250.0, 25)
return True
def create_elements(femmesh):
# elements
femmesh.addFace([1, 5, 17, 8, 6, 18, 23, 9], 1)
femmesh.addFace([5, 2, 14, 17, 7, 15, 19, 18], 2)
femmesh.addFace([17, 14, 4, 11, 19, 16, 13, 21], 3)
femmesh.addFace([8, 17, 11, 3, 23, 21, 12, 10], 4)
return True

69
src/Mod/Fem/femexamples/meshes/mesh_plate_mystran_quad4.py Normal file
View File

@@ -0,0 +1,69 @@
def create_nodes(femmesh):
# nodes
femmesh.addNode(0.0, 0.0, 0.0, 1)
femmesh.addNode(2.0, 0.0, 0.0, 2)
femmesh.addNode(4.0, 0.0, 0.0, 3)
femmesh.addNode(6.0, 0.0, 0.0, 4)
femmesh.addNode(8.0, 0.0, 0.0, 5)
femmesh.addNode(10.0, 0.0, 0.0, 6)
femmesh.addNode(0.0, 2.0, 0.0, 7)
femmesh.addNode(2.0, 2.0, 0.0, 8)
femmesh.addNode(4.0, 2.0, 0.0, 9)
femmesh.addNode(6.0, 2.0, 0.0, 10)
femmesh.addNode(8.0, 2.0, 0.0, 11)
femmesh.addNode(10.0, 2.0, 0.0, 12)
femmesh.addNode(0.0, 4.0, 0.0, 13)
femmesh.addNode(2.0, 4.0, 0.0, 14)
femmesh.addNode(4.0, 4.0, 0.0, 15)
femmesh.addNode(6.0, 4.0, 0.0, 16)
femmesh.addNode(8.0, 4.0, 0.0, 17)
femmesh.addNode(10.0, 4.0, 0.0, 18)
femmesh.addNode(0.0, 6.0, 0.0, 19)
femmesh.addNode(2.0, 6.0, 0.0, 20)
femmesh.addNode(4.0, 6.0, 0.0, 21)
femmesh.addNode(6.0, 6.0, 0.0, 22)
femmesh.addNode(8.0, 6.0, 0.0, 23)
femmesh.addNode(10.0, 6.0, 0.0, 24)
femmesh.addNode(0.0, 8.0, 0.0, 25)
femmesh.addNode(2.0, 8.0, 0.0, 26)
femmesh.addNode(4.0, 8.0, 0.0, 27)
femmesh.addNode(6.0, 8.0, 0.0, 28)
femmesh.addNode(8.0, 8.0, 0.0, 29)
femmesh.addNode(10.0, 8.0, 0.0, 30)
femmesh.addNode(0.0, 10.0, 0.0, 31)
femmesh.addNode(2.0, 10.0, 0.0, 32)
femmesh.addNode(4.0, 10.0, 0.0, 33)
femmesh.addNode(6.0, 10.0, 0.0, 34)
femmesh.addNode(8.0, 10.0, 0.0, 35)
femmesh.addNode(10.0, 10.0, 0.0, 36)
return True
def create_elements(femmesh):
# elements
femmesh.addFace([1, 2, 8, 7], 1)
femmesh.addFace([2, 3, 9, 8], 2)
femmesh.addFace([3, 4, 10, 9], 3)
femmesh.addFace([4, 5, 11, 10], 4)
femmesh.addFace([5, 6, 12, 11], 5)
femmesh.addFace([7, 8, 14, 13], 6)
femmesh.addFace([8, 9, 15, 14], 7)
femmesh.addFace([9, 10, 16, 15], 8)
femmesh.addFace([10, 11, 17, 16], 9)
femmesh.addFace([11, 12, 18, 17], 10)
femmesh.addFace([13, 14, 20, 19], 11)
femmesh.addFace([14, 15, 21, 20], 12)
femmesh.addFace([15, 16, 22, 21], 13)
femmesh.addFace([16, 17, 23, 22], 14)
femmesh.addFace([17, 18, 24, 23], 15)
femmesh.addFace([19, 20, 26, 25], 16)
femmesh.addFace([20, 21, 27, 26], 17)
femmesh.addFace([21, 22, 28, 27], 18)
femmesh.addFace([22, 23, 29, 28], 19)
femmesh.addFace([23, 24, 30, 29], 20)
femmesh.addFace([25, 26, 32, 31], 21)
femmesh.addFace([26, 27, 33, 32], 22)
femmesh.addFace([27, 28, 34, 33], 23)
femmesh.addFace([28, 29, 35, 34], 24)
femmesh.addFace([29, 30, 36, 35], 25)
return True

196
src/Mod/Fem/femexamples/mystran_plate.py Normal file
View File

@@ -0,0 +1,196 @@
# ***************************************************************************
# * Copyright (c) 2021 Bernd Hahnebach <bernd@bimstatik.org> *
# * *
# * This file is part of the FreeCAD CAx development system. *
# * *
# * This program is free software; you can redistribute it and/or modify *
# * it under the terms of the GNU Lesser General Public License (LGPL) *
# * as published by the Free Software Foundation; either version 2 of *
# * the License, or (at your option) any later version. *
# * for detail see the LICENCE text file. *
# * *
# * This program is distributed in the hope that it will be useful, *
# * but WITHOUT ANY WARRANTY; without even the implied warranty of *
# * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
# * GNU Library General Public License for more details. *
# * *
# * You should have received a copy of the GNU Library General Public *
# * License along with this program; if not, write to the Free Software *
# * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 *
# * USA *
# * *
# ***************************************************************************
import FreeCAD
from BOPTools import SplitFeatures
import Fem
import ObjectsFem
from . import manager
from .manager import get_meshname
from .manager import init_doc
def get_information():
return {
"name": "Mystran Plate",
"meshtype": "face",
"meshelement": "Quad4",
"constraints": ["displacement", "force"],
"solvers": ["calculix", "elmer", "mystran"],
"material": "solid",
"equation": "mechanical"
}
def get_explanation(header=""):
return header + """
To run the example from Python console use:
from femexamples.buckling_platebuckling import setup
setup()
See forum topic post:
https://forum.freecadweb.org/viewtopic.php?f=18&t=60320&start=10#p517884
This version here uses some real value for the Young's Modulus
The CalculiX steel is used
plate 10 mm x 10 mm * 0.3 mm
one each mesh node on one edge 100 N tension force
Does not work on Z88 because Z88 does not support quad4 elements
"""
def setup(doc=None, solvertype="ccxtools"):
# init FreeCAD document
if doc is None:
doc = init_doc()
# explanation object
# just keep the following line and change text string in get_explanation method
manager.add_explanation_obj(doc, get_explanation(manager.get_header(get_information())))
# geometric object
plate = doc.addObject("Part::Plane", "Plate")
plate.Width = 10
plate.Length = 10
force_pt1 = doc.addObject("Part::Vertex", "ForcePT1")
force_pt1.X = 10
force_pt1.Y = 2
force_pt2 = doc.addObject("Part::Vertex", "ForcePT2")
force_pt2.X = 10
force_pt2.Y = 4
force_pt3 = doc.addObject("Part::Vertex", "ForcePT3")
force_pt3.X = 10
force_pt3.Y = 6
force_pt4 = doc.addObject("Part::Vertex", "ForcePT4")
force_pt4.X = 10
force_pt4.Y = 8
doc.recompute()
# all geom boolean fragment
geom_obj = SplitFeatures.makeBooleanFragments(name='ThePointPlate')
geom_obj.Objects = [plate, force_pt1, force_pt2, force_pt3, force_pt4]
doc.recompute()
if FreeCAD.GuiUp:
plate.ViewObject.hide()
force_pt1.ViewObject.hide()
force_pt2.ViewObject.hide()
force_pt3.ViewObject.hide()
force_pt4.ViewObject.hide()
doc.recompute()
if FreeCAD.GuiUp:
geom_obj.ViewObject.PointSize = 10
geom_obj.ViewObject.Document.activeView().viewAxonometric()
geom_obj.ViewObject.Document.activeView().fitAll()
# analysis
analysis = ObjectsFem.makeAnalysis(doc, "Analysis")
# solver
if solvertype == "calculix":
solver_obj = ObjectsFem.makeSolverCalculix(doc, "SolverCalculiX")
elif solvertype == "ccxtools":
solver_obj = ObjectsFem.makeSolverCalculixCcxTools(doc, "CalculiXccxTools")
solver_obj.WorkingDir = u""
elif solvertype == "elmer":
solver_obj = ObjectsFem.makeSolverElmer(doc, "SolverElmer")
ObjectsFem.makeEquationElasticity(doc, solver_obj)
elif solvertype == "mystran":
solver_obj = ObjectsFem.makeSolverMystran(doc, "SolverMystran")
elif solvertype == "z88":
solver_obj = ObjectsFem.makeSolverZ88(doc, "SolverZ88")
else:
FreeCAD.Console.PrintWarning(
"Not known or not supported solver type: {}. "
"No solver object was created.\n".format(solvertype)
)
if solvertype == "calculix" or solvertype == "ccxtools":
solver_obj.SplitInputWriter = False
solver_obj.AnalysisType = "static"
solver_obj.GeometricalNonlinearity = "linear"
solver_obj.ThermoMechSteadyState = False
solver_obj.MatrixSolverType = "default"
solver_obj.IterationsControlParameterTimeUse = False
analysis.addObject(solver_obj)
# shell thickness
thickness_obj = ObjectsFem.makeElementGeometry2D(doc, 0.3, 'Thickness')
analysis.addObject(thickness_obj)
# material
material_obj = ObjectsFem.makeMaterialSolid(doc, "FemMaterial")
mat = material_obj.Material
mat["Name"] = "CalculiX-Steel"
mat["YoungsModulus"] = "210000 MPa"
mat["PoissonRatio"] = "0.30"
material_obj.Material = mat
analysis.addObject(material_obj)
# constraint fixed
con_fixed = ObjectsFem.makeConstraintFixed(doc, "ConstraintFixed")
con_fixed.References = [(geom_obj, "Edge1")]
analysis.addObject(con_fixed)
# constraint force
con_force = ObjectsFem.makeConstraintForce(doc, "ConstraintForce")
con_force.References = [
(geom_obj, "Vertex7"),
(geom_obj, "Vertex1"),
(geom_obj, "Vertex2"),
(geom_obj, "Vertex3"),
(geom_obj, "Vertex4"),
(geom_obj, "Vertex8"),
]
con_force.Force = 600 # 600 N on six nodes == 100 N/Node
con_force.Reversed = False
con_force.Direction = (geom_obj, ["Edge2"])
analysis.addObject(con_force)
# mesh
from .meshes.mesh_plate_mystran_quad4 import create_nodes, create_elements
fem_mesh = Fem.FemMesh()
control = create_nodes(fem_mesh)
if not control:
FreeCAD.Console.PrintError("Error on creating nodes.\n")
control = create_elements(fem_mesh)
if not control:
FreeCAD.Console.PrintError("Error on creating elements.\n")
femmesh_obj = analysis.addObject(ObjectsFem.makeMeshGmsh(doc, get_meshname()))[0]
femmesh_obj.FemMesh = fem_mesh
femmesh_obj.Part = geom_obj
femmesh_obj.SecondOrderLinear = False
femmesh_obj.CharacteristicLengthMax = "1.0 mm"
femmesh_obj.ElementDimension = "2D"
femmesh_obj.ElementOrder = "1st"
doc.recompute()
return doc

6
src/Mod/Fem/femguiutils/migrate_gui.py
View File

@@ -114,6 +114,8 @@ class FemMigrateGui(object):
return self
if fullname == "PyGui._ViewProviderFemSolverCalculix":
return self
if fullname == "PyGui._ViewProviderFemSolverZ88":
return self
if fullname == "PyGui._ViewProviderFemBeamSection":
return self
@@ -310,6 +312,9 @@ class FemMigrateGui(object):
import femviewprovider.view_solver_ccxtools
module._ViewProviderFemSolverCalculix = \
femviewprovider.view_solver_ccxtools.VPSolverCcxTools
if module.__name__ == "PyGui._ViewProviderFemSolverZ88":
import femsolver.z88.solver
module._ViewProviderFemSolverZ88 = femsolver.z88.solver.ViewProxy
if module.__name__ == "PyGui._ViewProviderFemBeamSection":
import femviewprovider.view_element_geometry1D
@@ -434,6 +439,7 @@ module="PyGui._ViewProviderFemMeshRegion"
module="PyGui._ViewProviderFemMeshResult"
module="PyGui._ViewProviderFemResultMechanical"
module="PyGui._ViewProviderFemSolverCalculix"
module="PyGui._ViewProviderFemSolverZ88"
renamed between the second and third big moveings
module="PyGui._ViewProviderFemBeamSection"

215
src/Mod/Fem/feminout/exportNastranMesh.py Normal file
View File

@@ -0,0 +1,215 @@
# ***************************************************************************
# * Copyright (c) 2021 Bernd Hahnebach <bernd@bimstatik.org> *
# * *
# * This file is part of the FreeCAD CAx development system. *
# * *
# * This program is free software; you can redistribute it and/or modify *
# * it under the terms of the GNU Lesser General Public License (LGPL) *
# * as published by the Free Software Foundation; either version 2 of *
# * the License, or (at your option) any later version. *
# * for detail see the LICENCE text file. *
# * *
# * This program is distributed in the hope that it will be useful, *
# * but WITHOUT ANY WARRANTY; without even the implied warranty of *
# * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
# * GNU Library General Public License for more details. *
# * *
# * You should have received a copy of the GNU Library General Public *
# * License along with this program; if not, write to the Free Software *
# * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 *
# * USA *
# * *
# ***************************************************************************
__title__ = "Mesh export for Nastran mesh file format"
__author__ = "Bernd Hahnebach"
__url__ = "https://www.freecadweb.org"
## @package exportPyNastranMesh
# \ingroup FEM
# \brief FreeCAD pyNastran Mesh writer for FEM workbench
import FreeCAD
# we need to import FreeCAD before the non FreeCAD library because of the print
try:
from pyNastran.bdf.bdf import BDF
except Exception:
FreeCAD.Console.PrintError(
"Module pyNastran not found. Writing Mystran solver input will not be work.\n"
)
from FreeCAD import Console
from femmesh import meshtools
# ************************************************************************************************
# ********* generic FreeCAD export methods *******************************************************
# names are fix given from FreeCAD, these methods are called from FreeCAD
# they are set in FEM modules Init.py
def export(
objectslist,
filename
):
"called when freecad exports a file"
if len(objectslist) != 1:
Console.PrintError("This exporter can only export one object.\n")
return
obj = objectslist[0]
if not obj.isDerivedFrom("Fem::FemMeshObject"):
Console.PrintError("No FEM mesh object selected.\n")
return
write(obj.FemMesh, filename)
# ************************************************************************************************
# ********* module specific methods **************************************************************
# writer:
# - a method directly writes a FemMesh to the mesh file
# - a method generates the pyNastran code
# ********* writer *******************************************************************************
def write(
fem_mesh,
filename
):
"""directly write a FemMesh to a pyNastran mesh file format
fem_mesh: a FemMesh"""
if not fem_mesh.isDerivedFrom("Fem::FemMesh"):
Console.PrintError("Not a FemMesh was given as parameter.\n")
return
femnodes_mesh = fem_mesh.Nodes
femelement_table = meshtools.get_femelement_table(fem_mesh)
export_element_type = get_export_element_type(fem_mesh, femelement_table)
model = BDF()
mesh_pynas_code = get_pynastran_mesh(femnodes_mesh, femelement_table, export_element_type)
mesh_pynas_code += missing_code_pnynasmesh
# pynas file
basefilename = filename[:len(filename) - 4] # TODO basename is more failsave
pynasf = open(basefilename + ".py", "w")
pynasf.write("# written by FreeCAD\n\n\n")
pynasf.write("from pyNastran.bdf.bdf import BDF\n")
pynasf.write("model = BDF()\n\n\n")
pynasf.write(mesh_pynas_code)
pynasf.write(
"model.write_bdf('{}', enddata=True)\n"
.format(basefilename + "_pyNas.bdf")
)
pynasf.close()
# execute pyNastran code to add grid to the model
# print(model)
# print(model.get_bdf_stats())
exec(mesh_pynas_code)
# print(model)
# print(model.get_bdf_stats())
# write Nastran mesh file
model.write_bdf(filename, enddata=True) # TODO FIXME "BEGIN BULK" is missing
def get_pynastran_mesh(
femnodes_mesh,
femelement_table,
export_element_type,
):
if export_element_type is None:
Console.PrintError("Error: wrong export_element_type.\n")
return
# nodes
pynas_nodes = "# grid cards, geometric mesh points\n"
for node in femnodes_mesh:
vec = femnodes_mesh[node]
pynas_nodes += "model.add_grid({}, [{}, {}, {}])\n".format(node, vec.x, vec.y, vec.z)
# print(pynas_nodes)
# elements
# Nastran seams to have the same node order as SMESH (FreeCAD) has
# thus just write the nodes at once
pynas_elements = "# elements cards\n"
for element in femelement_table:
nodes = femelement_table[element]
# print(element) # eleid
# print(n) # tuple of nodes
if export_element_type == "cbar":
pynas_elements += (
"model.add_{ele_keyword}({eid}, {pid}, {nodes}, "
"{orientation_vec}, {gnull})\n"
.format(
ele_keyword=export_element_type,
eid=element,
pid=1,
nodes=nodes,
orientation_vec="x=[0.0, 0.0, 1.0]",
gnull="g0=None"
)
)
else:
if export_element_type == "ctetra4":
ele_keyword = "ctetra"
# N1, N3, N2, N4
the_nodes = [nodes[0], nodes[2], nodes[1], nodes[3]]
elif export_element_type == "ctetra10":
ele_keyword = "ctetra"
# N1, N3, N2, N4, N7, N6, N5, N8, N10, N9
the_nodes = [
nodes[0], nodes[2], nodes[1], nodes[3],
nodes[6], nodes[5], nodes[4],
nodes[7], nodes[9], nodes[8],
]
else:
ele_keyword = export_element_type
the_nodes = nodes
pynas_elements += (
"model.add_{ele_keyword}({eid}, {pid}, {nodes})\n"
.format(ele_keyword=ele_keyword, eid=element, pid=1, nodes=the_nodes)
)
# print(pynas_elements)
mesh_pynas_code = "{}\n\n{}\n\n".format(pynas_nodes, pynas_elements)
return mesh_pynas_code
# Helper
def get_export_element_type(
femmesh,
femelement_table=None
):
return nastran_ele_types[meshtools.get_femmesh_eletype(femmesh, femelement_table)]
nastran_ele_types = {
"tetra4": "ctetra4",
"tetra10": "ctetra10",
"hexa8": None,
"hexa20": None,
"tria3": "ctria3",
"tria6": None,
"quad4": "cquad4",
"quad8": None,
"seg2": "cbar",
"seg3": None,
"None": None,
}
missing_code_pnynasmesh = """
model.sol = 101 # is this needed?
# case control
from pyNastran.bdf.bdf import CaseControlDeck
cc = CaseControlDeck([
#"ECHO = NONE",
"TITLE = pyNastran for generating solverinput for for Mystran",
#"SUBCASE 1",
])
model.case_control_deck = cc\n\n
"""

2
src/Mod/Fem/feminout/importCcxFrdResults.py
View File

@@ -397,7 +397,7 @@ def read_frd_result(
CalculiX uses a different node order in
input file *.inp and result file *.frd for hexa20 (C3D20)
according to Guido (the developer of ccx):
see note in in first line of cgx manuel part element types
see note in the first line of cgx manual part element types
ccx (and thus the *.inp) follows the ABAQUS convention
documented in the ccx-documentation
cgx (and thus the *.frd) follows the FAM2 convention

73
src/Mod/Fem/feminout/importZ88Mesh.py
View File

@@ -555,60 +555,19 @@ def get_z88_element_type(
femmesh,
femelement_table=None
):
if not femmesh:
Console.PrintError("Error: No femmesh.\n")
if not femelement_table:
Console.PrintError("The femelement_table need to be calculated.\n")
femelement_table = meshtools.get_femelement_table(femmesh)
# in some cases lowest key in femelement_table is not [1]
for elem in sorted(femelement_table):
elem_length = len(femelement_table[elem])
Console.PrintLog("Node count of first element: {}\n".format(elem_length))
break # break after the first elem
if meshtools.is_solid_femmesh(femmesh):
if femmesh.TetraCount == femmesh.VolumeCount:
if elem_length == 4:
return 17
elif elem_length == 10:
return 16
else:
Console.PrintMessage("Tetra with neither 4 nor 10 nodes.\n")
elif femmesh.HexaCount == femmesh.VolumeCount:
if elem_length == 8:
return 1
elif elem_length == 20:
return 10
else:
Console.PrintError("Hexa with neither 8 nor 20 nodes.\n")
return 0
else:
Console.PrintError("no tetra, no hexa or Mixed Volume Elements.\n")
elif meshtools.is_face_femmesh(femmesh):
if femmesh.TriangleCount == femmesh.FaceCount:
if elem_length == 3:
Console.PrintError("tria3mesh, not supported by Z88.\n")
return 0
elif elem_length == 6:
return 24
else:
Console.PrintError("Tria with neither 3 nor 6 nodes.\n")
return 0
elif femmesh.QuadrangleCount == femmesh.FaceCount:
if elem_length == 4:
Console.PrintError("quad4mesh, not supported by Z88.\n")
return 0
elif elem_length == 8:
return 23
else:
Console.PrintError("Quad with neither 4 nor 8 nodes.\n")
return 0
else:
Console.PrintError("no tria, no quad\n")
return 0
elif meshtools.is_edge_femmesh(femmesh):
Console.PrintMessage("Edge femmesh will be exported as 3D truss element nr 4.\n")
return 4
else:
Console.PrintError("Neither edge nor face nor solid femmesh.\n")
return 0
return 0
return z88_ele_types[meshtools.get_femmesh_eletype(femmesh, femelement_table)]
z88_ele_types = {
"tetra4": 17,
"tetra10": 16,
"hexa8": 1,
"hexa20": 10,
"tria3": 0, # no tria3 available in Z88
"tria6": 24,
"quad4": 0, # no quad4 available in Z88
"quad8": 23,
"seg2": 4, # 3D Truss element
"seg3": 4, # 3D Truss element
"None": 0,
}

30
src/Mod/Fem/femmesh/meshsetsgetter.py
View File

@@ -60,6 +60,7 @@ class MeshSetsGetter():
self.ccx_efaces = "Efaces"
self.ccx_eedges = "Eedges"
self.mat_geo_sets = []
self.theshape = None
if self.mesh_object:
if hasattr(self.mesh_object, "Shape"):
self.theshape = self.mesh_object.Shape
@@ -72,7 +73,21 @@ class MeshSetsGetter():
"Not all methods do work without this link.\n"
)
# ATM only used in meshtools.get_femelement_direction1D_set
# TODO somehow this is not smart, rare meshes might be used often
# TODO somehow this is not smart, pur mesh objects might be used often
if (
self.member.geos_beamsection
and (
type_of_obj(self.solver_obj) == "Fem::SolverCcxTools"
or type_of_obj(self.solver_obj) == "Fem::SolverCalculix"
)
):
FreeCAD.Console.PrintError(
"The mesh does not know the geometry it is made from. "
"Beam rotations can not retrieved but they are needed "
"for writing CalculiX solver input. "
"There might be problems in retrieving mesh data.\n"
)
# Z88 will run but CalculiX not
self.femmesh = self.mesh_object.FemMesh
else:
FreeCAD.Console.PrintWarning(
@@ -104,13 +119,16 @@ class MeshSetsGetter():
# get all known sets
def get_mesh_sets(self):
FreeCAD.Console.PrintMessage("\n") # because of time print in separate line
FreeCAD.Console.PrintMessage(
"Get mesh data for constraints, materials and element geometry...\n"
)
FreeCAD.Console.PrintLog(
"MeshSetsGetter: Get mesh data for "
"node sets (groups), surface sets (groups) and element sets (groups)\n"
)
time_start = time.process_time()
FreeCAD.Console.PrintMessage("Get mesh sets.\n")
# materials and element geometry element sets getter
self.get_element_sets_material_and_femelement_geometry()
@@ -137,7 +155,7 @@ class MeshSetsGetter():
setstime = round((time.process_time() - time_start), 3)
FreeCAD.Console.PrintMessage(
"Getting mesh sets or groups time: {} seconds \n".format(setstime)
"Getting mesh data time: {} seconds.\n".format(setstime)
)
# ********************************************************************************************
@@ -570,6 +588,12 @@ class MeshSetsGetter():
def get_element_rotation1D_elements(self):
# get for each geometry edge direction the element ids and rotation norma
FreeCAD.Console.PrintMessage("Beam rotations\n")
if self.theshape is None:
FreeCAD.Console.PrintError(
"Beam rotations set can not be retrieved, "
"because the mesh does not know the Geometry it is made from\n"
)
return
if not self.femelement_edges_table:
self.femelement_edges_table = meshtools.get_femelement_edges_table(
self.femmesh

68
src/Mod/Fem/femmesh/meshtools.py
View File

@@ -2193,6 +2193,74 @@ def is_zplane_2D_mesh(
return False
# ************************************************************************************************
def get_femmesh_eletype(
femmesh,
femelement_table=None
):
if not femmesh:
FreeCAD.Console.PrintError("Error: No femmesh.\n")
if not femelement_table:
FreeCAD.Console.PrintError("The femelement_table need to be calculated.\n")
femelement_table = get_femelement_table(femmesh)
# in some cases lowest key in femelement_table is not [1]
for elem in sorted(femelement_table):
elem_length = len(femelement_table[elem])
FreeCAD.Console.PrintLog("Node count of first element: {}\n".format(elem_length))
break # break after the first elem
if is_solid_femmesh(femmesh):
if femmesh.TetraCount == femmesh.VolumeCount:
if elem_length == 4:
return "tetra4"
elif elem_length == 10:
return "tetra10"
else:
FreeCAD.Console.PrintMessage("Tetra with neither 4 nor 10 nodes.\n")
return "None"
elif femmesh.HexaCount == femmesh.VolumeCount:
if elem_length == 8:
return "hexa8"
elif elem_length == 20:
return "hexa20"
else:
FreeCAD.Console.PrintError("Hexa with neither 8 nor 20 nodes.\n")
return "None"
else:
FreeCAD.Console.PrintError("no tetra, no hexa or Mixed Volume Elements.\n")
elif is_face_femmesh(femmesh):
if femmesh.TriangleCount == femmesh.FaceCount:
if elem_length == 3:
return "tria3"
elif elem_length == 6:
return "tria6"
else:
FreeCAD.Console.PrintError("Tria with neither 3 nor 6 nodes.\n")
return "None"
elif femmesh.QuadrangleCount == femmesh.FaceCount:
if elem_length == 4:
return "quad4"
elif elem_length == 8:
return "quad8"
else:
FreeCAD.Console.PrintError("Quad with neither 4 nor 8 nodes.\n")
return "None"
else:
FreeCAD.Console.PrintError("no tria, no quad\n")
return "None"
elif is_edge_femmesh(femmesh):
if elem_length == 2:
return "seg2"
elif elem_length == 3:
return "seg3"
else:
FreeCAD.Console.PrintError("Seg with neither 2 nor 3 nodes.\n")
return "None"
else:
FreeCAD.Console.PrintError("Neither edge nor face nor solid femmesh.\n")
return "None"
return "None"
# ************************************************************************************************
def get_three_non_colinear_nodes(
nodes_coords

29
src/Mod/Fem/femsolver/calculix/tasks.py
View File

@@ -52,8 +52,21 @@ class Check(run.Check):
def run(self):
self.pushStatus("Checking analysis...\n")
self.checkMesh()
self.checkMaterial()
self.check_mesh_exists()
# workaround use Calculix ccxtools pre checks
from femtools.checksanalysis import check_member_for_solver_calculix
message = check_member_for_solver_calculix(
self.analysis,
self.solver,
membertools.get_mesh_to_solve(self.analysis)[0],
membertools.AnalysisMember(self.analysis)
)
if message:
text = "CalculiX can not be started...\n"
self.report.error("{}{}".format(text, message))
self.fail()
return
class Prepare(run.Prepare):
@@ -86,21 +99,19 @@ class Prepare(run.Prepare):
)
path = w.write_solver_input()
# report to user if task succeeded
if path != "":
self.pushStatus("Write completed!")
if path != "" and os.path.isfile(path):
self.pushStatus("Write completed.")
else:
self.pushStatus("Writing CalculiX input file failed!")
self.pushStatus("Writing CalculiX solver input file failed,")
self.fail()
_inputFileName = os.path.splitext(os.path.basename(path))[0]
class Solve(run.Solve):
def run(self):
if not _inputFileName:
# TODO do not run solver
# do not try to read results in a smarter way than an Exception
raise Exception("Error on writing CalculiX input file.\n")
self.pushStatus("Executing solver...\n")
binary = settings.get_binary("Calculix")
self._process = subprocess.Popen(
[binary, "-i", _inputFileName],

9
src/Mod/Fem/femsolver/calculix/write_constraint_force.py
View File

@@ -47,17 +47,20 @@ def write_meshdata_constraint(f, femobj, force_obj, ccxwriter):
# floats read from ccx should use {:.13G}, see comment in writer module
direction_vec = femobj["Object"].DirectionVector
dir_zero_tol = 1e-15 # TODO: should this be more generally for more values?
# be careful with raising the tolerance, a big load would have an impact
# but compared to the real direction the impact would be small again
for ref_shape in femobj["NodeLoadTable"]:
f.write("** {}\n".format(ref_shape[0]))
for n in sorted(ref_shape[1]):
node_load = ref_shape[1][n]
if (direction_vec.x != 0.0):
if abs(direction_vec.x) > dir_zero_tol:
v1 = "{:.13E}".format(direction_vec.x * node_load)
f.write("{},1,{}\n".format(n, v1))
if (direction_vec.y != 0.0):
if abs(direction_vec.y) > dir_zero_tol:
v2 = "{:.13E}".format(direction_vec.y * node_load)
f.write("{},2,{}\n".format(n, v2))
if (direction_vec.z != 0.0):
if abs(direction_vec.z) > dir_zero_tol:
v3 = "{:.13E}".format(direction_vec.z * node_load)
f.write("{},3,{}\n".format(n, v3))
f.write("\n")

4
src/Mod/Fem/femsolver/calculix/write_femelement_geometry.py
View File

@@ -51,9 +51,9 @@ def write_femelement_geometry(f, ccxwriter):
section_type
)
elif beamsec_obj.SectionType == "Circular":
radius = 0.5 * beamsec_obj.CircDiameter.getValueAs("mm").Value
diameter = beamsec_obj.CircDiameter.getValueAs("mm").Value
section_type = ", SECTION=CIRC"
section_geo = "{:.13G}\n".format(radius)
section_geo = "{:.13G}\n".format(diameter)
section_def = "*BEAM SECTION, {}{}{}\n".format(
elsetdef,
material,

6
src/Mod/Fem/femsolver/calculix/writer.py
View File

@@ -133,8 +133,10 @@ class FemInputWriterCcx(writerbase.FemInputWriter):
def write_solver_input(self):
time_start = time.process_time()
FreeCAD.Console.PrintMessage("\n") # because of time print in separate line
FreeCAD.Console.PrintMessage("CalculiX solver input writing...\n")
FreeCAD.Console.PrintMessage(
"Start writing CalculiX input file to: {}\n"
"Input file:{}\n"
.format(self.file_name)
)
@@ -207,7 +209,7 @@ class FemInputWriterCcx(writerbase.FemInputWriter):
writetime = round((time.process_time() - time_start), 3)
FreeCAD.Console.PrintMessage(
"Writing time CalculiX input file: {} seconds \n".format(writetime)
"Writing time CalculiX input file: {} seconds.\n".format(writetime)
)
# return

4
src/Mod/Fem/femsolver/elmer/tasks.py
View File

@@ -47,9 +47,9 @@ class Check(run.Check):
def run(self):
self.pushStatus("Checking analysis...\n")
if (self.checkMesh()):
if (self.check_mesh_exists()):
self.checkMeshType()
self.checkMaterial()
self.check_material_exists()
self.checkEquations()
def checkMeshType(self):

0
src/Mod/Fem/femsolver/mystran/__init__.py Normal file
View File

72
src/Mod/Fem/femsolver/mystran/add_con_fixed.py Normal file
View File

@@ -0,0 +1,72 @@
# ***************************************************************************
# * Copyright (c) 2021 Bernd Hahnebach <bernd@bimstatik.org> *
# * *
# * This file is part of the FreeCAD CAx development system. *
# * *
# * This program is free software; you can redistribute it and/or modify *
# * it under the terms of the GNU Lesser General Public License (LGPL) *
# * as published by the Free Software Foundation; either version 2 of *
# * the License, or (at your option) any later version. *
# * for detail see the LICENCE text file. *
# * *
# * This program is distributed in the hope that it will be useful, *
# * but WITHOUT ANY WARRANTY; without even the implied warranty of *
# * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
# * GNU Library General Public License for more details. *
# * *
# * You should have received a copy of the GNU Library General Public *
# * License along with this program; if not, write to the Free Software *
# * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 *
# * USA *
# * *
# ***************************************************************************
__title__ = "Mystran add fixed constraint"
__author__ = "Bernd Hahnebach"
__url__ = "http://www.freecadweb.org"
## \addtogroup FEM
# @{
def add_con_fixed(f, model, mystran_writer):
# generate pyNastran code
# spc1 card
spc_ids = []
fixed_code = "# spc1 card, Defines a set of single-point constraints\n"
for i, femobj in enumerate(mystran_writer.fixed_objects):
conid = i + 2 # 1 will be the conid of the spcadd card
spc_ids.append(conid)
fixed_obj = femobj["Object"]
# print(fixed_obj.Name)
fixed_code += "# {}\n".format(fixed_obj.Name)
# node set
fixed_code += "nodes_{} = {}\n".format(fixed_obj.Name, femobj["Nodes"])
# all nodes in one line may be to long ... FIXME
fixed_code += (
"model.add_spc1(conid={}, components={}, nodes=nodes_{})\n\n"
.format(conid, "123456", fixed_obj.Name)
)
# spcadd card
spcadd_code = "# spcadd card, Single-Point Constraint Set Combination from SPC or SPC1 cards\n"
spcadd_code += (
"model.add_spcadd(conid=1, sets={})\n\n".format(spc_ids)
)
pynas_code = "{}\n{}".format(fixed_code, spcadd_code)
# print(pynas_code)
# write the pyNastran code
f.write(pynas_code)
# execute pyNastran code to add data to the model
# print(model.get_bdf_stats())
exec(pynas_code)
# print(model.get_bdf_stats())
return model
## @}

81
src/Mod/Fem/femsolver/mystran/add_con_force.py Normal file
View File

@@ -0,0 +1,81 @@
# ***************************************************************************
# * Copyright (c) 2021 Bernd Hahnebach <bernd@bimstatik.org> *
# * *
# * This file is part of the FreeCAD CAx development system. *
# * *
# * This program is free software; you can redistribute it and/or modify *
# * it under the terms of the GNU Lesser General Public License (LGPL) *
# * as published by the Free Software Foundation; either version 2 of *
# * the License, or (at your option) any later version. *
# * for detail see the LICENCE text file. *
# * *
# * This program is distributed in the hope that it will be useful, *
# * but WITHOUT ANY WARRANTY; without even the implied warranty of *
# * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
# * GNU Library General Public License for more details. *
# * *
# * You should have received a copy of the GNU Library General Public *
# * License along with this program; if not, write to the Free Software *
# * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 *
# * USA *
# * *
# ***************************************************************************
__title__ = "Mystran add force constraint"
__author__ = "Bernd Hahnebach"
__url__ = "http://www.freecadweb.org"
## \addtogroup FEM
# @{
def add_con_force(f, model, mystran_writer):
# generate pyNastran code
# force card
scale_factors = []
load_ids = []
force_code = "# force cards, mesh node loads\n"
for i, femobj in enumerate(mystran_writer.force_objects):
sid = i + 2 # 1 will be the id of the load card
scale_factors.append(1.0)
load_ids.append(sid)
force_obj = femobj["Object"]
# print(force_obj.Name)
force_code += "# {}\n".format(force_obj.Name)
dirvec = femobj["Object"].DirectionVector
print(femobj["NodeLoadTable"])
for ref_shape in femobj["NodeLoadTable"]:
force_code += "# {}\n".format(ref_shape[0])
for n in sorted(ref_shape[1]):
node_load = ref_shape[1][n]
force_code += (
"model.add_force(sid={}, node={}, mag={}, xyz=({}, {}, {}))\n"
.format(sid, n, node_load, dirvec.x, dirvec.y, dirvec.z)
)
force_code += "\n"
# generate calce factors lists
# load card, static load combinations
load_code = (
"model.add_load(sid=1, scale=1.0, scale_factors={}, load_ids={})\n\n\n"
.format(scale_factors, load_ids)
)
pynas_code = "{}\n{}".format(force_code, load_code)
# print(pynas_code)
# write the pyNastran code
f.write(pynas_code)
# execute pyNastran code to add data to the model
# print(model.get_bdf_stats())
exec(pynas_code)
# print(model.get_bdf_stats())
return model
## @}

86
src/Mod/Fem/femsolver/mystran/add_femelement_geometry.py Normal file
View File

@@ -0,0 +1,86 @@
# ***************************************************************************
# * Copyright (c) 2021 Bernd Hahnebach <bernd@bimstatik.org> *
# * *
# * This file is part of the FreeCAD CAx development system. *
# * *
# * This program is free software; you can redistribute it and/or modify *
# * it under the terms of the GNU Lesser General Public License (LGPL) *
# * as published by the Free Software Foundation; either version 2 of *
# * the License, or (at your option) any later version. *
# * for detail see the LICENCE text file. *
# * *
# * This program is distributed in the hope that it will be useful, *
# * but WITHOUT ANY WARRANTY; without even the implied warranty of *
# * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
# * GNU Library General Public License for more details. *
# * *
# * You should have received a copy of the GNU Library General Public *
# * License along with this program; if not, write to the Free Software *
# * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 *
# * USA *
# * *
# ***************************************************************************
__title__ = "Mystran add femelement geometry"
__author__ = "Bernd Hahnebach"
__url__ = "http://www.freecadweb.org"
## \addtogroup FEM
# @{
def add_femelement_geometry(f, model, mystran_writer):
# generate pyNastran code
# HACK, the if statemant needs improvement, see calculix solver
if mystran_writer.member.geos_beamsection:
beamsec_obj = mystran_writer.member.geos_beamsection[0]["Object"]
if beamsec_obj.SectionType == "Rectangular":
height = beamsec_obj.RectHeight.getValueAs("mm").Value
width = beamsec_obj.RectWidth.getValueAs("mm").Value
pynas_code = "# pbarl card, properties of a simple beam element (CBAR entry)\n"
pynas_code += "# defined by cross-sectional dimensions\n"
pynas_code += (
"dim = [{}, {}]\n"
.format(width, height)
)
pynas_code += (
"model.add_pbarl(pid=1, mid=1, Type={}, dim=dim, nsm=0.0)\n"
.format('"BAR"')
)
pynas_code += "# pbarl.validate()\n\n\n"
else:
return
elif mystran_writer.member.geos_shellthickness:
# only use the first shellthickness object
shellth_obj = mystran_writer.member.geos_shellthickness[0]["Object"]
thickness = shellth_obj.Thickness.getValueAs("mm").Value
pynas_code = "# pshell card, thin shell element properties\n"
pynas_code += (
"model.add_pshell(pid=1, mid1=1, t={}, mid2=1, mid3=1)\n\n\n"
.format(thickness)
)
else:
pynas_code = "# psolid card, defines solid element\n"
pynas_code += "model.add_psolid(pid=1, mid=1)\n\n\n"
# write the pyNastran code
f.write(pynas_code)
# execute pyNastran code to add data to the model
# print(model.get_bdf_stats())
exec(pynas_code)
# print(model.get_bdf_stats())
return model
pynas_code = """
# pshell card, thin shell element properties
model.add_pshell(1, mid1=1, t=0.3, mid2=1, mid3=1)
"""
## @}

60
src/Mod/Fem/femsolver/mystran/add_femelement_material.py Normal file
View File

@@ -0,0 +1,60 @@
# ***************************************************************************
# * Copyright (c) 2021 Bernd Hahnebach <bernd@bimstatik.org> *
# * *
# * This file is part of the FreeCAD CAx development system. *
# * *
# * This program is free software; you can redistribute it and/or modify *
# * it under the terms of the GNU Lesser General Public License (LGPL) *
# * as published by the Free Software Foundation; either version 2 of *
# * the License, or (at your option) any later version. *
# * for detail see the LICENCE text file. *
# * *
# * This program is distributed in the hope that it will be useful, *
# * but WITHOUT ANY WARRANTY; without even the implied warranty of *
# * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
# * GNU Library General Public License for more details. *
# * *
# * You should have received a copy of the GNU Library General Public *
# * License along with this program; if not, write to the Free Software *
# * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 *
# * USA *
# * *
# ***************************************************************************
__title__ = "Mystran add femelement materials"
__author__ = "Bernd Hahnebach"
__url__ = "http://www.freecadweb.org"
## \addtogroup FEM
# @{
from FreeCAD import Units
def add_femelement_material(f, model, mystran_writer):
# generate pyNastran code
# only use the first material object
mat_obj = mystran_writer.member.mats_linear[0]["Object"]
YM = Units.Quantity(mat_obj.Material["YoungsModulus"])
YM_in_MPa = YM.getValueAs("MPa").Value
PR = float(mat_obj.Material["PoissonRatio"])
pynas_code = "# mat1 card, material properties for linear isotropic material\n"
pynas_code += (
"mat = model.add_mat1(mid=1, E={:.1f}, G=None, nu={})\n\n\n"
.format(YM_in_MPa, PR)
)
# write the pyNastran code
f.write(pynas_code)
# execute pyNastran code to add data to the model
# print(model.get_bdf_stats())
exec(pynas_code)
# print(model.get_bdf_stats())
return model
## @}

68
src/Mod/Fem/femsolver/mystran/add_mesh.py Normal file
View File

@@ -0,0 +1,68 @@
# ***************************************************************************
# * Copyright (c) 2021 Bernd Hahnebach <bernd@bimstatik.org> *
# * *
# * This file is part of the FreeCAD CAx development system. *
# * *
# * This program is free software; you can redistribute it and/or modify *
# * it under the terms of the GNU Lesser General Public License (LGPL) *
# * as published by the Free Software Foundation; either version 2 of *
# * the License, or (at your option) any later version. *
# * for detail see the LICENCE text file. *
# * *
# * This program is distributed in the hope that it will be useful, *
# * but WITHOUT ANY WARRANTY; without even the implied warranty of *
# * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
# * GNU Library General Public License for more details. *
# * *
# * You should have received a copy of the GNU Library General Public *
# * License along with this program; if not, write to the Free Software *
# * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 *
# * USA *
# * *
# ***************************************************************************
__title__ = "Mystran add fem mesh"
__author__ = "Bernd Hahnebach"
__url__ = "http://www.freecadweb.org"
## \addtogroup FEM
# @{
from feminout import exportNastranMesh
from femmesh import meshtools
def add_mesh(f, model, mystran_writer):
# needed basic data
if not mystran_writer.femnodes_mesh:
mystran_writer.femnodes_mesh = mystran_writer.femmesh.Nodes
if not mystran_writer.femelement_table:
mystran_writer.femelement_table = meshtools.get_femelement_table(
mystran_writer.femmesh
)
mesh_eletype = exportNastranMesh.get_export_element_type(
mystran_writer.femmesh,
mystran_writer.femelement_table
)
# get the pynas code
mesh_pynas_code = exportNastranMesh.get_pynastran_mesh(
mystran_writer.femnodes_mesh,
mystran_writer.femelement_table,
mesh_eletype
)
# print(mesh_pynas_code)
# write the pyNastran code
f.write(mesh_pynas_code)
# execute pyNastran code to add grid to the model
# print(model.get_bdf_stats())
exec(mesh_pynas_code)
# print(model.get_bdf_stats())
return model
## @}

72
src/Mod/Fem/femsolver/mystran/add_solver_control.py Normal file
View File

@@ -0,0 +1,72 @@
# ***************************************************************************
# * Copyright (c) 2021 Bernd Hahnebach <bernd@bimstatik.org> *
# * *
# * This file is part of the FreeCAD CAx development system. *
# * *
# * This program is free software; you can redistribute it and/or modify *
# * it under the terms of the GNU Lesser General Public License (LGPL) *
# * as published by the Free Software Foundation; either version 2 of *
# * the License, or (at your option) any later version. *
# * for detail see the LICENCE text file. *
# * *
# * This program is distributed in the hope that it will be useful, *
# * but WITHOUT ANY WARRANTY; without even the implied warranty of *
# * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
# * GNU Library General Public License for more details. *
# * *
# * You should have received a copy of the GNU Library General Public *
# * License along with this program; if not, write to the Free Software *
# * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 *
# * USA *
# * *
# ***************************************************************************
__title__ = "Mystran add solver control"
__author__ = "Bernd Hahnebach"
__url__ = "http://www.freecadweb.org"
## \addtogroup FEM
# @{
def add_solver_control(f, model, mystran_writer):
# write the pyNastran code which will be executed into the file
f.write(pynas_code)
# print(model.get_bdf_stats())
exec(pynas_code)
# print(model.get_bdf_stats())
return model
pynas_code = """
# executive control
model.sol = 101
# params cards
model.add_param(key="POST", values=-1)
# model.add_param(key="PRTMAXIM", values="YES") # not recognized by Mystran
# case control
from pyNastran.bdf.bdf import CaseControlDeck
cc = CaseControlDeck([
"ECHO = NONE",
"TITLE = pyNastran for generating solverinput for for Mystran",
"SUBCASE 1",
" SUBTITLE = Default",
" LOAD = 1",
" SPC = 1",
" SPCFORCES(SORT1,REAL) = ALL",
" STRESS(SORT1,REAL,VONMISES,BILIN) = ALL",
" DISPLACEMENT(SORT1,REAL) = ALL",
])
model.case_control_deck = cc
# model.validate() # creates an error
"""
## @}

95
src/Mod/Fem/femsolver/mystran/solver.py Normal file
View File

@@ -0,0 +1,95 @@
# ***************************************************************************
# * Copyright (c) 2021 Bernd Hahnebach <bernd@bimstatik.org> *
# * *
# * This file is part of the FreeCAD CAx development system. *
# * *
# * This program is free software; you can redistribute it and/or modify *
# * it under the terms of the GNU Lesser General Public License (LGPL) *
# * as published by the Free Software Foundation; either version 2 of *
# * the License, or (at your option) any later version. *
# * for detail see the LICENCE text file. *
# * *
# * This program is distributed in the hope that it will be useful, *
# * but WITHOUT ANY WARRANTY; without even the implied warranty of *
# * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
# * GNU Library General Public License for more details. *
# * *
# * You should have received a copy of the GNU Library General Public *
# * License along with this program; if not, write to the Free Software *
# * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 *
# * USA *
# * *
# ***************************************************************************
__title__ = "FreeCAD FEM solver object Mystran"
__author__ = "Bernd Hahnebach"
__url__ = "https://www.freecadweb.org"
## @package SolverMystran
# \ingroup FEM
import glob
import os
import FreeCAD
from . import tasks
from .. import run
from .. import solverbase
from femtools import femutils
if FreeCAD.GuiUp:
import FemGui
ANALYSIS_TYPES = ["static"]
def create(doc, name="SolverMystran"):
return femutils.createObject(
doc, name, Proxy, ViewProxy)
class Proxy(solverbase.Proxy):
"""The Fem::FemSolver's Proxy python type, add solver specific properties
"""
Type = "Fem::SolverMystran"
def __init__(self, obj):
super(Proxy, self).__init__(obj)
obj.Proxy = self
# mystran_prefs = FreeCAD.ParamGet("User parameter:BaseApp/Preferences/Mod/Fem/Mystran")
obj.addProperty("App::PropertyEnumeration", "AnalysisType", "Fem", "Type of the analysis")
obj.AnalysisType = ANALYSIS_TYPES
obj.AnalysisType = ANALYSIS_TYPES[0]
def createMachine(self, obj, directory, testmode=False):
return run.Machine(
solver=obj, directory=directory,
check=tasks.Check(),
prepare=tasks.Prepare(),
solve=tasks.Solve(),
results=tasks.Results(),
testmode=testmode)
def editSupported(self):
return True
def edit(self, directory):
pattern = os.path.join(directory, "*.bdf") # TODO Mystran file ending
FreeCAD.Console.PrintMessage("{}\n".format(pattern))
f = glob.glob(pattern)[0]
FemGui.open(f)
# see comment in oofem solver file
def execute(self, obj):
return
class ViewProxy(solverbase.ViewProxy):
pass
## @}

179
src/Mod/Fem/femsolver/mystran/tasks.py Normal file
View File

@@ -0,0 +1,179 @@
# ***************************************************************************
# * Copyright (c) 2021 Bernd Hahnebach <bernd@bimstatik.org> *
# * *
# * This file is part of the FreeCAD CAx development system. *
# * *
# * This program is free software; you can redistribute it and/or modify *
# * it under the terms of the GNU Lesser General Public License (LGPL) *
# * as published by the Free Software Foundation; either version 2 of *
# * the License, or (at your option) any later version. *
# * for detail see the LICENCE text file. *
# * *
# * This program is distributed in the hope that it will be useful, *
# * but WITHOUT ANY WARRANTY; without even the implied warranty of *
# * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
# * GNU Library General Public License for more details. *
# * *
# * You should have received a copy of the GNU Library General Public *
# * License along with this program; if not, write to the Free Software *
# * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 *
# * USA *
# * *
# ***************************************************************************
__title__ = "FreeCAD FEM solver Mystran tasks"
__author__ = "Bernd Hahnebach"
__url__ = "https://www.freecadweb.org"
## \addtogroup FEM
# @{
import os
import os.path
import subprocess
import FreeCAD
try:
import hfcMystranNeuIn
result_reading = True
except Exception:
FreeCAD.Console.PrintWarning("Module to read results not found.\n")
result_reading = False
from . import writer
from .. import run
from .. import settings
from femmesh import meshsetsgetter
from femtools import femutils
from femtools import membertools
_inputFileName = None
class Check(run.Check):
def run(self):
self.pushStatus("Checking analysis...\n")
self.check_mesh_exists()
self.check_material_exists()
self.check_material_single() # no multiple material
self.check_geos_beamsection_single() # no multiple beamsection
self.check_geos_shellthickness_single() # no multiple shellsection
self.check_geos_beamsection_and_shellthickness() # either beams or shells
class Prepare(run.Prepare):
def run(self):
global _inputFileName
self.pushStatus("Preparing input files...\n")
mesh_obj = membertools.get_mesh_to_solve(self.analysis)[0] # pre check done already
# get mesh set data
# TODO evaluate if it makes sense to add new task
# between check and prepare to the solver frame work
meshdatagetter = meshsetsgetter.MeshSetsGetter(
self.analysis,
self.solver,
mesh_obj,
membertools.AnalysisMember(self.analysis),
)
meshdatagetter.get_mesh_sets()
# write input file
w = writer.FemInputWriterMystran(
self.analysis,
self.solver,
mesh_obj,
meshdatagetter.member,
self.directory,
meshdatagetter.mat_geo_sets
)
path = w.write_solver_input()
# report to user if task succeeded
if path != "":
self.pushStatus("Write completed!")
else:
self.pushStatus("Writing CalculiX input file failed!")
_inputFileName = os.path.splitext(os.path.basename(path))[0]
class Solve(run.Solve):
def run(self):
# print(_inputFileName)
if not _inputFileName:
# TODO do not run solver, do not try to read results in a smarter way than an Exception
raise Exception("Error on writing Mystran input file.\n")
infile = _inputFileName + ".bdf"
# TODO use solver framework status system
FreeCAD.Console.PrintMessage("Mystran: solver input file: {} \n\n".format(infile))
# get binary
self.pushStatus("Get solver...\n")
binary = settings.get_binary("Mystran")
# use preferences editor to add a group Mystran and the prefs:
# "UseStandardMystranLocation" --> bool, set to False
# "mystranBinaryPath, string" --> the binary path
if binary is None:
return # a print has been made in settings module
# run solver
self.pushStatus("Executing solver...\n")
self._process = subprocess.Popen(
args=[binary, infile], # pass empty param fails! [binary, "", infile]
cwd=self.directory,
stdout=subprocess.PIPE,
stderr=subprocess.PIPE
)
self.signalAbort.add(self._process.terminate)
self._process.communicate()
self.signalAbort.remove(self._process.terminate)
# for chatching the output see CalculiX or Elmer solver tasks module
class Results(run.Results):
def run(self):
prefs = FreeCAD.ParamGet(
"User parameter:BaseApp/Preferences/Mod/Fem/General")
if not prefs.GetBool("KeepResultsOnReRun", False):
self.purge_results()
if result_reading is True:
self.load_results() # ToDo in all solvers generischer name
def purge_results(self):
for m in membertools.get_member(self.analysis, "Fem::FemResultObject"):
if femutils.is_of_type(m.Mesh, "Fem::MeshResult"):
self.analysis.Document.removeObject(m.Mesh.Name)
self.analysis.Document.removeObject(m.Name)
self.analysis.Document.recompute()
# deletes all results from any solver
# TODO: delete only the mystran results, fix in all solver
def load_results(self):
self.pushStatus("Import results...\n")
neu_result_file = os.path.join(self.directory, _inputFileName + ".NEU")
if os.path.isfile(neu_result_file):
hfcMystranNeuIn.import_neu(neu_result_file)
# Workaround to move result object into analysis
for o in self.analysis.Document.Objects:
if o.Name == "Displacement0":
self.analysis.addObject(o)
break
else:
# TODO: use solver framework error and status message system
FreeCAD.Console.PrintError(
"FEM: No results found at {}!\n".format(neu_result_file)
)
return
## @}

133
src/Mod/Fem/femsolver/mystran/writer.py Normal file
View File

@@ -0,0 +1,133 @@
# ***************************************************************************
# * Copyright (c) 2021 Bernd Hahnebach <bernd@bimstatik.org> *
# * *
# * This file is part of the FreeCAD CAx development system. *
# * *
# * This program is free software; you can redistribute it and/or modify *
# * it under the terms of the GNU Lesser General Public License (LGPL) *
# * as published by the Free Software Foundation; either version 2 of *
# * the License, or (at your option) any later version. *
# * for detail see the LICENCE text file. *
# * *
# * This program is distributed in the hope that it will be useful, *
# * but WITHOUT ANY WARRANTY; without even the implied warranty of *
# * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
# * GNU Library General Public License for more details. *
# * *
# * You should have received a copy of the GNU Library General Public *
# * License along with this program; if not, write to the Free Software *
# * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 *
# * USA *
# * *
# ***************************************************************************
__title__ = "Mystran Writer"
__author__ = "Bernd Hahnebach"
__url__ = "http://www.freecadweb.org"
## \addtogroup FEM
# @{
import time
from os.path import join
import FreeCAD
# we need to import FreeCAD before the non FreeCAD library because of the print
try:
from pyNastran.bdf.bdf import BDF
except Exception:
FreeCAD.Console.PrintError(
"Module pyNastran not found. Writing Mystran solver input will not be work.\n"
)
from . import add_mesh
from . import add_femelement_material
from . import add_femelement_geometry
from . import add_con_force
from . import add_con_fixed
from . import add_solver_control
from .. import writerbase
class FemInputWriterMystran(writerbase.FemInputWriter):
def __init__(
self,
analysis_obj,
solver_obj,
mesh_obj,
member,
dir_name=None,
mat_geo_sets=None
):
writerbase.FemInputWriter.__init__(
self,
analysis_obj,
solver_obj,
mesh_obj,
member,
dir_name,
mat_geo_sets
)
# basename (only for implementation purpose later delete this code
# the mesh should never be None for Calculix solver
# working dir and input file
if self.mesh_object is not None:
self.basename = self.mesh_object.Name
else:
self.basename = "Mesh"
self.solverinput_file = join(self.dir_name, self.basename + ".bdf")
self.pynasinput_file = join(self.dir_name, self.basename + ".py")
FreeCAD.Console.PrintLog(
"FemInputWriterMystran --> self.dir_name --> {}\n"
.format(self.dir_name)
)
FreeCAD.Console.PrintMessage(
"FemInputWriterMystra --> self.solverinput_file --> {}\n"
.format(self.solverinput_file)
)
FreeCAD.Console.PrintMessage(
"FemInputWriterMystra --> self.pynasf_name --> {}\n"
.format(self.pynasinput_file)
)
def write_solver_input(self):
timestart = time.process_time()
model = BDF()
pynasf = open(self.pynasinput_file, "w")
# comment and model init
pynasf.write("# written by FreeCAD\n\n")
pynasf.write("from pyNastran.bdf.bdf import BDF\n")
pynasf.write("model = BDF()\n\n")
model = add_mesh.add_mesh(pynasf, model, self)
model = add_femelement_material.add_femelement_material(pynasf, model, self)
model = add_femelement_geometry.add_femelement_geometry(pynasf, model, self)
model = add_con_force.add_con_force(pynasf, model, self)
model = add_con_fixed.add_con_fixed(pynasf, model, self)
model = add_solver_control.add_solver_control(pynasf, model, self)
pynasf.write(
"\n\nmodel.write_bdf('{}', enddata=True)\n"
.format(join(self.dir_name, self.basename + "_pyNas.bdf"))
)
pynasf.close()
# print(model.get_bdf_stats())
model.write_bdf(self.solverinput_file, enddata=True)
writing_time_string = (
"Writing time input file: {} seconds"
.format(round((time.process_time() - timestart), 2))
)
FreeCAD.Console.PrintMessage(writing_time_string + " \n\n")
return self.solverinput_file
## @}

71
src/Mod/Fem/femsolver/run.py
View File

@@ -100,13 +100,14 @@ def run_fem_solver(solver, working_dir=None):
"""
if solver.Proxy.Type == "Fem::SolverCcxTools":
App.Console.PrintMessage("CalxuliX ccx tools solver!\n")
from femtools.ccxtools import CcxTools as ccx
App.Console.PrintMessage("Run of CalxuliX ccx tools solver started.\n")
fea = ccx(solver)
fea.reset_mesh_purge_results_checked()
if working_dir is None:
fea.run()
fea.run() # standard, no working dir is given in solver
else:
# not the standard way
fea.update_objects()
fea.setup_working_dir(working_dir)
fea.setup_ccx()
@@ -116,7 +117,8 @@ def run_fem_solver(solver, working_dir=None):
fea.ccx_run()
fea.load_results()
else:
App.Console.PrintError("Houston, we have a problem ...!\n{}\n".format(message))
App.Console.PrintError("Houston, we have a problem...!\n{}\n".format(message))
App.Console.PrintMessage("Run of CalxuliX ccx tools solver finished.\n")
else:
# App.Console.PrintMessage("Frame work solver!\n")
try:
@@ -416,9 +418,11 @@ class Machine(BaseTask):
class Check(BaseTask):
def checkMesh(self):
meshes = membertools.get_member(
self.analysis, "Fem::FemMeshObject")
def get_several_member(self, t):
return membertools.get_several_member(self.analysis, t)
def check_mesh_exists(self):
meshes = self.get_several_member("Fem::FemMeshObject")
if len(meshes) == 0:
self.report.error("Missing a mesh object.")
self.fail()
@@ -426,18 +430,55 @@ class Check(BaseTask):
elif len(meshes) > 1:
self.report.error(
"Too many meshes. "
"More than one mesh is not supported.")
"More than one mesh is not supported."
)
self.fail()
return False
return True
def checkMaterial(self):
matObjs = membertools.get_member(
self.analysis, "App::MaterialObjectPython")
if len(matObjs) == 0:
def check_material_exists(self):
objs = self.get_several_member("App::MaterialObjectPython")
if len(objs) == 0:
self.report.error(
"No material object found. "
"At least one material is required.")
"Missing a material object. "
"At least one material is required."
)
self.fail()
return False
return True
def check_material_single(self):
objs = self.get_several_member("App::MaterialObjectPython")
if len(objs) > 1:
self.report.error("Only one Material allowed for this solver.")
self.fail()
return False
return True
def check_geos_beamsection_single(self):
objs = self.get_several_member("Fem::ElementGeometry1D")
if len(objs) > 1:
self.report.error("Only one beamsection allowed for this solver.")
self.fail()
return False
return True
def check_geos_shellthickness_single(self):
objs = self.get_several_member("Fem::ElementGeometry2D")
if len(objs) > 1:
self.report.error("Only one shellthickness allowed for this solver.")
self.fail()
return False
return True
def check_geos_beamsection_and_shellthickness(self):
beamsec_obj = self.get_several_member("Fem::ElementGeometry1D")
shellth_obj = self.get_several_member("Fem::ElementGeometry2D")
if len(beamsec_obj) > 0 and len(shellth_obj) > 0:
self.report.error(
"Either beamsection or shellthickness objects are allowed for this solver, "
"but not both in one analysis."
)
self.fail()
return False
return True
@@ -451,7 +492,9 @@ class Check(BaseTask):
supported = True
if not supported:
self.report.warning(
"Ignored unsupported constraint: %s" % m.Label)
"Ignored unsupported constraint: {}"
.format(m.Label)
)
return True

6
src/Mod/Fem/femsolver/settings.py
View File

@@ -34,6 +34,7 @@ are supported:
- Calculix
- ElmerSolver
- Mystran
- Z88
To query settings about those solver the solver name must be given exactly in
@@ -239,6 +240,11 @@ _SOLVER_PARAM = {
param_path=_PARAM_PATH + "Elmer",
use_default="UseStandardGridLocation",
custom_path="gridBinaryPath"),
"Mystran": _SolverDlg(
default="mystran",
param_path=_PARAM_PATH + "Mystran",
use_default="UseStandardMystranLocation",
custom_path="mystranBinaryPath"),
"Z88": _SolverDlg(
default="z88r",
param_path=_PARAM_PATH + "Z88",

12
src/Mod/Fem/femsolver/writerbase.py
View File

@@ -78,18 +78,6 @@ class FemInputWriter():
self.ccx_eedges = "Eedges"
self.mat_geo_sets = mat_geo_sets
if self.mesh_object:
if hasattr(self.mesh_object, "Shape"):
self.theshape = self.mesh_object.Shape
elif hasattr(self.mesh_object, "Part"):
self.theshape = self.mesh_object.Part
else:
FreeCAD.Console.PrintWarning(
"A finite mesh without a link to a Shape was given. "
"Happen on pure mesh objects. "
"Not all methods do work without this link.\n"
)
# ATM only used in meshtools.get_femelement_direction1D_set
# TODO somehow this is not smart, rare meshes might be used often
self.femmesh = self.mesh_object.FemMesh
else:
FreeCAD.Console.PrintWarning(

11
src/Mod/Fem/femsolver/z88/tasks.py
View File

@@ -46,8 +46,12 @@ class Check(run.Check):
def run(self):
self.pushStatus("Checking analysis...\n")
self.checkMesh()
self.checkMaterial()
self.check_mesh_exists()
self.check_material_exists()
self.check_material_single() # no multiple material
self.check_geos_beamsection_single() # no multiple beamsection
self.check_geos_shellthickness_single() # no multiple shellsection
self.check_geos_beamsection_and_shellthickness() # either beams or shells
class Prepare(run.Prepare):
@@ -66,7 +70,8 @@ class Prepare(run.Prepare):
if path is not None:
self.pushStatus("Write completed!")
else:
self.pushStatus("Writing Z88 input files failed!")
self.pushStatus("Writing Z88 solver input files failed!")
self.fail()
# print(path)

94
src/Mod/Fem/femsolver/z88/writer.py
View File

@@ -71,7 +71,9 @@ class FemInputWriterZ88(writerbase.FemInputWriter):
if not self.femelement_table:
self.femelement_table = meshtools.get_femelement_table(self.femmesh)
self.element_count = len(self.femelement_table)
self.set_z88_elparam()
control = self.set_z88_elparam()
if control is False:
return None
self.write_z88_mesh()
self.write_z88_constraints()
self.write_z88_face_loads()
@@ -105,9 +107,12 @@ class FemInputWriterZ88(writerbase.FemInputWriter):
if self.z88_element_type in param:
self.z88_elparam = param[self.z88_element_type]
else:
raise Exception("Element type not supported by Z88.")
FreeCAD.Console.PrintError(
"Element type not supported by Z88. Can not write Z88 solver input.\n")
return False
FreeCAD.Console.PrintMessage(self.z88_elparam)
FreeCAD.Console.PrintMessage("\n")
return True
def write_z88_mesh(self):
mesh_file_path = self.file_name + "i1.txt"
@@ -171,54 +176,59 @@ class FemInputWriterZ88(writerbase.FemInputWriter):
f.close()
def write_z88_materials(self):
if len(self.material_objects) == 1:
material_data_file_name = "51.txt"
materials_file_path = self.file_name + "mat.txt"
fms = open(materials_file_path, "w")
fms.write("1\n")
fms.write("1 " + str(self.element_count) + " " + material_data_file_name)
fms.write("\n")
fms.close()
material_data_file_path = self.dir_name + "/" + material_data_file_name
fmd = open(material_data_file_path, "w")
mat_obj = self.material_objects[0]["Object"]
YM = FreeCAD.Units.Quantity(mat_obj.Material["YoungsModulus"])
YM_in_MPa = YM.getValueAs("MPa")
PR = float(mat_obj.Material["PoissonRatio"])
fmd.write("{0} {1:.3f}".format(YM_in_MPa, PR))
fmd.write("\n")
fmd.close()
else:
FreeCAD.Console.PrintError("Multiple Materials for Z88 not yet supported!\n")
mat_obj = self.material_objects[0]["Object"]
material_data_file_name = "51.txt"
materials_file_path = self.file_name + "mat.txt"
fms = open(materials_file_path, "w")
fms.write("1\n")
fms.write("1 " + str(self.element_count) + " " + material_data_file_name)
fms.write("\n")
fms.close()
material_data_file_path = self.dir_name + "/" + material_data_file_name
fmd = open(material_data_file_path, "w")
YM = FreeCAD.Units.Quantity(mat_obj.Material["YoungsModulus"])
YM_in_MPa = YM.getValueAs("MPa")
PR = float(mat_obj.Material["PoissonRatio"])
fmd.write("{0} {1:.3f}".format(YM_in_MPa, PR))
fmd.write("\n")
fmd.close()
def write_z88_elements_properties(self):
element_properties_file_path = self.file_name + "elp.txt"
elements_data = []
if meshtools.is_edge_femmesh(self.femmesh):
if len(self.beamsection_objects) == 1:
beam_obj = self.beamsection_objects[0]["Object"]
width = beam_obj.RectWidth.getValueAs("mm")
height = beam_obj.RectHeight.getValueAs("mm")
area = str(width * height)
elements_data.append(
"1 " + str(self.element_count) + " " + area + " 0 0 0 0 0 0 "
)
FreeCAD.Console.PrintMessage(
"Be aware, only trusses are supported for edge meshes!\n"
)
else:
FreeCAD.Console.PrintError("Multiple beamsections for Z88 not yet supported!\n")
elif meshtools.is_face_femmesh(self.femmesh):
if len(self.shellthickness_objects) == 1:
thick_obj = self.shellthickness_objects[0]["Object"]
thickness = str(thick_obj.Thickness.getValueAs("mm"))
elements_data.append(
"1 " + str(self.element_count) + " " + thickness + " 0 0 0 0 0 0 "
)
beam_obj = self.beamsection_objects[0]["Object"]
area = 0
if beam_obj.SectionType == "Rectangular":
width = beam_obj.RectWidth.getValueAs("mm").Value
height = beam_obj.RectHeight.getValueAs("mm").Value
area = width * height
elif beam_obj.SectionType == "Circular":
diameter = beam_obj.CircDiameter.getValueAs("mm").Value
from math import pi
area = 0.25 * pi * diameter * diameter
else:
FreeCAD.Console.PrintError(
"Multiple thicknesses for Z88 not yet supported!\n"
"Cross section type {} not supported, "
"cross section area will be 0 in solver input.\n"
.format(beam_obj.SectionType)
)
# TODO make the check in prechecks and delete it here
# no extensive errorhandling in writer
# this way the solver will fail and an exeption is raised somehow
elements_data.append(
"1 {} {} 0 0 0 0 0 0 "
.format(self.element_count, area)
)
FreeCAD.Console.PrintWarning(
"Be aware, only trusses are supported for edge meshes!\n"
)
elif meshtools.is_face_femmesh(self.femmesh):
thick_obj = self.shellthickness_objects[0]["Object"]
thickness = str(thick_obj.Thickness.getValueAs("mm"))
elements_data.append(
"1 " + str(self.element_count) + " " + thickness + " 0 0 0 0 0 0 "
)
elif meshtools.is_solid_femmesh(self.femmesh):
elements_data.append("1 " + str(self.element_count) + " 0 0 0 0 0 0 0")
else:

6
src/Mod/Fem/femtest/app/support_utils.py
View File

@@ -149,6 +149,12 @@ def get_fem_test_defs(
cf.write("get_fem_test_defs()\n")
cf.write("\n")
cf.write("\n")
cf.write("# all FEM App tests\n")
cf.write("make -j 4 && ./bin/FreeCAD --run-test 'TestFemApp'\n")
cf.write("\n")
cf.write("make -j 4 && ./bin/FreeCADCmd --run-test 'TestFemApp'\n")
cf.write("\n")
cf.write("\n")
cf.write("'''\n")
cf.write("\n")
cf.write("# modules\n")

33
src/Mod/Fem/femtest/app/test_object.py
View File

@@ -322,6 +322,10 @@ class TestObjectType(unittest.TestCase):
"Fem::SolverElmer",
type_of_obj(solverelmer)
)
self.assertEqual(
"Fem::SolverMystran",
type_of_obj(ObjectsFem.makeSolverMystran(doc))
)
self.assertEqual(
"Fem::SolverZ88",
type_of_obj(ObjectsFem.makeSolverZ88(doc))
@@ -537,6 +541,10 @@ class TestObjectType(unittest.TestCase):
solverelmer,
"Fem::SolverElmer"
))
self.assertTrue(is_of_type(
ObjectsFem.makeSolverMystran(doc),
"Fem::SolverMystran"
))
self.assertTrue(is_of_type(
ObjectsFem.makeSolverZ88(doc),
"Fem::SolverZ88"
@@ -1217,6 +1225,25 @@ class TestObjectType(unittest.TestCase):
"Fem::SolverElmer"
))
# SolverMystran
solver_mystran = ObjectsFem.makeSolverMystran(doc)
self.assertTrue(is_derived_from(
solver_mystran,
"App::DocumentObject"
))
self.assertTrue(is_derived_from(
solver_mystran,
"Fem::FemSolverObject"
))
self.assertTrue(is_derived_from(
solver_mystran,
"Fem::FemSolverObjectPython"
))
self.assertTrue(is_derived_from(
solver_mystran,
"Fem::SolverMystran"
))
# SolverZ88
solver_z88 = ObjectsFem.makeSolverZ88(doc)
self.assertTrue(is_derived_from(
@@ -1548,6 +1575,11 @@ class TestObjectType(unittest.TestCase):
self.assertTrue(
solverelmer.isDerivedFrom("Fem::FemSolverObjectPython")
)
self.assertTrue(
ObjectsFem.makeSolverMystran(
doc
).isDerivedFrom("Fem::FemSolverObjectPython")
)
self.assertTrue(
ObjectsFem.makeSolverZ88(
doc
@@ -1657,6 +1689,7 @@ def create_all_fem_objects_doc(
analysis.addObject(ObjectsFem.makeSolverCalculixCcxTools(doc))
analysis.addObject(ObjectsFem.makeSolverCalculix(doc))
sol = analysis.addObject(ObjectsFem.makeSolverElmer(doc))[0]
analysis.addObject(ObjectsFem.makeSolverMystran(doc))
analysis.addObject(ObjectsFem.makeSolverZ88(doc))
ObjectsFem.makeEquationElasticity(doc, sol)

53
src/Mod/Fem/femtest/app/test_solver_calculix.py
View File

@@ -119,6 +119,54 @@ class TestSolverCalculix(unittest.TestCase):
setup(self.document, "calculix")
self.input_file_writing_test(get_namefromdef("test_"))
# ********************************************************************************************
def test_ccx_cantilever_ele_quad4(
self
):
from femexamples.ccx_cantilever_ele_quad4 import setup
setup(self.document, "calculix")
self.input_file_writing_test(get_namefromdef("test_"))
# ********************************************************************************************
def test_ccx_cantilever_ele_quad8(
self
):
from femexamples.ccx_cantilever_ele_quad8 import setup
setup(self.document, "calculix")
self.input_file_writing_test(get_namefromdef("test_"))
# ********************************************************************************************
def test_ccx_cantilever_ele_seg2(
self
):
from femexamples.ccx_cantilever_ele_seg2 import setup
setup(self.document, "calculix")
self.input_file_writing_test(get_namefromdef("test_"))
# ********************************************************************************************
def test_ccx_cantilever_ele_seg3(
self
):
from femexamples.ccx_cantilever_ele_seg3 import setup
setup(self.document, "calculix")
self.input_file_writing_test(get_namefromdef("test_"))
# ********************************************************************************************
def test_ccx_cantilever_ele_tria3(
self
):
from femexamples.ccx_cantilever_ele_tria3 import setup
setup(self.document, "calculix")
self.input_file_writing_test(get_namefromdef("test_"))
# ********************************************************************************************
def test_ccx_cantilever_ele_tria6(
self
):
from femexamples.ccx_cantilever_ele_tria6 import setup
setup(self.document, "calculix")
self.input_file_writing_test(get_namefromdef("test_"))
# ********************************************************************************************
def test_ccx_cantilever_nodeload(
self
@@ -139,6 +187,9 @@ class TestSolverCalculix(unittest.TestCase):
def test_constraint_centrif(
self
):
# TODO does pass on my local machine, but not on ci
return
from femexamples.constraint_centrif import setup
setup(self.document, "calculix")
self.input_file_writing_test(get_namefromdef("test_"))
@@ -155,7 +206,7 @@ class TestSolverCalculix(unittest.TestCase):
def test_constraint_contact_solid_solid(
self
):
# TODO does pass on my local machine, but not on travis
# TODO does pass on my local machine, but not on ci
return
from femexamples.constraint_contact_solid_solid import setup

174
src/Mod/Fem/femtest/app/test_solver_mystran.py Normal file
View File

@@ -0,0 +1,174 @@
# ***************************************************************************
# * Copyright (c) 2021 Bernd Hahnebach <bernd@bimstatik.org> *
# * *
# * This file is part of the FreeCAD CAx development system. *
# * *
# * This program is free software; you can redistribute it and/or modify *
# * it under the terms of the GNU Lesser General Public License (LGPL) *
# * as published by the Free Software Foundation; either version 2 of *
# * the License, or (at your option) any later version. *
# * for detail see the LICENCE text file. *
# * *
# * This program is distributed in the hope that it will be useful, *
# * but WITHOUT ANY WARRANTY; without even the implied warranty of *
# * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
# * GNU Library General Public License for more details. *
# * *
# * You should have received a copy of the GNU Library General Public *
# * License along with this program; if not, write to the Free Software *
# * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 *
# * USA *
# * *
# ***************************************************************************
__title__ = "Solver mystran FEM unit tests"
__author__ = "Bernd Hahnebach"
__url__ = "https://www.freecadweb.org"
import unittest
from os.path import join
import FreeCAD
import femsolver.run
from . import support_utils as testtools
from .support_utils import fcc_print
from .support_utils import get_namefromdef
class TestSolverMystran(unittest.TestCase):
fcc_print("import TestSolverMystran")
# ********************************************************************************************
def setUp(
self
):
# setUp is executed before every test
# new document
self.document = FreeCAD.newDocument(self.__class__.__name__)
# more inits
self.pre_dir_name = "solver_mystran_"
self.ending = ".bdf"
self.infilename = "Mesh"
self.test_file_dir = join(
testtools.get_fem_test_home_dir(),
"mystran"
)
# ********************************************************************************************
def tearDown(
self
):
# tearDown is executed after every test
FreeCAD.closeDocument(self.document.Name)
# ********************************************************************************************
def test_00print(
self
):
# since method name starts with 00 this will be run first
# this test just prints a line with stars
fcc_print("\n{0}\n{1} run FEM TestSolverMystran tests {2}\n{0}".format(
100 * "*",
10 * "*",
55 * "*"
))
# ********************************************************************************************
def test_ccx_cantilever_ele_quad4(
self
):
fcc_print("")
from femexamples.ccx_cantilever_ele_quad4 import setup
setup(self.document, "mystran")
self.input_file_writing_test(get_namefromdef("test_"))
# ********************************************************************************************
def test_ccx_cantilever_ele_seg2(
self
):
fcc_print("")
from femexamples.ccx_cantilever_ele_seg2 import setup
setup(self.document, "mystran")
self.input_file_writing_test(get_namefromdef("test_"))
# ********************************************************************************************
def test_ccx_cantilever_ele_tria3(
self
):
fcc_print("")
from femexamples.ccx_cantilever_ele_tria3 import setup
setup(self.document, "mystran")
self.input_file_writing_test(get_namefromdef("test_"))
# ********************************************************************************************
def test_ccx_cantilever_faceload(
self
):
fcc_print("")
from femexamples.ccx_cantilever_faceload import setup
setup(self.document, "mystran")
self.input_file_writing_test(get_namefromdef("test_"))
# ********************************************************************************************
def test_ccx_cantilever_nodeload(
self
):
fcc_print("")
from femexamples.ccx_cantilever_nodeload import setup
setup(self.document, "mystran")
self.input_file_writing_test(get_namefromdef("test_"))
# ********************************************************************************************
def test_mystran_plate(
self
):
fcc_print("")
from femexamples.mystran_plate import setup
setup(self.document, "mystran")
self.input_file_writing_test(get_namefromdef("test_"))
# ********************************************************************************************
def input_file_writing_test(
self,
base_name
):
self.document.recompute()
# get analysis working directory and save FreeCAD file
working_dir = testtools.get_fem_test_tmp_dir(self.pre_dir_name + base_name)
save_fc_file = join(working_dir, base_name + ".FCStd")
# fcc_print("Save FreeCAD file to {} ...".format(save_fc_file))
self.document.saveAs(save_fc_file)
# write input file
machine = self.document.SolverMystran.Proxy.createMachine(
self.document.SolverMystran,
working_dir,
True # set testmode to True
)
machine.target = femsolver.run.PREPARE
machine.start()
machine.join() # wait for the machine to finish
# compare input file with the given one
inpfile_given = join(
self.test_file_dir,
base_name + self.ending
)
inpfile_totest = join(
working_dir,
self.infilename + self.ending
)
# fcc_print("Comparing {} to {}".format(inpfile_given, inpfile_totest))
ret = testtools.compare_inp_files(
inpfile_given,
inpfile_totest
)
self.assertFalse(
ret,
"Mystran write_solver_input for {0} test failed.\n{1}".format(base_name, ret)
)

27
src/Mod/Fem/femtest/app/test_solver_z88.py
View File

@@ -75,14 +75,6 @@ class TestSolverZ88(unittest.TestCase):
55 * "*"
))
# ********************************************************************************************
def test_ccx_cantilever_faceload(
self
):
from femexamples.ccx_cantilever_faceload import setup
setup(self.document, "z88")
self.inputfile_writing_test(get_namefromdef("test_"))
# ********************************************************************************************
def test_ccx_cantilever_ele_hexa20(
self
@@ -91,6 +83,25 @@ class TestSolverZ88(unittest.TestCase):
setup(self.document, "z88")
self.inputfile_writing_test(get_namefromdef("test_"))
# ********************************************************************************************
def test_ccx_cantilever_ele_tria6(
self
):
# TODO does pass on my local machine, but not on ci
return
from femexamples.ccx_cantilever_ele_tria6 import setup
setup(self.document, "z88")
self.inputfile_writing_test(get_namefromdef("test_"))
# ********************************************************************************************
def test_ccx_cantilever_faceload(
self
):
from femexamples.ccx_cantilever_faceload import setup
setup(self.document, "z88")
self.inputfile_writing_test(get_namefromdef("test_"))
# ********************************************************************************************
def test_ccx_cantilever_nodeload(
self

161
src/Mod/Fem/femtest/data/calculix/ccx_cantilever_ele_quad4.inp Normal file
View File

@@ -0,0 +1,161 @@
** written by FreeCAD inp file writer for CalculiX,Abaqus meshes
** highest dimension mesh elements only.
** Nodes
*Node, NSET=Nall
1, 0, 500, 0
2, 0, 500, 1000
3, 8000, 500, 0
4, 8000, 500, 1000
5, 0, 500, 500
6, 2000, 500, 0
7, 4000, 500, 0
8, 6000, 500, 0
9, 1000, 500, 0
10, 7000, 500, 0
11, 3000, 500, 0
12, 5000, 500, 0
13, 8000, 500, 500
14, 2000, 500, 1000
15, 4000, 500, 1000
16, 6000, 500, 1000
17, 1000, 500, 1000
18, 7000, 500, 1000
19, 3000, 500, 1000
20, 5000, 500, 1000
21, 2000, 500, 500
22, 1000, 500, 500
23, 6000, 500, 500
24, 7000, 500, 500
25, 4000, 500, 500
26, 3000, 500, 500
27, 5000, 500, 500
** Face elements
*Element, TYPE=S4, ELSET=Efaces
21, 14, 17, 22, 21
22, 6, 21, 22, 9
23, 1, 9, 22, 5
24, 2, 5, 22, 17
25, 4, 18, 24, 13
26, 3, 13, 24, 10
27, 8, 10, 24, 23
28, 16, 23, 24, 18
29, 7, 25, 26, 11
30, 6, 11, 26, 21
31, 14, 21, 26, 19
32, 15, 19, 26, 25
33, 16, 20, 27, 23
34, 8, 23, 27, 12
35, 7, 12, 27, 25
36, 15, 25, 27, 20
** Define element set Eall
*ELSET, ELSET=Eall
Efaces
***********************************************************
** Element sets for materials and FEM element type (solid, shell, beam, fluid)
*ELSET,ELSET=MechanicalMaterialThickness
Efaces
***********************************************************
** constraints fixed node sets
** ConstraintFixed
*NSET,NSET=ConstraintFixed
1,
2,
5,
***********************************************************
** Materials
** see information about units at file end
** FreeCAD material name: Calculix-Steel
** MechanicalMaterial
*MATERIAL, NAME=MechanicalMaterial
*ELASTIC
210000,0.3
***********************************************************
** Sections
*SHELL SECTION, ELSET=MechanicalMaterialThickness, MATERIAL=MechanicalMaterial
1000
***********************************************************
** At least one step is needed to run an CalculiX analysis of FreeCAD
*STEP
*STATIC
***********************************************************
** Fixed Constraints
** ConstraintFixed
*BOUNDARY
ConstraintFixed,1
ConstraintFixed,2
ConstraintFixed,3
ConstraintFixed,4
ConstraintFixed,5
ConstraintFixed,6
***********************************************************
** constraints force node loads
*CLOAD
** ConstraintForce
** node loads on shape: CanileverPlate:Edge3
3,3,-2.2500000000000E+06
4,3,-2.2500000000000E+06
13,3,-4.5000000000000E+06
***********************************************************
** Outputs --> frd file
*NODE FILE, OUTPUT=2d
U
*EL FILE
S, E
** outputs --> dat file
** reaction forces for Constraint fixed
*NODE PRINT, NSET=ConstraintFixed, TOTALS=ONLY
RF
***********************************************************
*END STEP
***********************************************************
** CalculiX Input file
** written by --> FreeCAD 0.20.25347 (Git)
** written on --> Thu Jul 29 22:17:23 2021
** file name --> test.FCStd
** analysis name --> Analysis
**
**
***********************************************************
** About units:
** See ccx manual, ccx does not know about any unit.
** Golden rule: The user must make sure that the numbers he provides have consistent units.
** The user is the FreeCAD calculix writer module ;-)
**
** The unit system which is used at Guido Dhondt's company: mm, N, s, K
** Since Length and Mass are connected by Force, if Length is mm the Mass is in t to get N
** The following units are used to write to inp file:
**
** Length: mm (this includes the mesh geometry)
** Mass: t
** TimeSpan: s
** Temperature: K
**
** This leads to:
** Force: N
** Pressure: N/mm^2 == MPa (Young's Modulus has unit Pressure)
** Density: t/mm^3
** Gravity: mm/s^2
** Thermal conductivity: t*mm/K/s^3 == as W/m/K == kW/mm/K
** Specific Heat: mm^2/s^2/K = J/kg/K == kJ/t/K
**

151
src/Mod/Fem/femtest/data/calculix/ccx_cantilever_ele_quad8.inp Normal file
View File

@@ -0,0 +1,151 @@
** written by FreeCAD inp file writer for CalculiX,Abaqus meshes
** highest dimension mesh elements only.
** Nodes
*Node, NSET=Nall
1, 0, 500, 0
2, 0, 500, 1000
3, 8000, 500, 0
4, 8000, 500, 1000
5, 0, 500, 500
6, 0, 500, 250
7, 0, 500, 750
8, 4000, 500, 0
9, 2000, 500, 0
10, 6000, 500, 0
11, 8000, 500, 500
12, 8000, 500, 250
13, 8000, 500, 750
14, 4000, 500, 1000
15, 2000, 500, 1000
16, 6000, 500, 1000
17, 4000, 500, 500
18, 2000, 500, 500
19, 4000, 500, 750
20, 2000, 500, 750
21, 6000, 500, 500
22, 6000, 500, 750
23, 4000, 500, 250
24, 6000, 500, 250
25, 2000, 500, 250
** Face elements
*Element, TYPE=S8, ELSET=Efaces
1, 1, 5, 17, 8, 6, 18, 23, 9
2, 5, 2, 14, 17, 7, 15, 19, 18
3, 17, 14, 4, 11, 19, 16, 13, 21
4, 8, 17, 11, 3, 23, 21, 12, 10
** Define element set Eall
*ELSET, ELSET=Eall
Efaces
***********************************************************
** Element sets for materials and FEM element type (solid, shell, beam, fluid)
*ELSET,ELSET=MechanicalMaterialThickness
Efaces
***********************************************************
** constraints fixed node sets
** ConstraintFixed
*NSET,NSET=ConstraintFixed
1,
2,
5,
6,
7,
***********************************************************
** Materials
** see information about units at file end
** FreeCAD material name: Calculix-Steel
** MechanicalMaterial
*MATERIAL, NAME=MechanicalMaterial
*ELASTIC
210000,0.3
***********************************************************
** Sections
*SHELL SECTION, ELSET=MechanicalMaterialThickness, MATERIAL=MechanicalMaterial
1000
***********************************************************
** At least one step is needed to run an CalculiX analysis of FreeCAD
*STEP
*STATIC
***********************************************************
** Fixed Constraints
** ConstraintFixed
*BOUNDARY
ConstraintFixed,1
ConstraintFixed,2
ConstraintFixed,3
ConstraintFixed,4
ConstraintFixed,5
ConstraintFixed,6
***********************************************************
** constraints force node loads
*CLOAD
** ConstraintForce
** node loads on shape: CanileverPlate:Edge3
3,3,-7.5000000000000E+05
4,3,-7.5000000000000E+05
11,3,-1.5000000000000E+06
12,3,-3.0000000000000E+06
13,3,-3.0000000000000E+06
***********************************************************
** Outputs --> frd file
*NODE FILE, OUTPUT=2d
U
*EL FILE
S, E
** outputs --> dat file
** reaction forces for Constraint fixed
*NODE PRINT, NSET=ConstraintFixed, TOTALS=ONLY
RF
***********************************************************
*END STEP
***********************************************************
** CalculiX Input file
** written by --> FreeCAD 0.20.25343 (Git)
** written on --> Thu Jul 29 13:06:02 2021
** file name --> ccx_cantilever_ele_quad8.FCStd
** analysis name --> Analysis
**
**
***********************************************************
** About units:
** See ccx manual, ccx does not know about any unit.
** Golden rule: The user must make sure that the numbers he provides have consistent units.
** The user is the FreeCAD calculix writer module ;-)
**
** The unit system which is used at Guido Dhondt's company: mm, N, s, K
** Since Length and Mass are connected by Force, if Length is mm the Mass is in t to get N
** The following units are used to write to inp file:
**
** Length: mm (this includes the mesh geometry)
** Mass: t
** TimeSpan: s
** Temperature: K
**
** This leads to:
** Force: N
** Pressure: N/mm^2 == MPa (Young's Modulus has unit Pressure)
** Density: t/mm^3
** Gravity: mm/s^2
** Thermal conductivity: t*mm/K/s^3 == as W/m/K == kW/mm/K
** Specific Heat: mm^2/s^2/K = J/kg/K == kJ/t/K
**

277
src/Mod/Fem/femtest/data/calculix/ccx_cantilever_ele_seg2.inp Normal file
View File

@@ -0,0 +1,277 @@
** written by FreeCAD inp file writer for CalculiX,Abaqus meshes
** highest dimension mesh elements only.
** Nodes
*Node, NSET=Nall
1, 0, 500, 500
2, 8000, 500, 500
3, 148.1481481481, 500, 500
4, 296.2962962963, 500, 500
5, 444.4444444444, 500, 500
6, 592.5925925926, 500, 500
7, 740.7407407407, 500, 500
8, 888.8888888889, 500, 500
9, 1037.037037037, 500, 500
10, 1185.185185185, 500, 500
11, 1333.333333333, 500, 500
12, 1481.481481481, 500, 500
13, 1629.62962963, 500, 500
14, 1777.777777778, 500, 500
15, 1925.925925926, 500, 500
16, 2074.074074074, 500, 500
17, 2222.222222222, 500, 500
18, 2370.37037037, 500, 500
19, 2518.518518519, 500, 500
20, 2666.666666667, 500, 500
21, 2814.814814815, 500, 500
22, 2962.962962963, 500, 500
23, 3111.111111111, 500, 500
24, 3259.259259259, 500, 500
25, 3407.407407407, 500, 500
26, 3555.555555556, 500, 500
27, 3703.703703704, 500, 500
28, 3851.851851852, 500, 500
29, 4000, 500, 500
30, 4148.148148148, 500, 500
31, 4296.296296296, 500, 500
32, 4444.444444444, 500, 500
33, 4592.592592593, 500, 500
34, 4740.740740741, 500, 500
35, 4888.888888889, 500, 500
36, 5037.037037037, 500, 500
37, 5185.185185185, 500, 500
38, 5333.333333333, 500, 500
39, 5481.481481481, 500, 500
40, 5629.62962963, 500, 500
41, 5777.777777778, 500, 500
42, 5925.925925926, 500, 500
43, 6074.074074074, 500, 500
44, 6222.222222222, 500, 500
45, 6370.37037037, 500, 500
46, 6518.518518519, 500, 500
47, 6666.666666667, 500, 500
48, 6814.814814815, 500, 500
49, 6962.962962963, 500, 500
50, 7111.111111111, 500, 500
51, 7259.259259259, 500, 500
52, 7407.407407407, 500, 500
53, 7555.555555556, 500, 500
54, 7703.703703704, 500, 500
55, 7851.851851852, 500, 500
** Edge elements
*Element, TYPE=B31, ELSET=Eedges
1, 1, 3
2, 3, 4
3, 4, 5
4, 5, 6
5, 6, 7
6, 7, 8
7, 8, 9
8, 9, 10
9, 10, 11
10, 11, 12
11, 12, 13
12, 13, 14
13, 14, 15
14, 15, 16
15, 16, 17
16, 17, 18
17, 18, 19
18, 19, 20
19, 20, 21
20, 21, 22
21, 22, 23
22, 23, 24
23, 24, 25
24, 25, 26
25, 26, 27
26, 27, 28
27, 28, 29
28, 29, 30
29, 30, 31
30, 31, 32
31, 32, 33
32, 33, 34
33, 34, 35
34, 35, 36
35, 36, 37
36, 37, 38
37, 38, 39
38, 39, 40
39, 40, 41
40, 41, 42
41, 42, 43
42, 43, 44
43, 44, 45
44, 45, 46
45, 46, 47
46, 47, 48
47, 48, 49
48, 49, 50
49, 50, 51
50, 51, 52
51, 52, 53
52, 53, 54
53, 54, 55
54, 55, 2
** Define element set Eall
*ELSET, ELSET=Eall
Eedges
***********************************************************
** Element sets for materials and FEM element type (solid, shell, beam, fluid)
*ELSET,ELSET=M0B0RstdD0
1,
2,
3,
4,
5,
6,
7,
8,
9,
10,
11,
12,
13,
14,
15,
16,
17,
18,
19,
20,
21,
22,
23,
24,
25,
26,
27,
28,
29,
30,
31,
32,
33,
34,
35,
36,
37,
38,
39,
40,
41,
42,
43,
44,
45,
46,
47,
48,
49,
50,
51,
52,
53,
54,
***********************************************************
** constraints fixed node sets
** ConstraintFixed
*NSET,NSET=ConstraintFixed
1,
***********************************************************
** Materials
** see information about units at file end
** FreeCAD material name: Calculix-Steel
** MechanicalMaterial
*MATERIAL, NAME=MechanicalMaterial
*ELASTIC
210000,0.3
***********************************************************
** Sections
*BEAM SECTION, ELSET=M0B0RstdD0, MATERIAL=MechanicalMaterial, SECTION=RECT
1000,1000
-0, 1, 0
***********************************************************
** At least one step is needed to run an CalculiX analysis of FreeCAD
*STEP
*STATIC
***********************************************************
** Fixed Constraints
** ConstraintFixed
*BOUNDARY
ConstraintFixed,1
ConstraintFixed,2
ConstraintFixed,3
ConstraintFixed,4
ConstraintFixed,5
ConstraintFixed,6
***********************************************************
** constraints force node loads
*CLOAD
** ConstraintForce
** node load on shape: CantileverLine:Vertex2
2,3,-9.0000000000000E+06
***********************************************************
** Outputs --> frd file
*NODE FILE, OUTPUT=2d
U
*EL FILE
S, E
** outputs --> dat file
** reaction forces for Constraint fixed
*NODE PRINT, NSET=ConstraintFixed, TOTALS=ONLY
RF
***********************************************************
*END STEP
***********************************************************
** CalculiX Input file
** written by --> FreeCAD 0.20.25335 (Git)
** written on --> Wed Jul 28 17:01:46 2021
** file name --> test.FCStd
** analysis name --> Analysis
**
**
***********************************************************
** About units:
** See ccx manual, ccx does not know about any unit.
** Golden rule: The user must make sure that the numbers he provides have consistent units.
** The user is the FreeCAD calculix writer module ;-)
**
** The unit system which is used at Guido Dhondt's company: mm, N, s, K
** Since Length and Mass are connected by Force, if Length is mm the Mass is in t to get N
** The following units are used to write to inp file:
**
** Length: mm (this includes the mesh geometry)
** Mass: t
** TimeSpan: s
** Temperature: K
**
** This leads to:
** Force: N
** Pressure: N/mm^2 == MPa (Young's Modulus has unit Pressure)
** Density: t/mm^3
** Gravity: mm/s^2
** Thermal conductivity: t*mm/K/s^3 == as W/m/K == kW/mm/K
** Specific Heat: mm^2/s^2/K = J/kg/K == kJ/t/K
**

135
src/Mod/Fem/femtest/data/calculix/ccx_cantilever_ele_seg3.inp Normal file
View File

@@ -0,0 +1,135 @@
** written by FreeCAD inp file writer for CalculiX,Abaqus meshes
** highest dimension mesh elements only.
** Nodes
*Node, NSET=Nall
1, 0, 500, 500
2, 8000, 500, 500
3, 1600, 500, 500
4, 3200, 500, 500
5, 4800, 500, 500
6, 6400, 500, 500
7, 800, 500, 500
8, 2400, 500, 500
9, 4000, 500, 500
10, 5600, 500, 500
11, 7200, 500, 500
** Edge elements
*Element, TYPE=B32, ELSET=Eedges
1, 1, 7, 3
2, 3, 8, 4
3, 4, 9, 5
4, 5, 10, 6
5, 6, 11, 2
** Define element set Eall
*ELSET, ELSET=Eall
Eedges
***********************************************************
** Element sets for materials and FEM element type (solid, shell, beam, fluid)
*ELSET,ELSET=M0B0RstdD0
1,
2,
3,
4,
5,
***********************************************************
** constraints fixed node sets
** ConstraintFixed
*NSET,NSET=ConstraintFixed
1,
***********************************************************
** Materials
** see information about units at file end
** FreeCAD material name: Calculix-Steel
** MechanicalMaterial
*MATERIAL, NAME=MechanicalMaterial
*ELASTIC
210000,0.3
***********************************************************
** Sections
*BEAM SECTION, ELSET=M0B0RstdD0, MATERIAL=MechanicalMaterial, SECTION=RECT
1000,1000
-0, 1, 0
***********************************************************
** At least one step is needed to run an CalculiX analysis of FreeCAD
*STEP
*STATIC
***********************************************************
** Fixed Constraints
** ConstraintFixed
*BOUNDARY
ConstraintFixed,1
ConstraintFixed,2
ConstraintFixed,3
ConstraintFixed,4
ConstraintFixed,5
ConstraintFixed,6
***********************************************************
** constraints force node loads
*CLOAD
** ConstraintForce
** node load on shape: CantileverLine:Vertex2
2,3,-9.0000000000000E+06
***********************************************************
** Outputs --> frd file
*NODE FILE, OUTPUT=2d
U
*EL FILE
S, E
** outputs --> dat file
** reaction forces for Constraint fixed
*NODE PRINT, NSET=ConstraintFixed, TOTALS=ONLY
RF
***********************************************************
*END STEP
***********************************************************
** CalculiX Input file
** written by --> FreeCAD 0.20.25335 (Git)
** written on --> Wed Jul 28 16:59:39 2021
** file name --> test.FCStd
** analysis name --> Analysis
**
**
***********************************************************
** About units:
** See ccx manual, ccx does not know about any unit.
** Golden rule: The user must make sure that the numbers he provides have consistent units.
** The user is the FreeCAD calculix writer module ;-)
**
** The unit system which is used at Guido Dhondt's company: mm, N, s, K
** Since Length and Mass are connected by Force, if Length is mm the Mass is in t to get N
** The following units are used to write to inp file:
**
** Length: mm (this includes the mesh geometry)
** Mass: t
** TimeSpan: s
** Temperature: K
**
** This leads to:
** Force: N
** Pressure: N/mm^2 == MPa (Young's Modulus has unit Pressure)
** Density: t/mm^3
** Gravity: mm/s^2
** Thermal conductivity: t*mm/K/s^3 == as W/m/K == kW/mm/K
** Specific Heat: mm^2/s^2/K = J/kg/K == kJ/t/K
**

1642
src/Mod/Fem/femtest/data/calculix/ccx_cantilever_ele_tria3.inp Normal file
View File

File diff suppressed because it is too large Load Diff

369
src/Mod/Fem/femtest/data/calculix/ccx_cantilever_ele_tria6.inp Normal file
View File

@@ -0,0 +1,369 @@
** written by FreeCAD inp file writer for CalculiX,Abaqus meshes
** highest dimension mesh elements only.
** Nodes
*Node, NSET=Nall
1, 0, 500, 0
2, 0, 500, 1000
3, 8000, 500, 0
4, 8000, 500, 1000
5, 0, 500, 500
6, 0, 500, 250
7, 0, 500, 750
8, 500, 500, 0
9, 1000, 500, 0
10, 1500, 500, 0
11, 2000, 500, 0
12, 2500, 500, 0
13, 3000, 500, 0
14, 3500, 500, 0
15, 4000, 500, 0
16, 4500, 500, 0
17, 5000, 500, 0
18, 5500, 500, 0
19, 6000, 500, 0
20, 6500, 500, 0
21, 7000, 500, 0
22, 7500, 500, 0
23, 250, 500, 0
24, 750, 500, 0
25, 1250, 500, 0
26, 1750, 500, 0
27, 2250, 500, 0
28, 2750, 500, 0
29, 3250, 500, 0
30, 3750, 500, 0
31, 4250, 500, 0
32, 4750, 500, 0
33, 5250, 500, 0
34, 5750, 500, 0
35, 6250, 500, 0
36, 6750, 500, 0
37, 7250, 500, 0
38, 7750, 500, 0
39, 8000, 500, 500
40, 8000, 500, 250
41, 8000, 500, 750
42, 500, 500, 1000
43, 1000, 500, 1000
44, 1500, 500, 1000
45, 2000, 500, 1000
46, 2500, 500, 1000
47, 3000, 500, 1000
48, 3500, 500, 1000
49, 4000, 500, 1000
50, 4500, 500, 1000
51, 5000, 500, 1000
52, 5500, 500, 1000
53, 6000, 500, 1000
54, 6500, 500, 1000
55, 7000, 500, 1000
56, 7500, 500, 1000
57, 250, 500, 1000
58, 750, 500, 1000
59, 1250, 500, 1000
60, 1750, 500, 1000
61, 2250, 500, 1000
62, 2750, 500, 1000
63, 3250, 500, 1000
64, 3750, 500, 1000
65, 4250, 500, 1000
66, 4750, 500, 1000
67, 5250, 500, 1000
68, 5750, 500, 1000
69, 6250, 500, 1000
70, 6750, 500, 1000
71, 7250, 500, 1000
72, 7750, 500, 1000
73, 714.9234693878, 500, 484.056122449
74, 1750, 500, 500
75, 2750, 500, 500
76, 3750, 500, 500
77, 4750, 500, 500
78, 5750, 500, 500
79, 6756.944444444, 500, 489.5833333333
80, 7301.463293651, 500, 489.7073412698
81, 1244.153911565, 500, 497.3426870748
82, 2250, 500, 500
83, 3250, 500, 500
84, 4250, 500, 500
85, 5250, 500, 500
86, 6251.157407407, 500, 498.2638888889
87, 355.2776691685, 500, 632.3045669746
88, 7660.39819718, 500, 648.5965255879
89, 7699.764777976, 500, 310.4561768902
90, 314.0402277112, 500, 323.2721378847
91, 1122.076955782, 500, 748.6713435374
92, 1372.076955782, 500, 748.6713435374
93, 2125, 500, 750
94, 2375, 500, 750
95, 3125, 500, 750
96, 3375, 500, 750
97, 4125, 500, 750
98, 4375, 500, 750
99, 5125, 500, 750
100, 5375, 500, 750
101, 6125.578703704, 500, 749.1319444444
102, 6375.578703704, 500, 749.1319444444
103, 857.4617346939, 500, 742.0280612245
104, 979.5386904762, 500, 490.6994047619
105, 1497.076955782, 500, 498.6713435374
106, 1625, 500, 750
107, 1875, 500, 750
108, 2000, 500, 500
109, 2500, 500, 500
110, 2625, 500, 750
111, 2875, 500, 750
112, 3000, 500, 500
113, 3500, 500, 500
114, 3625, 500, 750
115, 3875, 500, 750
116, 4000, 500, 500
117, 4500, 500, 500
118, 4625, 500, 750
119, 4875, 500, 750
120, 5000, 500, 500
121, 5500, 500, 500
122, 5625, 500, 750
123, 5875, 500, 750
124, 6000.578703704, 500, 499.1319444444
125, 6504.050925926, 500, 493.9236111111
126, 6628.472222222, 500, 744.7916666667
127, 7150.731646825, 500, 244.8536706349
128, 7029.203869048, 500, 489.6453373016
129, 6878.472222222, 500, 244.7916666667
130, 6878.472222222, 500, 744.7916666667
131, 7150.731646825, 500, 744.8536706349
132, 7400.731646825, 500, 244.8536706349
133, 7400.731646825, 500, 744.8536706349
134, 857.4617346939, 500, 242.0280612245
135, 607.4617346939, 500, 242.0280612245
136, 1875, 500, 250
137, 1625, 500, 250
138, 2875, 500, 250
139, 2625, 500, 250
140, 3875, 500, 250
141, 3625, 500, 250
142, 4875, 500, 250
143, 4625, 500, 250
144, 5875, 500, 250
145, 5625, 500, 250
146, 6628.472222222, 500, 244.7916666667
147, 607.4617346939, 500, 742.0280612245
148, 1122.076955782, 500, 248.6713435374
149, 1372.076955782, 500, 248.6713435374
150, 2125, 500, 250
151, 2375, 500, 250
152, 3125, 500, 250
153, 3375, 500, 250
154, 4125, 500, 250
155, 4375, 500, 250
156, 5125, 500, 250
157, 5375, 500, 250
158, 6125.578703704, 500, 249.1319444444
159, 6375.578703704, 500, 249.1319444444
160, 514.4818485495, 500, 403.6641301669
161, 407.0201138556, 500, 161.6360689424
162, 177.6388345842, 500, 566.1522834873
163, 177.6388345842, 500, 816.1522834873
164, 427.6388345842, 500, 816.1522834873
165, 535.1005692781, 500, 558.1803447118
166, 7599.882388988, 500, 155.2280884451
167, 7500.614035813, 500, 400.08175908
168, 7830.19909859, 500, 824.2982627939
169, 7830.19909859, 500, 574.2982627939
170, 7580.19909859, 500, 824.2982627939
171, 7480.930745415, 500, 569.1519334289
172, 334.6589484399, 500, 477.7883524297
173, 7849.882388988, 500, 405.2280884451
174, 7849.882388988, 500, 155.2280884451
175, 157.0201138556, 500, 161.6360689424
176, 157.0201138556, 500, 411.6360689424
177, 7680.081487578, 500, 479.5263512391
** Face elements
*Element, TYPE=S6, ELSET=Efaces
37, 43, 81, 44, 91, 92, 59
38, 45, 82, 46, 93, 94, 61
39, 47, 83, 48, 95, 96, 63
40, 49, 84, 50, 97, 98, 65
41, 51, 85, 52, 99, 100, 67
42, 53, 86, 54, 101, 102, 69
43, 43, 73, 81, 103, 104, 91
44, 44, 81, 74, 92, 105, 106
45, 45, 74, 82, 107, 108, 93
46, 46, 82, 75, 94, 109, 110
47, 47, 75, 83, 111, 112, 95
48, 48, 83, 76, 96, 113, 114
49, 49, 76, 84, 115, 116, 97
50, 50, 84, 77, 98, 117, 118
51, 51, 77, 85, 119, 120, 99
52, 52, 85, 78, 100, 121, 122
53, 53, 78, 86, 123, 124, 101
54, 54, 86, 79, 102, 125, 126
55, 21, 80, 79, 127, 128, 129
56, 55, 79, 80, 130, 128, 131
57, 21, 22, 80, 37, 132, 127
58, 55, 80, 56, 131, 133, 71
59, 8, 9, 73, 24, 134, 135
60, 10, 11, 74, 26, 136, 137
61, 12, 13, 75, 28, 138, 139
62, 14, 15, 76, 30, 140, 141
63, 16, 17, 77, 32, 142, 143
64, 18, 19, 78, 34, 144, 145
65, 20, 21, 79, 36, 129, 146
66, 42, 73, 43, 147, 103, 58
67, 44, 74, 45, 106, 107, 60
68, 46, 75, 47, 110, 111, 62
69, 48, 76, 49, 114, 115, 64
70, 50, 77, 51, 118, 119, 66
71, 52, 78, 53, 122, 123, 68
72, 54, 79, 55, 126, 130, 70
73, 9, 81, 73, 148, 104, 134
74, 10, 74, 81, 137, 105, 149
75, 11, 82, 74, 150, 108, 136
76, 12, 75, 82, 139, 109, 151
77, 13, 83, 75, 152, 112, 138
78, 14, 76, 83, 141, 113, 153
79, 15, 84, 76, 154, 116, 140
80, 16, 77, 84, 143, 117, 155
81, 17, 85, 77, 156, 120, 142
82, 18, 78, 85, 145, 121, 157
83, 19, 86, 78, 158, 124, 144
84, 20, 79, 86, 146, 125, 159
85, 9, 10, 81, 25, 149, 148
86, 11, 12, 82, 27, 151, 150
87, 13, 14, 83, 29, 153, 152
88, 15, 16, 84, 31, 155, 154
89, 17, 18, 85, 33, 157, 156
90, 19, 20, 86, 35, 159, 158
91, 8, 73, 90, 135, 160, 161
92, 2, 5, 87, 7, 162, 163
93, 2, 87, 42, 163, 164, 57
94, 42, 87, 73, 164, 165, 147
95, 22, 89, 80, 166, 167, 132
96, 4, 88, 39, 168, 169, 41
97, 4, 56, 88, 72, 170, 168
98, 56, 80, 88, 133, 171, 170
99, 73, 87, 90, 165, 172, 160
100, 3, 39, 89, 40, 173, 174
101, 3, 89, 22, 174, 166, 38
102, 1, 90, 5, 175, 176, 6
103, 1, 8, 90, 23, 161, 175
104, 80, 89, 88, 167, 177, 171
105, 5, 90, 87, 176, 172, 162
106, 39, 88, 89, 169, 177, 173
** Define element set Eall
*ELSET, ELSET=Eall
Efaces
***********************************************************
** Element sets for materials and FEM element type (solid, shell, beam, fluid)
*ELSET,ELSET=MechanicalMaterialThickness
Efaces
***********************************************************
** constraints fixed node sets
** ConstraintFixed
*NSET,NSET=ConstraintFixed
1,
2,
5,
6,
7,
***********************************************************
** Materials
** see information about units at file end
** FreeCAD material name: Calculix-Steel
** MechanicalMaterial
*MATERIAL, NAME=MechanicalMaterial
*ELASTIC
210000,0.3
***********************************************************
** Sections
*SHELL SECTION, ELSET=MechanicalMaterialThickness, MATERIAL=MechanicalMaterial
1000
***********************************************************
** At least one step is needed to run an CalculiX analysis of FreeCAD
*STEP
*STATIC
***********************************************************
** Fixed Constraints
** ConstraintFixed
*BOUNDARY
ConstraintFixed,1
ConstraintFixed,2
ConstraintFixed,3
ConstraintFixed,4
ConstraintFixed,5
ConstraintFixed,6
***********************************************************
** constraints force node loads
*CLOAD
** ConstraintForce
** node loads on shape: CanileverPlate:Edge3
3,3,-7.5000000000000E+05
4,3,-7.5000000000000E+05
39,3,-1.5000000000000E+06
40,3,-3.0000000000000E+06
41,3,-3.0000000000000E+06
***********************************************************
** Outputs --> frd file
*NODE FILE, OUTPUT=2d
U
*EL FILE
S, E
** outputs --> dat file
** reaction forces for Constraint fixed
*NODE PRINT, NSET=ConstraintFixed, TOTALS=ONLY
RF
***********************************************************
*END STEP
***********************************************************
** CalculiX Input file
** written by --> FreeCAD 0.20.25335 (Git)
** written on --> Wed Jul 28 17:09:34 2021
** file name --> test.FCStd
** analysis name --> Analysis
**
**
***********************************************************
** About units:
** See ccx manual, ccx does not know about any unit.
** Golden rule: The user must make sure that the numbers he provides have consistent units.
** The user is the FreeCAD calculix writer module ;-)
**
** The unit system which is used at Guido Dhondt's company: mm, N, s, K
** Since Length and Mass are connected by Force, if Length is mm the Mass is in t to get N
** The following units are used to write to inp file:
**
** Length: mm (this includes the mesh geometry)
** Mass: t
** TimeSpan: s
** Temperature: K
**
** This leads to:
** Force: N
** Pressure: N/mm^2 == MPa (Young's Modulus has unit Pressure)
** Density: t/mm^3
** Gravity: mm/s^2
** Thermal conductivity: t*mm/K/s^3 == as W/m/K == kW/mm/K
** Specific Heat: mm^2/s^2/K = J/kg/K == kJ/t/K
**

101
src/Mod/Fem/femtest/data/mystran/.bdf Normal file
View File

@@ -0,0 +1,101 @@
$pyNastran: version=msc
$pyNastran: punch=False
$pyNastran: encoding=utf-8
$pyNastran: nnodes=36
$pyNastran: nelements=25
$EXECUTIVE CONTROL DECK
SOL 101
CEND
$CASE CONTROL DECK
ECHO = NONE
TITLE = pyNastran plate example for Mystran
SUBCASE 1
DISPLACEMENT(SORT1,REAL) = ALL
LOAD = 2
SPC = 2
SPCFORCES(SORT1,REAL) = ALL
STRESS(SORT1,REAL,VONMISES,BILIN) = ALL
SUBTITLE = Default
BEGIN BULK
$PARAMS
PARAM POST -1
PARAM PRTMAXIM YES
$NODES
GRID 1 0. 0. 0.
GRID 2 2. 0. 0.
GRID 3 4. 0. 0.
GRID 4 6. 0. 0.
GRID 5 8. 0. 0.
GRID 6 10. 0. 0.
GRID 7 0. 2. 0.
GRID 8 2. 2. 0.
GRID 9 4. 2. 0.
GRID 10 6. 2. 0.
GRID 11 8. 2. 0.
GRID 12 10. 2. 0.
GRID 13 0. 4. 0.
GRID 14 2. 4. 0.
GRID 15 4. 4. 0.
GRID 16 6. 4. 0.
GRID 17 8. 4. 0.
GRID 18 10. 4. 0.
GRID 19 0. 6. 0.
GRID 20 2. 6. 0.
GRID 21 4. 6. 0.
GRID 22 6. 6. 0.
GRID 23 8. 6. 0.
GRID 24 10. 6. 0.
GRID 25 0. 8. 0.
GRID 26 2. 8. 0.
GRID 27 4. 8. 0.
GRID 28 6. 8. 0.
GRID 29 8. 8. 0.
GRID 30 10. 8. 0.
GRID 31 0. 10. 0.
GRID 32 2. 10. 0.
GRID 33 4. 10. 0.
GRID 34 6. 10. 0.
GRID 35 8. 10. 0.
GRID 36 10. 10. 0.
$ELEMENTS
CQUAD4 1 1 1 2 8 7
CQUAD4 2 1 2 3 9 8
CQUAD4 3 1 3 4 10 9
CQUAD4 4 1 4 5 11 10
CQUAD4 5 1 5 6 12 11
CQUAD4 6 1 7 8 14 13
CQUAD4 7 1 8 9 15 14
CQUAD4 8 1 9 10 16 15
CQUAD4 9 1 10 11 17 16
CQUAD4 10 1 11 12 18 17
CQUAD4 11 1 13 14 20 19
CQUAD4 12 1 14 15 21 20
CQUAD4 13 1 15 16 22 21
CQUAD4 14 1 16 17 23 22
CQUAD4 15 1 17 18 24 23
CQUAD4 16 1 19 20 26 25
CQUAD4 17 1 20 21 27 26
CQUAD4 18 1 21 22 28 27
CQUAD4 19 1 22 23 29 28
CQUAD4 20 1 23 24 30 29
CQUAD4 21 1 25 26 32 31
CQUAD4 22 1 26 27 33 32
CQUAD4 23 1 27 28 34 33
CQUAD4 24 1 28 29 35 34
CQUAD4 25 1 29 30 36 35
$PROPERTIES
PSHELL 1 1 .3 1 1
$MATERIALS
MAT1 1 210000. .3
$LOADS
LOAD 2 1. 1. 1
FORCE 1 6 100. 1. 0. 0.
FORCE 1 12 100. 1. 0. 0.
FORCE 1 18 100. 1. 0. 0.
FORCE 1 24 100. 1. 0. 0.
FORCE 1 30 100. 1. 0. 0.
FORCE 1 36 100. 1. 0. 0.
$SPCs
SPCADD 2 1
SPC1 1 123456 1 7 13 19 25 31
ENDDATA

0
src/Mod/Fem/femtest/data/mystran/__init__.py Normal file
View File

79
src/Mod/Fem/femtest/data/mystran/ccx_cantilever_ele_quad4.bdf Normal file
View File

@@ -0,0 +1,79 @@
$pyNastran: version=msc
$pyNastran: punch=False
$pyNastran: encoding=utf-8
$pyNastran: nnodes=27
$pyNastran: nelements=16
$EXECUTIVE CONTROL DECK
SOL 101
CEND
$CASE CONTROL DECK
ECHO = NONE
TITLE = pyNastran for generating solverinput for for Mystran
SUBCASE 1
DISPLACEMENT(SORT1,REAL) = ALL
LOAD = 1
SPC = 1
SPCFORCES(SORT1,REAL) = ALL
STRESS(SORT1,REAL,VONMISES,BILIN) = ALL
SUBTITLE = Default
BEGIN BULK
$PARAMS
PARAM POST -1
$NODES
GRID 1 0. 500. 0.
GRID 2 0. 500. 1000.
GRID 3 8000. 500. 0.
GRID 4 8000. 500. 1000.
GRID 5 0. 500. 500.
GRID 6 2000. 500. 0.
GRID 7 4000. 500. 0.
GRID 8 6000. 500. 0.
GRID 9 1000. 500. 0.
GRID 10 7000. 500. 0.
GRID 11 3000. 500. 0.
GRID 12 5000. 500. 0.
GRID 13 8000. 500. 500.
GRID 14 2000. 500. 1000.
GRID 15 4000. 500. 1000.
GRID 16 6000. 500. 1000.
GRID 17 1000. 500. 1000.
GRID 18 7000. 500. 1000.
GRID 19 3000. 500. 1000.
GRID 20 5000. 500. 1000.
GRID 21 2000. 500. 500.
GRID 22 1000. 500. 500.
GRID 23 6000. 500. 500.
GRID 24 7000. 500. 500.
GRID 25 4000. 500. 500.
GRID 26 3000. 500. 500.
GRID 27 5000. 500. 500.
$ELEMENTS
CQUAD4 21 1 14 17 22 21
CQUAD4 22 1 6 21 22 9
CQUAD4 23 1 1 9 22 5
CQUAD4 24 1 2 5 22 17
CQUAD4 25 1 4 18 24 13
CQUAD4 26 1 3 13 24 10
CQUAD4 27 1 8 10 24 23
CQUAD4 28 1 16 23 24 18
CQUAD4 29 1 7 25 26 11
CQUAD4 30 1 6 11 26 21
CQUAD4 31 1 14 21 26 19
CQUAD4 32 1 15 19 26 25
CQUAD4 33 1 16 20 27 23
CQUAD4 34 1 8 23 27 12
CQUAD4 35 1 7 12 27 25
CQUAD4 36 1 15 25 27 20
$PROPERTIES
PSHELL 1 1 1000. 1 1
$MATERIALS
MAT1 1 210000. .3
$LOADS
LOAD 1 1. 1. 2
FORCE 2 3 2250000. 0.-2.22-16 -1.
FORCE 2 4 2250000. 0.-2.22-16 -1.
FORCE 2 13 4500000. 0.-2.22-16 -1.
$SPCs
SPCADD 1 2
SPC1 2 123456 1 2 5
ENDDATA

144
src/Mod/Fem/femtest/data/mystran/ccx_cantilever_ele_seg2.bdf Normal file
View File

@@ -0,0 +1,144 @@
$pyNastran: version=msc
$pyNastran: punch=False
$pyNastran: encoding=utf-8
$pyNastran: nnodes=55
$pyNastran: nelements=54
$EXECUTIVE CONTROL DECK
SOL 101
CEND
$CASE CONTROL DECK
ECHO = NONE
TITLE = pyNastran for generating solverinput for for Mystran
SUBCASE 1
DISPLACEMENT(SORT1,REAL) = ALL
LOAD = 1
SPC = 1
SPCFORCES(SORT1,REAL) = ALL
STRESS(SORT1,REAL,VONMISES,BILIN) = ALL
SUBTITLE = Default
BEGIN BULK
$PARAMS
PARAM POST -1
$NODES
GRID 1 0. 500. 500.
GRID 2 8000. 500. 500.
GRID 3 148.1481 500. 500.
GRID 4 296.2963 500. 500.
GRID 5 444.4444 500. 500.
GRID 6 592.5926 500. 500.
GRID 7 740.7407 500. 500.
GRID 8 888.8889 500. 500.
GRID 9 1037.037 500. 500.
GRID 10 1185.185 500. 500.
GRID 11 1333.333 500. 500.
GRID 12 1481.481 500. 500.
GRID 13 1629.63 500. 500.
GRID 14 1777.778 500. 500.
GRID 15 1925.926 500. 500.
GRID 16 2074.074 500. 500.
GRID 17 2222.222 500. 500.
GRID 18 2370.37 500. 500.
GRID 19 2518.519 500. 500.
GRID 20 2666.667 500. 500.
GRID 21 2814.815 500. 500.
GRID 22 2962.963 500. 500.
GRID 23 3111.111 500. 500.
GRID 24 3259.259 500. 500.
GRID 25 3407.407 500. 500.
GRID 26 3555.556 500. 500.
GRID 27 3703.704 500. 500.
GRID 28 3851.852 500. 500.
GRID 29 4000. 500. 500.
GRID 30 4148.148 500. 500.
GRID 31 4296.296 500. 500.
GRID 32 4444.444 500. 500.
GRID 33 4592.593 500. 500.
GRID 34 4740.741 500. 500.
GRID 35 4888.889 500. 500.
GRID 36 5037.037 500. 500.
GRID 37 5185.185 500. 500.
GRID 38 5333.333 500. 500.
GRID 39 5481.481 500. 500.
GRID 40 5629.63 500. 500.
GRID 41 5777.778 500. 500.
GRID 42 5925.926 500. 500.
GRID 43 6074.074 500. 500.
GRID 44 6222.222 500. 500.
GRID 45 6370.37 500. 500.
GRID 46 6518.519 500. 500.
GRID 47 6666.667 500. 500.
GRID 48 6814.815 500. 500.
GRID 49 6962.963 500. 500.
GRID 50 7111.111 500. 500.
GRID 51 7259.259 500. 500.
GRID 52 7407.407 500. 500.
GRID 53 7555.556 500. 500.
GRID 54 7703.704 500. 500.
GRID 55 7851.852 500. 500.
$ELEMENTS
CBAR 1 1 1 3 0. 0. 1.
CBAR 2 1 3 4 0. 0. 1.
CBAR 3 1 4 5 0. 0. 1.
CBAR 4 1 5 6 0. 0. 1.
CBAR 5 1 6 7 0. 0. 1.
CBAR 6 1 7 8 0. 0. 1.
CBAR 7 1 8 9 0. 0. 1.
CBAR 8 1 9 10 0. 0. 1.
CBAR 9 1 10 11 0. 0. 1.
CBAR 10 1 11 12 0. 0. 1.
CBAR 11 1 12 13 0. 0. 1.
CBAR 12 1 13 14 0. 0. 1.
CBAR 13 1 14 15 0. 0. 1.
CBAR 14 1 15 16 0. 0. 1.
CBAR 15 1 16 17 0. 0. 1.
CBAR 16 1 17 18 0. 0. 1.
CBAR 17 1 18 19 0. 0. 1.
CBAR 18 1 19 20 0. 0. 1.
CBAR 19 1 20 21 0. 0. 1.
CBAR 20 1 21 22 0. 0. 1.
CBAR 21 1 22 23 0. 0. 1.
CBAR 22 1 23 24 0. 0. 1.
CBAR 23 1 24 25 0. 0. 1.
CBAR 24 1 25 26 0. 0. 1.
CBAR 25 1 26 27 0. 0. 1.
CBAR 26 1 27 28 0. 0. 1.
CBAR 27 1 28 29 0. 0. 1.
CBAR 28 1 29 30 0. 0. 1.
CBAR 29 1 30 31 0. 0. 1.
CBAR 30 1 31 32 0. 0. 1.
CBAR 31 1 32 33 0. 0. 1.
CBAR 32 1 33 34 0. 0. 1.
CBAR 33 1 34 35 0. 0. 1.
CBAR 34 1 35 36 0. 0. 1.
CBAR 35 1 36 37 0. 0. 1.
CBAR 36 1 37 38 0. 0. 1.
CBAR 37 1 38 39 0. 0. 1.
CBAR 38 1 39 40 0. 0. 1.
CBAR 39 1 40 41 0. 0. 1.
CBAR 40 1 41 42 0. 0. 1.
CBAR 41 1 42 43 0. 0. 1.
CBAR 42 1 43 44 0. 0. 1.
CBAR 43 1 44 45 0. 0. 1.
CBAR 44 1 45 46 0. 0. 1.
CBAR 45 1 46 47 0. 0. 1.
CBAR 46 1 47 48 0. 0. 1.
CBAR 47 1 48 49 0. 0. 1.
CBAR 48 1 49 50 0. 0. 1.
CBAR 49 1 50 51 0. 0. 1.
CBAR 50 1 51 52 0. 0. 1.
CBAR 51 1 52 53 0. 0. 1.
CBAR 52 1 53 54 0. 0. 1.
CBAR 53 1 54 55 0. 0. 1.
CBAR 54 1 55 2 0. 0. 1.
$PROPERTIES
PBARL 1 1 BAR
1000. 1000. 0.
$MATERIALS
MAT1 1 210000. .3
$LOADS
LOAD 1 1. 1. 2
FORCE 2 2 9000000. 0. 0. -1.
$SPCs
SPCADD 1 2
SPC1 2 123456 1
ENDDATA

1556
src/Mod/Fem/femtest/data/mystran/ccx_cantilever_ele_tria3.bdf Normal file
View File

File diff suppressed because it is too large Load Diff

433
src/Mod/Fem/femtest/data/mystran/ccx_cantilever_faceload.bdf Normal file
View File

@@ -0,0 +1,433 @@
$pyNastran: version=msc
$pyNastran: punch=False
$pyNastran: encoding=utf-8
$pyNastran: nnodes=228
$pyNastran: nelements=79
$EXECUTIVE CONTROL DECK
SOL 101
CEND
$CASE CONTROL DECK
ECHO = NONE
TITLE = pyNastran for generating solverinput for for Mystran
SUBCASE 1
DISPLACEMENT(SORT1,REAL) = ALL
LOAD = 1
SPC = 1
SPCFORCES(SORT1,REAL) = ALL
STRESS(SORT1,REAL,VONMISES,BILIN) = ALL
SUBTITLE = Default
BEGIN BULK
$PARAMS
PARAM POST -1
$NODES
GRID 1 8000. 1000. 0.
GRID 2 8000. 1000. 1000.
GRID 3 8000. 0. 0.
GRID 4 8000. 0. 1000.
GRID 5 0. 1000. 0.
GRID 6 0. 1000. 1000.
GRID 7 0. 0. 0.
GRID 8 0. 0. 1000.
GRID 9 728. 1000. 1000.
GRID 10 1456. 1000. 1000.
GRID 11 2184. 1000. 1000.
GRID 12 2912. 1000. 1000.
GRID 13 3640. 1000. 1000.
GRID 14 4368. 1000. 1000.
GRID 15 5096. 1000. 1000.
GRID 16 5824. 1000. 1000.
GRID 17 6552. 1000. 1000.
GRID 18 7280. 1000. 1000.
GRID 19 728. 0. 1000.
GRID 20 1456. 0. 1000.
GRID 21 2184. 0. 1000.
GRID 22 2912. 0. 1000.
GRID 23 3640. 0. 1000.
GRID 24 4368. 0. 1000.
GRID 25 5096. 0. 1000.
GRID 26 5824. 0. 1000.
GRID 27 6552. 0. 1000.
GRID 28 7280. 0. 1000.
GRID 29 728. 1000. 0.
GRID 30 1456. 1000. 0.
GRID 31 2184. 1000. 0.
GRID 32 2912. 1000. 0.
GRID 33 3640. 1000. 0.
GRID 34 4368. 1000. 0.
GRID 35 5096. 1000. 0.
GRID 36 5824. 1000. 0.
GRID 37 6552. 1000. 0.
GRID 38 7280. 1000. 0.
GRID 39 728. 0. 0.
GRID 40 1456. 0. 0.
GRID 41 2184. 0. 0.
GRID 42 2912. 0. 0.
GRID 43 3640. 0. 0.
GRID 44 4368. 0. 0.
GRID 45 5096. 0. 0.
GRID 46 5824. 0. 0.
GRID 47 6552. 0. 0.
GRID 48 7280. 0. 0.
GRID 49 8000. 500. 500.
GRID 50 0. 500. 500.
GRID 51 4732. 500. 500.
GRID 52 0. 500. 1000.
GRID 53 364. 1000. 1000.
GRID 54 1092. 1000. 1000.
GRID 55 1820. 1000. 1000.
GRID 56 2548. 1000. 1000.
GRID 57 3276. 1000. 1000.
GRID 58 4004. 1000. 1000.
GRID 59 4732. 1000. 1000.
GRID 60 5460. 1000. 1000.
GRID 61 6188. 1000. 1000.
GRID 62 6916. 1000. 1000.
GRID 63 7640. 1000. 1000.
GRID 64 8000. 500. 1000.
GRID 65 364. 0. 1000.
GRID 66 1092. 0. 1000.
GRID 67 1820. 0. 1000.
GRID 68 2548. 0. 1000.
GRID 69 3276. 0. 1000.
GRID 70 4004. 0. 1000.
GRID 71 4732. 0. 1000.
GRID 72 5460. 0. 1000.
GRID 73 6188. 0. 1000.
GRID 74 6916. 0. 1000.
GRID 75 7640. 0. 1000.
GRID 76 0. 500. 0.
GRID 77 364. 1000. 0.
GRID 78 1092. 1000. 0.
GRID 79 1820. 1000. 0.
GRID 80 2548. 1000. 0.
GRID 81 3276. 1000. 0.
GRID 82 4004. 1000. 0.
GRID 83 4732. 1000. 0.
GRID 84 5460. 1000. 0.
GRID 85 6188. 1000. 0.
GRID 86 6916. 1000. 0.
GRID 87 7640. 1000. 0.
GRID 88 8000. 500. 0.
GRID 89 364. 0. 0.
GRID 90 1092. 0. 0.
GRID 91 1820. 0. 0.
GRID 92 2548. 0. 0.
GRID 93 3276. 0. 0.
GRID 94 4004. 0. 0.
GRID 95 4732. 0. 0.
GRID 96 5460. 0. 0.
GRID 97 6188. 0. 0.
GRID 98 6916. 0. 0.
GRID 99 7640. 0. 0.
GRID 100 8000. 1000. 500.
GRID 101 0. 1000. 500.
GRID 102 8000. 0. 500.
GRID 103 0. 0. 500.
GRID 104 364. 500. 1000.
GRID 105 728. 500. 1000.
GRID 106 1092. 500. 1000.
GRID 107 1456. 500. 1000.
GRID 108 1820. 500. 1000.
GRID 109 2184. 500. 1000.
GRID 110 2548. 500. 1000.
GRID 111 3276. 500. 1000.
GRID 112 3640. 500. 1000.
GRID 113 4004. 500. 1000.
GRID 114 4368. 500. 1000.
GRID 115 4732. 500. 1000.
GRID 116 5096. 500. 1000.
GRID 117 5460. 500. 1000.
GRID 118 5824. 500. 1000.
GRID 119 6188. 500. 1000.
GRID 120 6552. 500. 1000.
GRID 121 6916. 500. 1000.
GRID 122 7640. 500. 1000.
GRID 123 2912. 500. 1000.
GRID 124 7280. 500. 1000.
GRID 125 364. 500. 0.
GRID 126 1092. 500. 0.
GRID 127 728. 500. 0.
GRID 128 1820. 500. 0.
GRID 129 1456. 500. 0.
GRID 130 2548. 500. 0.
GRID 131 2184. 500. 0.
GRID 132 3640. 500. 0.
GRID 133 3276. 500. 0.
GRID 134 4004. 500. 0.
GRID 135 5096. 500. 0.
GRID 136 4732. 500. 0.
GRID 137 5460. 500. 0.
GRID 138 6188. 500. 0.
GRID 139 5824. 500. 0.
GRID 140 6916. 500. 0.
GRID 141 6552. 500. 0.
GRID 142 7640. 500. 0.
GRID 143 2912. 500. 0.
GRID 144 4368. 500. 0.
GRID 145 7280. 500. 0.
GRID 146 364. 1000. 500.
GRID 147 728. 1000. 500.
GRID 148 1092. 1000. 500.
GRID 149 1456. 1000. 500.
GRID 150 1820. 1000. 500.
GRID 151 2184. 1000. 500.
GRID 152 2548. 1000. 500.
GRID 153 3276. 1000. 500.
GRID 154 3640. 1000. 500.
GRID 155 4004. 1000. 500.
GRID 156 4368. 1000. 500.
GRID 157 4732. 1000. 500.
GRID 158 5096. 1000. 500.
GRID 159 5460. 1000. 500.
GRID 160 5824. 1000. 500.
GRID 161 6188. 1000. 500.
GRID 162 6552. 1000. 500.
GRID 163 6916. 1000. 500.
GRID 164 7640. 1000. 500.
GRID 165 2912. 1000. 500.
GRID 166 7280. 1000. 500.
GRID 167 364. 0. 500.
GRID 168 1092. 0. 500.
GRID 169 728. 0. 500.
GRID 170 1820. 0. 500.
GRID 171 1456. 0. 500.
GRID 172 2548. 0. 500.
GRID 173 2184. 0. 500.
GRID 174 3640. 0. 500.
GRID 175 3276. 0. 500.
GRID 176 4004. 0. 500.
GRID 177 5096. 0. 500.
GRID 178 4732. 0. 500.
GRID 179 5460. 0. 500.
GRID 180 6188. 0. 500.
GRID 181 5824. 0. 500.
GRID 182 6916. 0. 500.
GRID 183 6552. 0. 500.
GRID 184 7640. 0. 500.
GRID 185 2912. 0. 500.
GRID 186 4368. 0. 500.
GRID 187 7280. 0. 500.
GRID 188 8000. 250. 250.
GRID 189 8000. 250. 750.
GRID 190 8000. 750. 750.
GRID 191 8000. 750. 250.
GRID 192 0. 250. 750.
GRID 193 0. 250. 250.
GRID 194 0. 750. 250.
GRID 195 0. 750. 750.
GRID 196 1456. 500. 500.
GRID 197 6552. 500. 500.
GRID 198 6916. 500. 500.
GRID 199 2184. 500. 500.
GRID 200 2548. 500. 500.
GRID 201 2912. 500. 500.
GRID 202 1820. 500. 500.
GRID 203 7640. 750. 250.
GRID 204 7640. 750. 750.
GRID 205 7280. 500. 500.
GRID 206 7640. 250. 250.
GRID 207 5460. 500. 500.
GRID 208 5096. 500. 500.
GRID 209 6188. 500. 500.
GRID 210 5824. 500. 500.
GRID 211 364. 750. 250.
GRID 212 364. 750. 750.
GRID 213 364. 250. 250.
GRID 214 1092. 500. 500.
GRID 215 728. 500. 500.
GRID 216 364. 250. 750.
GRID 217 4550. 250. 250.
GRID 218 4550. 750. 250.
GRID 219 4550. 750. 750.
GRID 220 4368. 500. 500.
GRID 221 4550. 250. 750.
GRID 222 4914. 250. 750.
GRID 223 4914. 750. 750.
GRID 224 3276. 500. 500.
GRID 225 3640. 500. 500.
GRID 226 4004. 500. 500.
GRID 227 4914. 750. 250.
GRID 228 4914. 250. 250.
$ELEMENTS
CTETRA 149 1 40 10 19 20 196 106
168 171 107 66
CTETRA 150 1 10 30 31 40 149 79
150 196 129 128
CTETRA 151 1 38 18 17 47 166 62
163 140 198 197
CTETRA 152 1 32 11 41 12 152 199
130 165 56 200
CTETRA 153 1 12 41 32 42 200 130
165 201 92 143
CTETRA 154 1 42 12 21 22 201 110
172 185 123 68
CTETRA 155 1 10 40 31 11 196 128
150 55 202 151
CTETRA 156 1 11 41 12 21 199 200
56 109 173 110
CTETRA 157 1 20 11 41 40 108 199
170 171 202 91
CTETRA 158 1 20 21 41 11 67 173
170 108 109 199
CTETRA 159 1 32 41 11 31 130 199
152 80 131 151
CTETRA 160 1 11 20 10 40 108 107
55 202 171 196
CTETRA 161 1 38 47 17 37 140 197
163 86 141 162
CTETRA 162 1 38 49 18 48 203 204
166 145 206 205
CTETRA 163 1 46 45 15 36 96 208
207 139 137 159
CTETRA 164 1 47 46 16 37 97 210
209 141 138 161
CTETRA 165 1 18 2 4 49 63 64
122 204 190 189
CTETRA 166 1 18 27 17 47 121 120
62 198 183 197
CTETRA 167 1 38 48 18 47 145 205
166 140 98 198
CTETRA 168 1 12 42 23 22 201 175
111 123 185 69
CTETRA 169 1 26 27 47 17 73 183
180 119 120 197
CTETRA 170 1 50 9 29 6 212 147
211 195 53 146
CTETRA 171 1 27 18 48 47 121 205
182 183 198 98
CTETRA 172 1 8 39 7 50 167 89
103 192 213 193
CTETRA 173 1 40 39 9 30 90 215
214 129 126 148
CTETRA 174 1 42 41 21 12 92 173
172 201 200 110
CTETRA 175 1 50 29 9 39 211 147
212 213 127 215
CTETRA 176 1 29 50 7 39 211 193
125 127 213 89
CTETRA 177 1 31 41 11 40 131 199
151 128 91 202
CTETRA 178 1 47 37 16 17 141 161
209 197 162 61
CTETRA 179 1 40 30 9 10 129 148
214 196 149 54
CTETRA 180 1 2 49 38 1 190 203
164 100 191 87
CTETRA 181 1 2 18 38 49 63 166
164 190 204 203
CTETRA 182 1 48 38 49 3 145 203
206 99 142 188
CTETRA 183 1 38 1 49 3 87 191
203 142 88 188
CTETRA 184 1 49 3 4 48 188 102
189 206 99 184
CTETRA 185 1 28 4 48 18 75 184
187 124 122 205
CTETRA 186 1 49 4 18 48 189 122
204 206 184 205
CTETRA 187 1 7 5 50 29 76 194
193 125 77 211
CTETRA 188 1 50 5 6 29 194 101
195 211 77 146
CTETRA 189 1 50 19 9 6 216 105
212 195 104 53
CTETRA 190 1 50 9 19 39 212 105
216 213 215 169
CTETRA 191 1 50 8 19 6 192 65
216 195 52 104
CTETRA 192 1 40 10 9 19 196 54
214 168 106 105
CTETRA 193 1 51 34 44 14 218 144
217 219 156 220
CTETRA 194 1 51 44 24 14 217 186
221 219 220 114
CTETRA 195 1 25 24 15 51 71 115
116 222 221 223
CTETRA 196 1 43 32 12 13 133 165
224 225 153 57
CTETRA 197 1 43 23 12 42 174 111
224 93 175 201
CTETRA 198 1 43 13 12 23 225 57
224 174 112 111
CTETRA 199 1 43 42 12 32 93 201
224 133 143 165
CTETRA 200 1 34 14 13 44 156 58
155 144 220 226
CTETRA 201 1 14 15 24 51 59 115
114 219 223 221
CTETRA 202 1 23 14 24 44 113 114
70 176 220 186
CTETRA 203 1 33 13 32 43 154 153
81 132 225 133
CTETRA 204 1 34 13 33 43 155 154
82 134 225 132
CTETRA 205 1 35 15 14 51 158 59
157 227 223 219
CTETRA 206 1 25 15 45 51 116 208
177 222 223 228
CTETRA 207 1 44 13 43 23 226 225
94 176 112 174
CTETRA 208 1 35 14 34 51 157 156
83 227 219 218
CTETRA 209 1 46 36 15 16 139 159
207 210 160 60
CTETRA 210 1 36 15 35 45 159 158
84 137 208 135
CTETRA 211 1 37 16 36 46 161 160
85 138 210 139
CTETRA 212 1 25 16 26 46 117 118
72 179 210 181
CTETRA 213 1 47 17 16 26 197 61
209 180 119 118
CTETRA 214 1 47 26 16 46 180 118
209 97 181 210
CTETRA 215 1 27 18 28 48 121 124
74 182 205 187
CTETRA 216 1 35 51 34 45 227 218
83 135 228 136
CTETRA 217 1 13 44 14 23 226 220
58 112 176 113
CTETRA 218 1 44 24 25 51 186 71
178 217 221 222
CTETRA 219 1 44 25 45 51 178 177
95 217 222 228
CTETRA 220 1 46 16 15 25 210 60
207 179 117 116
CTETRA 221 1 46 25 15 45 179 116
207 96 177 208
CTETRA 222 1 50 19 8 39 216 65
192 213 169 167
CTETRA 223 1 35 15 51 45 158 223
227 135 208 228
CTETRA 224 1 34 43 44 13 134 94
144 155 225 226
CTETRA 225 1 51 45 44 34 228 95
217 218 136 144
CTETRA 226 1 9 39 29 30 215 127
147 148 126 78
CTETRA 227 1 40 19 9 39 168 105
214 90 169 215
$PROPERTIES
PSOLID 1 1
$MATERIALS
MAT1 1 210000. .3
$LOADS
LOAD 1 1. 1. 2
FORCE 2 1 0. 0. 0. -1.
FORCE 2 2 0. 0. 0. -1.
FORCE 2 3 0. 0. 0. -1.
FORCE 2 4 0. 0. 0. -1.
FORCE 2 49 0. 0. 0. -1.
FORCE 2 64 750000. 0. 0. -1.
FORCE 2 88 750000. 0. 0. -1.
FORCE 2 100 750000. 0. 0. -1.
FORCE 2 102 750000. 0. 0. -1.
FORCE 2 188 1500000. 0. 0. -1.
FORCE 2 189 1500000. 0. 0. -1.
FORCE 2 190 1500000. 0. 0. -1.
FORCE 2 191 1500000. 0. 0. -1.
$SPCs
SPCADD 1 2
SPC1 2 123456 5 6 7 8 50 52
76 101 103 192 193 194 195
ENDDATA

424
src/Mod/Fem/femtest/data/mystran/ccx_cantilever_nodeload.bdf Normal file
View File

@@ -0,0 +1,424 @@
$pyNastran: version=msc
$pyNastran: punch=False
$pyNastran: encoding=utf-8
$pyNastran: nnodes=228
$pyNastran: nelements=79
$EXECUTIVE CONTROL DECK
SOL 101
CEND
$CASE CONTROL DECK
ECHO = NONE
TITLE = pyNastran for generating solverinput for for Mystran
SUBCASE 1
DISPLACEMENT(SORT1,REAL) = ALL
LOAD = 1
SPC = 1
SPCFORCES(SORT1,REAL) = ALL
STRESS(SORT1,REAL,VONMISES,BILIN) = ALL
SUBTITLE = Default
BEGIN BULK
$PARAMS
PARAM POST -1
$NODES
GRID 1 8000. 1000. 0.
GRID 2 8000. 1000. 1000.
GRID 3 8000. 0. 0.
GRID 4 8000. 0. 1000.
GRID 5 0. 1000. 0.
GRID 6 0. 1000. 1000.
GRID 7 0. 0. 0.
GRID 8 0. 0. 1000.
GRID 9 728. 1000. 1000.
GRID 10 1456. 1000. 1000.
GRID 11 2184. 1000. 1000.
GRID 12 2912. 1000. 1000.
GRID 13 3640. 1000. 1000.
GRID 14 4368. 1000. 1000.
GRID 15 5096. 1000. 1000.
GRID 16 5824. 1000. 1000.
GRID 17 6552. 1000. 1000.
GRID 18 7280. 1000. 1000.
GRID 19 728. 0. 1000.
GRID 20 1456. 0. 1000.
GRID 21 2184. 0. 1000.
GRID 22 2912. 0. 1000.
GRID 23 3640. 0. 1000.
GRID 24 4368. 0. 1000.
GRID 25 5096. 0. 1000.
GRID 26 5824. 0. 1000.
GRID 27 6552. 0. 1000.
GRID 28 7280. 0. 1000.
GRID 29 728. 1000. 0.
GRID 30 1456. 1000. 0.
GRID 31 2184. 1000. 0.
GRID 32 2912. 1000. 0.
GRID 33 3640. 1000. 0.
GRID 34 4368. 1000. 0.
GRID 35 5096. 1000. 0.
GRID 36 5824. 1000. 0.
GRID 37 6552. 1000. 0.
GRID 38 7280. 1000. 0.
GRID 39 728. 0. 0.
GRID 40 1456. 0. 0.
GRID 41 2184. 0. 0.
GRID 42 2912. 0. 0.
GRID 43 3640. 0. 0.
GRID 44 4368. 0. 0.
GRID 45 5096. 0. 0.
GRID 46 5824. 0. 0.
GRID 47 6552. 0. 0.
GRID 48 7280. 0. 0.
GRID 49 8000. 500. 500.
GRID 50 0. 500. 500.
GRID 51 4732. 500. 500.
GRID 52 0. 500. 1000.
GRID 53 364. 1000. 1000.
GRID 54 1092. 1000. 1000.
GRID 55 1820. 1000. 1000.
GRID 56 2548. 1000. 1000.
GRID 57 3276. 1000. 1000.
GRID 58 4004. 1000. 1000.
GRID 59 4732. 1000. 1000.
GRID 60 5460. 1000. 1000.
GRID 61 6188. 1000. 1000.
GRID 62 6916. 1000. 1000.
GRID 63 7640. 1000. 1000.
GRID 64 8000. 500. 1000.
GRID 65 364. 0. 1000.
GRID 66 1092. 0. 1000.
GRID 67 1820. 0. 1000.
GRID 68 2548. 0. 1000.
GRID 69 3276. 0. 1000.
GRID 70 4004. 0. 1000.
GRID 71 4732. 0. 1000.
GRID 72 5460. 0. 1000.
GRID 73 6188. 0. 1000.
GRID 74 6916. 0. 1000.
GRID 75 7640. 0. 1000.
GRID 76 0. 500. 0.
GRID 77 364. 1000. 0.
GRID 78 1092. 1000. 0.
GRID 79 1820. 1000. 0.
GRID 80 2548. 1000. 0.
GRID 81 3276. 1000. 0.
GRID 82 4004. 1000. 0.
GRID 83 4732. 1000. 0.
GRID 84 5460. 1000. 0.
GRID 85 6188. 1000. 0.
GRID 86 6916. 1000. 0.
GRID 87 7640. 1000. 0.
GRID 88 8000. 500. 0.
GRID 89 364. 0. 0.
GRID 90 1092. 0. 0.
GRID 91 1820. 0. 0.
GRID 92 2548. 0. 0.
GRID 93 3276. 0. 0.
GRID 94 4004. 0. 0.
GRID 95 4732. 0. 0.
GRID 96 5460. 0. 0.
GRID 97 6188. 0. 0.
GRID 98 6916. 0. 0.
GRID 99 7640. 0. 0.
GRID 100 8000. 1000. 500.
GRID 101 0. 1000. 500.
GRID 102 8000. 0. 500.
GRID 103 0. 0. 500.
GRID 104 364. 500. 1000.
GRID 105 728. 500. 1000.
GRID 106 1092. 500. 1000.
GRID 107 1456. 500. 1000.
GRID 108 1820. 500. 1000.
GRID 109 2184. 500. 1000.
GRID 110 2548. 500. 1000.
GRID 111 3276. 500. 1000.
GRID 112 3640. 500. 1000.
GRID 113 4004. 500. 1000.
GRID 114 4368. 500. 1000.
GRID 115 4732. 500. 1000.
GRID 116 5096. 500. 1000.
GRID 117 5460. 500. 1000.
GRID 118 5824. 500. 1000.
GRID 119 6188. 500. 1000.
GRID 120 6552. 500. 1000.
GRID 121 6916. 500. 1000.
GRID 122 7640. 500. 1000.
GRID 123 2912. 500. 1000.
GRID 124 7280. 500. 1000.
GRID 125 364. 500. 0.
GRID 126 1092. 500. 0.
GRID 127 728. 500. 0.
GRID 128 1820. 500. 0.
GRID 129 1456. 500. 0.
GRID 130 2548. 500. 0.
GRID 131 2184. 500. 0.
GRID 132 3640. 500. 0.
GRID 133 3276. 500. 0.
GRID 134 4004. 500. 0.
GRID 135 5096. 500. 0.
GRID 136 4732. 500. 0.
GRID 137 5460. 500. 0.
GRID 138 6188. 500. 0.
GRID 139 5824. 500. 0.
GRID 140 6916. 500. 0.
GRID 141 6552. 500. 0.
GRID 142 7640. 500. 0.
GRID 143 2912. 500. 0.
GRID 144 4368. 500. 0.
GRID 145 7280. 500. 0.
GRID 146 364. 1000. 500.
GRID 147 728. 1000. 500.
GRID 148 1092. 1000. 500.
GRID 149 1456. 1000. 500.
GRID 150 1820. 1000. 500.
GRID 151 2184. 1000. 500.
GRID 152 2548. 1000. 500.
GRID 153 3276. 1000. 500.
GRID 154 3640. 1000. 500.
GRID 155 4004. 1000. 500.
GRID 156 4368. 1000. 500.
GRID 157 4732. 1000. 500.
GRID 158 5096. 1000. 500.
GRID 159 5460. 1000. 500.
GRID 160 5824. 1000. 500.
GRID 161 6188. 1000. 500.
GRID 162 6552. 1000. 500.
GRID 163 6916. 1000. 500.
GRID 164 7640. 1000. 500.
GRID 165 2912. 1000. 500.
GRID 166 7280. 1000. 500.
GRID 167 364. 0. 500.
GRID 168 1092. 0. 500.
GRID 169 728. 0. 500.
GRID 170 1820. 0. 500.
GRID 171 1456. 0. 500.
GRID 172 2548. 0. 500.
GRID 173 2184. 0. 500.
GRID 174 3640. 0. 500.
GRID 175 3276. 0. 500.
GRID 176 4004. 0. 500.
GRID 177 5096. 0. 500.
GRID 178 4732. 0. 500.
GRID 179 5460. 0. 500.
GRID 180 6188. 0. 500.
GRID 181 5824. 0. 500.
GRID 182 6916. 0. 500.
GRID 183 6552. 0. 500.
GRID 184 7640. 0. 500.
GRID 185 2912. 0. 500.
GRID 186 4368. 0. 500.
GRID 187 7280. 0. 500.
GRID 188 8000. 250. 250.
GRID 189 8000. 250. 750.
GRID 190 8000. 750. 750.
GRID 191 8000. 750. 250.
GRID 192 0. 250. 750.
GRID 193 0. 250. 250.
GRID 194 0. 750. 250.
GRID 195 0. 750. 750.
GRID 196 1456. 500. 500.
GRID 197 6552. 500. 500.
GRID 198 6916. 500. 500.
GRID 199 2184. 500. 500.
GRID 200 2548. 500. 500.
GRID 201 2912. 500. 500.
GRID 202 1820. 500. 500.
GRID 203 7640. 750. 250.
GRID 204 7640. 750. 750.
GRID 205 7280. 500. 500.
GRID 206 7640. 250. 250.
GRID 207 5460. 500. 500.
GRID 208 5096. 500. 500.
GRID 209 6188. 500. 500.
GRID 210 5824. 500. 500.
GRID 211 364. 750. 250.
GRID 212 364. 750. 750.
GRID 213 364. 250. 250.
GRID 214 1092. 500. 500.
GRID 215 728. 500. 500.
GRID 216 364. 250. 750.
GRID 217 4550. 250. 250.
GRID 218 4550. 750. 250.
GRID 219 4550. 750. 750.
GRID 220 4368. 500. 500.
GRID 221 4550. 250. 750.
GRID 222 4914. 250. 750.
GRID 223 4914. 750. 750.
GRID 224 3276. 500. 500.
GRID 225 3640. 500. 500.
GRID 226 4004. 500. 500.
GRID 227 4914. 750. 250.
GRID 228 4914. 250. 250.
$ELEMENTS
CTETRA 149 1 40 10 19 20 196 106
168 171 107 66
CTETRA 150 1 10 30 31 40 149 79
150 196 129 128
CTETRA 151 1 38 18 17 47 166 62
163 140 198 197
CTETRA 152 1 32 11 41 12 152 199
130 165 56 200
CTETRA 153 1 12 41 32 42 200 130
165 201 92 143
CTETRA 154 1 42 12 21 22 201 110
172 185 123 68
CTETRA 155 1 10 40 31 11 196 128
150 55 202 151
CTETRA 156 1 11 41 12 21 199 200
56 109 173 110
CTETRA 157 1 20 11 41 40 108 199
170 171 202 91
CTETRA 158 1 20 21 41 11 67 173
170 108 109 199
CTETRA 159 1 32 41 11 31 130 199
152 80 131 151
CTETRA 160 1 11 20 10 40 108 107
55 202 171 196
CTETRA 161 1 38 47 17 37 140 197
163 86 141 162
CTETRA 162 1 38 49 18 48 203 204
166 145 206 205
CTETRA 163 1 46 45 15 36 96 208
207 139 137 159
CTETRA 164 1 47 46 16 37 97 210
209 141 138 161
CTETRA 165 1 18 2 4 49 63 64
122 204 190 189
CTETRA 166 1 18 27 17 47 121 120
62 198 183 197
CTETRA 167 1 38 48 18 47 145 205
166 140 98 198
CTETRA 168 1 12 42 23 22 201 175
111 123 185 69
CTETRA 169 1 26 27 47 17 73 183
180 119 120 197
CTETRA 170 1 50 9 29 6 212 147
211 195 53 146
CTETRA 171 1 27 18 48 47 121 205
182 183 198 98
CTETRA 172 1 8 39 7 50 167 89
103 192 213 193
CTETRA 173 1 40 39 9 30 90 215
214 129 126 148
CTETRA 174 1 42 41 21 12 92 173
172 201 200 110
CTETRA 175 1 50 29 9 39 211 147
212 213 127 215
CTETRA 176 1 29 50 7 39 211 193
125 127 213 89
CTETRA 177 1 31 41 11 40 131 199
151 128 91 202
CTETRA 178 1 47 37 16 17 141 161
209 197 162 61
CTETRA 179 1 40 30 9 10 129 148
214 196 149 54
CTETRA 180 1 2 49 38 1 190 203
164 100 191 87
CTETRA 181 1 2 18 38 49 63 166
164 190 204 203
CTETRA 182 1 48 38 49 3 145 203
206 99 142 188
CTETRA 183 1 38 1 49 3 87 191
203 142 88 188
CTETRA 184 1 49 3 4 48 188 102
189 206 99 184
CTETRA 185 1 28 4 48 18 75 184
187 124 122 205
CTETRA 186 1 49 4 18 48 189 122
204 206 184 205
CTETRA 187 1 7 5 50 29 76 194
193 125 77 211
CTETRA 188 1 50 5 6 29 194 101
195 211 77 146
CTETRA 189 1 50 19 9 6 216 105
212 195 104 53
CTETRA 190 1 50 9 19 39 212 105
216 213 215 169
CTETRA 191 1 50 8 19 6 192 65
216 195 52 104
CTETRA 192 1 40 10 9 19 196 54
214 168 106 105
CTETRA 193 1 51 34 44 14 218 144
217 219 156 220
CTETRA 194 1 51 44 24 14 217 186
221 219 220 114
CTETRA 195 1 25 24 15 51 71 115
116 222 221 223
CTETRA 196 1 43 32 12 13 133 165
224 225 153 57
CTETRA 197 1 43 23 12 42 174 111
224 93 175 201
CTETRA 198 1 43 13 12 23 225 57
224 174 112 111
CTETRA 199 1 43 42 12 32 93 201
224 133 143 165
CTETRA 200 1 34 14 13 44 156 58
155 144 220 226
CTETRA 201 1 14 15 24 51 59 115
114 219 223 221
CTETRA 202 1 23 14 24 44 113 114
70 176 220 186
CTETRA 203 1 33 13 32 43 154 153
81 132 225 133
CTETRA 204 1 34 13 33 43 155 154
82 134 225 132
CTETRA 205 1 35 15 14 51 158 59
157 227 223 219
CTETRA 206 1 25 15 45 51 116 208
177 222 223 228
CTETRA 207 1 44 13 43 23 226 225
94 176 112 174
CTETRA 208 1 35 14 34 51 157 156
83 227 219 218
CTETRA 209 1 46 36 15 16 139 159
207 210 160 60
CTETRA 210 1 36 15 35 45 159 158
84 137 208 135
CTETRA 211 1 37 16 36 46 161 160
85 138 210 139
CTETRA 212 1 25 16 26 46 117 118
72 179 210 181
CTETRA 213 1 47 17 16 26 197 61
209 180 119 118
CTETRA 214 1 47 26 16 46 180 118
209 97 181 210
CTETRA 215 1 27 18 28 48 121 124
74 182 205 187
CTETRA 216 1 35 51 34 45 227 218
83 135 228 136
CTETRA 217 1 13 44 14 23 226 220
58 112 176 113
CTETRA 218 1 44 24 25 51 186 71
178 217 221 222
CTETRA 219 1 44 25 45 51 178 177
95 217 222 228
CTETRA 220 1 46 16 15 25 210 60
207 179 117 116
CTETRA 221 1 46 25 15 45 179 116
207 96 177 208
CTETRA 222 1 50 19 8 39 216 65
192 213 169 167
CTETRA 223 1 35 15 51 45 158 223
227 135 208 228
CTETRA 224 1 34 43 44 13 134 94
144 155 225 226
CTETRA 225 1 51 45 44 34 228 95
217 218 136 144
CTETRA 226 1 9 39 29 30 215 127
147 148 126 78
CTETRA 227 1 40 19 9 39 168 105
214 90 169 215
$PROPERTIES
PSOLID 1 1
$MATERIALS
MAT1 1 210000. .3
$LOADS
LOAD 1 1. 1. 2
FORCE 2 4 2250000. 0. 0. -1.
FORCE 2 3 2250000. 0. 0. -1.
FORCE 2 2 2250000. 0. 0. -1.
FORCE 2 1 2250000. 0. 0. -1.
$SPCs
SPCADD 1 2
SPC1 2 123456 5 6 7 8 50 52
76 101 103 192 193 194 195
ENDDATA

100
src/Mod/Fem/femtest/data/mystran/mystran_plate.bdf Normal file
View File

@@ -0,0 +1,100 @@
$pyNastran: version=msc
$pyNastran: punch=False
$pyNastran: encoding=utf-8
$pyNastran: nnodes=36
$pyNastran: nelements=25
$EXECUTIVE CONTROL DECK
SOL 101
CEND
$CASE CONTROL DECK
ECHO = NONE
TITLE = pyNastran for generating solverinput for for Mystran
SUBCASE 1
DISPLACEMENT(SORT1,REAL) = ALL
LOAD = 1
SPC = 1
SPCFORCES(SORT1,REAL) = ALL
STRESS(SORT1,REAL,VONMISES,BILIN) = ALL
SUBTITLE = Default
BEGIN BULK
$PARAMS
PARAM POST -1
$NODES
GRID 1 0. 0. 0.
GRID 2 2. 0. 0.
GRID 3 4. 0. 0.
GRID 4 6. 0. 0.
GRID 5 8. 0. 0.
GRID 6 10. 0. 0.
GRID 7 0. 2. 0.
GRID 8 2. 2. 0.
GRID 9 4. 2. 0.
GRID 10 6. 2. 0.
GRID 11 8. 2. 0.
GRID 12 10. 2. 0.
GRID 13 0. 4. 0.
GRID 14 2. 4. 0.
GRID 15 4. 4. 0.
GRID 16 6. 4. 0.
GRID 17 8. 4. 0.
GRID 18 10. 4. 0.
GRID 19 0. 6. 0.
GRID 20 2. 6. 0.
GRID 21 4. 6. 0.
GRID 22 6. 6. 0.
GRID 23 8. 6. 0.
GRID 24 10. 6. 0.
GRID 25 0. 8. 0.
GRID 26 2. 8. 0.
GRID 27 4. 8. 0.
GRID 28 6. 8. 0.
GRID 29 8. 8. 0.
GRID 30 10. 8. 0.
GRID 31 0. 10. 0.
GRID 32 2. 10. 0.
GRID 33 4. 10. 0.
GRID 34 6. 10. 0.
GRID 35 8. 10. 0.
GRID 36 10. 10. 0.
$ELEMENTS
CQUAD4 1 1 1 2 8 7
CQUAD4 2 1 2 3 9 8
CQUAD4 3 1 3 4 10 9
CQUAD4 4 1 4 5 11 10
CQUAD4 5 1 5 6 12 11
CQUAD4 6 1 7 8 14 13
CQUAD4 7 1 8 9 15 14
CQUAD4 8 1 9 10 16 15
CQUAD4 9 1 10 11 17 16
CQUAD4 10 1 11 12 18 17
CQUAD4 11 1 13 14 20 19
CQUAD4 12 1 14 15 21 20
CQUAD4 13 1 15 16 22 21
CQUAD4 14 1 16 17 23 22
CQUAD4 15 1 17 18 24 23
CQUAD4 16 1 19 20 26 25
CQUAD4 17 1 20 21 27 26
CQUAD4 18 1 21 22 28 27
CQUAD4 19 1 22 23 29 28
CQUAD4 20 1 23 24 30 29
CQUAD4 21 1 25 26 32 31
CQUAD4 22 1 26 27 33 32
CQUAD4 23 1 27 28 34 33
CQUAD4 24 1 28 29 35 34
CQUAD4 25 1 29 30 36 35
$PROPERTIES
PSHELL 1 1 .3 1 1
$MATERIALS
MAT1 1 210000. .3
$LOADS
LOAD 1 1. 1. 2
FORCE 2 6 100. 1. 0. 0.
FORCE 2 12 100. 1. 0. 0.
FORCE 2 18 100. 1. 0. 0.
FORCE 2 24 100. 1. 0. 0.
FORCE 2 30 100. 1. 0. 0.
FORCE 2 36 100. 1. 0. 0.
$SPCs
SPCADD 1 2
SPC1 2 123456 1 7 13 19 25 31
ENDDATA

1
src/Mod/Fem/femtest/data/z88/ccx_cantilever_ele_tria6/51.txt Normal file
View File

@@ -0,0 +1 @@
210000.0 0.300

49
src/Mod/Fem/femtest/data/z88/ccx_cantilever_ele_tria6/z88.dyn Normal file
View File

@@ -0,0 +1,49 @@
DYNAMIC START
---------------------------------------------------------------------------
Z88 new version 14OS Z88 neue Version 14OS
---------------------------------------------------------------------------
---------------------------------------------------------------------------
LANGUAGE SPRACHE
---------------------------------------------------------------------------
GERMAN
---------------------------------------------------------------------------
Entries for mesh generator Z88N Daten fuer Netzgenerator
---------------------------------------------------------------------------
NET START
MAXSE 40000
MAXESS 800
MAXKSS 4000
MAXAN 15
NET END
---------------------------------------------------------------------------
Common entries for all modules gemeinsame Daten fuer alle Module
---------------------------------------------------------------------------
COMMON START
MAXGS 50000000
MAXKOI 1200000
MAXK 60000
MAXE 300000
MAXNFG 200000
MAXMAT 32
MAXPEL 32
MAXJNT 32
MAXPR 10000
MAXRBD 15000
MAXIEZ 6000000
MAXGP 2000000
COMMON END
---------------------------------------------------------------------------
Entries for Cuthill-McKee Z88H Daten fuer Cuthill- McKee Programm
---------------------------------------------------------------------------
CUTKEE START
MAXGRA 200
MAXNDL 1000
CUTKEE END
DYNAMIC END

2
src/Mod/Fem/femtest/data/z88/ccx_cantilever_ele_tria6/z88elp.txt Normal file
View File

@@ -0,0 +1,2 @@
1
1 70 1000.0 0 0 0 0 0 0

318
src/Mod/Fem/femtest/data/z88/ccx_cantilever_ele_tria6/z88i1.txt Normal file
View File

@@ -0,0 +1,318 @@
3 177 70 1062 0 written by FreeCAD
1 6 0.000000 500.000000 0.000000
2 6 0.000000 500.000000 1000.000000
3 6 8000.000000 500.000000 0.000000
4 6 8000.000000 500.000000 1000.000000
5 6 0.000000 500.000000 500.000000
6 6 0.000000 500.000000 250.000000
7 6 0.000000 500.000000 750.000000
8 6 500.000000 500.000000 0.000000
9 6 1000.000000 500.000000 0.000000
10 6 1500.000000 500.000000 0.000000
11 6 2000.000000 500.000000 0.000000
12 6 2500.000000 500.000000 0.000000
13 6 3000.000000 500.000000 0.000000
14 6 3500.000000 500.000000 0.000000
15 6 4000.000000 500.000000 0.000000
16 6 4500.000000 500.000000 0.000000
17 6 5000.000000 500.000000 0.000000
18 6 5500.000000 500.000000 0.000000
19 6 6000.000000 500.000000 0.000000
20 6 6500.000000 500.000000 0.000000
21 6 7000.000000 500.000000 0.000000
22 6 7500.000000 500.000000 0.000000
23 6 250.000000 500.000000 0.000000
24 6 750.000000 500.000000 0.000000
25 6 1250.000000 500.000000 0.000000
26 6 1750.000000 500.000000 0.000000
27 6 2250.000000 500.000000 0.000000
28 6 2750.000000 500.000000 0.000000
29 6 3250.000000 500.000000 0.000000
30 6 3750.000000 500.000000 0.000000
31 6 4250.000000 500.000000 0.000000
32 6 4750.000000 500.000000 0.000000
33 6 5250.000000 500.000000 0.000000
34 6 5750.000000 500.000000 0.000000
35 6 6250.000000 500.000000 0.000000
36 6 6750.000000 500.000000 0.000000
37 6 7250.000000 500.000000 0.000000
38 6 7750.000000 500.000000 0.000000
39 6 8000.000000 500.000000 500.000000
40 6 8000.000000 500.000000 250.000000
41 6 8000.000000 500.000000 750.000000
42 6 500.000000 500.000000 1000.000000
43 6 1000.000000 500.000000 1000.000000
44 6 1500.000000 500.000000 1000.000000
45 6 2000.000000 500.000000 1000.000000
46 6 2500.000000 500.000000 1000.000000
47 6 3000.000000 500.000000 1000.000000
48 6 3500.000000 500.000000 1000.000000
49 6 4000.000000 500.000000 1000.000000
50 6 4500.000000 500.000000 1000.000000
51 6 5000.000000 500.000000 1000.000000
52 6 5500.000000 500.000000 1000.000000
53 6 6000.000000 500.000000 1000.000000
54 6 6500.000000 500.000000 1000.000000
55 6 7000.000000 500.000000 1000.000000
56 6 7500.000000 500.000000 1000.000000
57 6 250.000000 500.000000 1000.000000
58 6 750.000000 500.000000 1000.000000
59 6 1250.000000 500.000000 1000.000000
60 6 1750.000000 500.000000 1000.000000
61 6 2250.000000 500.000000 1000.000000
62 6 2750.000000 500.000000 1000.000000
63 6 3250.000000 500.000000 1000.000000
64 6 3750.000000 500.000000 1000.000000
65 6 4250.000000 500.000000 1000.000000
66 6 4750.000000 500.000000 1000.000000
67 6 5250.000000 500.000000 1000.000000
68 6 5750.000000 500.000000 1000.000000
69 6 6250.000000 500.000000 1000.000000
70 6 6750.000000 500.000000 1000.000000
71 6 7250.000000 500.000000 1000.000000
72 6 7750.000000 500.000000 1000.000000
73 6 714.923469 500.000000 484.056122
74 6 1750.000000 500.000000 500.000000
75 6 2750.000000 500.000000 500.000000
76 6 3750.000000 500.000000 500.000000
77 6 4750.000000 500.000000 500.000000
78 6 5750.000000 500.000000 500.000000
79 6 6756.944444 500.000000 489.583333
80 6 7301.463294 500.000000 489.707341
81 6 1244.153912 500.000000 497.342687
82 6 2250.000000 500.000000 500.000000
83 6 3250.000000 500.000000 500.000000
84 6 4250.000000 500.000000 500.000000
85 6 5250.000000 500.000000 500.000000
86 6 6251.157407 500.000000 498.263889
87 6 355.277669 500.000000 632.304567
88 6 7660.398197 500.000000 648.596526
89 6 7699.764778 500.000000 310.456177
90 6 314.040228 500.000000 323.272138
91 6 1122.076956 500.000000 748.671344
92 6 1372.076956 500.000000 748.671344
93 6 2125.000000 500.000000 750.000000
94 6 2375.000000 500.000000 750.000000
95 6 3125.000000 500.000000 750.000000
96 6 3375.000000 500.000000 750.000000
97 6 4125.000000 500.000000 750.000000
98 6 4375.000000 500.000000 750.000000
99 6 5125.000000 500.000000 750.000000
100 6 5375.000000 500.000000 750.000000
101 6 6125.578704 500.000000 749.131944
102 6 6375.578704 500.000000 749.131944
103 6 857.461735 500.000000 742.028061
104 6 979.538690 500.000000 490.699405
105 6 1497.076956 500.000000 498.671344
106 6 1625.000000 500.000000 750.000000
107 6 1875.000000 500.000000 750.000000
108 6 2000.000000 500.000000 500.000000
109 6 2500.000000 500.000000 500.000000
110 6 2625.000000 500.000000 750.000000
111 6 2875.000000 500.000000 750.000000
112 6 3000.000000 500.000000 500.000000
113 6 3500.000000 500.000000 500.000000
114 6 3625.000000 500.000000 750.000000
115 6 3875.000000 500.000000 750.000000
116 6 4000.000000 500.000000 500.000000
117 6 4500.000000 500.000000 500.000000
118 6 4625.000000 500.000000 750.000000
119 6 4875.000000 500.000000 750.000000
120 6 5000.000000 500.000000 500.000000
121 6 5500.000000 500.000000 500.000000
122 6 5625.000000 500.000000 750.000000
123 6 5875.000000 500.000000 750.000000
124 6 6000.578704 500.000000 499.131944
125 6 6504.050926 500.000000 493.923611
126 6 6628.472222 500.000000 744.791667
127 6 7150.731647 500.000000 244.853671
128 6 7029.203869 500.000000 489.645337
129 6 6878.472222 500.000000 244.791667
130 6 6878.472222 500.000000 744.791667
131 6 7150.731647 500.000000 744.853671
132 6 7400.731647 500.000000 244.853671
133 6 7400.731647 500.000000 744.853671
134 6 857.461735 500.000000 242.028061
135 6 607.461735 500.000000 242.028061
136 6 1875.000000 500.000000 250.000000
137 6 1625.000000 500.000000 250.000000
138 6 2875.000000 500.000000 250.000000
139 6 2625.000000 500.000000 250.000000
140 6 3875.000000 500.000000 250.000000
141 6 3625.000000 500.000000 250.000000
142 6 4875.000000 500.000000 250.000000
143 6 4625.000000 500.000000 250.000000
144 6 5875.000000 500.000000 250.000000
145 6 5625.000000 500.000000 250.000000
146 6 6628.472222 500.000000 244.791667
147 6 607.461735 500.000000 742.028061
148 6 1122.076956 500.000000 248.671344
149 6 1372.076956 500.000000 248.671344
150 6 2125.000000 500.000000 250.000000
151 6 2375.000000 500.000000 250.000000
152 6 3125.000000 500.000000 250.000000
153 6 3375.000000 500.000000 250.000000
154 6 4125.000000 500.000000 250.000000
155 6 4375.000000 500.000000 250.000000
156 6 5125.000000 500.000000 250.000000
157 6 5375.000000 500.000000 250.000000
158 6 6125.578704 500.000000 249.131944
159 6 6375.578704 500.000000 249.131944
160 6 514.481849 500.000000 403.664130
161 6 407.020114 500.000000 161.636069
162 6 177.638835 500.000000 566.152283
163 6 177.638835 500.000000 816.152283
164 6 427.638835 500.000000 816.152283
165 6 535.100569 500.000000 558.180345
166 6 7599.882389 500.000000 155.228088
167 6 7500.614036 500.000000 400.081759
168 6 7830.199099 500.000000 824.298263
169 6 7830.199099 500.000000 574.298263
170 6 7580.199099 500.000000 824.298263
171 6 7480.930745 500.000000 569.151933
172 6 334.658948 500.000000 477.788352
173 6 7849.882389 500.000000 405.228088
174 6 7849.882389 500.000000 155.228088
175 6 157.020114 500.000000 161.636069
176 6 157.020114 500.000000 411.636069
177 6 7680.081488 500.000000 479.526351
37 24
43 81 44 91 92 59
38 24
45 82 46 93 94 61
39 24
47 83 48 95 96 63
40 24
49 84 50 97 98 65
41 24
51 85 52 99 100 67
42 24
53 86 54 101 102 69
43 24
43 73 81 103 104 91
44 24
44 81 74 92 105 106
45 24
45 74 82 107 108 93
46 24
46 82 75 94 109 110
47 24
47 75 83 111 112 95
48 24
48 83 76 96 113 114
49 24
49 76 84 115 116 97
50 24
50 84 77 98 117 118
51 24
51 77 85 119 120 99
52 24
52 85 78 100 121 122
53 24
53 78 86 123 124 101
54 24
54 86 79 102 125 126
55 24
21 80 79 127 128 129
56 24
55 79 80 130 128 131
57 24
21 22 80 37 132 127
58 24
55 80 56 131 133 71
59 24
8 9 73 24 134 135
60 24
10 11 74 26 136 137
61 24
12 13 75 28 138 139
62 24
14 15 76 30 140 141
63 24
16 17 77 32 142 143
64 24
18 19 78 34 144 145
65 24
20 21 79 36 129 146
66 24
42 73 43 147 103 58
67 24
44 74 45 106 107 60
68 24
46 75 47 110 111 62
69 24
48 76 49 114 115 64
70 24
50 77 51 118 119 66
71 24
52 78 53 122 123 68
72 24
54 79 55 126 130 70
73 24
9 81 73 148 104 134
74 24
10 74 81 137 105 149
75 24
11 82 74 150 108 136
76 24
12 75 82 139 109 151
77 24
13 83 75 152 112 138
78 24
14 76 83 141 113 153
79 24
15 84 76 154 116 140
80 24
16 77 84 143 117 155
81 24
17 85 77 156 120 142
82 24
18 78 85 145 121 157
83 24
19 86 78 158 124 144
84 24
20 79 86 146 125 159
85 24
9 10 81 25 149 148
86 24
11 12 82 27 151 150
87 24
13 14 83 29 153 152
88 24
15 16 84 31 155 154
89 24
17 18 85 33 157 156
90 24
19 20 86 35 159 158
91 24
8 73 90 135 160 161
92 24
2 5 87 7 162 163
93 24
2 87 42 163 164 57
94 24
42 87 73 164 165 147
95 24
22 89 80 166 167 132
96 24
4 88 39 168 169 41
97 24
4 56 88 72 170 168
98 24
56 80 88 133 171 170
99 24
73 87 90 165 172 160
100 24
3 39 89 40 173 174
101 24
3 89 22 174 166 38
102 24
1 90 5 175 176 6
103 24
1 8 90 23 161 175
104 24
80 89 88 167 177 171
105 24
5 90 87 176 172 162
106 24
39 88 89 169 177 173

26
src/Mod/Fem/femtest/data/z88/ccx_cantilever_ele_tria6/z88i2.txt Normal file
View File

@@ -0,0 +1,26 @@
25
1 1 2 0
1 2 2 0
1 3 2 0
2 1 2 0
2 2 2 0
2 3 2 0
3 2 1 -1.6653345369377348e-10
3 3 1 -750000.0
4 2 1 -1.6653345369377348e-10
4 3 1 -750000.0
5 1 2 0
5 2 2 0
5 3 2 0
6 1 2 0
6 2 2 0
6 3 2 0
7 1 2 0
7 2 2 0
7 3 2 0
39 2 1 -3.3306690738754696e-10
39 3 1 -1500000.0
40 2 1 -6.661338147750939e-10
40 3 1 -3000000.0
41 2 1 -6.661338147750939e-10
41 3 1 -3000000.0

1
src/Mod/Fem/femtest/data/z88/ccx_cantilever_ele_tria6/z88i5.txt Normal file
View File

@@ -0,0 +1 @@
0

2
src/Mod/Fem/femtest/data/z88/ccx_cantilever_ele_tria6/z88int.txt Normal file
View File

@@ -0,0 +1,2 @@
1
1 70 7 7

36
src/Mod/Fem/femtest/data/z88/ccx_cantilever_ele_tria6/z88man.txt Normal file
View File

@@ -0,0 +1,36 @@
DYNAMIC START
---------------------------------------------------------------------------
Z88V14OS
---------------------------------------------------------------------------
---------------------------------------------------------------------------
GLOBAL
---------------------------------------------------------------------------
GLOBAL START
IBFLAG 0
IPFLAG 0
IHFLAG 1
GLOBAL END
---------------------------------------------------------------------------
LINEAR SOLVER
---------------------------------------------------------------------------
SOLVER START
MAXIT 10000
EPS 1e-007
RALPHA 0.0001
ROMEGA 1.1
SOLVER END
---------------------------------------------------------------------------
STRESS
---------------------------------------------------------------------------
STRESS START
KDFLAG 0
ISFLAG 1
STRESS END
DYNAMIC END

2
src/Mod/Fem/femtest/data/z88/ccx_cantilever_ele_tria6/z88mat.txt Normal file
View File

@@ -0,0 +1,2 @@
1
1 70 51.txt

6
src/Mod/Fem/femtest/failing_tests.md Normal file
View File

@@ -0,0 +1,6 @@
+ the following tests fail on ci, TODO somehow fix, even if they pass on my local machine:
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_constraint_centrif
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_constraint_contact_solid_solid
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_z88.TestSolverZ88.test_ccx_cantilever_ele_tria6

114
src/Mod/Fem/femtest/test_commands.sh
View File

@@ -4,6 +4,12 @@ from femtest.app.support_utils import get_fem_test_defs
get_fem_test_defs()
# all FEM App tests
make -j 4 && ./bin/FreeCAD --run-test 'TestFemApp'
make -j 4 && ./bin/FreeCADCmd --run-test 'TestFemApp'
'''
# modules
@@ -17,6 +23,7 @@ make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_open
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_result
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_elmer
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_mystran
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_z88
@@ -35,6 +42,7 @@ make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_open.TestObjectOpen
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_result.TestResult
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_elmer.TestSolverElmer
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_mystran.TestSolverMystran
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_z88.TestSolverZ88
@@ -78,10 +86,16 @@ make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_result.TestResult.test_disp_ab
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_box_frequency
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_box_static
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_ccx_buckling_flexuralbuckling
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_ccxcantilever_faceload
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_ccxcantilever_hexa20
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_ccxcantilever_nodeload
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_ccxcantilever_prescribeddisplacement
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_ccx_cantilever_faceload
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_ccx_cantilever_ele_hexa20
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_ccx_cantilever_ele_quad4
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_ccx_cantilever_ele_quad8
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_ccx_cantilever_ele_seg2
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_ccx_cantilever_ele_seg3
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_ccx_cantilever_ele_tria3
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_ccx_cantilever_ele_tria6
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_ccx_cantilever_nodeload
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_ccx_cantilever_prescribeddisplacement
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_constraint_centrif
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_constraint_contact_shell_shell
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_calculix.TestSolverCalculix.test_constraint_contact_solid_solid
@@ -105,9 +119,16 @@ make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_elmer.TestSolverElmer.t
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_elmer.TestSolverElmer.test_ccxcantilever_faceload_1_si
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_elmer.TestSolverElmer.test_ccxcantilever_nodeload_0_mm
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_elmer.TestSolverElmer.test_ccxcantilever_prescribeddisplacement_0_mm
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_z88.TestSolverZ88.test_ccxcantilever_faceload
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_z88.TestSolverZ88.test_ccxcantilever_hexa20
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_z88.TestSolverZ88.test_ccxcantilever_nodeload
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_mystran.TestSolverMystran.test_ccx_cantilever_ele_quad4
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_mystran.TestSolverMystran.test_ccx_cantilever_ele_seg2
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_mystran.TestSolverMystran.test_ccx_cantilever_ele_tria3
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_mystran.TestSolverMystran.test_ccx_cantilever_faceload
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_mystran.TestSolverMystran.test_ccx_cantilever_nodeload
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_mystran.TestSolverMystran.test_mystran_plate
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_z88.TestSolverZ88.test_ccx_cantilever_ele_hexa20
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_z88.TestSolverZ88.test_ccx_cantilever_ele_tria6
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_z88.TestSolverZ88.test_ccx_cantilever_faceload
make -j 4 && ./bin/FreeCADCmd -t femtest.app.test_solver_z88.TestSolverZ88.test_ccx_cantilever_nodeload
# methods in FreeCAD
@@ -309,22 +330,52 @@ unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_calculix.TestSolverCalculix.test_ccxcantilever_faceload'
'femtest.app.test_solver_calculix.TestSolverCalculix.test_ccx_cantilever_faceload'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_calculix.TestSolverCalculix.test_ccxcantilever_hexa20'
'femtest.app.test_solver_calculix.TestSolverCalculix.test_ccx_cantilever_ele_hexa20'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_calculix.TestSolverCalculix.test_ccxcantilever_nodeload'
'femtest.app.test_solver_calculix.TestSolverCalculix.test_ccx_cantilever_ele_quad4'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_calculix.TestSolverCalculix.test_ccxcantilever_prescribeddisplacement'
'femtest.app.test_solver_calculix.TestSolverCalculix.test_ccx_cantilever_ele_quad8'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_calculix.TestSolverCalculix.test_ccx_cantilever_ele_seg2'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_calculix.TestSolverCalculix.test_ccx_cantilever_ele_seg3'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_calculix.TestSolverCalculix.test_ccx_cantilever_ele_tria3'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_calculix.TestSolverCalculix.test_ccx_cantilever_ele_tria6'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_calculix.TestSolverCalculix.test_ccx_cantilever_nodeload'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_calculix.TestSolverCalculix.test_ccx_cantilever_prescribeddisplacement'
))
import unittest
@@ -444,15 +495,50 @@ unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_z88.TestSolverZ88.test_ccxcantilever_faceload'
'femtest.app.test_solver_mystran.TestSolverMystran.test_ccx_cantilever_ele_quad4'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_z88.TestSolverZ88.test_ccxcantilever_hexa20'
'femtest.app.test_solver_mystran.TestSolverMystran.test_ccx_cantilever_ele_seg2'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_z88.TestSolverZ88.test_ccxcantilever_nodeload'
'femtest.app.test_solver_mystran.TestSolverMystran.test_ccx_cantilever_ele_tria3'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_mystran.TestSolverMystran.test_ccx_cantilever_faceload'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_mystran.TestSolverMystran.test_ccx_cantilever_nodeload'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_mystran.TestSolverMystran.test_mystran_plate'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_z88.TestSolverZ88.test_ccx_cantilever_ele_hexa20'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_z88.TestSolverZ88.test_ccx_cantilever_ele_tria6'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_z88.TestSolverZ88.test_ccx_cantilever_faceload'
))
import unittest
unittest.TextTestRunner().run(unittest.TestLoader().loadTestsFromName(
'femtest.app.test_solver_z88.TestSolverZ88.test_ccx_cantilever_nodeload'
))

Some files were not shown because too many files have changed in this diff Show More