Nazwa	wskaźniki konsolidacyjne				g(ro)	y	gs	ys	wn	ysr	y'	gsr	gd	yd	e	n	Mo	M	B(betta)	v	Ko
gruntu.	Id	Sr	I L	gr.kon.		gamma
	-----------	---	0,42	C	2,05	20,5	2,68	26,8	20
Gpz	-------------	0,720840480274442	0,42	C	2,05	20,5	2,68	26,8	20	20,7089552238806	10,7089552238806	2,07089552238806	1,70833333333333	3,50208333333333	0,568780487804878	0,362562189054726	18000	30000	0,6	0,32	0,470588235294118
Pr	0,36	0,713706407137064	--------	---------------	2	20	2,65	26,5	22	20,2072378595731	10,2072378595732	2,02072378595731	1,63934426229508	3,27868852459016	0,6165	0,381379523662233	78000	86666,6666666667	0,9	0,25	0,333333333333333
pył	-------------	0,699087823504455	0,56	C	1,95	19,5	2,67	26,7	26	19,6798644551454	9,67986445514536	1,96798644551454	1,54761904761905	3,01785714285714	0,725230769230769	0,420367397895488	13500	22500	0,6	0,32	0,470588235294118
y=	g*g		g=	10	g-przyśpieszenie ziemskie
				y-ciężar objętościowy gruntu							yd=	g*gd		yd-ciężar objętościowy szkieletu gruntowego
ys=	gs*g			ys-ciężar właściwy szkieletu gruntowego
											y'=	ysr - yw		y'-ciężar objętościowy gruntu z uwzględnieniem wyporu wody
Sr=	wn*g / 100*e*gw			Sr-stopień wilgotności gruntu niespoistego
											ysr=	g*gsr		ysr-ciężar objętościowy gruntu przy całkowitym nasyceniu porów wodą
gw=	1			gw-gęstość wody
											yw=	10		yw-ciężar objętościowy wody
e=	gs-gd/gd			e-wskaźnik porowatości
											B=	Mo/M		B-wskaźnik skonsolidowania gruntu(z normy)
n=	gs-gd/gs			n-porowatość
											M=	Mo/B		M-edometryczny moduł ściślwości pierwotnej
gd=	100*g/100+wn			gd-gęstość objętościowa szkieletu gruntowego										Mo-edometryczny moduł ściślwości wtórnej (z normy)
											K=	v/1-v		K-granica parcia bocznego
gsr=	Vs*gs+Vp*gw / V													v-współczynnik Poisson'a(z normy)
gsr=	gd+n*gw			gsr-gęstość objętościowa gruntu przy całkowitym nasyceniu porów wodą
	4.	Rozkład naprężeń pierwotnych w gruncie.
												n
	4,1	Nzprężenia całowite									s=	(suma) yi*hi
												i=1			yi	ciężar objętościowy
		sigma zg1=	22,55	kPa	y1*h1		h1=	1,1							hi	grubość poszczególnych warstw
		sigma zg2=	41	kPa	sigma zg1+y1*h2		h2=	0,9
		sigma zg3=	82	kPa	sigma zg2+y1*h3		h3=	2
		sigma zg4=	142	kPa	sigma zg3+y2*h4		h4=	3
		sigma zg5=	259		sigma zg4+y3*h5		h5=	6
	4.2		Ciśniene porowe U
							Ui=	U(i-1)+hw*yw			hw	wysokość słupa wody
											yw=	10 kN/m^3		ciężar objętościowy wody
			U1(z=1,1)=		0	kPa	hw=	0			U(n-1)	ciśniene porowe warstwy wcześniejszej
			U2(z=0,9)=		0	kPa	hw=	0
			U3(z=2)=		20	kPa	hw=	2
			U4(z=3)=		50	kPa	hw=	3
			U5(z=2)=		70	kPa	hw=	2
	4.3		Naprężenia efektywne
							sig. zgi=		sig. zgi - Ui
			sig'.zg1=	22,55
			sig'.zg2=	41
			sig'.zg3=	62
			sig'.zg4=	92
			sig'.zg5=	189
	5.		Podział warstw na warstewki obliczeniowe
		h1=	0,9
		h2=	1
		h3=	1
		h4=	1
		h5=	1
		h6=	1
		h7=	2
		h8=	2
		h9=	2
		6.		Wyznaczanie rozkładu napężeń w odłożu metodą punktów narożnych.
				prostokąt	I	II	III	IV			warunek
				L	7,5	12,5	12,5	9			L.>=B
				B	5	5	9	7,5
				L / B	1,5	2,5	1,38888888888889	1,2
			Tabela - odciążenie podłoża wykopem.
			I		II		III		IV		odciążenie podłoża
	zx=teta		L/B=1,5		L/B=2,5		L/B=1,39		L/B=1,
	b-sigma	z	Z / B	zx nI	Z / B	zx nII	Z / B	zx nIII	Z / B	zx nIV	suma zx ni	b""zgi		b""zg=	qw*zxn
		0	0	0,25	0	0,25	0	0,25	0	0,25	1	22,55
		0,9	0,18	0,2493	0,18	0,2494	0,1	0,2499	0,12	0,2497	0,9931	22,394405		qw=	yd*D
		1,9	0,38	0,2441	0,38	0,2449	0,211111111111111	0,2488	0,253333333333333	0,2477	0,9738	21,95919		qw=	22,55
		2,9	0,58	0,2324	0,58	0,2349	0,322222222222222	0,2461	0,386666666666667	0,2427	0,9391	20,6602
		3,9	0,78	0,2154	0,78	0,2206	0,433333333333333	0,2412	0,52	0,2341	0,8916	19,6152
		4,9	0,98	0,1957	0,98	0,2041	0,544444444444444	0,2342	0,653333333333333	0,2226	0,8363	18,3986
		5,9	1,18	0,1754	1,18	0,187	0,655555555555556	0,2253	0,786666666666667	0,2089	0,7594	16,7068
		7,9	1,58	0,1382	1,58	0,1553	0,877777777777778	0,2039	1,05333333333333	0,1791	0,6467	13,871715
		9,9	1,98	0,1084	1,98	0,1286	1,1	0,1808	1,32	0,1503	0,5454	11,69883
		11,9	2,38	0,0857	2,38	0,1069	1,32222222222222	0,1579	1,58666666666667	0,1252	0,4757	10,203765
				A
			7.	Wyznaczanie naprężeń od obciążeń zewnętrznych.
			7.1.	Metoda punktów środkowych.
			Prostokąty
											L	B	L/B=
										I	5,5	2	2,75
					I					II	2	1,5	1,33333333333333
							A	II		III =III-VI+V	7	5,5	1,27272727272727
										IV	2	1,5	1,33333333333333
							V	IV		V	2	1,5	1,33333333333333
					III					VI=V+VI	7	1,5	4,66666666666667
							VI
				I		II		III		IV		V		VI		obciążenie podłoża
				L/B=2,75		L/B=1,33		L/B=1,273		L/B=1,33		L/B=1,33		L/B=4,667					b1zq=	Q*zxn		D-	głębokość wyopu
			z	Z / B	zx nI	Z / B	zx nII	Z / B	zx nIII	Z / B	zx nIV	Z / B	zx nV	Z / B	zx nVI	suma zx ni	b1 zqi		Q=	400		D=	1,1
	zx=teta		0	0	0,25	0	0,25	0	0,25	0	0,25	0	0,25	0	0,25	1	400					zxn-	współ. zanikania naprężeń(z normy)
	b-sigma		0,9	0,45	0,2421	0,6	0,2294	0,163636363636364	0,2494	0,6	0,2294	0,6	0,2294	0,6	0,2342	0,9551	382,04					qw-	ciężar wykopu
			1,9	0,95	0,2072	1,26666666666667	0,1614	0,345454545454545	0,2449	1,26666666666667	0,1614	1,26666666666667	0,1614	1,26666666666667	0,1833	0,7968	318,72					yw-	ciężar objętościowy wykopanego gruntu
			2,9	1,45	0,1665	1,93333333333333	0,1051	0,527272727272727	0,2344	1,93333333333333	0,1051	1,93333333333333	0,1051	1,93333333333333	0,1396	0,6456	258,24
			3,9	1,95	0,1326	2,6	0,0701	0,709090909090909	0,2186	2,6	0,0701	2,6	0,0701	2,6	0,109	0,5303	212,12
			4,9	2,45	0,1064	3,26666666666667	0,049	0,890909090909091	0,1997	3,26666666666667	0,049	3,26666666666667	0,049	3,26666666666667	0,0874	0,4425	177
			5,9	2,95	0,0863	3,93333333333333	0,0358	1,07272727272727	0,1797	3,93333333333333	0,0358	3,93333333333333	0,0358	3,93333333333333	0,0715	0,3733	149,32		qw=	yd*D
			7,9	3,95	0,0588	5,26666666666667	0,0212	1,43636363636364	0,1423	5,26666666666667	0,0212	5,26666666666667	0,0212	5,26666666666667	0,05	0,2723	108,92		qw=	22,55
			9,9	4,95	0,0419	6,6	0,0139	1,8	0,1117	6,6	0,0139	6,6	0,0139	6,6	0,0365	0,204	81,6
			11,9	5,95	0,031	7,93333333333333	0,0098	2,16363636363636	0,0884	7,93333333333333	0,0098	7,93333333333333	0,0098	7,93333333333333	0,0276	0,1568	62,72
	7.1	Metoda punktów śrdkowych.
							Q=	400		b zqi=	q * zxni
			L=4	A
							L/B=	1,33333333333333		zxn=	f(L/B;Z/B)
										zx ni=	4 * zx n
			Z [m2]	Z/B	ex ni	b'zqi
			0	0	0,25	100
			0,9	0,3	0,2467	98,68
			1,9	0,633333333333333	0,2455	98,2
			2,9	0,966666666666667	0,2455	98,2
			3,9	1,3	0,2455	98,2
			4,9	1,63333333333333	0,2455	98,2
			5,9	1,96666666666667	0,2455	98,2
			7,9	2,63333333333333	0,2236	89,44
			9,9	3,3	0,2236	89,44
			11,9	3,96666666666667	0,2236	89,44
	7.2	Wyznaczane naprężeń spowodowanych obciążeniem od sąsiada.
														b zqo= 3*Q*z^3/2pi*Ro^5
														Ro=	pierw.(r^2+z^2)
														r=	8,732
														q=	340
														R=	2
														Q=	4272,576
			r	z	Ro	b2 zqo
			8,732	0	0	0
			8,732	0,9	8,77825859723898	0,028530924414516
			8,732	1,9	8,93632049559549	0,245525887658629
			8,732	2,9	9,2009686446591	0,754495072159113
			8,732	3,9	9,56335840591578	1,51276625369578
			8,732	4,9	10,0128829015424	2,38463943380126
			8,732	5,9	10,5383976011536	3,22339332739155
			8,732	7,9	11,7753056860533	4,44273582921008
			8,732	9,9	13,2006751342498	4,93804490264605
			8,732	11,9	14,7600075880739	4,90724959470127
	8.		Obliczanie rozkładu naprężeń wtórnych b zs i dodatkowych b zd.
			b"""zg	naprężenia od odciążenia podłoża
			b zs=	b"""zgi						z	b1 zqi	b2 zqo	b zqc	b"""zgi	b zd
			b zqc=	b1zqi+b2zqo						0	400	0	400	22,55	377,45
			b zd=	bzqc - bzg						0,9	382,04	0,0285	382,0685	22,3944	359,6741
			b zd-	naprężenia dodatkowe						1,9	318,72	0,2455	318,9655	21,9592	297,0063
										2,9	258,24	0,7545	258,9945	20,6602	238,3343
										3,9	212,12	1,5128	213,6328	19,6152	194,0176
										4,9	177	2,3846	179,3846	18,3986	160,986
										5,9	149,32	3,2234	152,5434	16,7068	135,8366
										7,9	108,92	4,4427	113,3627	13,8717	99,491
										9,9	81,6	4,938	86,538	11,6988	74,8392
										11,9	62,72	4,9072	67,6272	10,2038	57,4234
	9.		Obliczenie osiadania pod punktem A
			Określenie głębokości strefy aktywnej
			b z max.d	<=	b"""z max.g
				S=	S' + S"
				S'=	b zd*hi/Mo		osiadanie pierwotne warstwy
				S"=	lamb.*(b zs*hi/M)		osiadanie wturne warstwy
				lambda= 0 lub 1