Learning Forum

[Dadojin]

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]

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]

!

!

---

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]

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]

!

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]

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]

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]

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]

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]

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

[jpasquerell]

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

[Dadojin]

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]

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

[John Doyle]

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]

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]

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), SOLID187, CPT212, CPT213, CPT215, CPT216, CPT217"

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.