.N IăDUNGăPH NGăPHÁPăNGHIểNă CU

Một phần của tài liệu Phân tích ứng xử động của tấm nindlin chịu tác dụng của hệ khối lượng lò xo cản di chuyển sử dụng phương pháp phần tử chuyển động mem (Trang 42 - 46)

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 n1.

 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 n1 là 1

n

d  .

 Tìm các giá tr v n t c và chuy n v t i th i đi m n1 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 mn1.

 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 n1 là dn1.

 Tìm các giá tr v n t c và gia t c t i th i đi m n1theo 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

Một phần của tài liệu Phân tích ứng xử động của tấm nindlin chịu tác dụng của hệ khối lượng lò xo cản di chuyển sử dụng phương pháp phần tử chuyển động mem (Trang 42 - 46)