Ðề: Implicit integration

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.

Ðề: Implicit integration

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ể.

Ðề: 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.

Ðề: Implicit integration

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

Ðề: 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)

Ðề: Implicit integration

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.
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.
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.

Ðề: 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.

Ðề: 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.

Ðề: 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??? :-)