~~~~~~~~~~~~~~~~~~~~~~
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 20
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 1
~~~~~~~~~~~~~~~~
# define units------------------------------------------------
set mm 1.;
set Kg 1.;
set sec 1.;
set LunitTXT "mm";
set FunitTXT "Kg";
set TunitTXT "sec";
set m [expr 1000.*$mm];
set N [expr $Kg*$m/pow($sec,2)];
set KN [expr 1000.*$N];
set KNm [expr $KN*$m];
set Pa [expr $N/pow($m,2)];
set MPa [expr 1.e6*$Pa];
set GPa [expr 1.e9*$Pa]
set ru [expr $Kg/pow($m,3)]; # density
set mm2 [expr $mm*$mm];
set mm4 [expr $mm*$mm*$mm*$mm];
set PI [expr 2*asin(1.0)];
source LibUnits.tcl
wipe
model basic -ndm 2 -ndf 3
node 1 0 0
node 2 500 0
node 3 1000 0
node 4 1500 0
node 5 2000 0
node 6 2500 0
node 7 3000 0
fix 1 0 1 0
fix 7 1 1 0
set fpc [expr -34.2*$MPa];
set E [expr 36.3*$GPa];
set epsc0 [expr $fpc/$E];
set fpcu [expr $fpc*0.05];
set epsU -0.02
set lambda 0.1
set ft [expr 3*$MPa]
set Ets [expr $ft/0.002];
uniaxialMaterial Concrete01 4 $fpc $epsc0 $fpcu $epsU
set Fy [expr 335*$MPa]
set E [expr 200*$GPa]
#set E [expr 180*$GPa]
set bate 0.00390
set R0 18.5
set CR1 0.925
set CR2 0.15
set Fu [expr 530*$MPa]
#set Fu [expr 510*$MPa]
set Esh [expr $E*0.043]
set esh 0.005
set eult 0.12
uniaxialMaterial Steel02 1 $Fy $E $bate $R0 $CR1 $CR2
set Fu [expr 1863.3*$MPa]
set Fy [expr 0.8*$Fu*$MPa]
set E [expr 195*$GPa]
set Fp [expr 1197*$MPa]
uniaxialMaterial Steel02 5 $Fy $E 0.02 18 0.925 0.15 0 1 0 1 $Fp;
set S 200
set d1 40
set d2 40
section Fiber 1 {
patch rect 4 3 10 -150 100 150 150
patch rect 4 3 6 -100 -150 100 -100
fiber -50 -125 $S 1
fiber 50 -125 $S 1
}
geomTransf Linear 1
set eletype forceBeamColumn;
set num 3
element $eletype 1 1 2 $num 1 1
element $eletype 2 2 3 $num 1 1
element $eletype 3 3 4 $num 1 1
element $eletype 4 4 5 $num 1 1
element $eletype 5 5 6 $num 1 1
element $eletype 6 6 7 $num 1 1
node 101 0 -50
node 102 500 -50
node 103 1000 -50
node 104 1500 -50
node 105 2000 -50
node 106 2500 -50
node 107 3000 -50
section Fiber 4 {
patch circ 5 20 1 0 0 0 5.92 0 360
fiber 50 -125 $S 1
}
element truss 11 101 102 4
element truss 12 102 103 4
element truss 13 103 104 4
element truss 14 104 105 4
element truss 15 105 106 4
element truss 16 106 107 4
rigidLink beam 1 101
rigidLink beam 4 104
rigidLink beam 7 107
recorder Node -file node4.out -time -node 4 -dof 2 disp
set W [expr -0.63*$KN/$m]
pattern Plain 1 Linear {
for {set level 1} {$level <=6} {incr level 1} {
set eleID $level
eleLoad -ele $eleID -type -beamUniform $W
}
}
constraints Transformation
numberer Plain
system BandGeneral
test EnergyIncr 1.0e-1 200
algorithm Newton
integrator LoadControl 1;
analysis Static
analyze 1
loadConst -time 0.0
set F [expr 0.5*$KN]
pattern Plain 2 Linear {
load 3 0 [expr -$F] 0
load 5 0 [expr -$F] 0
}
constraints Transformation
numberer Plain
system BandGeneral
test EnergyIncr 1.0e-6 200
algorithm Newton
analysis Static
integrator DisplacementControl 4 2 -1
analyze 60
~~~~~~~~~~~~~~~~~~
help ,plz. I change most parameter ,it make no difference
.