Download Python Interface

Python Interface, Release 4.0
•B = MdBrick(’qu term’, MdBrick pb[, int bnum[, int nfem]]) Update the tangent matrix with a int (Qu).v term.
The Q(x) parameter is a matrix field of size qdim x qdim. An example of use is for the “iku” part of Robin
boundary conditions partial_n u + iku = ...
•B = MdBrick(’mass matrix’, MeshIm mim, MeshFem mf_u[,
’real’|’complex’]) Build a mass-matrix brick.
•B = MdBrick(’generic elliptic’, MeshIm mim, MeshFem mfu[,
’scalar’|’matrix’|’tensor’][, ’real’|’complex’]) Setup a generic elliptic problem.
a(x)*grad(U).grad(V)
The brick parameter a may be a scalar field, a matrix field, or a tensor field (default is scalar).
•B = MdBrick(’helmholtz’, MeshIm mim, MeshFem mfu[, ’real’|’complex’])
Setup a Helmholtz problem.
The brick has one parameter, ‘wave_number’.
•B = MdBrick(’isotropic linearized elasticity’, MeshIm mim, MeshFem
mfu) Setup a linear elasticity problem.
The brick has two scalar parameter, ‘lambda’ and ‘mu’ (the Lame coefficients).
•B = MdBrick(’linear incompressibility term’, MdBrick pb, MeshFem mfp[,
int nfem]) Add an incompressibily constraint (div u = 0).
•B = MdBrick(’nonlinear elasticity’, MeshIm mim, MeshFem mfu, string
law) Setup a nonlinear elasticity (large deformations) problem.
The material law can be chosen among: – ‘SaintVenant Kirchhoff’ Linearized material law.
– ‘Mooney Rivlin’ To be used with the nonlinear incompressibily term.
– ‘Ciarlet Geymonat’
•B = MdBrick(’nonlinear elasticity incompressibility term’, MdBrick pb,
MeshFem mfp[, int nfem]) Add an incompressibily constraint to a large strain elasticity problem.
•B = MdBrick(’small deformations plasticity’, MeshIm mim, MeshFem mfu,
scalar THRESHOLD) Setup a plasticity problem (with small deformations).
The THRESHOLD parameter is the maximum value of the Von Mises stress before ‘plastification’ of the
material.
•B = MdBrick(’dynamic’, MdBrick pb, scalar rho[, int numfem])
brick. This brick is not fully working.
Dynamic
•B = MdBrick(’bilaplacian’, MeshIm mim, MeshFem mfu[,
’Kirchhoff-Love’]) Setup a bilaplacian problem.
If the ‘Kirchhoff-Love’ option is specified, the Kirchhoff-Love plate model is used.
•B = MdBrick(’navier stokes’, MeshIm mim, MeshFem mfu, MeshFem mfp) Setup
a Navier-Stokes problem (this brick is not ready, do not use it).
•B = MdBrick(’isotropic_linearized_plate’, MeshIm mim, MeshIm mims,
MeshFem mfut, MeshFem mfu3, MeshFem mftheta, scalar eps) Setup a linear plate
model brick.
For moderately thick plates, using the Reissner-Mindlin model. eps is the plate thinkness, the MeshFem
mfut and mfu3 are used respectively for the membrane displacement and the transverse displacement of
the plate. The MeshFem mftheta is the rotation of the normal (“section rotations”).
The second integration method mims can be chosen equal to mim, or different if you want to perform
sub-integration on the transverse shear term (mitc4 projection).
This brick has two parameters “lambda” and “mu” (the Lame coefficients)
•B = MdBrick(’mixed_isotropic_linearized_plate’, MeshIm mim, MeshFem
mfut, MeshFem mfu3, MeshFem mftheta, scalar eps) Setup a mixed linear plate model
brick.
For thin plates, using Kirchhoff-Love model. For a non-mixed version, use the bilaplacian brick.
7.8. MdBrick
35