A perpendicular load is put at the node of 100.
I tried to run the problem but it showed as follow. Not quite not where the problem is, and hopefully got some help.
Thanks
System
Geometry In
Beam Elements In
analysis
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 1171
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 0.1
OpenSees > analyze failed, returned: -3 error flag
Static Analysis Done
# build the model
wipe
puts "System"
model basic -ndm 2 -ndf 3
set NBoxes 8.
set NSpans [expr ($NBoxes -1)]
set B1 48
set B2 54
set B3 54
set B4 54
set B5 54
set B6 54
set B7 48
set C1 12
set C2 12
set KS 50
set BL [expr ($C1 + $B1 + $B2 + $B3 + $B4 + $B5 + $B6 + $B7 + $C2)]
set NCBNds [expr ($BL + 1)]
set FXNd [expr ($C1 + $B1 + $B2 + $B3 + $B4 + 1)]
# add nodes - command: node nodeId xCrd yCrd
for {set i 0} {$i <= [expr ($NCBNds - 1)]} {incr i 1} {
node $i [expr ($i - 1)] 0.0
}
fix $FXNd 1 0 0
set xb 0
set ni [expr ($NCBNds + 1)]
set xB1 [expr ($xb + $C1)]; set na1 [expr ($ni + 0)];
node $na1 $xB1 0.00;
fix $na1 1 1 1
set xB2 [expr ($xB1 + $B1)]; set na2 [expr ($ni + 1)];
node $na2 $xB2 0.00;
fix $na2 1 1 1
set xB3 [expr ($xB2 + $B2)]; set na3 [expr ($ni + 2)];
node $na3 $xB3 0.00;
fix $na3 1 1 1
set xB4 [expr ($xB3 + $B3)]; set na4 [expr ($ni + 3)];
node $na4 $xB4 0.00;
fix $na4 1 1 1
set xB5 [expr ($xB4 + $B4)]; set na5 [expr ($ni + 4)];
node $na5 $xB5 0.00;
fix $na5 1 1 1
set xB6 [expr ($xB5 + $B5)]; set na6 [expr ($ni + 5)];
node $na6 $xB6 0.00;
fix $na6 1 1 1
set xB7 [expr ($xB6 + $B6)]; set na7 [expr ($ni + 6)];
node $na7 $xB7 0.00;
fix $na7 1 1 1
set xB8 [expr ($xB7 + $B7)]; set na8 [expr ($ni + 7)];
node $na8 $xB8 0.00;
fix $na8 1 1 1
puts "Geometry In"
# add material - command: material <matType> matID <matArgs>
uniaxialMaterial Elastic 1 $KS
#uniaxialMaterial Elastic 1 100
set E 29000
set A 9.855
set Iz 21.2085
geomTransf Linear 1
set a 0
for {set n1 $a} {$n1 <= [expr ($BL - 1)]} {incr n1 1} {
element elasticBeamColumn $n1 $n1 [expr ($n1 + 1)] $A $E $Iz 1
}
puts "Beam Elements In"
element zeroLength [expr ($BL + 1)] [expr ($xB1 + 1)] $na1 -mat 1 -dir 1
element zeroLength [expr ($BL + 2)] [expr ($xB2 + 1)] $na2 -mat 1 -dir 1
element zeroLength [expr ($BL + 3)] [expr ($xB3 + 1)] $na3 -mat 1 -dir 1
element zeroLength [expr ($BL + 4)] [expr ($xB4 + 1)] $na4 -mat 1 -dir 1
element zeroLength [expr ($BL + 5)] [expr ($xB5 + 1)] $na5 -mat 1 -dir 1
element zeroLength [expr ($BL + 6)] [expr ($xB6 + 1)] $na6 -mat 1 -dir 1
element zeroLength [expr ($BL + 7)] [expr ($xB7 + 1)] $na7 -mat 1 -dir 1
element zeroLength [expr ($BL +
] [expr ($xB8 + 1)] $na8 -mat 1 -dir 1#set Ld 20
#set a -1
timeSeries Linear 1
#command pattern Plain $tag $timeSeriesTag { $loads }
#For uniformly-distributed load:
#eleLoad -ele $eleTag1 <$eleTag2 ....> -type -beamUniform $Wz <$Wx>
pattern Plain 1 1 {
load 100 0.0 -1.0 0.0
}
recorder display frame 200 0 500 500 -wipe
vup 0 1 0 # dirn defining up direction of view plane(ghab bordar movazi)
vpn 0 0 1 # direction of outward normal to view plane(amud bar safhe ghab)
prp 0 0 1 # eye location in local coord sys defined by viewing system
viewWindow -5 5 -5 5
display 1 2 1
#
puts "analysis"
# Create the system of equation
system BandGeneral
# Create the DOF numberer, the reverse Cuthill-McKee algorithm
numberer Plain
# Create the constraint handler, a Plain handler is used as homo constraints
constraints Plain
test NormDispIncr 1e-6 50
# Create the integration scheme, the LoadControl scheme using steps of 1.0
integrator LoadControl 0.1
# Create the solution algorithm, a Linear algorithm is created
algorithm Newton
# create the analysis object
analysis Static
recorder Element -file eleGlobal.out -time -ele all globalForces
analyze 10
puts "Static Analysis Done"