Zerolength orientation

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

Moderators: silvia, selimgunay, Moderators

Post Reply
ehsanomranimemphis
Posts: 22
Joined: Thu Jun 08, 2017 2:31 pm
Location: Babol Noshirvani University of Technology

Zerolength orientation

Post by ehsanomranimemphis » Fri Mar 16, 2018 9:14 am

I used zero length element to model abutment. Since the deck in my model is not straight, I have to rotate the direction of zerolength element relative to deck axis to be able to consider local coordinate of abutment. But the results are the same with original direction (original direction means zero length element in global direction ). Is there some one to help me?
Thank you all.

ehsanomranimemphis
Posts: 22
Joined: Thu Jun 08, 2017 2:31 pm
Location: Babol Noshirvani University of Technology

Re: Zerolength orientation

Post by ehsanomranimemphis » Sat Mar 17, 2018 8:08 am

This is my OpenSees code.
wipe;
wipeAnalysis;

file mkdir Data;
model BasicBuilder -ndm 2 -ndf 3;

node 1 0.0 0.0;
node 2 0.0 0.0;


equalDOF 1 2 1 2;

fix 1 1 1 1;


set eleID 1 ;
set K 26.43e6 ;
set as 0.035 ;
set My 3.46e5 ;
set LS 3.6 ;
set LD 3.6;
set LA 3.6 ;
set LK 3.6 ;
set c 1. ;
set th_p 0.06 ;
set th_pc 0.225 ;
set Res 0.0 ;
set th_u 0.3 ;
set D 1. ;

uniaxialMaterial ModIMKPeakOriented $eleID $K $as $as $My [expr -$My] $LS $LD $LA $LK $c $c $c $c $th_p $th_p $th_pc $th_pc $Res $Res $th_u $th_u $D $D

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



recorder Node -file data/node1reac.out -node 1 -dof 3 reaction
recorder Node -file data/node2disp.out -node 2 -dof 3 disp

pattern Plain 20 Linear {
load 2 0.0 0.0 1.0e10;
};

constraints Transformation
numberer RCM
system BandGeneral
test EnergyIncr 1.e-5 200
algorithm KrylovNewton

set disploop [list 8.225 16.45 16.45 16.45 24.675 32.9 32.9 32.9 49.35 65.8 65.8 65.8 82.25 98.7 98.7 98.7 115.15 131.6 131.6 131.6 148.05 164.5 164.5 164.5 180.95 194.4 194.4 194.4]

for {set i 1} {$i<=[llength $disploop]} {incr i 1} {
set dispincr [expr 1.0*pow((-1.0),[expr $i-1])*[lindex $disploop [expr $i-1]]/100000]
integrator DisplacementControl 2 1 $dispincr

analysis Static
analyze 100
}
puts "Run lateral completed!"

nilz
Posts: 34
Joined: Sat Nov 12, 2016 7:03 am
Location: Ferdowsi university of mashhad

Re: Zerolength orientation

Post by nilz » Mon Mar 19, 2018 1:32 am

Hello Dears
How must calculate the parameters of orient for zero-length element?

Post Reply