hi all :) i've tried using limesh for a simple contact problem of two
elastic bodies, like Mengda describes, but i have not checked for the issues
that Roy mentioned above ...and the most important part for which i have not
yet come into a clear conclusion is the decision of a proper "penalty"
para
unction agnostic
> code... You never know when it might come in handy! (and libMesh
> makes it so dang easy).
>
> Derek
>
> Sent from my iPad
>
> On Jul 7, 2011, at 5:25 PM, Vasilis Vavourakis wrote:
>
> > how about the following code, would that do for you?
>
how about the following code, would that do for you?
// access one node of the mesh through this pointer
Node* node;
// set the solution vector (ie. displacement) to this node
Real soln[3] ={ u, v, w };
// loop for all system variables :
for (unsigned int ivar=0; ivarn_vars(system.number()); ivar+
to my mind, Blitz++ is indeed the best choice for implementing tensor
calculus!
2011/6/10 Boyce Griffith
> On 6/10/11 12:05 PM, John Peterson wrote:
> > On Fri, Jun 10, 2011 at 9:50 AM, Ataollah Mesgarnejad
> > wrote:
> >> Dear all,
> >>
> >> Is there a way to implement a fourth order tensor
hi all :)
it happened for me to need to evaluate the "volume" of Face elements in my
code through these functions:
Real Tri3::volume () const
Real Quad4::volume () const
...however, i noticed that the same function is not implemented for Tri6,
Quad8, Quad9 and so on... due to the fact that thes
2011/5/5 Roy Stogner
>
> On Thu, 5 May 2011, robert wrote:
>
> > is there a possibility of mesh-deformation in libmesh? I would need it
> > to model tectonic emplacement of lithospheric plates.
>
> Yes - you should be able to deform the mesh manually and still get
> valid results, just by changin
2011/1/10 Roy Stogner
>
> On Mon, 10 Jan 2011, Vasilis Vavourakis wrote:
>
> (1) in BEM mostly boundary integrations occur, which means that in 2D
>> problems only 1D discretization is needed. however, since only Edge
>> elements
>> are "loaded" to the Me
hi all and happy new year :)
i'd like to ask if anyone has implemented libmesh in Boundary Element
Methods.
i'm currently trying to built a simple 2D BEM solver for the Laplace PDE
"Grad^2(u) = 0". since i'm not that a libmesh expert i face some
difficulties in the implementation:
(1) in BEM mos
hello everybody :)
i know that this is the users area and not the developers one...but i have
built two System classes that solve steady-state and transient (Newmark)
non-linear systems through a Newton-Raphson iterative procedure.
my question is: how is it possible to attach some few source file
2010/3/23 John Peterson
> On Tue, Mar 23, 2010 at 2:41 PM, Roy Stogner
> wrote:
> >
> > On Tue, 23 Mar 2010, John Peterson wrote:
> >
> >> Roy: Is my memory correct in that this is the same "perfectly
> >> unnested" element (with respect to h-refinement) we looked at for
> >> shallow water equat
hi all :)
not sure whether this topic is already discussed in this list...
has anyone utilized LibMesh in multi-region problems (i.e. fluid-solid
interaction, different material solid-solid etc)???
i'm really really interested in that matter...
cheers,
Vasilis
---
hi Ted
i installed Libmesh+Petsc in my ubuntu x86 machine ...here are the steps i
followed in order to install all these stuff. just copy-paste each line to a
terminal window:
PETSC:
---
export PETSC_DIR=/home/Xuser/Programs/petsc-2.3.3-p15
config/./configure.py --with-cc=gcc --wi
hi Mengda
i was thinking the same thing too...i was to use libmesh so as to solve
incompressible N-S problems in 3d...
so far i have found this piece of open source code:
http://www.vmtk.org/Main/Download
go to "gnuid CFD solver" ;) and download it... it works...
however, i do want to use some o
yeah, i've been through that process some time ago and i've coded using
lbMesh for 2d elastostatics problems (as a starting case)...
i've attached a file with the code... hope you find some of the comments
helpful ;)
Vasilis
2009/5/8
> Hi all,
>
> I was wondering whether any of you has a simp
John you are right, the matrix is not DENSE, it's banded ;) i just
typed the class of matrix i declared in my assembly function!
2009/3/19 Roy Stogner :
>
> On Thu, 19 Mar 2009, Vasilis Vavourakis wrote:
>
>> i don't know if this might sound a silly question BUT i want to
i don't know if this might sound a silly question BUT i want to
perform the inversion of a "DenseMatrix" of a system...
to be more specific, i want to implement the "arc-length" procedure
for non-linear structural problems.
for this topic see: Crisfield M.A. - Non-Linear Finite Element
Analysis of
thnx for the info, i did in the past "make doc" and as u said they are not
present in class FEBase.
since that it seems that i need to use these functions i will add them in my
libmesh fe_base.h header file so as to have access to the corresponding
vectors of shape functions derivatives!
Vasilis
You are right, although the URL appointed doesn't work! Paraview calculator
filter does work for the job!
many thnx
Vasilis
2009/1/7 Norbert Stoop
> vas...@mech.upatras.gr wrote:
> > Mostly i use Paraview to visualise results through VTK or EXODUS2
> > formats...but since the velocity field v
in case the attached file is not present, i copied-pasted the c++ code
below:
** SOURCE CODE STARTS HERE
***
//=
#include
#include
Hello all :)
I attached a source file that implements a displacement-based formulation
for steady-state 2D elasticity (very simple example for someone to
understand it) and this piece of code works FINE. For my practice i also
implemented a mixed formulation (displacements-stresses) and i did mana
20 matches
Mail list logo
