Here I am using same elements in X and Z directions, just with different transformations so as to orient members have stronger axis normal to the plane of moment resisting frame. Thus the stiffness of frames in X and Z directions should be of the same order.
The bug is that I am getting very very less deformations in Z direction as compared to X direction (about 500 times less when same load is applied in X and Z direction at the top floor). Is there any problem with fix command for center of mass nodes. Or is there any problem with the node slavings.
The geometry is Los Angeles Building with pre-Northridge design. Its figure and table for sections is attached herewith.
Input code for geometry is-
Code: Select all
#-----------------------------------------------------------------------------------------------------
wipe;
source DisplayPlane.tcl; # PROCEDURE FOR DISPLAYING A PLANE IN MODEL
source DisplayModel3D.tcl; # PROCEDURE FOR DISPLAYING 3D PERSPECTIVES
source DisplayEigen.tcl; # PROCEDURE FOR DISPLAYING MODE SHAPES
source Wsection.tcl; # PROCEDURE TO DEFINE FIBER W SECTION
puts "\n\nCreating Geometry";
#
model BasicBuilder -ndm 3 -ndf 6;
#-----------------------------------------------------------------------------------------------------
puts "Creating Nodes...";
#-----------------
# FOUNDATION (BOTTOM NODES)
# tag x y z
node 111 0 0 0;
node 112 360 0 0;
node 113 720 0 0;
node 114 1080 0 0;
node 115 1440 0 0;
node 116 1800 0 0;
node 117 2160 0 0;
#
node 121 0 0 360;
node 127 2160 0 360;
#
node 131 0 0 720;
node 137 2160 0 720;
#
node 141 0 0 1080;
node 147 2160 0 1080;
#
node 151 0 0 1440;
node 152 360 0 1440;
node 153 720 0 1440
node 154 1080 0 1440;
node 155 1440 0 1440;
node 156 1800 0 1440;
node 157 2160 0 1440;
#-----------------
# FIRST STOREY
node 211 0 156 0;
node 212 360 156 0;
node 213 720 156 0;
node 214 1080 156 0;
node 215 1440 156 0;
node 216 1800 156 0;
node 217 2160 156 0;
node 511 0 156 0;
node 512 360 156 0;
node 513 1440 156 0;
node 514 1800 156 0;
node 515 1800 156 0;
node 516 2160 156 0;
#
node 221 0 156 360;
node 227 2160 156 360;
#
node 231 0 156 720;
node 237 2160 156 720;
#
node 241 0 156 1080;
node 247 2160 156 1080;
node 541 0 156 1080;
node 547 2160 156 1080;
#
node 251 0 156 1440;
node 252 360 156 1440;
node 253 720 156 1440;
node 254 1080 156 1440;
node 255 1440 156 1440;
node 256 1800 156 1440;
node 257 2160 156 1440;
node 551 0 156 1440;
node 552 360 156 1440;
node 553 1440 156 1440;
node 554 1800 156 1440;
node 555 1800 156 1440;
node 556 2160 156 1440;
node 557 0 156 1440;
node 558 2160 156 1440;
#-----------------
# SECOND STOREY
node 311 0 312 0;
node 312 360 312 0;
node 313 720 312 0;
node 314 1080 312 0;
node 315 1440 312 0;
node 316 1800 312 0;
node 317 2160 312 0;
node 611 0 312 0;
node 612 360 312 0;
node 613 1440 312 0;
node 614 1800 312 0;
node 615 1800 312 0;
node 616 2160 312 0;
#
node 321 0 312 360;
node 327 2160 312 360;
#
node 331 0 312 720;
node 337 2160 312 720;
#
node 341 0 312 1080;
node 347 2160 312 1080;
node 641 0 312 1080;
node 647 2160 312 1080;
#
node 351 0 312 1440;
node 352 360 312 1440;
node 353 720 312 1440;
node 354 1080 312 1440;
node 355 1440 312 1440;
node 356 1800 312 1440;
node 357 2160 312 1440;
node 651 0 312 1440;
node 652 360 312 1440;
node 653 1440 312 1440;
node 654 1800 312 1440;
node 655 1800 312 1440;
node 656 2160 312 1440;
node 657 0 312 1440;
node 658 2160 312 1440;
#-----------------
# THIRD STOREY
node 411 0 468 0;
node 412 360 468 0;
node 413 720 468 0;
node 414 1080 468 0;
node 415 1440 468 0;
node 416 1800 468 0;
node 417 2160 468 0;
node 711 0 468 0;
node 712 360 468 0;
node 713 1440 468 0;
node 714 1800 468 0;
node 715 1800 468 0;
node 716 2160 468 0;
#
node 421 0 468 360;
node 427 2160 468 360;
#
node 431 0 468 720;
node 437 2160 468 720;
#
node 441 0 468 1080;
node 447 2160 468 1080;
node 741 0 468 1080;
node 747 2160 468 1080;
#
node 451 0 468 1440;
node 452 360 468 1440;
node 453 720 468 1440;
node 454 1080 468 1440;
node 455 1440 468 1440;
node 456 1800 468 1440;
node 457 2160 468 1440;
node 751 0 468 1440;
node 752 360 468 1440;
node 753 1440 468 1440;
node 754 1800 468 1440;
node 755 1800 468 1440;
node 756 2160 468 1440;
node 757 0 468 1440;
node 758 2160 468 1440;
#------------------------
# EXTRA NODES FOR CONCENTRATED MASS AT GIVEN STOREY
set Xecc 0; # % ECCENTRICITY TO BE INPUTTED BY USER
set Zecc 0;
set Xcm [expr 1080.0 + 2160.0*$Xecc/100.0]; # CENTRE OF MASS COORDINATES
set Zcm [expr 720.0 + 1440.0*$Zecc/100.0];
#set Xcm 0;
#set Zcm 0;
node 200 $Xcm 156 $Zcm;
node 300 $Xcm 312 $Zcm;
node 400 $Xcm 468 $Zcm;
#
#-----------------------------------------------------------------------------------------------------
# NODAL MASSES
puts "Assigning Masses...";
mass 200 5.460 1e-6 5.460 1e-6 1e-6 1e-6;
mass 300 5.460 1e-6 5.460 1e-6 1e-6 1e-6;
mass 400 5.908 1e-6 5.908 1e-6 1e-6 1e-6;
#
#-----------------------------------------------------------------------------------------------------
puts "Assigning Rigid Diaphragm...";
# DX DY DZ RX RY RZ
fix 200 0 1 0 1 0 1;
fix 300 0 1 0 1 0 1;
fix 400 0 1 0 1 0 1;
# Perp Mst Slv
rigidDiaphragm 2 200 211 212 213 214 215 216 217 221 227 231 237 241 247 251 252 253 254 255 256 257;
rigidDiaphragm 2 300 311 312 313 314 315 316 317 321 327 331 337 341 347 351 352 353 354 355 356 357;
rigidDiaphragm 2 400 411 412 413 414 415 416 417 421 427 431 437 441 447 451 452 453 454 455 456 457;
#
#-----------------------------------------------------------------------------------------------------
# SETTING NODE SLAVING FOR THE PINNED BEAM
puts "Setting Node Slaving..."
#----------------
# FIRST STOREY
equalDOF 211 511 1 2 3 4 5;
equalDOF 212 512 1 2 3 4 5;
equalDOF 215 513 1 2 3 4 5;
equalDOF 216 514 1 2 3 4 5;
equalDOF 216 515 1 2 3 4 5;
equalDOF 217 516 1 2 3 4 5;
#
equalDOF 241 541 1 2 3 5 6;
equalDOF 247 547 1 2 3 5 6;
#
equalDOF 251 551 1 2 3 4 5;
equalDOF 252 552 1 2 3 4 5;
equalDOF 255 553 1 2 3 4 5;
equalDOF 256 554 1 2 3 4 5;
equalDOF 256 555 1 2 3 4 5;
equalDOF 257 556 1 2 3 4 5;
#
equalDOF 251 557 1 2 3 5 6;
equalDOF 257 558 1 2 3 5 6;
#----------------
# SECOND STOREY
equalDOF 311 611 1 2 3 4 5;
equalDOF 312 612 1 2 3 4 5;
equalDOF 315 613 1 2 3 4 5;
equalDOF 316 614 1 2 3 4 5;
equalDOF 316 615 1 2 3 4 5;
equalDOF 317 616 1 2 3 4 5;
#
equalDOF 341 641 1 2 3 5 6;
equalDOF 347 647 1 2 3 5 6;
#
equalDOF 351 651 1 2 3 4 5;
equalDOF 352 652 1 2 3 4 5;
equalDOF 355 653 1 2 3 4 5;
equalDOF 356 654 1 2 3 4 5;
equalDOF 356 655 1 2 3 4 5;
equalDOF 357 656 1 2 3 4 5;
#
equalDOF 351 657 1 2 3 5 6;
equalDOF 357 658 1 2 3 5 6;
#----------------
# THIRD STOREY
equalDOF 411 711 1 2 3 4 5;
equalDOF 412 712 1 2 3 4 5;
equalDOF 415 713 1 2 3 4 5;
equalDOF 416 714 1 2 3 4 5;
equalDOF 416 715 1 2 3 4 5;
equalDOF 417 716 1 2 3 4 5;
#
equalDOF 441 741 1 2 3 5 6;
equalDOF 447 747 1 2 3 5 6;
#
equalDOF 451 751 1 2 3 4 5;
equalDOF 452 752 1 2 3 4 5;
equalDOF 455 753 1 2 3 4 5;
equalDOF 456 754 1 2 3 4 5;
equalDOF 456 755 1 2 3 4 5;
equalDOF 457 756 1 2 3 4 5;
#
equalDOF 451 757 1 2 3 5 6;
equalDOF 457 758 1 2 3 5 6;
#
#-----------------------------------------------------------------------------------------------------
puts "Assigning Boundary Conditions...";
#fixY 0.0 1 1 1 1 1 1; # ALL BOTTOM NODES ARE FIXED
# tag x y z
fix 111 1 1 1 1 1 1;
fix 112 1 1 1 1 1 1;
fix 113 1 1 1 1 1 1;
fix 114 1 1 1 1 1 1;
fix 115 1 1 1 1 1 1;
fix 116 1 1 1 1 1 1;
fix 117 1 1 1 1 1 1;
#
fix 121 1 1 1 1 1 1;
fix 127 1 1 1 1 1 1;
#
fix 131 1 1 1 1 1 1;
fix 137 1 1 1 1 1 1;
#
fix 141 1 1 1 1 1 1;
fix 147 1 1 1 1 1 1;
#
fix 151 1 1 1 1 1 1;
fix 152 1 1 1 1 1 1;
fix 153 1 1 1 1 1 1;
fix 154 1 1 1 1 1 1;
fix 155 1 1 1 1 1 1;
fix 156 1 1 1 1 1 1;
#
#-----------------------------------------------------------------------------------------------------
puts "Defining Materials...";
set matTag 1;
set E 29000.0;
set epsyP 0.002;
puts " INFO: Yield stress is [expr $E*$epsyP]";
uniaxialMaterial ElasticPP $matTag $E $epsyP
set nu 0.3;
set Gs [expr $E/2.0/[expr 1+$nu]]; # FOR TORSIONAL STIFFNESS
#
#-----------------------------------------------------------------------------------------------------
puts "Defining Sections...";
# INPUT PARAMETERS
# secID - section ID number
# matID - material ID number
# d = nominal depth
# tw = web thickness
# bf = flange width
# tf = flange thickness
# nfdw = number of fibers along web depth
# nftw = number of fibers along web thickness
# nfbf = number of fibers along flange width
# nftf = number of fibers along flange thickness
#
# SETTING THE SECTION ID's
# FINAL SECTIONS
# MOMENT RESISTING MEMBERS
set ColExtMom 1;
set ColIntMom 2;
set Bem1stMom 3;
set Bem2ndMom 4;
set Bem3rdMom 5;
# GRAVITY FRAME MEMBERS
set ColAllGrv 6;
set Bem1stGrv 7;
set Bem2ndGrv 8;
set Bem3rdGrv 9;
# FIBER SECTIONS
# MOMENT RESISTING MEMBERS
set ColExtMomF 10;
set ColIntMomF 20;
set Bem1stMomF 30;
set Bem2ndMomF 40;
set Bem3rdMomF 50;
# GRAVITY FRAME MEMBERS
set ColAllGrvF 60;
set Bem1stGrvF 70;
set Bem2ndGrvF 80;
set Bem3rdGrvF 90;
# TORSIONAL PROPERTIES
set TorsionSec1 101;
set TorsionSec2 102;
set TorsionSec3 103;
set TorsionSec4 104;
set TorsionSec5 105;
set TorsionSec6 106;
set TorsionSec7 107;
set TorsionSec8 108;
set TorsionSec9 109;
# SOME COMMON VALUES
set matTag 1;
set nfdw 16;
set nftw 2;
set nfbf 16;
set nftf 4;
#-------------------
# W14X257
set d 16.38;
set tw 1.175;
set bf 15.995;
set tf 1.890;
Wsection $ColExtMomF $matTag $d $bf $tf $tw $nfdw $nftw $nfbf $nftf;
#
set GJ [expr 4690*$Gs];
uniaxialMaterial Elastic $TorsionSec1 $GJ;
section Aggregator $ColExtMom $TorsionSec1 T -section $ColExtMomF;
#-------------------
# W14X311
set d 17.12;
set tw 1.410;
set bf 16.230;
set tf 2.260;
Wsection $ColIntMomF $matTag $d $bf $tf $tw $nfdw $nftw $nfbf $nftf;
#
set GJ [expr 5940*$Gs];
uniaxialMaterial Elastic $TorsionSec2 $GJ;
section Aggregator $ColIntMom $TorsionSec2 T -section $ColIntMomF;
#-------------------
# W33X118
set d 32.86;
set tw 0.550;
set bf 11.480;
set tf 0.740;
Wsection $Bem1stMomF $matTag $d $bf $tf $tw $nfdw $nftw $nfbf $nftf;
#
set GJ [expr 6087*$Gs];
uniaxialMaterial Elastic $TorsionSec3 $GJ;
section Aggregator $Bem1stMom $TorsionSec3 T -section $Bem1stMomF;
#-------------------
# W30X116
set d 30.01;
set tw 0.565;
set bf 10.495;
set tf 0.850;
Wsection $Bem2ndMomF $matTag $d $bf $tf $tw $nfdw $nftw $nfbf $nftf;
#
set GJ [expr 5094*$Gs];
uniaxialMaterial Elastic $TorsionSec4 $GJ;
section Aggregator $Bem2ndMom $TorsionSec4 T -section $Bem2ndMomF;
#-------------------
# W24X68
set d 23.73;
set tw 0.415;
set bf 8.965;
set tf 0.585;
Wsection $Bem3rdMomF $matTag $d $bf $tf $tw $nfdw $nftw $nfbf $nftf;
#
set GJ [expr 1900*$Gs];
uniaxialMaterial Elastic $TorsionSec5 $GJ;
section Aggregator $Bem3rdMom $TorsionSec5 T -section $Bem3rdMomF;
#-------------------
# W14X68
set d 14.04;
set tw 0.415;
set bf 10.035;
set tf 0.720;
Wsection $ColAllGrvF $matTag $d $bf $tf $tw $nfdw $nftw $nfbf $nftf;
#
set GJ [expr 841*$Gs];
uniaxialMaterial Elastic $TorsionSec6 $GJ;
section Aggregator $ColAllGrv $TorsionSec6 T -section $ColAllGrvF;
#-------------------
# W18X35
set d 17.70;
set tw 0.300;
set bf 6.000;
set tf 0.425;
Wsection $Bem1stGrvF $matTag $d $bf $tf $tw $nfdw $nftw $nfbf $nftf;
#
set GJ [expr 525*$Gs];
uniaxialMaterial Elastic $TorsionSec7 $GJ;
section Aggregator $Bem1stGrv $TorsionSec7 T -section $Bem1stGrvF;
#-------------------
# W18X35
set d 17.70;
set tw 0.300;
set bf 6.000;
set tf 0.425;
Wsection $Bem2ndGrvF $matTag $d $bf $tf $tw $nfdw $nftw $nfbf $nftf;
#
set GJ [expr 525*$Gs];
uniaxialMaterial Elastic $TorsionSec8 $GJ;
section Aggregator $Bem2ndGrv $TorsionSec8 T -section $Bem2ndGrvF;
#-------------------
# W16X26
set d 15.69;
set tw 0.250;
set bf 5.500;
set tf 0.345;
Wsection $Bem3rdGrvF $matTag $d $bf $tf $tw $nfdw $nftw $nfbf $nftf;
#
set GJ [expr 311*$Gs];
uniaxialMaterial Elastic $TorsionSec9 $GJ;
section Aggregator $Bem3rdGrv $TorsionSec9 T -section $Bem3rdGrvF;
#
#-----------------------------------------------------------------------------------------------------
# DEFINING TRANSFORMATION FOR ELEMENTS
puts "Defining Transformations...";
set TransColmNS 1;
set TransColmEW 2;
set TransBeamNS 3;
set TransBeamEW 4;
set ColmTransType Linear; # OPTIONS, LINEAR, PDELTA, COROTATIONAL FOR COLUMN
set BeamTransType Linear; # OPTIONS, LINEAR, PDELTA, COROTATIONAL FOR BEAM
geomTransf $ColmTransType $TransColmNS 1.0 0.0 0.0;
geomTransf $ColmTransType $TransColmEW 0.0 0.0 1.0;
geomTransf $BeamTransType $TransBeamNS 1.0 0.0 0.0;
geomTransf $BeamTransType $TransBeamEW 0.0 0.0 1.0;
#
if {$ColmTransType == "Linear"} {;
puts " INFO: P-Delta is NOT INCLUDED for columns";
} elseif {$ColmTransType == "PDelta"} {;
puts " INFO: P-Delta is INCLUDED for columns";
} elseif {$ColmTransType == "Corotational"} {;
puts " INFO: COROTATIONAL transformation for columns";
} else {;
puts "\aInvalid Transformation option... Exiting";
return -1;
}
#
#-----------------------------------------------------------------------------------------------------
# DEFINE ELEMENTS
puts "Defining Elements..."
set np 10;
puts " INFO: Number of integration points used is $np";
#_____________________________
# FIRST STOREY COLUMNS
#
# $eleTag $iNode $jNode $numIntgrPts $secTag $transfTag
element nonlinearBeamColumn 111 111 211 $np $ColExtMom $TransColmNS;
element nonlinearBeamColumn 112 112 212 $np $ColExtMom $TransColmEW;
element nonlinearBeamColumn 113 113 213 $np $ColIntMom $TransColmEW;
element nonlinearBeamColumn 114 114 214 $np $ColIntMom $TransColmEW;
element nonlinearBeamColumn 115 115 215 $np $ColExtMom $TransColmEW;
element nonlinearBeamColumn 116 116 216 $np $ColAllGrv $TransColmEW;
element nonlinearBeamColumn 117 117 217 $np $ColExtMom $TransColmNS;
#
element nonlinearBeamColumn 121 121 221 $np $ColIntMom $TransColmNS;
element nonlinearBeamColumn 127 127 227 $np $ColIntMom $TransColmNS;
#
element nonlinearBeamColumn 131 131 231 $np $ColIntMom $TransColmNS;
element nonlinearBeamColumn 137 137 237 $np $ColIntMom $TransColmNS;
#
element nonlinearBeamColumn 141 141 241 $np $ColExtMom $TransColmNS;
element nonlinearBeamColumn 147 147 247 $np $ColExtMom $TransColmNS;
#
element nonlinearBeamColumn 151 151 251 $np $ColAllGrv $TransColmNS;
element nonlinearBeamColumn 152 152 252 $np $ColExtMom $TransColmEW;
element nonlinearBeamColumn 153 153 253 $np $ColIntMom $TransColmEW;
element nonlinearBeamColumn 154 154 254 $np $ColIntMom $TransColmEW;
element nonlinearBeamColumn 155 155 255 $np $ColExtMom $TransColmEW;
element nonlinearBeamColumn 156 156 256 $np $ColAllGrv $TransColmEW;
element nonlinearBeamColumn 157 157 257 $np $ColAllGrv $TransColmNS;
#_____________________________
# SECOND STOREY COLUMNS
#
element nonlinearBeamColumn 211 211 311 $np $ColExtMom $TransColmNS;
element nonlinearBeamColumn 212 212 312 $np $ColExtMom $TransColmEW;
element nonlinearBeamColumn 213 213 313 $np $ColIntMom $TransColmEW;
element nonlinearBeamColumn 214 214 314 $np $ColIntMom $TransColmEW;
element nonlinearBeamColumn 215 215 315 $np $ColExtMom $TransColmEW;
element nonlinearBeamColumn 216 216 316 $np $ColAllGrv $TransColmEW;
element nonlinearBeamColumn 217 217 317 $np $ColExtMom $TransColmNS;
#
element nonlinearBeamColumn 221 221 321 $np $ColIntMom $TransColmNS;
element nonlinearBeamColumn 227 227 327 $np $ColIntMom $TransColmNS;
#
element nonlinearBeamColumn 231 231 331 $np $ColIntMom $TransColmNS;
element nonlinearBeamColumn 237 237 337 $np $ColIntMom $TransColmNS;
#
element nonlinearBeamColumn 241 241 341 $np $ColExtMom $TransColmNS;
element nonlinearBeamColumn 247 247 347 $np $ColExtMom $TransColmNS;
#
element nonlinearBeamColumn 251 251 351 $np $ColAllGrv $TransColmNS;
element nonlinearBeamColumn 252 252 352 $np $ColExtMom $TransColmEW;
element nonlinearBeamColumn 253 253 353 $np $ColIntMom $TransColmEW;
element nonlinearBeamColumn 254 254 354 $np $ColIntMom $TransColmEW;
element nonlinearBeamColumn 255 255 355 $np $ColExtMom $TransColmEW;
element nonlinearBeamColumn 256 256 356 $np $ColAllGrv $TransColmEW;
element nonlinearBeamColumn 257 257 357 $np $ColAllGrv $TransColmNS;
#__________________________________
# THIRD STOREY COLUMNS
#
element nonlinearBeamColumn 311 311 411 $np $ColExtMom $TransColmNS;
element nonlinearBeamColumn 312 312 412 $np $ColExtMom $TransColmEW;
element nonlinearBeamColumn 313 313 413 $np $ColIntMom $TransColmEW;
element nonlinearBeamColumn 314 314 414 $np $ColIntMom $TransColmEW;
element nonlinearBeamColumn 315 315 415 $np $ColExtMom $TransColmEW;
element nonlinearBeamColumn 316 316 416 $np $ColAllGrv $TransColmEW;
element nonlinearBeamColumn 317 317 417 $np $ColExtMom $TransColmNS;
#
element nonlinearBeamColumn 321 321 421 $np $ColIntMom $TransColmNS;
element nonlinearBeamColumn 327 327 427 $np $ColIntMom $TransColmNS;
#
element nonlinearBeamColumn 331 331 431 $np $ColIntMom $TransColmNS;
element nonlinearBeamColumn 337 337 437 $np $ColIntMom $TransColmNS;
#
element nonlinearBeamColumn 341 341 441 $np $ColExtMom $TransColmNS;
element nonlinearBeamColumn 347 347 447 $np $ColExtMom $TransColmNS;
#
element nonlinearBeamColumn 351 351 451 $np $ColAllGrv $TransColmNS;
element nonlinearBeamColumn 352 352 452 $np $ColExtMom $TransColmEW;
element nonlinearBeamColumn 353 353 453 $np $ColIntMom $TransColmEW;
element nonlinearBeamColumn 354 354 454 $np $ColIntMom $TransColmEW;
element nonlinearBeamColumn 355 355 455 $np $ColExtMom $TransColmEW;
element nonlinearBeamColumn 356 356 456 $np $ColAllGrv $TransColmEW;
element nonlinearBeamColumn 357 357 457 $np $ColAllGrv $TransColmNS;
#____________________________
# FIRST STOREY BEAMS
#
element nonlinearBeamColumn 2001 211 221 $np $Bem1stMom $TransBeamNS;
element nonlinearBeamColumn 2002 221 231 $np $Bem1stMom $TransBeamNS;
element nonlinearBeamColumn 2003 231 241 $np $Bem1stMom $TransBeamNS;
element nonlinearBeamColumn 2004 541 557 $np $Bem1stGrv $TransBeamNS;
element nonlinearBeamColumn 2005 551 552 $np $Bem1stGrv $TransBeamEW;
element nonlinearBeamColumn 2006 252 253 $np $Bem1stMom $TransBeamEW;
element nonlinearBeamColumn 2007 253 254 $np $Bem1stMom $TransBeamEW;
element nonlinearBeamColumn 2008 254 255 $np $Bem1stMom $TransBeamEW;
element nonlinearBeamColumn 2009 553 554 $np $Bem1stGrv $TransBeamEW;
element nonlinearBeamColumn 2010 555 556 $np $Bem1stGrv $TransBeamEW;
element nonlinearBeamColumn 2011 558 547 $np $Bem1stGrv $TransBeamNS;
element nonlinearBeamColumn 2012 247 237 $np $Bem1stMom $TransBeamNS;
element nonlinearBeamColumn 2013 237 227 $np $Bem1stMom $TransBeamNS;
element nonlinearBeamColumn 2014 227 217 $np $Bem1stMom $TransBeamNS;
element nonlinearBeamColumn 2015 516 515 $np $Bem1stGrv $TransBeamEW;
element nonlinearBeamColumn 2016 514 513 $np $Bem1stGrv $TransBeamEW;
element nonlinearBeamColumn 2017 215 214 $np $Bem1stMom $TransBeamEW;
element nonlinearBeamColumn 2018 214 213 $np $Bem1stMom $TransBeamEW;
element nonlinearBeamColumn 2019 213 212 $np $Bem1stMom $TransBeamEW;
element nonlinearBeamColumn 2020 512 511 $np $Bem1stGrv $TransBeamEW;
#____________________________
# SECOND STOREY BEAMS
#
element nonlinearBeamColumn 3001 311 321 $np $Bem2ndMom $TransBeamNS;
element nonlinearBeamColumn 3002 321 331 $np $Bem2ndMom $TransBeamNS;
element nonlinearBeamColumn 3003 331 341 $np $Bem2ndMom $TransBeamNS;
element nonlinearBeamColumn 3004 641 657 $np $Bem2ndGrv $TransBeamNS;
element nonlinearBeamColumn 3005 651 652 $np $Bem2ndGrv $TransBeamEW;
element nonlinearBeamColumn 3006 352 353 $np $Bem2ndMom $TransBeamEW;
element nonlinearBeamColumn 3007 353 354 $np $Bem2ndMom $TransBeamEW;
element nonlinearBeamColumn 3008 354 355 $np $Bem2ndMom $TransBeamEW;
element nonlinearBeamColumn 3009 653 654 $np $Bem2ndGrv $TransBeamEW;
element nonlinearBeamColumn 3010 655 656 $np $Bem2ndGrv $TransBeamEW;
element nonlinearBeamColumn 3011 658 647 $np $Bem2ndGrv $TransBeamNS;
element nonlinearBeamColumn 3012 347 337 $np $Bem2ndMom $TransBeamNS;
element nonlinearBeamColumn 3013 337 327 $np $Bem2ndMom $TransBeamNS;
element nonlinearBeamColumn 3014 327 317 $np $Bem2ndMom $TransBeamNS;
element nonlinearBeamColumn 3015 616 615 $np $Bem2ndGrv $TransBeamEW;
element nonlinearBeamColumn 3016 614 613 $np $Bem2ndGrv $TransBeamEW;
element nonlinearBeamColumn 3017 315 314 $np $Bem2ndMom $TransBeamEW;
element nonlinearBeamColumn 3018 314 313 $np $Bem2ndMom $TransBeamEW;
element nonlinearBeamColumn 3019 313 312 $np $Bem2ndMom $TransBeamEW;
element nonlinearBeamColumn 3020 612 611 $np $Bem2ndGrv $TransBeamEW;
#____________________________
# THIRD STOREY BEAMS
#
element nonlinearBeamColumn 4001 411 421 $np $Bem3rdMom $TransBeamNS;
element nonlinearBeamColumn 4002 421 431 $np $Bem3rdMom $TransBeamNS;
element nonlinearBeamColumn 4003 431 441 $np $Bem3rdMom $TransBeamNS;
element nonlinearBeamColumn 4004 741 757 $np $Bem3rdGrv $TransBeamNS;
element nonlinearBeamColumn 4005 751 752 $np $Bem3rdGrv $TransBeamEW;
element nonlinearBeamColumn 4006 452 453 $np $Bem3rdMom $TransBeamEW;
element nonlinearBeamColumn 4007 453 454 $np $Bem3rdMom $TransBeamEW;
element nonlinearBeamColumn 4008 454 455 $np $Bem3rdMom $TransBeamEW;
element nonlinearBeamColumn 4009 753 754 $np $Bem3rdGrv $TransBeamEW;
element nonlinearBeamColumn 4010 755 756 $np $Bem3rdGrv $TransBeamEW;
element nonlinearBeamColumn 4011 758 747 $np $Bem3rdGrv $TransBeamNS;
element nonlinearBeamColumn 4012 447 437 $np $Bem3rdMom $TransBeamNS;
element nonlinearBeamColumn 4013 437 427 $np $Bem3rdMom $TransBeamNS;
element nonlinearBeamColumn 4014 427 417 $np $Bem3rdMom $TransBeamNS;
element nonlinearBeamColumn 4015 716 715 $np $Bem3rdGrv $TransBeamEW;
element nonlinearBeamColumn 4016 714 713 $np $Bem3rdGrv $TransBeamEW;
element nonlinearBeamColumn 4017 415 414 $np $Bem3rdMom $TransBeamEW;
element nonlinearBeamColumn 4018 414 413 $np $Bem3rdMom $TransBeamEW;
element nonlinearBeamColumn 4019 413 412 $np $Bem3rdMom $TransBeamEW;
element nonlinearBeamColumn 4020 712 711 $np $Bem3rdGrv $TransBeamEW;
#
#-----------------------------------------------------------------------------------------------------
puts "Defining Regions...";
#
set RegColm1st 1;
set RegColm2nd 2;
set RegColm3rd 3;
set RegBeam1st 6;
set RegBeam2nd 7;
set RegBeam3rd 8;
#
set RegNode0th 10;
set RegNode1st 11;
set RegNode2nd 12;
set RegNode3rd 13;
#
region $RegColm1st -eleRange 111 157 ;
region $RegColm2nd -eleRange 211 257 ;
region $RegColm3rd -eleRange 311 357 ;
region $RegBeam1st -eleRange 2001 2020;
region $RegBeam2nd -eleRange 3001 3020;
region $RegBeam3rd -eleRange 4001 4020;
#
region $RegNode0th -node 111 112 113 114 115 116 117 121 127 131 137 141 147 151 152 153 154 155 156 157;
region $RegNode1st -node 211 212 213 214 215 216 217 221 227 231 237 241 247 251 252 253 254 255 256 257 200;
region $RegNode2nd -node 311 312 313 314 315 316 317 321 327 331 337 341 347 351 352 353 354 355 356 357 300;
region $RegNode3rd -node 411 412 413 414 415 416 417 421 427 431 437 441 447 451 452 453 454 455 456 457 400;
#
#-----------------------------------------------------------------------------------------------------
puts "COMPLETE GEOMETRY BUILT\n";
Code: Select all
#-----------------------------------------------------------------------------------------------------
# DEFINING GRAVITY LOADS
puts "Defining Gravity Loads...";
pattern Plain 1 "Linear -factor 1.0" {
load 400 100 0 100 0 0 0;
}
#
#-----------------------------------------------------------------------------------------------------
#DEFINING PARAMETERS FOR ANALYSIS
set Tol 1.0e-8;
#
constraints Transformation;
numberer RCM;
system UmfPack;
test EnergyIncr $Tol 10;
algorithm Newton;
set NstepGravity 10;
set DGravity [expr 1.0/$NstepGravity];
integrator LoadControl $DGravity;
#integrator DisplacementControl 200 3 0.1;
analysis Static;
#
#-----------------------------------------------------------------------------------------------------
#CARRYING OUT ANALYSIS
puts "Carrying out Gravity Analysis..."
puts " INFO: Gravity is applied in $NstepGravity steps";
analyze $NstepGravity;
#analyze 60;
#
loadConst -time 0.0;
#
puts "GRAVITY ANALYSIS DONE\n";
What can be the reason for very high stiffness in Z direction?
Thanks a lot for your help
Thanks