How to add temperature load

[J.Zhang]

Hi,
I tried to add temperature load into a beam, but got the information:
"WARNING: CTestEnergyIncr::test() - failed to converge after: 6 iterations
NewtnRaphson::solveCurrentStep() -the ConvergenceTest object failed in test()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at load factor 0.1"

I can not find why. The code is following:



wipe
set dataDir Data
file mkdir $dataDir;
model BasicBuilder -ndm 2 -ndf 3

# Define Nodal Coordinates

node 1 0.0 0.0 0.0
node 2 3.0 0.0 0.0

puts "Nodes defined"

# Define Nodal Restraints

fix 1 1 1 0
fix 2 0 1 0
puts "Nodal Restr defined"

# Define Elastic Material Properties

uniaxialMaterial Elastic 1 29000.0
puts "ElasticMat defined"

# Define Geometric Transformations

geomTransf Linear 1
puts "Linear GeomTransf defined"

# Define Elastic Beam Column Elements

element elasticBeamColumn 1 1 2 50.0 29000.0 0.0050 1


# Define Node Recorder

recorder Node -file $dataDir/Node_2_disp.out -time -node 2 -dof 1 2 3 disp
puts "Node recorder defined"

# Define Element Recorder

recorder Element -file $dataDir/Element_1_globalForce.out -time -ele 1 globalForce
puts "Element recorder defined"



# Define Gravity

pattern Plain 100 Linear {
eleLoad -ele 1 -type -beamUniform 100.0
}
constraints Plain;

numberer Plain;

system BandGeneral;

test EnergyIncr 1.0e-8 6 ;

algorithm Newton;

integrator LoadControl 0.1;

analysis Static

analyze 10;

loadConst -time 0.0;
puts "gravity defined"

# `Define TempLoad
pattern Plain 101 Linear {
eleLoad -ele 1 -type -beamTemp 10. 20.
}
integrator LoadControl 0.1;
analyze 10;

puts "done!"