hi, I'm trying to build a very simple pile model for the analysis lateral behavior.
I wanted to make a spring using zerolength element, but the my coding has some error below.
' warning SuperLU:: solve<void>- Error 14 returned in factorization dgstrf<>
warning Linear:: solveCurrentStep<>-the LinearSOE failed in solve<>
StaticAnalysis::analyze<>-the Algorithm failed at iteration: 0 with domain at load factor 1'
I attaced my coding, is there something wrong?
==========================================================================
# Date: 2012. 09. 4 Tuesday
# Written by Young-Eun Jang for rigid model
# Unit N, m
wipe all
# pile properties, ElasticBeamColumn, D =3.5m, E = 210,000kPa, Iz = 0.6507 m4
set E_beam 210000000
set I_beam 0.6507
set A_beam 3.9
########################################
# Beam's dimension & DOF #
########################################
model BasicBuilder -ndm 2 -ndf 3
########################################
# Create beam node #
########################################
node 1 0 0
node 2 0 1
node 3 0 2
node 4 0 3
node 5 0 4
node 6 0 5
puts "Finished create all nodes"
########################################
# Create beam Element #
########################################
geomTransf Linear 1
element elasticBeamColumn 1 1 2 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 2 2 3 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 3 3 4 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 4 4 5 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 5 5 6 $A_beam $E_beam $I_beam 1
puts "Finished create all elements"
################################################
############# Add Node for 6x6 #################
################################################
node 7 0 0
################################################
######### 6x6 Springs with pile node 1 #########
################################################
uniaxialMaterial Elastic 1 50000
uniaxialMaterial Elastic 2 30000
uniaxialMaterial Elastic 3 80000
# lateral,horizontal is x & z , vertical is y
# x y z(translation) x y z(rotation)
# 1 2 3 4 5 6
element zeroLength 6 1 7 -mat 1 2 3 -dir 1 3 5
########################################
# Boundary Condition #
########################################
fix 1 0 1 0 1 0 1
puts "Finished create boundary condition"
########################################
# Define load #
########################################
timeSeries Linear 1
pattern Plain 1 1 {
#Create the nodal load - command: load nodeID xForce yForce
load 6 100.0 0.0 0.0
}
puts "Finished apply load"
########################################
# recorder #
########################################
# nodal displacement
recorder Node -file Data/disp1.out -time -node 1 -dof 1 disp
recorder Node -file Data/disp2.out -time -node 2 -dof 1 disp
recorder Node -file Data/disp3.out -time -node 3 -dof 1 disp
recorder Node -file Data/disp4.out -time -node 4 -dof 1 disp
recorder Node -file Data/disp5.out -time -node 5 -dof 1 disp
recorder Node -file Data/disp6.out -time -node 6 -dof 1 disp
# nodel bending moment
#########################################################################
# #
# | → end node 2'shear, axial, and bending moment #
# | #
# | #
# | #
# | → end node 1's shear, axial, and bending moment #
# #
#########################################################################
recorder Element -file Data/Fele1.out -time -ele 1 force
recorder Element -file Data/Fele2.out -time -ele 2 force
recorder Element -file Data/Fele3.out -time -ele 3 force
recorder Element -file Data/Fele4.out -time -ele 4 force
recorder Element -file Data/Fele5.out -time -ele 5 force
puts "Finished create recorder"
########################################
# analysis #
########################################
integrator LoadControl 1.0
#integrator Newmark 0.5 0.25
test EnergyIncr 1.0e-12 10 0
algorithm Linear
numberer RCM
constraints Plain
system SparseGeneral
#system ProfileSPD
analysis Static
analyze 1
puts "Finished analysis"
wipe all
==========================================================================
zerolength element
Moderator: Moderators