The OpenSees Community

Dear fmk
thank you

Dear fmk
thank you

Dear fmk
best regards
i want to update parameter damping coefficient viscous material through dynamic analysis in my project. because updateparameter do not work for this problem i use remove element through dynamic analysis and again added with new parameter through analysis.
can you offer a solution for this problem.????

I need to use remove element command in my project. I encountered problem in use remove element that explain in the following two example below.

############################################
example 1

wipe all
model basic -ndm 2 -ndf 3
file mkdir rem
node 1 0. 0.
node 2 0. 0.
fix 1 1 1 1
equalDOF 1 2 2 3
mass 2 11025. 0. 0
uniaxialMaterial Elastic 1 290000.
uniaxialMaterial Viscous 2 290. 1
element zeroLength 1 1 2 -mat 1 -dir 1
element zeroLength 2 1 2 -mat 2 -dir 1
recorder Node -file rem/node2disp.txt -time -node 2 -dof 1 disp

set accel "Series -dt 0.01 -filePath elcentro.txt -factor [expr 9.81]"
pattern UniformExcitation 2 1 -accel $accel

constraints Plain
numberer Plain
system BandGeneral
test EnergyIncr 1.e-8 8
algorithm Newton
integrator Newmark 0.5 0.25
analysis VariableTransient
analyze 4000 0.01

#######################################################################
example 2
wipe all

model basic -ndm 2 -ndf 3
file mkdir rem
node 1 0. 0.
node 2 0. 0.
fix 1 1 1 1
equalDOF 1 2 2 3
mass 2 11025. 0. 0
uniaxialMaterial Elastic 1 290000.
uniaxialMaterial Viscous 2 290. 1
recorder Node -file rem/node2disp.txt -time -node 2 -dof 1 disp

set accel "Series -dt 0.01 -filePath elcentro.txt -factor [expr 9.81]"
pattern UniformExcitation 2 1 -accel $accel

constraints Plain
numberer Plain
system BandGeneral
test EnergyIncr 1.e-8 8
algorithm Newton
integrator Newmark 0.5 0.25
analysis VariableTransient


for {set i 1} {$i<4000} {incr i} {
element zeroLength 1 1 2 -mat 1 -dir 1
element zeroLength 2 1 2 -mat 2 -dir 1
analyze 1 0.01
remove element 1
remove element 2
}

################################################################
the problem is that the response for node 2 is very different in the two example so that in the two example response for node 2 should be equal.

would you like for solution this problem what should I do???????????

I need to use remove element command in my project. I encountered problem in use remove element that explain in the following two example below.

############################################
example 1

wipe all
model basic -ndm 2 -ndf 3
file mkdir rem
node 1 0. 0.
node 2 0. 0.
fix 1 1 1 1
equalDOF 1 2 2 3
mass 2 11025. 0. 0
uniaxialMaterial Elastic 1 290000.
uniaxialMaterial Viscous 2 290. 1
element zeroLength 1 1 2 -mat 1 -dir 1
element zeroLength 2 1 2 -mat 2 -dir 1
recorder Node -file rem/node2disp.txt -time -node 2 -dof 1 disp

set accel "Series -dt 0.01 -filePath elcentro.txt -factor [expr 9.81]"
pattern UniformExcitation 2 1 -accel $accel

constraints Plain
numberer Plain
system BandGeneral
test EnergyIncr 1.e-8 8
algorithm Newton
integrator Newmark 0.5 0.25
analysis VariableTransient
analyze 4000 0.01

#######################################################################
example 2
wipe all

model basic -ndm 2 -ndf 3
file mkdir rem
node 1 0. 0.
node 2 0. 0.
fix 1 1 1 1
equalDOF 1 2 2 3
mass 2 11025. 0. 0
uniaxialMaterial Elastic 1 290000.
uniaxialMaterial Viscous 2 290. 1
recorder Node -file rem/node2disp.txt -time -node 2 -dof 1 disp

set accel "Series -dt 0.01 -filePath elcentro.txt -factor [expr 9.81]"
pattern UniformExcitation 2 1 -accel $accel

constraints Plain
numberer Plain
system BandGeneral
test EnergyIncr 1.e-8 8
algorithm Newton
integrator Newmark 0.5 0.25
analysis VariableTransient


for {set i 1} {$i<4000} {incr i} {
element zeroLength 1 1 2 -mat 1 -dir 1
element zeroLength 2 1 2 -mat 2 -dir 1
analyze 1 0.01
remove element 1
remove element 2
}

################################################################
the problem is that the response for node 2 is very different in the two example so that in the two example response for node 2 should be equal.

would you like for solution this problem what should I do???????????

I need to update Parameter damping coefficient of uniaxialMaterial Viscous through dynamic analysis....updateParameter command not work for this Viscous Material.

what way suggest me for this problem in the OpenSees ???????????????
the matlab program not work for my master thesis...


the code below is simple example for my problem that Viscous material can not update through dynamic analysis..

wipe all
model basic -ndm 2 -ndf 3
file mkdir rem
node 1 0. 0.
node 2 0. 0.
fix 1 1 1 1
equalDOF 1 2 2 3
mass 2 11025. 0. 0
uniaxialMaterial Elastic 1 290000.
uniaxialMaterial Viscous 2 290. 1
element zeroLength 1 1 2 -mat 1 -dir 1
element zeroLength 2 1 2 -mat 2 -dir 1
recorder Node -file rem/node22disp.txt -time -node 2 -dof 1 disp


set accel "Series -dt 0.01 -filePath elcentro.txt -factor [expr 9.81]"
pattern UniformExcitation 2 1 -accel $accel

constraints Plain
numberer Plain
system BandGeneral
test EnergyIncr 1.e-8 8
algorithm Newton
integrator Newmark 0.5 0.25
analysis VariableTransient

for {set i 1} {$i<4000} {incr i} {
set f 200.
parameter 1 element 2 C
updateParameter 1 $f
remove parameter 1
}

would you like to explain the best idea about this problem solution??? (except link OpenSees with matlab program)!

I appreciate if somebody help me to find out a solution for my problem.

I need to update Parameter damping coefficient of uniaxialMaterial Viscous through dynamic analysis....updateParameter command not work for this Viscous Material.

what way suggest me for this problem in the OpenSees ???????????????
the matlab program not work for my master thesis...


the code below is simple example for my problem that Viscous material can not update through dynamic analysis..

wipe all
model basic -ndm 2 -ndf 3
file mkdir rem
node 1 0. 0.
node 2 0. 0.
fix 1 1 1 1
equalDOF 1 2 2 3
mass 2 11025. 0. 0
uniaxialMaterial Elastic 1 290000.
uniaxialMaterial Viscous 2 290. 1
element zeroLength 1 1 2 -mat 1 -dir 1
element zeroLength 2 1 2 -mat 2 -dir 1
recorder Node -file rem/node22disp.txt -time -node 2 -dof 1 disp


set accel "Series -dt 0.01 -filePath elcentro.txt -factor [expr 9.81]"
pattern UniformExcitation 2 1 -accel $accel

constraints Plain
numberer Plain
system BandGeneral
test EnergyIncr 1.e-8 8
algorithm Newton
integrator Newmark 0.5 0.25
analysis VariableTransient

for {set i 1} {$i<4000} {incr i} {
set f 200.
parameter 1 element 2 C
updateParameter 1 $f
remove parameter 1
}

would you like to explain the best idea about this problem solution??? (except link OpenSees with matlab program)!

I appreciate if somebody help me to find out a solution for my problem.