Dear Thomas,
this is a very interesting topic. I am thinking about using xfem for 3d
crack problems, but I have not figured out completely how to realize
especially enriched dofs in deal.II. If I can help, let me know.
Best,
Martin
----- Original Message -----
From: <[email protected]>
To: "Xujun Zhao" <[email protected]>
Cc: <[email protected]>
Sent: Saturday, August 18, 2012 11:09 AM
Subject: Re: [deal.II] xfem in deal. II
Hi Xujun,
I am working together with Wolfgang to implement XFEM in deal.ii.
Soon there will be a kind of tutorial available from the website of deal.
You need a special contruction to manage the enriched dofs. In the
tutorial we will explain how to use them for different kind of
discontinuities.
Which kind of application do you have in mind?
Thomas
Hi there,
Is there anyone working on xfem? I browsed the online tutorial, but
unfortunately couldn't find any implementation of XFEM in deal II. I have
written XFEM, Fast marching, level set codes with Matlab, but not
familiar
with deal ii library yet. anyone can give me some hints? I read some
examples recently, but still have some points unclear
1. By justifying the signs of level set function in each element, I can
easily identify the standard element and enriched element. But how can I
change the settings of dof's in these elements. For example, in the
enriched elements, additional dof's exist due to the enrichment. Also,
the
local_dof_indices are changed in this process.
2. We can achieve loops over each active cell to calculate local
stiffness
matrix using DoFHandler<dim>:: active_cell_iterator. My question is can I
do the loop operation over each nodes and also find the corresponding
supporting elements around it? Any connection functions from nodes to
elements? for example, 3x3 Q1 element, totally 4x4 nodes, I would like to
do the following loops
for i=1:total_num_nodes
support_element_number = some_function( i );
n_support_element = size( support_element_number );
for j=1:n_support_element
... operations in current j-th element
...
end
end
3. If I assembled two sets of global matrix, say, [ A ], [ B ] and rhs
{f},
respectively
[ A ]*{u} + [ B ]*{a} = {f}
[ B ]_T *{a} = {g}
how can I combine them as a big matrix K ? for example,
K = [ A B ]
[ BT 0 ]
F = {f}
{g}
Thank you very much.
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
