zerolength element

A forum dedicated to users with questions regarding soil materials and elements.

forum currently locked

Moderator: Moderators

Locked
gokelly
Posts: 1
Joined: Sun Mar 06, 2011 8:23 pm
Location: UST

zerolength element

Post by gokelly » Tue Sep 04, 2012 2:03 am

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

==========================================================================

Locked