MFE_HAL.py

import os
import sys
baseName = os.path.basename(__file__)
dirName = os.path.dirname(__file__)
print('basename:    ', baseName)
print('dirname:     ', dirName)
sys.path.append(dirName + r'/../..')

sys.stdout.write("\r                                                                       ")
sys.stdout.write("\rIMPORTEREN BIBLIOTHEKEN\n")
sys.stdout.flush()

# dict = {'A':'1','B':'2','C':'3','D':'4','E':'5','F':'6','G':'7','H':'8',\
#         'I':'9','J':'10','K':'11','L':'12','M':'13','N':'14','O':'15','P':'16',\
#         'Q':'17','R':'18','S':'19','T':'20','U':'21','V':'22','W':'23','X':'24',\
#         'Y':'25','Z':'26','AA':'27','AB':'28','AC':'29','AD':'30','AE':'31','AF':'32',\
#         'AG':'33','AH':'34','AI':'35','AJ':'36','AK':'37','AL':'38','AM':'39','AN':'40'}

import time
time1 = time.time()

from RFEM.enums import *
from RFEM.initModel import Model, Calculate_all, CalculateSelectedCases, ConvertToDlString, SetAddonStatus, FirstFreeIdNumber, insertSpaces
from RFEM.baseSettings import BaseSettings
from RFEM.BasicObjects.material import Material
from RFEM.BasicObjects.section import Section
from RFEM.BasicObjects.node import Node
from RFEM.BasicObjects.line import Line
from RFEM.BasicObjects.member import Member
from RFEM.BasicObjects.memberSet import MemberSet
from RFEM.BasicObjects.surface import Surface
from RFEM.TypesForNodes.nodalSupport import NodalSupport
from RFEM.TypesForMembers.memberHinge import MemberHinge
# from RFEM.TypesForMembers.memberNonlinearity import MemberNonlinearity
from RFEM.LoadCasesAndCombinations.staticAnalysisSettings import StaticAnalysisSettings
from RFEM.LoadCasesAndCombinations.designSituation import clearAttributes, DesignSituation
from RFEM.LoadCasesAndCombinations.loadCase import LoadCase
from RFEM.LoadCasesAndCombinations.loadCombination import LoadCombination
from RFEM.LoadCasesAndCombinations.loadCasesAndCombinations import LoadCasesAndCombinations
from RFEM.LoadCasesAndCombinations.combinationWizard import CombinationWizard
from RFEM.Loads.memberLoad import MemberLoad
from RFEM.Calculate.meshSettings import GetModelInfo
# from RFEM.ImportExport.exports import ExportDetailsOfDesignToCSV
from RFEM.dataTypes import inf
from RFEM.SteelDesign.steelUltimateConfigurations import SteelDesignUltimateConfigurations
from RFEM.SteelDesign.steelServiceabilityConfiguration import SteelDesignServiceabilityConfigurations
from RFEM.TypesForSteelDesign.steelEffectiveLengths import SteelEffectiveLengths
from RFEM.Results.resultTables import ResultTables, GetMaxValue, GetMinValue
from RFEM.Results.designOverview import GetDesignOverview, GetPartialDesignOverview

from MFE_ALG import *
# import MFE_KolomLetter
import MFE_ZoekNode
import MFE_getMembers
import MFE_getNodes
from MFE_Wind_Walls_and_Flat_Roof import Wind
from MFE_Sneeuw import Sneeuw

import math

# print(ord("a"))
EffLengthMembers = [[],[],[],[],[]]
EffLengthSets = [[],[],[],[],[]]
DesignPropsViaParentSet = bool


def BeamLimitCharacteristic(value: int = 100, sls_no = 1):

    with open(dirName+r"./sls_char.js", "w") as std:

        std.write("STEEL_DESIGN.steel_design_sls_configurations[{}].settings_ec3.property_sl_beam_limit_characteristic = {}".format(sls_no, value))

    Model.clientModel.service.run_script(dirName+r"./sls_char.js")

    os.remove(dirName+r"./sls_char.js")

def SnowWizardMonopitch(nodeList= [1,2,3,4]):

    with open(dirName+r"./SnowWizard.js", "w") as snowwiz:

        snowwiz.write("var snowLoadWizard = new SnowLoadWizard();\n")
        snowwiz.write('snowLoadWizard.SetMonoPitchRoofType(undefined, [4, 9, 57, 52], load_cases[300])\n')

    Model.clientModel.service.run_script(dirName+r"./SnowWizard.js")\

    os.remove(dirName+r"./SnowWizard.js")

def MemberToLineList(ListIn = [1,2,3,4], members=[{"no":1}]):
    listOut = []
    for m in ListIn:
        member = next(item for item in members if item["no"] == int(m))
        listOut.append(member["line"])
    return listOut


if __name__ == '__main__':
    dy = 5 #float(input("H.o.h. afstand tussen assen // x-as [m]: "))
    dx = 5 #float(input("H.o.h. afstand tussen assen // y-as [m]: "))
    nx = 8 #int(input("Aantal assen in x-richting: "))
    ny = 5 #int(input("Aantal assen in y-richting: "))
    h = 14.4 #float(input("Hoogte hal incl. dakrand [m]: "))
    h_dr = 1.2 #float(input("Hoogte dakrand t.o.v. hart dakliggers [m]: "))
    dhr = 0.8 #float(input("Hoogteverschil bovenkant dakrand en bovenkant dakafwerking [m]: "))
    kst_kol = 2 #int(input("Aantal kniksteunen van de kolommen: "))
    nodes_frame = 2*(3+kst_kol)+ny-2
    members_frame = 2*(2+kst_kol)+ny-1
    sup_nodes = []
    b = (ny-1)*dy
    L = (nx-1)*dx
    GewichtGevel = 0.3 # Gewicht gevel in kN/m²
    PB_Dak = 0.5 # Gewicht dakplaten en dakbedekking in kN/m²

    qpz = 0.97 # "Winddruk in kN/m2 : "

    nodes=[]

    Model(True, "MFE_HAL")
    Model.clientModel.service.delete_all()
    Model.clientModel.service.begin_modification()
    BaseSettings(9.81, GlobalAxesOrientationType.E_GLOBAL_AXES_ORIENTATION_ZUP, LocalAxesOrientationType.E_LOCAL_AXES_ORIENTATION_ZDOWN,[0.01,0.01,0.01,0.01])

    SetAddonStatus(Model.clientModel, AddOn.steel_design_active, True)

    Model.clientModel.service.finish_modification(); Model.clientModel.service.begin_modification()
    # Model.clientModel.service.begin_modification()

    Material(1, 'S235')
    # Section(1, 'HEA 240 | Euronorm 53-62; ... | SZS', 1, "Profiel Kolommen 1")
    # Section(2, 'HEA 200 | Euronorm 53-62; ... | SZS', 1, "Profiel Kolommen 2")
    Section(1, 'IPE 330 | Euronorm 19-57; ... | SZS', 1, "Profiel Kolommen 1")
    Section(2, 'IPE 300 | Euronorm 19-57; ... | SZS', 1, "Profiel Kolommen 2")
    # Section(2, 'HEA 200 | Euronorm 53-62; ... | SZS', 1, "Profiel Kolommen 2")
    Section(3, 'IPE 400 | Euronorm 19-57; ... | SZS', 1, "Profiel Hoofdliggers")
    Section(4, 'HEA 160 | Euronorm 53-62; ... | SZS', 1, "Profiel Dakliggers KOPSE GEVELS")
    Section(5, 'HEA 140 | Euronorm 53-62; ... | SZS', 1, "Profiel Dakliggers LANGSGEVELS")
    Section(6, 'SHS 100x5 | EN 10210-2:2006-04', 1, "Profiel Kipsteunen dak")
    Section(7, 'SHS 80x4 | EN 10210-2:2006-04', 1, "Profiel Horizontale gevelliggers KOPSE GEVELS")
    Section(8, 'HEA 100 | Euronorm 53-62; ... | SZS', 1, "Profiel Horizontale gevelliggers LANGSGEVELS")

    Section(9, 'HEA 100 | Euronorm 53-62; ... | SZS', 1, "Profiel Dakrand Verticaal")
    Section(10, 'HEA 100 | Euronorm 53-62; ... | SZS', 1, "Profiel Dakrand Horizontaal")
    Section(11, 'L 90x90x9', 1, "Schoren Dak 1")
    Section(12, 'L 90x90x9', 1, "Schoren Dak 2")
    Section(13, 'FL 150x10', 1, "Schoren Gevel 1")
    Section(14, 'FL 100x8', 1, "Schoren Gevel 2")

    MemberHinge(1,"Local","",inf,inf,inf,inf,0,0)
    MemberHinge(2,"Local","",inf,inf,inf,0,0,0)

    #MemberNonlinearity(1,"")

#----------------------------------------------------------------------------------------------
#
#                       GEOMETRIE
#
#----------------------------------------------------------------------------------------------
    sys.stdout.write("\r                                                                       ")
    sys.stdout.write("\rGEOMETRIE\n")
    sys.stdout.flush()

    # Lijsten aanmaken voor staven en lijnen waarmee uiteindelijk ook de Load Transfer 2D elementen worden gedefinieerd:
    MemListRandLiggerX0 = []
    MemListRandLiggerXL = []
    MemListRandLiggerY0 = []
    MemListRandLiggerYB = []
    MemListKol00 = []
    MemListKol0B = []
    MemListKolL0 = []
    MemListKolLB = []
    MemListDak = []
    MemListX0 = []
    MemListXL = []
    MemListY0 = []
    MemListYB = []

    MemListDakRandX0 = []
    MemListDakRandXL = []
    MemListDakRandY0 = []
    MemListDakRandYB = []

    LinListRandLiggerX0 = []
    LinListRandLiggerXL = []
    LinListRandLiggerY0 = []
    LinListRandLiggerYB = []
    LinListKol00 = []
    LinListKol0B = []
    LinListKolL0 = []
    LinListKolLB = []

    linesMV_Y0 = [] #Lijst met lijnen grondvlak tbv loadtransfer surfactelement op Y=0
    linesMV_YB = [] #Lijst met lijnen grondvlak tbv loadtransfer surfactelement op Y=B
    linesMV_X0 = [] #Lijst met lijnen grondvlak tbv loadtransfer surfactelement op X=0
    linesMV_XL = [] #Lijst met lijnen grondvlak tbv loadtransfer surfactelement op X=L

    exclParMembersDak = []
    exclParMembersX0 = []
    exclParMembersXL = []
    exclParMembersY0 = []
    exclParMembersYB = []
    exclParDakRandMembersX0 = []
    exclParDakRandMembersXL = []
    exclParDakRandMembersY0 = []
    exclParDakRandMembersYB = []

    #Beginnen met het invoeren van de geometrie

    for i in range(nx):
        n = i*nodes_frame
        m = i*members_frame

        #Model.clientModel.service.finish_modification(); Model.clientModel.service.begin_modification()

    #Kolommen op x-as maken
        Node(n+1, i*dx, 0.0, 0.0)
        sup_nodes.append(n+1)
        Members_in_Set = []

        if kst_kol==0:
            if i==0 or i==nx-1: EffLengthMembers[4].append(m+kst_kol+1)
            else: EffLengthMembers[1].append(m+kst_kol+1)
            DesignPropsViaParentSet = False
        else: DesignPropsViaParentSet = True

        if kst_kol >= 1:
            for k in range(kst_kol):
                Node(n+1+k+1, i*dx, 0.0, (k+1)*(h-h_dr)/(kst_kol+1))
                Member(m+k+1,n+k+1,n+1+k+1, math.radians(90) ,1,1,params={'design_properties_via_member': False, 'design_properties_via_parent_member_set': True})
                Members_in_Set.append(m+k+1)
                if i==0: MemListKol00.append(m+k+1)
                if i==nx-1: MemListKolL0.append(m+k+1)

                #Gevelligers op x-as
                if i>0:
                    EffLengthMembers[2].append(i+nx*members_frame+(ny-2+2)*(nx-1)+k*(nx-1))
                    Member(i+nx*members_frame+(ny-2+2)*(nx-1)+k*(nx-1),n+1+k+1-nodes_frame,n+1+k+1,math.radians(0) ,8,8,1,2,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False})

        Node(n+kst_kol+2, i*dx, 0.0, h-h_dr)
        Member(m+kst_kol+1,n+kst_kol+1,n+kst_kol+2,math.radians(90),1,1,params={'design_properties_via_member': not DesignPropsViaParentSet, 'design_properties_via_parent_member_set': DesignPropsViaParentSet})
        Members_in_Set.append(m+kst_kol+1)
        if i==0: MemListKol00.append(m+kst_kol+1)
        if i==nx-1: MemListKolL0.append(m+kst_kol+1)

        #Randbalk x-as maken
        if i>0:
            EffLengthMembers[2].append(i+nx*members_frame)
            Member(i+nx*members_frame,n-nodes_frame+kst_kol+2,n+kst_kol+2,math.radians(0),5,5,1,1,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False})
            MemListRandLiggerY0.append(i+nx*members_frame)

        if kst_kol >= 1:
            if i==0 or i==nx-1: EffLengthSets[4].append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER_SET))
            else: EffLengthSets[1].append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER_SET))
            MemberSet(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER_SET), insertSpaces(Members_in_Set),SetType.SET_TYPE_CONTINUOUS,"Kolom A" + str(i+1),params={'design_properties_activated': True, 'user_defined_name_enabled': True, 'name': "Kolom A" + str(i+1)})


    #Kolommen op y = b maken

        Node(n+nodes_frame, i*dx, b, 0.0)
        sup_nodes.append(n+nodes_frame)
        Members_in_Set = []

        if kst_kol==0:
            if i==0 or i==nx-1: EffLengthMembers[4].append(m+members_frame-kst_kol)
            else: EffLengthMembers[1].append(m+members_frame-kst_kol)
            DesignPropsViaParentSet = False
        else: DesignPropsViaParentSet = True

        if kst_kol >= 1:
            for k in range(kst_kol):
                Node(n+nodes_frame-k-1, i*dx, b, (k+1)*(h-h_dr)/(kst_kol+1))
                Member(m+members_frame-k,n+nodes_frame-k,n+nodes_frame-k-1, math.radians(90) ,1,1,params={'design_properties_via_member': False, 'design_properties_via_parent_member_set': True})
                Members_in_Set.append(m+members_frame-k)
                if i==0: MemListKol0B.append(m+members_frame-k)
                if i==nx-1: MemListKolLB.append(m+members_frame-k)

                if i>0:
                    EffLengthMembers[2].append(i+nx*members_frame+(ny+kst_kol+1)*(nx-1)+k*(nx-1))
                    Member(i+nx*members_frame+(ny+kst_kol+1)*(nx-1)+k*(nx-1),n-k-1,n+nodes_frame-k-1,math.radians(0) ,8,8,1,2,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False})

        Node(n+nodes_frame-kst_kol-1, i*dx, b, h-h_dr)
        Member(m+members_frame-kst_kol,n+nodes_frame-kst_kol,n+nodes_frame-kst_kol-1,math.radians(90),1,1,params={'design_properties_via_member': not DesignPropsViaParentSet, 'design_properties_via_parent_member_set': DesignPropsViaParentSet})
        Members_in_Set.append(m+members_frame-kst_kol)
        if i==0: MemListKol0B.append(m+members_frame-kst_kol)
        if i==nx-1: MemListKolLB.append(m+members_frame-kst_kol)

        #Randbalk maken
        if i>0:
            EffLengthMembers[2].append(i+nx*members_frame+(ny-1)*(nx-1))
            Member(i+nx*members_frame+(ny-1)*(nx-1),n-kst_kol-1,n+nodes_frame-kst_kol-1,math.radians(0),5,5,1,1,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False})
            MemListRandLiggerYB.append(i+nx*members_frame+(ny-1)*(nx-1))

        if kst_kol >= 1:
            if i==0 or i==nx-1: EffLengthSets[4].append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER_SET))
            else: EffLengthSets[1].append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER_SET))
            MemberSet(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER_SET), insertSpaces(Members_in_Set),SetType.SET_TYPE_CONTINUOUS,"Kolom " + KolomLetter(ny) + str(i+1),params={'design_properties_activated': True, 'user_defined_name_enabled': True, 'name': "Kolom " + KolomLetter(ny) + str(i+1)})

    #Hoofdliggers maken
        Members_in_Set = []

        if ny==2:
            EffLengthMembers[0].append(m+kst_kol+ny+1)
            DesignPropsViaParentSet = False
        else: DesignPropsViaParentSet = True

        if ny >= 3:
            Node(n+kst_kol+4, i*dx, dy, h-h_dr)
            if i==0 or i==nx-1:
                EffLengthMembers[2].append(m+kst_kol+3)
                Member(m+kst_kol+3,n+kst_kol+2,n+kst_kol+4, math.radians(0) ,4,4,1,1,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False})
            else:
                Member(m+kst_kol+3,n+kst_kol+2,n+kst_kol+4, math.radians(0) ,3,3,1,0,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False})
                Members_in_Set.append(m+kst_kol+3)
            if i==0: MemListRandLiggerX0.append(m+kst_kol+3)
            if i==nx-1: MemListRandLiggerXL.append(m+kst_kol+3)

            #Kipsteun maken
            if i>0:
                EffLengthMembers[2].append(i+nx*members_frame+(0+1)*(nx-1))
                Member(i+nx*members_frame+(0+1)*(nx-1),n+kst_kol-nodes_frame+4,n+kst_kol+4,math.radians(0),6,6,1,2)

            for d in range(1,ny-2):
                Node(n+kst_kol+d+4, i*dx, (d+1)*dy, h-h_dr)
                if i==0 or i==nx-1:
                    EffLengthMembers[2].append(m+kst_kol+d+3)
                    Member(m+kst_kol+d+3,n+kst_kol+d+3,n+kst_kol+d+4, math.radians(0) ,4,4,1,1,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False})
                else:
                    Member(m+kst_kol+d+3,n+kst_kol+d+3,n+kst_kol+d+4, math.radians(0) ,3,3)
                    Members_in_Set.append(m+kst_kol+d+3)
                if i==0: MemListRandLiggerX0.append(m+kst_kol+d+3)
                if i==nx-1: MemListRandLiggerXL.append(m+kst_kol+d+3)

                #Kipsteun maken
                if i>0:
                    EffLengthMembers[2].append(i+nx*members_frame+(d+1)*(nx-1))
                    Member(i+nx*members_frame+(d+1)*(nx-1),n+kst_kol-nodes_frame+d+4,n+kst_kol+d+4,math.radians(0),6,6,1,2)




        #Model.clientModel.service.finish_modification(); Model.clientModel.service.begin_modification()

        if ny==2:
            if i==0 or i==nx-1:
                EffLengthMembers[2].append(m+kst_kol+ny+1)
                Member(m+kst_kol+ny+1,n+kst_kol+2,n+nodes_frame-kst_kol-1, math.radians(0) ,4,4,1,1,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False})
            else:
                Member(m+kst_kol+ny+1,n+kst_kol+2,n+nodes_frame-kst_kol-1, math.radians(0) ,3,3,0,1,params={'design_properties_via_member': not DesignPropsViaParentSet, 'design_properties_via_parent_member_set': DesignPropsViaParentSet})
        else:
            if i==0 or i==nx-1:
                EffLengthMembers[2].append(m+kst_kol+ny+1)
                Member(m+kst_kol+ny+1,n+kst_kol+ny+1,n+nodes_frame-kst_kol-1, math.radians(0) ,4,4,1,1,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False})
            else:
                Member(m+kst_kol+ny+1,n+kst_kol+ny+1,n+nodes_frame-kst_kol-1, math.radians(0) ,3,3,0,1,params={'design_properties_via_member': not DesignPropsViaParentSet, 'design_properties_via_parent_member_set': DesignPropsViaParentSet})
        if i==0: MemListRandLiggerX0.append(m+kst_kol+ny+1)
        if i==nx-1: MemListRandLiggerXL.append(m+kst_kol+ny+1)

        #Model.clientModel.service.finish_modification(); Model.clientModel.service.begin_modification()

        Members_in_Set.append(m+kst_kol+ny+1)
        if ny >= 3 and i>0 and i<nx-1:
            EffLengthSets[0].append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER_SET))
            MemberSet(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER_SET), insertSpaces(Members_in_Set),SetType.SET_TYPE_CONTINUOUS,"Hoofdligger as " + str(i+1),params={'design_properties_activated': True, 'user_defined_name_enabled': True, 'name': "Hoofdligger as " + str(i+1)})

    #Dakrand verticale en horizontale profielen maken
        Node(n+kst_kol+3, i*dx, 0, h)
        EffLengthMembers[3].append(m+kst_kol+2)
        Member(m+kst_kol+2,n+kst_kol+2,n+kst_kol+3, math.radians(90) ,9,9)
        if i==0:
            MemListDakRandY0.append(m+kst_kol+2)
            MemListDakRandX0.append(m+kst_kol+2)
        if i>0:
            EffLengthMembers[2].append(i+nx*members_frame+(ny+kst_kol)*(nx-1))
            Member(i+nx*members_frame+(ny+kst_kol)*(nx-1),n+kst_kol+3-nodes_frame,n+kst_kol+3,math.radians(0),10,10,1,1)
            MemListDakRandY0.append(i+nx*members_frame+(ny+kst_kol)*(nx-1))
            if i==1: exclParDakRandMembersY0.append(i+nx*members_frame+(ny+kst_kol)*(nx-1))
        if i == nx-1:
            MemListDakRandY0.append(m+kst_kol+2)
            MemListDakRandXL.append(m+kst_kol+2)

        Node(n+nodes_frame-kst_kol-2, i*dx, b, h)
        EffLengthMembers[3].append(m+members_frame-kst_kol-1)
        Member(m+members_frame-kst_kol-1,n+nodes_frame-kst_kol-1,n+nodes_frame-kst_kol-2,math.radians(90),9,9)
        if i==0:
            MemListDakRandYB.append(m+members_frame-kst_kol-1)
            MemListDakRandX0.append(m+members_frame-kst_kol-1)
        if i>0:
            EffLengthMembers[2].append(i+nx*members_frame+(ny+2*kst_kol+1)*(nx-1))
            Member(i+nx*members_frame+(ny+2*kst_kol+1)*(nx-1),n-kst_kol-2,n+nodes_frame-kst_kol-2,math.radians(0),10,10,1,1)
            MemListDakRandYB.append(i+nx*members_frame+(ny+2*kst_kol+1)*(nx-1))
            if i==1: exclParDakRandMembersYB.append(i+nx*members_frame+(ny+2*kst_kol+1)*(nx-1))
        if i == nx-1:
            MemListDakRandYB.append(m+members_frame-kst_kol-1)
            MemListDakRandXL.append(m+members_frame-kst_kol-1)


    # Lijnen op grondvlak langsgevels tbv load transfer 2D elements
        if i>0:
            Model.clientModel.service.finish_modification(); Model.clientModel.service.begin_modification()
            lineNodes = str(n+1-nodes_frame) + " " + str(n+1)
            linesMV_Y0.append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_LINE))
            Line(linesMV_Y0[i-1],lineNodes)
            lineNodes = str(n) + " " + str(n+nodes_frame)
            linesMV_YB.append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_LINE))
            Line(linesMV_YB[i-1],lineNodes)

        Model.clientModel.service.finish_modification(); Model.clientModel.service.begin_modification()

    #Kopse gevels
    #Op x=0 goldt eerder n = 0 en m = 0
    #Op x=L goldt eerder n = (nx-1)*nodes_frame en m = (nx-1)*members_frame

    n = nx*nodes_frame
    m = nx*members_frame+(ny+2*(kst_kol+1))*(nx-1)

    if ny >= 3:
        for d in range(ny-2):

            #Kolommen as x=0 (kopse gevel)
            Members_in_Set = []
            Node(n+d*(kst_kol+2)+1, 0.0, (d+1)*dy, 0.0)
            sup_nodes.append(n+d*(kst_kol+2)+1)

            if kst_kol==0:
                EffLengthMembers[1].append(m+d*(kst_kol+2)+kst_kol+1)
                DesignPropsViaParentSet = False
            else: DesignPropsViaParentSet = True

            if kst_kol >= 1:
                for k in range(kst_kol):
                    Node(n+d*(kst_kol+2)+k+2, 0.0, (d+1)*dy, (k+1)*(h-h_dr)/(kst_kol+1))
                    Member(m+d*(kst_kol+2)+k+1,n+d*(kst_kol+2)+k+1,n+d*(kst_kol+2)+k+2, math.radians(0) ,2,2,params={'design_properties_via_member': False, 'design_properties_via_parent_member_set': True})
                    Members_in_Set.append(m+d*(kst_kol+2)+k+1)


                    #Horizontale gevelliggers as x=0:
                    EffLengthMembers[2].append(m+d+1+2*(kst_kol+2)*(ny-2)+k*(ny-1))
                    if d==0:
                        Member(m+d+1+2*(kst_kol+2)*(ny-2)+k*(ny-1),0+1+k+1,n+d*(kst_kol+2)+k+2,math.radians(0) ,8,8,1,2,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False})
                    else:
                        Member(m+d+1+2*(kst_kol+2)*(ny-2)+k*(ny-1),n+(d-1)*(kst_kol+2)+k+2,n+d*(kst_kol+2)+k+2,math.radians(0) ,8,8,1,2,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False})


            Member(m+d*(kst_kol+2)+kst_kol+1,n+d*(kst_kol+2)+kst_kol+1,0+kst_kol+d+4,math.radians(0),2,2,params={'design_properties_via_member': not DesignPropsViaParentSet, 'design_properties_via_parent_member_set': DesignPropsViaParentSet})
            Members_in_Set.append(m+d*(kst_kol+2)+kst_kol+1)
            if kst_kol >= 1:
                EffLengthSets[1].append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER_SET))
                MemberSet(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER_SET), insertSpaces(Members_in_Set),SetType.SET_TYPE_CONTINUOUS,"Kolom " + KolomLetter(d+2) + "1",params={'design_properties_activated': True, 'user_defined_name_enabled': True, 'name':"Kolom " + KolomLetter(d+2) + "1"})

            #Verticaal profiel dakrand x=0
            Node(n+d*(kst_kol+2)+kst_kol+2, 0.0, (d+1)*dy, h)
            EffLengthMembers[3].append(m+d*(kst_kol+2)+kst_kol+2)
            Member(m+d*(kst_kol+2)+kst_kol+2,0+kst_kol+d+4,n+d*(kst_kol+2)+kst_kol+2,math.radians(0),9,9,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False})

            #Horizontale profielen dakrand x=0
            EffLengthMembers[2].append(m+d+1+2*(kst_kol+2)*(ny-2)+kst_kol*(ny-1))
            if d==0:
                Member(m+d+1+2*(kst_kol+2)*(ny-2)+kst_kol*(ny-1),0+1+kst_kol+2,n+d*(kst_kol+2)+kst_kol+2,math.radians(0) ,10,10,1,1,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False})
                exclParDakRandMembersX0.append(m+d+1+2*(kst_kol+2)*(ny-2)+kst_kol*(ny-1))
            else:
                Member(m+d+1+2*(kst_kol+2)*(ny-2)+kst_kol*(ny-1),n+d*(kst_kol+2),n+d*(kst_kol+2)+kst_kol+2,math.radians(0) ,10,10,1,1,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False})
            MemListDakRandX0.insert(len(MemListDakRandX0)-1,m+d+1+2*(kst_kol+2)*(ny-2)+kst_kol*(ny-1))

            #Kolommen as x=L (kopse gevel)
            Members_in_Set = []
            Node(n+(ny-2+d)*(kst_kol+2)+1, (nx-1)*dx, (d+1)*dy, 0.0)
            sup_nodes.append(n+(ny-2+d)*(kst_kol+2)+1)

            if kst_kol==0:
                EffLengthMembers[1].append(m+(ny-2+d)*(kst_kol+2)+kst_kol+1)
                DesignPropsViaParentSet = False
            else: DesignPropsViaParentSet = True

            if kst_kol >= 1:
                for k in range(kst_kol):
                    Node(n+(ny-2+d)*(kst_kol+2)+k+2, (nx-1)*dx, (d+1)*dy, (k+1)*(h-h_dr)/(kst_kol+1))
                    Member(m+(ny-2+d)*(kst_kol+2)+k+1,n+(ny-2+d)*(kst_kol+2)+k+1,n+(ny-2+d)*(kst_kol+2)+k+2, math.radians(0) ,2,2,params={'design_properties_via_member': False, 'design_properties_via_parent_member_set': True})
                    Members_in_Set.append(m+(ny-2+d)*(kst_kol+2)+k+1)

                    #Horizontale gevelliggers as x=L:
                    EffLengthMembers[2].append(m+d+1+2*(kst_kol+2)*(ny-2)+(kst_kol+k+1)*(ny-1))
                    if d==0:
                        Member(m+d+1+2*(kst_kol+2)*(ny-2)+(kst_kol+k+1)*(ny-1),(nx-1)*nodes_frame+1+k+1,n+(ny-2+d)*(kst_kol+2)+k+2,math.radians(0) ,8,8,1,2,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False})
                    else:
                        Member(m+d+1+2*(kst_kol+2)*(ny-2)+(kst_kol+k+1)*(ny-1),n+(ny-2+d-1)*(kst_kol+2)+k+2,n+(ny-2+d)*(kst_kol+2)+k+2,math.radians(0) ,8,8,1,2,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False})

            Member(m+(ny-2+d)*(kst_kol+2)+kst_kol+1,n+(ny-2+d)*(kst_kol+2)+kst_kol+1,(nx-1)*nodes_frame+kst_kol+d+4,math.radians(0),2,2,params={'design_properties_via_member': not DesignPropsViaParentSet, 'design_properties_via_parent_member_set': DesignPropsViaParentSet})
            Members_in_Set.append(m+(ny-2+d)*(kst_kol+2)+kst_kol+1)
            if kst_kol >= 1:
                EffLengthSets[1].append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER_SET))
                MemberSet(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER_SET), insertSpaces(Members_in_Set),SetType.SET_TYPE_CONTINUOUS,"Kolom " + KolomLetter(d+2) + str(i+1),params={'design_properties_activated': True,  'user_defined_name_enabled': True, 'name':"Kolom " + KolomLetter(d+2) + str(i+1)})

            #Verticaal profiel dakrand x=L
            Node(n+(ny-2+d)*(kst_kol+2)+kst_kol+2, (nx-1)*dx, (d+1)*dy, h)
            EffLengthMembers[3].append(m+(ny-2+d)*(kst_kol+2)+kst_kol+2)
            Member(m+(ny-2+d)*(kst_kol+2)+kst_kol+2,(nx-1)*nodes_frame+kst_kol+d+4,n+(ny-2+d)*(kst_kol+2)+kst_kol+2,math.radians(0),9,9,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False})

            #Horizontale profielen dakrand x=L
            EffLengthMembers[2].append(m+d+1+2*(kst_kol+2)*(ny-2)+(2*kst_kol+1)*(ny-1))
            if d==0:
                Member(m+d+1+2*(kst_kol+2)*(ny-2)+(2*kst_kol+1)*(ny-1),(nx-1)*nodes_frame+1+kst_kol+2,n+(ny-2+d)*(kst_kol+2)+kst_kol+2,math.radians(0) ,10,10,1,1,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False})
                exclParDakRandMembersX0.append(m+d+1+2*(kst_kol+2)*(ny-2)+(2*kst_kol+1)*(ny-1))
            else:
                Member(m+d+1+2*(kst_kol+2)*(ny-2)+(2*kst_kol+1)*(ny-1),n+(ny-2+d)*(kst_kol+2),n+(ny-2+d)*(kst_kol+2)+kst_kol+2,math.radians(0) ,10,10,1,1,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False})
            MemListDakRandXL.insert(len(MemListDakRandXL)-1,m+d+1+2*(kst_kol+2)*(ny-2)+(2*kst_kol+1)*(ny-1))

            Model.clientModel.service.finish_modification(); Model.clientModel.service.begin_modification()

            # Lijnen op grondvlak kopse gevels tbv load transfer 2D elements
            if d==0:
                lineNodes = "1" + " " + str(n+d*(kst_kol+2)+1)
                linesMV_X0.append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_LINE))
                Line(linesMV_X0[d],lineNodes)
                lineNodes = str((nx-1)*nodes_frame+1) + " " + str(n+(ny-2+d)*(kst_kol+2)+1)
                linesMV_XL.append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_LINE))
                Line(linesMV_XL[d],lineNodes)
            else:
                lineNodes = str(n+(d-1)*(kst_kol+2)+1) + " " + str(n+d*(kst_kol+2)+1)
                linesMV_X0.append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_LINE))
                Line(linesMV_X0[d],lineNodes)
                lineNodes = str(n+(ny-2+d-1)*(kst_kol+2)+1) + " " + str(n+(ny-2+d)*(kst_kol+2)+1)
                linesMV_XL.append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_LINE))
                Line(linesMV_XL[d],lineNodes)

        #Horizontale gevelliggers in de kopgevels die aansluiten bij y = b:
        for k in range(kst_kol):
            EffLengthMembers[2].append(m+(ny-2)+1+2*(kst_kol+2)*(ny-2)+k*(ny-1))
            Member(m+(ny-2)+1+2*(kst_kol+2)*(ny-2)+k*(ny-1),n+((ny-2)-1)*(kst_kol+2)+k+2,nodes_frame-k-1,math.radians(0) ,8,8,1,2,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False}) #x=0
            EffLengthMembers[2].append(m+(ny-2)+1+2*(kst_kol+2)*(ny-2)+(kst_kol+k+1)*(ny-1))
            Member(m+(ny-2)+1+2*(kst_kol+2)*(ny-2)+(kst_kol+k+1)*(ny-1),n+(ny-2+(ny-2)-1)*(kst_kol+2)+k+2,(nx-1)*nodes_frame+nodes_frame-k-1,math.radians(0) ,8,8,1,2,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False}) #x=L

        #Horizontale profielen dakrand in de kopgevels die aansluiten bij y=b
        EffLengthMembers[2].append(m+(ny-2)+1+2*(kst_kol+2)*(ny-2)+kst_kol*(ny-1))
        Member(m+(ny-2)+1+2*(kst_kol+2)*(ny-2)+kst_kol*(ny-1),n+((ny-2)-1)*(kst_kol+2)+kst_kol+2,nodes_frame-kst_kol-2,math.radians(0) ,10,10,1,1,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False}) #x=0
        MemListDakRandX0.insert(len(MemListDakRandX0)-1,m+(ny-2)+1+2*(kst_kol+2)*(ny-2)+kst_kol*(ny-1))
        EffLengthMembers[2].append(m+(ny-2)+1+2*(kst_kol+2)*(ny-2)+(2*kst_kol+1)*(ny-1))
        Member(m+(ny-2)+1+2*(kst_kol+2)*(ny-2)+(2*kst_kol+1)*(ny-1),n+(ny-2+(ny-2))*(kst_kol+2),nx*nodes_frame-kst_kol-2,math.radians(0) ,10,10,1,1,params={'design_properties_via_member': True, 'design_properties_via_parent_member_set': False}) #x=L
        MemListDakRandXL.insert(len(MemListDakRandXL)-1,m+(ny-2)+1+2*(kst_kol+2)*(ny-2)+(2*kst_kol+1)*(ny-1))

        lineNodes = str(n+((ny-2)-1)*(kst_kol+2)+1) + " " + str(nodes_frame)
        linesMV_X0.append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_LINE))
        Line(linesMV_X0[d+1],lineNodes)
        lineNodes = str(n+(ny-2+(ny-2)-1)*(kst_kol+2)+1) + " " + str((nx-1)*nodes_frame+nodes_frame)
        linesMV_XL.append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_LINE))
        Line(linesMV_XL[d+1],lineNodes)

    NodalSupport(1,insertSpaces(sup_nodes), [inf, inf, inf, 0.0, 0.0, inf])

    Model.clientModel.service.finish_modification(); Model.clientModel.service.begin_modification()

    #----------------------------------------------------------------------------
    #-                      WINDVERBANDEN TOEVOEGEN                             -
    #----------------------------------------------------------------------------

    nodes = MFE_getNodes.getNodes()             # met deze lijst kan je knopen opzoeken op basis van hun coördinaten
    members = MFE_getMembers.getMembers()       # met deze lijst kan je staven opzoeken op basis van hun begin- en eindknoop
    TensionMembers = []

    #Schoren in langsgevels maken:
    xyzxyzList = []

    for k in range(kst_kol+1):
        xyzxyzList.append([1*dx,0.0,k*(h-h_dr)/(kst_kol+1),2*dx,0.0,(k+1)*(h-h_dr)/(kst_kol+1)])
        xyzxyzList.append([2*dx,0.0,k*(h-h_dr)/(kst_kol+1),1*dx,0.0,(k+1)*(h-h_dr)/(kst_kol+1)])
        xyzxyzList.append([1*dx,b,k*(h-h_dr)/(kst_kol+1),2*dx,b,(k+1)*(h-h_dr)/(kst_kol+1)])
        xyzxyzList.append([2*dx,b,k*(h-h_dr)/(kst_kol+1),1*dx,b,(k+1)*(h-h_dr)/(kst_kol+1)])
        if nx>=10:
            xyzxyzList.append([(nx-2)*dx,0.0,k*(h-h_dr)/(kst_kol+1),(nx-3)*dx,0.0,(k+1)*(h-h_dr)/(kst_kol+1)])
            xyzxyzList.append([(nx-3)*dx,0.0,k*(h-h_dr)/(kst_kol+1),(nx-2)*dx,0.0,(k+1)*(h-h_dr)/(kst_kol+1)])
            xyzxyzList.append([(nx-2)*dx,b,k*(h-h_dr)/(kst_kol+1),(nx-3)*dx,b,(k+1)*(h-h_dr)/(kst_kol+1)])
            xyzxyzList.append([(nx-3)*dx,b,k*(h-h_dr)/(kst_kol+1),(nx-2)*dx,b,(k+1)*(h-h_dr)/(kst_kol+1)])

    t = 0 #teller om ervoor te zorgen dat er maar twee schoor staven worden gebruikt als excluded members in load transfer surfaces
    exclParMembersY0.append(1+nx*members_frame)
    exclParMembersYB.append(1+nx*members_frame+(ny-1)*(nx-1))
    for m in xyzxyzList:
        if t<2:
            exclParMembersY0.append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER))
            t+=1
        elif t<4:
            exclParMembersYB.append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER))
            t+=1
        Node1 = MFE_ZoekNode.ZoekNode(m[0],m[1],m[2],Model,nodes)
        Node2 = MFE_ZoekNode.ZoekNode(m[3],m[4],m[5],Model,nodes)
        TensionMembers.append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER))
        Member.Tension(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER),Node1["no"],Node2["no"],MemberRotationSpecificationType.COORDINATE_SYSTEM_ROTATION_VIA_ANGLE, [math.radians(90)], 14)

    #Schoren in dakvlak maken:

    if ny>2:
        xyzxyzList = []

        for d in range(ny-1):
            xyzxyzList.append([1*dx,d*dy,(h-h_dr),2*dx,(d+1)*dy,(h-h_dr)])
            xyzxyzList.append([2*dx,d*dy,(h-h_dr),1*dx,(d+1)*dy,(h-h_dr)])
            if nx>=10:
                xyzxyzList.append([(nx-2)*dx,d*dy,(h-h_dr),(nx-3)*dx,(d+1)*dy,(h-h_dr)])
                xyzxyzList.append([(nx-3)*dx,d*dy,(h-h_dr),(nx-2)*dx,(d+1)*dy,(h-h_dr)])

        for m in xyzxyzList:
            Node1 = MFE_ZoekNode.ZoekNode(m[0],m[1],m[2],Model,nodes)
            Node2 = MFE_ZoekNode.ZoekNode(m[3],m[4],m[5],Model,nodes)
            TensionMembers.append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER))
            Member.Tension(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER),Node1["no"],Node2["no"],MemberRotationSpecificationType.COORDINATE_SYSTEM_ROTATION_VIA_ANGLE, [math.radians(0)], 12)

        Model.clientModel.service.finish_modification(); Model.clientModel.service.begin_modification()

    #Schoren in kopgevels maken:
    xyzxyzList = []

    for k in range(kst_kol+1):
        xyzxyzList.append([0.0,1*dy,k*(h-h_dr)/(kst_kol+1),0.0,2*dy,(k+1)*(h-h_dr)/(kst_kol+1)])
        xyzxyzList.append([0.0,2*dy,k*(h-h_dr)/(kst_kol+1),0.0,1*dy,(k+1)*(h-h_dr)/(kst_kol+1)])
        xyzxyzList.append([L,1*dy,k*(h-h_dr)/(kst_kol+1),L,2*dy,(k+1)*(h-h_dr)/(kst_kol+1)])
        xyzxyzList.append([L,2*dy,k*(h-h_dr)/(kst_kol+1),L,1*dy,(k+1)*(h-h_dr)/(kst_kol+1)])

    t = 0 #teller om ervoor te zorgen dat er maar twee schoor staven worden gebruikt als excluded members in load transfer surfaces
    exclParMembersX0.append(kst_kol+3)
    exclParMembersXL.append((nx-1)*members_frame+kst_kol+3)
    for m in xyzxyzList:
        if t<2:
            exclParMembersX0.append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER))
            t+=1
        elif t<4:
            exclParMembersXL.append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER))
            t+=1
        Node1 = MFE_ZoekNode.ZoekNode(m[0],m[1],m[2],Model,nodes)
        Node2 = MFE_ZoekNode.ZoekNode(m[3],m[4],m[5],Model,nodes)
        TensionMembers.append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER))
        Member.Tension(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER),Node1["no"],Node2["no"],MemberRotationSpecificationType.COORDINATE_SYSTEM_ROTATION_VIA_ANGLE, [math.radians(90)], 13)

    #Schoren in dakvlak (in y-richting) maken:

    if ny>2:
        xyzxyzList = []

        for d in range(nx-1):
            if d !=1 and not(nx >= 10 and d == nx-3):
                xyzxyzList.append([d*dx,1*dy,(h-h_dr),(d+1)*dx,2*dy,(h-h_dr)])
                xyzxyzList.append([d*dx,2*dy,(h-h_dr),(d+1)*dx,1*dy,(h-h_dr)])

        t = 0 #teller om ervoor te zorgen dat er maar twee schoor staven worden gebruikt als excluded members in load transfer surfaces
        exclParMembersDak.append(1+nx*members_frame+(ny-1)*(nx-1))
        for m in xyzxyzList:
            if t<2:
                exclParMembersDak.append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER))
                t+=1
            Node1 = MFE_ZoekNode.ZoekNode(m[0],m[1],m[2],Model,nodes)
            Node2 = MFE_ZoekNode.ZoekNode(m[3],m[4],m[5],Model,nodes)
            TensionMembers.append(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER))
            Member.Tension(FirstFreeIdNumber(ObjectTypes.E_OBJECT_TYPE_MEMBER),Node1["no"],Node2["no"],MemberRotationSpecificationType.COORDINATE_SYSTEM_ROTATION_VIA_ANGLE, [math.radians(0)], 11)

        Model.clientModel.service.finish_modification(); Model.clientModel.service.begin_modification()


#TODO: staafeigenschappen van schoren aanpassen (profielen aanmaken, scharnieren toevoegen, tension only)


#TODO: elif toevoegen dat er nog wel horizontale gevelliggers en dakrandligger moeten worden tegevoegd als ny=2\

    #Load transfer surfaces maken voor gevels en dak
    MemListKol0Br = MemListKol0B.copy(); MemListKol0Br.reverse()
    MemListKolLBr = MemListKolLB.copy(); MemListKolLBr.reverse()
    MemListKolL0r = MemListKolL0.copy(); MemListKolL0r.reverse()
    MemListKolLBr = MemListKolLB.copy(); MemListKolLBr.reverse()
    MemListRandLiggerX0r = MemListRandLiggerX0.copy(); MemListRandLiggerX0r.reverse()
    MemListRandLiggerXLr = MemListRandLiggerXL.copy(); MemListRandLiggerXLr.reverse()
    MemListRandLiggerY0r = MemListRandLiggerY0.copy(); MemListRandLiggerY0r.reverse()
    MemListRandLiggerYBr = MemListRandLiggerYB.copy(); MemListRandLiggerYBr.reverse()

    MemListDak = MemListRandLiggerY0 + MemListRandLiggerXL + MemListRandLiggerYBr + MemListRandLiggerX0r
    MemListX0 = MemListKol00 + MemListRandLiggerX0 + MemListKol0Br
    MemListXL = MemListKolL0 + MemListRandLiggerXL + MemListKolLBr
    MemListY0 = MemListKol00 + MemListRandLiggerY0 + MemListKolL0r
    MemListYB = MemListKol0B + MemListRandLiggerYB + MemListKolLBr

    MemListDakRandX0 += MemListRandLiggerX0r
    MemListDakRandXL += MemListRandLiggerXLr
    MemListDakRandY0 += MemListRandLiggerY0r
    MemListDakRandYB += MemListRandLiggerYBr

    LineListDak = MemberToLineList(MemListDak,members)

    linesMV_X0r = linesMV_X0.copy(); linesMV_X0r.reverse()
    linesMV_XLr = linesMV_XL.copy(); linesMV_XLr.reverse()
    linesMV_Y0r = linesMV_Y0.copy(); linesMV_Y0r.reverse()
    linesMV_YBr = linesMV_YB.copy(); linesMV_YBr.reverse()

    LineListX0 = MemberToLineList(MemListX0,members) + linesMV_X0r
    LineListXL = MemberToLineList(MemListXL,members) + linesMV_XLr
    LineListY0 = MemberToLineList(MemListY0,members) + linesMV_Y0r
    LineListYB = MemberToLineList(MemListYB,members) + linesMV_YBr
    LineListDakRandX0 = MemberToLineList(MemListDakRandX0,members)
    LineListDakRandXL = MemberToLineList(MemListDakRandXL,members)
    LineListDakRandY0 = MemberToLineList(MemListDakRandY0,members)
    LineListDakRandYB = MemberToLineList(MemListDakRandYB,members)

    LineListXL.reverse()
    LineListY0.reverse()
    LineListDakRandXL.reverse()
    LineListDakRandY0.reverse()

    Surface(1,insertSpaces(LineListDak),params={'geometry':'GEOMETRY_PLANE','type':'TYPE_LOAD_TRANSFER','load_transfer_direction':"LOAD_TRANSFER_DIRECTION_IN_X",'load_distribution':"LOAD_DISTRIBUTION_VARYING",'excluded_parallel_to_members':insertSpaces(exclParMembersDak)})
    Surface(2,insertSpaces(LineListX0),params={'geometry':'GEOMETRY_PLANE','type':'TYPE_LOAD_TRANSFER','load_transfer_direction':"LOAD_TRANSFER_DIRECTION_IN_X",'load_distribution':"LOAD_DISTRIBUTION_VARYING",'excluded_parallel_to_members':insertSpaces(exclParMembersX0)})
    Surface(3,insertSpaces(LineListXL),params={'geometry':'GEOMETRY_PLANE','type':'TYPE_LOAD_TRANSFER','load_transfer_direction':"LOAD_TRANSFER_DIRECTION_IN_X",'load_distribution':"LOAD_DISTRIBUTION_VARYING",'excluded_parallel_to_members':insertSpaces(exclParMembersXL)})
    Surface(4,insertSpaces(LineListY0),params={'geometry':'GEOMETRY_PLANE','type':'TYPE_LOAD_TRANSFER','load_transfer_direction':"LOAD_TRANSFER_DIRECTION_IN_X",'load_distribution':"LOAD_DISTRIBUTION_VARYING",'excluded_parallel_to_members':insertSpaces(exclParMembersY0)})
    Surface(5,insertSpaces(LineListYB),params={'geometry':'GEOMETRY_PLANE','type':'TYPE_LOAD_TRANSFER','load_transfer_direction':"LOAD_TRANSFER_DIRECTION_IN_X",'load_distribution':"LOAD_DISTRIBUTION_VARYING",'excluded_parallel_to_members':insertSpaces(exclParMembersYB)})
    Surface(6,insertSpaces(LineListDakRandX0),params={'geometry':'GEOMETRY_PLANE','type':'TYPE_LOAD_TRANSFER','load_transfer_direction':"LOAD_TRANSFER_DIRECTION_IN_X",'load_distribution':"LOAD_DISTRIBUTION_VARYING",'excluded_parallel_to_members':insertSpaces(exclParDakRandMembersX0)})
    Surface(7,insertSpaces(LineListDakRandXL),params={'geometry':'GEOMETRY_PLANE','type':'TYPE_LOAD_TRANSFER','load_transfer_direction':"LOAD_TRANSFER_DIRECTION_IN_X",'load_distribution':"LOAD_DISTRIBUTION_VARYING",'excluded_parallel_to_members':insertSpaces(exclParDakRandMembersXL)})
    Surface(8,insertSpaces(LineListDakRandY0),params={'geometry':'GEOMETRY_PLANE','type':'TYPE_LOAD_TRANSFER','load_transfer_direction':"LOAD_TRANSFER_DIRECTION_IN_X",'load_distribution':"LOAD_DISTRIBUTION_VARYING",'excluded_parallel_to_members':insertSpaces(exclParDakRandMembersY0)})
    Surface(9,insertSpaces(LineListDakRandYB),params={'geometry':'GEOMETRY_PLANE','type':'TYPE_LOAD_TRANSFER','load_transfer_direction':"LOAD_TRANSFER_DIRECTION_IN_X",'load_distribution':"LOAD_DISTRIBUTION_VARYING",'excluded_parallel_to_members':insertSpaces(exclParDakRandMembersYB)})

#----------------------------------------------------------------------------------------------
#
#                       INSTELLINGEN EC3 TOETSINGEN
#
#----------------------------------------------------------------------------------------------
    sys.stdout.write("\r                                                                       ")
    sys.stdout.write("\rINSTELLINGEN EC3 TOETSINGEN\n")
    sys.stdout.flush()

    #Instellingen Steel Effective Lengths voor hoofdliggers maken:
    NodalSupportsList = []
    # factors = []
    TussenPunten = 5
    # for f in range(TussenPunten):
    #     factors.append([1,1,1,1,1,1,1,1,1,1])
    NodalSupportsList.append([SteelEffectiveLengthsSupportType.SUPPORT_TYPE_FIXED_IN_Z_Y_AND_TORSION, True, 0.0, \
                                SteelEffectiveLengthsEccentricityType.ECCENTRICITY_TYPE_NONE, 0,0,0,0, \
                                SteelEffectiveLengthsSupportTypeInY.SUPPORT_STATUS_YES, SteelEffectiveLengthsRestraintTypeAboutX.SUPPORT_STATUS_YES, \
                                SteelEffectiveLengthsRestraintTypeAboutZ.SUPPORT_STATUS_NO, SteelEffectiveLengthsRestraintTypeWarping.SUPPORT_STATUS_NO, \
                                ""])
    for i in range(TussenPunten): #TODO: Range afhankelijk maken van ny voor het geval de hoofdliggers meer kipsteunen hebben.
        NodalSupportsList.append([SteelEffectiveLengthsSupportType.SUPPORT_TYPE_FIXED_IN_Y, False, 0.0, \
                                    SteelEffectiveLengthsEccentricityType.ECCENTRICITY_TYPE_AT_UPPER_FLANGE, 1,4,2,3, \
                                    SteelEffectiveLengthsSupportTypeInY.SUPPORT_STATUS_YES, SteelEffectiveLengthsRestraintTypeAboutX.SUPPORT_STATUS_NO, \
                                    SteelEffectiveLengthsRestraintTypeAboutZ.SUPPORT_STATUS_NO, SteelEffectiveLengthsRestraintTypeWarping.SUPPORT_STATUS_NO, \
                                    ""],)

    NodalSupportsList.append([SteelEffectiveLengthsSupportType.SUPPORT_TYPE_FIXED_IN_Z_Y_AND_TORSION, True, 0.0, \
                                SteelEffectiveLengthsEccentricityType.ECCENTRICITY_TYPE_NONE, 0,0,0,0, \
                                SteelEffectiveLengthsSupportTypeInY.SUPPORT_STATUS_YES, SteelEffectiveLengthsRestraintTypeAboutX.SUPPORT_STATUS_YES, \
                                SteelEffectiveLengthsRestraintTypeAboutZ.SUPPORT_STATUS_NO, SteelEffectiveLengthsRestraintTypeWarping.SUPPORT_STATUS_NO, \
                                ""])

    SteelEffectiveLengths(1, insertSpaces(EffLengthMembers[0]), insertSpaces(EffLengthSets[0]), False, False, False, True, True, True, 'SEL1', NodalSupportsList, intermediate_nodes=True,
                          different_properties=True, determination_of_mcr=SteelEffectiveLengthsDeterminationMcrEurope.DETERMINATION_EUROPE_EIGENVALUE)

    #Instellingen Steel Effective Lengths voor kolommen maken:
    NodalSupportsList = []
    NodalSupportsList.append([SteelEffectiveLengthsSupportType.SUPPORT_TYPE_FIXED_IN_Z_Y_AND_TORSION, True, 0.0, \
                                SteelEffectiveLengthsEccentricityType.ECCENTRICITY_TYPE_NONE, 0,0,0,0, \
                                SteelEffectiveLengthsSupportTypeInY.SUPPORT_STATUS_YES, SteelEffectiveLengthsRestraintTypeAboutX.SUPPORT_STATUS_YES, \
                                SteelEffectiveLengthsRestraintTypeAboutZ.SUPPORT_STATUS_NO, SteelEffectiveLengthsRestraintTypeWarping.SUPPORT_STATUS_NO, \
                                ""])
    for i in range(5): #TODO: Range afhankelijk maken van ny voor het geval de kolommen meer kniksteunen hebben.
        NodalSupportsList.append([SteelEffectiveLengthsSupportType.SUPPORT_TYPE_FIXED_IN_Y, False, 0.0, \
                                    SteelEffectiveLengthsEccentricityType.ECCENTRICITY_TYPE_NONE, 0,0,0,0, \
                                    SteelEffectiveLengthsSupportTypeInY.SUPPORT_STATUS_YES, SteelEffectiveLengthsRestraintTypeAboutX.SUPPORT_STATUS_NO, \
                                    SteelEffectiveLengthsRestraintTypeAboutZ.SUPPORT_STATUS_NO, SteelEffectiveLengthsRestraintTypeWarping.SUPPORT_STATUS_NO, \
                                    ""],)

    NodalSupportsList.append([SteelEffectiveLengthsSupportType.SUPPORT_TYPE_FIXED_IN_Z_Y_AND_TORSION, True, 0.0, \
                                SteelEffectiveLengthsEccentricityType.ECCENTRICITY_TYPE_NONE, 0,0,0,0, \
                                SteelEffectiveLengthsSupportTypeInY.SUPPORT_STATUS_YES, SteelEffectiveLengthsRestraintTypeAboutX.SUPPORT_STATUS_YES, \
                                SteelEffectiveLengthsRestraintTypeAboutZ.SUPPORT_STATUS_NO, SteelEffectiveLengthsRestraintTypeWarping.SUPPORT_STATUS_NO, \
                                ""])

    SteelEffectiveLengths(2, insertSpaces(EffLengthMembers[1]), insertSpaces(EffLengthSets[1]), True, True, False, True, True, True, 'SEL2', NodalSupportsList, intermediate_nodes=True,
                          different_properties=True, determination_of_mcr=SteelEffectiveLengthsDeterminationMcrEurope.DETERMINATION_EUROPE_EIGENVALUE)

    #Instellingen Steel Effective Lengths voor andere liggers maken:
    NodalSupportsList = []
    NodalSupportsList.append([SteelEffectiveLengthsSupportType.SUPPORT_TYPE_FIXED_IN_Z_Y_AND_TORSION, True, 0.0, \
                                SteelEffectiveLengthsEccentricityType.ECCENTRICITY_TYPE_NONE, 0,0,0,0, \
                                SteelEffectiveLengthsSupportTypeInY.SUPPORT_STATUS_YES, SteelEffectiveLengthsRestraintTypeAboutX.SUPPORT_STATUS_YES, \
                                SteelEffectiveLengthsRestraintTypeAboutZ.SUPPORT_STATUS_NO, SteelEffectiveLengthsRestraintTypeWarping.SUPPORT_STATUS_NO, \
                                ""])

    NodalSupportsList.append([SteelEffectiveLengthsSupportType.SUPPORT_TYPE_FIXED_IN_Z_Y_AND_TORSION, True, 0.0, \
                                SteelEffectiveLengthsEccentricityType.ECCENTRICITY_TYPE_NONE, 0,0,0,0, \
                                SteelEffectiveLengthsSupportTypeInY.SUPPORT_STATUS_YES, SteelEffectiveLengthsRestraintTypeAboutX.SUPPORT_STATUS_YES, \
                                SteelEffectiveLengthsRestraintTypeAboutZ.SUPPORT_STATUS_NO, SteelEffectiveLengthsRestraintTypeWarping.SUPPORT_STATUS_NO, \
                                ""])

    SteelEffectiveLengths(3, insertSpaces(EffLengthMembers[2]), insertSpaces(EffLengthSets[2]), False, False, False, True, True, True, 'SEL3', NodalSupportsList, intermediate_nodes=False,
                          different_properties=True, determination_of_mcr=SteelEffectiveLengthsDeterminationMcrEurope.DETERMINATION_EUROPE_EIGENVALUE)

    #Instellingen Steel Effective Lengths voor verticale dakrandprofielen maken:
    NodalSupportsList = []
    NodalSupportsList.append([SteelEffectiveLengthsSupportType.SUPPORT_TYPE_FIXED_IN_Z_Y_AND_TORSION, True, 0.0, \
                                SteelEffectiveLengthsEccentricityType.ECCENTRICITY_TYPE_NONE, 0,0,0,0, \
                                SteelEffectiveLengthsSupportTypeInY.SUPPORT_STATUS_YES, SteelEffectiveLengthsRestraintTypeAboutX.SUPPORT_STATUS_YES, \
                                SteelEffectiveLengthsRestraintTypeAboutZ.SUPPORT_STATUS_YES, SteelEffectiveLengthsRestraintTypeWarping.SUPPORT_STATUS_YES, \
                                ""])

    NodalSupportsList.append([SteelEffectiveLengthsSupportType.SUPPORT_TYPE_NONE, True, 0.0, \
                                SteelEffectiveLengthsEccentricityType.ECCENTRICITY_TYPE_NONE, 0,0,0,0, \
                                SteelEffectiveLengthsSupportTypeInY.SUPPORT_STATUS_NO, SteelEffectiveLengthsRestraintTypeAboutX.SUPPORT_STATUS_NO, \
                                SteelEffectiveLengthsRestraintTypeAboutZ.SUPPORT_STATUS_NO, SteelEffectiveLengthsRestraintTypeWarping.SUPPORT_STATUS_NO, \
                                ""])

    SteelEffectiveLengths(4, insertSpaces(EffLengthMembers[3]), insertSpaces(EffLengthSets[3]), False, False, False, True, True, True, 'SEL4', NodalSupportsList, intermediate_nodes=False,
                          different_properties=True, determination_of_mcr=SteelEffectiveLengthsDeterminationMcrEurope.DETERMINATION_EUROPE_EIGENVALUE)

    #Instellingen Steel Effective Lengths voor hoekkolommen maken:
    NodalSupportsList = []
    NodalSupportsList.append([SteelEffectiveLengthsSupportType.SUPPORT_TYPE_FIXED_IN_Z_Y_AND_TORSION, True, 0.0, \
                                SteelEffectiveLengthsEccentricityType.ECCENTRICITY_TYPE_NONE, 0,0,0,0, \
                                SteelEffectiveLengthsSupportTypeInY.SUPPORT_STATUS_YES, SteelEffectiveLengthsRestraintTypeAboutX.SUPPORT_STATUS_YES, \
                                SteelEffectiveLengthsRestraintTypeAboutZ.SUPPORT_STATUS_NO, SteelEffectiveLengthsRestraintTypeWarping.SUPPORT_STATUS_NO, \
                                ""])
    for i in range(5): #TODO: Range afhankelijk maken van ny voor het geval de kolommen meer kniksteunen hebben.
        NodalSupportsList.append([SteelEffectiveLengthsSupportType.SUPPORT_TYPE_INDIVIDUALLY, True, 0.0, \
                                    SteelEffectiveLengthsEccentricityType.ECCENTRICITY_TYPE_NONE, 0,0,0,0, \
                                    SteelEffectiveLengthsSupportTypeInY.SUPPORT_STATUS_YES, SteelEffectiveLengthsRestraintTypeAboutX.SUPPORT_STATUS_NO, \
                                    SteelEffectiveLengthsRestraintTypeAboutZ.SUPPORT_STATUS_NO, SteelEffectiveLengthsRestraintTypeWarping.SUPPORT_STATUS_NO, \
                                    ""])

    NodalSupportsList.append([SteelEffectiveLengthsSupportType.SUPPORT_TYPE_FIXED_IN_Z_Y_AND_TORSION, True, 0.0, \
                                SteelEffectiveLengthsEccentricityType.ECCENTRICITY_TYPE_NONE, 0,0,0,0, \
                                SteelEffectiveLengthsSupportTypeInY.SUPPORT_STATUS_YES, SteelEffectiveLengthsRestraintTypeAboutX.SUPPORT_STATUS_YES, \
                                SteelEffectiveLengthsRestraintTypeAboutZ.SUPPORT_STATUS_NO, SteelEffectiveLengthsRestraintTypeWarping.SUPPORT_STATUS_NO, \
                                ""])

    SteelEffectiveLengths(5, insertSpaces(EffLengthMembers[4]), insertSpaces(EffLengthSets[4]), True, True, False, True, True, True, 'SEL5', NodalSupportsList, intermediate_nodes=True,
                          different_properties=True, determination_of_mcr=SteelEffectiveLengthsDeterminationMcrEurope.DETERMINATION_EUROPE_EIGENVALUE)


    SteelDesignUltimateConfigurations(1, name="EC3 checks UGT")
    SteelDesignServiceabilityConfigurations(1,"EC3 checks BGT")
    #BeamLimitCharacteristic(134,1)

#----------------------------------------------------------------------------------------------
#
#                       STATIC ANALYSIS SETTINGS
#
#----------------------------------------------------------------------------------------------

    StaticAnalysisSettings.GeometricallyLinear(1, "Linear")
    StaticAnalysisSettings.SecondOrderPDelta(2, "SecondOrder")
    StaticAnalysisSettings.LargeDeformation(3, "LargeDeformation")

#----------------------------------------------------------------------------------------------
#
#                       BELASTINGEN
#
#----------------------------------------------------------------------------------------------
    sys.stdout.write("\r                                                                       ")
    sys.stdout.write("\rBELASTINGEN\n")
    sys.stdout.flush()

    LoadCasesAndCombinations({
                    "current_standard_for_combination_wizard": 6047,
                    "activate_combination_wizard_and_classification": True,
                    "activate_combination_wizard": True,
                    "result_combinations_active": False,
                    "result_combinations_parentheses_active": False,
                    "result_combinations_consider_sub_results": False,
                    "combination_name_according_to_action_category": False},
                 model= Model)

    LoadCase.StaticAnalysis(10, 'PB: Eigen Gewicht',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_PERMANENT_G, self_weight=[True, 0.0, 0.0, -1.0])
    LoadCase.StaticAnalysis(20, 'PB: Vloerafwerking en dakbedekking',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_PERMANENT_G, self_weight=[False])
    LoadCase.StaticAnalysis(30, 'PB: Gevelbekleding',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_PERMANENT_G, self_weight=[False])
    # LoadCase.StaticAnalysis(40, 'PB: Trappen, leuningen en bordessen',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_PERMANENT_G, self_weight=[False])
    # LoadCase.StaticAnalysis(50, 'PB: Klein installatie werk (gelijkmatig verdeeld)',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_PERMANENT_IMPOSED_GQ, self_weight=[False])
    # LoadCase.StaticAnalysis(60, 'PB: Grote leidingen en leidingsupports',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_PERMANENT_IMPOSED_GQ, self_weight=[False])
    # LoadCase.StaticAnalysis(70, 'PB: Overige supports',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_PERMANENT_IMPOSED_GQ, self_weight=[False])
    # LoadCase.StaticAnalysis(80, 'PB: Gronddrukken',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_PERMANENT_G, self_weight=[False])
    # LoadCase.StaticAnalysis(90, 'PB: Reserve',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_PERMANENT_G, self_weight=[False])
    # LoadCase.StaticAnalysis(100, 'VB: Opgelegde belasting op vloeren (momentaan)',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_IMPOSED_LOADS_CATEGORY_E_STORAGE_AREAS_QI_E, self_weight=[False])
    LoadCase.StaticAnalysis(300, 'VB: Sneeuw en/of aanvriezing',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_SNOW_ICE_LOADS_H_LESS_OR_EQUAL_TO_1000_M_QS, self_weight=[False])
    LoadCase.StaticAnalysis(401, 'VB: Wind (+X richting; overdruk)',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_WIND_QW, self_weight=[False])
    LoadCase.StaticAnalysis(402, 'VB: Wind (-X richting; overdruk)',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_WIND_QW, self_weight=[False])
    LoadCase.StaticAnalysis(403, 'VB: Wind (+Y richting; overdruk)',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_WIND_QW, self_weight=[False])
    LoadCase.StaticAnalysis(404, 'VB: Wind (-Y richting; overdruk)',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_WIND_QW, self_weight=[False])
    LoadCase.StaticAnalysis(411, 'VB: Wind (+X richting; onderdruk)',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_WIND_QW, self_weight=[False])
    LoadCase.StaticAnalysis(412, 'VB: Wind (-X richting; onderdruk)',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_WIND_QW, self_weight=[False])
    LoadCase.StaticAnalysis(413, 'VB: Wind (+Y richting; onderdruk)',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_WIND_QW, self_weight=[False])
    LoadCase.StaticAnalysis(414, 'VB: Wind (-Y richting; onderdruk)',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_WIND_QW, self_weight=[False])
    # LoadCase.StaticAnalysis(510, 'VB: Thermische belasting (geval 1)',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_TEMPERATURE_NON_FIRE_QT, self_weight=[False])
    # LoadCase.StaticAnalysis(520, 'VB: Thermische belasting (geval 2)',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_TEMPERATURE_NON_FIRE_QT, self_weight=[False])
    # LoadCase.StaticAnalysis(611, 'Hijslasten (geval 1 met hor. kracht in +X richting)',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_IMPOSED_LOADS_CATEGORY_E_STORAGE_AREAS_QI_E, self_weight=[False])
    # LoadCase.StaticAnalysis(612, 'Hijslasten (geval 1 met hor. kracht in -X richting)',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_IMPOSED_LOADS_CATEGORY_E_STORAGE_AREAS_QI_E, self_weight=[False])
    # LoadCase.StaticAnalysis(613, 'Hijslasten (geval 1 met hor. kracht in +Y richting)',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_IMPOSED_LOADS_CATEGORY_E_STORAGE_AREAS_QI_E, self_weight=[False])
    # LoadCase.StaticAnalysis(614, 'Hijslasten (geval 1 met hor. kracht in -Y richting)',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_IMPOSED_LOADS_CATEGORY_E_STORAGE_AREAS_QI_E, self_weight=[False])
    # LoadCase.StaticAnalysis(710, 'BB: Stootbelastingen of ontploffingen (geval 1)',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_ACCIDENTAL_ACTIONS_A, self_weight=[False])
    # LoadCase.StaticAnalysis(720, 'BB: Stootbelastingen of ontploffingen (geval 2)',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_ACCIDENTAL_ACTIONS_A, self_weight=[False])
    # LoadCase.StaticAnalysis(900, 'PB: Equipement - lege gewicht',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_PERMANENT_G, self_weight=[False])
    # LoadCase.StaticAnalysis(910, 'PB: Equipement - operationele gewicht',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_PERMANENT_IMPOSED_GQ, self_weight=[False])
    # LoadCase.StaticAnalysis(920, 'PB: Equipement - test gewicht',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_PERMANENT_IMPOSED_GQ, self_weight=[False])
    LoadCase.StaticAnalysis(990, 'Voorspanning trekschoren',analysis_settings_no=1,action_category=ActionCategoryType.ACTION_CATEGORY_PRESTRESS_P, self_weight=[False])

    Sneeuw(L,b,dhr,1)

    #SnowWizardMonopitch([4,9,57,52])

#----------------------------------------------------------------------------------------------
#
#                       WINDBELASTING
#
#----------------------------------------------------------------------------------------------

    Wind(qpz,L,b,h,h_dr,dhr,1)

# COMBINATIONWIZARD:
w = CombinationWizard(1,consider_imperfection_case=False, params={'consider_initial_state': False,'structure_modification_enabled': False})
DesignSituation(1,DesignSituationType.DESIGN_SITUATION_TYPE_STR_PERMANENT_AND_TRANSIENT_6_10A_6_10B,True,'ULS (STR/GEO) - Permanent and transient - Eq. 6.10a and 6.10b',params={'combination_wizard': 1})
DesignSituation(2,DesignSituationType.DESIGN_SITUATION_TYPE_SLS_CHARACTERISTIC,True,'SLS - Characteristic',params={'combination_wizard': 1})


#VOORSPANNING

MemberLoad.InitialPrestress(1,990,insertSpaces(TensionMembers),MemberLoadDirection.LOAD_DIRECTION_LOCAL_X,100)




Model.clientModel.service.finish_modification()
Model.clientModel.service.generate_load_cases_and_combinations()
Calculate_all()


MIF = ResultTables.MembersInternalForces(CaseObjectType.E_OBJECT_TYPE_DESIGN_SITUATION,1,17)
MIFS = ResultTables.MembersInternalForcesBySection(CaseObjectType.E_OBJECT_TYPE_DESIGN_SITUATION,1,2,False)

print("test")
print("T1 = " + str(time.time()-time1) + "s")