TÍNH TOÁN ỨNG DỤNG TRONG KTDK Báo cáo tính toán ud bao đẹp

t i ng dụng EXCEL và MATLAB kiểm tra toán tính hệ số an to n mái dốc GVHD h m n ng h m 1 n t 31300783 2 gu n nh o 31301017 3 gu n nh uốc 31203005 4 gu n ao 1511675 5 gu n n h nh 1513055 6 gu n u n r c[.]

t i ng dụng EXCEL và MATLAB kiểm tra toán tính hệ số an

to n mái dốc

GVHD : h m n ng

h m : 1 n t 31300783

2 gu n nh o 31301017

3 gu n nh uốc 31203005

4 gu n ao 1511675

5 gu n n h nh 1513055

6 gu n u n r c 1513804

hí inh – 2017

MỤC LỤC

1 Đặt vấn đề………3

2 Tính toán bằng Excel……… 5

3 Tính toán bằng Mablat……….8

4 Kết luận………9

Báo cáo n o n n n on ậ

2

1 ĐẶT VẤN ĐỀ

1.1 Yêu cầu

- Hiện nay các ph n mềm hỗ trợ cho việc tính toán và mô phỏng ngày càng phổ biến giúp cho các k sư có ể giải quyết công việc một cách nhanh

c ón N ưn sự chính xác của các ph n mềm hỗ trợ này vẫn c ưa được

đ n i một cách chi tiết và bài bản

- Nay chúng tôi tiến hành kiểm tra kết quả tính toán của một ph n mềm c thể và so sánh chúng với kết quả từ Excel, Mablat

1.2 Bài toán

Hình 1: Mái dốc tự nhiên

- Cho mái dốc và các thông số cơ lý n ư ìn vẽ

- Mô phỏng và tính toán Hệ số an toàn mái dốc bằng ph n mềm Geostudio

Hình 2: Kết quả tính toán trên Geostudio

- Công th c hệ số an toàn

W *sin

resi ce mobilized

S





- Với bài toán trên kết quả hệ số an toàn

0.688

Báo cáo n o n n n on ậ

4

2 TÍNH TOÀN BẰNG EXCEL

2.1 Thu thập số liệu

Bảng 1

Slice

#

Width (m) Mid-height (m)

Weight (kN) Alpha (degrees) beta(m)

30 0.48012 0.18128 1.3055 -13.76 0.49431

- Số liệ có được do ph n mềm cung cấp

2.2 Tính toán

- Với số liệu từ bẳng trên chúng tôi dùng các hàm trong Excel và các phép toán cộng, trừ, n ân, c ia để tính toán

- Các hàm sử d ng:

Radians(x): chuyển từ đơn vị độ sang radian của góc

Cos(x): tính cosin của góc x

Tan(x): tính tan của góc x

Sin(x): tính sin của góc x

Sum(x): tính tổng nhiều số hạng

ABS(x): lấy giá trị tuyệ đối của x

- Từ số liệu và các hàm sử d ng chúng tôi tính toán theo m c đ c của mình

Bảng 2

5.75550 1.30699 0.47571 -3.21294 6.23121

5.08550 4.23462 1.54127 -8.98192 6.62677

4.61620 5.83967 2.12547 -14.43772 6.74167

4.26465 10.48367 3.81574 -17.75024 8.08039

4.40190 14.39304 5.23864 -22.12189 9.64054

4.13450 16.26698 5.92070 -22.80968 10.05520

3.91665 17.90443 6.51668 -23.08059 10.43333

3.73520 19.31103 7.02864 -23.01448 10.76384

3.58145 20.48428 7.45567 -22.67836 11.03712

3.44935 21.43604 7.80208 -22.11642 11.25143

3.33460 22.16854 8.06869 -21.37205 11.40329

3.23395 22.68699 8.25739 -20.47770 11.49134

3.14500 22.99739 8.37036 -19.46321 11.51536

3.06590 23.10396 8.40915 -18.35403 11.47505

Báo cáo n o n n n on ậ

6

2.99520 23.01249 8.37586 -17.16986 11.37106 2.93165 22.72832 8.27243 -15.93169 11.20408 2.87440 22.25682 8.10082 -14.65617 10.97522 2.81260 21.60183 7.86242 -13.35928 10.67502 2.77565 20.76965 7.55953 -12.05379 10.33518 2.73300 19.76382 7.19344 -10.75367 9.92644

2.69420 18.59062 6.76643 -9.47076 9.46063

2.65890 17.25371 6.27984 -8.21593 8.93874

2.62680 15.75850 5.73563 -6.99934 8.36243

2.59755 14.11068 5.13587 -5.83099 7.73342

2.57095 12.31352 4.48175 -4.71957 7.05270

2.54680 10.37251 3.77528 -3.67412 6.32208

2.52495 8.29305 3.01842 -2.70290 5.54337

2.50520 6.06948 2.20911 -1.81066 4.71431

2.48745 3.73582 1.35973 -1.01367 3.84718

2.47155 1.26803 0.46153 -0.31052 2.93308

ổn -388.54413 266.14149

- Sa i có được các số liệ n ư Bảng 2

- Tính hệ số an toàn

 * * tan  266.14149

0.685

c N Fs

ABS

 











3 TÍNH TOÁN TRONG MABLAT

- Viết code trong Editor của Mablat

- Code:

%Chuong trinh tinh he so an toan mai doc

disp( 'nhap chi tieu co ly dat' )

phi=input( 'phi=' );

c=input( 'c=' );

unit_weight=input( 'unit_weight=' );

Width=xlsread( 'D:\HK6\TIN HOC UNG DUNG\bao cao

14-4-2016' , 'B2:B31' );

Mid_height=xlsread( 'D:\HK6\TIN HOC UNG DUNG\bao cao 14-4-2016' , 'C2:C32' );

Weight=xlsread( 'D:\HK6\TIN HOC UNG DUNG\bao cao 14-4-2016' , 'D2:D32' );

Alpha=xlsread( 'D:\HK6\TIN HOC UNG DUNG\bao cao

14-4-2016' , 'E2:E32' );

Beta=xlsread( 'D:\HK6\TIN HOC UNG DUNG\bao cao

14-4-2016' , 'F2:F32' );

a=sum(c*Beta);

b=sum(Weight.*cos(Alpha*pi/180));

c=b*tan(phi*pi/180);

d=abs(sum(Weight.*sin(Alpha*pi/180)));

disp( 'He so an toan Fs' )

Fs=(a+c)/d

- Kết quả sau khi chạy code:

nhap chi tieu co ly dat

phi=20

c=5

unit_weight=15

He so an toan Fs

Fs = 0.6850

- Kết quả khi tính toán trong Mablat

Báo cáo n o n n n on ậ

8

4 K N

Sau khi tính toán lại bằng Excel và Mabtal có các kết quả sau:

Ph n mềm

Hệ số an toàn

Đ n i : các kết quả có sự chệnh lệch không lớn

Kết luận: Chúng ta có thể in ưởng vào kết quả của ph n mềm Geostudio tính toán