## Point load 3.
cLC = loadCaseManager.setCurrentLoadCase('p3')
pLoad = xc.Vector([4.45e3, -890e3, 0.0])
p1.getNode().newLoad(pLoad)

## Control sections
topNode = p1.getNode()
bottomNode = p0.getNode()
bottomElement = preprocessor.getSets.defSet('bottomElement')
bottomElement.nodes.append(bottomNode)
steelColumn.errFileName = "/tmp/erase.log"  # Don't print error messages.
bottomElement.fillUpwards()
steelColumn.errFileName = "cerr"  # Display errors if any.

## Load combinations
combContainer = combs.CombContainer()

combContainer.SLS.qp.add('c0', '1.0*p0')
combContainer.SLS.qp.add('c1', '1.0*p1')
combContainer.SLS.qp.add('c2', '1.0*p2')
combContainer.SLS.qp.add('c3', '1.0*p3')

### Dump combination definition into XC.
combContainer.dumpCombinations(preprocessor)

# Reference values
refMom = [38.0e3, 53.1e3, 67.7e3, 96.2e3]
refDisp = [22.9e-3, 33.9e-3, 44.6e-3, 65.4e-3]

# Solution
analysis = predefined_solutions.penalty_newton_raphson(steelColumn)
Ejemplo n.º 2
0
# -*- coding: utf-8 -*-
from actions import combinations as cc
#Combinaciones determinantes en el diseño de armaduras (para tanteo rápido)

combContainer = cc.CombContainer()
#Tablero (tensiones normales)
combContainer.ULS.perm.add("ELU007",
                           "1.00*G12 + 1.00*G3 + 1.35*Q1a1 + 0.90*Q33")
combContainer.ULS.perm.add("ELU010",
                           "1.00*G12 + 1.00*G3 + 1.35*Q1a1 + 0.90*Q22")
combContainer.ULS.perm.add("ELU034",
                           "1.00*G12 + 1.00*G3 + 1.35*Q1c + 0.90*Q22")
combContainer.ULS.perm.add("ELU076", "1.00*G12 + 1.00*G3 + 0.6*Q1e + 1.50*Q33")
combContainer.ULS.perm.add("ELU120",
                           "1.00*G12 + 1.35*G3 + 1.35*Q1c + 0.90*Q22")
combContainer.ULS.perm.add("ELU179",
                           "1.35*G12 + 1.00*G3 + 1.35*Q1a1 + 0.90*Q33")
combContainer.ULS.perm.add("ELU182",
                           "1.35*G12 + 1.00*G3 + 1.35*Q1a1 + 0.90*Q22")
combContainer.ULS.perm.add("ELU190",
                           "1.35*G12 + 1.00*G3 + 1.35*Q1b1 + 0.90*Q34")
combContainer.ULS.perm.add("ELU192",
                           "1.35*G12 + 1.00*G3 + 1.35*Q1b1 + 0.90*Q32")
combContainer.ULS.perm.add("ELU198",
                           "1.35*G12 + 1.00*G3 + 1.35*Q1b2 + 0.90*Q32")
combContainer.ULS.perm.add("ELU215",
                           "1.35*G12 + 1.00*G3 + 1.35*Q1e + 0.90*Q33")
combContainer.ULS.perm.add("ELU268",
                           "1.35*G12 + 1.35*G3 + 1.35*Q1a1 + 0.90*Q22")
combContainer.ULS.perm.add("ELU286",
                           "1.35*G12 + 1.35*G3 + 1.35*Q1b2 + 0.90*Q22")
Ejemplo n.º 3
0
QtrafSit1unif.addLstLoads([
    trafLoadSit1Lane1Hw, trafLoadSit1Lane2Hw, trafLoadSit1Lane3Hw,
    trafLoadSit1RestHw, trafLoadSit12Lane1Sr
])  #defined only for the purpose of displaying
QtrafSit2unif = lcases.LoadCase(preprocessor=prep,
                                name='QtrafSit2unif',
                                loadPType="default",
                                timeSType="constant_ts")
QtrafSit2unif.create()
QtrafSit2unif.addLstLoads([
    trafLoadSit2Lane1Hw, trafLoadSit2Lane2Hw, trafLoadSit2Lane3Hw,
    trafLoadSit2RestHw, trafLoadSit12Lane1Sr
])  #defined only for the purpose of displaying

#LOAD COMBINATIONS
combContainer = cc.CombContainer()  #Container of load combinations

# COMBINATIONS OF ACTIONS FOR SERVICEABILITY LIMIT STATES
# name:        name to identify the combination
# rare:        combination for a rare design situation
# freq:        combination for a frequent design situation
# qp:          combination for a quasi-permanent design situation
# earthquake:  combination for a seismic design situation
#Characteristic combinations.
combContainer.SLS.rare.add(
    'ELSR01', '1.0*GselfWeight+1.0*GdeadLoad+1.0*GearthPress+1.0*QtrafSit1a')
combContainer.SLS.rare.add(
    'ELSR02', '1.0*GselfWeight+1.0*GdeadLoad+1.0*GearthPress+1.0*QtrafSit1b')
combContainer.SLS.rare.add(
    'ELSR03', '1.0*GselfWeight+1.0*GdeadLoad+1.0*GearthPress+1.0*QtrafSit2a')
combContainer.SLS.rare.add(
Ejemplo n.º 4
0
    gammaWater=1000 * gravity,
    qUnif=0)
wall.createBackFillPressures(backFillPressureModelAcc, Delta=backFillDelta)
zGroundFrontFill = zGroundBackFill - wall.stemHeight + frontFillDepth  #Front fill
frontFillPressureModel = earth_pressure.EarthPressureModel(
    zGround=zGroundFrontFill,
    zBottomSoils=[-1e3],
    KSoils=[Ka],
    gammaSoils=[gSoil],
    zWater=-1e3,
    gammaWater=1000 * gravity,
    qUnif=0)
wall.createFrontFillPressures(frontFillPressureModel)

#Load combinations
combContainer = combinations.CombContainer()

## Quasi-permanent situations.
combContainer.SLS.qp.add('SLS01', '1.0*selfWeight+1.0*earthPress')
slsCombinations = ['SLS01']

## Stability ultimate states. (type 1)
combContainer.ULS.perm.add('ULS01', '0.9*selfWeight+1.00*earthPress')
combContainer.ULS.perm.add('ULS02', '1.1*selfWeight+1.00*earthPress')
combContainer.ULS.perm.add('ULS03', '0.9*selfWeight+1.50*earthPress')
combContainer.ULS.perm.add('ULS04', '1.1*selfWeight+1.50*earthPress')
combContainer.ULS.perm.add('ULS05', '0.9*selfWeight+1.00*earthPressAcc')
combContainer.ULS.perm.add('ULS06', '1.1*selfWeight+1.00*earthPressAcc')
stabilityULSCombinations = [
    'ULS01', 'ULS02', 'ULS03', 'ULS04', 'ULS05', 'ULS06'
]