factorization failed, matrix singular

[cexuhan]

Hi all,

I used VS2010 to create a new dynamic link library(DLL) file for the uniaxial material class and successfully invoked it in opensees. Now I want to draw a simple load-displacement curve of the column structure with this new uniaxial material, and the solver reported an error during the analysis. The errors are as follows:

WARNING BandGenLinLapackSolver::solve() -factorization failed, matrix singular U(i,i) = 0, i= 1
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 68 with domain at load factor -1.#IND
OpenSees > analyze failed, returned: -3 error flag

When I used an existing uniaxialmaterial instead, the model operated successfully. Can anyone tell me the reason why this kind of mistake happened? And which part of my developed model can cause this sort of problem?

Here are whole TCL commands in opensees （SteelConfinedConcrete is my created material）

wipe;

model basic -ndm 2 -ndf 3
set e0 40;
set gap 9.145;
node 1 [expr (0*$gap)] 0
node 2 [expr (1*$gap)] 0
node 3 [expr (2*$gap)] 0
node 4 [expr (3*$gap)] 0
node 5 [expr (4*$gap)] 0
node 6 [expr (5*$gap)] 0
node 7 [expr (6*$gap)] 0
node 8 [expr (7*$gap)] 0
node 9 [expr (8*$gap)] 0
node 10 [expr (9*$gap)] 0
node 11 [expr (10*$gap)] 0
node 12 [expr (11*$gap)] 0
node 13 [expr (12*$gap)] 0
node 14 [expr (13*$gap)] 0
node 15 [expr (14*$gap)] 0
node 16 [expr (15*$gap)] 0
node 17 [expr (16*$gap)] 0
node 18 [expr (17*$gap)] 0
node 19 [expr (18*$gap)] 0
node 20 [expr (19*$gap)] 0
node 21 [expr (20*$gap)] 0
node 22 [expr (21*$gap)] 0
node 23 [expr (22*$gap)] 0
node 24 [expr (23*$gap)] 0
node 25 [expr (24*$gap)] 0
node 26 [expr (25*$gap)] 0
node 27 [expr (26*$gap)] 0
node 28 [expr (27*$gap)] 0
node 29 [expr (28*$gap)] 0
node 30 [expr (29*$gap)] 0
node 31 [expr (30*$gap)] 0
node 32 [expr (31*$gap)] 0
node 33 [expr (32*$gap)] 0
node 34 [expr (33*$gap)] 0
node 35 [expr (34*$gap)] 0
node 36 [expr (35*$gap)] 0
node 37 [expr (36*$gap)] 0
node 38 [expr (37*$gap)] 0
node 39 [expr (38*$gap)] 0
node 40 [expr (39*$gap)] 0
node 41 [expr (40*$gap)] 0
node 42 [expr (41*$gap)] 0
node 43 [expr (42*$gap)] 0
node 44 [expr (43*$gap)] 0
node 45 [expr (44*$gap)] 0
node 46 [expr (45*$gap)] 0
node 47 [expr (46*$gap)] 0
node 48 [expr (47*$gap)] 0
node 49 [expr (48*$gap)] 0
node 50 [expr (49*$gap)] 0
node 51 [expr (50*$gap)] 0
node 52 [expr (51*$gap)] 0
node 53 [expr (52*$gap)] 0
node 54 [expr (53*$gap)] 0
node 55 [expr (54*$gap)] 0
node 56 [expr (55*$gap)] 0
node 57 [expr (56*$gap)] 0
node 58 [expr (57*$gap)] 0
node 59 [expr (58*$gap)] 0
node 60 [expr (59*$gap)] 0
node 61 [expr (60*$gap)] 0
node 62 [expr (61*$gap)] 0
node 63 [expr (62*$gap)] 0
node 64 [expr (63*$gap)] 0
node 65 [expr (64*$gap)] 0
node 66 [expr (65*$gap)] 0
node 67 [expr (66*$gap)] 0
node 68 [expr (67*$gap)] 0
node 69 [expr (68*$gap)] 0
node 70 [expr (69*$gap)] 0
node 71 [expr (70*$gap)] 0
node 72 [expr (71*$gap)] 0
node 73 [expr (72*$gap)] 0
node 74 [expr (73*$gap)] 0
node 75 [expr (74*$gap)] 0
node 76 [expr (75*$gap)] 0
node 77 [expr (76*$gap)] 0
node 78 [expr (77*$gap)] 0
node 79 [expr (78*$gap)] 0
node 80 [expr (79*$gap)] 0
node 81 [expr (80*$gap)] 0
node 82 [expr (81*$gap)] 0
node 83 [expr (82*$gap)] 0
node 84 [expr (83*$gap)] 0
node 85 [expr (84*$gap)] 0
node 86 [expr (85*$gap)] 0
node 87 [expr (86*$gap)] 0
node 88 [expr (87*$gap)] 0
node 89 [expr (88*$gap)] 0
node 90 [expr (89*$gap)] 0
node 91 [expr (90*$gap)] 0
node 92 [expr (91*$gap)] 0
node 93 [expr (92*$gap)] 0
node 94 [expr (93*$gap)] 0
node 95 [expr (94*$gap)] 0
node 96 [expr (95*$gap)] 0
node 97 [expr (96*$gap)] 0
node 98 [expr (97*$gap)] 0
node 99 [expr (98*$gap)] 0
node 100 [expr (99*$gap)] 0
node 101 [expr (100*$gap)] 0
node 102 [expr (101*$gap)] 0
node 103 [expr (102*$gap)] 0
node 104 [expr (103*$gap)] 0
node 105 [expr (104*$gap)] 0
node 106 [expr (105*$gap)] 0
node 107 [expr (106*$gap)] 0
node 108 [expr (107*$gap)] 0
node 109 [expr (108*$gap)] 0
node 110 [expr (109*$gap)] 0
node 111 [expr (110*$gap)] 0
node 112 [expr (111*$gap)] 0
node 113 [expr (112*$gap)] 0
node 114 [expr (113*$gap)] 0
node 115 [expr (114*$gap)] 0
node 116 [expr (115*$gap)] 0
node 117 [expr (116*$gap)] 0
node 118 [expr (117*$gap)] 0
node 119 [expr (118*$gap)] 0
node 120 [expr (119*$gap)] 0
node 121 [expr (120*$gap)] 0
node 122 [expr (121*$gap)] 0
node 123 [expr (122*$gap)] 0
node 124 [expr (123*$gap)] 0
node 125 [expr (124*$gap)] 0
node 126 [expr (125*$gap)] 0
node 127 [expr (126*$gap)] 0
node 128 [expr (127*$gap)] 0
node 129 [expr (128*$gap)] 0
node 130 [expr (129*$gap)] 0
node 131 [expr (130*$gap)] 0
node 132 [expr (131*$gap)] 0
node 133 [expr (132*$gap)] 0
node 134 [expr (133*$gap)] 0
node 135 [expr (134*$gap)] 0
node 136 [expr (135*$gap)] 0
node 137 [expr (136*$gap)] 0
node 138 [expr (137*$gap)] 0
node 139 [expr (138*$gap)] 0
node 140 [expr (139*$gap)] 0
node 141 [expr (140*$gap)] 0
node 142 [expr (141*$gap)] 0
node 143 [expr (142*$gap)] 0
node 144 [expr (143*$gap)] 0
node 145 [expr (144*$gap)] 0
node 146 [expr (145*$gap)] 0
node 147 [expr (146*$gap)] 0
node 148 [expr (147*$gap)] 0
node 149 [expr (148*$gap)] 0
node 150 [expr (149*$gap)] 0
node 151 [expr (150*$gap)] 0
node 152 [expr (151*$gap)] 0
node 153 [expr (152*$gap)] 0
node 154 [expr (153*$gap)] 0
node 155 [expr (154*$gap)] 0
node 156 [expr (155*$gap)] 0
node 157 [expr (156*$gap)] 0
node 158 [expr (157*$gap)] 0
node 159 [expr (158*$gap)] 0
node 160 [expr (159*$gap)] 0
node 161 [expr (160*$gap)] 0
node 162 [expr (161*$gap)] 0
node 163 [expr (162*$gap)] 0
node 164 [expr (163*$gap)] 0
node 165 [expr (164*$gap)] 0
node 166 [expr (165*$gap)] 0
node 167 [expr (166*$gap)] 0
node 168 [expr (167*$gap)] 0
node 169 [expr (168*$gap)] 0
node 170 [expr (169*$gap)] 0
node 171 [expr (170*$gap)] 0
node 172 [expr (171*$gap)] 0
node 173 [expr (172*$gap)] 0
node 174 [expr (173*$gap)] 0
node 175 [expr (174*$gap)] 0
node 176 [expr (175*$gap)] 0
node 177 [expr (176*$gap)] 0
node 178 [expr (177*$gap)] 0
node 179 [expr (178*$gap)] 0
node 180 [expr (179*$gap)] 0
node 181 [expr (180*$gap)] 0
node 182 [expr (181*$gap)] 0
node 183 [expr (182*$gap)] 0
node 184 [expr (183*$gap)] 0
node 185 [expr (184*$gap)] 0
node 186 [expr (185*$gap)] 0
node 187 [expr (186*$gap)] 0
node 188 [expr (187*$gap)] 0
node 189 [expr (188*$gap)] 0
node 190 [expr (189*$gap)] 0
node 191 [expr (190*$gap)] 0
node 192 [expr (191*$gap)] 0
node 193 [expr (192*$gap)] 0
node 194 [expr (193*$gap)] 0
node 195 [expr (194*$gap)] 0
node 196 [expr (195*$gap)] 0
node 197 [expr (196*$gap)] 0
node 198 [expr (197*$gap)] 0
node 199 [expr (198*$gap)] 0
node 200 [expr (199*$gap)] 0
node 201 [expr (200*$gap)] 0

node 202 0 [expr $e0]
node 203 [expr (200*$gap)] [expr $e0]
fix 1 0 1 0
fix 201 0 1 0
set CSec 1;
set Esec 2;

#define MATERIAL TAG------------------------------------------------------------------

#This is assigning a number to the different material models used
set CSteel 1
set CConc 2

#define GEOMETRY----------------------------------------------------------------------

set CODSteel 152.4; #Outernal diameter of steel
set CThSteel 3.12; #steel thickness

set CODConc [expr $CODSteel - 2*$CThSteel]; #Internal diameter of concrete



#define NUMBER OF SUBDIVISION FOR EACH MATERIAL--------------------------------------

#This is defining the discretization for the patch

set CnumSubCircSteel 16;
set CnumSubRadSteel 2;

set CnumSubCircConc 16
set CnumSubRadConc 20

set PI 3.1415926535
set Cec0 -0.002;
set Cfc0 -55; #Unconfined concrete strength
set Cfy -382.0; # Yield strength of the steel material
set Cvst 0.25; #Poisson ratios of the steel tube
set CR [expr 0.5*$CODSteel]; #-Outer radius of circular column section
set CEs 198000; #Elastic modulus of steel
set Cm 4.0;
set Cfcc -87.5893; #-115(e=0), -90(e=100000) #Compressive strength of FRP confined concrete (Max)
set Cecc -0.0107; #Empirical coefficient to evaluate the strength of confined concrete
set CEc [expr 4730.0*sqrt(-0.85*$Cfc0)]; #Compressive elastic modulus of unconfined concrete
# set CEc 48000; #Compressive elastic modulus of unconfined concrete
# set Cft [expr 0.632*sqrt(-$Cfc0)]; #Tensile strength
set Cft 4.70;
set CEts [expr 0.05*$CEc];
#This is the Steel tube confined concrete model. This is not available in OpenSees. To use this model you have to have the SteelConfinedConcrete.dll in the same folder as this file.
# uniaxialMaterial FRPConfinedConcrete02 $CConc $Cfc0 $CEc $Cec0 -Ultimate $Cfcc $Cecc $Cft $CEts [expr 1]
# uniaxialMaterial SteelConfinedConcrete $CConc $Cfc0 $CEc $Cec0 -SteelC $CThSteel $Cfy $Cvst $CR $CEs $Cm $Cft [expr 1]
# uniaxialMaterial Concrete04 $CConc $Cfc0 $Cecc -0.025 $CEc $Cft 0.002
uniaxialMaterial SteelConfinedConcrete $CConc $Cfc0 $CEc $Cec0 -Ultimate $Cfcc $Cecc $Cft [expr 1]
uniaxialMaterial Steel02 $CSteel [expr -$Cfy] $CEs 0.01 18 0.925 0.15
section Fiber 1 {
patch circ $CSteel $CnumSubCircSteel $CnumSubRadSteel 0 0 [expr 0.5*$CODConc] [expr 0.5*$CODSteel] 0 360
patch circ $CConc $CnumSubCircConc $CnumSubRadConc 0 0 0 [expr 0.5*$CODConc] 0 360
}
set nP 8;
geomTransf Corotational 1
element dispBeamColumn 1 1 2 $nP $CSec 1
element dispBeamColumn 2 2 3 $nP $CSec 1
element dispBeamColumn 3 3 4 $nP $CSec 1
element dispBeamColumn 4 4 5 $nP $CSec 1
element dispBeamColumn 5 5 6 $nP $CSec 1
element dispBeamColumn 6 6 7 $nP $CSec 1
element dispBeamColumn 7 7 8 $nP $CSec 1
element dispBeamColumn 8 8 9 $nP $CSec 1
element dispBeamColumn 9 9 10 $nP $CSec 1
element dispBeamColumn 10 10 11 $nP $CSec 1
element dispBeamColumn 11 11 12 $nP $CSec 1
element dispBeamColumn 12 12 13 $nP $CSec 1
element dispBeamColumn 13 13 14 $nP $CSec 1
element dispBeamColumn 14 14 15 $nP $CSec 1
element dispBeamColumn 15 15 16 $nP $CSec 1
element dispBeamColumn 16 16 17 $nP $CSec 1
element dispBeamColumn 17 17 18 $nP $CSec 1
element dispBeamColumn 18 18 19 $nP $CSec 1
element dispBeamColumn 19 19 20 $nP $CSec 1
element dispBeamColumn 20 20 21 $nP $CSec 1
element dispBeamColumn 21 21 22 $nP $CSec 1
element dispBeamColumn 22 22 23 $nP $CSec 1
element dispBeamColumn 23 23 24 $nP $CSec 1
element dispBeamColumn 24 24 25 $nP $CSec 1
element dispBeamColumn 25 25 26 $nP $CSec 1
element dispBeamColumn 26 26 27 $nP $CSec 1
element dispBeamColumn 27 27 28 $nP $CSec 1
element dispBeamColumn 28 28 29 $nP $CSec 1
element dispBeamColumn 29 29 30 $nP $CSec 1
element dispBeamColumn 30 30 31 $nP $CSec 1
element dispBeamColumn 31 31 32 $nP $CSec 1
element dispBeamColumn 32 32 33 $nP $CSec 1
element dispBeamColumn 33 33 34 $nP $CSec 1
element dispBeamColumn 34 34 35 $nP $CSec 1
element dispBeamColumn 35 35 36 $nP $CSec 1
element dispBeamColumn 36 36 37 $nP $CSec 1
element dispBeamColumn 37 37 38 $nP $CSec 1
element dispBeamColumn 38 38 39 $nP $CSec 1
element dispBeamColumn 39 39 40 $nP $CSec 1
element dispBeamColumn 40 40 41 $nP $CSec 1
element dispBeamColumn 41 41 42 $nP $CSec 1
element dispBeamColumn 42 42 43 $nP $CSec 1
element dispBeamColumn 43 43 44 $nP $CSec 1
element dispBeamColumn 44 44 45 $nP $CSec 1
element dispBeamColumn 45 45 46 $nP $CSec 1
element dispBeamColumn 46 46 47 $nP $CSec 1
element dispBeamColumn 47 47 48 $nP $CSec 1
element dispBeamColumn 48 48 49 $nP $CSec 1
element dispBeamColumn 49 49 50 $nP $CSec 1
element dispBeamColumn 50 50 51 $nP $CSec 1
element dispBeamColumn 51 51 52 $nP $CSec 1
element dispBeamColumn 52 52 53 $nP $CSec 1
element dispBeamColumn 53 53 54 $nP $CSec 1
element dispBeamColumn 54 54 55 $nP $CSec 1
element dispBeamColumn 55 55 56 $nP $CSec 1
element dispBeamColumn 56 56 57 $nP $CSec 1
element dispBeamColumn 57 57 58 $nP $CSec 1
element dispBeamColumn 58 58 59 $nP $CSec 1
element dispBeamColumn 59 59 60 $nP $CSec 1
element dispBeamColumn 60 60 61 $nP $CSec 1
element dispBeamColumn 61 61 62 $nP $CSec 1
element dispBeamColumn 62 62 63 $nP $CSec 1
element dispBeamColumn 63 63 64 $nP $CSec 1
element dispBeamColumn 64 64 65 $nP $CSec 1
element dispBeamColumn 65 65 66 $nP $CSec 1
element dispBeamColumn 66 66 67 $nP $CSec 1
element dispBeamColumn 67 67 68 $nP $CSec 1
element dispBeamColumn 68 68 69 $nP $CSec 1
element dispBeamColumn 69 69 70 $nP $CSec 1
element dispBeamColumn 70 70 71 $nP $CSec 1
element dispBeamColumn 71 71 72 $nP $CSec 1
element dispBeamColumn 72 72 73 $nP $CSec 1
element dispBeamColumn 73 73 74 $nP $CSec 1
element dispBeamColumn 74 74 75 $nP $CSec 1
element dispBeamColumn 75 75 76 $nP $CSec 1
element dispBeamColumn 76 76 77 $nP $CSec 1
element dispBeamColumn 77 77 78 $nP $CSec 1
element dispBeamColumn 78 78 79 $nP $CSec 1
element dispBeamColumn 79 79 80 $nP $CSec 1
element dispBeamColumn 80 80 81 $nP $CSec 1
element dispBeamColumn 81 81 82 $nP $CSec 1
element dispBeamColumn 82 82 83 $nP $CSec 1
element dispBeamColumn 83 83 84 $nP $CSec 1
element dispBeamColumn 84 84 85 $nP $CSec 1
element dispBeamColumn 85 85 86 $nP $CSec 1
element dispBeamColumn 86 86 87 $nP $CSec 1
element dispBeamColumn 87 87 88 $nP $CSec 1
element dispBeamColumn 88 88 89 $nP $CSec 1
element dispBeamColumn 89 89 90 $nP $CSec 1
element dispBeamColumn 90 90 91 $nP $CSec 1
element dispBeamColumn 91 91 92 $nP $CSec 1
element dispBeamColumn 92 92 93 $nP $CSec 1
element dispBeamColumn 93 93 94 $nP $CSec 1
element dispBeamColumn 94 94 95 $nP $CSec 1
element dispBeamColumn 95 95 96 $nP $CSec 1
element dispBeamColumn 96 96 97 $nP $CSec 1
element dispBeamColumn 97 97 98 $nP $CSec 1
element dispBeamColumn 98 98 99 $nP $CSec 1
element dispBeamColumn 99 99 100 $nP $CSec 1
element dispBeamColumn 100 100 101 $nP $CSec 1
element dispBeamColumn 101 101 102 $nP $CSec 1
element dispBeamColumn 102 102 103 $nP $CSec 1
element dispBeamColumn 103 103 104 $nP $CSec 1
element dispBeamColumn 104 104 105 $nP $CSec 1
element dispBeamColumn 105 105 106 $nP $CSec 1
element dispBeamColumn 106 106 107 $nP $CSec 1
element dispBeamColumn 107 107 108 $nP $CSec 1
element dispBeamColumn 108 108 109 $nP $CSec 1
element dispBeamColumn 109 109 110 $nP $CSec 1
element dispBeamColumn 110 110 111 $nP $CSec 1
element dispBeamColumn 111 111 112 $nP $CSec 1
element dispBeamColumn 112 112 113 $nP $CSec 1
element dispBeamColumn 113 113 114 $nP $CSec 1
element dispBeamColumn 114 114 115 $nP $CSec 1
element dispBeamColumn 115 115 116 $nP $CSec 1
element dispBeamColumn 116 116 117 $nP $CSec 1
element dispBeamColumn 117 117 118 $nP $CSec 1
element dispBeamColumn 118 118 119 $nP $CSec 1
element dispBeamColumn 119 119 120 $nP $CSec 1
element dispBeamColumn 120 120 121 $nP $CSec 1
element dispBeamColumn 121 121 122 $nP $CSec 1
element dispBeamColumn 122 122 123 $nP $CSec 1
element dispBeamColumn 123 123 124 $nP $CSec 1
element dispBeamColumn 124 124 125 $nP $CSec 1
element dispBeamColumn 125 125 126 $nP $CSec 1
element dispBeamColumn 126 126 127 $nP $CSec 1
element dispBeamColumn 127 127 128 $nP $CSec 1
element dispBeamColumn 128 128 129 $nP $CSec 1
element dispBeamColumn 129 129 130 $nP $CSec 1
element dispBeamColumn 130 130 131 $nP $CSec 1
element dispBeamColumn 131 131 132 $nP $CSec 1
element dispBeamColumn 132 132 133 $nP $CSec 1
element dispBeamColumn 133 133 134 $nP $CSec 1
element dispBeamColumn 134 134 135 $nP $CSec 1
element dispBeamColumn 135 135 136 $nP $CSec 1
element dispBeamColumn 136 136 137 $nP $CSec 1
element dispBeamColumn 137 137 138 $nP $CSec 1
element dispBeamColumn 138 138 139 $nP $CSec 1
element dispBeamColumn 139 139 140 $nP $CSec 1
element dispBeamColumn 140 140 141 $nP $CSec 1
element dispBeamColumn 141 141 142 $nP $CSec 1
element dispBeamColumn 142 142 143 $nP $CSec 1
element dispBeamColumn 143 143 144 $nP $CSec 1
element dispBeamColumn 144 144 145 $nP $CSec 1
element dispBeamColumn 145 145 146 $nP $CSec 1
element dispBeamColumn 146 146 147 $nP $CSec 1
element dispBeamColumn 147 147 148 $nP $CSec 1
element dispBeamColumn 148 148 149 $nP $CSec 1
element dispBeamColumn 149 149 150 $nP $CSec 1
element dispBeamColumn 150 150 151 $nP $CSec 1
element dispBeamColumn 151 151 152 $nP $CSec 1
element dispBeamColumn 152 152 153 $nP $CSec 1
element dispBeamColumn 153 153 154 $nP $CSec 1
element dispBeamColumn 154 154 155 $nP $CSec 1
element dispBeamColumn 155 155 156 $nP $CSec 1
element dispBeamColumn 156 156 157 $nP $CSec 1
element dispBeamColumn 157 157 158 $nP $CSec 1
element dispBeamColumn 158 158 159 $nP $CSec 1
element dispBeamColumn 159 159 160 $nP $CSec 1
element dispBeamColumn 160 160 161 $nP $CSec 1
element dispBeamColumn 161 161 162 $nP $CSec 1
element dispBeamColumn 162 162 163 $nP $CSec 1
element dispBeamColumn 163 163 164 $nP $CSec 1
element dispBeamColumn 164 164 165 $nP $CSec 1
element dispBeamColumn 165 165 166 $nP $CSec 1
element dispBeamColumn 166 166 167 $nP $CSec 1
element dispBeamColumn 167 167 168 $nP $CSec 1
element dispBeamColumn 168 168 169 $nP $CSec 1
element dispBeamColumn 169 169 170 $nP $CSec 1
element dispBeamColumn 170 170 171 $nP $CSec 1
element dispBeamColumn 171 171 172 $nP $CSec 1
element dispBeamColumn 172 172 173 $nP $CSec 1
element dispBeamColumn 173 173 174 $nP $CSec 1
element dispBeamColumn 174 174 175 $nP $CSec 1
element dispBeamColumn 175 175 176 $nP $CSec 1
element dispBeamColumn 176 176 177 $nP $CSec 1
element dispBeamColumn 177 177 178 $nP $CSec 1
element dispBeamColumn 178 178 179 $nP $CSec 1
element dispBeamColumn 179 179 180 $nP $CSec 1
element dispBeamColumn 180 180 181 $nP $CSec 1
element dispBeamColumn 181 181 182 $nP $CSec 1
element dispBeamColumn 182 182 183 $nP $CSec 1
element dispBeamColumn 183 183 184 $nP $CSec 1
element dispBeamColumn 184 184 185 $nP $CSec 1
element dispBeamColumn 185 185 186 $nP $CSec 1
element dispBeamColumn 186 186 187 $nP $CSec 1
element dispBeamColumn 187 187 188 $nP $CSec 1
element dispBeamColumn 188 188 189 $nP $CSec 1
element dispBeamColumn 189 189 190 $nP $CSec 1
element dispBeamColumn 190 190 191 $nP $CSec 1
element dispBeamColumn 191 191 192 $nP $CSec 1
element dispBeamColumn 192 192 193 $nP $CSec 1
element dispBeamColumn 193 193 194 $nP $CSec 1
element dispBeamColumn 194 194 195 $nP $CSec 1
element dispBeamColumn 195 195 196 $nP $CSec 1
element dispBeamColumn 196 196 197 $nP $CSec 1
element dispBeamColumn 197 197 198 $nP $CSec 1
element dispBeamColumn 198 198 199 $nP $CSec 1
element dispBeamColumn 199 199 200 $nP $CSec 1
element dispBeamColumn 200 200 201 $nP $CSec 1
# element dispBeamColumn 201 1 202 $nP $CSec 1
# element dispBeamColumn 202 201 203 $nP $CSec 1
element elasticBeamColumn 201 1 202 10000 1000000 1000000000 1
element elasticBeamColumn 202 201 203 10000 1000000 1000000000 1

recorder Node -file CFSTdisplacement.out -time -node 101 -dof 2 disp
recorder Element -file Concrete_StressStrainT.out -time -ele 100 101 section 1 fiber [expr ($CODSteel/2-$CThSteel)] 0 $CConc stressStrain
recorder Element -file Concrete_StressStrainM.out -time -ele 100 101 section 1 fiber 0 0 $CConc stressStrain
recorder Element -file Concrete_StressStrainB.out -time -ele 100 101 section 1 fiber [expr -($CODSteel/2-$CThSteel)] 0 $CConc stressStrain
recorder Element -file SteelTube_StressStrainT.out -time -ele 100 101 section 1 fiber [expr $CODSteel/2] 0 $CSteel stressStrain
# pattern Plain 101 Constant {
# load 1 100000 0 0
# load 201 -100000 0 0

# }
# integrator LoadControl 1
# system SparseGeneral -piv;
# test EnergyIncr 1e-10 200
# numberer Plain
# constraints Plain
# algorithm Newton
# analysis Static
# analyze 5
pattern Plain 102 Linear {
load 1 1 0 -$e0
load 201 -1 0 $e0
# load 202 1 0 0
# load 203 -1 0 0
# load 101 0 -1 0;
}
puts "Computing"
constraints Plain
numberer Plain
system BandGeneral
test EnergyIncr 1.0e-5 200
algorithm Newton
integrator DisplacementControl 101 2 -0.3
# integrator LoadControl 10000
analysis Static
analyze 200


Any tiny help or suggestion is needed. Thanks a lot.

[selimgunay]

It looks like the stiffness matrix cannot be constructed properly with the new material

[cexuhan]

Hi Semiguna

It does like that. I found that there are some problems occurred in my reloading and unloading part so that the overall stiffness matrix cannot be decomposed by PLU. Now I had solved it already. Thanks anyway.

Best
Xu

[Feliciasatriya]

Dear cecuhan,

How did you solve it? What changes did you make to your model ?

[sinceyourethere]

Hey OP,

Managed to get past that error? If yes, please let me know, as I've did all I could to make it work but I always get a similar error with yours.

Thanks for your help.

[aminulislam]

Hi cexuhan, How you solve it? It would be great if you kindly share.