beamY_rg=gm.IJKRange((0,0,0),(0,lastYpos,0)) beamY=gridGeom.genLinOneRegion(ijkRange=beamY_rg,setName='beamY') # *** MATERIALS *** S235JR= EC3_materials.S235JR S235JR.gammaM= 1.00 # Steel material-section appropiate for 3D beam analysis, including shear # deformations. # Attributes: # steel: steel material ( # name: name of the standard steel profile. Types: IPEShape, HEShape, # UPNShape, AUShape, CHSShape # (defined in materials.sections.structural_shapes.arcelor_metric_shapes) beamY_mat= EC3_materials.IPEShape(steel=S235JR,name='IPE_A_450') beamY_mat.defElasticShearSection3d(preprocessor) # ***FE model - MESH*** beamY_mesh=fem.LinSetToMesh(linSet=beamY,matSect=beamY_mat,elemSize=eSize,vDirLAxZ=xc.Vector([1,0,0]),elemType='ElasticBeam3d',coordTransfType='linear') beamY_mesh.generateMesh(preprocessor) # EC3beam definition lstLines=gridGeom.getLstLinRange(beamY_rg) from materials.ec3 import EC3Beam as ec3b ec3beam=ec3b.EC3Beam(name='ec3bm',ec3Shape=None,lstLines=lstLines) ec3beam.setControlPoints()
from __future__ import division __author__= "Luis C. Pérez Tato (LCPT)" __copyright__= "Copyright 2014, LCPT" __license__= "GPL" __version__= "3.0" __email__= "*****@*****.**" from geom_utils import interpolation as intp from materials.ec3 import EC3_limit_state_checking as EC3lsc from materials.ec3 import EC3_materials S355JR= EC3_materials.S355JR gammaM0= 1.05 S355JR.gammaM= gammaM0 IPE400= EC3_materials.IPEShape(S355JR,"IPE_400") # Geometry k1= 1.0; k2= 1.0 #Check results page 32 L= 6.0 # Bar length (m) x= [0.0,0.25*L,0.5*L,0.75*L,1.0*L] M= [-93.7,0,114.3,0,111.4] Mi=intp.interpEquidistPoints(xi=x,yi=M,nDiv=4) mgf= EC3lsc.MomentGradientFactorC1(Mi) Mcr1= IPE400.getMcr(L,Mi) Mcr1Teor= 164.7e3 ratio1= abs(Mcr1-Mcr1Teor)/Mcr1Teor #NOTE: Here there is a big difference between the results
beamY_mat.setupElasticShear3DSection(preprocessor=prep)
columnZconcr_mat = tm.BeamMaterialData(name='columnZconcr_mat',
section=geomSectColumnZ,
material=concrProp)
columnZconcr_mat.setupElasticShear3DSection(preprocessor=prep)
# Steel material-section appropiate for 3D beam analysis, including shear
# deformations.
# Attributes:
# steel: steel material (
# name: name of the standard steel profile. Types: IPEShape, HEShape,
# UPNShape, AUShape, CHSShape
# (defined in materials.sections.structural_shapes.arcelor_metric_shapes)
columnZsteel_mat = EC3_materials.HEShape(steel=S235JR, name='HE_200_A')
columnZsteel_mat.defElasticShearSection3d(prep)
beamXsteel_mat = EC3_materials.IPEShape(steel=S235JR, name='IPE_A_300')
beamXsteel_mat.defElasticShearSection3d(prep)
# ***FE model - MESH***
# IMPORTANT: it's convenient to generate the mesh of surfaces before meshing
# the lines, otherwise, sets of shells can take also beam elements touched by
# them
beamXconcr_mesh = fem.LinSetToMesh(linSet=beamXconcr,
matSect=beamXconcr_mat,
elemSize=eSize,
vDirLAxZ=xc.Vector([0, 1, 0]),
elemType='ElasticBeam3d',
dimElemSpace=3,
coordTransfType='linear')
beamY_mesh = fem.LinSetToMesh(linSet=beamY,
#Coefficients alpha, beta for biaxial bending (clause 6.2.9 of EC3.1.1)
__author__ = "Ana Ortega (AO_O)"
__copyright__ = "Copyright 2018, AO_O"
__license__ = "GPL"
__version__ = "3.0"
__email__ = "*****@*****.**"
from materials.ec3 import EC3_materials
S235JR = EC3_materials.S235JR
S235JR.gammaM = 1.00
sctClass = 1 #section class
Nd = 7e5 #axial force [N]
IPE300 = EC3_materials.IPEShape(steel=S235JR, name='IPE_A_300')
alpha_IPE300, beta_IPE300 = IPE300.getBiaxBendCoeffs(Nd,
IPE300.getNcRd(sctClass))
alpha_IPE300_comp = 2
beta_IPE300_comp = 5 * Nd / 1093455.0
IPN300 = EC3_materials.IPNShape(steel=S235JR, name='IPN_300')
alpha_IPN300, beta_IPN300 = IPN300.getBiaxBendCoeffs(Nd,
IPN300.getNcRd(sctClass))
alpha_IPN300_comp = 2
beta_IPN300_comp = 5 * Nd / 1621500.0
HE300 = EC3_materials.HEShape(steel=S235JR, name='HE_300_A')
alpha_HE300, beta_HE300 = HE300.getBiaxBendCoeffs(Nd, HE300.getNcRd(sctClass))
alpha_HE300_comp = 2
beta_HE300_comp = 5 * Nd / 2643750.0
#Lateral-torsional buckling (see Aitziber López, Danny J. Yong and Miguel A. Serna article)
from __future__ import division
import math
import xc_base
import geom
import xc
import scipy.interpolate
from materials.ec3 import EC3_materials
from rough_calculations import ng_simple_beam as sb
S235JR= EC3_materials.S235JR
S235JR.gammaM= 1.00
IPE450A= EC3_materials.IPEShape(S235JR,'IPE_A_450')
# Geometry
# k1=lateral bending and warping coefficient at first end (free:1, prevented:0.5)
# k2=lateral bending and warping coefficient at last end (free:1, prevented:0.5)
k1= 1.0; k2= 1.0
#Check results pages 34 and 35
L= 3.2 # Bar length (m)
x= [0.0,L]
M= [0.0,-1296e3] #values of the moment at sections in x abcissae, each of them
#with the corresponding sign.
overlineLambdaLT= IPE450A.getLateralBucklingNonDimensionalBeamSlenderness(sectionClass=1,xi=x,Mi=M) #xi: abcissae for the moment diagram, ordinates for the moment diagram
alphaLT= IPE450A.getLateralBucklingImperfectionFactor()
# phiLT= IPE450A.getLateralBucklingIntermediateFactor(1,x,M)