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

Fem: Update writer of the constraint rigid body

Browse Source
This commit is contained in:
FEA-eng
2024-03-15 19:12:12 +01:00
committed by Chris Hennes
parent ab6afeb261
commit 3739d06528
6 changed files with 105 additions and 19 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
src/Mod/Fem/App/FemConstraintRigidBody.cpp
View File

@@ -77,7 +77,7 @@ void ConstraintRigidBody::onChanged(const App::Property* prop)
if (prop == &References) {
std::vector<Base::Vector3d> points;
std::vector<Base::Vector3d> normals;
int scale = 1; // OvG: Enforce use of scale
double scale = 1; // OvG: Enforce use of scale
if (getPoints(points, normals, &scale)) {
Points.setValues(points);
Normals.setValues(normals);

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

@@ -232,6 +232,7 @@ SET(FemSolverCalculix_SRCS
femsolver/calculix/write_constraint_planerotation.py
femsolver/calculix/write_constraint_pressure.py
femsolver/calculix/write_constraint_rigidbody.py
femsolver/calculix/write_constraint_rigidbody_step.py
femsolver/calculix/write_constraint_sectionprint.py
femsolver/calculix/write_constraint_selfweight.py
femsolver/calculix/write_constraint_temperature.py

28
src/Mod/Fem/femsolver/calculix/write_constraint_rigidbody.py
View File

@@ -66,9 +66,9 @@ def write_constraint(f, femobj, rb_obj, ccxwriter):
is_ref_bc_defined = any((rb_obj.xDisplacement,
rb_obj.yDisplacement,
rb_obj.zDisplacement,
rb_obj.xLoad,
rb_obj.yLoad,
rb_obj.zLoad))
rb_obj.xForce,
rb_obj.yForce,
rb_obj.zForce))
is_rot_bc_defined = any((rb_obj.xRotation,
rb_obj.yRotation,
@@ -78,10 +78,14 @@ def write_constraint(f, femobj, rb_obj, ccxwriter):
rb_obj.zMoment))
# FIXME: This needs to be implemented
ref_node_idx = -1
rot_node_idx = -1
ref_node_idx = 10000000
rot_node_idx = 20000000
kw_line = "*RIGID BODY,NSET={}".format(rb_obj.Name)
f.write("*NODE\n")
f.write("{},{},{},{}\n".format(ref_node_idx, rb_obj.xRefNode, rb_obj.yRefNode, rb_obj.zRefNode))
f.write("{},{},{},{}\n".format(rot_node_idx, rb_obj.xRefNode, rb_obj.yRefNode, rb_obj.zRefNode))
kw_line = "*RIGID BODY, NSET={}, REF NODE={}, ROT NODE={}".format(rb_obj.Name, ref_node_idx, rot_node_idx)
if is_ref_bc_defined:
kw_line = kw_line + ",REF NODE={}".format(ref_node_idx)
@@ -91,15 +95,3 @@ def write_constraint(f, femobj, rb_obj, ccxwriter):
f.write(kw_line + "\n")
# TODO: Displacement definitions need fixing
if is_ref_bc_defined:
f.write("*CLOAD\n")
f.write("{},1,{}\n".format(ref_node_idx, rb_obj.xLoad))
f.write("{},2,{}\n".format(ref_node_idx, rb_obj.yLoad))
f.write("{},3,{}\n".format(ref_node_idx, rb_obj.zLoad))
if is_rot_bc_defined:
f.write("*CLOAD\n")
f.write("{},1,{}\n".format(rot_node_idx, rb_obj.xMoment))
f.write("{},2,{}\n".format(rot_node_idx, rb_obj.yMoment))
f.write("{},3,{}\n".format(rot_node_idx, rb_obj.zMoment))

88
src/Mod/Fem/femsolver/calculix/write_constraint_rigidbody_step.py Normal file
View File

@@ -0,0 +1,88 @@
# ***************************************************************************
# * Copyright (c) 2022 Ajinkya Dahale <dahale.a.p@gmail.com> *
# * *
# * 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 calculix constraint rigid body"
__author__ = "Ajinkya Dahale"
__url__ = "https://www.freecadweb.org"
def get_analysis_types():
return "all" # write for all analysis types
def get_sets_name():
return "constraints_rigidbody_node_sets"
def get_constraint_title():
return "Rigid Body Constraints"
def get_before_write_meshdata_constraint():
return ""
def get_after_write_meshdata_constraint():
return ""
def get_before_write_constraint():
return ""
def get_after_write_constraint():
return ""
def write_constraint(f, femobj, rb_obj, ccxwriter):
# floats read from ccx should use {:.13G}, see comment in writer module
is_ref_bc_defined = any((rb_obj.xDisplacement,
rb_obj.yDisplacement,
rb_obj.zDisplacement,
rb_obj.xForce,
rb_obj.yForce,
rb_obj.zForce))
is_rot_bc_defined = any((rb_obj.xRotation,
rb_obj.yRotation,
rb_obj.zRotation,
rb_obj.xMoment,
rb_obj.yMoment,
rb_obj.zMoment))
# FIXME: This needs to be implemented
ref_node_idx = 10000000
rot_node_idx = 20000000
# TODO: Displacement definitions need fixing
f.write("*CLOAD\n")
f.write("{},1,{}\n".format(ref_node_idx, rb_obj.xForce))
f.write("{},2,{}\n".format(ref_node_idx, rb_obj.yForce))
f.write("{},3,{}\n".format(ref_node_idx, rb_obj.zForce))
f.write("*CLOAD\n")
f.write("{},1,{}\n".format(rot_node_idx, rb_obj.xMoment))
f.write("{},2,{}\n".format(rot_node_idx, rb_obj.yMoment))
f.write("{},3,{}\n".format(rot_node_idx, rb_obj.zMoment))

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

@@ -46,6 +46,7 @@ from . import write_constraint_initialtemperature as con_itemp
from . import write_constraint_planerotation as con_planerotation
from . import write_constraint_pressure as con_pressure
from . import write_constraint_rigidbody as con_rigidbody
from . import write_constraint_rigidbody_step as con_rigidbody_step
from . import write_constraint_sectionprint as con_sectionprint
from . import write_constraint_selfweight as con_selfweight
from . import write_constraint_temperature as con_temperature
@@ -190,6 +191,7 @@ class FemInputWriterCcx(writerbase.FemInputWriter):
# constraints dependent from steps
self.write_constraints_propdata(inpfile, self.member.cons_fixed, con_fixed)
self.write_constraints_propdata(inpfile, self.member.cons_rigidbody_step, con_rigidbody_step)
self.write_constraints_propdata(inpfile, self.member.cons_displacement, con_displacement)
self.write_constraints_propdata(inpfile, self.member.cons_sectionprint, con_sectionprint)
self.write_constraints_propdata(inpfile, self.member.cons_selfweight, con_selfweight)

3
src/Mod/Fem/femtools/membertools.py
View File

@@ -301,6 +301,9 @@ class AnalysisMember():
self.cons_rigidbody = self.get_several_member(
"Fem::ConstraintRigidBody"
)
self.cons_rigidbody_step = self.get_several_member(
"Fem::ConstraintRigidBody"
)
self.cons_force = self.get_several_member(
"Fem::ConstraintForce"
)