fire resistance analysis of concrete slab modelled with solid226 and reinf264 elements

I have some problems with elements solid226 (for thermal + structural analysis) and reinf264 to make them work correctly together in a different temperatures then a normal one (room - approx.20°C).
 
It seems like the material model "MISO" defined in dependance on temperature for the elements REINF264 applied into base elements solid226 is not working the correct way. The tamperature of the reinforcement elements is shown correctly, something around the temperature of the base ones, but then it achieves stress which should not be achieved in that temperature.
If the elements solid185 are used as the base ones, the stress - strain dependance of the reinf elements seems to be reasonable there. (although so far I have not achieved succesfully some higher temperatures the structure should withstand, but at least it behaves not the same way as if there is normal room temperature)
 
any experience with using these elements together?

Comments

  • edited July 2018

    Hello Dadojin,

      When you say it doesn't work correctly, can you please elaborate and show some pictures? It would in general help if you could post snapshots about to model/settings?

    Regards,

    Sandeep

  • DadojinDadojin Member
    edited July 2018

    !

    ! ------- GEOMETRIA -------

    !

    xxxx=2.000! 1/2 rozpatia

    yyyy=0.075! 1/2 vysky

    zzzz=0.075! 1/2 sirky

    rozmerelementu=0.025

    !

    /PREP7

    !

    K,1,0,0,0

    K,2,0,-1*yyyy,0

    K,3,0,yyyy,0

    K,7,0,-1*yyyy,zzzz

    K,8,0,0,zzzz

    K,9,0,yyyy,zzzz

    !

    A,8,1,3,9

    A,7,2,1,8

    !

    KSEL,S,LOC,X,0

    LSLK,S,1

    ASLL,S,1

    VEXT,ALL,,,xxxx-yyyy,

    !

    KSEL,S,LOC,X,xxxx-yyyy

    LSLK,S,1

    ASLL,S,1

    VEXT,ALL,,,yyyy,

    !

    KSEL,S,LOC,X,xxxx

    LSLK,S,1

    ASLL,S,1

    VEXT,ALL,,,yyyy

    !

    ! velkosti

    !

    ALLS

    KSEL,S,LOC,X,0

    LSLK,S,1

    KSEL,S,LOC,X,xxxx-yyyy

    LSLK,A,1

    KSEL,S,LOC,X,xxxx

    LSLK,A,1

    KSEL,S,LOC,X,xxxx+yyyy

    LSLK,A,1

    CM,ciary_z_y,lines

    LESIZE,ALL,rozmerelementu, , , ,1, , ,1,

    !

    ALLS

    KSEL,S,LOC,X,0

    KSEL,A,LOC,X,xxxx-yyyy

    LSLK,S,1

    CMSEL,U,ciary_z_y

    CM,ciary_nosnik,lines

    LESIZE,ALL,2*rozmerelementu, , ,3,1, , ,1,

    !

    ALLS

    KSEL,S,LOC,X,xxxx

    LSLK,S,0

    CMSEL,U,ciary_z_y

    CMSEL,U,ciary_nosnik

    CM,ciary_podpora,lines

    LESIZE,ALL,yyyy, , , ,1, , ,1,

    !

    ALLS

    /REPLO

    !

    ! ------- DEFINOVANIE BETONU -------

    !

    ! - ELASTIC CONSTANTS - C25/30

    beton=25! OZNACENIE ELEMENTU <-----

    joungov_bet =31E9 ! [Pa]

    poisson_bet =0.2! [-]

    epsilon_c1 =0.0021! [-]

    epsilon_cu1 =0.0035! [-]

    hustota_bet =2300! [kg/m^3]

    !

    ET,beton,SOLID226! element

    KEYOPT,beton,1,11! structural + thermal

    ! KEYOPT,beton,2,1 ! 0 = strong / 1 = weak coupling

    !

    ! Joungov modul pruznosti a poissonov koeficient

    !

    MPTEMP,,,,,,,,

    MPTEMP,1,0

    MPTEMP,2,20

    MPTEMP,3,57

    MPTEMP,4,100

    MPTEMP,5,200

    MPTEMP,6,300

    MPTEMP,7,400

    MPTEMP,8,500

    MPTEMP,9,600

    MPTEMP,10,700

    MPTEMP,11,800

    MPTEMP,12,900

    MPTEMP,13,1000

    MPTEMP,14,1100

    MPTEMP,15,1200

    !

    MPDATA,EX,beton,,31E9

    MPDATA,EX,beton,,31E9

    MPDATA,EX,beton,,31E9

    MPDATA,EX,beton,,29.13421E9

    MPDATA,EX,beton,,24.79518E9

    MPDATA,EX,beton,,20.45614E9

    MPDATA,EX,beton,,16.11711E9

    MPDATA,EX,beton,,11.77807E9

    MPDATA,EX,beton,,3.9E9

    MPDATA,EX,beton,,3.1E9

    MPDATA,EX,beton,,2.48E9

    MPDATA,EX,beton,,1.86E9

    MPDATA,EX,beton,,1.24E9

    MPDATA,EX,beton,,0.62E9

    MPDATA,EX,beton,,1

    !

    MPDATA,NUXY,beton,,poisson_beton

    MPDATA,NUXY,beton,,poisson_beton

    MPDATA,NUXY,beton,,poisson_beton

    MPDATA,NUXY,beton,,poisson_beton

    MPDATA,NUXY,beton,,poisson_beton

    MPDATA,NUXY,beton,,poisson_beton

    MPDATA,NUXY,beton,,poisson_beton

    MPDATA,NUXY,beton,,poisson_beton

    MPDATA,NUXY,beton,,poisson_beton

    MPDATA,NUXY,beton,,poisson_beton

    MPDATA,NUXY,beton,,poisson_beton

    MPDATA,NUXY,beton,,poisson_beton

    MPDATA,NUXY,beton,,poisson_beton

    MPDATA,NUXY,beton,,poisson_beton

    MPDATA,NUXY,beton,,poisson_beton

    !

    ! dlzkova teplotna roztaznost

    !

    ! Kremicite kamenivo

    !

    MPTEMP,,,,,,,,

    MPTEMP,1,0

    MPTEMP,2,20

    MPTEMP,3,100

    MPTEMP,4,200

    MPTEMP,5,300

    MPTEMP,6,400

    MPTEMP,7,500

    MPTEMP,8,600

    MPTEMP,9,700

    MPTEMP,10,1200

    !

    MPDATA,THSX,beton,,0

    MPDATA,THSX,beton,,0.000000184

    MPDATA,THSX,beton,,0.000743

    MPDATA,THSX,beton,,0.001804

    MPDATA,THSX,beton,,0.003141

    MPDATA,THSX,beton,,0.004892

    MPDATA,THSX,beton,,0.007195

    MPDATA,THSX,beton,,0.010188

    MPDATA,THSX,beton,,0.014

    MPDATA,THSX,beton,,0.0141

    !

    !  Menetrey-Willam beton, exponencialne zmakcenie, kamenivo kremicite

    Rc =33E6! [Pa] - uniaxial compressive strength

    Rt =2.6E6! [Pa] - uniaxial tensile strength

    psi =11.5! [°] - dilatancy angle

    !

    TB,CONCRETE,beton,15,,MW

    TBTEMP,0

    TBDATA,1,33E6,2.6E6,39.6E6

    TBTEMP,20

    TBDATA,1,33E6,2.6E6,39.6E6

    TBTEMP,57

    TBDATA,1,33E6,2.6E6,39.6E6

    TBTEMP,100

    TBDATA,1,33E6,2.6E6,39.6E6

    TBTEMP,200

    TBDATA,1,31.35E6,2.08E6,37.62E6

    TBTEMP,300

    TBDATA,1,28.05E6,1.56E6,33.66E6

    TBTEMP,400

    TBDATA,1,24.75E6,1.04E6,29.7E6

    TBTEMP,500

    TBDATA,1,19.8E6,0.52E6,23.76E6

    TBTEMP,600

    TBDATA,1,14.85E6,3.5  ,17.82E6

    TBTEMP,700

    TBDATA,1,9.9E6,3  ,11.88E6

    TBTEMP,800

    TBDATA,1,4.95E6,2.5  ,5.94E6

    TBTEMP,900

    TBDATA,1,2.64E6,2  ,3.168E6

    TBTEMP,1000

    TBDATA,1,1.32E6,1.5  ,1.584E6

    TBTEMP,1100

    TBDATA,1,0.33E6,1  ,0.396E6

    TBTEMP,1200

    TBDATA,1,1  ,0.5  ,1.2 

    !

    TB,CONCRETE,beton,2,,DILA

    TBTEMP,0

    TBDATA,1,psi

    TBTEMP,1200

    TBDATA,1,psi

    !

    TB,CONCRETE,beton,15,,HSD2

    TBTEMP,0

    TBDATA,1,0.00143548,0.01935831,0.4,0.6028,0.2,65

    TBDATA,7,0.1

    TBTEMP,20

    TBDATA,1,0.00143548,0.01935831,0.4,0.6028,0.2,65

    TBDATA,7,0.1

    TBTEMP,57

    TBDATA,1,0.00212923,0.02051456,0.4,0.6028,0.2,65

    TBDATA,7,0.1

    TBTEMP,100

    TBDATA,1,0.00286731,0.02181722,0.4,0.6028,0.2,65

    TBDATA,7,0.1

    TBTEMP,200

    TBDATA,1,0.00423564,0.02423785,0.4,0.6028,0.2,65

    TBDATA,7,0.1

    TBTEMP,300

    TBDATA,1,0.00562877,0.02667343,0.4,0.6028,0.2,65

    TBDATA,7,0.1

    TBTEMP,400

    TBDATA,1,0.00846436,0.02907432,0.4,0.6028,0.2,65

    TBDATA,7,0.1

    TBTEMP,500

    TBDATA,1,0.01331891,0.03148664,0.4,0.6028,0.2,65

    TBDATA,7,0.1

    TBTEMP,600

    TBDATA,1,0.02119231,0.03270473,0.4,0.6028,0.2,65

    TBDATA,7,0.1

    TBTEMP,700

    TBDATA,1,0.02180645,0.03557493,0.4,0.6028,0.2,65

    TBDATA,7,0.1

    TBTEMP,800

    TBDATA,1,0.02300403,0.03879683,0.4,0.6028,0.2,65

    TBDATA,7,0.1

    TBTEMP,900

    TBDATA,1,0.02358065,0.04164442,0.4,0.6028,0.2,65

    TBDATA,7,0.1

    TBTEMP,1000

    TBDATA,1,0.02393548,0.04435831,0.4,0.6028,0.2,65

    TBDATA,7,0.1

    TBTEMP,1100

    TBDATA,1,0.02446774,0.04717916,0.4,0.6028,0.2,65

    TBDATA,7,0.1

    TBTEMP,1200

    TBDATA,1,0.025,0.048,0.4,0.6028,0.2,65

    TBDATA,7,0.1

    !

    ! Definicia hustoty betonu [kg / m^3 ]

    !

    MPTEMP,,,,,,,,

    MPTEMP,1,0

    MPTEMP,2,20

    MPTEMP,3,115

    MPTEMP,4,200

    MPTEMP,5,400

    MPTEMP,6,1200

    !

    MPDATA,DENS,beton,,hustota_bet

    MPDATA,DENS,beton,,hustota_bet

    MPDATA,DENS,beton,,hustota_bet

    MPDATA,DENS,beton,,0.98*hustota_bet

    MPDATA,DENS,beton,,0.95*hustota_bet

    MPDATA,DENS,beton,,0.88*hustota_bet

    !

    ! Definicia mernej tepelnej kapacity betonu c [J/ (kg*K) ]

    ! pre vlhkost 1,5% hmotnosti betonu

    !

    MPTEMP,,,,,,,,

    MPTEMP,1,0

    MPTEMP,2,20

    MPTEMP,3,100

    MPTEMP,4,101

    MPTEMP,5,115

    MPTEMP,6,200

    MPTEMP,7,400

    MPTEMP,8,1200

    !

    MPDATA,C,beton,,900

    MPDATA,C,beton,,900

    MPDATA,C,beton,,900

    MPDATA,C,beton,,1470

    MPDATA,C,beton,,1470

    MPDATA,C,beton,,1000

    MPDATA,C,beton,,1100

    MPDATA,C,beton,,1100

    !

    ! Definicia tepelnej vodivosti lambda [W/ (m*K) ]

    ! Dolna medz

    !

    MPTEMP,,,,,,,,

    MPTEMP,1,0

    MPTEMP,2,20

    MPTEMP,3,100

    MPTEMP,4,200

    MPTEMP,5,300

    MPTEMP,6,400

    MPTEMP,7,500

    MPTEMP,8,600

    MPTEMP,9,700

    MPTEMP,10,800

    MPTEMP,11,900

    MPTEMP,12,1000

    MPTEMP,13,1100

    MPTEMP,14,1200

    !

    MPDATA,KXX,beton,,1.36

    MPDATA,KXX,beton,,1.333028

    MPDATA,KXX,beton,,1.2297

    MPDATA,KXX,beton,,1.1108

    MPDATA,KXX,beton,,1.0033

    MPDATA,KXX,beton,,0.9072

    MPDATA,KXX,beton,,0.8225

    MPDATA,KXX,beton,,0.7492

    MPDATA,KXX,beton,,0.6873

    MPDATA,KXX,beton,,0.6368

    MPDATA,KXX,beton,,0.5977

    MPDATA,KXX,beton,,0.57

    MPDATA,KXX,beton,,0.5537

    MPDATA,KXX,beton,,0.5488

    !

    ! ------- DEFINOVANIE PODLOZKY -------

    !

    podlozka = 500! OZNACENIE ELEMENTU <-----

    !

    ET,podlozka,SOLID226! element

    KEYOPT,podlozka,1,11! structural + thermal

    ! KEYOPT,podlozka,2,1 ! 0 = strong / 1 = weak coupling

    !

    MPTEMP,,,,,,,,

    MPTEMP,1,0

    MPTEMP,2,20

    MPTEMP,3,57

    MPTEMP,4,100

    MPTEMP,5,200

    MPTEMP,6,300

    MPTEMP,7,400

    MPTEMP,8,500

    MPTEMP,9,600

    MPTEMP,10,700

    MPTEMP,11,800

    MPTEMP,12,900

    MPTEMP,13,1000

    MPTEMP,14,1100

    MPTEMP,15,1200

    !

    MPDATA,EX,podlozka,,31E9

    MPDATA,EX,podlozka,,31E9

    MPDATA,EX,podlozka,,31E9

    MPDATA,EX,podlozka,,29.13421E9

    MPDATA,EX,podlozka,,24.79518E9

    MPDATA,EX,podlozka,,20.45614E9

    MPDATA,EX,podlozka,,16.11711E9

    MPDATA,EX,podlozka,,3.9E9

    MPDATA,EX,podlozka,,3.1E9

    MPDATA,EX,podlozka,,2.48E9

    MPDATA,EX,podlozka,,1.86E9

    MPDATA,EX,podlozka,,1.24E9

    MPDATA,EX,podlozka,,0.62E9

    MPDATA,EX,podlozka,,1

    !

    MPDATA,NUXY,podlozka,,poisson_beton

    MPDATA,NUXY,podlozka,,poisson_beton

    MPDATA,NUXY,podlozka,,poisson_beton

    MPDATA,NUXY,podlozka,,poisson_beton

    MPDATA,NUXY,podlozka,,poisson_beton

    MPDATA,NUXY,podlozka,,poisson_beton

    MPDATA,NUXY,podlozka,,poisson_beton

    MPDATA,NUXY,podlozka,,poisson_beton

    MPDATA,NUXY,podlozka,,poisson_beton

    MPDATA,NUXY,podlozka,,poisson_beton

    MPDATA,NUXY,podlozka,,poisson_beton

    MPDATA,NUXY,podlozka,,poisson_beton

    MPDATA,NUXY,podlozka,,poisson_beton

    MPDATA,NUXY,podlozka,,poisson_beton

    MPDATA,NUXY,podlozka,,poisson_beton

    !

    MPTEMP,,,,,,,,

    MPTEMP,1,0

    MPTEMP,2,20

    MPTEMP,3,100

    MPTEMP,4,200

    MPTEMP,5,300

    MPTEMP,6,400

    MPTEMP,7,500

    MPTEMP,8,600

    MPTEMP,9,700

    MPTEMP,10,1200

    !

    MPDATA,THSX,podlozka,,0

    MPDATA,THSX,podlozka,,0.000000184

    MPDATA,THSX,podlozka,,0.000743

    MPDATA,THSX,podlozka,,0.001804

    MPDATA,THSX,podlozka,,0.003141

    MPDATA,THSX,podlozka,,0.004892

    MPDATA,THSX,podlozka,,0.007195

    MPDATA,THSX,podlozka,,0.010188

    MPDATA,THSX,podlozka,,0.014

    MPDATA,THSX,podlozka,,0.0154

    !

    MPTEMP,,,,,,,,

    MPTEMP,1,0

    MPTEMP,2,20

    MPTEMP,3,115

    MPTEMP,4,200

    MPTEMP,5,400

    MPTEMP,6,1200

    !

    MPDATA,DENS,podlozka,,hustota_bet

    MPDATA,DENS,podlozka,,hustota_bet

    MPDATA,DENS,podlozka,,hustota_bet

    MPDATA,DENS,podlozka,,0.98*hustota_bet

    MPDATA,DENS,podlozka,,0.95*hustota_bet

    MPDATA,DENS,podlozka,,0.88*hustota_bet

    !

    MPTEMP,,,,,,,,

    MPTEMP,1,0

    MPTEMP,2,20

    MPTEMP,3,100

    MPTEMP,4,101

    MPTEMP,5,115

    MPTEMP,6,200

    MPTEMP,7,400

    MPTEMP,8,1200

    !

    MPDATA,C,podlozka,,900

    MPDATA,C,podlozka,,900

    MPDATA,C,podlozka,,900

    MPDATA,C,podlozka,,1470

    MPDATA,C,podlozka,,1470

    MPDATA,C,podlozka,,1000

    MPDATA,C,podlozka,,1100

    MPDATA,C,podlozka,,1100

    !

    MPTEMP,,,,,,,,

    MPTEMP,1,0

    MPTEMP,2,20

    MPTEMP,3,100

    MPTEMP,4,200

    MPTEMP,5,300

    MPTEMP,6,400

    MPTEMP,7,500

    MPTEMP,8,600

    MPTEMP,9,700

    MPTEMP,10,800

    MPTEMP,11,900

    MPTEMP,12,1000

    MPTEMP,13,1100

    MPTEMP,14,1200

    !

    MPDATA,KXX,podlozka,,1.36

    MPDATA,KXX,podlozka,,1.333028

    MPDATA,KXX,podlozka,,1.2297

    MPDATA,KXX,podlozka,,1.1108

    MPDATA,KXX,podlozka,,1.0033

    MPDATA,KXX,podlozka,,0.9072

    MPDATA,KXX,podlozka,,0.8225

    MPDATA,KXX,podlozka,,0.7492

    MPDATA,KXX,podlozka,,0.6873

    MPDATA,KXX,podlozka,,0.6368

    MPDATA,KXX,podlozka,,0.5977

    MPDATA,KXX,podlozka,,0.57

    MPDATA,KXX,podlozka,,0.5537

    MPDATA,KXX,podlozka,,0.5488

    !

    ! ------- MESHING --------

    !

    KSEL,S,LOC,X,0,xxxx-yyyy

    LSLK,S,1

    ASLL,S,1

    VSLA,S,1

    VATT,beton,beton,beton,,beton

    !

    KSEL,S,LOC,x,xxxx-yyyy,xxxx+yyyy

    LSLK,S,1

    ASLL,S,1

    VSLA,S,1

    VATT,podlozka,podlozka,podlozka,,podlozka

    !

    ALLS

    VMESH,ALL

    /PNUM,MAT,1

    /NUMBER,1

    EPLO

    !

    ! ------- DEFINOVANIE VYSTUZE -------

    !

    ocel=632 ! OZNACENIE ELEMENTU <-----

    hustota_ocel =7850!

    poisson_ocel=0.3

    priemer=12! [mm]

    pi=acos(-1)

    plochavystuze=pi*(priemer/1000)*priemer/4000

    !

    ET,ocel,REINF264

    !

    priemer_horna=8! [mm]

    priemer_strm=6! [mm]

    !

    plocha_hore=pi*(priemer_horna/1000)*priemer_horna/4000

    plocha_strm=pi*(priemer_strm/1000)*priemer_strm/4000

    !

    v_hore = ocel+1

    v_strm_z= ocel+2

    v_strm_y= ocel+3

    !

    SECT,ocel,REINF,DISCRETE,prut

    SECDATA,ocel,plochavystuze,EDGo,3, 0.5,0.5,0.5,0.5

    !

    SECT,v_hore,REINF,DISCRETE,vyst_hore

    SECDATA,ocel,plocha_hore,EDGo,3, 0.5,0.5,0.5,0.5

    !

    SECT,v_strm_z,REINF,DISCRETE,strm_z

    SECDATA,ocel,plocha_strm,EDGo,1, 0.5,0.5,0.5,0.5

    !

    SECT,v_strm_y,REINF,DISCRETE,strm_y

    SECDATA,ocel,plocha_strm,EDGo,2, 0.5,0.5,0.5,0.5

    !

    MPTEMP,,,,,,,,

    MPTEMP,1,0

    MPTEMP,2,20

    MPTEMP,3,100

    MPTEMP,4,200

    MPTEMP,5,300

    MPTEMP,6,400

    MPTEMP,7,500

    MPTEMP,8,600

    MPTEMP,9,700

    MPTEMP,10,800

    MPTEMP,11,900

    MPTEMP,12,1000

    MPTEMP,13,1100

    MPTEMP,14,1200

    !

    MPDATA,EX,ocel,,200000E6

    MPDATA,EX,ocel,,200000E6

    MPDATA,EX,ocel,,200000E6

    MPDATA,EX,ocel,,180000E6

    MPDATA,EX,ocel,,160000E6

    MPDATA,EX,ocel,,140000E6

    MPDATA,EX,ocel,,120000E6

    MPDATA,EX,ocel,,62000E6

    MPDATA,EX,ocel,,26000E6

    MPDATA,EX,ocel,,18000E6

    MPDATA,EX,ocel,,14000E6

    MPDATA,EX,ocel,,8000E6

    MPDATA,EX,ocel,,4000E6

    MPDATA,EX,ocel,,1

    !

    MPTEMP,,,,,,,,

    MPTEMP,1,0

    MPTEMP,2,1200

    MPDATA,NUXY,ocel,,poisson_ocel

    MPDATA,NUXY,ocel,,poisson_ocel

    !

    !  Multilinear isotropic hardening; valcované za tepla 

    !

    TB,MISO,ocel,14,14,0

    TBTEMP,0

    TBPT,,0,0

    TBPT,,0.00316,632E6

    TBPT,,0.15,632.01E6

    TBTEMP,20

    TBPT,,0,0E6

    TBPT,,0.00316,632E6

    TBPT,,0.15,632.01E6

    TBTEMP,100

    TBPT,,0,0E6

    TBPT,,0.00316,632E6

    TBPT,,0.15,632.01E6

    TBTEMP,200

    TBPT,,0,0E6

    TBPT,,0.002844,511.92E6

    TBPT,,0.0044,559.14E6

    TBPT,,0.0059,578.27E6

    TBPT,,0.0074,591.87E6

    TBPT,,0.0089,602.33E6

    TBPT,,0.0104,610.61E6

    TBPT,,0.0119,617.2E6

    TBPT,,0.0134,622.39E6

    TBPT,,0.0149,626.35E6

    TBPT,,0.0164,629.22E6

    TBPT,,0.0179,631.06E6

    TBPT,,0.02,632E6

    TBPT,,0.15,632.01E6

    TBTEMP,300

    TBPT,,0,0E6

    TBPT,,0.0024095,385.52E6

    TBPT,,0.004,475.12E6

    TBPT,,0.0055,515.21E6

    TBPT,,0.007,544.04E6

    TBPT,,0.0085,566.37E6

    TBPT,,0.01,584.15E6

    TBPT,,0.0115,598.39E6

    TBPT,,0.013,609.7E6

    TBPT,,0.0145,618.46E6

    TBPT,,0.016,624.92E6

    TBPT,,0.0175,629.26E6

    TBPT,,0.02,632E6

    TBPT,,0.15,632.01E6

    TBTEMP,400

    TBPT,,0,0E6

    TBPT,,0.001896,265.44E6

    TBPT,,0.0035,384.7E6

    TBPT,,0.0051,448.84E6

    TBPT,,0.0067,495.12E6

    TBPT,,0.0083,530.93E6

    TBPT,,0.0099,559.33E6

    TBPT,,0.0115,581.92E6

    TBPT,,0.0131,599.69E6

    TBPT,,0.0147,613.24E6

    TBPT,,0.0163,622.96E6

    TBPT,,0.0179,629.11E6

    TBPT,,0.02,632E6

    TBPT,,0.15,632.01E6

    TBTEMP,500

    TBPT,,0,0E6

    TBPT,,0.001896,227.52E6

    TBPT,,0.0035,317.39E6

    TBPT,,0.0051,363.28E6

    TBPT,,0.0067,396.17E6

    TBPT,,0.0083,421.54E6

    TBPT,,0.0099,441.63E6

    TBPT,,0.0115,457.6E6

    TBPT,,0.0131,470.15E6

    TBPT,,0.0147,479.72E6

    TBPT,,0.0163,486.58E6

    TBPT,,0.0179,490.92E6

    TBPT,,0.02,492.96E6

    TBPT,,0.15,492.97E6

    TBTEMP,600

    TBPT,,0,0E6

    TBPT,,0.00183483870967742,113.76E6

    TBPT,,0.0034,170.38E6

    TBPT,,0.005,202.95E6

    TBPT,,0.0066,226.56E6

    TBPT,,0.0082,244.87E6

    TBPT,,0.0098,259.42E6

    TBPT,,0.0114,271.02E6

    TBPT,,0.013,280.16E6

    TBPT,,0.0146,287.15E6

    TBPT,,0.0162,292.2E6

    TBPT,,0.0178,295.43E6

    TBPT,,0.02,297.04E6

    TBPT,,0.15,297.05E6

    TBTEMP,700

    TBPT,,0,0E6

    TBPT,,0.00170153846153846,44.24E6

    TBPT,,0.0033,72.96E6

    TBPT,,0.0049,90.95E6

    TBPT,,0.0065,104.34E6

    TBPT,,0.0081,114.84E6

    TBPT,,0.0097,123.25E6

    TBPT,,0.0113,129.99E6

    TBPT,,0.0129,135.33E6

    TBPT,,0.0145,139.43E6

    TBPT,,0.0161,142.41E6

    TBPT,,0.0177,144.34E6

    TBPT,,0.02,145.36E6

    TBPT,,0.15,145.37E6

    TBTEMP,800

    TBPT,,0,0E6

    TBPT,,0.00175555555555556,31.6E6

    TBPT,,0.0034,44.73E6

    TBPT,,0.005,51.17E6

    TBPT,,0.0066,55.8E6

    TBPT,,0.0082,59.37E6

    TBPT,,0.0098,62.21E6

    TBPT,,0.0114,64.46E6

    TBPT,,0.013,66.24E6

    TBPT,,0.0146,67.6E6

    TBPT,,0.0162,68.58E6

    TBPT,,0.0178,69.21E6

    TBPT,,0.02,69.52E6

    TBPT,,0.15,69.53E6

    TBTEMP,900

    TBPT,,0,0E6

    TBPT,,0.00180571428571429,25.28E6

    TBPT,,0.0034,30.08E6

    TBPT,,0.005,32.15E6

    TBPT,,0.0066,33.62E6

    TBPT,,0.0082,34.74E6

    TBPT,,0.0098,35.63E6

    TBPT,,0.0114,36.34E6

    TBPT,,0.013,36.9E6

    TBPT,,0.0146,37.32E6

    TBPT,,0.0162,37.63E6

    TBPT,,0.0178,37.82E6

    TBPT,,0.02,37.92E6

    TBPT,,0.15,37.93E6

    TBTEMP,1000

    TBPT,,0,0E6

    TBPT,,0.00158,12.64E6

    TBPT,,0.0032,17.17E6

    TBPT,,0.0048,19.28E6

    TBPT,,0.0064,20.79E6

    TBPT,,0.008,21.95E6

    TBPT,,0.0096,22.87E6

    TBPT,,0.0112,23.6E6

    TBPT,,0.0128,24.18E6

    TBPT,,0.0144,24.63E6

    TBPT,,0.016,24.95E6

    TBPT,,0.0176,25.16E6

    TBPT,,0.02,25.28E6

    TBPT,,0.15,25.29E6

    TBTEMP,1100

    TBPT,,0,0E6

    TBPT,,0.00158,6.32E6

    TBPT,,0.0032,8.59E6

    TBPT,,0.0048,9.64E6

    TBPT,,0.0064,10.39E6

    TBPT,,0.008,10.97E6

    TBPT,,0.0096,11.43E6

    TBPT,,0.0112,11.8E6

    TBPT,,0.0128,12.09E6

    TBPT,,0.0144,12.31E6

    TBPT,,0.016,12.48E6

    TBPT,,0.0176,12.58E6

    TBPT,,0.02,12.64E6

    TBPT,,0.15,12.65E6

    TBTEMP,1200

    TBPT,,0,0

    TBPT,,0.15,0.15

    !

    ! Sucinitel dlzkovej teplotnej roztaznosti ocele

    ! Betonarska ocel

    !

    MPTEMP,,,,,,,,

    MPTEMP,1,0

    MPTEMP,2,20

    MPTEMP,3,100

    MPTEMP,4,200

    MPTEMP,5,300

    MPTEMP,6,400

    MPTEMP,7,500

    MPTEMP,8,600

    MPTEMP,9,750

    MPTEMP,10,860

    MPTEMP,11,1000

    MPTEMP,12,1100

    MPTEMP,13,1200

    !

    MPDATA,THSX,ocel,,0

    MPDATA,THSX,ocel,,0

    MPDATA,THSX,ocel,,0.0009984

    MPDATA,THSX,ocel,,0.0023184

    MPDATA,THSX,ocel,,0.0037184

    MPDATA,THSX,ocel,,0.0051984

    MPDATA,THSX,ocel,,0.0067584

    MPDATA,THSX,ocel,,0.0083984

    MPDATA,THSX,ocel,,0.011

    MPDATA,THSX,ocel,,0.011

    MPDATA,THSX,ocel,,0.0138

    MPDATA,THSX,ocel,,0.0158

    MPDATA,THSX,ocel,,0.0178

    !

    ! hustota ocele

    !

    MPTEMP,,,,,,,,

    MPTEMP,1,0

    MPTEMP,2,1200

    !

    MPDATA,DENS,ocel,,hustota_ocel

    MPDATA,DENS,ocel,,hustota_ocel

    !

    ! Definicia mernej tepelnej kapacity oceli c [J/ (kg*K) ]

    !

    MPTEMP,,,,,,,,

    MPTEMP,1,0

    MPTEMP,2,20

    MPTEMP,3,100

    MPTEMP,4,200

    MPTEMP,5,300

    MPTEMP,6,400

    MPTEMP,7,500

    MPTEMP,8,600

    MPTEMP,9,690

    MPTEMP,10,700

    MPTEMP,11,710

    MPTEMP,12,720

    MPTEMP,13,730

    MPTEMP,14,735

    MPTEMP,15,740

    MPTEMP,16,750

    MPTEMP,17,760

    MPTEMP,18,770

    MPTEMP,19,780

    MPTEMP,20,790

    MPTEMP,21,800

    MPTEMP,22,900

    MPTEMP,23,1200

    !

    MPDATA,C,ocel,,425

    MPDATA,C,ocel,,440

    MPDATA,C,ocel,,488

    MPDATA,C,ocel,,530

    MPDATA,C,ocel,,565

    MPDATA,C,ocel,,606

    MPDATA,C,ocel,,667

    MPDATA,C,ocel,,760

    MPDATA,C,ocel,,937

    MPDATA,C,ocel,,1008

    MPDATA,C,ocel,,1130

    MPDATA,C,ocel,,1388

    MPDATA,C,ocel,,2291

    MPDATA,C,ocel,,5000

    MPDATA,C,ocel,,2525

    MPDATA,C,ocel,,1483

    MPDATA,C,ocel,,1159

    MPDATA,C,ocel,,1002

    MPDATA,C,ocel,,909

    MPDATA,C,ocel,,847

    MPDATA,C,ocel,,803

    MPDATA,C,ocel,,650

    MPDATA,C,ocel,,650

    !

    ! Definicia tepelnej vodivosti lambda [W/ (m*K) ]

    !

    MPTEMP,,,,,,,,

    MPTEMP,1,0

    MPTEMP,2,800

    MPTEMP,3,1200

    !

    MPDATA,KXX,ocel,,54

    MPDATA,KXX,ocel,,27.3

    MPDATA,KXX,ocel,,27.3

    !

    ! ------- OKRAJOVE PODMIENKY -------

    !

    NSEL,S,LOC,X,xxxx

    NSEL,R,LOC,Y,-1*yyyy

    D,ALL, , , , , ,UY,UZ

    !

    NSEL,S,LOC,X,0

    D,ALL, , , , , ,UX

    !

    NSEL,S,LOC,Z,0

    NSEL,A,LOC,Z,zzzz

    D,ALL, , , , , ,UZ

    !

    ALLS

    ! EPLO

  • DadojinDadojin Member
    edited July 2018

    hello Sandeep,

    those above are the setting for material models. After that, elements along the lenght of the "beam" - actually a piece of "slab" as some kind of plane stress is modelled are chosen and the reinf 264 is applied.

    the temperature load is here:

     

  • DadojinDadojin Member
    edited July 2018

    !

    tepl_refer = 20

    cas_horenia=240 ! minut

    !

    ! ------- APLIKOVANIE ZATAZENIA -------

    ! ------- NASTAVENIE RIESICA -------

    !

    /SOLU

    ANTYPE,4! transient analysis

    TRNOPT,FULL! full transient analysis

    KBC,0! ramp loads up and down          

    ! Time stepping        

    TIME,60*cas_horenia! end time for load step      

    AUTOTS,ON! use automatic time stepping  

    ! substep size (seconds) ! minimum value shorter than smallest time change in the table arrays below

    DELTIM,10,5,12.2

    OUTRES,ALL,ALL

    ! LNSRCH,1

    ! NEQIT,100 

    ! LUMPM,0

    TREF,tepl_refer

    !

    ! K?ivka normová : 20 + 345 log10 (8t + 1)

    !

    *DIM,TEPLOTA,TABLE,37,1,1,TIME,,,0

    *SET,TEPLOTA(0,1,1),0

    *SET, TEPLOTA (1 , 0,1), 0

    *SET, TEPLOTA (2 , 0,1), 15

    *SET, TEPLOTA (3 , 0,1), 30

    *SET, TEPLOTA (4 , 0,1), 45

    *SET, TEPLOTA (5 , 0,1), 60

    *SET, TEPLOTA (6 , 0,1), 120

    *SET, TEPLOTA (7 , 0,1), 180

    *SET, TEPLOTA (8 , 0,1), 240

    *SET, TEPLOTA (9 , 0,1), 300

    *SET, TEPLOTA (10 , 0,1), 360

    *SET, TEPLOTA (11 , 0,1), 420

    *SET, TEPLOTA (12 , 0,1), 480

    *SET, TEPLOTA (13 , 0,1), 540

    *SET, TEPLOTA (14 , 0,1), 600

    *SET, TEPLOTA (15 , 0,1), 1200

    *SET, TEPLOTA (16 , 0,1), 1800

    *SET, TEPLOTA (17 , 0,1), 2400

    *SET, TEPLOTA (18 , 0,1), 3000

    *SET, TEPLOTA (19 , 0,1), 3600

    *SET, TEPLOTA (20 , 0,1), 4200

    *SET, TEPLOTA (21 , 0,1), 4800

    *SET, TEPLOTA (22 , 0,1), 5400

    *SET, TEPLOTA (23 , 0,1), 6000

    *SET, TEPLOTA (24 , 0,1), 6600

    *SET, TEPLOTA (25 , 0,1), 7200

    *SET, TEPLOTA (26 , 0,1), 7800

    *SET, TEPLOTA (27 , 0,1), 8400

    *SET, TEPLOTA (28 , 0,1), 9000

    *SET, TEPLOTA (29 , 0,1), 9600

    *SET, TEPLOTA (30 , 0,1), 10200

    *SET, TEPLOTA (31 , 0,1), 10800

    *SET, TEPLOTA (32 , 0,1), 11400

    *SET, TEPLOTA (33 , 0,1), 12000

    *SET, TEPLOTA (34 , 0,1), 12600

    *SET, TEPLOTA (35 , 0,1), 13200

    *SET, TEPLOTA (36 , 0,1), 13800

    *SET, TEPLOTA (37 , 0,1), 14400

    !

    *SET, TEPLOTA (1 , 1,1), 20

    *SET, TEPLOTA (2 , 1,1), 185

    *SET, TEPLOTA (3 , 1,1), 261

    *SET, TEPLOTA (4 , 1,1), 312

    *SET, TEPLOTA (5 , 1,1), 349

    *SET, TEPLOTA (6 , 1,1), 445

    *SET, TEPLOTA (7 , 1,1), 502

    *SET, TEPLOTA (8 , 1,1), 544

    *SET, TEPLOTA (9 , 1,1), 576

    *SET, TEPLOTA (10 , 1,1), 603

    *SET, TEPLOTA (11 , 1,1), 626

    *SET, TEPLOTA (12 , 1,1), 645

    *SET, TEPLOTA (13 , 1,1), 663

    *SET, TEPLOTA (14 , 1,1), 678

    *SET, TEPLOTA (15 , 1,1), 781

    *SET, TEPLOTA (16 , 1,1), 842

    *SET, TEPLOTA (17 , 1,1), 885

    *SET, TEPLOTA (18 , 1,1), 918

    *SET, TEPLOTA (19 , 1,1), 945

    *SET, TEPLOTA (20 , 1,1), 968

    *SET, TEPLOTA (21 , 1,1), 988

    *SET, TEPLOTA (22 , 1,1), 1006

    *SET, TEPLOTA (23 , 1,1), 1022

    *SET, TEPLOTA (24 , 1,1), 1036

    *SET, TEPLOTA (25 , 1,1), 1049

    *SET, TEPLOTA (26 , 1,1), 1061

    *SET, TEPLOTA (27 , 1,1), 1072

    *SET, TEPLOTA (28 , 1,1), 1082

    *SET, TEPLOTA (29 , 1,1), 1092

    *SET, TEPLOTA (30 , 1,1), 1101

    *SET, TEPLOTA (31 , 1,1), 1110

    *SET, TEPLOTA (32 , 1,1), 1118

    *SET, TEPLOTA (33 , 1,1), 1126

    *SET, TEPLOTA (34 , 1,1), 1133

    *SET, TEPLOTA (35 , 1,1), 1140

    *SET, TEPLOTA (36 , 1,1), 1146

    *SET, TEPLOTA (37 , 1,1), 1153

    !

    ! /title,TEPLOTA vs CAS normova krivka

    ! *VPLOT,TEPLOTA(1,0),TEPLOTA(1,1)

    ! /ui,copy,save,png,graph,color,reverse,portrait,yes

    ! /ui,copy,save,png,graph,mono,reverse,portrait,yes

    ! /title

    !

    *DIM,ERKO,TABLE,2,1,1,TIME,,,0

    *SET,ERKO(0,1,1),0

    *SET, ERKO (1 , 0,1), 0

    *SET, ERKO (2 , 0,1), 14400

    !

    *SET, ERKO (1 , 1,1), 1000

    *SET, ERKO (2 , 1,1), 1000

    !

    NSEL,S,LOC,Y,-1*yyyy

    SF,ALL,CONV, %ERKO% , %TEPLOTA%

    ALLS

    !

    *DIM,HORE,TABLE,2,1,1,TIME,,,0

    *SET,HORE(0,1,1),0

    *SET, HORE (1 , 0,1), 0

    *SET, HORE (2 , 0,1), 14400

    !

    *SET, HORE (1 , 1,1), 20

    *SET, HORE (2 , 1,1), 20

    !

    *DIM,ERKO_HORE,TABLE,2,1,1,TIME,,,0

    *SET,ERKO_HORE(0,1,1),0

    *SET, ERKO_HORE (1 , 0,1), 0

    *SET, ERKO_HORE (2 , 0,1), 14400

    !

    *SET, ERKO_HORE (1 , 1,1), 1

    *SET, ERKO_HORE (2 , 1,1), 1

    !

    !NSEL,S,LOC,Y,yyyy

    !SF,ALL,CONV, %ERKO_HORE% , %HORE%

    ALLS

    !

  • DadojinDadojin Member
    edited July 2018

    actually in this case it was not working as well - i mean with the "tref,20" - I needed to put tref,15 - but it does not really matter, as +/- 5°C is almost nothing there.

  • jpasquerelljpasquerell Canonsburg, PAMember
    edited July 2018

    Dadojin:

     

    Try adding an ERESX,NO command prior to the SOLVE. Often the extrapolation of linear results appears as though the stresses are off the stress strain curve. 

     

     

  • jjdoylejjdoyle PittsburghMember
    edited July 2018
    Lots of temperature dependent materials.  What if you take out all temperature dependence (as a test)? 
     aAlso SOLID226 is higher order, SOLID185 is lower order, ...and both element types have different KEYOPT switches.  You should refer to the elements manual docs for each of these element types and go thru each element specification carefully to discover possible differences.  What about SOLID186?
    Also how localized is the result in question?  Perhaps there is enough of a difference in the integration point locations and densities that is causing some interpolation/extrapolation issues. 
    Also, compare the two solver outputs side by side to consider the all the solver specifications and convergence tolerances that might be different.
    Regards,
    John
  • DadojinDadojin Member
    edited July 2018

    thanks for the advices.

    in one temperature only (the normal one = 20°C), solid226 with reinf264 is working as well as solid185 with reinf264.

    however when it comes to temperature loads:

    in case of solid185 I insert load step from previous thermal transient analysis by solid70 elements. That is working for 2 cases - R30 and R60 (30 and 60 minutes of fire) fire resistance according to European standard. The resistance of the structure is a bit higher then what are the caculations by hand according to mentioned standard. In ansys I put the thermal field from the results file as 1st load step, and then putting pressure loads in next steps till it does not convergate.

    However when it comes to R90, the solution falls at approx. time 0,9 (time 1 is the temperature of R90 for the structure) - not talking about R120,180,240 - the structure should be definitelly resisting also these much higher temperatures.

    there is only little difference if i first put load and then temperature, or if i put it in once in one load step, there is smth wrong with the temperature which occurs in the structure in that particular moment somewhere.

    So I tried to use solid226 as well and to put there the temperature load. also working only till some point, however this stresses in reinforcement are weird. When I remove HSD2 from MW material model of concrete, in case of solid226 the solution goes till some 9875 seconds (0second is the initiation of the fire and 240*60 is 240 minutes of fire). after that it was possible to restart the solution form that step, it reached the 240th minute, but the results were this way wrong (stress in reinf.)

    in normal temperature, if in case of solid 185 the enhanced strain is used, it was not working, so full integration is being used since then.

    so far i have not used solid 186 - could give it a try as well

  • jpasquerelljpasquerell Canonsburg, PAMember
    edited July 2018

    The TB,CONC nonlinear material model is only used by SOLID65 (legacy) so it is not used by SOLID226 elements.  

  • jpasquerelljpasquerell Canonsburg, PAMember
    edited July 2018

    See Chapter 4 of the Mechanical APDL Element Reference for a table of Element Support for Material Models.

  • DadojinDadojin Member
    edited July 2018

    sure, but I am using menetrey willam - mat model, which is supported by solid 185 and also 186 and also solid 226 - it is working very well in case of 20°C

  • DadojinDadojin Member
    edited July 2018

    it is defined through concrete,,,MW or smth like that like above

  • jjdoylejjdoyle PittsburghMember
    edited July 2018

    How did you conclude that SOLID226 supports Menetrey-Willam?

    According to the Elements Manual, I only see Extended Drucker Prager concrete option listed under the SOLID226 element chapter. 

     

  • DadojinDadojin Member
    edited July 2018

    well, the same way it is not mentioned in case of solid185, which is perfectly working with MW mat model (see here: https://ansyshelp.ansys.com/account/secured?returnurl=/Views/Secured/corp/v190/ans_elem/elem_matsupp.html )

    so I just tried it, first with temperature 20°C, and the results were almost identical as when I used elements solid185.

  • jjdoylejjdoyle PittsburghMember
    edited July 2018

    Thanks for this feedback.

    There appears to be an inconsistency in our documentation. 

    The Table “Material Model Elements Support” (Chapter 2 of Material Reference Manual) lists the following elements as supporting 'CONCR - Menetry-Willam';

    "PLANE182 (excluding plane stress), PLANE183 (excluding plane stress), SOLID185 (excluding layered structural solid), SOLID186 (excluding layered structural solid), SOLID187CPT212CPT213CPT215CPT216CPT217"

    However the Table “Elements Support for Material Models” (Chapter 4 of Material Reference Manual) does not list Menetry-Willam under SOLID18x and CPTx elements.  It should.

    We will try to fix this inconsistency in future release.

    Thanks again for the feedback.

Sign In or Register to comment.