Re: [deal.II] Fwd: Regarding Step 60 tutorial of deal.ii

2019-07-25 Thread Ramprasad R
Hello all,

I would like to know if I can use one name space (step-x) in another name 
space (step-y). So that I can couple the 2 or more steps. 

Thank you.

On Wednesday, June 26, 2019 at 11:30:21 AM UTC+2, Ramprasad R wrote:
>
> Hi Daniel,
>
> Thank you so much for your detailed Explanation. As you said earlier, I 
> have to do the eigen value analysis.  You explanation has given me an 
> approach to start with. I can not thank you enough.
>
> Warm regards
>
> Ramprasad
>
> On Tue, Jun 25, 2019 at 6:27 PM Daniel Garcia-Sanchez <
> danielgarsanc...@gmail.com> wrote:
>
>>
>>
>> On Monday, June 24, 2019 at 11:36:32 PM UTC+2, Ramprasad R wrote:
>>>
>>> Hello Bangerth,
>>>
>>> The terms N_* are the resultant forces in the * direction and these 
>>> forces are calculated using the strains which in turn are calculated using 
>>> the displacements u. So the terms N_* are not constants, rather change with 
>>> each element. And these values directly depend on u.
>>>
>>>
>>>>
>> Hi Ramprasad,
>>
>> I think that you want to do an eigenvalue calculation (step-36)
>>
>> I think that first you have to do an static calculation before your 
>> eigenvalue calculation in order to obtain the values for N_* (step-8)
>>
>> As discussed in your paper, the typical eigenvalue problem for the 
>> elastic equation takes this form, where omega^2 is your eigenvalue and Phi 
>> the eigenvector
>> (K - omega^2 M)*Phi = 0
>>
>> For the buckling case lambda is your eigenvalue.
>> (K - lambda G)*Phi = 0
>>
>> The calculation of K can be found in step-8, step-62 (or other tutorials).
>>
>> I think that in order to calculate G you need the resulting strain of an 
>> static calculation. You can do the static calculation, store the strain in 
>> a temporary buffer and use that data to calculate G. step-18 shows you how 
>> to do this.
>>
>> Once you have K and G, you can do an eigenvalue calculation. step-36 
>> shows you how to do an eigenvalue calculation. Note that in step36 the 
>> stiffness matrix is called A.
>>
>> The equation in step36 is 
>> (D-epsilon M)*Phi=0
>> which is very similar to the buckling equation. 
>>
>> Best,
>> Daniel 
>>
>> -- 
>> The deal.II project is located at http://www.dealii.org/
>> For mailing list/forum options, see 
>> https://groups.google.com/d/forum/dealii?hl=en
>> --- 
>> You received this message because you are subscribed to the Google Groups 
>> "deal.II User Group" group.
>> To unsubscribe from this group and stop receiving emails from it, send an 
>> email to dealii+unsubscr...@googlegroups.com.
>> To view this discussion on the web visit 
>> https://groups.google.com/d/msgid/dealii/d026c101-ffdd-4325-a1bb-f40c8283bbc9%40googlegroups.com
>>  
>> <https://groups.google.com/d/msgid/dealii/d026c101-ffdd-4325-a1bb-f40c8283bbc9%40googlegroups.com?utm_medium=email_source=footer>
>> .
>> For more options, visit https://groups.google.com/d/optout.
>>
>

-- 
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- 
You received this message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/ab345e27-7cef-487a-a128-864b8c346ccc%40googlegroups.com.


Re: [deal.II] Re: Compatibility of Petsc with step 18

2019-07-19 Thread Ramprasad R
Hi Daniel,
The problem is now solved. The issue was that, the bash rc did not have the
location of PETSc. This line in the bash file solved the problem "export
PETSC_DIR=/path/to/petsc-x-y-z""

Regards,
Ramprasad

On Wed, Jul 17, 2019 at 5:31 PM Bruno Turcksin 
wrote:

> Ramprasad,
>
> Le mer. 17 juil. 2019 à 11:24, Ramprasad R
>  a écrit :
> > Yes I did recompile it. The details of the config files are in the
> attachment.
> There is something strange in you installation: detailed.log says that
> you are using PETSc but config.h says that you are not. It looks like
> you compiled deal.II without PETSc and then in the same directory
> recompiled deal.II with PETSc supports. This can cause the strange
> behavior that you are seeing. Remove your build directory and
> reinstall deal.II from scratch.
>
> Best,
>
> Bruno
>
> --
> The deal.II project is located at http://www.dealii.org/
> For mailing list/forum options, see
> https://groups.google.com/d/forum/dealii?hl=en
> ---
> You received this message because you are subscribed to the Google Groups
> "deal.II User Group" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to dealii+unsubscr...@googlegroups.com.
> To view this discussion on the web visit
> https://groups.google.com/d/msgid/dealii/CAGVt9eP5bqzDKi5CHFW0%3DeMF1bv%3DALRVJg2z61T2bZZEUronjA%40mail.gmail.com
> .
> For more options, visit https://groups.google.com/d/optout.
>

-- 
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- 
You received this message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/CAKS2Wq5OCjHG3m39AsTsVyz8hBC%3DNV-ddZfhz-%3DHZ0DKew0XhQ%40mail.gmail.com.


Re: [deal.II] Re: Compatibility of Petsc with step 18

2019-07-17 Thread Ramprasad R
Hello,
Yes I did recompile it. The details of the config files are in the
attachment.
Regards,
Ramprasad Raghunath

On Wed, Jul 17, 2019 at 2:14 PM Daniel Arndt  wrote:

> Ramprasad,
>
> Did you also recompile, i.e. run 'make', afterwards? Can you also show us
> the file 'detailed.log' in your build directory and the the file
> 'base/config.h' in your installation directory?
>
> Best,
> Daniel
>
> --
> The deal.II project is located at http://www.dealii.org/
> For mailing list/forum options, see
> https://groups.google.com/d/forum/dealii?hl=en
> ---
> You received this message because you are subscribed to the Google Groups
> "deal.II User Group" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to dealii+unsubscr...@googlegroups.com.
> To view this discussion on the web visit
> https://groups.google.com/d/msgid/dealii/827f9523-4b62-4557-ae4f-af916a39d7af%40googlegroups.com
> 
> .
> For more options, visit https://groups.google.com/d/optout.
>

-- 
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- 
You received this message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/CAKS2Wq56Dh6TNLz--npx7Bgxo5E2S2ugK2Z1EvkqRzr_ru83kQ%40mail.gmail.com.
For more options, visit https://groups.google.com/d/optout.
###
#
#  deal.II configuration:
#CMAKE_BUILD_TYPE:   DebugRelease
#BUILD_SHARED_LIBS:  ON
#CMAKE_INSTALL_PREFIX:   /home/ramprasad/Documents/Studienarbeit/bin
#CMAKE_SOURCE_DIR:   /home/ramprasad/Documents/Studienarbeit/dealii-9.0.1
#(version 9.0.1)
#CMAKE_BINARY_DIR:   /home/ramprasad/Documents/Studienarbeit/build
#CMAKE_CXX_COMPILER: GNU 7.4.0 on platform Linux x86_64
#/usr/bin/c++
#CMAKE_C_COMPILER:   /usr/bin/cc
#CMAKE_Fortran_COMPILER: /usr/bin/gfortran
#CMAKE_GENERATOR:Unix Makefiles
#
#  Base configuration (prior to feature configuration):
#DEAL_II_CXX_FLAGS:-pedantic -fPIC -Wall -Wextra -Wpointer-arith -Wwrite-strings -Wsynth -Wsign-compare -Wswitch -Woverloaded-virtual -Wno-placement-new -Wno-deprecated-declarations -Wno-literal-suffix -fopenmp-simd -std=c++17
#DEAL_II_CXX_FLAGS_RELEASE:-O2 -funroll-loops -funroll-all-loops -fstrict-aliasing -Wno-unused-local-typedefs
#DEAL_II_CXX_FLAGS_DEBUG:  -Og -ggdb -Wa,--compress-debug-sections
#DEAL_II_LINKER_FLAGS: -Wl,--as-needed -rdynamic
#DEAL_II_LINKER_FLAGS_RELEASE: 
#DEAL_II_LINKER_FLAGS_DEBUG:   -ggdb
#DEAL_II_DEFINITIONS:  
#DEAL_II_DEFINITIONS_RELEASE:  
#DEAL_II_DEFINITIONS_DEBUG:DEBUG
#DEAL_II_USER_DEFINITIONS: 
#DEAL_II_USER_DEFINITIONS_REL: 
#DEAL_II_USER_DEFINITIONS_DEB: DEBUG
#DEAL_II_INCLUDE_DIRS  
#DEAL_II_USER_INCLUDE_DIRS:
#DEAL_II_BUNDLED_INCLUDE_DIRS: 
#DEAL_II_LIBRARIES:m
#DEAL_II_LIBRARIES_RELEASE:
#DEAL_II_LIBRARIES_DEBUG:  
#DEAL_II_COMPILER_VECTORIZATION_LEVEL: 1
#
#  Configured Features (DEAL_II_ALLOW_BUNDLED = ON, DEAL_II_ALLOW_AUTODETECTION = ON):
#  ( DEAL_II_WITH_64BIT_INDICES = OFF )
#  ( DEAL_II_WITH_ADOLC = OFF )
#DEAL_II_WITH_ARPACK set up with external dependencies
#ARPACK_LINKER_FLAGS = 
#ARPACK_LIBRARIES = /usr/lib/x86_64-linux-gnu/libparpack.so;/usr/lib/x86_64-linux-gnu/libarpack.so;/usr/lib/x86_64-linux-gnu/liblapack.so;/usr/lib/x86_64-linux-gnu/libf77blas.so;/usr/lib/x86_64-linux-gnu/libatlas.so;gfortran;quadmath;m;c;gcc_s;gcc;/usr/lib/x86_64-linux-gnu/openmpi/lib/libmpi.so
#  ( DEAL_II_WITH_ASSIMP = OFF )
#DEAL_II_WITH_BOOST set up with external dependencies
#BOOST_VERSION = 1.65.1
#BOOST_CXX_FLAGS = -Wno-unused-local-typedefs
#BOOST_DEFINITIONS = BOOST_NO_AUTO_PTR
#BOOST_USER_DEFINITIONS = BOOST_NO_AUTO_PTR
#BOOST_INCLUDE_DIRS = /usr/include
#BOOST_USER_INCLUDE_DIRS = /usr/include
#BOOST_LIBRARIES = /usr/lib/x86_64-linux-gnu/libboost_iostreams.so;/usr/lib/x86_64-linux-gnu/libboost_serialization.so;/usr/lib/x86_64-linux-gnu/libboost_system.so;/usr/lib/x86_64-linux-gnu/libboost_thread.so;/usr/lib/x86_64-linux-gnu/libboost_regex.so;/usr/lib/x86_64-linux-gnu/libboost_chrono.so;/usr/lib/x86_64-linux-gnu/libboost_date_time.so;/usr/lib/x86_64-linux-gnu/libboost_atomic.so;pthread
#  ( DEAL_II_WITH_CUDA = OFF )
#  

Re: [deal.II] Re: Compatibility of Petsc with step 18

2019-07-17 Thread Ramprasad R
Hello Bruno,
Yes I have used the option -DDEAL_II_WITH_PETSC=ON. and I have also used
the command to specify the location of petsc.
"-DPETSC_DIR=/home/ramprasad/petsc-3.10.5". The error still pertains.

Thank you.

Regards,
Ramprasad

On Tue, Jul 16, 2019 at 9:13 PM Bruno Turcksin 
wrote:

> Ramprasad,
>
> On Tuesday, July 16, 2019 at 11:51:00 AM UTC-4, Ramprasad R wrote:
>>
>> The version of deal.ii is 9.0.1 and the version of PETSc is 3.10.5. I
>> also changed the config file of deal.ii to have MPI PETSc on.
>>
> What do you mean by "I changed the config file of deal.II"? You should not
> change any file generated by deal.II. When you compiled deal.II, did you
> use the option -DDEAL_II_WITH_PETSC=ON? If not, deal.II was not compiled
> with PETSc supports.
>
> Best,
>
> Bruno
>
> --
> The deal.II project is located at http://www.dealii.org/
> For mailing list/forum options, see
> https://groups.google.com/d/forum/dealii?hl=en
> ---
> You received this message because you are subscribed to the Google Groups
> "deal.II User Group" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to dealii+unsubscr...@googlegroups.com.
> To view this discussion on the web visit
> https://groups.google.com/d/msgid/dealii/0372d95b-ae52-457e-9400-9e8159469a62%40googlegroups.com
> <https://groups.google.com/d/msgid/dealii/0372d95b-ae52-457e-9400-9e8159469a62%40googlegroups.com?utm_medium=email_source=footer>
> .
> For more options, visit https://groups.google.com/d/optout.
>

-- 
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- 
You received this message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/CAKS2Wq74jFVqjSuJq_cJLysdd5oR9kHBRxqYaSY1hzC0EPJtUA%40mail.gmail.com.
For more options, visit https://groups.google.com/d/optout.


Re: [deal.II] Re: Compatibility of Petsc with step 18

2019-07-16 Thread Ramprasad R
Dear Daniel,
The version of deal.ii is 9.0.1 and the version of PETSc is 3.10.5. I also
changed the config file of deal.ii to have MPI PETSc on.

Thank you.

Best regards,
Ramprasad

On Tue, Jul 16, 2019 at 5:12 PM Daniel Arndt  wrote:

> Ramprasad,
>
> Let's first rule out some simple problems. Which version of deal.II are
> you using? Did you compile deal.II with MPI and PETSc support?
> We only define the namecpacePETScWrappers when DEAL_II_WITH_PETSC=ON. Can
> you provide us with the content of the detailed.log file in your build
> directory?
>
> Best,
> Daniel
>
> --
> The deal.II project is located at http://www.dealii.org/
> For mailing list/forum options, see
> https://groups.google.com/d/forum/dealii?hl=en
> ---
> You received this message because you are subscribed to the Google Groups
> "deal.II User Group" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to dealii+unsubscr...@googlegroups.com.
> To view this discussion on the web visit
> https://groups.google.com/d/msgid/dealii/1b94eabe-7cd4-4b79-8ffb-207234fb026c%40googlegroups.com
> 
> .
> For more options, visit https://groups.google.com/d/optout.
>

-- 
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- 
You received this message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/CAKS2Wq6k9FjcM4CV34raLPsshwgXui%3DTCq9n7xQ2s_VVU3Agmg%40mail.gmail.com.
For more options, visit https://groups.google.com/d/optout.


[deal.II] Compatibility of Petsc with step 18

2019-07-16 Thread Ramprasad R
Hello All,

I am getting an error when I try to compile the step-18. It states 
"PETScWrappers does not name a type". I wanted to know if anyone else has 
faced the same problem. And if there is a solution to this.

Thank you. 

-- 
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- 
You received this message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/2d25a910-8d32-4268-a8aa-818ff5416670%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.


Re: [deal.II] Fwd: Regarding Step 60 tutorial of deal.ii

2019-06-26 Thread Ramprasad R
Hi Daniel,

Thank you so much for your detailed Explanation. As you said earlier, I
have to do the eigen value analysis.  You explanation has given me an
approach to start with. I can not thank you enough.

Warm regards

Ramprasad

On Tue, Jun 25, 2019 at 6:27 PM Daniel Garcia-Sanchez <
danielgarsanc...@gmail.com> wrote:

>
>
> On Monday, June 24, 2019 at 11:36:32 PM UTC+2, Ramprasad R wrote:
>>
>> Hello Bangerth,
>>
>> The terms N_* are the resultant forces in the * direction and these
>> forces are calculated using the strains which in turn are calculated using
>> the displacements u. So the terms N_* are not constants, rather change with
>> each element. And these values directly depend on u.
>>
>>
>>>
> Hi Ramprasad,
>
> I think that you want to do an eigenvalue calculation (step-36)
>
> I think that first you have to do an static calculation before your
> eigenvalue calculation in order to obtain the values for N_* (step-8)
>
> As discussed in your paper, the typical eigenvalue problem for the elastic
> equation takes this form, where omega^2 is your eigenvalue and Phi the
> eigenvector
> (K - omega^2 M)*Phi = 0
>
> For the buckling case lambda is your eigenvalue.
> (K - lambda G)*Phi = 0
>
> The calculation of K can be found in step-8, step-62 (or other tutorials).
>
> I think that in order to calculate G you need the resulting strain of an
> static calculation. You can do the static calculation, store the strain in
> a temporary buffer and use that data to calculate G. step-18 shows you how
> to do this.
>
> Once you have K and G, you can do an eigenvalue calculation. step-36 shows
> you how to do an eigenvalue calculation. Note that in step36 the stiffness
> matrix is called A.
>
> The equation in step36 is
> (D-epsilon M)*Phi=0
> which is very similar to the buckling equation.
>
> Best,
> Daniel
>
> --
> The deal.II project is located at http://www.dealii.org/
> For mailing list/forum options, see
> https://groups.google.com/d/forum/dealii?hl=en
> ---
> You received this message because you are subscribed to the Google Groups
> "deal.II User Group" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to dealii+unsubscr...@googlegroups.com.
> To view this discussion on the web visit
> https://groups.google.com/d/msgid/dealii/d026c101-ffdd-4325-a1bb-f40c8283bbc9%40googlegroups.com
> <https://groups.google.com/d/msgid/dealii/d026c101-ffdd-4325-a1bb-f40c8283bbc9%40googlegroups.com?utm_medium=email_source=footer>
> .
> For more options, visit https://groups.google.com/d/optout.
>

-- 
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- 
You received this message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/CAKS2Wq6yWDpn6sRzNbRfeFyeSJPky8pUSRKdrX_ARQaEAULqyA%40mail.gmail.com.
For more options, visit https://groups.google.com/d/optout.


Re: [deal.II] Fwd: Regarding Step 60 tutorial of deal.ii

2019-06-24 Thread Ramprasad R
Hello Bangerth,

The terms N_* are the resultant forces in the * direction and these forces
are calculated using the strains which in turn are calculated using the
displacements u. So the terms N_* are not constants, rather change with
each element. And these values directly depend on u.

Thank you.

Regards,
Ramprasad

On Mon, Jun 24, 2019 at 10:56 PM Wolfgang Bangerth 
wrote:

> On 6/24/19 9:24 AM, Ramprasad R wrote:
> >
> > The following image has the exact formulation of the geometric stiffness
> > matrix. Which is denoted by G in the image.
>
> Ramprasad,
> I don't think I quite understand the connection step-60. In the context
> of that problem, for which of the two variables (u and lambda) do you
> want to compute that matrix? Assuming you mean that it is supposed to be
> for u, it is not difficult to write down what you need to do based on
> what you already learn in step-3 or step-4, for example. In fact, these
> two early programs already implement the first two terms of G_ij if you
> assume that N_x=N_y=N_xy=1.
>
> I am not familiar with the terms used in the page you show. Are the N_*
> factors constants? Are they functions of space?
>
> Best
>   W.
>
>
> --
> 
> Wolfgang Bangerth  email: bange...@colostate.edu
> www: http://www.math.colostate.edu/~bangerth/
>
> --
> The deal.II project is located at http://www.dealii.org/
> For mailing list/forum options, see
> https://groups.google.com/d/forum/dealii?hl=en
> ---
> You received this message because you are subscribed to the Google Groups
> "deal.II User Group" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to dealii+unsubscr...@googlegroups.com.
> To view this discussion on the web visit
> https://groups.google.com/d/msgid/dealii/53a1ebd3-fdd4-485d-2486-b70913095c55%40colostate.edu
> .
> For more options, visit https://groups.google.com/d/optout.
>

-- 
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- 
You received this message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/CAKS2Wq58q8AU%2BJ%3DKR2Y7aExCd%2BFa1onWQMHR6EJv6-ADST8UKQ%40mail.gmail.com.
For more options, visit https://groups.google.com/d/optout.


[deal.II] Fwd: Regarding Step 60 tutorial of deal.ii

2019-06-22 Thread Ramprasad R
Dear Sir/Madam,

I am a student in Technical University of Braunschweig, Germany. I am
trying to work on a project for the aerospace department of my university.
In this project I am trying to extend the step 60 tutorial into buckling
analysis. I am trying to find a function to calculate the geometric
stiffness matrix. It would be of great help if I could find some help
regarding this issue.

Hoping to hear from you soon.

Thank you.

Best regards,

Ramprasad Raghunath

-- 
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- 
You received this message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/CAKS2Wq4aScmkeXcTbd_%2BuBtCig_ZVfO%3DGTBrHwrAsczXY7%3DLRg%40mail.gmail.com.
For more options, visit https://groups.google.com/d/optout.