hi,everyone.I am a new guy studying opensees.
there are several questions about opensees .
I study the example of Solid-Fluid Fully Coupled elements,using the block3d to mesh
the model.
but,uh, there are two question,
1 how to define the nodes of the new model meshing by "block3d"?
2 After that,there are several elements ,if I want to 'recorder' the data of the
element, like this "eval "recorder Element -ele 1 -time -file PRESS6_2_2.out -dT
0.01 material 6 pressure",how to do this ? maybe this is the same question as the
question 1, firstly how to define the elements?define by myself,or following the
opensees results,but how to do this?how does the nodes and elements define of?
thx very much!
question about block3d
Moderator: Moderators
-
- Posts: 5
- Joined: Fri Apr 03, 2009 2:28 am
- Location: wuhan
Re: question about block3d
anyone here ?
please help me
can you tell me some rules about block3d?
please help me
can you tell me some rules about block3d?
-
- Posts: 73
- Joined: Tue Feb 20, 2007 6:48 am
- Location: IIEES
- Contact:
this is an example:
wipe all
model BasicBuilder -ndm 3 -ndf 3
set nx 1
set ny 4
set nz 20
set e1 1
set n1 1
set dx 1.0
set dy 120.0
set dz 60.0
nDMaterial PressureIndependMultiYield 1 3 1.9 182590.0 852086.666667 10.0 0.1 21.4 50 0 -9 1.e-06 1 3.16e-06 0.913 1.e-05 0.761 3.16e-05 0.565 1.e-04 0.4 3.16e-04 0.261 1.e-03 0.152 3.16e-03 0.076 1.e-02 0.037
# block3D $nx $ny $nz $e1 $n1 stdBrick 1 {
# 1 $dx $dy 0
# 2 0 $dy 0
# 3 0 0 0
# 4 $dx 0 0
# 5 $dx $dy $dz
# 6 0 $dy $dz
# 7 0 0 $dz
# 8 $dx 0 $dz
# }
block3D $nx $ny $nz $e1 $n1 stdBrick 1 {
1 1 60 0
2 0 60 0
3 0 0 0
4 1 0 0
5 1 60 120
6 0 60 120
7 0 0 120
8 1 0 120
}
fixZ 0.0 1 1 1
fixX 0.0 1 0 0
fixX $dx 1 0 0
updateMaterialStage -material 1 -stage 0
#############################################################
# GRAVITY APPLICATION (elastic behavior)
system ProfileSPD
test NormDispIncr 1.e-5 10 0
algorithm ModifiedNewton
constraints Transformation
integrator LoadControl 1 1 1 1
numberer RCM
analysis Static
analyze 1
# switch the material to plastic
updateMaterialStage -material 1 -stage 1
# NOW APPLY LOADING SEQUENCE AND ANALYZE
# rezero time
loadConst -time 0.0
wipeAnalysis
source inputMotion.tcl
set saveFolder data
file mkdir ./$saveFolder
set nodeRec [expr ($nx+1)*($ny+1)*(2*$nz+1)/2]
recorder Node -file $saveFolder/acc -time -node $nodeRec -dof 2 -dT 0.033 accel
recorder Node -file $saveFolder/disp -time -node $nodeRec -dof 3 -dT 0.033 disp
source display.tcl
# create the Analysis
constraints Transformation
test NormDispIncr 1.e-5 10 0
numberer RCM
algorithm Newton
system ProfileSPD
integrator Newmark 0.55 0.275625
analysis VariableTransient
#analyze
set deltaT 0.033
set startT [clock seconds]
analyze 1362 $deltaT [expr $deltaT/100] $deltaT 5
set endT [clock seconds]
puts "Execution time: [expr $endT-$startT] seconds."
wipe #flush ouput stream
wipe all
model BasicBuilder -ndm 3 -ndf 3
set nx 1
set ny 4
set nz 20
set e1 1
set n1 1
set dx 1.0
set dy 120.0
set dz 60.0
nDMaterial PressureIndependMultiYield 1 3 1.9 182590.0 852086.666667 10.0 0.1 21.4 50 0 -9 1.e-06 1 3.16e-06 0.913 1.e-05 0.761 3.16e-05 0.565 1.e-04 0.4 3.16e-04 0.261 1.e-03 0.152 3.16e-03 0.076 1.e-02 0.037
# block3D $nx $ny $nz $e1 $n1 stdBrick 1 {
# 1 $dx $dy 0
# 2 0 $dy 0
# 3 0 0 0
# 4 $dx 0 0
# 5 $dx $dy $dz
# 6 0 $dy $dz
# 7 0 0 $dz
# 8 $dx 0 $dz
# }
block3D $nx $ny $nz $e1 $n1 stdBrick 1 {
1 1 60 0
2 0 60 0
3 0 0 0
4 1 0 0
5 1 60 120
6 0 60 120
7 0 0 120
8 1 0 120
}
fixZ 0.0 1 1 1
fixX 0.0 1 0 0
fixX $dx 1 0 0
updateMaterialStage -material 1 -stage 0
#############################################################
# GRAVITY APPLICATION (elastic behavior)
system ProfileSPD
test NormDispIncr 1.e-5 10 0
algorithm ModifiedNewton
constraints Transformation
integrator LoadControl 1 1 1 1
numberer RCM
analysis Static
analyze 1
# switch the material to plastic
updateMaterialStage -material 1 -stage 1
# NOW APPLY LOADING SEQUENCE AND ANALYZE
# rezero time
loadConst -time 0.0
wipeAnalysis
source inputMotion.tcl
set saveFolder data
file mkdir ./$saveFolder
set nodeRec [expr ($nx+1)*($ny+1)*(2*$nz+1)/2]
recorder Node -file $saveFolder/acc -time -node $nodeRec -dof 2 -dT 0.033 accel
recorder Node -file $saveFolder/disp -time -node $nodeRec -dof 3 -dT 0.033 disp
source display.tcl
# create the Analysis
constraints Transformation
test NormDispIncr 1.e-5 10 0
numberer RCM
algorithm Newton
system ProfileSPD
integrator Newmark 0.55 0.275625
analysis VariableTransient
#analyze
set deltaT 0.033
set startT [clock seconds]
analyze 1362 $deltaT [expr $deltaT/100] $deltaT 5
set endT [clock seconds]
puts "Execution time: [expr $endT-$startT] seconds."
wipe #flush ouput stream
operating system: Debian Gnu/Linux lenny 64Bit debian (DOT) org/
email : e.roknabadi (AT) iiees.ac.ir
email : e.roknabadi (AT) iiees.ac.ir