I am getting some problems running the eigenvalue analysis when using Rigid Diaphragms with equalDOF to pin the beams at the ends.
Error Message:
ProfileSPDLinDirectSolver::solve() - aii < 0 (i, aii): (4, 0)
ProfileSPDLinDirectSolver::solve() - aii < 0 (i, aii): (9, 0)
ProfileSPDLinDirectSolver::solve() - aii < 0 (i, aii): (15, 0)
ProfileSPDLinDirectSolver::solve() - aii < 0 (i, aii): (22, 0)
T1 = 1.5161886558381732 s
T2 = 1.5161886558381732 s
T3 = 1.5161886558381732 s
The actual period should be around 0.281 seconds.
My script is below:
set startTime [clock clicks -milliseconds]
############################################################################
# File Setup
############################################################################
puts "Importing TCL File Section: File Setup"
wipe all;
model BasicBuilder -ndm 3 -ndf 6;
logFile log.dat;
############################################################################
# Define Geometric Transformations
############################################################################
puts "Importing TCL File Section: Define Geometric Transformations"
geomTransf Linear 1 0 1 0
geomTransf Linear 2 0 0 1
geomTransf Linear 3 0 0 1
############################################################################
# Defining Nodes
############################################################################
puts "Importing TCL File Section: Defining Nodes"
node 1000 0.000000 0.000000 0.000000;
node 1200 360.000000 0.000000 0.000000;
node 1020 0.000000 360.000000 0.000000;
node 1220 360.000000 360.000000 0.000000;
node 1001 0.000000 0.000000 144.000000;
node 1201 360.000000 0.000000 144.000000;
node 1021 0.000000 360.000000 144.000000;
node 1221 360.000000 360.000000 144.000000;
node 2001021 0.000000 20.000000 144.000000;
node 2001029 0.000000 340.000000 144.000000;
node 2001022 0.000000 20.000000 144.000000;
node 2001028 0.000000 340.000000 144.000000;
node 2201221 360.000000 20.000000 144.000000;
node 2201229 360.000000 340.000000 144.000000;
node 2201222 360.000000 20.000000 144.000000;
node 2201228 360.000000 340.000000 144.000000;
node 2001201 20.000000 0.000000 144.000000;
node 2001209 340.000000 0.000000 144.000000;
node 2001202 20.000000 0.000000 144.000000;
node 2001208 340.000000 0.000000 144.000000;
node 2021221 20.000000 360.000000 144.000000;
node 2021229 340.000000 360.000000 144.000000;
node 2021222 20.000000 360.000000 144.000000;
node 2021228 340.000000 360.000000 144.000000;
############################################################################
# Defining Diaphragms
############################################################################
puts "Importing TCL File Section: Defining Diaphragms"
############################################################################
# Defining Fixity
############################################################################
puts "Importing TCL File Section: Defining Fixity"
equalDOF 2021229 2021228 1 2 3
fix 1200 1 1 1 1 1 1;
equalDOF 2001201 2001202 1 2 3
fix 1220 1 1 1 1 1 1;
equalDOF 2201221 2201222 1 2 3
equalDOF 2001209 2001208 1 2 3
fix 1000 1 1 1 1 1 1;
fix 1020 1 1 1 1 1 1;
equalDOF 2001021 2001022 1 2 3
equalDOF 2021221 2021222 1 2 3
equalDOF 2201229 2201228 1 2 3
equalDOF 2001029 2001028 1 2 3
############################################################################
# Defining Masses
############################################################################
puts "Importing TCL File Section: Defining Masses"
mass 1201 0.058230 0.058230 0.000000 0.000000 0.000000 0.000000
mass 1221 0.058230 0.058230 0.000000 0.000000 0.000000 0.000000
mass 1001 0.058230 0.058230 0.000000 0.000000 0.000000 0.000000
mass 1021 0.058230 0.058230 0.000000 0.000000 0.000000 0.000000
############################################################################
# Writing Elements
############################################################################
puts "Importing TCL File Section: Writing Elements"
element elasticBeamColumn 2000001 1000 1001 1050.000000 29000.000000 11500.000000 1000.000000 1000.000000 1000.000000 1; #Plot_Name: Columns - Gravity
element elasticBeamColumn 2020021 1020 1021 1050.000000 29000.000000 11500.000000 1000.000000 1000.000000 1000.000000 1; #Plot_Name: Columns - Gravity
element elasticBeamColumn 2200201 1200 1201 1050.000000 29000.000000 11500.000000 1000.000000 1000.000000 1000.000000 1; #Plot_Name: Columns - Gravity
element elasticBeamColumn 2220221 1220 1221 1050.000000 29000.000000 11500.000000 1000.000000 1000.000000 1000.000000 1; #Plot_Name: Columns - Gravity
element elasticBeamColumn 9100021 1001 2001021 15.600000 2900.000000 111538.461538 1.940000 541.000000 57.700000 2; #Plot_Name:
element elasticBeamColumn 300021 2001022 2001028 15.600000 290.000000 11153.846154 1.940000 541.000000 57.700000 2; #Plot_Name: Beam - W24x131
element elasticBeamColumn 9200021 2001029 1021 15.600000 2900.000000 111538.461538 1.940000 541.000000 57.700000 2; #Plot_Name:
element elasticBeamColumn 9120221 1201 2201221 15.600000 2900.000000 111538.461538 1.940000 541.000000 57.700000 2; #Plot_Name:
element elasticBeamColumn 320221 2201222 2201228 15.600000 290.000000 11153.846154 1.940000 541.000000 57.700000 2; #Plot_Name: Beam - W24x131
element elasticBeamColumn 9220221 2201229 1221 15.600000 2900.000000 111538.461538 1.940000 541.000000 57.700000 2; #Plot_Name:
element elasticBeamColumn 9100201 1001 2001201 15.600000 2900.000000 111538.461538 1.940000 541.000000 57.700000 2; #Plot_Name:
element elasticBeamColumn 300201 2001202 2001208 15.600000 290.000000 11153.846154 1.940000 541.000000 57.700000 2; #Plot_Name: Beam - W24x131
element elasticBeamColumn 9200201 2001209 1201 15.600000 2900.000000 111538.461538 1.940000 541.000000 57.700000 2; #Plot_Name:
element elasticBeamColumn 9102221 1021 2021221 15.600000 2900.000000 111538.461538 1.940000 541.000000 57.700000 2; #Plot_Name:
element elasticBeamColumn 302221 2021222 2021228 15.600000 290.000000 11153.846154 1.940000 541.000000 57.700000 2; #Plot_Name: Beam - W24x131
element elasticBeamColumn 9202221 2021229 1221 15.600000 2900.000000 111538.461538 1.940000 541.000000 57.700000 2; #Plot_Name:
############################################################################
# Eigenvalue Analysis
############################################################################
puts "Importing TCL File Section: Eigenvalue Analysis"
set pi [expr 2.0*asin(1.0)];
set nModes 3;
set lambdaN [eigen 3];
set lambdaN1 [lindex $lambdaN [expr 0]];
set lambdaN2 [lindex $lambdaN [expr 1]];
set lambdaN3 [lindex $lambdaN [expr 2]];
set w1 [expr pow($lambdaN1,0.5)];
set w2 [expr pow($lambdaN2,0.5)];
set w3 [expr pow($lambdaN3,0.5)];
set T1 [expr 2.0*$pi/$w1];
set T2 [expr 2.0*$pi/$w2];
set T3 [expr 2.0*$pi/$w3];
puts "T1 = $T1 s"
puts "T2 = $T2 s"
puts "T3 = $T3 s"
3D Model - Rigid Diaphragm with Pinned Beams
Moderators: silvia, selimgunay, Moderators
Re: 3D Model - Rigid Diaphragm with Pinned Beams
your beams are spinning around their axis because of your constraints!
use these:
equalDOF 2021229 2021228 1 2 3 4
fix 1200 1 1 1 1 1 1;
equalDOF 2001201 2001202 1 2 3 4
fix 1220 1 1 1 1 1 1;
equalDOF 2201221 2201222 1 2 3 5
equalDOF 2001209 2001208 1 2 3 4
fix 1000 1 1 1 1 1 1;
fix 1020 1 1 1 1 1 1;
equalDOF 2001021 2001022 1 2 3 5
equalDOF 2021221 2021222 1 2 3 5
equalDOF 2201229 2201228 1 2 3 4
equalDOF 2001029 2001028 1 2 3 5
use these:
equalDOF 2021229 2021228 1 2 3 4
fix 1200 1 1 1 1 1 1;
equalDOF 2001201 2001202 1 2 3 4
fix 1220 1 1 1 1 1 1;
equalDOF 2201221 2201222 1 2 3 5
equalDOF 2001209 2001208 1 2 3 4
fix 1000 1 1 1 1 1 1;
fix 1020 1 1 1 1 1 1;
equalDOF 2001021 2001022 1 2 3 5
equalDOF 2021221 2021222 1 2 3 5
equalDOF 2201229 2201228 1 2 3 4
equalDOF 2001029 2001028 1 2 3 5