Problem in Bilin usage.

[alima2065]

Hi dears
I want to use modified Ibarra-Krawinkler deterioration model for my steel frames. For testing I have made a simple cantilever beam with concentrated spring at an end. But I don’t know why i don't have any reaction in node 1 that has been fixed!! I always use these recorders and i think i made another wrong. I have checked it exactly. Would please inform me what I do wrong?


wipe;
model BasicBuilder -ndm 2 -ndf 3;

node 1 0 0
node 2 0 0
node 3 3 0

fix 1 1 1 1

equalDOF 1 2 1 2

set E 2e8;

set AW30 0.005709666; # W14X30
set IW30 0.000121123; # W14X30

#----- Units are KN and m and Rad.
set eleID 30
set K 8259
set asPos 0.03
set asNeg 0.03
set MyPos 267.375
set MyNeg -267.375
set LS 1000
set LK 1000
set LA 1000
set LD 1000
set cS 1
set cK 1
set cA 1
set cD 1
set th_pP 0.032
set th_pN 0.032
set th_pcP 0.173
set th_pcN 0.173
set ResP 0.4
set ResN 0.4
set th_uP 0.4
set th_uN 0.4
set DP 1
set DN 1

uniaxialMaterial Bilin $eleID $K $asPos $asNeg $MyPos $MyNeg $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
#uniaxialMaterial Elastic 1 2e8
geomTransf PDelta 1


element elasticBeamColumn 1 2 3 $AW30 $E $IW30 1

element zeroLength 2 1 2 -mat $eleID -dir 3

#source DisplayModel2D.tcl

#source DisplayPlane.tcl


recorder Node -file node2.txt -node 2 -dof 3 disp;
recorder Node -file Rx.txt -node 1 -dof 1 2 3 reaction;

#DisplayModel2D DeformedShape

# ---------- Load Pattern for Gravity Loads ------------------

pattern Plain 100 Linear {
load 3 0 1 0
}

constraints Lagrange
numberer Plain;
system BandGeneral
test NormDispIncr 1e-5 1000 ;
algorithm Newton;

integrator DisplacementControl 3 2 0.018
analysis Static;
analyze 10

integrator DisplacementControl 3 2 -0.018
analysis Static;
analyze 20

integrator DisplacementControl 3 2 0.018
analysis Static;
analyze 20

integrator DisplacementControl 3 2 -0.018
analysis Static;
analyze 20

integrator DisplacementControl 3 2 0.018
analysis Static;
analyze 10

[johnyjohn]

nice............

[bailiagng]

Hello !
Maybe you should change "element zeroLength 2 1 2 -mat $eleID -dir 3" into "element zeroLength 2 1 2 -mat $eleID -dir 6". Try it !

[bailiagng]

I have the similiar problem.But I get a very large peroid when I do the eigen analysis. When i change the "equalDOF 1 2 1 2" into "equalDOF 1 2 1 2 3", the peroid seems reasonable. But then i do the pushover analysis , the base force is 0. Who can help me ?