máy bơm nước hỏa tiễn công ty in ấn https://duan-knparadise.com/ Bơm nước hỏa tiễn bóp da nam Túi nilon PP cong ty in bieu mau
Ma trận khối lượng, độ cứng, damping - Diễn đàn của các kỹ sư kết cấu Việt Nam










































































Features
Go Back   Diễn đàn của các kỹ sư kết cấu Việt Nam > THIẾT KẾ KẾT CẤU CÔNG TRÌNH BTCT VÀ GẠCH ĐÁ > Kết cấu khung, khung-vách, vách cứng, lõi bê tông cốt thép, kết cấu gạch đá
Tên thành viên
Mật mã
Tin HOT Thông tin BQT diễn đàn

Result  Reload Thống Kê - Diễn đàn của các kỹ sư kết cấu Việt Nam
All Forums | Thông báo | Văn Bản Pháp Luật | Văn Hóa Văn Nghệ Giao Lưu
Please wait...
Please wait...
Loading...
Trả lời
 
Ðiều Chỉnh Xếp Bài
Old 30-05-2017   #1
dslbkxd
Thành viên
 
dslbkxd's Avatar
 
Tham gia ngày: Apr 2010
Nơi Cư Ngụ: Yokohama
Bài gởi: 50
Thanks: 60
Thanked 603 Times in 34 Posts
Default Ma trận khối lượng, độ cứng, damping

Chào mọi người,

Trước khi có các chuơng trình tính toán hiện đại như ngày nay. Các nhà cao tầng đã có trước. Và các kĩ sư đã tính toán động học bằng tay. Tôi đang cố gắng đi tìm phuơng pháp của các kĩ sư ngày trước.

Hiện tại tôi đang cố gắng tính toán Natural Frequency của một Shear Building (X,Y, RZ) 60 tầng với kích thước kết cấu cho trước. Việc lấy natural frequency từ các chuơng trình tính toán là không quá khó. Nhưng tôi muốn xây dựng nên ma trận độ cứng, khối lượng và damping, dùng Mathlab để triển khai giải cho MDOF (mỗi tầng 3 DOF: X, Y, RZ = 180 DOF). Tôi có 2 thắc mắc như sau:

1. Có những phuơng pháp nào để xây dựng nên các ma trận tổng thể trong không gian 3D (X,Y,RZ) (độ cứng, khối lượng, damping mỗi tầng đã biết)? Tôi đã thử qua phuơng pháp Stability analysis of building based on individual bracing units, nhưng chỉ lấy được kết quả tại mode dao động đầu tiên vì không thông qua giải các ma trận tổng thể.

2. Tôi đã thử giải trong từng phương riêng biệt (60 DOF cho X) và không đi đến kết quả chính xác. Vậy điều kiện nào có thể decouple các ma trận để đưa về giải hệ 2D theo mỗi phương?

Cám ơn rất nhiều,
__________________
Doan Sy Long

Think out of the box
dslbkxd vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
Old 01-06-2017   #2
umy
Administrator
 
umy's Avatar
 
Tham gia ngày: Feb 2009
Bài gởi: 6,246
Thanks: 7,059
Thanked 14,283 Times in 4,568 Posts
Default Ðề: Ma trận khối lượng, độ cứng, damping

Thử tìm xem thêm

1) Vibration Simulation Using MATLAB and ANSYS 656Page

Michael R. Hatch - 2000 - ‎Technology & Engineering
Rayleigh form of damping matrix, 198–199 Reduction of model size ... system definition and equation, 15 frequency response derivation, 17–21 MATLAB code, ... matrix, 353–354 stiffness matrix, 11–12 stiffness matrix, finite elements global matrix ... eigenvalues/vectors, 113–115 equations of motion, 108–109 input matrix ...


2) a308123.pdf 98Page
MATLAB', Non-Uniform Beam Dynamics, Frequency Response, 91. Eigen Analysis ...
2.2.4 ELEMENT AL MASS AND STIFFNESS MATRIX FORMULATION . ...
3.1.1 EIGENVECTORS OF REPEATED EIGENVALUES . .....
The Euler-Bernoulli beam element used does not model shear effects, nor the cross-sectional rotational.

http://www.dtic.mil/dtic/tr/fulltext/u2/a308123.pdf

3) Structural Dynamics and Vibration in Practice Douglas Thorby (419 trang)

http://www.mediafire.com/view/73j450...n_Practice.pdf

thay đổi nội dung bởi: umy, 01-06-2017 lúc 09:24 PM
umy vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
The Following 2 Users Say Thank You to umy For This Useful Post:
detuphuho151 (25-04-2018), dslbkxd (02-06-2017)
Old 02-06-2017   #3
umy
Administrator
 
umy's Avatar
 
Tham gia ngày: Feb 2009
Bài gởi: 6,246
Thanks: 7,059
Thanked 14,283 Times in 4,568 Posts
Default Ðề: Ma trận khối lượng, độ cứng, damping

Trích:
Nguyên văn bởi dslbkxd View Post
Chào mọi người,

Trước khi có các chuơng trình tính toán hiện đại như ngày nay. Các nhà cao tầng đã có trước. Và các kĩ sư đã tính toán động học bằng tay. Tôi đang cố gắng đi tìm phuơng pháp của các kĩ sư ngày trước.

Hiện tại tôi đang cố gắng tính toán Natural Frequency của một Shear Building (X,Y, RZ) 60 tầng với kích thước kết cấu cho trước. Việc lấy natural frequency từ các chuơng trình tính toán là không quá khó. Nhưng tôi muốn xây dựng nên ma trận độ cứng, khối lượng và damping, dùng Mathlab để triển khai giải cho MDOF (mỗi tầng 3 DOF: X, Y, RZ = 180 DOF). Tôi có 2 thắc mắc như sau:

1. Có những phuơng pháp nào để xây dựng nên các ma trận tổng thể trong không gian 3D (X,Y,RZ) (độ cứng, khối lượng, damping mỗi tầng đã biết)? Tôi đã thử qua phuơng pháp Stability analysis of building based on individual bracing units, nhưng chỉ lấy được kết quả tại mode dao động đầu tiên vì không thông qua giải các ma trận tổng thể.

2. Tôi đã thử giải trong từng phương riêng biệt (60 DOF cho X) và không đi đến kết quả chính xác. Vậy điều kiện nào có thể decouple các ma trận để đưa về giải hệ 2D theo mỗi phương?

Cám ơn rất nhiều,
0- Giãi tay hệ thống phương trình ngày xưa khoãng 1974 (chưa có máy tính cá nhân), dùng Fourier Reihe hoặc Ritz- Ansazt (tiếng Đức) hệ thống tổng số của Sin, Cos hoặc Sinh Cosh với nhiều lời giãi !
Tìm xem trong mạng, các ĐH ngoại vẫn còn bài giãng chỉ day cho các nghiên cứu sinh về động học. Tự học thêm !

1- Lâu quá trên 40 năm rồi, tôi không còn nhớ chi tiết để hướng dẩn. Trong nhà KC đã có nhiều thành viên tuổi trung lưu khoãng 40 đến 50. Đã xong TS ở ngoài nước về dạy ĐH bên nhà hoặc công tác tiếp. (Phuho, Hiennghiem, Kienjapan, Structural ...) Có đủ trình để giúp cho cậu Long.

2- theo ngu ý, khi tách ra 2D thì ko chỉ có kết quả eigenvalues đầu tạm đúng ... còn các phần sau lệch lạc nhiều lắm.

Trong các lý thuyết của mềm, SAP, ANSYS, ABAQUS, NX ... trong Help đều có hướng dẩn (xem modal analysis) , có nhiều phương pháp đặc biệt rất hay để giãi các ma trận thật lớn

3- Cậu dslbkxd đã cày MatLab và lập trình nhuyễn tay rồi chứ ?, đưa vài bài đã làm, cho các cao thủ còn sống sót trong Dđ lấy hứng thú vào trao đổi bổ ích trở lại.

thay đổi nội dung bởi: umy, 03-06-2017 lúc 06:38 AM
umy vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
The Following User Says Thank You to umy For This Useful Post:
dslbkxd (02-06-2017)
Old 02-06-2017   #4
dslbkxd
Thành viên
 
dslbkxd's Avatar
 
Tham gia ngày: Apr 2010
Nơi Cư Ngụ: Yokohama
Bài gởi: 50
Thanks: 60
Thanked 603 Times in 34 Posts
Default Ðề: Ma trận khối lượng, độ cứng, damping

Gởi Bác Umy,

Trích:
Nguyên văn bởi umy View Post
1- Lâu quá trên 40 năm rồi, tôi không còn nhớ chi tiết để hướng dẩn. Trong nhà KC đã có nhiều thành viên tuổi trung lưu khoãng 40 đến 50. Đã xong TS ở ngoài nước về dạy ĐH bên nhà hoặc công tác tiếp. (Phuho, Hiennghiem, Kienjapan, Structure ...) Có đủ trình để giúp cho cậu Long.
.
Cảm ơn bác nhiều về các câu trả lời và giải thích. Khi cháu post bài, người cháu mong chờ có câu trả lời là bác Umy. Các đàn anh mà bác giới thiệu, cháu cũng muốn thọ giáo nhiều. Mọi người ai cũng bận rộn với công việc nên cháu không muốn gởi tin nhắn riêng mà post bài lên hi vọng mọi người có dịp sẽ trao đổi. Riêng anh Kienjapan là đàn anh xa xôi, hiện cháu đang học với thầy của ảnh tại Nhật.

Cháu cũng đang đọc về "Stick model" của thầy Nagarajaiah Satish (http://rice.academia.edu/SatishNagarajaiah). Hy vọng rằng sẽ áp dụng được.

1. Các phần mềm cháu cũng chưa thử hết. Cháu mới chạy thử bằng ETABS thì phần Story stiffness không hiện lên, cháu phải thông qua việc tính tay (lực đơn vị / chuyển vị tại mỗi tầng).

2. Việc decoupled đưa về 2D cháu nghĩ chỉ có hệ khung dầm cột có thể làm được, nhưng nếu có sự tham gia của shear wall thì việc giải trong 2D có vẻ không còn hợp lí.

Trích:
Nguyên văn bởi umy View Post
3- Cậu dslbkxd đã cày MatLab và lập trình nhuyễn tay rồi chứ ?, đưa vài bài đã làm, cho các cao thủ còn sống sót trong Dđ lấy hứng thú vào trao đổi bổ ích trở lại.
3. Mathlab cháu cũng chỉ mới sử dụng để giải hệ MDOF proportional damping, MDOF non proportional damping nhưng là thông qua code của thầy. Riêng phần cháu tự giải và nghiên cứu chỉ có phần xây dựng response spectum theo từng phuơng của trận động đất lớn Kumamoto tại Nhật bằng phuơng pháp Newmark Beta thông qua cả Excel và Mathlab. Vì vậy cháu có thể chia sẻ lên diễn đàn cùng mọi người.
https://www.dropbox.com/s/jmew019jfx...quake.rar?dl=0
__________________
Doan Sy Long

Think out of the box
dslbkxd vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
The Following 2 Users Say Thank You to dslbkxd For This Useful Post:
detuphuho151 (20-03-2018), umy (03-06-2017)
Old 03-06-2017   #5
umy
Administrator
 
umy's Avatar
 
Tham gia ngày: Feb 2009
Bài gởi: 6,246
Thanks: 7,059
Thanked 14,283 Times in 4,568 Posts
Default Ðề: Ma trận khối lượng, độ cứng, damping

1) Một lối giãi cổ điển mà Etabs và các cố nhân thường dùng cho nhà cao tầng, tháp nước, tháp TV BTCT ... tất cả đều đưa về dạng dầm consol đứng với vài mass-points. Giãi tay dể dàng, Kết quả chấp nhận được mode đầu tiên khá đúng, các mode sau lêch sai ...

Thí dụ theo bài giãng căn bản của University of Maryland
(cho mọi người dể theo dỏi, tìm hiểu)

Modal Analysis of a 4 Story Shear Building

https://www.isr.umd.edu/~austin/alad...-building.html

(lập trình ADLADIN tợ như Matlab, C++ ... có thể chuyển qua Matlab được)

- Problem Description
Lập ma trận với độ cưng côt cho mỗi tầng, Mass của sàn vào trọng điểm.

- Input File
Lựa chọn phương pháp Subsspace, ở các mềm khác như ansys có nhiều phương pháp khác nhau để lựa chọn (hiểu được, phải có căn bản! biết tại sao ?)

- Output File

- Analysis Results
Xem kết quả ở Figure 2 : Modal Shapes and Natural Periods of Vibration: chỉ có 4 Modes cho 4 mass-point

Câu hỏi Kết quả sẻ ra sao ? nếu chia mạng lưới mịn, các điễm mass cho từng phân tử cột dầm ??

2) Ma trận cho khung (với cột dầm) tương tác với vách (shearwall) khá phức tạp -dù cho 2D-
Xem thêm một TL đã được phổ biến:

*** 5_vol1_174.pdf 10Page FEM for the Interaction of frames and shear walls
http://www.iitk.ac.in/nicee/wcee/article/5_vol1_174.pdf

còn tiếp để thảo luận ....

Những TL có thể tìm thấy nhiều trong mạng, khi trích dẩn nhớ ghi nguồn. Tôn trọng bản quyền tinh thần của người tìm thấy, phổ biến.
umy vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
The Following User Says Thank You to umy For This Useful Post:
dslbkxd (03-06-2017)
Old 03-06-2017   #6
umy
Administrator
 
umy's Avatar
 
Tham gia ngày: Feb 2009
Bài gởi: 6,246
Thanks: 7,059
Thanked 14,283 Times in 4,568 Posts
Default Ðề: Ma trận khối lượng, độ cứng, damping

Vài TL lớn để tra cứu thêm cho nghiên cứu làm LV:

1) *** basic_dynamics.pdf 404Page
Fundamentals of Dynamic Analysis _Siemens Nastran NX
https://docs.plm.automation.siemens....c_dynamics.pdf

2) DIANA FEA
Tangent and Secant Stiffness Matrix
https://dianafea.com/manuals/d94/MatLib/node264.html

3) E-Wilson3D_Static_And_Dynamic_Analysis_Of_Structures 423Page
Căn bản của SAP và ETABS
Đã có cho lên Dđ, nhưng ko biết links ở đâu ? Tự tìm xem thêm.


4) Vibration_Simulation_Using_MATLAB_and_ANSYS_-_Michael_R_Hatch 656Page
Đã có cho lên Dđ, nhưng ko biết links ở đâu ? Tự tìm xem thêm.

Cậu dslbkxd cho mail lên đây, tôi gởi TL 3) và 4) ! Up lên lại cho Dđ !

còn tiếp ...

thay đổi nội dung bởi: umy, 03-06-2017 lúc 02:25 AM Lý do: Thêm
umy vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
The Following 3 Users Say Thank You to umy For This Useful Post:
detuphuho151 (25-04-2018), dslbkxd (03-06-2017), the546 (03-06-2017)
Old 03-06-2017   #7
the546
Thành viên
 
Tham gia ngày: Dec 2011
Bài gởi: 58
Thanks: 33
Thanked 10 Times in 8 Posts
Default Ðề: Ma trận khối lượng, độ cứng, damping

Em đang có 1 số tài liệu về môn động lực học mà thầy giáo em cho,nếu a cần dùng để tham khảo thì để lại email em gửi cho,em không up lên dđ được
the546 vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
The Following User Says Thank You to the546 For This Useful Post:
dslbkxd (03-06-2017)
Old 03-06-2017   #8
atomdie
Thành viên rất tích cực
 
Tham gia ngày: Sep 2010
Bài gởi: 347
Thanks: 35
Thanked 171 Times in 104 Posts
Default Ðề: Ma trận khối lượng, độ cứng, damping

Em hiểu câu hỏi của bác như này đúng không ạ. Để lấy được natural frequency, tức giải phương trình như sau:

[M]u''+[C]u'+[K]u=0 trong đó u là chuyển vị, [M] là ma trận khối lượng, [C] là ma trận damping, [K] là ma trận độ cứng.

- Với ma trận khối lượng [M] tính thông qua density của từng phần tử, có hai cách xây dựng. Cách 1, direct mass lumping. Ví dụ dầm dài L, diện tích A, density v gồm 2 nút thì mỗi nút chịu 1/2 trọng lượng. Một cách khác là Variational Mass Lumping, tính kinetic energy của từng phần tử thông qua shape function và gán lại vào ma trận global mass.
atomdie vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
The Following 2 Users Say Thank You to atomdie For This Useful Post:
dslbkxd (03-06-2017), umy (03-06-2017)
Old 03-06-2017   #9
dslbkxd
Thành viên
 
dslbkxd's Avatar
 
Tham gia ngày: Apr 2010
Nơi Cư Ngụ: Yokohama
Bài gởi: 50
Thanks: 60
Thanked 603 Times in 34 Posts
Default Ðề: Ma trận khối lượng, độ cứng, damping

Gởi bác Umy,

1. Cảm ơn bác nhiều vì các tài liệu quí giá.

2. Bài giảng 4 tầng 2D cháu có từng xem qua. Với bài toán Undamp và độ cứng cho trước thì rất tốt cho việc thực hành các chương trình.

Ở đây cháu dùng tạm Matlab giải với non proportional damping. Chú ý ma trận trong matlab được viết từ tầng 1 đến tầng 4 và đơn vị của ma trận độ cứng phải là N/m

Mode 1:

Mode 2:


Mode 3:


Mode 4:


Kết quả ra tương tự với bài giảng của thầy Mark Austin. Mode 1 là 10.85 Rad/s tương đương 0.58s

3. Với việc chia nhỏ mass tại mỗi tầng thành nhiều lump mass của dầm và nối với nhau bởi độ cứng sàn. Bài toán sẽ phức tạp hơn. Ngay cả ma trận độ cứng cũng sẽ thay đổi. Các mode dao động cũng nhiều hơn và mode shape cũng đa dạng hơn.

Ở những mode đầu tiên khi mà các lump mass mode shape cùng chiều nhau thì bài toán có dạng mode shape giống 4 lump mass. Nhưng tiến đến mode lớn hơn (phải là rất lớn vì độ cứng ngang của sàn khá lớn), các mass mode shape của các dầm di chuyển ngược chiều nhau. Thì 2 bài toán đã khác nhau.

Đó là những gì cháu có thể hình dung ra được.



Trích:
Nguyên văn bởi umy View Post

3) E-Wilson3D_Static_And_Dynamic_Analysis_Of_Structures 423Page
Căn bản của SAP và ETABS
Đã có cho lên Dđ, nhưng ko biết links ở đâu ? Tự tìm xem thêm.


4) Vibration_Simulation_Using_MATLAB_and_ANSYS_-_Michael_R_Hatch 656Page
Đã có cho lên Dđ, nhưng ko biết links ở đâu ? Tự tìm xem thêm.

Cậu dslbkxd cho mail lên đây, tôi gởi TL 3) và 4) ! Up lên lại cho Dđ !

còn tiếp ...
4. Tình cờ cháu cũng có 2 cuốn sách số (3) và số (4) bác nêu. Cháu xin chia sẻ ở bên dưới

(3) Hatch, M. R. (2000). Vibration simulation using MATLAB and ANSYS: CRC Press.: https://www.dropbox.com/s/p3zn9xx9az...ANSYS.pdf?dl=0

(4) Wilson, E. L. (1996). Three-dimensional static and dynamic analysis of structures. Computers and Structures, Inc., Berkeley, CA. : https://www.dropbox.com/s/x9necvao8l...tures.pdf?dl=0
__________________
Doan Sy Long

Think out of the box

thay đổi nội dung bởi: dslbkxd, 03-06-2017 lúc 03:46 AM Lý do: Hình ảnh không xuất hiện
dslbkxd vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
The Following 2 Users Say Thank You to dslbkxd For This Useful Post:
detuphuho151 (25-04-2018), umy (03-06-2017)
Old 03-06-2017   #10
dslbkxd
Thành viên
 
dslbkxd's Avatar
 
Tham gia ngày: Apr 2010
Nơi Cư Ngụ: Yokohama
Bài gởi: 50
Thanks: 60
Thanked 603 Times in 34 Posts
Default Ðề: Ma trận khối lượng, độ cứng, damping

Trích:
Nguyên văn bởi the546 View Post
Em đang có 1 số tài liệu về môn động lực học mà thầy giáo em cho,nếu a cần dùng để tham khảo thì để lại email em gửi cho,em không up lên dđ được
Cảm ơn The564, bạn có thể gởi cho tôi theo địa chỉ mail: dslbkxd@gmail.com được ko?
__________________
Doan Sy Long

Think out of the box
dslbkxd vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
Old 03-06-2017   #11
umy
Administrator
 
umy's Avatar
 
Tham gia ngày: Feb 2009
Bài gởi: 6,246
Thanks: 7,059
Thanked 14,283 Times in 4,568 Posts
Default Ðề: Ma trận khối lượng, độ cứng, damping

Trích:
Nguyên văn bởi atomdie View Post
Em hiểu câu hỏi của bác như này đúng không ạ. Để lấy được natural frequency, tức giải phương trình như sau:

[M]u''+[C]u'+[K]u=0 trong đó u là chuyển vị, [M] là ma trận khối lượng, [C] là ma trận damping, [K] là ma trận độ cứng.

- Với ma trận khối lượng [M] tính thông qua density của từng phần tử, có hai cách xây dựng. Cách 1, direct mass lumping. Ví dụ dầm dài L, diện tích A, density v gồm 2 nút thì mỗi nút chịu 1/2 trọng lượng. Một cách khác là Variational Mass Lumping, tính kinetic energy của từng phần tử thông qua shape function và gán lại vào ma trận global mass.
Đố vui để học thêm cho cậu atomdie, đã dùng nhiều Ansys modal analysis.
Tìm hiểu và giãi thích các Method sau đây cho MODOPT.

Mode extraction method to be used for the modal analysis.

LANB — Block Lanczos

LANPCG — PCG Lanczos

SUBSP — Subspace iteration

REDUC — Householder (reduced)

UNSYM — Unsymmetric matrix. This option cannot be followed by a subsequent spectrum analysis.

DAMP — Damped system. This option cannot be followed by a subsequent spectrum analysis.

QRDAMP — Damped system using QR algorithm. This option cannot be followed by a subsequent spectrum analysis.

VT — Solve with the Variational Technology method of ANSYS DesignXplorer.

thay đổi nội dung bởi: umy, 03-06-2017 lúc 04:54 AM
umy vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
The Following User Says Thank You to umy For This Useful Post:
dslbkxd (03-06-2017)
Old 03-06-2017   #12
umy
Administrator
 
umy's Avatar
 
Tham gia ngày: Feb 2009
Bài gởi: 6,246
Thanks: 7,059
Thanked 14,283 Times in 4,568 Posts
Default Ðề: Ma trận khối lượng, độ cứng, damping

Có gởi theo Mail đến cậu dslksxd, thưởng thêm 4 Lecture nhỏ, Up lên cho ACE xem thêm.

1-Structural_analysis_for performance based earthquake engineering 85 Trang

- Xem kỹ trang 36 -45, OpenSees (PEER) Rất tiếc! Tôi không rành !!
- Trang 46 các mềm thương mãi tính động.


2-E.Wilson dynamic_analysis_basics 96 trang
Hướng dẩn của Thầy viết SAP, ETABS về tính động cho Tower.

3-Dynamic Reponse of the Transmission Tower-Line
Bài tính động (động đất , gió động) cho tháp thép trong thực tế

4-Modal and Harmonic Analysis with ANSYS 58Trang
TL hướng dẩn chính thức của Ansys về modal và harmonic analysis >> Harmonic analysis thường được dùng để phân tích cộng hưởng (Resonance)

Căn dặn của Admin umy:
4 Lecture nầy các cậu làm nghiên cứu nên giử lấy học thêm

thay đổi nội dung bởi: umy, 03-06-2017 lúc 04:57 AM
umy vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
The Following 2 Users Say Thank You to umy For This Useful Post:
detuphuho151 (25-04-2018), dslbkxd (03-06-2017)
Old 03-06-2017   #13
atomdie
Thành viên rất tích cực
 
Tham gia ngày: Sep 2010
Bài gởi: 347
Thanks: 35
Thanked 171 Times in 104 Posts
Default Ðề: Ma trận khối lượng, độ cứng, damping

Trích:
Nguyên văn bởi umy View Post
Đố vui để học thêm cho cậu atomdie, đã dùng nhiều Ansys modal analysis.
Tìm hiểu và giãi thích các Method sau đây cho MODOPT.

Mode extraction method to be used for the modal analysis.

LANB — Block Lanczos

LANPCG — PCG Lanczos

SUBSP — Subspace iteration

REDUC — Householder (reduced)

UNSYM — Unsymmetric matrix. This option cannot be followed by a subsequent spectrum analysis.

DAMP — Damped system. This option cannot be followed by a subsequent spectrum analysis.

QRDAMP — Damped system using QR algorithm. This option cannot be followed by a subsequent spectrum analysis.

VT — Solve with the Variational Technology method of ANSYS DesignXplorer.
Cháu xin trả lời thế này. Công thức tổng quát để tìm Eigenvalues và Eigenvectors như sau:

L*L[M]{E}+[K]{E}+L[C]{E}=0 trong đó L là Eigenvalues, {E} là Eigenvector, [M] ma trận khối lượng, [K] ma trận độ cứng, và C là ma trận damped.

Có hai cách giải chủ yếu cho hệ phương trình tuyến tính direct solver và lặp.

1. Với phương pháp giải trực tiếp direct solver, yêu cầu bộ nhớ RAM lớn, giải 1 lần là xong, không cần lặp.

2. Giải lặp, có nhiều phương pháp như PCG, Jacobi, bicgstab... Thường đi kèm với phương pháp precondition matrix để giúp hội tụ nhanh hơn.

Tiếp tục, do đặc điểm của hệ phương trình tuyến tính, có thể chia ra:
1. Hệ đối xứng, SYM
2. Hệ không đối xứng, UNSYM

Với hệ phương trình gốc, có thể kể đến:
1. Có damped matrix [C] hay không


Kết hợp tất cả các yếu tố trên, ANSYS đưa ra lựa chọn để extract Eigenvector.

Hầu hết các phần mềm đều sử dụng UMFPACK làm core thư viện để giải hệ phương trình tuyến tính. Ansys sử dụng thư viện Intel MKL được Intel tối ưu hóa trên nền UMFPACK nhằm tối ưu tốc độ tính toán multicore của chip Intel.

Với mỗi dạng ma trận, Intel MKl có một cách tiếp cận riêng nhằm tối ưu hóa việc giải hệ phương trình, do đó ANSYS cũng phân ra nhiều option cho users lựa chọn.

Ngoài ra trong version mới, Ansys tích hợp thêm CUDA solver với thư viện cuBlas và cuSparse cho matrix-vector operation.

//Bản thân cháu đang dùng cả Intel MKL, cuSparse, cuBlas, cuSolver để lập trình trên nền C++.

Cháu có cái note về GPU computing, https://www.researchgate.net/publica...puting_for_FEM

thay đổi nội dung bởi: atomdie, 03-06-2017 lúc 05:33 AM
atomdie vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
The Following 2 Users Say Thank You to atomdie For This Useful Post:
dslbkxd (03-06-2017), umy (03-06-2017)
Old 03-06-2017   #14
dslbkxd
Thành viên
 
dslbkxd's Avatar
 
Tham gia ngày: Apr 2010
Nơi Cư Ngụ: Yokohama
Bài gởi: 50
Thanks: 60
Thanked 603 Times in 34 Posts
Default Ðề: Ma trận khối lượng, độ cứng, damping

Trích:
Nguyên văn bởi umy View Post
Có gởi theo Mail đến cậu dslksxd, thưởng thêm 4 Lecture nhỏ, Up lên cho ACE xem thêm.

1-Structural_analysis_for performance based earthquake engineering 85 Trang

- Xem kỹ trang 36 -45, OpenSees (PEER) Rất tiếc! Tôi không rành !!
- Trang 46 các mềm thương mãi tính động.


2-E.Wilson dynamic_analysis_basics 96 trang
Hướng dẩn của Thầy viết SAP, ETABS về tính động cho Tower.

3-Dynamic Reponse of the Transmission Tower-Line
Bài tính động (động đất , gió động) cho tháp thép trong thực tế

4-Modal and Harmonic Analysis with ANSYS 58Trang
TL hướng dẩn chính thức của Ansys về modal và harmonic analysis >> Harmonic analysis thường được dùng để phân tích cộng hưởng (Resonance)

Căn dặn của Admin umy:
4 Lecture nầy các cậu làm nghiên cứu nên giử lấy học thêm
Bác admin Umy có gởi đến anh em diễn đàn 4 bài giảng quan trọng (Cảm ơn Bác):

1. Structural analysis for performance based earthquake engineering - FEMA
2. Dynamic analysis - Ed Wilson
3. Dynamic response of transmission tower line - Bochen et.al.
4. Modal and harmonic analysis with ANSYS - David Herrin


Tôi cũng gởi 1 vài bài giảng hữu ích về chủ đề động lực học thông qua các chuơng trình summer school (Smart structure) giữa các trường tại Mỹ và Nhật:
1.
SDOF system - Narutoshi Nakata
2.
MDOF structural dynamic - Erik Johnson
3.
Adaptive negative stiffness structural system - Satish Snagarajaiah
4.
Strength Deformability Integrity and strong column in seismic design of multi story structures - Akira Wada
5.
Structural control for civil engineering applications - Hyung Yo Yung
6.
Real time hybrid simulation: Design and execution
__________________
Doan Sy Long

Think out of the box

thay đổi nội dung bởi: dslbkxd, 03-06-2017 lúc 06:04 AM
dslbkxd vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
The Following User Says Thank You to dslbkxd For This Useful Post:
umy (03-06-2017)
Old 07-06-2017   #15
umy
Administrator
 
umy's Avatar
 
Tham gia ngày: Feb 2009
Bài gởi: 6,246
Thanks: 7,059
Thanked 14,283 Times in 4,568 Posts
Default Ðề: Ma trận khối lượng, độ cứng, damping

Xem thêm #42 và 44, có nhắc đến cách giãi ma trận

http://www.ketcau.com/forum/showthread.php?t=606&page=3
(thời vàng son của D đ kết cấu, có nhiều cao thủ hăng say đóng góp)


Trong các ngôn ngữ lập trình, thường có sẳn Matrix operation để giãi hệ thống phương trình ! ko cần lăng tăng tự viết mất thời giờ nhiêu.
umy vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
Old 07-06-2017   #16
atomdie
Thành viên rất tích cực
 
Tham gia ngày: Sep 2010
Bài gởi: 347
Thanks: 35
Thanked 171 Times in 104 Posts
Default Ðề: Ma trận khối lượng, độ cứng, damping

Trích:
Nguyên văn bởi umy View Post
Xem thêm #42 và 44, có nhắc đến cách giãi ma trận

http://www.ketcau.com/forum/showthread.php?t=606&page=3
(thời vàng son của D đ kết cấu, có nhiều cao thủ hăng say đóng góp)


Trong các ngôn ngữ lập trình, thường có sẳn Matrix operation để giãi hệ thống phương trình ! ko cần lăng tăng tự viết mất thời giờ nhiêu.
Vừa tuần trước cháu có nói chuyện với mode tnlinh, và bác ấy đề cập tới cái bộ code của bác ấy mua hồi năm 2005, hình như đúng là cái Mtxvec (Intel MKL based) bác umy ạ.
atomdie vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
The Following User Says Thank You to atomdie For This Useful Post:
umy (07-06-2017)
Old 07-06-2017   #17
umy
Administrator
 
umy's Avatar
 
Tham gia ngày: Feb 2009
Bài gởi: 6,246
Thanks: 7,059
Thanked 14,283 Times in 4,568 Posts
Default Ðề: Ma trận khối lượng, độ cứng, damping

TU Darmstadt có nhiều nhân tài gốc Việt đấy nhé !

Gần đây có xem trên TV thấy có Anh Trần Quốc Khánh về ĐH dạy về đèn LED

http://www.lichttechnik.tu-darmstadt...q_34952.de.jsp

Cậu atomdie cho tôi gửi lời thăm anh tnlinh và gia quyến.
Cóc già umy xin cuối đầu thán phục các bạn trẻ tài ba.

thay đổi nội dung bởi: umy, 07-06-2017 lúc 09:16 PM
umy vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
The Following 2 Users Say Thank You to umy For This Useful Post:
atomdie (07-06-2017), detuphuho151 (20-03-2018)
Loading...
Trả lời

Ðiều Chỉnh
Xếp Bài

Quyền sử dụng ở Diễn Ðàn
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is Mở
Smilies đang Mở
[IMG] đang Mở
HTML đang Mở

Chuyển đến



Múi giờ GMT. Hiện tại là 07:46 AM.


DIỄN ĐÀN http://ketcau.com/forum NƠI HỘI TỤ CỦA CÁC KỸ SƯ KẾT CÂU VIỆT NAM
WWW.KETCAU.COM - CẦU NỐI CỦA CÁC KỸ SƯ KẾT CẤU CÔNG TRÌNH, ĐỊA KỸ THUẬT VIỆT NAM. DIỄN ĐÀN ĐƯỢC MUA BẢN QUYỀN CỦA JELSOFT ENTERPRISES Ltd. Thiết kế website