import xc_base import geom import xc from materials.sections.structural_shapes import arcelor_metric_shapes from materials.ec3 import EC3_materials from model import predefined_spaces from solution import predefined_solutions test = xc.FEProblem() preprocessor = test.getPreprocessor S355JR = EC3_materials.S355JR S355JR.gammaM = 1.05 epp = S355JR.getDesignElasticPerfectlyPlasticMaterial(preprocessor, "epp") IPE200 = arcelor_metric_shapes.IPEShape(S355JR, 'IPE_200') fs3d = IPE200.getFiberSection3d(preprocessor, 'epp') L = 1.0 nodes = preprocessor.getNodeHandler modelSpace = predefined_spaces.StructuralMechanics2D(nodes) nodes.defaultTag = 10 # First node number. nod1 = nodes.newNodeXY(0, 0.0) nod2 = nodes.newNodeXY(L, 0.0) # Geometric transformations lin = modelSpace.newLinearCrdTransf("lin") # Elements definition elements = preprocessor.getElementHandler elements.defaultTransformation = lin.name # Coordinate transformation for the new elements elements.defaultMaterial = IPE200.fiberSection3dName
# -*- coding: utf-8 -*- import math from rough_calculations import ng_simple_beam as sb from rough_calculations import ng_pinned_fixed_beam as pfb from materials.sections.structural_shapes import arcelor_metric_shapes from materials.ec3 import EC3_materials S355JR = EC3_materials.S355JR S355JR.gammaM = 1.05 IPE450A = arcelor_metric_shapes.IPEShape(S355JR, 'IPE_A_450') #beam= sb.SimpleBeam(205e9,IPE450A.Iz()) beam = pfb.PinnedFixedBeam(205e9, IPE450A.Iz()) beam.l = math.sqrt(15.5**2 + 1**2) b = 15.5 / 5.0 x3 = beam.l / 2.0 x2 = x3 - b x1 = x2 - b x4 = x3 + b x5 = x4 + b pos = [x1, x2, x3, x4, x5] w = IPE450A.getRho() * 9.81 P0 = 14.5063917522e3 defl0 = 0.0 P1 = 19.9114175252e3 defl1 = 0.0