[deal.II] Re: Step-42

Dear All, Thanks very much for your help. I really appreciate your kindness if you let me know why this step cannot be run as other tutorials? Bests, Ben On Monday, May 15, 2017 at 4:12:17 PM UTC-5, Bruno Turcksin wrote: > > Benhour, > > do the following three steps: > cmake . > make >

[deal.II] Re: Step-42

Dear Daniel, According to your noteworthy comments, I tried exactly what you said, however, I got exactly the same as before, as what you got. Thanks, Benhour On Monday, May 15, 2017 at 12:23:44 PM UTC-5, Daniel Arndt wrote: > > Benhour, > > `make run` gives me: > > Scanning dependencies of

[deal.II] Re: Step-42

Dear All, Unfortunately I could not run this step. In fact, none of my friends also could run this step. I was wondering if you let me know it works in your system. Thanks, Benhour On Monday, May 15, 2017 at 6:25:51 AM UTC-5, Jean-Paul Pelteret wrote: > > Dear Benhour, > > I would suggest

[deal.II] Re: Step-42

Hi Daniel, I tried with version 8.4.1 and also the new version, 8.5, however it does not run at all. The errors that I got are given as below: make[3]: *** [CMakeFiles/run] Error 1 make[2]: *** [CMakeFiles/run.dir/all] Error 2 make[1]: *** [CMakeFiles/run.dir/rule] Error 2 make: *** [run] Error 2

[deal.II] Step-42

Dear All, I really appreciate it if you let me know how I can run step 42 because whenever I ran it, I got two error without any information about that. This is exactly the same with "goal oriented mesh adaptivity in elastoplasticity proiblems" example which is written by Shahram. Looking

[deal.II] distributed body force

Dear All, It would be very kind of you if you let me know how I can imply distributed load on the inner radius of a ring instead of using nodal load. Thanks, Benhour -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see

[deal.II] Re: Neumann boundary condition for the curved domain

Dear Wolfgang, It should not be zero. In addition, the local quadrature point is for the cell. Thanks, Benhour On Thursday, December 8, 2016 at 8:16:03 PM UTC-6, benhour@gmail.com wrote: > > Dear All, > I still have problem with implying boundary condition on the curved > domain. consider

[deal.II] general question

Dear All, I hope you have had a great time. I have a general question. I define total strain and deviatoric part of the elastic strain tensor in two different ways as below: 1- Total strain tensor: const SymmetricTensor<2, dim> total_strain = 0.5*(Grad_u + transpose(Grad_u)); (Tensor<2, dim>

Re: [deal.II] Re: Apply Neumann Boundary condition

Dear Jean, Thanks very much for your help and support. Most of my problem is related to applying axisymmetry boundary condition for a curved domain. I attached my code for your consideration. According to the fact that this is my first project with dealii, Please accept my lack of knowledge in

[deal.II] Re: Apply Neumann Boundary condition

Dear Jean, Thanks very much for your help and support. I do really appreciate your time and know you are really busy. The fact is that my code ran just for the first time step and can not solve for the other time step. when I refine the global mesh twice, this problem happened, however, with

[deal.II] Re: Apply Neumann Boundary condition

Dear Jean, I know that this website is not the code reviewer. I want to apply Neumann boundary condition on the curved side of my domain. It should be note that the domain is symmetry. I have considered half of a circle. It would be very kind of you if you let me know whether it is applied

[deal.II] Apply Neumann Boundary condition

Dear All, I want to apply non-zero neumann boundary condition on the curved domain of half of the circle. The way that I define my geometry, the boundary id and apply the Neumann B.C is written as follow: *Defining boundary id for different sections of domain (half_hyper_ball):* const

[deal.II] Error with the content "dealii::numbers::is_finite(sum)"

Dear All, I am solving a thermoelastic problem in cylindrical coordinate system. I used CG solver for thermal problem and Newton-Raphson method for elasticity section. For solving linear system, I used direct solver SparseDirectUMFPACK with SSOR as a preconditioner. Unfortunately, I got some

Re: [deal.II] Re: Defining tensor in different dim and spacedim

Dear All, I want to solve my thermoelasticity problem in cylindrical coordinate systems. my domain includes two quarter of a circle in 2D and two quarter of two spheres in 3D. They are bulk. I attached my domain that was created by quarter_hyper_shell for 2D. But, I really do not kno how I can

Re: [deal.II] Re: Defining tensor in different dim and spacedim

Dear All, I still have some problem in solving my equations on a curved domain (sphere in 3D or circle in 2D). I got error with this content: invalid initialization of reference of type ‘dealii::Triangulation<2, 3>&’ from expression of type ‘dealii::Triangulation<3, 3>’ Does it mean that I

Re: [deal.II] Re: Defining tensor in different dim and spacedim

Yes exactly. The only difference is in weak form of my equation which should be multiplied by 2*pi*r which I define r_value = fe_values_eta.quadrature_point(q_point)[0]; In fact, I should solve my problem on a sphere (circle). According to your noteworthy comments, the only thing that I have to

Re: [deal.II] Re: Defining tensor in different dim and spacedim

Yes You are right. I mean curved domain. My question is should i use codim instead of mapping because of the fact that there is only two tutorials in dealii which solved a differential equation on c curved domain according to the fact that my derived weak expressions are in cylindrical

[deal.II] Re: Defining tensor in different dim and spacedim

Dear Jean-Paul, Let me put it in a such way. I want to solve my equation on a curved edge. According to dealii library, I should use MappingQ whenever I have nonlinear domain. Am I correct? In addition, Should I replace all dim in defining my tensors to spacedim? It should be noted that dim =

[deal.II] Re: Defining tensor in different dim and spacedim

Dear All, I want to solve a thermoelastic problem in cylindrical coordinate system. My domain is a circle in 2D or Sphere in 3D. I am following step 44 with considering the fact that I am solving for small strain. Because I have curved domain and I should Impose a boundary condition that change

[deal.II] Defining tensor in different dim and spacedim

Dear All, I am solving an elastic problem coupled with heat equation transfer (thermal-elasticity). I should solve the problem in a curved domain (circle in 2d and sphere in 3d) . I am developing step 44 in this case. I have to define some symmetric tensors. I have attached some parts of my

[deal.II] Tensor with variable components

Dear All, I have a question. I want to define a tensor with components like the below expression, where eta is a scalar function in cylindrical coordinate system. I defined this component in dealii like this: FE_Qfe const QGauss quadrature_formula(3)

[deal.II] Derivative of a scalar function in the code

Dear All, I want to extract the first derivative and the gradient of a scalar function in dealii. Maybe I should use the void FEValuesViews::Scalar < dim, spacedim >::get_function_gradients. However, I really do

[deal.II] Re: geometry

Daniel, In that code, I had a half circle in the right part of the coordinate axis. By writing v(1) <=0, I wanted to choose the part with negative y components from the geometry, however, The attached file has been obtained. I really do not know how I should modify the source you sent to me

[deal.II] Geometry and boundary conditions

Dear All, I want to define a new geometry (a quarter of a circle) and apply symmetry boundary condition on two sides of the quarter. It would be very kind of you if you help me for that. Best, Benhour -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options,

[deal.II] Re: Elasticity

Dear all, Thanks very much for your response and help. If I can make it simple, in my equation, the lame constants in stress-strain tensor are related to another function that changes with solution vector. For instance, consider a situation in which you have a constant rank-2 tensor for initial

[deal.II] Elasticity

Dear all, In step 18, the elasticity problem has been solved with this assumption that there was not initial strain in the equation. I want to add the initial strain with consideration of the fact that this tensor is a function of the solution. In other words, this tensor should be multiplied