bandgenlinklapacksolver warning error

Forum for OpenSees users to post questions, comments, etc. on the use of the OpenSees interpreter, OpenSees.exe

Moderators: silvia, selimgunay, Moderators

Post Reply
mijin
Posts: 8
Joined: Tue Sep 15, 2015 5:16 am
Location: PNU

bandgenlinklapacksolver warning error

Post by mijin »

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
fmk
Site Admin
Posts: 5884
Joined: Fri Jun 11, 2004 2:33 pm
Location: UC Berkeley
Contact:

Re: bandgenlinklapacksolver warning error

Post by fmk »

look at the element connectivity a.k.a the nodes at ends of the 2 elements of yours
mijin
Posts: 8
Joined: Tue Sep 15, 2015 5:16 am
Location: PNU

Re: bandgenlinklapacksolver warning error

Post by mijin »

i don't know the problem of connecting node..
what is problem nodes at ends of the 2 elements of mine?
fmk
Site Admin
Posts: 5884
Joined: Fri Jun 11, 2004 2:33 pm
Location: UC Berkeley
Contact:

Re: bandgenlinklapacksolver warning error

Post by fmk »

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.
Post Reply