Hi there
to compute gradient of solutions of elasticity equation (as tutorial s.8) i
used fe_values[velocities].symmetric_gradient as :
for (unsigned int icell = 0; cell!=endc; ++cell, ++icell){
for (unsigned int i=0; i<dofs_per_cell; ++i)
for (unsigned int j=0; j<dofs_per_cell; ++j) {
for (unsigned int q_point=0; q_point<n_q_points; ++q_point){
for (unsigned int k=0; k<dofs_per_cell; ++k){
epsilon[k] =
fe_values[velocities].symmetric_gradient (k, q_point);
}
eps(i,j) += ((epsilon[i] * epsilon[j])
+(( i == j) ?
(epsilon[i] * epsilon[j]) :
0)) * fe_values.JxW(q_point);} }
in this code some parameters are defined as :
const FEValuesExtractors::Vector velocities(1) ;
std::vector<SymmetricTensor<2,dim> > epsilon(dofs_per_cell);
after checking ( aps ) calculated in code i found all zero .
is there any problem in using symmetry.gradient here ? how can i fix problem ?
any help would be appreciated.
regards,
S.M.Mohseni
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
