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.
Zerolength orientation
Moderators: silvia, selimgunay, Moderators
-
- Posts: 22
- Joined: Thu Jun 08, 2017 2:31 pm
- Location: Babol Noshirvani University of Technology
-
- Posts: 22
- Joined: Thu Jun 08, 2017 2:31 pm
- Location: Babol Noshirvani University of Technology
Re: Zerolength orientation
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!"
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!"
Re: Zerolength orientation
Hello Dears
How must calculate the parameters of orient for zero-length element?
How must calculate the parameters of orient for zero-length element?