joining 2dof with 3dof elements

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

Moderators: silvia, selimgunay, Moderators

partla
Posts: 35
Joined: Wed Feb 01, 2012 3:04 am
Location: Urmia university

Re: joining 2dof with 3dof elements

Post by partla » Tue Nov 27, 2012 11:09 am

Thank you vesna,

I'll check it

partla
Posts: 35
Joined: Wed Feb 01, 2012 3:04 am
Location: Urmia university

Re: joining 2dof with 3dof elements

Post by partla » Mon Dec 03, 2012 11:43 am

hi vesna

I wrote a smaller script for my problem and in that case my analyses completed. but again in big script the analyses have errors:(
the big script is exactly same as the small one just number of elements is more. what is the problem in your idea? here is the small script. please take a look at it.
thank you

wipe

#-----------------------------------------------------------------------------------------
# 1. CREATE SOIL NODES
#-----------------------------------------------------------------------------------------
model BasicBuilder -ndm 2 -ndf 2

# define soil nodes
node 1 0 0
node 2 0.5 0
node 3 1 0
node 4 1.5 0
node 5 2 0
node 6 2.5 0
node 7 3 0
node 8 0 0.5
node 9 0.5 0.5
node 10 1 0.5
node 11 1.5 0.5
node 12 2 0.5
node 13 2.5 0.5
node 14 3 0.5
node 15 0 1
node 16 0.5 1
node 17 1 1
node 18 1.5 1
node 19 2 1
node 20 2.5 1
node 21 3 1
node 22 0 1.5
node 23 0.5 1.5
node 24 1 1.5
node 25 2 1.5
node 26 2.5 1.5
node 27 3 1.5
node 28 0 2
node 29 0.5 2
node 30 1 2
node 31 2 2
node 32 2.5 2
node 33 3 2

puts "Finshed creating all soil nodes..."

#-----------------------------------------------------------------------------------------
# 2. Define SOIL materials
#-----------------------------------------------------------------------------------------
# Material "Nevada dry sand (Dr=74%)": matTag nd rho Gr Br phi gammamax pr d phipt contrac dilat1 dilat2 liqf1 liqf2 liqf3
# noYS e cs1 cs2 cs3 pa
nDMaterial PressureDependMultiYield 1 2 +1.692 +5.300000E+004 +1.150000E+005 +3.500000E+001 +1.000000E-001 +5.400000E+001 +5.000000E-001 +2.700000E+001 +5.000000E-002 +6.000000E-001 +3.000000E+000 +0.000000E+000 +0.000000E+000 +0.000000E+000 \
11 +5.660000E-001 +9.000000E-001 +2.000000E-002 +7.000000E-001 +1.010000E+002

updateMaterialStage -material 1 -stage 0
# body force in x-direction
set xWgt1 0.00
# body force in y-direction
set yWgt1 [expr -9.81*1.692]

puts "Finished creating all soil materials..."

#-----------------------------------------------------------------------------------------
# 3. define quadElement for soil
#-----------------------------------------------------------------------------------------

element quad 1 1 2 9 8 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 2 2 3 10 9 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 3 3 4 11 10 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 4 4 5 12 11 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 5 5 6 13 12 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 6 6 7 14 13 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 7 8 9 16 15 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 8 9 10 17 16 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 9 10 11 18 17 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 10 11 12 19 18 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 11 12 13 20 19 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 12 13 14 21 20 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 13 15 16 23 22 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 14 16 17 24 23 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 15 19 20 26 25 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 16 20 21 27 26 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 17 22 23 29 28 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 18 23 24 30 29 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 19 25 26 32 31 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 20 26 27 33 32 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1

puts "Finished creating all soil elements..."

#-----------------------------------------------------------------------------------------
# 4. define fixities for Bottom boundary
#-----------------------------------------------------------------------------------------
# SPC tag Dx Dy Rz
fix 1 1 1 1
fix 2 1 1 1
fix 3 1 1 1
fix 4 1 1 1
fix 5 1 1 1
fix 6 1 1 1
fix 7 1 1 1

puts "Finished creating all fixities..."

#-----------------------------------------------------------------------------------------
# 5. define equal DOF (soil to soil)
#-----------------------------------------------------------------------------------------
equalDOF 1 7 1 2
equalDOF 8 14 1 2
equalDOF 15 21 1 2
equalDOF 22 28 1 2
equalDOF 29 35 1 2

puts "Finished creating equal DOF (soil to soil)..."

#-----------------------------------------------------------------------------------------
# 6. creat zerolength elements nodes
#-----------------------------------------------------------------------------------------
node 34 1 1
node 35 1 1.5
node 36 1 2
node 37 1.5 1
node 38 2 1
node 39 2 1.5
node 40 2 2
puts "Finished creating zerolength elements nodes..."

#-----------------------------------------------------------------------------------------
# 7. Define zerolength elements materials
#-----------------------------------------------------------------------------------------

# Material: matTag E
uniaxialMaterial ENT 2 +4.000000E+004
uniaxialMaterial Viscous 3 1200 1
uniaxialMaterial Parallel 4 2 3

puts "Finished creating all Zerolength elements materials..."

element zeroLength 21 17 34 -mat 4 -dir 1
element zeroLength 22 24 35 -mat 4 -dir 1
element zeroLength 23 30 36 -mat 4 -dir 1
element zeroLength 24 18 37 -mat 4 -dir 2
element zeroLength 25 19 38 -mat 4 -dir 1
element zeroLength 26 25 39 -mat 4 -dir 1
element zeroLength 27 31 40 -mat 4 -dir 1

puts "Finished creating zerolength elements ..."

#-----------------------------------------------------------------------------------------
# 9. create the 1st analysis method & analyze
#-----------------------------------------------------------------------------------------
system BandGeneral
test NormDispIncr 1e-5 20
algorithm NewtonLineSearch
constraints Transformation
integrator LoadControl 1 1 1 1
numberer RCM
analysis Static
analyze 1

puts "finished 1st analysis (static)"

#-----------------------------------------------------------------------------------------
# 10. update Material Stage
#-----------------------------------------------------------------------------------------
updateMaterialStage -material 1 -stage 1

#-----------------------------------------------------------------------------------------
# 11. Creat Wall and base nodes
#-----------------------------------------------------------------------------------------
model BasicBuilder -ndm 2 -ndf 3

node 41 1 1
node 42 1 1.5
node 43 1 2
node 44 1.5 1
node 45 2 1
node 46 2 1.5
node 47 2 2

puts "Finished creating wall and base nodes ..."

#-----------------------------------------------------------------------------------------
# 12. define equal DOF (soil to wall)
#-----------------------------------------------------------------------------------------

equalDOF 41 34 1 2
equalDOF 42 35 1 2
equalDOF 43 36 1 2
equalDOF 44 37 1 2
equalDOF 45 38 1 2
equalDOF 46 39 1 2
equalDOF 47 40 1 2

puts "Finished defining equal DOF (soil to wall) ..."

#-----------------------------------------------------------------------------------------
# 13. define wall and base elements
#-----------------------------------------------------------------------------------------
geomTransf Linear 1

element elasticBeamColumn 28 41 42 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 29 42 43 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 30 41 44 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 31 44 45 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 32 45 46 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 33 46 47 1.40E-01 70000000 2.43E-03 1


puts "Finished creating wall and base elements ..."

#-----------------------------------------------------------------------------------------
# 14. settime and wipeAnalysis
#-----------------------------------------------------------------------------------------
wipeAnalysis
loadConst -time 0.0

#-----------------------------------------------------------------------------------------
# 15. apply the weight of structure
#-----------------------------------------------------------------------------------------
set wstiffwall -2.38
set wbase -19.62
set wbs -22
set wstiffwallh -1.19
set wflexiblewallh -1.015

pattern Plain 5 Linear {
load 43 0.00E+00 $wstiffwall 0.00E+00
load 42 0.00E+00 $wstiffwall 0.00E+00
load 41 0.00E+00 $wbs 0.00E+00
load 44 0.00E+00 $wbase 0.00E+00
load 45 0.00E+00 $wbs 0.00E+00
load 46 0.00E+00 $wstiffwall 0.00E+00
load 47 0.00E+00 $wstiffwall 0.00E+00
}
puts "Finished applying the structures weight ..."

#-----------------------------------------------------------------------------------------
# 16. Recorder
#-----------------------------------------------------------------------------------------

# Define damping parameters
rayleigh 0 0.0044 0 0
#-----------------------------------------------------------------------------------------
# 17. create the 2nd analysis method & analyze
#-----------------------------------------------------------------------------------------

constraints Transformation
test NormDispIncr +1.000000E-005 40 1
integrator LoadControl 0.05
algorithm NewtonLineSearch
numberer RCM
system BandGeneral
analysis Static

set startT [clock seconds]

analyze 40

set endT [clock seconds]
puts "finished 2nd analysis (static)"


#-----------------------------------------------------------------------------------------
# 18. earthquake excitation
#-----------------------------------------------------------------------------------------

wipeAnalysis
loadConst -time 0.0

timeSeries Path 3 -dt +0.005 -filePath loma.txt -factor +9.810000E+000
# LoadPattern "loma": patternTag dir tsTag
pattern UniformExcitation 3 1 -accel 3

#-----------------------------------------------------------------------------------------
# 19. Define recorder
#-----------------------------------------------------------------------------------------

#-----------------------------------------------------------------------------------------
# 20. create the 3rd analysis method & analyze
#-----------------------------------------------------------------------------------------

constraints Penalty 1e12 1e12
test NormDispIncr +1.000000E-004 10 1
integrator Newmark 0.6 0.3025
algorithm NewtonLineSearch
numberer RCM
system BandGeneral
analysis Transient

set startT [clock seconds]

analyze 7989 0.005

set endT [clock seconds]

puts "[expr $endT-$startT] seconds"

puts "finished 3rd analysis (dynamic)"

vesna
Posts: 3033
Joined: Tue May 23, 2006 11:23 am
Location: UC Berkeley

Re: joining 2dof with 3dof elements

Post by vesna » Mon Dec 03, 2012 12:24 pm

Post the error messages.

Expand your model slowly. This way you will be able to capture the errors more successfully.

partla
Posts: 35
Joined: Wed Feb 01, 2012 3:04 am
Location: Urmia university

Re: joining 2dof with 3dof elements

Post by partla » Mon Dec 03, 2012 12:27 pm

thank you for reply,
here is the warnings.

Finshed creating all soil nodes...
WARNING: updateMaterialStage - MaterialStageParameter::setDomain() - no effect w
ith material tag 1
Finished creating all soil materials...
Finished creating all soil elements...
Finished creating all fixities...
Finished creating equal DOF (soil to soil)...
Finished creating zerolength elements nodes...
Finished creating all Zerolength elements materials...
Finished creating zerolength elements ...
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 2
WARNING NewtonLineSearch::solveCurrentStep() -the LinearSysOfEqn failed in solve
()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 1
OpenSees > analyze failed, returned: -3 error flag
finished 1st analysis (static)
Finished creating wall and base nodes ...
Finished defining equal DOF (soil to wall) ...
Finished creating wall and base elements ...
Finished applying the structures weight ...
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 2
WARNING NewtonLineSearch::solveCurrentStep() -the LinearSysOfEqn failed in solve
()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 0.05
OpenSees > analyze failed, returned: -3 error flag
finished 2nd analysis (static)
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 17
WARNING NewtonLineSearch::solveCurrentStep() -the LinearSysOfEqn failed in solve
()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 0.005
OpenSees > analyze failed, returned: -3 error flag
0 seconds
finished 3rd analysis (dynamic)

vesna
Posts: 3033
Joined: Tue May 23, 2006 11:23 am
Location: UC Berkeley

Re: joining 2dof with 3dof elements

Post by vesna » Mon Dec 03, 2012 2:05 pm

the messages indicate the problem with the model.

partla
Posts: 35
Joined: Wed Feb 01, 2012 3:04 am
Location: Urmia university

Re: joining 2dof with 3dof elements

Post by partla » Wed Dec 12, 2012 7:06 am

thankyou vesna you were right. the problem was solved and now I can run my sccript. but again I have some questions. I use zerolength elements as interface elements of soil and structure and assign springs for these elements. but I'm not sure wich materials is good for a nonlinear springs. could you please help me?

thanks again

vesna
Posts: 3033
Joined: Tue May 23, 2006 11:23 am
Location: UC Berkeley

Re: joining 2dof with 3dof elements

Post by vesna » Wed Dec 12, 2012 1:30 pm

If you show me a backbone curve of the material that you are trying to model with your spring I can tell you what OpenSees material to use.

partla
Posts: 35
Joined: Wed Feb 01, 2012 3:04 am
Location: Urmia university

Re: joining 2dof with 3dof elements

Post by partla » Sat Dec 15, 2012 7:40 am

I'm not sure about the backbone curve of materials. As I mentioned before I'm working on report peer 2008/104 and here is the explanation of that report about springs. whats your idea about it? I need the parametrs of two material that they used.

"The soil-structure interaction was simulated by zero-length nonlinear springs. Each nonlinear spring consisted of an elastic-no-tension component in parallel with a viscous component or a dashpot representing radiation damping. The properties of the nonlinear springs were generally selected to be proportional to the strength and the damping of the adjacent soil elements"

Post Reply