General Mechanical

General Mechanical

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

    • Dadojin
      Subscriber
      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?
    • Sandeep Medikonda
      Ansys Employee

      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

    • Dadojin
      Subscriber

      !


      !


      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

    • Dadojin
      Subscriber

      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:


       

    • Dadojin
      Subscriber

      !


      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


      !

    • Dadojin
      Subscriber

      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.

    • jpasquerell
      Ansys Employee

      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. 


       


       

    • John Doyle
      Ansys Employee
      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
    • Dadojin
      Subscriber

      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

    • jpasquerell
      Ansys Employee

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

    • jpasquerell
      Ansys Employee

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

    • Dadojin
      Subscriber

      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

    • Dadojin
      Subscriber

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

    • John Doyle
      Ansys Employee

      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. 


       

    • Dadojin
      Subscriber

      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.

    • John Doyle
      Ansys Employee

      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.

Viewing 15 reply threads
  • You must be logged in to reply to this topic.