Hi developers and users,
I made a 3D finite deformation elasticity input code and got the executable
file, while the logical error displayed
in the screen when running it. The Newmark solver is used for time steps and
no Newton solver used because
I would make the code simpler at beginning of trail. This local errors makes
me wonder, I have no any clue to
modify the code, does it lie in the main or assemble part, or even in the
boundary condition? The assemble part
is quit straightforward based on the Zienkiewicz's, the [K] martix follows
the classic steps involving the displacement, deformation gradient, push
forward the dphi[][], neoHookean algorithm and elastic tensor c_ijkl. The
mass matrix [M] is set as lumped, which
contains only diagonal components. No damping and body force here. The whole
process looks fine but still get
this error, does any one could give me some advices to solve it? Thanks in
advance!
Liang
************************** error
**************************************************
li...@b356:~/Desktop/libmesh/libmesh_liang/dld3d_v5$ ./dld3d_v5-dbg
Running ./dld3d_v5-dbg
Mesh Information:
mesh_dimension()=3
spatial_dimension()=3
n_nodes()=2013
n_local_nodes()=2013
n_elem()=1200
n_local_elem()=1200
n_active_elem()=1200
n_subdomains()=1
n_processors()=1
processor_id()=0
EquationSystems
n_systems()=1
System "Disp_finite_deform_3d"
Type "Newmark"Error in libMesh internal logic
Variables="u" "v" "w"
Finite Element Types="LAGRANGE", "JACOBI_20_00" "LAGRANGE",
"JACOBI_20_00" "LAGRANGE", "JACOBI_20_00"
Infinite Element Mapping="CARTESIAN" "CARTESIAN" "CARTESIAN"
Approximation Orders="FIRST", "THIRD" "FIRST", "THIRD" "FIRST", "THIRD"
n_dofs()=6039
n_local_dofs()=6039
n_constrained_dofs()=0
n_vectors()=9
Assertion `!_points.empty()' failed.
[0] /home/liang/Desktop/libmesh/include/quadrature/quadrature.h, line 116,
compiled Apr 28 2010 at 12:08:50
terminate called after throwing an instance of 'libMesh::LogicError'
what(): Error in libMesh internal logic
Aborted (core dumped)
------------------------------------------------------------------------------
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
