1 / 18

Internship presentation

Internship presentation. Numerical simulation of Multi-scale dissipative processes in geomaterials. 27 April 2010. Given by: Ghassen MOALLA ( Master student at TPS) Supervised by: Ali KARRECH Thomas POULET. Framework.

crwys
Download Presentation

Internship presentation

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Internship presentation Numerical simulation of Multi-scale dissipative processes in geomaterials 27 April 2010 Given by: Ghassen MOALLA ( Master student at TPS) Supervised by: Ali KARRECH Thomas POULET

  2. Framework • Effects of dissipative processes on the behavior of geomaterials • Behavior of poro-elasto-plastic materials : • the thermal effects • fluid flow • chemical reactions • important for the analysis of geo-materials • Escript : Fluid + T° + C° • mechanical description to be finalized

  3. Purpose of the work • Implementation of a finite element code in Pyhton • programming standing alone FEM • testing Abaqus subroutines in python • including the effects of thermal expansion and gravity • starting point to describe geomaterial problems within • a unified system

  4. Purpose of the work • Combine the code with escript • powerful library of solvers (finely ) • wide range of elements • efficient tool for solving fluid transport problem • Handle nonlinear mechanic problems

  5. Outline • Python FEM • Numerical methods • Code architecture • Results & discussion • Escript • Current achievements • Future developments

  6. Numerical methods • Newton Raphson method • Return mapping algorithm • Von mises criterion

  7. Code Architecture • Composed by key components of FEA • Element: evaluating the contribution of the element • Assembly : adding the element part to the global system • Total program : solving the system of equations

  8. Code Architecture Total Incremental Method Cond Incr dF_g Total F Increment Loading U(BC) dU_g(BC) statev_g Statev_g One Increment u_g du_g

  9. Code Architecture One Increment error>tol dF_g statev_g* dU_g(BC) One Iteration statev_g du_g* du_g Calcul error

  10. Code Architecture One Iteration statev_g* statev_g ASTIF Calcul ASTIF du_g rhs_g du_g* dF_g Solver dU_g(BC)

  11. Code Architecture Calcul ASTIF Loop nelt statev statev_g Isolate du_g du ASTIF ESM statev* Treat Element rhs statev_g* statev_g* Integrate rhs_g rhs_g

  12. Code Architecture Treat Element Loop nip 12 statev* statev* du ESMi+ ESM TreatIntegration point statev rhs rhsi+

  13. TreatIntegration point dε du Dstrain σ_ip statev Statevar_in stv_ip stv_ip* σ_ip* CalculD DDSDE statev* Statevar_out ESMi ESMCPE4 rhsi

  14. Results & Discussion Example:

  15. Escript • Mechanical description integrated • Solve the non-linearities • Von Mises • DruckerPrager • Unified framework • deformation • temperature

  16. Current achievements

  17. Future developments • Enlarge the code • plasticity criteria • return mapping algorithm • programming interface • Benchmark the results

  18. Internship presentation Numerical simulation of Multi-scale dissipative processes in geomaterials 27 April 2010 Given by: Ghassen MOALLA ( Master student at TPS) Supervised by: Ali KARRECH Thomas POULET

More Related