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
Implicit integration - 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 > ỨNG DỤNG TIN HỌC TRONG THIẾT KẾ XÂY DỰNG > Phương pháp tính
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...
Trả lời
Ðiều Chỉnh Xếp Bài
Unread 14-06-2006   #1
kiendh
Thành viên
 
kiendh's Avatar
 
Tham gia ngày: Oct 2004
Bài gởi: 37
Thanks: 0
Thanked 8 Times in 5 Posts
Default Implicit integration

Chào các bác, hiện nay em thấy trong PTHH hay nói đến implicit integration and explicit integration để hội tụ trong cac lần lặp.
Bài toắn của em là em có cái constitutive model cho soil, đại khái là nó xét đến các ứng suất và biến dạng trên các mặt phẳng bất kỳ.sau đấy sẽ tính tổng các biến dạng trên các mặt phẳng đấy để có được global strain, ưu điểm của nó la constitutive đơn giản, chỉ xét đến us tiếp và pháp, và no có thể xét được đến soil-anisotropy (chính là cái multilaminate model mà bây h đã được cho đá) Nói chung tất cả các constitutive model đều có thể đưa về được cái framework đấy, việc lập ra constitutive model ko có j khó khăn. Nhưng vấn đề khó khăn là ở chỗ nếu cái constitutive model đấy mà nó có hardening hoặc softening thi nếu dùng substep và explicit để giải nó thì substep fai nhỏ khủng khiếp mới hội tụ được, mà tất nhiên ko fải lúc nào cũng hội tụ. Thế nên em mới định dùng implicit (backward euler, return mapping) nhưng khổ cái em tài hèn sức mọn nên làm mãi mà nó vẫn không hội tụ (cho element test chứ đừng nói đến benchmark). Có bác nào đã làm về cái đấy rồi có thể cho em xin cái code (cho implicit integration) được ko a.? Em xin cảm ơn các bác rất nhiều.
kiendh vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
Unread 14-06-2006   #2
hien nghiem
Thành viên nhiệt huyết
 
hien nghiem's Avatar
 
Tham gia ngày: Apr 2005
Nơi Cư Ngụ: Từ Liêm - Hà Nội
Bài gởi: 515
Thanks: 6
Thanked 369 Times in 123 Posts
Default Ðề: Implicit integration

Trích:
Nguyên văn bởi kiendh
Chào các bác, hiện nay em thấy trong PTHH hay nói đến implicit integration and explicit integration để hội tụ trong cac lần lặp.
Bài toắn của em là em có cái constitutive model cho soil, đại khái là nó xét đến các ứng suất và biến dạng trên các mặt phẳng bất kỳ.sau đấy sẽ tính tổng các biến dạng trên các mặt phẳng đấy để có được global strain, ưu điểm của nó la constitutive đơn giản, chỉ xét đến us tiếp và pháp, và no có thể xét được đến soil-anisotropy (chính là cái multilaminate model mà bây h đã được cho đá) Nói chung tất cả các constitutive model đều có thể đưa về được cái framework đấy, việc lập ra constitutive model ko có j khó khăn. Nhưng vấn đề khó khăn là ở chỗ nếu cái constitutive model đấy mà nó có hardening hoặc softening thi nếu dùng substep và explicit để giải nó thì substep fai nhỏ khủng khiếp mới hội tụ được, mà tất nhiên ko fải lúc nào cũng hội tụ. Thế nên em mới định dùng implicit (backward euler, return mapping) nhưng khổ cái em tài hèn sức mọn nên làm mãi mà nó vẫn không hội tụ (cho element test chứ đừng nói đến benchmark). Có bác nào đã làm về cái đấy rồi có thể cho em xin cái code (cho implicit integration) được ko a.? Em xin cảm ơn các bác rất nhiều.
Bác vào website của Plaxis download phần manual của nó, trong đó có nói đến thuật toán Backward euler và constitutive model có hardening và softening.
hien nghiem vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
Unread 14-06-2006   #3
kiendh
Thành viên
 
kiendh's Avatar
 
Tham gia ngày: Oct 2004
Bài gởi: 37
Thanks: 0
Thanked 8 Times in 5 Posts
Default Ðề: Implicit integration

Trích:
Nguyên văn bởi hien nghiem
Bác vào website của Plaxis download phần manual của nó, trong đó có nói đến thuật toán Backward euler và constitutive model có hardening và softening.
Cảm ơn bác hien nghiem nhưng em xem trong Plaxis no chỉ nói vài dòng về implicit integration ( em xem trong fan scientific) mà thực ra không fải về implicit mà nói là nó có dùng implicit integration. Về mặt lý thuyết thì em có vài chục article về các kỹ thuật implicit integration, cả flow chart nữa, nhưng về mặt thực hành thì em chịu, kiểm di kiểm lại mà vẫn không hiểu em sai hoặc hiểu nhầm ở chỗ nào. Bởi vậy em mới xin mấy bác cái code cho 1 cái constitutive model cụ thể.
kiendh vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
Unread 17-06-2006   #4
gdn
Thành viên mới
 
Tham gia ngày: Jun 2006
Bài gởi: 4
Thanks: 0
Thanked 0 Times in 0 Posts
Default Ðề: Implicit integration

Mi`nh co' code, nhung la`m nonlocal constitutive model (do co' softening) nen hoi phuc tap. Do do mi`nh chi du`ng duoc implicit o muc don gian nhat. Ba.n xem Crisfield (Non-linear Finite Element Analysis of Solids and Structures, Vol. 1, Chap 6), co noi ca'ch na`y. Dai khai la` du`ng khai trien Taylor cua yield function o 'elastic trial point' roi tu ddo' ti'nh ra plastic multiplier. Cach na`y cung duoc xem la` implicit do du`ng 'elastic trial point'; ket qua mi`nh thay tot hon han (o^?n dinh) so voi khi du`ng explicit.
gdn vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
Unread 18-06-2006   #5
kiendh
Thành viên
 
kiendh's Avatar
 
Tham gia ngày: Oct 2004
Bài gởi: 37
Thanks: 0
Thanked 8 Times in 5 Posts
Default Ðề: Implicit integration

Trích:
Nguyên văn bởi gdn
Mi`nh co' code, nhung la`m nonlocal constitutive model (do co' softening) nen hoi phuc tap. Do do mi`nh chi du`ng duoc implicit o muc don gian nhat. Ba.n xem Crisfield (Non-linear Finite Element Analysis of Solids and Structures, Vol. 1, Chap 6), co noi ca'ch na`y. Dai khai la` du`ng khai trien Taylor cua yield function o 'elastic trial point' roi tu ddo' ti'nh ra plastic multiplier. Cach na`y cung duoc xem la` implicit do du`ng 'elastic trial point'; ket qua mi`nh thay tot hon han (o^?n dinh) so voi khi du`ng explicit.
Oẹ, nonlocal constitute là cái gì thế hả bác ?
Em dùng khai triển Taylor order bậc 1 cho cái đại lượng tạm
r(ij)=sig(ij) - (sig(ij)(trial)-lamda*dg/dsig*epsilon(ij))
trong đấy cái trial tính dực trên cái predict elastic strain,cái này giữ nguyên trong các vòng lặp, sau đấy tình ra cái multiplier.
Còn cái model của em thi trong cái hardening law no có tính đến degradation (H=-Hd+Hh trong đó Hh là hardening law dựa theo Modifeed camclay và cái Hd là cái degradation law cho đất dính)nên đến 1 lúc nào đấy nếu cài Hd>Hh thì có thể xem nó là softening.
Nếu bác có cái code về cái softening này nữa thì quí hóa qua', em đang lo có khi nó ko converge la vì cái quả softening, bác có thể gửi cho em xem 1 ti được ko a.? Không hiểu được nhiều thì cũng được 1 ít.
email của em: dhk.pkl@gmail.com

thay đổi nội dung bởi: kiendh, 18-06-2006 lúc 12:44 AM
kiendh vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
Unread 18-06-2006   #6
gdn
Thành viên mới
 
Tham gia ngày: Jun 2006
Bài gởi: 4
Thanks: 0
Thanked 0 Times in 0 Posts
Default Ðề: Implicit integration

Neu ban tinh tat ca (stress, flow vector,...) tu elastic trial point va giu nguyen trong vong lap thi` mi`nh nghi chang can phai lap nhu the nua??? Ban thu chi du`ng 1 lan lap ket hop voi sub-incrementation xem sao, trong do so sub-increment duoc tinh tu elastic strain increment. Minh dung cahc na`y voi nonlocal model, chang thay co van de gi` ca?, va` cung khong den noi cham du` mi`nh phai giai ca he phuong tr`inh khi cap nhat ung suat (do ung xu tai cac diem Gauss co anh huong den nhau trong nonlocal model).

Ban khoan la`m gi den softening da (do no la`m loi giai phu thuoc luoi chia phan tu), du` theo mi`nh thi` khi cap nhat ung suat hardening hay softening cung chang khac gi` (vi` chi la` giai 1 phuong tri`nh vi phan ODE). Co`n lam cho softening thi phai du`ng nonlocal hay gradient, hay don gian nhat la` ***** band (mi`nh khong biet '***** band' cho soil thi` the na`o)

thay đổi nội dung bởi: gdn, 18-06-2006 lúc 06:03 AM
gdn vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
Unread 18-06-2006   #7
kiendh
Thành viên
 
kiendh's Avatar
 
Tham gia ngày: Oct 2004
Bài gởi: 37
Thanks: 0
Thanked 8 Times in 5 Posts
Default Ðề: Implicit integration

Trích:
Nguyên văn bởi gdn
Neu ban tinh tat ca (stress, flow vector,...) tu elastic trial point va giu nguyen trong vong lap thi` mi`nh nghi chang can phai lap nhu the nua??? Ban thu chi du`ng 1 lan lap ket hop voi sub-incrementation xem sao, trong do so sub-increment duoc tinh tu elastic strain increment.
Em chỉ giử lại cái trial stress predict (sigma(ij)(trial) thôi con multiplier (lamda)và flow vector (dg(or df)/dsignma(ij) thì fải lặp chứ anh, cón nếu giữ nguyên thì lặp cái j nữa. Initial value thì em tính theo semi-backwar euler, sau đó cứ tiép tực lặp dựa theo cái khai triển taylor bậc 1 của cái r(ij)=0.
Trích:
Nguyên văn bởi gdn
Ban khoan la`m gi den softening da (do no la`m loi giai phu thuoc luoi chia phan tu), du` theo mi`nh thi` khi cap nhat ung suat hardening hay softening cung chang khac gi` (vi` chi la` giai 1 phuong tri`nh vi phan ODE). Co`n lam cho softening thi phai du`ng nonlocal hay gradient, hay don gian nhat la` ***** band (mi`nh khong biet '***** band' cho soil thi` the na`o)
Hi`, anh nói rất đúng, softening fụ thuộc rất nhiều vào mesh nhưng em đã đọc 1 số article đề xuất 1 số những improve cho FE dể loặi bỏ cái mesh depend này và thấy rất có lý mặc dù về mặt physical nó ko co lý lắm, trong đáy cái framework của em có 1 số điểm khá vượt trội trong vấn đề xử lý numerical cho strain softening.
Trích:
Nguyên văn bởi gdn
Minh dung cahc na`y voi nonlocal model, chang thay co van de gi` ca?, va` cung khong den noi cham du` mi`nh phai giai ca he phuong tr`inh khi cap nhat ung suat (do ung xu tai cac diem Gauss co anh huong den nhau trong nonlocal model).)
Về mặt time-consuming, có lẽ em fải nói sơ qua về cái frame work của em, thực tề với 1 stress point, thay vì tính 1 lần, em fải tính nó với n lần tùy thuộc vào integration law mà em chọn (nó cũng giống nhu là cách mình tình tọa độ điểm gauss và các weight coefficient nhưng tính cho 1 cái spatial côrdinate, trong th của em là 1 cái hình cầu, ví dụ em chọn 64 điểm integrate trên cái hình cầu đấy thì em fải tính với 64 lần, sau đó cũng cộng tổng các strain sinh ra trện từng cái plane ay lại (tat nhiên là có nhân vói weight coefficient) Tức là riêng về mặt lý thuyết đã chậm hơn 64(n) lần (n là số mặt phẳng mà mình tình trên đấy). Có thể hiểu nôm na là em integrate 2 lần: lần đầu cho chính các điểm Gauss và sau đấy là cho cái global.
Cái thứ 2 mà theo ngu ý của em ảnh hưởng đến convergence là ví các plastic strain generate trên từng mf đã ko chính xác nên em integrate 1 lần nữa sẽ càng ko chính xác=> nếu dùng sub step thì step sẽ cực bé' so với các frameworh thông thường, do đó em muốn tính chính xác ngay trên từng plane.
Chính ví thế mà implicit integration có vai trò cực kỳ quan trọng trong cai framwork của em, nếu ko với 1 bài toắn khoẳng 5p (với frame work thông thường tức là chỉ integrate 1 lần) thì em fải lảm tới 10 tiếng. Lí do:
1. Integrate 2 lần (lâu hơn n lần).
2. Converge cực lâu ( do integrate 2 lần nên kết quả ko còn hội tụ nhanh nữa) (ko kiểm soat được, nếu ko converge thì tiếp tục chia nho them substep)
3. Vì ko đảm bảo converge nen fải chia ra lam nhiều sub step nhỏ, và dĩ nhiên, nhỏ hơn nhiều lần so với các substep của các framework chi integrate cac diem Gauss
Nói tóm lại là bây h em cần 1 cái ví dụ cụ thể cho implicit integration vì đấy là cách duy nhất để em củng cố cái stable cua cái code(cu chuối) của em.
Cái nữa là anh co thể giải thích cho em quả nonlocal constitutive dc ko a.? nếu có article nào thì mail cho em vao cai dc email cua em, em rất chi là cảm ơn anh.

thay đổi nội dung bởi: kiendh, 18-06-2006 lúc 02:38 PM
kiendh vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
Unread 19-06-2006   #8
gdn
Thành viên mới
 
Tham gia ngày: Jun 2006
Bài gởi: 4
Thanks: 0
Thanked 0 Times in 0 Posts
Default Ðề: Implicit integration

Mi`nh khong biet la cai framework cua ban la.i phuc tap nhu vay. Co cam tuong la` no giong giong microplane model??? Ban thu ti`m nhung ba`i bao tuong tu lien quan den microplane model (hau het cua Bazant) hay local-global constitutive law (e.g. multiple fixed *****) xem.

Ban co the email cho mi`nh cai framework ban dang du`ng duoc khong? Mi`nh se goi truoc nay mai cho may ba`i bao ve nonlocal. Dai khai la` ung xu cua 1 diem trong nonlocal model se phu thuoc va`o ung xu cua cac diem xung quanh hay cua toan ket cau. Theo mi`nh (va` ban cung da thay) thi` nonlocal hay gradient, hay may tro` gi` khac de khac phu.c nhung van de do softening gay ra deu co it y nghia vat ly, ma` chi hoac la` dam bao (ve mat toa'n ho.c) phuong tri`nh vi phan tung phan PDE khong bi. mat tinh elliptic (nonlocal hay gradient) hay it ra dam bao ve mat tieu tan nang luong (***** band hay smeared *****: dissipation la` khong doi khi thay doi luoi chia phan tu).

Ve cap nhat ung suat thi` mi`nh chi co the co chut y kien nhu vay (do framework cua ban phuc tap qua :-)) . Mi`nh khong co san vi. du cu the, nhung sao ban khong thu test cai implicit integration cua ban truoc voi constitutive law don gian thoi (local law thong thuong chang han, vi du co san trong cuon sach cua Crisfield). Ban co thu du`ng Runge-Kutta de tang do chinh xac chua?

Mi`nh khong biet bai toan 10 tieng ma` ban noi la` ba`i toan the na`o: BVP voi nonlinear law tren framework cua ban hay chi la bai toan cap nhat ung suat (voi local-global law) tai 1 material point? Neu la` truong hop thu 2 (material poiint) thi` qua that khong the chiu dung duoc :-); minh nghi la` co bug cho na`o ddo'. Mi`nh du`ng nonlocal model, moi lan cap nhat ung suat phai giai nhieu lan (20-30 lan) 1 he khoang 1000 an so (tu`y vao mesh) ma` cung chi cap nhat mat it phut.
gdn vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
Unread 19-06-2006   #9
kiendh
Thành viên
 
kiendh's Avatar
 
Tham gia ngày: Oct 2004
Bài gởi: 37
Thanks: 0
Thanked 8 Times in 5 Posts
Default Ðề: Implicit integration

Hehehe, bác nói đúng rồi a., nguồn gốc chính xác là microplane bặt nguồn từ concrete của bác bazant sau đấy được 1 số các bác trong lĩnh vực soil ứng dụng vào geomechanics. Cái này lúc chuyển từ global constitutive về local nó đơn giản lắm, có fức tạo j đâu bác, chắc tại em nói năng ko rõ nên đam ra thành fức tạp.
10 tiếng của em tất nhiên là cho toằn bộ mesh, khoẳng 2000 phần tử 6 nodes.
Cái non local của bác có vẻ thích hợp để dùng giải quyết cho các bài toắn softening và shear band bác nhỉ. Hấp dẫn qua'. Bác gửi cho em với. Tiên đây em gửi cho bác quả framework microplane mà cho rock với soil.
QUyển sách của Crisfield mãi hôm qua em mới kiếm ra và có đọc sơ qua, hôm nay đang ngồi xem lại vì từ trước đến h em chỉ xem có 2 cuôn 1 của Bathe với 1 của Smith.
P/S. Bác có quả messenger ko? Cho em cái có j em hỏi bác trực tiếp có khi vỡ ra được nhiều cái.
kiendh vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
Unread 19-06-2006   #10
gdn
Thành viên mới
 
Tham gia ngày: Jun 2006
Bài gởi: 4
Thanks: 0
Thanked 0 Times in 0 Posts
Default Ðề: Implicit integration

ma~i moi thay 1 nguoi di ga^`n duong voi mi`nh :-)
Nonlocal cung chang phai cua mi`nh :-) (la.i Bazant :-)), nhung co 1001 cach de dua nonlocal va`o mo hi`nh. Mi`nh chi du`ng no' cho mo hi`nh damage thoi, co`n shear band chac de bac nao la`m soil.
10 tieng voi mesh nhu the thi` cung khong te nhi, mi`nh co khi phai mat ca nga`y cho cai nonlocal voi may ba`i phuc tap.
Mi`nh khong du`ng messenger, nen co gi` email duoc roi. De mi`nh tim la.i may ba`i bao ve nonlocal roi mai goi cho ba.n.
Co ve ban dang o chau Au nhi, co khi o Phap la`m voi Zhao??? :-)
gdn vẫn chưa có mặt trong diễn đàn   Trả Lời Với Trích Dẫn
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à 02:40 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.