### ===================================================================================================================
### Tool for creating IFC of a model
### ===================================================================================================================
# Copyright ©VIIA 2025
### ===================================================================================================================
### 1. Import modules
### ===================================================================================================================
# General references
from pathlib import Path
from typing import List, Optional
from math import pi
# References for functions and classes in the rhdhv_fem package
from rhdhv_fem.materials import DirectStiffnessMatrixModel, LinearElasticOrthotropicModel, SpringBehaviour, NoCapacity,\
Masonry, Concrete
from rhdhv_fem.geometries import TrussProfile, Rectangle
from rhdhv_fem.loads import ModalPushOverLoad
from rhdhv_fem.supports import SurfaceSupport
### ===================================================================================================================
### 2. viia_create_scia_model_for_ifc_export
### ===================================================================================================================
[docs]def viia_create_scia_model_for_ifc_export(
json: Path = None, diana_dat_file: Path = None, items_to_remove: Optional[List[str]] = None):
"""
Create a scia xml file that can be used to export an IFC model in SCIA Engineering.
Input:
- json (obj): Path of a json file with the model. Default is None.
- diana_dat_file (obj): Path of a diana dat file with the model. Default is None.
- items_to_remove (list of str): List with name of items that should be removed to be able to create the model
in SCIA. Can be used if the automatic update of the model is not sufficient.
Output:
- A xml will be created with the model that can be exported as IFC.
"""
from viiapackage.viiaStatus import viia_create_project
if (json is None and diana_dat_file is None) or (json is not None and diana_dat_file is not None):
raise ValueError(
f"ERROR: One of the arguments 'json' and 'diana_dat_file' should be not None. Provided was {json=} and "
f"{diana_dat_file=}.")
project = viia_create_project(f'test-ifc_from_json')
if json is not None:
project.viia_read_dump(json)
elif diana_dat_file is not None:
project.from_diana(diana_dat_file)
else:
raise FileNotFoundError(f"ERROR: File not found. Provided was {json=} and {diana_dat_file=}.")
# Remove items that cannot be made in scia
removed_items = []
for analysis in reversed(project.collections.analyses):
removed_items.append(analysis.name)
analysis.remove_analysis()
for spectrum in reversed(project.collections.response_spectra):
removed_items.append(spectrum.name)
project.collections.response_spectra = []
for shape in reversed(project.collections.line_masses):
removed_items.append(shape.name)
shape.remove_shape()
for shape in reversed(project.collections.springs):
removed_items.append(shape.name)
shape.remove_connection()
for shape in reversed(project.collections.interfaces + project.collections.boundary_interfaces):
removed_items.append(shape.name)
shape.remove_connection()
for load in reversed(project.collections.loads):
if isinstance(load, ModalPushOverLoad):
removed_items.append(load.name)
load.remove_load()
for support in reversed(project.collections.supports):
if isinstance(support, SurfaceSupport):
removed_items.append(support.name)
support.remove_support()
# fix for line masses from dat file
LIJNMASSA = [shape for shape in project.collections.shapes if 'LIJNMASSA' in shape.name]
for line_mass in reversed(LIJNMASSA):
removed_items.append(line_mass.name)
line_mass.remove_shape()
concrete = [material for material in project.collections.materials if isinstance(material, Concrete)]
if concrete:
concrete = concrete[0]
for line in project.collections.lines:
if isinstance(line.geometry.geometry_model, Rectangle) and isinstance(line.material, Masonry):
if not concrete:
concrete = project.create_material('c20/25')
line.material = concrete
project.write_log(
f"The material of line {line} is changed from masonry to concrete so it can be made in SCIA.")
for geometry in reversed(project.collections.geometries):
if isinstance(geometry.geometry_model, TrussProfile):
diameter = 2 * ((geometry.geometry_model.cross_sectional_area/pi)**(1/2))
geometry.geometry_model = project.create_circle_profile(
profile_name=geometry.geometry_model.name,
diameter=diameter)
project.write_log(
f"The geometry model for geometry {geometry} is changed to circular so it can be made "
f"in SCIA.")
for material in reversed(project.collections.materials):
if isinstance(material.material_model, DirectStiffnessMatrixModel):
membrane_stiffness_matrix = list(material.material_model.membrane_stiffness_matrix.values())
membrane_stiffness_matrix = sum(membrane_stiffness_matrix)/len(membrane_stiffness_matrix)
shear_stiffness_matrix = list(material.material_model.shear_stiffness_matrix.values())
shear_stiffness_matrix = sum(shear_stiffness_matrix) / len(shear_stiffness_matrix)
poissons_ratio = membrane_stiffness_matrix / (2*shear_stiffness_matrix) - 1
poissons_ratio = min(poissons_ratio, 0)
material.material_model = project.create_linear_elastic_isotropic_model(
youngs_modulus=membrane_stiffness_matrix, poissons_ratio=poissons_ratio)
project.write_log(
f"The material model for material {material} is changed to linear elastic isotropic so it can be made "
f"in SCIA.")
elif isinstance(material.material_model, LinearElasticOrthotropicModel):
youngs_moduli = list(material.material_model.youngs_moduli.values())
youngs_moduli = sum(youngs_moduli)/len(youngs_moduli)
poissons_ratios = list(material.material_model.poissons_ratios.values())
poissons_ratios = sum(poissons_ratios) / len(poissons_ratios)
material.material_model = project.create_linear_elastic_isotropic_model(
youngs_modulus=youngs_moduli, poissons_ratio=poissons_ratios)
project.write_log(
f"The material model for material {material} is changed to linear elastic isotropic so it can be made "
f"in SCIA.")
elif isinstance(material, SpringBehaviour):
removed_items.append(material.name)
material.remove_material()
# Fix for material without capacities
concrete_without_capacity = [
material
for material in project.collections.materials
if (material.capacity is None or isinstance(material.capacity, NoCapacity)) and
('BETON' in material.name or 'KPV' in material.name)
]
if concrete_without_capacity:
capacity = project.create_EC_concrete_capacity(strength_class='C20/25')
for material in concrete_without_capacity:
project.write_log(f"Concrete capacity added for material {material}.")
material.capacity = capacity
if items_to_remove:
removed_items = []
for item in items_to_remove:
_item = project.find(description=item, collection='all_collections')
if _item is not None:
for method_name in [
'remove_shape', 'remove_connection', 'remove_load', 'remove_material', 'remove_support']:
if hasattr(_item, method_name):
removed_items.append(item)
func = getattr(_item, method_name)
func()
break
else:
project.write_log(f"For item {_item} no remove method is found.")
not_removed_items = [item for item in items_to_remove if item not in removed_items]
if not_removed_items:
project.write_log(f"The following items could not be removed:")
for item in not_removed_items:
project.write_log(f"\t{item}")
project.purge_shape_geometries()
project.write_dump()
return project.create_model(software='scia')
### ===================================================================================================================
### 3. End of script
### ===================================================================================================================