Hello. i'am opensees first users.. i got this error when gravity analysis..
↓
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 2
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
I guess that its problem is modling. but i don't know where is error and why because my model looks like different other conventional model..
finally i wanna pushover analysis structure shaped angle..
shape of my model looks like example below
Example)
■ fix (node 2)
│
│
│
│
■ ─ ─ ─ ─ ─ ─ ● roller(node 3)
please help me
my script file is below..
###################################################################################################
wipe all; # clear memory of past model definitions
model BasicBuilder -ndm 2 -ndf 3; # Define the model builder, ndm = #dimension, ndf = #dofs
#source DisplayModel2D.tcl; # procedure for displaying a 2D perspective of model
#source DisplayPlane.tcl; # procedure for displaying a plane in the model
###################################################################################################
# Define Analysis Type
###################################################################################################
# Define type of analysis: "pushover" = pushover
set analysisType "pushover";
if {$analysisType == "pushover"} {
set dataDir Anglefibermodel-Pushover-Output; # name of output folder
file mkdir $dataDir; # create output folder
}
###################################################################################################
# Define Building Geometry, Nodes, and Constraints
###################################################################################################
# input data
set t 15.9 ;
set width 254 ;
set g2 114 ;
set L2 38 ;
set y 279 ;
set Es 200000 ;
set a 0.015;
set b 1.5 ;
set Weight1 80.80 ;
set Weight2 14.76 ;
set g 9810;
#nodal coordinates: unit N,mm
node 1 0 0 ;
node 2 0 137.8125;
node 3 61.8125 0 ;
# nodal masses:
mass 1 [expr $Weight1/$g] 1e-9 0.; # node#, Mx My Mz, Mass=Weight/g, neglect rotational inertia at nodes
# assign boundary condidtions
# command: fix nodeID dxFixity dLdyFixity rzFixity
# fixity values: 1 = constrained; 0 = unconstrained
# fix the base of the building; pin P-delta column at base
fix 2 1 1 1; # fix
fix 3 0 1 0; # roller
rigidLink bar 3 1;
###################################################################################################
# Define Section Properties and Elements
###################################################################################################
# define Hysterical Materials
uniaxialMaterial Hysteretic 1 279 0.001395 418.5 0.047895 -279 -0.001395 -418.5 -0.047895 1 1 0 0;
# define angle element [spring]
set nfdv 1; # number of fibers along vertical depth
set nftv 5; # number of fibers along vertical thickness
set nfbh 1; # number of fibers along horizontal width (you want this many in a bi-directional loading)
set nfth 5; # number of fibers along horizontal thickness
set yz1_v [expr -$t/2];
set yz2_v [expr $t/2];
section fiberSec 1 {
# nfIJ nfJK yI zI yJ zJ yK zK yL zL
patch quadr 1 $nfdv $nftv $yz1_v $yz1_v $yz2_v $yz1_v $yz2_v $yz2_v $yz1_v $yz2_v
;}
# define transformation
geomTransf Linear 1 ;
#element nonlinearBeamColumn $eleTag $iNode $jNode $numlntgrPts $secTag $transfTag
element nonlinearBeamColumn 1 1 2 3 1 1 ;
element nonlinearBeamColumn 2 1 3 3 1 1 ;
# Define RECORDERS -------------------------------------------------------------
recorder Node -file $dataDir/DRollder.out -time -node 3 -dof 1 disp; # displacements of roller nodes
recorder Node -file $dataDir/DFree.out -time -node 1 -dof 1 disp; # displacements of free nodes
recorder Element -file $dataDir/Force.out -time -ele 1 2 globalForce; # element forces
# define GRAVITY -------------------------------------------------------------
pattern Plain 1 Linear {
load 1 0 -80.55 152.0099
}
pattern Plain 2 Linear {
load 2 0 -80.55 152.0099; # distributed superstructure-weight on beam
}
# Gravity-analysis parameters -- load-controlled static analysis
set Tol 1.0e-8; # convergence tolerance for test
constraints Plain; # how it handles boundary conditions
numberer Plain; # renumber dof's to minimize band-width (optimization), if you want to
system BandGeneral; # how to store and solve the system of equations in the analysis
test NormDispIncr $Tol 10 ; # determine if convergence has been achieved at the end of an iteration step
algorithm Newton; # use Newton's solution algorithm: updates tangent stiffness at every iteration
set NstepGravity 10; # apply gravity in 10 steps
set DGravity [expr 1./$NstepGravity]; # first load increment;
integrator LoadControl $DGravity; # determine the next time step for an analysis
analysis Static; # define type of analysis static or transient
analyze $NstepGravity; # apply gravity
# ------------------------------------------------- maintain constant gravity loads and reset time to zero
loadConst -time 0.0
bandgenlinklapacksolver warning error
Moderators: silvia, selimgunay, Moderators
Re: bandgenlinklapacksolver warning error
look at the element connectivity a.k.a the nodes at ends of the 2 elements of yours
Re: bandgenlinklapacksolver warning error
i don't know the problem of connecting node..
what is problem nodes at ends of the 2 elements of mine?
what is problem nodes at ends of the 2 elements of mine?
Re: bandgenlinklapacksolver warning error
look at element 2 nodes 1 and 3 .. node 3 is only node fixed and it is a pin! .. so either you goofed on your boundary conditions or ele 2 was supposed to go between nodes 2 and 3.