Computational-Micromechanics.md

@@ -32,7 +32,19 @@ from dG3Dpy import*
 ```
 ### define Material law and Domain
 The basic principle of the FE implementation is to divide the whole domain of interest by domains. Each domain possesses its own behavior govering by a material law. Mechanical, thermal, or multi-physical coupled domain-material law can be used.
-
+```python
+lawnum1 = 11 # unique number of law
+E = 210E3
+nu = 0.3
+K = E/3./(1.-2.*nu)	# Bulk mudulus
+mu =E/2./(1.+nu)	  # Shear mudulus
+rho = 2.7 # Bulk mass
+sy0 = 507.
+h = E/100.
+# creation of material law
+harden = LinearExponentialJ2IsotropicHardening(lawnum1, sy0, h, 0., 10.)
+law1   = J2LinearDG3DMaterialLaw(lawnum1,rho,E,nu,harden)
+```
 ### define Domain
 ### define solver
 ### define solver options