Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

[Problem] FEM ConstraintForce may generate values in an incorrect format causing a CalculiX error #12805

Closed
2 tasks done
FEA-eng opened this issue Mar 8, 2024 · 2 comments
Closed
2 tasks done

[Problem] FEM ConstraintForce may generate values in an incorrect format causing a CalculiX error #12805

FEA-eng opened this issue Mar 8, 2024 · 2 comments
Labels
Bug This issue or PR is related to a bug WB FEM Related to the FEM Workbench

Comments

@FEA-eng
Copy link
Contributor

FEA-eng commented Mar 8, 2024

Is there an existing issue for this?

  • I have searched the existing issues

Problem description

While helping a user (forum thread), I've noticed that FEM ConstraintForce may generate values which can't be correctly read by CalculiX. The following error appears when submitting an analysis:

 *ERROR reading *CLOAD. Card image:
        6833,1,1.7302535485473438E-06

This is caused by the fact that CalculiX uses the F20.0 Fortran format to read the values and can only read 20 characters so the following value:

1.7302535485473438E-06

is read as:

1.7302535485473438E-

which is obviously incorrect.

Thus, it would be necessary to ensure that FEM ConstraintForce doesn't generate so many decimals and follows the F20.0 format required by CalculiX.

Model (just change the extension to .FCStd and regenerate the mesh, then try running this):
Number format error.zip

@marioalexis84 Do you know how to fix it ?

Full version info

OS: Windows 10 build 19045
Word size of FreeCAD: 64-bit
Version: 0.22.0dev.36277 (Git)
Build type: Release
Branch: main
Hash: 9e1903d46112b3660bf10c6a4537d728101d560b
Python 3.10.13, Qt 5.15.8, Coin 4.0.2, Vtk 9.2.6, OCC 7.6.3
Locale: Polish/Poland (pl_PL)

Subproject(s) affected?

FEM

Anything else?

No response

Code of Conduct

  • I agree to follow this project's Code of Conduct
@maxwxyz maxwxyz added Bug This issue or PR is related to a bug WB FEM Related to the FEM Workbench labels Mar 8, 2024
UR-0 pushed a commit to UR-0/FreeCAD that referenced this issue Mar 10, 2024
FEM ConstraintForce may generate values in an incorrect format causin…
2e2154b 
…g a CalculiX error

issue FreeCAD#12805
fixed by forcing to format {:.13G}
@UR-0 UR-0 mentioned this issue Mar 10, 2024
Closed
@marioalexis84
Copy link
Member

Model (just change the extension to .FCStd and regenerate the mesh, then try running this):
Number format error.zip

For me, gmsh fails to mesh the compound shape.

@FEA-eng
Copy link
Contributor Author

FEA-eng commented Mar 13, 2024

For me, gmsh fails to mesh the compound shape.

Strange. But maybe you have different Gmsh settings. Mine are:

gmsh

Especially the Second Order Linear property has to be true.

I only get a warning (2 ill-shaped tets) but no errors and the analysis works if I replace the force with pressure.

chennes added a commit that referenced this issue Mar 13, 2024
@chennes
Merge pull request #12932 from marioalexis84/fem-ccx_force_format
5264d6c 
Fem: Format force constraint in CalculiX writer - fixes #12805
davesrocketshop pushed a commit to davesrocketshop/FreeCAD that referenced this issue Mar 16, 2024
@marioalexis84 @davesrocketshop
Fem: Format force constraint in CalculiX writer - fixes FreeCAD#12805
3c4cbac
davesrocketshop added a commit to davesrocketshop/FreeCAD that referenced this issue Mar 17, 2024
@davesrocketshop
Material: Material appearance
e6c1747 
Uses new materail system for appearance

Each feature object now has a property called ShapeMaterial that
describes its physical properties. If it has a shape, it has a
material.

The ShapeColor attribute is replaced by a ShapeAppearance attribute.
This is a material list that describes all appearance properties, not
just diffuse color. As a list in can be used for all elements of a
shape, such as edges and faces.

A new widget is provided to allow the user to select materials in a
consistent fashion. It can also launch the material editor with its
more advanced capabilities.

commit 0c77489
Author: David Carter <dcarter@davidcarter.ca>
Date:   Sun Mar 17 03:07:28 2024 -0400

    Migrate view parameters

commit 25b0740
Author: David Carter <dcarter@davidcarter.ca>
Date:   Sat Mar 16 20:18:00 2024 -0400

    Build fixes

commit 73161b2
Author: David Carter <dcarter@davidcarter.ca>
Date:   Sat Mar 16 17:57:43 2024 -0400

    Merge fixes

commit e950dea
Author: David Carter <dcarter@davidcarter.ca>
Date:   Sat Mar 16 12:23:16 2024 -0400

    Rename Material library to Materials

commit cb20474
Author: David Carter <dcarter@davidcarter.ca>
Date:   Sat Mar 16 01:13:20 2024 -0400

    Build fixes

commit 49fa5ff
Author: David Carter <dcarter@davidcarter.ca>
Date:   Tue Mar 12 21:37:45 2024 -0400

    Compatibility with perevious versions

commit 1ea547c
Author: David Carter <dcarter@davidcarter.ca>
Date:   Sat Mar 2 01:19:47 2024 -0500

    ShapeMaterial becomes part of Part::Feature

commit 721bfe8
Author: David Carter <dcarter@davidcarter.ca>
Date:   Mon Feb 26 23:29:01 2024 -0500

    Set material option

commit 8535cef
Author: David Carter <dcarter@davidcarter.ca>
Date:   Fri Feb 23 13:54:13 2024 -0500

    Use materials to set appearance

commit 8bbc04e
Author: David Carter <dcarter@davidcarter.ca>
Date:   Wed Feb 14 10:35:11 2024 -0500

    ShapeColor to ShapeMaterial

commit f86dcde
Author: David Carter <dcarter@davidcarter.ca>
Date:   Wed Feb 14 03:42:30 2024 -0500

    Pass more unit tests

commit 53b1126
Author: David Carter <dcarter@davidcarter.ca>
Date:   Tue Feb 13 21:38:33 2024 -0500

    Creaate and use ShapeMaterial

commit d107c63
Author: Max Wilfinger <6246609+maxwxyz@users.noreply.github.com>
Date:   Sat Mar 16 14:41:46 2024 +0100

    Update labels.yml for correct GH label

commit 10a6538
Author: wmayer <wmayer@users.sourceforge.net>
Date:   Fri Mar 15 15:09:13 2024 +0100

    Gui: fixes FreeCAD#12891: Random color not working for App::Part

commit 93cb44f
Author: wmayer <wmayer@users.sourceforge.net>
Date:   Fri Mar 15 15:08:50 2024 +0100

    App: add overloaded template method getExtension()

    This allows it to write more elegant code like:
    auto ext = obj->getExtension<>(GroupExtension);
    instead of
    auto ext = dynamic_cast<GroupExtension*>(obj->getExtension(GroupExtension::getExtensionClassTypeId(), true, true));

commit c19dedf
Author: wmayer <wmayer@users.sourceforge.net>
Date:   Sat Mar 16 01:45:37 2024 +0100

    TD: fixes FreeCAD#12984: Doing Ctrl-Z in TechDraw reliably causes a segfault

commit 0e123b3
Author: wmayer <wmayer@users.sourceforge.net>
Date:   Fri Mar 15 11:45:11 2024 +0100

    Drawing: Fixes FreeCAD#12953

    The deprecated Drawing module still uses some classes that have been removed with Qt6.
    These are the classes:
    * QRegExp that is replaced with QRegularExpression
    * QGLWidget that is replaced with QOpenGLWidget

commit 286aae6
Author: wmayer <wmayer@users.sourceforge.net>
Date:   Wed Mar 13 12:39:23 2024 +0100

    Mesh: fixes FreeCAD#12798: Macro not recording mesh decimation

commit dab9a95
Author: wmayer <wmayer@users.sourceforge.net>
Date:   Wed Mar 13 11:54:51 2024 +0100

    Mesh: add unit tests for mesh export

commit 2065c6e
Author: wmayer <wmayer@users.sourceforge.net>
Date:   Mon Mar 11 07:59:25 2024 +0100

    App: Fix segmentation fault while dragging object

    For more details see: https://forum.freecad.org/viewtopic.php?t=86001

    The reason of the crash is that the GIL is not locked while accessing certain Python functions

commit 2c33a24
Author: bgbsww <bgbsww@gmail.com>
Date:   Thu Mar 14 10:24:32 2024 -0400

    Toponaming/Part: Review and lint cleanups

commit 9830a62
Author: bgbsww <bgbsww@gmail.com>
Date:   Tue Mar 12 11:20:56 2024 -0400

    Toponaming/Part: clean and test attacher

commit 21d6a67
Author: Zheng, Lei <realthunder.dev@gmail.com>
Date:   Mon Mar 11 11:13:27 2024 -0400

    Toponaming/Part:  Bring over attacher differences

commit 21729eb
Author: FEA-eng <59876896+FEA-eng@users.noreply.github.com>
Date:   Thu Mar 14 19:31:17 2024 +0100

    DRAFT: Update preferences-draftsnap.ui

commit 2e18429
Author: Chris Hennes <chennes@pioneerlibrarysystem.org>
Date:   Thu Mar 14 22:23:02 2024 -0500

    FEM: Add two missing PreCompiled.h includes

commit 05f2f97
Author: FEA-eng <59876896+FEA-eng@users.noreply.github.com>
Date:   Fri Mar 15 04:15:29 2024 +0100

    Start wb FEM examples update (FreeCAD#12871)

    * FEM: Delete data/examples/FemCalculixCantilever3D_newSolver.FCStd

    * FEM: Delete data/examples/FemCalculixCantilever3D.FCStd

    * FEM: Delete data/examples/FemCalculixCantilever2D.FCStd

    * FEM: add updated cantilever 3D

    * FEM: add updated cantilever 2D

    * FEM: add cantilever 1D

    * FEM: Update CMakeLists.txt

    * FEM: Delete data/examples/FemCalculixCantilever1D.FCStd

    * FEM: Delete data/examples/FemCalculixCantilever2D.FCStd

    * FEM: Delete data/examples/FemCalculixCantilever3D.FCStd

    * FEM: reupload the 1D example with a new license

    * FEM: reupload the 2D example with a new license

    * FEM: reupload the 3D example with a new license

commit 9a79aa8
Author: David Carter <dcarter@davidcarter.ca>
Date:   Wed Mar 13 17:42:53 2024 -0400

    Material: Material editor enhancements

    Removes localization from the material card dictionary.

    Fixes FreeCAD#12935

    The previous card implementation transferred strings directly to the
    card dictionary without interpretation. The new material system parses
    the file data, producing a dictionary for compatibility. The new
    dictionaries localized the strings which is incorrect. This PR
    removes that localization.

commit 46b753a
Author: FEA-eng <59876896+FEA-eng@users.noreply.github.com>
Date:   Thu Mar 14 20:39:16 2024 +0100

    FEM: Update solver.py

commit 1b1b232
Author: wandererfan <wandererfan@gmail.com>
Date:   Thu Mar 14 14:22:29 2024 -0400

    [TD]fix over/under tolerance

    - left justify tolerances
    - use format spec for tolerance == zero
    - prevent uncommanded tolerance format changes
    - use 'w' format for HoleShaftFit
    - allow zero tolerance values for HoleShaftFit

commit 1411e70
Author: wandererfan <wandererfan@gmail.com>
Date:   Thu Mar 14 12:45:19 2024 -0400

    [TD]fix thread line alignment

commit f98a91f
Author: wandererfan <wandererfan@gmail.com>
Date:   Thu Mar 14 12:42:17 2024 -0400

    [TD]prevent fail on delete with open dialog

commit d2e27b3
Author: sliptonic <shopinthewoods@gmail.com>
Date:   Thu Mar 14 11:37:32 2024 -0500

    mark dynamic fields to not be translated

commit 755f2ec
Author: hoshengwei <hoshengwei@gmail.com>
Date:   Wed Mar 13 21:15:57 2024 +0800

    Arch: Fixed bug in BuildingPart obj

    Fixes FreeCAD#12921.

commit 7e3fc8d
Author: marioalexis <mario.passaglia@gmail.com>
Date:   Wed Mar 13 16:41:58 2024 -0300

    Fem: Update force constraint .inp files

commit 3c4cbac
Author: marioalexis <mario.passaglia@gmail.com>
Date:   Wed Mar 13 14:50:00 2024 -0300

    Fem: Format force constraint in CalculiX writer - fixes FreeCAD#12805

commit 50b2e09
Author: luzpaz <luzpaz@users.noreply.github.com>
Date:   Wed Mar 13 18:10:04 2024 -0400

    labels.yml: Add Materials label

    Auto-assign 'Materials' label to Material related PRs

commit ecfc88b
Author: marioalexis <mario.passaglia@gmail.com>
Date:   Sat Mar 9 15:30:11 2024 -0300

    Fem: Partial transparency fix in FemPostObject display modes

commit 9d7cf7e
Author: marioalexis <mario.passaglia@gmail.com>
Date:   Sat Mar 9 15:05:53 2024 -0300

    Fem: Export/import mesh elements to Vtk using proper order

commit 53d4069
Author: David Carter <dcarter@davidcarter.ca>
Date:   Sat Mar 16 01:13:20 2024 -0400

    Build fixes

commit 0ba77b1
Merge: 1b22ed6 cc2412e
Author: David Carter <dcarter@davidcarter.ca>
Date:   Fri Mar 15 22:04:02 2024 -0400

    Merge branch 'materials_04022024' into materials_15032024

commit 1b22ed6
Author: David Carter <dcarter@davidcarter.ca>
Date:   Thu Mar 14 00:45:02 2024 -0400

    Material: Call to import_materials only returning solids

    Fixes FreeCAD#12933

    The previous card implementation used a parameter to select either
    solids or fluids. The current material dictionary ignored the
    parameter and only returned solids. This PR restores the use of the
    material type filter.

commit cc2412e
Author: David Carter <dcarter@davidcarter.ca>
Date:   Tue Mar 12 21:37:45 2024 -0400

    Compatibility with perevious versions

commit e8ed41d
Author: David Carter <dcarter@davidcarter.ca>
Date:   Sun Mar 3 21:01:34 2024 -0500

    Revert "Ondsel solver update"

    This reverts commit 6db79f0.

commit b0e6e96
Author: David Carter <dcarter@davidcarter.ca>
Date:   Sat Mar 2 01:19:47 2024 -0500

    ShapeMaterial becomes part of Part::Feature

commit 20444f1
Author: David Carter <dcarter@davidcarter.ca>
Date:   Mon Feb 26 23:29:01 2024 -0500

    Set material option

commit f6ed92d
Author: David Carter <dcarter@davidcarter.ca>
Date:   Mon Feb 26 08:21:48 2024 -0500

    Add uuid to App/Material

commit a85f500
Author: David Carter <dcarter@davidcarter.ca>
Date:   Fri Feb 23 13:54:13 2024 -0500

    Use materials to set appearance

commit 9ae37a6
Merge: 957c531 1303d70
Author: David Carter <dcarter@davidcarter.ca>
Date:   Thu Feb 22 15:12:37 2024 -0500

    Merge branch 'materials_04022024' of https://github.com/davesrocketshop/FreeCAD into materials_04022024

commit 957c531
Author: David Carter <dcarter@davidcarter.ca>
Date:   Thu Feb 22 15:12:31 2024 -0500

    Remove unused Python code

commit 1303d70
Author: David Carter <dcarter@davidcarter.ca>
Date:   Thu Feb 22 15:10:47 2024 -0500

    Move appearance settings to materials

commit 45bd641
Author: David Carter <dcarter@davidcarter.ca>
Date:   Thu Feb 22 10:25:43 2024 -0500

    Move Std_SetAppearnce

commit e3660e3
Author: David Carter <dcarter@davidcarter.ca>
Date:   Tue Feb 20 17:25:53 2024 -0500

    Add custom attributes

commit bfa694f
Author: David Carter <dcarter@davidcarter.ca>
Date:   Sun Feb 18 11:01:34 2024 -0500

    ShapeAppearance to PropertyMaterialList

commit ec4d37d
Author: David Carter <dcarter@davidcarter.ca>
Date:   Wed Feb 14 10:35:11 2024 -0500

    ShapeColor to ShapeMaterial

commit 17f567f
Author: David Carter <dcarter@davidcarter.ca>
Date:   Wed Feb 14 06:08:08 2024 -0500

    Fix remaining test cases

commit 9297f89
Author: David Carter <dcarter@davidcarter.ca>
Date:   Wed Feb 14 03:42:30 2024 -0500

    Pass more unit tests

commit 0767454
Author: David Carter <dcarter@davidcarter.ca>
Date:   Tue Feb 13 21:40:38 2024 -0500

    Filter changes

commit ee72cb0
Author: David Carter <dcarter@davidcarter.ca>
Date:   Tue Feb 13 21:38:33 2024 -0500

    Creaate and use ShapeMaterial

commit 4fdcf80
Author: David Carter <dcarter@davidcarter.ca>
Date:   Sun Feb 4 17:40:55 2024 -0500

    FIlter materials by model
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Assignees
No one assigned
Labels
Bug This issue or PR is related to a bug WB FEM Related to the FEM Workbench
Projects
FEM Roadmap
Status: Done
Milestone
No milestone
Development

No branches or pull requests
3 participants
@maxwxyz @marioalexis84 @FEA-eng