3.1. Thông s đ u vào và thu t toán s d ng.
Lu n v n s d ng ph ng pháp Newmark gia t c trung bình tìm nghi m d ng chuy n v. Các b c đ gi i ph ng trình chuy n đ ng trong lu n v n theo ph ng pháp Newmark đ c trình bƠy nh sau:
3.1.1. Thơng s đ u vào
Xác đnh các d li u c a bài tốn g m: các thơng s c a h kh i l ng – lò xo – c n g m kh i l ng m, đ c ng lò xo k, h s c n c và di chuy n v i v n t i V; các thông s k t c u t m Mindlin bên trên g m tr ng l ng riêng , module đƠn h i E , h s Poission c a t m, kích th c L B, , h; các thông s c a đ t n n l n l t là: kwf , ksf, c . Các thông s lf n l t đ c li t kê trong các B ng 3.1, B ng 3.2 và B ng 3.3. B ng 3.1: T a đ và tr ng s các đi m Gauss Module đƠn h i (N/m2) H s Poisson Kích th c t m (m) Tr ng l ng riêng (kg/m3) E L B h
B ng 3.2: Thông s n n hai thông s Pasternak
H s đ c ng (N/m3) H s kháng c t (N/m) H s c n n n (N.s/m3) kwf ksf cf B ng 3.3: Thông s t i tr ng Kh i l ng c ng lò xo H s c n V n t c (kg) (N/m) (N.s/m) (m/s) m k c V
30
Thi t l p các ma tr n kh i l ng M, các ma tr n đ c ng K, ma tr n c n C c a
k t c u t m và n n b ng cách ghép n i ma tr n.
Xác đnh ma tr n t i tr ng tác d ng lên t m c n kh o sát. Sau đó thi t l p ph ng trình chuy n đ ng và ch n b c th i gian t. Nh p đi u ki n ban đ u d ,0 d0 và 1 0 0 0 0 d M P Cd Kd . R i r c hóa vect t i tr ng theo bi n th i gian.
Gi i ph ng trình chuy n đ ng b ng ph ng pháp tích phơn Newmark. Gi i bài tốn theo d ng tìm chuy n v , xu t các k t qu , v bi u đ và l p b ng th ng kê đ phân tích các y u t nh h ng đ n chuy n v c a t m Mindlin trên n n Pasternak. T đó rút ra nh n xét, đánh giá vƠ k t lu n.
3.1.2. Gi i bài toán theo d ng chuy n v
Xác đnh ma tr n kh i l ng hi u d ng.
Tính vect t i tr ng hi u d ng t i th i đi m n1.
Gi i h ph ng trình đ i s tuy n tính (2.89) đ tìm gia t c t i th i đi m n1 là 1
n
d .
Tìm các giá tr v n t c và chuy n v t i th i đi m n1 theo các ph ng trình (2.86) và (2.87).
3.1.3. Gi i bài toán theo d ng gia t c
Xác đnh ma tr n đ c ng hi u d ng.
Tính vect t i tr ng hi u d ng t i th i đi mn1.
Gi i h ph ng trình đ i s tuy n tính (2.94) đ tìm chuy n v t i th i đi m n1 là dn1.
Tìm các giá tr v n t c và gia t c t i th i đi m n1theo các ph ng trình (2.92) và (2.93).
3.1.4. năđ nh và h i t c aăph ngăphápăNewmark
Nh đư đ c p trong m c 2.4, ph ng pháp Newmark v i 1 2
và 1 4
còn g i lƠ ph ng pháp gia t c trung bình cho s n đ nh không đi u ki n vƠ đ chính xác t t.
Ch ng 3: N i dung ph ng pháp nghiên c u
31
Do đó Lu n v n nƠy s d ng ph ng pháp Newmark gia t c trung bình v i 1 2
và 1
4
đ gi i bài toán. S h i t s đ c ti n hành ki m tra trong Lu n v n.
Trong Lu n v n nƠy s d ng ph ng pháp Newmark vƠ ngôn ng l p trình Matlab phiên b n (R2019a) đ tìm nghi m.
32