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
