Re: [deal.II] Re: Function to output stress and strain

2018-11-02 Thread Yuxiang Wang
Hi Denis, Thank you for the pointer! This is very helpful. The below is my understanding (and please let me know if I am wrong): 1. If we do the L2 projection over a domain of elements, then it smooths the stress jumps in the sense of integration. The stress would be continuous over the

Re: [deal.II] Gradient of the shape function wrt natural coordinate system?

2018-10-14 Thread Yuxiang Wang
Exactly what I was looking for! Thank you so much Professor Bangerth. Have a great one, Shawn On Sun, Oct 14, 2018 at 1:47 PM Wolfgang Bangerth wrote: > On 10/14/2018 11:57 AM, Yuxiang Wang wrote: > > > > While I can use the shape_grad function to get the gradient of the sha

[deal.II] Gradient of the shape function wrt natural coordinate system?

2018-10-14 Thread Yuxiang Wang
Dear all, While I can use the shape_grad function to get the gradient of the shape function with respect to the real cell coordinates, would there be a way to get the shape_grad with respect to the reference natural coordinates? For example, for shell element (that I'm working on), I have the

Re: [deal.II] Is "-- Configuring incomplete, errors occurred!" OK?

2018-10-12 Thread Yuxiang Wang
n. > > Best, > Matthias > > > > On Fri, Oct 12, 2018, at 17:36 CDT, Yuxiang Wang < > yuxiang.shawn.w...@gmail.com> wrote: > > > Thank you Prof. Bangerth for the response! > > > > Actually I kept that window open so I still have the original message.

Re: [deal.II] Can GridIn read 3D quad mesh?

2018-10-04 Thread Yuxiang Wang
ciate this > contribution! > > Best, > Jean-Paul > > On 04 Oct 2018, at 08:54, Yuxiang Wang > wrote: > > Hi Jean-Paul, > > Thank you again for the fix and sorry to bother again. I wonder that, > would you think it'd be easy to have spacedim = 3 and dim = 1 im

Re: [deal.II] Can GridIn read 3D quad mesh?

2018-10-04 Thread Yuxiang Wang
t; box (I’d forgotten that someone else had already done much of the work >> necessary to support codim-1). >> >> I’ve made a patch to add this functionality. See here: >> https://github.com/dealii/dealii/pull/7217 >> Thank a lot for the grids and the minimal test case!

Re: [deal.II] Can GridIn read 3D quad mesh?

2018-09-20 Thread Yuxiang Wang
gt; > Best, > J-P > > On 20 Sep 2018, at 08:11, Yuxiang Wang > wrote: > > Oh I see! Thank you for explaining. My guess would be that the ABAQUS > equivalent is "node set". Please let me know if you'd like me to try that! > > Best, > Shawn > > On Wed

Re: [deal.II] Can GridIn read 3D quad mesh?

2018-09-20 Thread Yuxiang Wang
ay you please construct and give us a fully structured >> grid in codimension 1 (so something like a rectangle warped in the third >> dimension) that also has some boundary colouring? That would help a great >> deal with implementing this feature, and would also ensure that the >> assig

Re: [deal.II] Can GridIn read 3D quad mesh?

2018-09-19 Thread Yuxiang Wang
> deal with implementing this feature, and would also ensure that the > assignment of boundary ID’s is possible. > > Best, > Jean-Paul > > On 19 Sep 2018, at 01:38, Yuxiang Wang > wrote: > > Prof. Bangerth, > > Thank you for the prompt response! Yes please find attached the

Re: [deal.II] Re: Can GridIn read 3D quad mesh?

2018-09-18 Thread Yuxiang Wang
3> gridin; gridin.attach_triangulation(triangulation); std::ifstream f("mesh.inp"); gridin.read_abaqus(f); GridOut gridout; std::ofstream fout("mesh.vtu"); gridout.write_vtu(triangulation, fout); } ``` Shawn On Tue, Sep 18, 2018 at 4:28 PM Wolfgang Bangerth wrote: &

Re: [deal.II] Re: Can GridIn read 3D quad mesh?

2018-09-18 Thread Yuxiang Wang
18, 2018 at 1:56:19 PM UTC-7, Wolfgang Bangerth wrote: > > On 09/18/2018 02:40 PM, Yuxiang Wang wrote: > > > > Thank you for your help! I just tried a code snippet (modified step-49) > below > > and it does not work, either when I set the dim=2 or dim=3. > > If

[deal.II] Re: Can GridIn read 3D quad mesh?

2018-09-18 Thread Yuxiang Wang
Oh and another detail - the mesh is not really 2D mesh because the four nodes of a quad element may not be co-planar. Therefore, there does not exist a transform that can reduce the mesh to 2D. Shawn On Monday, September 17, 2018 at 12:07:43 PM UTC-7, Daniel Arndt wrote: > > Shawn, > > If I

Re: [deal.II] Re: Can GridIn read 3D quad mesh?

2018-09-18 Thread Yuxiang Wang
Hi Daniel, Thank you for your help! I just tried a code snippet (modified step-49) below and it does not work, either when I set the dim=2 or dim=3. ``` int main () { Triangulation<2> triangulation; GridIn<2> gridin; gridin.attach_triangulation(triangulation); std::ifstream

[deal.II] Can GridIn read 3D quad mesh?

2018-09-16 Thread Yuxiang Wang
Dear all, Sorry for the beginner question. I am trying to load a 3D quad mesh (i.e. not all nodes are co-planar for one element) to implement a Mindlin–Reissner plate. I wonder, can GridIn class read such meshes? Upon reading the document, it seems like if the input element is quad then it is

Re: [deal.II] Re: Help needed - remote debug in VS2017 hangs when MPI_InitFinalize is called

2018-08-02 Thread Yuxiang Wang
to everyone who helped! Shawn On Friday, July 13, 2018 at 4:59:56 PM UTC-7, Wolfgang Bangerth wrote: > > On 07/12/2018 04:36 PM, Yuxiang Wang wrote: > > > > I tried and I can neither step over (n) nor step into (s) that line of > code - > > either way, it hangs rig

Re: [deal.II] Re: Interior node locations for higher order elements

2018-04-20 Thread Yuxiang Wang
Hi Bruno, Exactly what I was looking for. Thank you so much! Best, Shawn On Fri, Apr 20, 2018 at 1:29 PM, Bruno Turcksin <bruno.turck...@gmail.com> wrote: > Shawn, > > On Friday, April 20, 2018 at 2:31:28 PM UTC-4, Yuxiang Wang wrote: >> >> Could you pl

[deal.II] Interior node locations for higher order elements

2018-04-20 Thread Yuxiang Wang
Hello dear deal.II community, Previously, I thought that for a high order element in 1D (cubic, for example), the interior nodes would be equally distributed at: Node 0 - 0 Node 1 - 0.01 Node 2 - 0.0033 Node 3 - 0.0067 However, in deal.II I found them to be at: Node 0 - 0 Node 1 -

Re: [deal.II] Question wrt step-3: shape functions and test functions

2018-04-07 Thread Yuxiang Wang
Thank you Professor Bangerth for the answer! I appreciate your link to the video a lot. And thank you for the kind words - I'll try to work harder to match your expectations and one day contribute to deal.II :) Have a great one, Shawn On Monday, April 2, 2018 at 1:22:42 PM UTC-7, Wolfgang

[deal.II] Question wrt step-3: shape functions and test functions

2018-03-31 Thread Yuxiang Wang
Dear all, Sorry for the spam - a basic mathematical question. In the step-3 example (https://www.dealii.org/8.4.1/doxygen/deal.II/step_3.html), I noticed that at first the variable phi is used as a test function. It is in the same space with u, except for that the value is 0 on the

Re: [deal.II] Are shell elements available in deal.II?

2018-03-30 Thread Yuxiang Wang
on why you can’t implement Mindlin-Reissner plate theory > in deal.II. I would start with Timoshenko beam theory and then build up to > plates. > > A > > > > On 26 Mar 2018, at 14:40, Wolfgang Bangerth <bang...@colostate.edu > > wrote: > > On 03/25/2018 11:44 A

Re: [deal.II] Are shell elements available in deal.II?

2018-03-25 Thread Yuxiang Wang
Hi Professor Bangerth, Thank you for the pointer! Those are very helpful. I just realized that Mindlin-Reissner plate theory can be expressed in forms of mixed elements. My current need for shell element is for solid mechanics problem (Mindlin-Reissner plates). A simple example could be that a

[deal.II] Re: For docker, cmake runs examples fine when I run it in the pre-installed src example folder, but not after I copied it to ~/shared

2018-03-23 Thread Yuxiang Wang
uot;/home/dealii/shared/examples/step-1/CMakeFiles/CMakeError.log". ``` On Thursday, March 22, 2018 at 11:58:34 AM UTC-7, Yuxiang Wang wrote: > > Short summary: cmake runs examples fine when I run it in the pre-installed > src example folder, but not after I copied it to ~/shared. Co

[deal.II] For docker, cmake runs examples fine when I run it in the pre-installed src example folder, but not after I copied it to ~/shared

2018-03-22 Thread Yuxiang Wang
Short summary: cmake runs examples fine when I run it in the pre-installed src example folder, but not after I copied it to ~/shared. Could you please help see why? Full story: first, when I run it in the pre-installed src example folder, it works: ``` E:\dealii>docker pull

Re: [deal.II] Are shell elements available in deal.II?

2018-03-06 Thread Yuxiang Wang
Thank you Professor Bangerth for your response! Yes that's exactly what I am interested in. I am excited to learn that there are arguably better methods (mixed methods) to solve shell problems. I am not familiar with those; could you please kindly provide some pointers for me to start

[deal.II] Are shell elements available in deal.II?

2018-03-05 Thread Yuxiang Wang
Hi, Sorry for the spam. I tried to search but did not find an implementation of shell elements in deal.II. Since this is a commonly used element, I'd like to make sure that it's not me missing it. Could you please help confirm? Best, Shawn -- The deal.II project is located at