Hi Developers,
Questions again, I made a small deformation elastic problem with Great
libmesh,
then extended it to the Large Deformation. A couple of arrays are
introduced to
accommodate the deformation gradient, elastic module and so on, which are
defined as following,
std::vector<Real> F; // deformation gradient
//spatial elastic modulus c_ijkl 6x6 = aa[][]
std::vector<std::vector<Real> > aa;
aa.resize(6);
for (int i = 0; i < 6; ++i)
aa[i].resize(6);
for(int i=0;i<6;i++){
for (int j=0;j<6;j++){
aa[j][i] = 0.0;
}
}
Now the errors below bother me, which looks like container problem.
did I make mistake on declaring the container? My data structure
background is pretty weak. Thanks for any gives.
Bests,
Liang
++++++++++++++++ output in screen +++++++++++++++++
li...@liang-laptop:~/lib/libmesh/libmesh_liang/ldm3d_v7$ mpirun -np 2
./ldm3d_v7-dbg
Running ./ldm3d_v7-dbg
Mesh Information:
mesh_dimension()=3
spatial_dimension()=3
n_nodes()=1331
n_local_nodes()=735
n_elem()=1000
n_local_elem()=500
n_active_elem()=1000
n_subdomains()=1
n_processors()=2
processor_id()=0
EquationSystems
n_systems()=1
System "largeDeform_disp_3d"
Type "TransientLinearImplicit"
Variables="u0" "u1" "u2" "v0" "v1" "v2"
Finite Element Types="LAGRANGE" "LAGRANGE" "LAGRANGE" "LAGRANGE"
"LAGRANGE" "LAGRANGE"
Approximation Orders="FIRST" "FIRST" "FIRST" "FIRST" "FIRST" "FIRST"
n_dofs()=7986
n_local_dofs()=4410
n_constrained_dofs()=0
n_vectors()=3
Solving time step 0, time=0.00100...
/usr/include/c++/4.4/debug/vector:265:error: attempt to subscript
/usr/include/c++/4.4/debug/vector:265:error: attempt to subscript container
with out-of-bounds index 0, but container only holds 0 elements.
Objects involved in the operation:
sequence "this" @ 0x0xbfc77c88 {
type = NSt7__debug6vectorIdSaIdEEE;
}
container
with out-of-bounds index 0, but container only holds 0 elements.
Objects involved in the operation:
sequence "this" @ 0x0xbfe1ccd8 {
type = NSt7__debug6vectorIdSaIdEEE;
}
rank 1 in job 4 liang-laptop_34550 caused collective abort of all ranks
exit status of rank 1: killed by signal 9
------------------------------------------------------------------------------
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
