STRANGE BEHAVIOUR

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

Moderators: silvia, selimgunay, Moderators

Post Reply
sedatacemogluone
Posts: 4
Joined: Sun Apr 18, 2021 8:15 am

STRANGE BEHAVIOUR

Post by sedatacemogluone » Sat Jan 22, 2022 6:59 am

Hello,
When I put the beam in model, it gives error. My model below consists of two columns and one beam. Beam is having same fibers as columns. Their size also same. I dont understand whats the problem in this code. it does not simply run.

Thanks for any help.

Sedat Acemoglu



////////CODE
wipe
model BasicBuilder -ndm 3 -ndf 6
node 0 0 0 0 -mass 0 0 0.0 0.0 0.0 0.0
node 1 0 0 5 -mass 4776.88 4776.88 0.0 0.0 0.0 0.0
node 2 10 0 0 -mass 0 0 0.0 0.0 0.0 0.0
node 3 10 0 5 -mass 4776.88 4776.88 0.0 0.0 0.0 0.0
fix 0 1 1 1 1 1 1
fix 2 1 1 1 1 1 1
uniaxialMaterial ReinforcingSteel 1000001000 420000000 504000000 200000000000 1078305519.8973 0.008 0.08
uniaxialMaterial Concrete07 1000002000 -30000000 -0.002 7083631551.92269 1250000 0.000352926317761605 2 2.3 3.86923076923077
uniaxialMaterial Concrete07 1000003000 -26134372.9539653 -0.000711457651321774 7083631551.92269 1250000 0.000352926317761605 2 30 -2.92980199358746
section Fiber 100000 -GJ 12720393247.5672 {
patch quadr 1000003000 20 20 -0.196 0.196 -0.196 -0.196 0.196 -0.196 0.196 0.196
patch quadr 1000002000 2 20 -0.25 0.25 -0.196 0.196 0.196 0.196 0.25 0.25
patch quadr 1000002000 2 20 -0.196 -0.196 -0.25 -0.25 0.25 -0.25 0.196 -0.196
patch quadr 1000002000 20 2 0.25 0.25 0.196 0.196 0.196 -0.196 0.25 -0.25
patch quadr 1000002000 20 2 -0.25 0.25 -0.196 0.196 -0.196 -0.196 -0.25 -0.25
fiber 0.19 0.19 0.000314159265358979 1000001000
fiber 0.19 0.095 0.000314159265358979 1000001000
fiber 0.19 0 0.000314159265358979 1000001000
fiber 0.19 -0.095 0.000314159265358979 1000001000
fiber 0.19 -0.19 0.000314159265358979 1000001000
fiber 0.095 0.19 0.000314159265358979 1000001000
fiber 0.095 -0.19 0.000314159265358979 1000001000
fiber 0 0.19 0.000314159265358979 1000001000
fiber 0 -0.19 0.000314159265358979 1000001000
fiber -0.095 0.19 0.000314159265358979 1000001000
fiber -0.095 -0.19 0.000314159265358979 1000001000
fiber -0.19 0.19 0.000314159265358979 1000001000
fiber -0.19 0.095 0.000314159265358979 1000001000
fiber -0.19 0 0.000314159265358979 1000001000
fiber -0.19 -0.095 0.000314159265358979 1000001000
fiber -0.19 -0.19 0.000314159265358979 1000001000
}

geomTransf PDelta 1 -1 0 0
uniaxialMaterial Elastic 133333 16000000000
section Aggregator 144444 133333 T -section 100000
element dispBeamColumn 1 0 1 5 144444 1
uniaxialMaterial ReinforcingSteel 1000011000 420000000 504000000 200000000000 1078305519.8973 0.008 0.08
uniaxialMaterial Concrete07 1000012000 -30000000 -0.002 7083631551.92269 1250000 0.000352926317761605 2 2.3 3.86923076923077
uniaxialMaterial Concrete07 1000013000 -26134372.9539653 -0.000711457651321774 7083631551.92269 1250000 0.000352926317761605 2 30 -2.92980199358746
section Fiber 100001 -GJ 12720393247.5672 {
patch quadr 1000013000 20 20 -0.196 0.196 -0.196 -0.196 0.196 -0.196 0.196 0.196
patch quadr 1000012000 2 20 -0.25 0.25 -0.196 0.196 0.196 0.196 0.25 0.25
patch quadr 1000012000 2 20 -0.196 -0.196 -0.25 -0.25 0.25 -0.25 0.196 -0.196
patch quadr 1000012000 20 2 0.25 0.25 0.196 0.196 0.196 -0.196 0.25 -0.25
patch quadr 1000012000 20 2 -0.25 0.25 -0.196 0.196 -0.196 -0.196 -0.25 -0.25
fiber 0.19 0.19 0.000314159265358979 1000011000
fiber 0.19 0.095 0.000314159265358979 1000011000
fiber 0.19 0 0.000314159265358979 1000011000
fiber 0.19 -0.095 0.000314159265358979 1000011000
fiber 0.19 -0.19 0.000314159265358979 1000011000
fiber 0.095 0.19 0.000314159265358979 1000011000
fiber 0.095 -0.19 0.000314159265358979 1000011000
fiber 0 0.19 0.000314159265358979 1000011000
fiber 0 -0.19 0.000314159265358979 1000011000
fiber -0.095 0.19 0.000314159265358979 1000011000
fiber -0.095 -0.19 0.000314159265358979 1000011000
fiber -0.19 0.19 0.000314159265358979 1000011000
fiber -0.19 0.095 0.000314159265358979 1000011000
fiber -0.19 0 0.000314159265358979 1000011000
fiber -0.19 -0.095 0.000314159265358979 1000011000
fiber -0.19 -0.19 0.000314159265358979 1000011000
}

geomTransf PDelta 2 -1 0 0
uniaxialMaterial Elastic 133334 16000000000
section Aggregator 144445 133334 T -section 100001
element dispBeamColumn 2 2 3 5 144445 2
uniaxialMaterial ReinforcingSteel 1000021000 420000000 504000000 200000000000 1078305519.8973 0.008 0.08
uniaxialMaterial Concrete07 1000022000 -30000000 -0.002 5477225575.05166 1250000 0.000456435464587638 2 2.3 3.86923076923077
uniaxialMaterial Concrete07 1000023000 -26134372.9539653 -0.000711457651321774 5477225575.05166 1250000 0.000456435464587638 2 30 -2.92980199358746
section Fiber 100002 -GJ 12720393247.5672 {
patch quadr 1000023000 20 20 -0.196 0.196 -0.196 -0.196 0.196 -0.196 0.196 0.196
patch quadr 1000022000 2 20 -0.25 0.25 -0.196 0.196 0.196 0.196 0.25 0.25
patch quadr 1000022000 2 20 -0.196 -0.196 -0.25 -0.25 0.25 -0.25 0.196 -0.196
patch quadr 1000022000 20 2 0.25 0.25 0.196 0.196 0.196 -0.196 0.25 -0.25
patch quadr 1000022000 20 2 -0.25 0.25 -0.196 0.196 -0.196 -0.196 -0.25 -0.25
fiber 0.19 0.19 0.000314159265358979 1000021000
fiber 0.19 0.095 0.000314159265358979 1000021000
fiber 0.19 0 0.000314159265358979 1000021000
fiber 0.19 -0.095 0.000314159265358979 1000021000
fiber 0.19 -0.19 0.000314159265358979 1000021000
fiber 0.095 0.19 0.000314159265358979 1000021000
fiber 0.095 -0.19 0.000314159265358979 1000021000
fiber 0 0.19 0.000314159265358979 1000021000
fiber 0 -0.19 0.000314159265358979 1000021000
fiber -0.095 0.19 0.000314159265358979 1000021000
fiber -0.095 -0.19 0.000314159265358979 1000021000
fiber -0.19 0.19 0.000314159265358979 1000021000
fiber -0.19 0.095 0.000314159265358979 1000021000
fiber -0.19 0 0.000314159265358979 1000021000
fiber -0.19 -0.095 0.000314159265358979 1000021000
fiber -0.19 -0.19 0.000314159265358979 1000021000
}

geomTransf PDelta 3 0 0 1
uniaxialMaterial Elastic 133335 16000000000
section Aggregator 144446 133335 T -section 100002
element dispBeamColumn 3 1 3 5 144446 3
pattern Plain 1 "Linear" {
load 1 0 0 -46861.15875 0 0 0
load 3 0 0 -46861.15875 0 0 0
}
test NormDispIncr 1.0e-6 100
system BandGeneral
numberer RCM
constraints Plain
integrator LoadControl 0.1
algorithm Newton
analysis Static
analyze 10
loadConst -time 0.0
recorder Node -precision 10 -xml "C:\\Users\\90543\\Documents\\NTHA\\yksl\\recs\\RSN4352_UBMARCHE.P_A-NCR270.AT2!nodesOutX.txt" -dof 1 2 3 4 5 6 disp
timeSeries Path 2 -dt 0.005 -filePath ivmeopensees.txt -factor 4.905;
pattern UniformExcitation 2 1 -accel 2;
pattern UniformExcitation 3 2 -accel 2;
source SmartAnalyze.tcl;
wipeAnalysis;
constraints Transformation;
numberer Plain;
system BandGeneral;
integrator Newmark 0.5 0.25;
SmartAnalyzeTransient 0.005 1000 ;
wipe;

mhscott
Posts: 872
Joined: Tue Jul 06, 2004 3:38 pm
Location: Corvallis, Oregon USA
Contact:

Re: STRANGE BEHAVIOUR

Post by mhscott » Sat Jan 22, 2022 2:45 pm

Backpedal to a linear-elastic model where everything works, then re-introduce nonlinearity to find the problem.
https://portwooddigital.com/2020/08/21/ ... -a-sprint/

sedatacemogluone
Posts: 4
Joined: Sun Apr 18, 2021 8:15 am

Re: STRANGE BEHAVIOUR

Post by sedatacemogluone » Sun Jan 23, 2022 12:52 am

Thank you Professor, I later changed number of integration points from 5 to 4 and it worked properly.

Post Reply