| ... | @@ -31,8 +31,8 @@ mu =E/2./(1.+nu) # Shear mudulus |
... | @@ -31,8 +31,8 @@ mu =E/2./(1.+nu) # Shear mudulus |
|
|
rho = 2.7 # Bulk mass
|
|
rho = 2.7 # Bulk mass
|
|
|
sy0 = 507.
|
|
sy0 = 507.
|
|
|
h = E/100.
|
|
h = E/100.
|
|
|
# creation of material law
|
|
|
|
|
harden = LinearExponentialJ2IsotropicHardening(lawnum1, sy0, h, 0., 10.)
|
|
harden = LinearExponentialJ2IsotropicHardening(lawnum1, sy0, h, 0., 10.)
|
|
|
|
# creation of material law
|
|
|
law1 = J2LinearDG3DMaterialLaw(lawnum1,rho,E,nu,harden)
|
|
law1 = J2LinearDG3DMaterialLaw(lawnum1,rho,E,nu,harden)
|
|
|
```
|
|
```
|
|
|
### *Domain*
|
|
### *Domain*
|
| ... | @@ -41,7 +41,8 @@ A domain is defined as |
... | @@ -41,7 +41,8 @@ A domain is defined as |
|
|
domain = dG3DDomain(tag,phys,sp,lnum, fdg, dim, nonLocalVar, nonConstitutiveExtraDOFVar, nonCurlData)
|
|
domain = dG3DDomain(tag,phys,sp,lnum, fdg, dim, nonLocalVar, nonConstitutiveExtraDOFVar, nonCurlData)
|
|
|
```
|
|
```
|
|
|
where
|
|
where
|
|
|
- *tag*: an integer associated to this domain.
|
|
- *tag*: a positive integer associated to this domain. If only one domain is available, this number can be choose freely as a positive integer. When multiple domains are presents, this number provides different possibility:
|
|
|
|
* If *fdg=False*, and *tag* is the same for all domains, a conventional constitous Galerkin FE is found.
|
|
|
- *phys*
|
|
- *phys*
|
|
|
- *sp*
|
|
- *sp*
|
|
|
|
|
|
| ... | | ... | |
