Thermal analysis in opensees

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

Moderators: silvia, selimgunay, Moderators

Post Reply
Deyanj
Posts: 1
Joined: Wed Aug 19, 2015 10:16 am

Thermal analysis in opensees

Post by Deyanj » Wed Aug 19, 2015 10:36 am

I am trying to analyze a 5 storey steel structure with some elements subjected to fire but it keeps failing at iteration 1781. Please what do i do? The code is posted below:

##########################################################################################
# 2D Frame with two bays subject to uniform temperature increment
# 6 elements per Beam
# 4 elements per Column
# Beam section is 356x171x51UB, 8 fibers for web;4 fibers for each flange
# Column section is 305x305x137UC, 8 fibers for web;4 fibers for each flange
# Plastic material class Steel01Thermal is used with fy=3.08e8
# Units are in N,m,s
# Geometric model #-------------------------------------------------------------

wipe;
file mkdir Data; # create data directory
model BasicBuilder -ndm 2 -ndf 3;
source DisplayPlane.tcl; # procedure for displaying a plane in model
source DisplayModel2D.tcl; # procedure for displaying 2D perspective of model

# nodal coordinates:
node 1 0 0;
node 2 6 0;
node 3 12 0;
node 4 18 0;
node 5 0 4;
node 6 6 4;
node 7 12 4;
node 8 18 4;
node 9 0 8;
node 10 6 8;
node 11 12 8;
node 12 18 8;
node 13 0 12;
node 14 6 12;
node 15 12 12;
node 16 18 12;
node 17 0 16;
node 18 6 16;
node 19 12 16;
node 20 18 16;
node 21 0 20;
node 22 6 20;
node 23 12 20;
node 24 18 20;
node 25 0 1;
node 26 0 2;
node 27 0 3;
node 28 0 5;
node 29 0 6;
node 30 0 7;
node 31 0 9;
node 32 0 10;
node 33 0 11;
node 34 0 13;
node 35 0 14;
node 36 0 15;
node 37 0 17;
node 38 0 18;
node 39 0 19;
node 40 6 1;
node 41 6 2;
node 42 6 3;
node 43 6 5;
node 44 6 6;
node 45 6 7;
node 46 6 9;
node 47 6 10;
node 48 6 11;
node 49 6 13;
node 50 6 14;
node 51 6 15;
node 52 6 17;
node 53 6 18;
node 54 6 19;
node 55 12 1;
node 56 12 2;
node 57 12 3;
node 58 12 5;
node 59 12 6;
node 60 12 7;
node 61 12 9;
node 62 12 10;
node 63 12 11;
node 64 12 13;
node 65 12 14;
node 66 12 15;
node 67 12 17;
node 68 12 18;
node 69 12 19;
node 70 18 1;
node 71 18 2;
node 72 18 3;
node 73 18 5;
node 74 18 6;
node 75 18 7;
node 76 18 9;
node 77 18 10;
node 78 18 11;
node 79 18 13;
node 80 18 14;
node 81 18 15;
node 82 18 17;
node 83 18 18;
node 84 18 19;

node 85 1 4;
node 86 2 4;
node 87 3 4;
node 88 4 4;
node 89 5 4;
node 90 7 4;
node 91 8 4;
node 92 9 4;
node 93 10 4;
node 94 11 4;
node 95 13 4;
node 96 14 4;
node 97 15 4;
node 98 16 4;
node 99 17 4;
node 100 1 8;
node 101 2 8;
node 102 3 8;
node 103 4 8;
node 104 5 8;
node 105 7 8;
node 106 8 8;
node 107 9 8;
node 108 10 8;
node 109 11 8;
node 110 13 8;
node 111 14 8;
node 112 15 8;
node 113 16 8;
node 114 17 8;
node 115 1 12;
node 116 2 12;
node 117 3 12;
node 118 4 12;
node 119 5 12;
node 120 7 12;
node 121 8 12;
node 122 9 12;
node 123 10 12;
node 124 11 12;
node 125 13 12;
node 126 14 12;
node 127 15 12;
node 128 16 12;
node 129 17 12;
node 130 1 16;
node 131 2 16;
node 132 3 16;
node 133 4 16;
node 134 5 16;
node 135 7 16;
node 136 8 16;
node 137 9 16;
node 138 10 16;
node 139 11 16;
node 140 13 16;
node 141 14 16;
node 142 15 16;
node 143 16 16;
node 144 17 16;
node 145 1 20;
node 146 2 20;
node 147 3 20;
node 148 4 20;
node 149 5 20;
node 150 7 20;
node 151 8 20;
node 152 9 20;
node 153 10 20;
node 154 11 20;
node 155 13 20;
node 156 14 20;
node 157 15 20;
node 158 16 20;
node 159 17 20;


# Single point constraints -- Boundary Conditions
fix 1 1 1 1;
fix 2 1 1 1;
fix 3 1 1 1;
fix 4 1 1 1;


#define steel01 material: $matTag $yieldStress $E $rat
uniaxialMaterial Steel01Thermal 1 3.08e8 2.1e11 0.01;

#define BEAM fibred section:fiber $yLoc $zLoc $A $matTag
# origin of section is the center of rectangular
section FiberThermal 1 {
# 8fibres in web
fiber 0.02075 0 0.0003071 1;
fiber 0.06225 0 0.0003071 1;
fiber 0.10375 0 0.0003071 1;
fiber 0.14525 0 0.0003071 1;
fiber -0.02075 0 0.0003071 1;
fiber -0.06225 0 0.0003071 1;
fiber -0.10375 0 0.0003071 1;
fiber -0.14525 0 0.0003071 1;
# 4 fibres in top flange
fiber 0.16743 0 0.000493 1;
fiber 0.17031 0 0.000493 1;
fiber 0.17319 0 0.000493 1;
fiber 0.17606 0 0.000493 1;
# 4 fibres in bottom flange
fiber -0.16743 0 0.000493 1;
fiber -0.17031 0 0.000493 1;
fiber -0.17319 0 0.000493 1;
fiber -0.17606 0 0.000493 1;

};

#define COLUMN fibred section:fiber $yLoc $zLoc $A $matTag
# origin of section is the center of rectangular
section FiberThermal 2 {
# 8fibres in web
fiber 0.01732 0 0.000478 1;
fiber 0.05195 0 0.000478 1;
fiber 0.08659 0 0.000478 1;
fiber 0.12123 0 0.000478 1;
fiber -0.01732 0 0.000478 1;
fiber -0.05195 0 0.000478 1;
fiber -0.08659 0 0.000478 1;
fiber -0.12123 0 0.000478 1;
# 4 fibres in top flange
fiber 0.14126 0 0.001677 1;
fiber 0.14669 0 0.001677 1;
fiber 0.15211 0 0.001677 1;
fiber 0.15754 0 0.001677 1;
# 4 fibres in bottom flange
fiber -0.14126 0 0.001677 1;
fiber -0.14669 0 0.001677 1;
fiber -0.15211 0 0.001677 1;
fiber -0.15754 0 0.001677 1;

};

#geomTransf choice is Linear, Corotational, PDelta ;
geomTransf Linear 1 ;

#define beam element: dispBeamColumnThermal $eleTag $iNode $jNode $numIntgrPts $secTag $TransfTag;
#"numIntgrPts" is the number of integration points along the element;
#"TransfTag" is pre-defined coordinate-transformation;
#element choice is dispBeamColumnTemperature or forceBeamColumnTemperature
element dispBeamColumnThermal 1 1 25 5 2 1;
element dispBeamColumnThermal 2 25 26 5 2 1;
element dispBeamColumnThermal 3 26 27 5 2 1;
element dispBeamColumnThermal 4 27 5 5 2 1;
element dispBeamColumnThermal 5 5 28 5 2 1;
element dispBeamColumnThermal 6 28 29 5 2 1;
element dispBeamColumnThermal 7 29 30 5 2 1;
element dispBeamColumnThermal 8 30 9 5 2 1;
element dispBeamColumnThermal 9 9 31 5 2 1;
element dispBeamColumnThermal 10 31 32 5 2 1;
element dispBeamColumnThermal 11 32 33 5 2 1;
element dispBeamColumnThermal 12 33 13 5 2 1;
element dispBeamColumnThermal 13 13 34 5 2 1;
element dispBeamColumnThermal 14 34 35 5 2 1;
element dispBeamColumnThermal 15 35 36 5 2 1;
element dispBeamColumnThermal 16 36 17 5 2 1;
element dispBeamColumnThermal 17 17 37 5 2 1;
element dispBeamColumnThermal 18 37 38 5 2 1;
element dispBeamColumnThermal 19 38 39 5 2 1;
element dispBeamColumnThermal 20 39 21 5 2 1;

element dispBeamColumnThermal 21 2 40 5 2 1;
element dispBeamColumnThermal 22 40 41 5 2 1;
element dispBeamColumnThermal 23 41 42 5 2 1;
element dispBeamColumnThermal 24 42 6 5 2 1;
element dispBeamColumnThermal 25 6 43 5 2 1;
element dispBeamColumnThermal 26 43 44 5 2 1;
element dispBeamColumnThermal 27 44 45 5 2 1;
element dispBeamColumnThermal 28 45 10 5 2 1;
element dispBeamColumnThermal 29 10 46 5 2 1;
element dispBeamColumnThermal 30 46 47 5 2 1;
element dispBeamColumnThermal 31 47 48 5 2 1;
element dispBeamColumnThermal 32 48 14 5 2 1;
element dispBeamColumnThermal 33 14 49 5 2 1;
element dispBeamColumnThermal 34 49 50 5 2 1;
element dispBeamColumnThermal 35 50 51 5 2 1;
element dispBeamColumnThermal 36 51 18 5 2 1;
element dispBeamColumnThermal 37 18 52 5 2 1;
element dispBeamColumnThermal 38 52 53 5 2 1;
element dispBeamColumnThermal 39 53 54 5 2 1;
element dispBeamColumnThermal 40 54 22 5 2 1;

element dispBeamColumnThermal 41 3 55 5 2 1;
element dispBeamColumnThermal 42 55 56 5 2 1;
element dispBeamColumnThermal 43 56 57 5 2 1;
element dispBeamColumnThermal 44 57 7 5 2 1;
element dispBeamColumnThermal 45 7 58 5 2 1;
element dispBeamColumnThermal 46 58 59 5 2 1;
element dispBeamColumnThermal 47 59 60 5 2 1;
element dispBeamColumnThermal 48 60 11 5 2 1;
element dispBeamColumnThermal 49 11 61 5 2 1;
element dispBeamColumnThermal 50 61 62 5 2 1;
element dispBeamColumnThermal 51 62 63 5 2 1;
element dispBeamColumnThermal 52 63 15 5 2 1;
element dispBeamColumnThermal 53 15 64 5 2 1;
element dispBeamColumnThermal 54 64 65 5 2 1;
element dispBeamColumnThermal 55 65 66 5 2 1;
element dispBeamColumnThermal 56 66 19 5 2 1;
element dispBeamColumnThermal 57 19 67 5 2 1;
element dispBeamColumnThermal 58 67 68 5 2 1;
element dispBeamColumnThermal 59 68 69 5 2 1;
element dispBeamColumnThermal 60 69 23 5 2 1;

element dispBeamColumnThermal 61 4 70 5 2 1;
element dispBeamColumnThermal 62 70 71 5 2 1;
element dispBeamColumnThermal 63 71 72 5 2 1;
element dispBeamColumnThermal 64 72 8 5 2 1;
element dispBeamColumnThermal 65 8 73 5 2 1;
element dispBeamColumnThermal 66 73 74 5 2 1;
element dispBeamColumnThermal 67 74 75 5 2 1;
element dispBeamColumnThermal 68 75 12 5 2 1;
element dispBeamColumnThermal 69 12 76 5 2 1;
element dispBeamColumnThermal 70 76 77 5 2 1;
element dispBeamColumnThermal 71 77 78 5 2 1;
element dispBeamColumnThermal 72 78 16 5 2 1;
element dispBeamColumnThermal 73 16 79 5 2 1;
element dispBeamColumnThermal 74 79 80 5 2 1;
element dispBeamColumnThermal 75 80 81 5 2 1;
element dispBeamColumnThermal 76 81 20 5 2 1;
element dispBeamColumnThermal 77 20 82 5 2 1;
element dispBeamColumnThermal 78 82 83 5 2 1;
element dispBeamColumnThermal 79 83 84 5 2 1;
element dispBeamColumnThermal 80 84 24 5 2 1;



element dispBeamColumnThermal 81 5 85 5 1 1;
element dispBeamColumnThermal 82 85 86 5 1 1;
element dispBeamColumnThermal 83 86 87 5 1 1;
element dispBeamColumnThermal 84 87 88 5 1 1;
element dispBeamColumnThermal 85 88 89 5 1 1;
element dispBeamColumnThermal 86 89 6 5 1 1;
element dispBeamColumnThermal 87 6 90 5 1 1;
element dispBeamColumnThermal 88 90 91 5 1 1;
element dispBeamColumnThermal 89 91 92 5 1 1;
element dispBeamColumnThermal 90 92 93 5 1 1;
element dispBeamColumnThermal 91 93 94 5 1 1;
element dispBeamColumnThermal 92 94 7 5 1 1;
element dispBeamColumnThermal 93 7 95 5 1 1;
element dispBeamColumnThermal 94 95 96 5 1 1;
element dispBeamColumnThermal 95 96 97 5 1 1;
element dispBeamColumnThermal 96 97 98 5 1 1;
element dispBeamColumnThermal 97 98 99 5 1 1;
element dispBeamColumnThermal 98 99 8 5 1 1;

element dispBeamColumnThermal 99 9 100 5 1 1;
element dispBeamColumnThermal 100 100 101 5 1 1;
element dispBeamColumnThermal 101 101 102 5 1 1;
element dispBeamColumnThermal 102 102 103 5 1 1;
element dispBeamColumnThermal 103 103 104 5 1 1;
element dispBeamColumnThermal 104 104 10 5 1 1;
element dispBeamColumnThermal 105 10 105 5 1 1;
element dispBeamColumnThermal 106 105 106 5 1 1;
element dispBeamColumnThermal 107 106 107 5 1 1;
element dispBeamColumnThermal 108 107 108 5 1 1;
element dispBeamColumnThermal 109 108 109 5 1 1;
element dispBeamColumnThermal 110 109 11 5 1 1;
element dispBeamColumnThermal 111 11 110 5 1 1;
element dispBeamColumnThermal 112 110 111 5 1 1;
element dispBeamColumnThermal 113 111 112 5 1 1;
element dispBeamColumnThermal 114 112 113 5 1 1;
element dispBeamColumnThermal 115 113 114 5 1 1;
element dispBeamColumnThermal 116 114 12 5 1 1;

element dispBeamColumnThermal 117 13 115 5 1 1;
element dispBeamColumnThermal 118 115 116 5 1 1;
element dispBeamColumnThermal 119 116 117 5 1 1;
element dispBeamColumnThermal 120 117 118 5 1 1;
element dispBeamColumnThermal 121 118 119 5 1 1;
element dispBeamColumnThermal 122 119 14 5 1 1;
element dispBeamColumnThermal 123 14 120 5 1 1;
element dispBeamColumnThermal 124 120 121 5 1 1;
element dispBeamColumnThermal 125 121 122 5 1 1;
element dispBeamColumnThermal 126 122 123 5 1 1;
element dispBeamColumnThermal 127 123 124 5 1 1;
element dispBeamColumnThermal 128 124 15 5 1 1;
element dispBeamColumnThermal 129 15 125 5 1 1;
element dispBeamColumnThermal 130 125 126 5 1 1;
element dispBeamColumnThermal 131 126 127 5 1 1;
element dispBeamColumnThermal 132 127 128 5 1 1;
element dispBeamColumnThermal 133 128 129 5 1 1;
element dispBeamColumnThermal 134 129 16 5 1 1;

element dispBeamColumnThermal 135 17 130 5 1 1;
element dispBeamColumnThermal 136 130 131 5 1 1;
element dispBeamColumnThermal 137 131 132 5 1 1;
element dispBeamColumnThermal 138 132 133 5 1 1;
element dispBeamColumnThermal 139 133 134 5 1 1;
element dispBeamColumnThermal 140 134 18 5 1 1;
element dispBeamColumnThermal 141 18 135 5 1 1;
element dispBeamColumnThermal 142 135 136 5 1 1;
element dispBeamColumnThermal 143 136 137 5 1 1;
element dispBeamColumnThermal 144 137 138 5 1 1;
element dispBeamColumnThermal 145 138 139 5 1 1;
element dispBeamColumnThermal 146 139 19 5 1 1;
element dispBeamColumnThermal 147 19 140 5 1 1;
element dispBeamColumnThermal 148 140 141 5 1 1;
element dispBeamColumnThermal 149 141 142 5 1 1;
element dispBeamColumnThermal 150 142 143 5 1 1;
element dispBeamColumnThermal 151 143 144 5 1 1;
element dispBeamColumnThermal 152 144 20 5 1 1;

element dispBeamColumnThermal 153 21 145 5 1 1;
element dispBeamColumnThermal 154 145 146 5 1 1;
element dispBeamColumnThermal 155 146 147 5 1 1;
element dispBeamColumnThermal 156 147 148 5 1 1;
element dispBeamColumnThermal 157 148 149 5 1 1;
element dispBeamColumnThermal 158 149 22 5 1 1;
element dispBeamColumnThermal 159 22 150 5 1 1;
element dispBeamColumnThermal 160 150 151 5 1 1;
element dispBeamColumnThermal 161 151 152 5 1 1;
element dispBeamColumnThermal 162 152 153 5 1 1;
element dispBeamColumnThermal 163 153 154 5 1 1;
element dispBeamColumnThermal 164 154 23 5 1 1;
element dispBeamColumnThermal 165 23 155 5 1 1;
element dispBeamColumnThermal 166 155 156 5 1 1;
element dispBeamColumnThermal 167 156 157 5 1 1;
element dispBeamColumnThermal 168 157 158 5 1 1;
element dispBeamColumnThermal 169 158 159 5 1 1;
element dispBeamColumnThermal 170 159 24 5 1 1;



#define output
recorder Node -file Data/YDispBeam.out -time -node 117 122 -dof 2 disp; # displacements of free nodes

recorder Node -file Data/XDispCol.out -time -node 32 47 62 -dof 1 disp; # displacements of free nodes

recorder Node -file Data/DispLeftJoint.out -time -node 13 -dof 1 2 3 disp; # displacements of free nodes

recorder Node -file Data/YTopJoints.out -time -node 21 22 -dof 2 disp; # displacements of free nodes

# view the deformed shape
set ViewScale 6;
DisplayModel2D DeformedShape $ViewScale ; # display deformed shape, the scaling factor needs to be adjusted for each model


#apply load

pattern Plain 1 Linear {

eleLoad -ele 81 -type -beamUniform -30600 0. 0.;
eleLoad -ele 82 -type -beamUniform -30600 0. 0.;
eleLoad -ele 83 -type -beamUniform -30600 0. 0.;
eleLoad -ele 84 -type -beamUniform -30600 0. 0.;
eleLoad -ele 85 -type -beamUniform -30600 0. 0.;
eleLoad -ele 86 -type -beamUniform -30600 0. 0.;
eleLoad -ele 87 -type -beamUniform -30600 0. 0.;
eleLoad -ele 88 -type -beamUniform -30600 0. 0.;
eleLoad -ele 89 -type -beamUniform -30600 0. 0.;
eleLoad -ele 90 -type -beamUniform -30600 0. 0.;
eleLoad -ele 91 -type -beamUniform -30600 0. 0.;
eleLoad -ele 92 -type -beamUniform -30600 0. 0.;
eleLoad -ele 93 -type -beamUniform -30600 0. 0.;
eleLoad -ele 94 -type -beamUniform -30600 0. 0.;
eleLoad -ele 95 -type -beamUniform -30600 0. 0.;
eleLoad -ele 96 -type -beamUniform -30600 0. 0.;
eleLoad -ele 97 -type -beamUniform -30600 0. 0.;
eleLoad -ele 98 -type -beamUniform -30600 0. 0.;

eleLoad -ele 99 -type -beamUniform -30600 0. 0.;
eleLoad -ele 100 -type -beamUniform -30600 0. 0.;
eleLoad -ele 101 -type -beamUniform -30600 0. 0.;
eleLoad -ele 102 -type -beamUniform -30600 0. 0.;
eleLoad -ele 103 -type -beamUniform -30600 0. 0.;
eleLoad -ele 104 -type -beamUniform -30600 0. 0.;
eleLoad -ele 105 -type -beamUniform -30600 0. 0.;
eleLoad -ele 106 -type -beamUniform -30600 0. 0.;
eleLoad -ele 107 -type -beamUniform -30600 0. 0.;
eleLoad -ele 108 -type -beamUniform -30600 0. 0.;
eleLoad -ele 109 -type -beamUniform -30600 0. 0.;
eleLoad -ele 110 -type -beamUniform -30600 0. 0.;
eleLoad -ele 111 -type -beamUniform -30600 0. 0.;
eleLoad -ele 112 -type -beamUniform -30600 0. 0.;
eleLoad -ele 113 -type -beamUniform -30600 0. 0.;
eleLoad -ele 114 -type -beamUniform -30600 0. 0.;
eleLoad -ele 115 -type -beamUniform -30600 0. 0.;
eleLoad -ele 116 -type -beamUniform -30600 0. 0.;

eleLoad -ele 117 -type -beamUniform -30600 0. 0.;
eleLoad -ele 118 -type -beamUniform -30600 0. 0.;
eleLoad -ele 119 -type -beamUniform -30600 0. 0.;
eleLoad -ele 120 -type -beamUniform -30600 0. 0.;
eleLoad -ele 121 -type -beamUniform -30600 0. 0.;
eleLoad -ele 122 -type -beamUniform -30600 0. 0.;
eleLoad -ele 123 -type -beamUniform -30600 0. 0.;
eleLoad -ele 124 -type -beamUniform -30600 0. 0.;
eleLoad -ele 125 -type -beamUniform -30600 0. 0.;
eleLoad -ele 126 -type -beamUniform -30600 0. 0.;
eleLoad -ele 127 -type -beamUniform -30600 0. 0.;
eleLoad -ele 128 -type -beamUniform -30600 0. 0.;
eleLoad -ele 129 -type -beamUniform -30600 0. 0.;
eleLoad -ele 130 -type -beamUniform -30600 0. 0.;
eleLoad -ele 131 -type -beamUniform -30600 0. 0.;
eleLoad -ele 132 -type -beamUniform -30600 0. 0.;
eleLoad -ele 133 -type -beamUniform -30600 0. 0.;
eleLoad -ele 134 -type -beamUniform -30600 0. 0.;

eleLoad -ele 135 -type -beamUniform -30600 0. 0.;
eleLoad -ele 136 -type -beamUniform -30600 0. 0.;
eleLoad -ele 137 -type -beamUniform -30600 0. 0.;
eleLoad -ele 138 -type -beamUniform -30600 0. 0.;
eleLoad -ele 139 -type -beamUniform -30600 0. 0.;
eleLoad -ele 140 -type -beamUniform -30600 0. 0.;
eleLoad -ele 141 -type -beamUniform -30600 0. 0.;
eleLoad -ele 142 -type -beamUniform -30600 0. 0.;
eleLoad -ele 143 -type -beamUniform -30600 0. 0.;
eleLoad -ele 144 -type -beamUniform -30600 0. 0.;
eleLoad -ele 145 -type -beamUniform -30600 0. 0.;
eleLoad -ele 146 -type -beamUniform -30600 0. 0.;
eleLoad -ele 147 -type -beamUniform -30600 0. 0.;
eleLoad -ele 148 -type -beamUniform -30600 0. 0.;
eleLoad -ele 149 -type -beamUniform -30600 0. 0.;
eleLoad -ele 150 -type -beamUniform -30600 0. 0.;
eleLoad -ele 151 -type -beamUniform -30600 0. 0.;
eleLoad -ele 152 -type -beamUniform -30600 0. 0.;

eleLoad -ele 153 -type -beamUniform -30600 0. 0.;
eleLoad -ele 154 -type -beamUniform -30600 0. 0.;
eleLoad -ele 155 -type -beamUniform -30600 0. 0.;
eleLoad -ele 156 -type -beamUniform -30600 0. 0.;
eleLoad -ele 157 -type -beamUniform -30600 0. 0.;
eleLoad -ele 158 -type -beamUniform -30600 0. 0.;
eleLoad -ele 159 -type -beamUniform -30600 0. 0.;
eleLoad -ele 160 -type -beamUniform -30600 0. 0.;
eleLoad -ele 161 -type -beamUniform -30600 0. 0.;
eleLoad -ele 162 -type -beamUniform -30600 0. 0.;
eleLoad -ele 163 -type -beamUniform -30600 0. 0.;
eleLoad -ele 164 -type -beamUniform -30600 0. 0.;
eleLoad -ele 165 -type -beamUniform -30600 0. 0.;
eleLoad -ele 166 -type -beamUniform -30600 0. 0.;
eleLoad -ele 167 -type -beamUniform -30600 0. 0.;
eleLoad -ele 168 -type -beamUniform -30600 0. 0.;
eleLoad -ele 169 -type -beamUniform -30600 0. 0.;
eleLoad -ele 170 -type -beamUniform -30600 0. 0.;
};

constraints Plain;
numberer Plain;
system BandGeneral;
test NormDispIncr 1e-8 100 ;
algorithm Newton;
integrator LoadControl 0.1
analysis Static;
analyze 10;
loadConst -time 0.0

#define thermal load (i.e. temperature distribution in section)
#-beamThermal $T1 $LocY1 $T2 LocY2....; two temperature means uniform or linear temperature distribution
# T1 is the bottom temp of beam and T2 is the top
# Y1 is the coordinate of bottom of beam section; Y2 is for top
# the temperature will be interpolated along the section

#PK create timeseries for fire
timeSeries Path 2 -time {0.0 50 100 150 200 250 300 350 400 450 500 550 600 650 700 750 800 850 900 950 1000 1050 1100 1150 1200 1250 1300 1350 1400 1800 2200 2600 3000 3300 4500 5000} -values {0.0215 0.3488 0.4493 0.5107 0.5550 0.5897 0.6182 0.6425 0.6635 0.6821 0.6988 0.7139 0.7277 0.7404 0.7522 0.7631 0.7734 0.7831 0.7922 0.8008 0.8090 0.8167 0.8242 0.8313 0.8381 0.8446 0.8509 0.8569 0.8627 0.9029 0.9350 0.9618 0.9847 1.0 0.0215 0.0215}

timeSeries Path 3 -time {600 650 700 750 800 850 900 950 1000 1050 1100 1150 1200 1250 1300 1350 1400 1800 2200 2600 3000 3300 4500 5000} -values {0.0222 0.3604 0.4643 0.5277 0.5735 0.6093 0.6388 0.6638 0.6855 0.7048 0.7220 0.7376 0.7519 0.7650 0.7772 0.7885 0.7991 0.8659 0.9134 0.9503 0.9805 1.0 0.0222 0.0222}

#PK create fire pattern
pattern Plain 2 2 {

eleLoad -ele 29 -type -beamThermal 932.3313 -0.16025 932.3313 0.16025;
eleLoad -ele 30 -type -beamThermal 932.3313 -0.16025 932.3313 0.16025;
eleLoad -ele 31 -type -beamThermal 932.3313 -0.16025 932.3313 0.16025;
eleLoad -ele 32 -type -beamThermal 932.3313 -0.16025 932.3313 0.16025;
eleLoad -ele 49 -type -beamThermal 932.3313 -0.16025 932.3313 0.16025;
eleLoad -ele 50 -type -beamThermal 932.3313 -0.16025 932.3313 0.16025;
eleLoad -ele 51 -type -beamThermal 932.3313 -0.16025 932.3313 0.16025;
eleLoad -ele 52 -type -beamThermal 932.3313 -0.16025 932.3313 0.16025;


eleLoad -ele 123 -type -beamThermal 932.3313 -0.1775 932.3313 0.1775;
eleLoad -ele 124 -type -beamThermal 932.3313 -0.1775 932.3313 0.1775;
eleLoad -ele 125 -type -beamThermal 932.3313 -0.1775 932.3313 0.1775;
eleLoad -ele 126 -type -beamThermal 932.3313 -0.1775 932.3313 0.1775;
eleLoad -ele 127 -type -beamThermal 932.3313 -0.1775 932.3313 0.1775;
eleLoad -ele 128 -type -beamThermal 932.3313 -0.1775 932.3313 0.1775;
};

pattern Plain 3 3 {
eleLoad -ele 9 -type -beamThermal 902.34 -0.16025 902.34 0.16025;
eleLoad -ele 10 -type -beamThermal 902.34 -0.16025 902.34 0.16025;
eleLoad -ele 11 -type -beamThermal 902.34 -0.16025 902.34 0.16025;
eleLoad -ele 12 -type -beamThermal 902.34 -0.16025 902.34 0.16025;

eleLoad -ele 117 -type -beamThermal 902.34 -0.1775 902.34 0.1775;
eleLoad -ele 118 -type -beamThermal 902.34 -0.1775 902.34 0.1775;
eleLoad -ele 119 -type -beamThermal 902.34 -0.1775 902.34 0.1775;
eleLoad -ele 120 -type -beamThermal 902.34 -0.1775 902.34 0.1775;
eleLoad -ele 121 -type -beamThermal 902.34 -0.1775 902.34 0.1775;
eleLoad -ele 122 -type -beamThermal 902.34 -0.1775 902.34 0.1775;
};


set Nstep 5000;
set Factor [expr 5000./$Nstep]; # first load increment;
integrator LoadControl $Factor; # determine the next time step for an analysis
analysis Static; # define type of analysis static or transient
analyze $Nstep ; # apply fire load
Top
Post Reply

Return to “OpenSees.exe Users”