search

OpenSeesforFire / openseesforfire.github.io

This a repository for OpenSees for fire
http://openseesforfire.github.io/
18 stars 14 forks source link

Beam3dThermalAction problem #4

Open reza378salami opened 6 years ago

reza378salami commented 6 years ago

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

reza378salami commented 6 years ago

I also tried this but not working:

----------------------------------------------------

pattern Plain 2 Linear {

eleLoad -ele $eleTag -type -beamThermal $T1 $Y1 ... $T5 $Y5 $T6 $T7 $Z1 ...$T14 $T15 $Z5

eleLoad -ele 2 -type -beamThermal 300 -0.0348 300 -0.0174 300 0.0 300 0.0174 300 0.0348 300 300 0.023 300 300 0.015 300 300 0.0 300 300 -0.015 300 300 -0.023

};

li-ming-jiang commented 6 years ago

Hi,

I will resolve your problem soon and let you know.

Regards, Liming

Sent from my iPhone

On 17 Nov 2017, at 3:18 AM, reza378salami notifications@github.com<mailto:notifications@github.com> wrote:

I also tried this but not working:

pattern Plain 2 Linear {

eleLoad -ele $eleTag -type -beamThermal $T1 $Y1 ... $T5 $Y5 $T6 $T7 $Z1 ...$T14 $T15 $Z5

eleLoad -ele 2 -type -beamThermal 300 -0.0348 300 -0.0174 300 0.0 300 0.0174 300 0.0348 300 300 0.023 300 300 0.015 300 300 0.0 300 300 -0.015 300 300 -0.023

};

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHubhttps://github.com/OpenSeesforFire/openseesforfire.github.io/issues/4#issuecomment-345033174, or mute the threadhttps://github.com/notifications/unsubscribe-auth/AJE4vcbiZC1psaN5-z9K-jnYLhkokvMzks5s3IqHgaJpZM4Qg7bo.

li-ming-jiang commented 6 years ago

Hi，

I have updated the executable file of OpenSees for fire on openseesforfire.github.io website.

Would please try it for your case?

Thanks!

Best regards, Liming

Sent from Mailhttps://go.microsoft.com/fwlink/?LinkId=550986 for Windows 10

From: reza378salami notifications@github.com Sent: Friday, November 17, 2017 3:18:31 AM To: OpenSeesforFire/openseesforfire.github.io Cc: Subscribed Subject: Re: [OpenSeesforFire/openseesforfire.github.io] Beam3dThermalAction problem (#4)

I also tried this but not working:

pattern Plain 2 Linear {

eleLoad -ele $eleTag -type -beamThermal $T1 $Y1 ... $T5 $Y5 $T6 $T7 $Z1 ...$T14 $T15 $Z5

eleLoad -ele 2 -type -beamThermal 300 -0.0348 300 -0.0174 300 0.0 300 0.0174 300 0.0348 300 300 0.023 300 300 0.015 300 300 0.0 300 300 -0.015 300 300 -0.023

};

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHubhttps://github.com/OpenSeesforFire/openseesforfire.github.io/issues/4#issuecomment-345033174, or mute the threadhttps://github.com/notifications/unsubscribe-auth/AJE4vcbiZC1psaN5-z9K-jnYLhkokvMzks5s3IqHgaJpZM4Qg7bo.

hexiaohuangcan commented 5 years ago

Hello Liming,

I face the same problem with Reza by adopting the .exe file downloaded in this website. It says "SectionForceDeformation::getTemperatureStress(double *dataMixed) - should not be calle".

Thank you for your help. Regards Shaw

hexiaohuangcan commented 5 years ago

Hello Liming,

I found the trick. Replace "beamThermal" with "beam3dThermal", then the program works well.

Regards Shaw

hexiaohuangcan commented 5 years ago

Hello,

I'm modelling a 3D two-beam substructure under fire condition. However, after the gravitational load was added, the displacement hardly changed under thermal action. By contrast, for 2D structure it worked quite well.

I wonder how thermal action should be correctly added? I've tried the following three cases, but all failed:

# eleLoad -range 1 6 -type -beam3dThermal -source Tdata.dat -$HalfD  $HalfD -$HalfB $HalfB; 
# eleLoad -range 11 16 -type -beam3dThermal -source Tdata.dat -$HalfD  $HalfD -$HalfB $HalfB; 

# eleLoad -range 1 6 -type -beam3dThermal 1000 -$HalfD  1000 $HalfD 
# eleLoad -range 11 16 -type -beam3dThermal  1000 -$HalfD 1000 $HalfD 

eleLoad -range 1 6 -type -beam3dThermalAction 1000 -$HalfD  1000 $HalfD 
eleLoad -range 11 16 -type -beam3dThermalAction  1000 -$HalfD 1000 $HalfD

Could you help me?

Attached is the source file.

Regards Shaw

example-SteelBeam-uniform_3d.txt