Dear getfem users,
I'd like to compute the divergence of stress tensor of a linear elasticity
problem which is necessary in an a posteriori error estimate.
Do you have any idea how to compute it using generic assembly of GetFEM?
Thank you
in advance for your support.
Best regards,
Phuoc
t-edges-faces
>
> Best regards,
>
> Yves
>
>
>
> Le 16/01/2018 à 06:08, Huu Phuoc BUI a écrit :
>
> Dear getfem users,
>
> I'd like to compute the divergence of stress tensor of a linear elasticity
> problem which is necessary in an a posteriori error e
Dear Kostas,
I'm interested in computation of stress of a linear elasticity problem.
I found this discussion in the getfem-users archive.
Could you please explain me what is the "mass matrix" in GetFEM++, and
how do you need the line "sigmaij =
np.transpose(gf.linsolve_mumps(mass_mat_t, sigm
Dear GetFEM++ users,
I am working on adaptive refinement of a linear elasticity problem using a
posteriori error estimate.
At each adaptive refinement step, I need to compute some quantity of
interest defined on a subdomain $\Omega$. Let's call this quantity of
interest J(u) = \int_{\Omega} div(u
,"t",True,mfef, np.zeros(mfef.nbdof()))
>
> However, if mfef is a Lagrange finite element, what you wrote will also do
> the job, may be except that it sums all the components of course.
>
> Best regards,
>
> Yves.
>
> - Original Me
Konstantinos Poulios
wrote:
> Dear Phuoc
>
> Which linear solver are you using?
>
> Best regards
> Kostas
>
> On Thu, Apr 19, 2018 at 3:30 PM, Huu Phuoc BUI > wrote:
>
>> Dear Yves,
>>
>> Thank you very much for your reply. I checked with what you p
Dear Getfem-users,
I have a problem with computing the outer faces when using a mixed mesh
(having tetras, hexa, prisms and pyramids).
The error is:
face_left = m.outer_faces_with_direction([0, 0., -1.0], 0.01)
File "/usr/local/lib/python2.7/dist-packages/getfem/getfem.py", line
1372, in o
ramid element, the transformation is singular
on the geometrical nodes.
I replaced the mean of normal vector by the normal vector on the mean
of geometrical nodes of the face. It fixes the problem. (I push the
change on the master branch).
Best regards,
Yves
Le 14/06/2018 à 15:40, Huu Phu
ingular at the pyramid
tip is a difficulty.
Best regards,
Yves
Le 15/06/2018 à 11:53, Huu Phuoc BUI a écrit :
Dear Yves,
Thank you very much for your help.
Now it worked when computing the outer faces. But I faced another
problem of this mixed mesh. In fact, I'd like to compute the
Phuoc
On 15/06/2018 16:43, Huu Phuoc BUI wrote:
Dear Yves,
That's great! It's indeed working now. Thank you very much for your help.
Best regards,
Phuoc
On 15/06/2018 16:38, Yves Renard wrote:
Dear Phuoc,
The problem is indeed the same, to compute the location of gauss
points on
a_para*Trace(Grad_u)*Id(qdim(u)) + mu_para*(Grad_u + Grad_u'));
sig_u ..."
you can also define macros with parameters, see
http://getfem.org/userdoc/gasm_high.html#macro-definition
Best regards,
Yves
Le 18/06/2018 à 13:35, Huu Phuoc BUI a écrit :
Dear Yves,
Please find enclosed
Dear GetFEM-users,
As Grad(expression) is available in the new version of GetFEM++, I
compute the divergence of the first Piola-Kirchhoff stress tensor as
(python code):
/div_sigma = "( Trace(Grad(first_Piola_Kirchhoff_stress)(1,:,:)) +
Trace(Grad(first_Piola_Kirchhoff_stress)(2,:,:)) +
Tra
respect to the given indices, so that
Contract(A, 2, 3) = A_{i,j,j} with summation on j.
Can your try again with the fix I commited to the master branch ?
Best regards,
Yves
- Original Message -----
From: "Huu Phuoc BUI"
To: "yves renard"
Cc: "franz chouly&quo
13 matches
Mail list logo