107 lines
4.1 KiB
Python
107 lines
4.1 KiB
Python
# ***************************************************************************
|
|
# * 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 calculix constraint planerotation"
|
|
__author__ = "Bernd Hahnebach"
|
|
__url__ = "https://www.freecad.org"
|
|
|
|
|
|
from femmesh import meshtools
|
|
|
|
|
|
def get_analysis_types():
|
|
return "all" # write for all analysis types
|
|
|
|
|
|
def get_sets_name():
|
|
return "constraints_planerotation_node_sets"
|
|
|
|
|
|
def get_constraint_title():
|
|
return "PlaneRotation 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, fric_obj, ccxwriter):
|
|
# write nodes to file
|
|
if not ccxwriter.femnodes_mesh:
|
|
ccxwriter.femnodes_mesh = ccxwriter.femmesh.Nodes
|
|
# info about ccxwriter.constraint_conflict_nodes:
|
|
# is used to check if MPC and constraint fixed and
|
|
# constraint displacement share same nodes
|
|
# because MPC"s and constraints fixed and
|
|
# constraints displacement can't share same nodes.
|
|
# Thus call write_node_sets_constraints_planerotation has to be
|
|
# after constraint fixed and constraint displacement
|
|
l_nodes = femobj["Nodes"]
|
|
f.write(f"*NSET,NSET={fric_obj.Name}\n")
|
|
# Code to extract nodes and coordinates on the PlaneRotation support face
|
|
nodes_coords = []
|
|
for node in l_nodes:
|
|
nodes_coords.append(
|
|
(
|
|
node,
|
|
ccxwriter.femnodes_mesh[node].x,
|
|
ccxwriter.femnodes_mesh[node].y,
|
|
ccxwriter.femnodes_mesh[node].z,
|
|
)
|
|
)
|
|
node_planerotation = meshtools.get_three_non_colinear_nodes(nodes_coords)
|
|
for i in range(len(l_nodes)):
|
|
if l_nodes[i] not in node_planerotation:
|
|
node_planerotation.append(l_nodes[i])
|
|
MPC_nodes = []
|
|
for i in range(len(node_planerotation)):
|
|
cnt = 0
|
|
for j in range(len(ccxwriter.constraint_conflict_nodes)):
|
|
if node_planerotation[i] == ccxwriter.constraint_conflict_nodes[j]:
|
|
cnt = cnt + 1
|
|
if cnt == 0:
|
|
MPC = node_planerotation[i]
|
|
MPC_nodes.append(MPC)
|
|
for i in range(len(MPC_nodes)):
|
|
f.write(f"{MPC_nodes[i]},\n")
|
|
|
|
|
|
def write_constraint(f, femobj, fric_obj, ccxwriter):
|
|
|
|
# floats read from ccx should use {:.13G}, see comment in writer module
|
|
|
|
f.write("*MPC\n")
|
|
f.write(f"PLANE,{fric_obj.Name}\n")
|