94 lines
4.0 KiB
Python
94 lines
4.0 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 write inpfile mesh"
|
|
__author__ = "Bernd Hahnebach"
|
|
__url__ = "https://www.freecad.org"
|
|
|
|
|
|
import codecs
|
|
from os.path import join
|
|
|
|
from femmesh import meshtools
|
|
|
|
|
|
def write_mesh(ccxwriter):
|
|
|
|
element_param = 1 # highest element order only
|
|
group_param = False # do not write mesh group data
|
|
|
|
# Use reduced integration beam elements if this option is enabled in ccx solver settings
|
|
vol_variant = "standard"
|
|
edge_variant = "beam"
|
|
if ccxwriter.solver_obj.BeamReducedIntegration:
|
|
edge_variant = "beam reduced"
|
|
# Check to see if fluid sections are in analysis and use D network element type
|
|
if ccxwriter.member.geos_fluidsection:
|
|
edge_variant = "network"
|
|
|
|
# Use 2D elements if model space is not set to 3D
|
|
if ccxwriter.solver_obj.ModelSpace == "3D":
|
|
face_variant = "shell"
|
|
elif ccxwriter.solver_obj.ModelSpace == "plane stress":
|
|
face_variant = "stress"
|
|
elif ccxwriter.solver_obj.ModelSpace == "plane strain":
|
|
face_variant = "strain"
|
|
elif ccxwriter.solver_obj.ModelSpace == "axisymmetric":
|
|
face_variant = "axisymmetric"
|
|
|
|
if ccxwriter.split_inpfile:
|
|
write_name = "femesh"
|
|
file_name_split = ccxwriter.mesh_name + "_" + write_name + ".inp"
|
|
ccxwriter.femmesh_file = join(ccxwriter.dir_name, file_name_split)
|
|
|
|
ccxwriter.femmesh.writeABAQUS(
|
|
ccxwriter.femmesh_file,
|
|
element_param,
|
|
group_param,
|
|
volVariant=vol_variant,
|
|
faceVariant=face_variant,
|
|
edgeVariant=edge_variant,
|
|
)
|
|
|
|
inpfile = codecs.open(ccxwriter.file_name, "w", encoding="utf-8")
|
|
inpfile.write("{}\n".format(59 * "*"))
|
|
inpfile.write(f"** {write_name}\n")
|
|
inpfile.write(f"*INCLUDE,INPUT={file_name_split}\n")
|
|
|
|
else:
|
|
ccxwriter.femmesh_file = ccxwriter.file_name
|
|
ccxwriter.femmesh.writeABAQUS(
|
|
ccxwriter.femmesh_file,
|
|
element_param,
|
|
group_param,
|
|
volVariant=vol_variant,
|
|
faceVariant=face_variant,
|
|
edgeVariant=edge_variant,
|
|
)
|
|
|
|
# reopen file with "append" to add all the rest
|
|
inpfile = codecs.open(ccxwriter.femmesh_file, "a", encoding="utf-8")
|
|
inpfile.write("\n\n")
|
|
|
|
return inpfile
|