forked from FreeCAD/FreeCAD
-
Notifications
You must be signed in to change notification settings - Fork 1
/
write_constraint_rigidbody.py
105 lines (77 loc) · 3.95 KB
/
write_constraint_rigidbody.py
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
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
# ***************************************************************************
# * 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_meshdata_constraint(f, femobj, rb_obj, ccxwriter):
f.write("*NSET,NSET=" + rb_obj.Name + "\n")
for n in femobj["Nodes"]:
f.write("{},\n".format(n))
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.xLoad,
rb_obj.yLoad,
rb_obj.zLoad))
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 = -1
rot_node_idx = -1
kw_line = "*RIGID BODY,NSET={}".format(rb_obj.Name)
if is_ref_bc_defined:
kw_line = kw_line + ",REF NODE={}".format(ref_node_idx)
if is_rot_bc_defined:
kw_line = kw_line + ",ROT NODE={}".format(rot_node_idx)
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))