perete 1 b-c

10

Click here to load reader

Upload: maria-acceleanu

Post on 04-Jan-2016

216 views

Category:

Documents


2 download

TRANSCRIPT

Page 1: PERETE 1 B-C

1.Calculul rezistentei la Forta axiala

280.00 cm A totala= 2.00 seism pozitiv --> ETAJ 1-seism Y negativ zona comprimata stanga

25.00 cm A totala= 2.00 seism negativ <-- 3.20 1.68 9.25 300 1750 615.75

11.20 73.75 KNm289.83 KN

coef de echivalenta 5.49

latimea talpii echivalente 137.30 cm coordonata y a C.G. al peretelui 169.13 cm

seism pozitiv

l talpa parter etaj 1

comprimata 625 625 aria comprimata Ac= 2024.55 cm^2intinsa 0 0

0.93 cmseism negativ

l talpa parter etaj 1 lungimea comprimata de pe talpa Xc= 80.98 cm

comprimata 0 0 lung. Totala a zonei comprimate x= 80.98 cmintinsa 625 625 y`=X/2

coordonata y a C.G. al zonei comprimate y`= 40.49 cm

PARTER-seism Y negativ zona comprimata stanga distanta C.G.zonei comprimate-C.G.perete 128.64 cm

seism pozitiv --> seism <-- negativN (kN) M(kNm) N (kN) M(kNm) 181.33 KNm Mcap. Al sect de beton 372.83 Knm

ETAJ 1

27.38 4.27 27.12 4.29 499.48 KN72.61 3.36 71.95 3.36 Mcap.al armaturii din stalpisor 323.27 Knm

117.84 2.46 116.78 2.42 coordonata y a C.G. al peretelui 165.62 cm

PARTER

150.69 2.14 149.59 2.22 Mcap 696.10 Knm202.12 4.34 201.90 4.46 aria comprimata Ac= 3489.01 cm^2253.54 6.54 254.21 6.69 59.21 kN

*pentru valorile de la mijlocul nivelului am mediat extremitatiile lungimea comprimata de pe talpa Xc= 139.56 cm <25cmca valoare de calcul pentru forta taietoare se va folosi maximul dintre cele doua valori

*M si N sunt trecute in modul lung. Totala a zonei comprimate x= 139.56 cmy`=X/2

seism --> seism <-- seism --> seism <-- coordonata y a C.G. al zonei comprimate y`= 69.78 cm

ehi ei

ETAJ10.16 0.16 0.17 0.17 distanta C.G.zonei comprimate-C.G.perete 95.84 cm

0.05 0.05 0.06 0.06 3.1Rezistenta de proiectare la lunecare in rost orizontal

0.02 0.02 0.03 0.03 Mcap. Al sect de beton 478.70 Knm

PARTER0.01 0.01 0.02 0.02 seism negativ seism pozitiv

0.02 0.02 0.03 0.03 Mcap.al armaturii din stalpisor 323.27 Knm Ned= 108.73 Ned= 109.14 kN

0.03 0.03 0.04 0.04 Med= 56.13 Med= 53.60 kNm

em/t Mcap 801.97 Knm Ved= 31.29 Ved= 26.50 kN

seism --> seism <-- seism -->seism <--etaj 1 0.22 0.22 0.49 0.49 59.21 kN e=Med/Ned Aper= 2139889.56 mmpParter 0.12 0.22 0.69 0.49 ca valoare de calcul pentru forta taietoare se va folosi maximul dintre cele doua valori en= 0.52 m

ep= 0.49 m

lcn=1.5lw-3e= 1.08 m

hef/texcentricitatea relativa emk/t lcp=1.5lw-3e= 1.15 m

0.05 0.10 0.15 0.20 0.25 0.30 3.1Rezistenta de proiectare la lunecare in rost orizontal lad=lcn+lcp-lw= 0.48

5.00 0.89 0.79 0.69 0.59 0.49 0.39 59.21 kN6.00 0.88 0.78 0.68 0.58 0.48 0.38 seism negativ seism pozitiv aria totala a zonei comprimate Azc= 80.98 cm^2

7.00 0.88 0.77 0.67 0.57 0.47 0.37 Ned= 197.09 Ned= 195.70 kN calc.rezistentei unitare caract la forfecare 0.20 N/mm^2

8.00 0.86 0.76 0.66 0.56 0.45 0.35 Med= 157.37 Med= 145.51 kNm calc.rezistentei unitare de proiectare la forfecare 0.11 N/mm^29.00 0.85 0.75 0.65 0.54 0.44 0.34 Ved= 64.00 Ved= 55.07 kN factor de participare la armaturii 0.35

10.00 0.84 0.73 0.63 0.53 0.42 0.32 108.73 kN

11.00 0.82 0.72 0.61 0.51 0.40 0.30 e=Med/Ned Aper= 2000485 mmp11.20 0.85 0.74 0.64 0.54 0.43 0.33 en= 0.80 m rez.de pr. La forfecare a ZNA Vrdl*= 56.07 KN

12.00 0.80 0.70 0.59 0.49 0.38 0.28 ep= 0.74 m

13.00 0.79 0.68 0.57 0.47 0.36 0.26 lcn=1.5lw-3e= 0.60 m rez.de pr. La forfecare a arm stalpisor comp 129.31 KN

14.00 0.77 0.66 0.55 0.45 0.34 0.24 lcp=1.5lw-3e= 0.77 m

15.00 0.75 0.64 0.53 0.42 0.32 0.22 lad=lcn+lcp-lw= -0.63 m deci lad=0 rez.de pr. a stalpisorului comprimat 96.25

16.00 0.72 0.61 0.51 0.40 0.30 0.20 59.21 kN rez. de proiectare la lunecare in rost orizontal Vrdl= 281.63 KN17.00 0.70 0.59 0.48 0.38 0.28 0.18 aria totala a zonei comprimate Azc= 3489.01 cm^2

18.00 0.68 0.57 0.46 0.35 0.25 0.16 calc.rezistentei unitare caract la forfecare 0.20 N/mm^2 3.2Rezistenta de proiectare la cedare pe sectiune inclinata a peretilor din zidarie confinata

19.00 0.65 0.54 0.44 0.33 0.23 0.14 calc.rezistentei unitare de proiectare la forfecare 0.11 N/mm^2 108.73 kN

20.00 0.63 0.52 0.41 0.31 0.21 0.13 Tabel CR6 factor de participare la armaturii 0.35 0.05 N/mm^2260.81 kN rez. unitara caract. de cedare pe sectiune inclinata fvk,i= 0.08 N/mm^2

Determinarea factorului de reducere Comparatie Nrd-Nedcapacitate-cerinta rez.de pr. La forfecare a ZNA Vrdl*= 104.32 KN rez.de pr. La forfecare a ZNA Vrdi*= 200.75 KN

seism --> seism<-- seism --> seism<-- seism --> seism<-- rez.de pr. La forfecare a arm stalpisor comp 129.31 KN rez.de pr. la forfecare a arm stalpisor comprimat 64.65 KN

Φ Nrd Ned

ETAJ1-0.32 0.99 -237.48 726.97 27.38 27.12 rez.de pr. a stalpisorului comprimat 192.50 rez.de pr. a stalpisorului comprimat 192.50

0.49 0.49 359.77 358.03 72.61 71.95 rez. de proiectare la lunecare in rost orizontal Vrdl= 426.13 KN Vrdi= 457.91 KN

0.76 1.00 558.77 735.07 117.84 116.78

PARTER0.81 1.00 598.11 735.42 150.69 149.59 3.2Rezistenta de proiectare la cedare pe sectiune inclinata a peretilor din zidarie confinata Vrd=min (Vrdl ,Vrdi) 281.63 kN

0.69 0.49 510.98 358.03 202.12 201.90 260.81 kN daca Ved>Vrd se va propune armare transversala

0.72 1.00 529.77 734.74 253.54 254.21 0.13 N/mm^2 Vrd3=0.8*lw*Asw/s*fydrez. unitara caract. de cedare pe sectiune inclinata fvk,i= 0.11 N/mm^2

3.3Rezistenta la lunecare verticala asociata incovoierii peretelui rez.de pr. La forfecare a ZNA Vrdi*= 141.06 KN

73.68 kN

rez.de pr. la forfecare a arm stalpisor comprimat 258.62 KN

3.4Rezistenta de proiectare a peretelui supus la inc. perpendicualara pe planul median

rez.de pr. a stalpisorului comprimat 192.50

Vrdi= 592.18 KN

Vrd=min (Vrdl ,Vrdi) 426.13 kN 10416666.67daca Ved>Vrd se va propune armare transversala fxk1 =0.24/1.9= 0.13Vrd3=0.8*lw*Asw/s*fyd fxk2 =0.48/1.9= 0.25

1315789.47

2631578.9511

ETAJ 1-seism Y pozitiv zona comprimata dreapta

70.47 KNm238.59 KN

coordonata y a C.G. al peretelui 80.86 cm

aria comprimata Ac= 1666.62 cm^2

lungimea comprimata de pe talpa Xc= 2.67 cmlung. Totala a zonei comprimate x= 2.67 cmy`=X/2coordonata y a C.G. al zonei comprimate y`= 1.33 cm

distanta C.G.zonei comprimate-C.G.perete 79.53 cm

5Mcap. Al sect de beton 189.74 Knm

PARTER-seism Y pozitiv zona comprimata dreapta Mcap.al armaturii din stalpisor 646.54 Knm

169.19 KNm Mcap 836.28 Knm497.03 KN

59.21 kNcoordonata y a C.G. al peretelui 34.38 cm

ca valoare de calcul pentru forta taietoare se va folosi maximul dintre cele doua valori

aria comprimata Ac= 3471.90 cm^2

lungimea comprimata de pe talpa Xc= 5.56 cm <25cm3.00 lung. Totala a zonei comprimate x= 5.56 cm

y`=X/2 3.1Rezistenta de proiectare la lunecare in rost orizontalcoordonata y a C.G. al zonei comprimate y`= 2.78 cm

seism negativ seism pozitiv

distanta C.G.zonei comprimate-C.G.perete 31.60 cm Ned= 108.73 Ned= 109.14 kN

Med= 56.13 Med= 53.60 kNm

Mcap. Al sect de beton 157.07 Knm Ved= 31.29 Ved= 26.50 kN

Mcap.al armaturii din stalpisor 323.27 Knm e=Med/Ned Aper= 2280625 mmp

en= 0.52 mMcap 480.34 Knm ep= 0.49 m

lcn=1.5lw-3e= 1.08 m59.21 kN lcp=1.5lw-3e= 1.15 m

ca valoare de calcul pentru forta taietoare se va folosi maximul dintre cele doua valorilad=lcn+lcp-lw= 0.48

59.21 kNaria totala a zonei comprimate Azc= 2.67 cm^2

calc.rezistentei unitare caract la forfecare 0.20 N/mm^2

3.1Rezistenta de proiectare la lunecare in rost orizontal calc.rezistentei unitare de proiectare la forfecare 0.11 N/mm^2factor de participare la armaturii 0.35

seism negativ seism pozitiv 126.10 kNNed= 197.09 Ned= 195.70 kN

Med= 157.37 Med= 145.51 kNm rez.de pr. La forfecare a ZNA Vrdl*= 63.02 KNVed= 64.00 Ved= 55.07 kN

rez.de pr. La forfecare a arm stalpisor comp 64.65 KNe=Med/Ned Aper= 2155625 mmp

en= 0.80 m rez.de pr. a stalpisorului comprimat 96.25

ep= 0.74 m rez. de proiectare la lunecare in rost orizontal Vrdl= 223.92 KNlcn=1.5lw-3e= 0.23 mlcp=1.5lw-3e= 0.39 m 3.2Rezistenta de proiectare la cedare pe sectiune inclinata a peretilor din zidarie confinata

lad=lcn+lcp-lw= -2.88 m lad=0 126.10 kN

59.21 kN 0.06 N/mm^2aria totala a zonei comprimate Azc= 3471.90 cm^2 rez. unitara caract. de cedare pe sectiune inclinata fvk,i= 0.08 N/mm^2

calc.rezistentei unitare caract la forfecare 0.20 N/mm^2

calc.rezistentei unitare de proiectare la forfecare 0.11 N/mm^2 rez.de pr. La forfecare a ZNA Vrdi*= 209.46 KN

factor de participare la armaturii 0.35

195.70 kN rez.de pr. la forfecare a arm stalpisor comprimat 64.65 KN

rez.de pr. La forfecare a ZNA Vrdl*= 78.28 KN rez.de pr. a stalpisorului comprimat 192.50

Vrdi= 466.61 KN

rez.de pr. La forfecare a arm stalpisor comp 64.65 KN

Vrd=min (Vrdl ,Vrdi) 223.92 kN

rez.de pr. a stalpisorului comprimat 96.25 daca Ved>Vrd se va propune armare transversalarez. de proiectare la lunecare in rost orizontal Vrdl= 239.18 KN Vrd3=0.8*lw*Asw/s*fyd

3.2Rezistenta de proiectare la cedare pe sectiune inclinata a peretilor din zidarie confinata 3.3Rezistenta la lunecare verticala asociata incovoierii peretelui 195.70 kN

0.09 N/mm^2 73.68 kNrez. unitara caract. de cedare pe sectiune inclinata fvk,i= 0.10 N/mm^2

3.4Rezistenta de proiectare a peretelui supus la inc. perpendicualara pe planul medianrez.de pr. La forfecare a ZNA Vrdi*= 171.72 KN

rez.de pr. la forfecare a arm stalpisor comprimat 64.65 KN

10416666.67

rez.de pr. a stalpisorului comprimat 192.50 fxk1 =0.24/1.9= 0.13Vrdi= 428.87 KN fxk2 =0.48/1.9= 0.25

1315789.47

Vrd=min (Vrdl ,Vrdi) 239.18 kN 2631578.95daca Ved>Vrd se va propune armare transversala 11Vrd3=0.8*lw*Asw/s*fyd

8

10 17

h etaj m2 fk fd fcd fyd linima Aa

t perete m2

het/t= bech=25*nech Med=

Nrd=Φi(m)*Ai*fd Nrd=rezistenta la compresiune axiala a peretilor din ZC nech= fcd/fd Ned=

Φi(m)=1-2ei/t nech=

ei=e0i+ehi+ea ≥0.05*t=0.0125m bech= yCG=

e0i =excentricitate incarcarii variabile Ac=Ned/(0.85*fd)

ea=excentricitatea accidentala

ehi=excentricitatea datorata fortelor Xc=Ac/Ltalpa

ea=max(t/30;het/300)=

ehi=Mi/Ni

em=2/3*ei0+ehm+ea 2.Calculul rezistentei de proiectare la compresiune si incovoiere in planul peretilor structurali din zidarie confinata Mrdehm=excentricitatea datorata efectului incarcarilor orizontale, in sect. de la mijl. peretelui yzc=YCG-y`

ehm=Mm/Nm yzc=

Mrd(ZNA)=Ned*yzc

Med= Mrd(ZNA)=

Nsus Msus Nsus Msus Ned= Mrd(As)=As*fyd*(Larm1+Larm

2)

Nmijl Mmijl Nmijl Mmijl Mrd(As)=

Njos Mjos Njos Mjos yCG= Mrd=Mrd(ZNA)+Mrd(As)

Nsus Msus Nsus Msus Ac=Ned/(0.85*fd) Mrd=

Nmijl Mmijl Nmijl Mmijl

Njos Mjos Njos Mjos Xc=Ac/Ltalpa Ved≤lpan*t*fvd=

yzc=YCG-y`3.Calculul rezistentei de proiectare la forta taietoare a peretilor

structurali din zidarie confinata Vrdehi sus ehi sus ei sus ei sus yzc=

ehm mijl ehi mijl em mijl em mijl Mrd(ZNA)=Ned*yzc

ehi jos ehi jos ei jos ei jos Mrd(ZNA)=

ehi sus ehi sus ei sus ei sus Mrd(As)=As*fyd*(Larm1+Larm

2)

ehm mijl ehm mijl em mijl em mijl Mrd(As)=

ehi jos ehi jos ei jos ei jos Mrd=Mrd(ZNA)+Mrd(As)

Φi mijl Mrd=

Ved≤lpan*t*fvd=

Determinarea lui Φi mijl 3.Calculul rezistentei de proiectare la forta taietoare a peretilor structurali din zidarie confinata Vrd

Vedcalcul

fvk0=

fvd0=λc=

Ned*=Ned+0.8Ved*hpan/lpan Ned*=

Vrdl*=1/γM*fvk0*t*lad+0.4*Ned*

Vrd2=λc*Aasc*fyd

Vrd2=

Vrsc=Astalpisor*fcvd

Vrsc=

Vedcalcul

fvk0=

fvd0= Ned*=Ned+0.8Ved*hpan/lpan Ned*=

λc= Ϭ0d*=Ned*/Aw Ϭ0d*=

Ned*=Ned+0.8Ved*hpan/lpan Ned*=

Vrdl*=1/γM*fvk0*t*lad+0.4*Ned* Vrdi*=Aw/b*fvd

Vrd2=λc*Aasc*fyd Vrd2=λc*Aasc*fyd

Vrd2= Vrd2=

Vrsc=Astalpisor*fcvd Vrsc=Astalpisor*fcvd

Φi sus Ni sus Vrsc= Vrsc=

Φi mijl Ni mijl

Φi jos Ni jos

Φi sus Ni sus

Φi mijl Ni mijl Ned*=Ned+0.8Ved*hpan/lpan Ned*=

Φi jos Ni jos Ϭ0d*=Ned*/Aw Ϭ0d*=

Vrdi*=Aw/b*fvd

Vrd2=λc*Aasc*fyd VLdh=hetaj*t*fvk0/γM=

Vrd2=

Vrsc=Astalpisor*fcvd

Vrsc=

MRxd1=ww*fxd1

MRxd2=ww*fxd2

Ww=1000*t2/6 mm3

MRxd1=

MRxd2=

Med=Ned=

yCG=

Ac=Ned/(0.85*fd)

Xc=Ac/Ltalpa

yzc=YCG-y`yzc=

Mrd(ZNA)=Ned*yzc

2.Calculul rezistentei de proiectare la compresiune si incovoiere in planul peretilor structurali din zidarie confinata Mrd

Mrd(ZNA)=

Mrd(As)=As*fyd*(Larm1+Larm

2)Mrd(As)=

Mrd=Mrd(ZNA)+Mrd(As)

Med= Mrd=Ned=

Ved≤lpan*t*fvd=yCG=

Ac=Ned/(0.85*fd)

Xc=Ac/Ltalpa

3.Calculul rezistentei de proiectare la forta taietoare a peretilor structurali din zidarie confinata Vrd

yzc=YCG-y`yzc=

Mrd(ZNA)=Ned*yzc

Mrd(ZNA)=

Mrd(As)=As*fyd*(Larm1+Larm

2)Mrd(As)=

Mrd=Mrd(ZNA)+Mrd(As)

Mrd=

Ved≤lpan*t*fvd=

Vedcalcul

3.Calculul rezistentei de proiectare la forta taietoare a peretilor structurali din zidarie confinata Vrd fvk0=

fvd0=λc=

Ned*=Ned+0.8Ved*hpan/lpan Ned*=

Vrdl*=1/γM*fvk0*t*lad+0.4*Ned*

Vrd2=λc*Aasc*fyd

Vrd2=

Vrsc=Astalpisor*fcvd

Vrsc=

Ned*=Ned+0.8Ved*hpan/lpan Ned*=

Vedcalcul Ϭ0d*=Ned*/Aw Ϭ0d*=

fvk0= Vrdi*=Aw/b*fvd

fvd0=

λc= Vrd2=λc*Aasc*fyd

Ned*=Ned+0.8Ved*hpan/lpan Ned*= Vrd2=

Vrdl*=1/γM*fvk0*t*lad+0.4*Ned* Vrsc=Astalpisor*fcvd

Vrsc=

Vrd2=λc*Aasc*fyd

Vrd2=

Vrsc=Astalpisor*fcvd

Vrsc=

Ned*=Ned+0.8Ved*hpan/lpan Ned*=

Ϭ0d*=Ned*/Aw Ϭ0d*= VLdh=hetaj*t*fvk0/γM=

Vrdi*=Aw/b*fvd

Vrd2=λc*Aasc*fyd MRxd1=ww*fxd1

Vrd2= MRxd2=ww*fxd2

Vrsc=Astalpisor*fcvd Ww=1000*t2/6 mm3

Vrsc=

MRxd1=

MRxd2=

Page 2: PERETE 1 B-C

Story Pier Load Loc P V2 V3 T M2 M3STORY4 TALPA MA GSYPMP Bottom -277.94 0.06 1.64 7.442 -0.514 2.413STORY3 TALPA MA GSYPMP Bottom -571.49 81.01 3.2 9.647 0.419 143.438STORY2 TALPA MA GSYPMP Bottom -880.03 136.2 4.22 9.167 2.192 337.481STORY1 TALPA MA GSYPMP Bottom -1194.94 152.28 5.91 1.524 5.687 601.353

Page 3: PERETE 1 B-C

Story Pier Load Loc P V2 V3 T M2 M3STORY4 TMICA GSYNMP Bottom -418.51 -6.37 3.69 3.629 7.033 16.559STORY3 TMICA GSYNMP Bottom -835.73 -84.44 0.15 2.183 11.257 -162.97STORY2 TMICA GSYNMP Bottom -1261.8 -140.51 -2.91 0.945 15.983 -554.007STORY1 TMICA GSYNMP Bottom -1664.95 -160.25 -6.22 6.495 14.373 -946.428

Page 4: PERETE 1 B-C

Predimensionare Fundatiepef<pconvpconv 260 kPA

Perete ax A L 12.4 m Perete ax 1 L 11.7 m

Aaf 14.28 Aaf 35

N supras 638.3342 kN N supras 898.6365 kNN 702.1676 kN N 988.5002 kNB 0.217794 m B 0.324951 m

Perete ax B L 15.1 m Perete ax 2 L 11.7 m

Aaf 31 Aaf 55.55kN

N supras 964.6959 kN N supras 1181.536 kNN 1061.165 kN N 1299.69 kNB 0.270292 m B 0.427248 m

Perete ax C L 15.1 m Perete ax 3 L 11.7 m

Aaf 30 Aaf 39.6

N supras 950.9295 kN N supras 2062.809 kNN 1046.022 kN N 2269.09 kNB 0.266435 m B 0.74592 m

Perete ax D L 15.1 m Perete ax 4 L 8.85 m

Aaf 42.5 Aaf 12.35N subsol 245.1798 kN

N supras 1123.01 kN N supras 1013.907 kNN 1235.31 kN N 1384.995 kNB 0.314649 m B 0.60191 m

Perete ax E L 15.1 m

Aaf 28.75

N supras 933.7215 kNN 1027.094 kNB 0.261613 m

Bfundatie= 0.74592 mBefectiv 70cm

m2 m2

m2 m2

m2 m2

m2 m2

m2

Page 5: PERETE 1 B-C

Calcul armatura necesara scarab d1 fyd z1 d2 z2

1000 85 300 76.5 115 103.5

Nr.crt Med Asnec Asef Mcap diametru1 2.61 84.05797 251.3274 7.803716 5fi82 5.55 178.744 251.3274 7.803716 5fi83 5.25 169.0821 251.3274 7.803716 5fi84 0.6 19.32367 251.3274 7.803716 5fi85 2 64.41224 251.3274 7.803716 5fi86 0.5 21.78649 251.3274 5.767964 5fi87 0.3 13.0719 251.3274 5.767964 5fi88 1.2 52.28758 251.3274 5.767964 5fi8

Asnec=Med/(z*fyd)

Page 6: PERETE 1 B-C

Verificarea deplasarilor laterale

Tc c1 0.136 2

q= 3.125Calcul drift directia Y

Story

STORY4 0.0000365 0.000073 0.00011 0.00046 0.0015 0.004STORY3 0.000044 0.000088 0.00014 0.00055 0.0015 0.004STORY2 0.0000495 0.000099 0.00015 0.00062 0.0015 0.004STORY1 0.0000395 0.000079 0.00012 0.00049 0.0015 0.004

Calcul drift directia XStory

STORY4 0.000022 0.000044 0.00007 0.00028 0.0015 0.004STORY3 0.0000345 0.000069 0.00011 0.00043 0.0015 0.004STORY2 0.0000415 0.000083 0.00013 0.00052 0.0015 0.004STORY1 0.000035 0.00007 0.00011 0.00044 0.0015 0.004

T SLU

dr SLS = g*q*dre SLS

dr SLU = c*q*dre SLU

dre Y SLS dre Y SLU dr Y SLS dr Y SLU dr,a SLS dr,a SLU

dre X SLS dre X SLU dr X SLS dr X SLU dr,a SLS dr,a SLU

STORY4

STORY3

STORY2

STORY1

0.00000

0.00002

0.00004

0.00006

0.00008

0.00010

0.00012

0.00014

0.00016

0.00018

dr X SLSdr Y SLS

STORY4

STORY3

STORY2

STORY1

0.00000 0.00010 0.00020 0.00030 0.00040 0.00050 0.00060 0.00070

dr X SLUdr Y SLU