BeamThermal(3D) not working!!!

[reza378salami]

Hi all,

I'm working on a 3d model with temperature variation. The 2D model works perfectly fine but in the 3D cant apply the temperature load. I was wondering it is working on the latest version, or maybe not compiled yet.

The error I get is:
"SectionForceDeformation::getTemperatureStress<double*dataMix>-sould not be called
Thermal section not defined, defult 0 return"


here is a simple code for 1 3D column.
Thank you for your help
Cheers
Reza

#-----------------------
wipe;
model BasicBuilder -ndm 3 -ndf 6;

# node#, X, Y Z
node 1 0 0 0;
node 2 0 1 0;
node 3 0 2 0;
node 4 0 3 0;
node 5 0 4 0;

fix 1 1 1 1 1 1 1;
fix 5 1 0 1 1 1 1;

set matID 1
set secID 2;
set ColTransfTag 3;
set SecTagTorsion 4;
set ColSecTag3D 5;

uniaxialMaterial Steel01Thermal $matID 3.55e8 2e11 0.0;

# ELEMENT properties
# Structural-Steel W-section properties
set d 0.08; # depth
set bf 0.046; # flange width
set tf 0.0052; # flange thickness
set tw 0.0038; # web thickness
set nfdw 16; # number of fibers along dw
set nftw 2; # number of fibers along tw
set nfbf 16; # number of fibers along bf
set nftf 4; # number of fibers along tf

set dw [expr $d - 2 * $tf]
set y1 [expr -$d/2]
set y2 [expr -$dw/2]
set y3 [expr $dw/2]
set y4 [expr $d/2]
set z1 [expr -$bf/2]
set z2 [expr -$tw/2]
set z3 [expr $tw/2]
set z4 [expr $bf/2]

section FiberThermal $secID {
# nfIJ nfJK yI zI yJ zJ yK zK yL zL
patch quadr $matID $nfbf $nftf $y1 $z4 $y1 $z1 $y2 $z1 $y2 $z4
patch quadr $matID $nftw $nfdw $y2 $z3 $y2 $z2 $y3 $z2 $y3 $z3
patch quadr $matID $nfbf $nftf $y3 $z4 $y3 $z1 $y4 $z1 $y4 $z4
}

# associate a tag to column transformation
geomTransf Linear $ColTransfTag 0 0 1;
set Ubig 9.E10;
uniaxialMaterial Elastic $SecTagTorsion $Ubig
section Aggregator $ColSecTag3D $SecTagTorsion T -section $secID

set numIntgrPts 3;
element dispBeamColumnThermal 1 1 2 $numIntgrPts $ColSecTag3D $ColTransfTag;
element dispBeamColumnThermal 2 2 3 $numIntgrPts $ColSecTag3D $ColTransfTag;
element dispBeamColumnThermal 3 3 4 $numIntgrPts $ColSecTag3D $ColTransfTag;
element dispBeamColumnThermal 4 4 5 $numIntgrPts $ColSecTag3D $ColTransfTag;

recorder Node -file TDFree5_3D.out -time -node 5 -dof 2 disp;

pattern Plain 1 Linear {
load 5 0 -100000 0 0 0 0;
}

# Gravity-analysis parameters -- load-controlled static analysis
set Tol 1.0e-8;
constraints Plain;
numberer Plain;
system BandGeneral;
test NormDispIncr $Tol 100 ;
algorithm Newton;
set NstepGravity 100;
set DGravity [expr 1./$NstepGravity];
integrator LoadControl $DGravity;
analysis Static;
analyze $NstepGravity;
loadConst -time 0.0

# ----------------------------------------------------
# Beam3dThermalAction
#eleLoad -ele $eleTag -type -beamThermal $T1 $Y1 $T2 $Y2
pattern Plain 2 Linear {
eleLoad -ele 1 -type -beamThermal 500 -0.04 500 0.04
eleLoad -ele 2 -type -beamThermal 500 -0.04 500 0.04
eleLoad -ele 3 -type -beamThermal 500 -0.04 500 0.04
eleLoad -ele 4 -type -beamThermal 500 -0.04 500 0.04
};
integrator LoadControl 0.01
analysis Static;
analyze 100;
loadConst -time 0.0