Bond_SP01, multilinear
Moderators: silvia, selimgunay, Moderators
Bond_SP01, multilinear
Hello guys.
I used "uniaxialmaterial Bond_SP01" command to make yield strength. but error comes out.
So I tried "uniaxialMaterial Mulilinear" command, but same error comes out.
error is
WARNING SuperLU::solve(void)- Error 1 returned in factorization dgstrf()
WARNING ModifiedNewton::solveCurrentStep() -the LinearSysOfEqn failed in solve()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 0.01
OpenSees > analyze failed, returned: -3 error flag
What is problem?
Thanks for answering
---------------------------------------------------------------------------------------------
# units: kg, mm, sec, N
wipe all
# Define Acc file
source model.tcl
set Accfile "$acc.acc"
set outputfile "$output-non(3)"
set dt $dt_acc
set Nsteps $Nsteps_acc
# Define Save
file mkdir $outputfile
# Define model
model BasicBuilder -ndm 2 -ndf 3
set h 3500
set m 5.908
set g 9800
set W [expr $m*$g]
node 1 0.0 0.0
node 2 0.0 $h -mass $m 1e-9 0.0
fix 1 1 1 1
fix 2 0 0 1
set Tn 1; # sec (Natural period)
set pi [expr acos(-1.0)];
# Columns and Beam Properties
set K [expr pow(2*$pi/$Tn,2)*$m]; # KN/mm
set alpha 0
set Fy 0.255
set E $K
set b 0
set Su 0.00328
set Sy [expr $Su*(1/3)]
set Fu $Fy
# Define Material
set materialTag 1
uniaxialMaterial Bond_SP01 $materialTag $Fy $Sy $Fu $Su $b 1.0
#uniaxialMaterial MultiLinear $materialTag $Sy $Fy $Su $Fu
# Define Element
set elementTag 1
element twoNodeLink $elementTag 1 2 -mat $materialTag -dir 2 -doRayleigh;
puts "Model Built"
# Define timeSeries
timeSeries Path 1001 -dt $dt -filePath $Accfile -factor 1;
# Define RECORDERS
recorder Node -file $outputfile/Aabs.out -timeSeries 1001 -time -node 2 -dof 1 accel;
recorder Node -file $outputfile/Afree.out -time -time -node 2 -dof 1 accel; # displacements of free nodes
recorder Element -file $outputfile/Deform.out -time -ele 1 deformations ;
recorder Element -file $outputfile/Force.out -time -ele 1 globalForce ;
# Define Damping
#set damping based on first eigen mode
#set lam [eigen 1]
#set omega [expr sqrt($lam)]
#puts "DoneDone"
#set omega [expr sqrt($K/$m)]
set omega [expr sqrt($K/$m)]
set period [expr 2*($pi/$omega)]
set freq [expr 1/$period]
set damp 0.05;
rayleigh 0 0 0 [expr 2.*$damp/($omega)]
#pattern UniformExcitation $patternTag $dir -accel $tsTag <-vel0 $ver0>
pattern UniformExcitation 1 1 -accel 1001; # define where and how (pattern tag, dof) acceleration is applied
# create the analysis
#wipeAnalysis; # clear previously-define analysis parameters
constraints Transformation; # how it handles boundary conditions
numberer Plain; # renumber dof's to minimize band-width (optimization), if you want to
system SparseGeneral -piv; # how to store and solve the system of equations in the analysis (large model: try UmfPack)
test EnergyIncr 1.0e-8 100 0; # test Eneregy incerment
algorithm ModifiedNewton; # use Kyrlow-Newton algorithm
integrator Newmark 0.5 0.25 ; # determine the next time step for an analysis
analysis Transient; # define type of analysis: time-dependent
analyze $Nsteps $dt; # apply 10*4096 steps for 0.001-sec time steps in analysis
wipe
puts "You are the best!"
I used "uniaxialmaterial Bond_SP01" command to make yield strength. but error comes out.
So I tried "uniaxialMaterial Mulilinear" command, but same error comes out.
error is
WARNING SuperLU::solve(void)- Error 1 returned in factorization dgstrf()
WARNING ModifiedNewton::solveCurrentStep() -the LinearSysOfEqn failed in solve()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 0.01
OpenSees > analyze failed, returned: -3 error flag
What is problem?
Thanks for answering
---------------------------------------------------------------------------------------------
# units: kg, mm, sec, N
wipe all
# Define Acc file
source model.tcl
set Accfile "$acc.acc"
set outputfile "$output-non(3)"
set dt $dt_acc
set Nsteps $Nsteps_acc
# Define Save
file mkdir $outputfile
# Define model
model BasicBuilder -ndm 2 -ndf 3
set h 3500
set m 5.908
set g 9800
set W [expr $m*$g]
node 1 0.0 0.0
node 2 0.0 $h -mass $m 1e-9 0.0
fix 1 1 1 1
fix 2 0 0 1
set Tn 1; # sec (Natural period)
set pi [expr acos(-1.0)];
# Columns and Beam Properties
set K [expr pow(2*$pi/$Tn,2)*$m]; # KN/mm
set alpha 0
set Fy 0.255
set E $K
set b 0
set Su 0.00328
set Sy [expr $Su*(1/3)]
set Fu $Fy
# Define Material
set materialTag 1
uniaxialMaterial Bond_SP01 $materialTag $Fy $Sy $Fu $Su $b 1.0
#uniaxialMaterial MultiLinear $materialTag $Sy $Fy $Su $Fu
# Define Element
set elementTag 1
element twoNodeLink $elementTag 1 2 -mat $materialTag -dir 2 -doRayleigh;
puts "Model Built"
# Define timeSeries
timeSeries Path 1001 -dt $dt -filePath $Accfile -factor 1;
# Define RECORDERS
recorder Node -file $outputfile/Aabs.out -timeSeries 1001 -time -node 2 -dof 1 accel;
recorder Node -file $outputfile/Afree.out -time -time -node 2 -dof 1 accel; # displacements of free nodes
recorder Element -file $outputfile/Deform.out -time -ele 1 deformations ;
recorder Element -file $outputfile/Force.out -time -ele 1 globalForce ;
# Define Damping
#set damping based on first eigen mode
#set lam [eigen 1]
#set omega [expr sqrt($lam)]
#puts "DoneDone"
#set omega [expr sqrt($K/$m)]
set omega [expr sqrt($K/$m)]
set period [expr 2*($pi/$omega)]
set freq [expr 1/$period]
set damp 0.05;
rayleigh 0 0 0 [expr 2.*$damp/($omega)]
#pattern UniformExcitation $patternTag $dir -accel $tsTag <-vel0 $ver0>
pattern UniformExcitation 1 1 -accel 1001; # define where and how (pattern tag, dof) acceleration is applied
# create the analysis
#wipeAnalysis; # clear previously-define analysis parameters
constraints Transformation; # how it handles boundary conditions
numberer Plain; # renumber dof's to minimize band-width (optimization), if you want to
system SparseGeneral -piv; # how to store and solve the system of equations in the analysis (large model: try UmfPack)
test EnergyIncr 1.0e-8 100 0; # test Eneregy incerment
algorithm ModifiedNewton; # use Kyrlow-Newton algorithm
integrator Newmark 0.5 0.25 ; # determine the next time step for an analysis
analysis Transient; # define type of analysis: time-dependent
analyze $Nsteps $dt; # apply 10*4096 steps for 0.001-sec time steps in analysis
wipe
puts "You are the best!"
-
- Posts: 913
- Joined: Mon Sep 09, 2013 8:50 pm
- Location: University of California, Berkeley
Re: Bond_SP01, multilinear
It looks like you need to fix the vertical DOF at node 2 as well
fix 2 1 0 1
fix 2 1 0 1
Re: Bond_SP01, multilinear
I put earthquake acceleration on direction 1 and i want to know deformation-x.selimgunay wrote: ↑Mon Aug 03, 2020 3:37 pmIt looks like you need to fix the vertical DOF at node 2 as well
fix 2 1 0 1
So i cant do 'fix 2 1 0 1'
Thank you for answering
-
- Posts: 913
- Joined: Mon Sep 09, 2013 8:50 pm
- Location: University of California, Berkeley
Re: Bond_SP01, multilinear
Then it should be 'fix 2 0 1 1' to fix the vertical dof