[vtk-developers] Pentagonal prism bug resolved!

Tue Jun 16 20:20:05 EDT 2015

Brilliant work Bill,
   That is a really interesting paper. I'm glad you found a solution.
Your work here just underscores the fundamental importance of good testing.
I'm looking forward to seeing it in VTK.

I had never heard of SymPy, I had always used Maple but it is too expensive
now that I have retired. I'll give it a SymPy a go.

Nothing wrong with punch cards ... except when you drop a whole box. I used
them to do Tempromandibular joint (TMJ) analysis using FORTRAN, way back in
the 70's. Then I graduated to using a Teletype and APL!

Regards
    Andrew

On Tue, Jun 16, 2015 at 6:26 PM, <vtk-developers-request at vtk.org> wrote:

> Send vtk-developers mailing list submissions to
>         vtk-developers at vtk.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
>         http://public.kitware.com/mailman/listinfo/vtk-developers
> or, via email, send a message with subject or body 'help' to
>         vtk-developers-request at vtk.org
>
> You can reach the person managing the list at
>         vtk-developers-owner at vtk.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of vtk-developers digest..."
>
> Today's Topics:
>
>    1. Pentagonal prism bug resolved! (Bill Lorensen)
>    2. Re: Pentagonal prism bug resolved! (Berk Geveci)
>    3. Re: Pentagonal prism bug resolved! (Bill Lorensen)
>    4. Re: Pentagonal prism bug resolved! (Berk Geveci)
>    5. Re: VTK for android example: missing      vtkRenderingOpenGL2
>       (Lonni Besan?on)
>    6. android_native_app_glue.c doesn't exist (Lonni Besan?on)
>
>
> ---------- Forwarded message ----------
> From: Bill Lorensen <bill.lorensen at gmail.com>
> To: VTK Developers <vtk-developers at vtk.org>
> Cc:
> Date: Mon, 15 Jun 2015 13:31:58 -0400
> Subject: [vtk-developers] Pentagonal prism bug resolved!
> Folks,
>
> I reported recently of issues with the shape functions and derivatives
> for the pentagonal prism. Since I could not find any references for
> the current implementation I used the formulation in this paper:
> http://dilbert.engr.ucdavis.edu/~suku/nem/papers/polyelas.pdf
> Since this paper did not have the shape function derivatives, I used
> SymPy ( http://www.sympy.org/en/index.html ) to compute the
> derivatives.
>
> I'm happy to report that the new formulation does a great job and
> passes the new unit test I created for cells.
>
> I have not submitted a patch yet since I want the code to document the
> process. I wrote c++ code and SymPy code to generate the statements
> for the new formulation.
>
> I wanted to give a heads up in case someone else is looking into this.
>
> Bill
>
>
>
> ---------- Forwarded message ----------
> From: Berk Geveci <berk.geveci at kitware.com>
> To: Bill Lorensen <bill.lorensen at gmail.com>
> Cc: VTK Developers <vtk-developers at vtk.org>
> Date: Mon, 15 Jun 2015 14:17:00 -0400
> Subject: Re: [vtk-developers] Pentagonal prism bug resolved!
> You rock Bill.
>
> On Mon, Jun 15, 2015 at 1:31 PM, Bill Lorensen <bill.lorensen at gmail.com>
> wrote:
>
>> Folks,
>>
>> I reported recently of issues with the shape functions and derivatives
>> for the pentagonal prism. Since I could not find any references for
>> the current implementation I used the formulation in this paper:
>> http://dilbert.engr.ucdavis.edu/~suku/nem/papers/polyelas.pdf
>> Since this paper did not have the shape function derivatives, I used
>> SymPy ( http://www.sympy.org/en/index.html ) to compute the
>> derivatives.
>>
>> I'm happy to report that the new formulation does a great job and
>> passes the new unit test I created for cells.
>>
>> I have not submitted a patch yet since I want the code to document the
>> process. I wrote c++ code and SymPy code to generate the statements
>> for the new formulation.
>>
>> I wanted to give a heads up in case someone else is looking into this.
>>
>> Bill
>> _______________________________________________
>> Powered by www.kitware.com
>>
>> Visit other Kitware open-source projects at
>> http://www.kitware.com/opensource/opensource.html
>>
>> Search the list archives at: http://markmail.org/search/?q=vtk-developers
>>
>> Follow this link to subscribe/unsubscribe:
>> http://public.kitware.com/mailman/listinfo/vtk-developers
>>
>>
>
>
> ---------- Forwarded message ----------
> From: Bill Lorensen <bill.lorensen at gmail.com>
> To: Berk Geveci <berk.geveci at kitware.com>
> Cc: VTK Developers <vtk-developers at vtk.org>
> Date: Mon, 15 Jun 2015 14:27:34 -0400
> Subject: Re: [vtk-developers] Pentagonal prism bug resolved!
> Brings back memories of my finite element days, 1968-1978, at
> Watervliet Arsenal. At that time we worked with NASTRAN and we added a
> few elements for crack tip analysis....
>
>
> Pu, S. L., M. A. Hussain, and W. E. Lorensen, The Collapsed Cubic
> Isoparametric Element as a Singular Element for Crack Problems,
> International Journal for Numerical Methods in Engineering, vol. 12,
> no. 11, pp. 1727-1742, 1978.
>
> Hussain, M. A. and W. E. Lorensen, Isoparametric Elements as Singular
> Elements for Crack Problems, AR-LCB-TR-77041, Benet Weapons
> Laboratory, Watervliet, NY, November 1977.
>
> Hussain, M. A.; Pu, S. L.; Lorensen, W. E., Singular Plastic Element:
> NASTRAN Implementation and Application, Sixth NASTRAN Users' Colloq.,
> pp. 257-274, October, 1977, (NASA CP-2018).
>
> Pu, S. L., M. A. Hussain, and W. E. Lorensen, Collapsed 12-Node
> Triangular Elements as Crack Tip Elements for Elastic Fracture,
> ARLCB-TR-77047, Watervliet Arsenal, Watervliet, NY, December 1977.
>
> Lorensen, W. E. and M. A. Hussain, "The Stability of Isoparametric
> Elements as a Singular Element for Crack Problems," in Developments in
> Mechanics, ed. T. C. T. Ting, 1977.
>
> Pu, S. L., M. A. Hussain, and W. E. Lorensen, The Collapsed Cubic
> Isoparametric Element as a Singular Element for Crack Problems,
> AR-LCB-TR77047, Watervliet Arsenal, 1977.
>
> Hussain, M. A., W. E. Lorensen, and G. Pflegl, The Quarter-Point
> Quadratic Isoparametric Element as A Singular Element for Crack
> Problems, TM-X-3428, p. 419, NASA, 1976.
>
>
> On Mon, Jun 15, 2015 at 2:17 PM, Berk Geveci <berk.geveci at kitware.com>
> wrote:
> > You rock Bill.
> >
> > On Mon, Jun 15, 2015 at 1:31 PM, Bill Lorensen <bill.lorensen at gmail.com>
> > wrote:
> >>
> >> Folks,
> >>
> >> I reported recently of issues with the shape functions and derivatives
> >> for the pentagonal prism. Since I could not find any references for
> >> the current implementation I used the formulation in this paper:
> >> http://dilbert.engr.ucdavis.edu/~suku/nem/papers/polyelas.pdf
> >> Since this paper did not have the shape function derivatives, I used
> >> SymPy ( http://www.sympy.org/en/index.html ) to compute the
> >> derivatives.
> >>
> >> I'm happy to report that the new formulation does a great job and
> >> passes the new unit test I created for cells.
> >>
> >> I have not submitted a patch yet since I want the code to document the
> >> process. I wrote c++ code and SymPy code to generate the statements
> >> for the new formulation.
> >>
> >> I wanted to give a heads up in case someone else is looking into this.
> >>
> >> Bill
> >> _______________________________________________
> >> Powered by www.kitware.com
> >>
> >> Visit other Kitware open-source projects at
> >> http://www.kitware.com/opensource/opensource.html
> >>
> >> Search the list archives at:
> http://markmail.org/search/?q=vtk-developers
> >>
> >> Follow this link to subscribe/unsubscribe:
> >> http://public.kitware.com/mailman/listinfo/vtk-developers
> >>
> >
>
>
>
> --
> Unpaid intern in BillsBasement at noware dot com
>
>
>
> ---------- Forwarded message ----------
> From: Berk Geveci <berk.geveci at kitware.com>
> To: Bill Lorensen <bill.lorensen at gmail.com>
> Cc: VTK Developers <vtk-developers at vtk.org>
> Date: Mon, 15 Jun 2015 15:03:09 -0400
> Subject: Re: [vtk-developers] Pentagonal prism bug resolved!
> Finite element analysis with punch cards. That's hard core!
>
> -berk
>
> On Mon, Jun 15, 2015 at 2:27 PM, Bill Lorensen <bill.lorensen at gmail.com>
> wrote:
>
>> Brings back memories of my finite element days, 1968-1978, at
>> Watervliet Arsenal. At that time we worked with NASTRAN and we added a
>> few elements for crack tip analysis....
>>
>>
>> Pu, S. L., M. A. Hussain, and W. E. Lorensen, The Collapsed Cubic
>> Isoparametric Element as a Singular Element for Crack Problems,
>> International Journal for Numerical Methods in Engineering, vol. 12,
>> no. 11, pp. 1727-1742, 1978.
>>
>> Hussain, M. A. and W. E. Lorensen, Isoparametric Elements as Singular
>> Elements for Crack Problems, AR-LCB-TR-77041, Benet Weapons
>> Laboratory, Watervliet, NY, November 1977.
>>
>> Hussain, M. A.; Pu, S. L.; Lorensen, W. E., Singular Plastic Element:
>> NASTRAN Implementation and Application, Sixth NASTRAN Users' Colloq.,
>> pp. 257-274, October, 1977, (NASA CP-2018).
>>
>> Pu, S. L., M. A. Hussain, and W. E. Lorensen, Collapsed 12-Node
>> Triangular Elements as Crack Tip Elements for Elastic Fracture,
>> ARLCB-TR-77047, Watervliet Arsenal, Watervliet, NY, December 1977.
>>
>> Lorensen, W. E. and M. A. Hussain, "The Stability of Isoparametric
>> Elements as a Singular Element for Crack Problems," in Developments in
>> Mechanics, ed. T. C. T. Ting, 1977.
>>
>> Pu, S. L., M. A. Hussain, and W. E. Lorensen, The Collapsed Cubic
>> Isoparametric Element as a Singular Element for Crack Problems,
>> AR-LCB-TR77047, Watervliet Arsenal, 1977.
>>
>> Hussain, M. A., W. E. Lorensen, and G. Pflegl, The Quarter-Point
>> Quadratic Isoparametric Element as A Singular Element for Crack
>> Problems, TM-X-3428, p. 419, NASA, 1976.
>>
>>
>> On Mon, Jun 15, 2015 at 2:17 PM, Berk Geveci <berk.geveci at kitware.com>
>> wrote:
>> > You rock Bill.
>> >
>> > On Mon, Jun 15, 2015 at 1:31 PM, Bill Lorensen <bill.lorensen at gmail.com
>> >
>> > wrote:
>> >>
>> >> Folks,
>> >>
>> >> I reported recently of issues with the shape functions and derivatives
>> >> for the pentagonal prism. Since I could not find any references for
>> >> the current implementation I used the formulation in this paper:
>> >> http://dilbert.engr.ucdavis.edu/~suku/nem/papers/polyelas.pdf
>> >> Since this paper did not have the shape function derivatives, I used
>> >> SymPy ( http://www.sympy.org/en/index.html ) to compute the
>> >> derivatives.
>> >>
>> >> I'm happy to report that the new formulation does a great job and
>> >> passes the new unit test I created for cells.
>> >>
>> >> I have not submitted a patch yet since I want the code to document the
>> >> process. I wrote c++ code and SymPy code to generate the statements
>> >> for the new formulation.
>> >>
>> >> I wanted to give a heads up in case someone else is looking into this.
>> >>
>> >> Bill
>> >> _______________________________________________
>> >> Powered by www.kitware.com
>> >>
>> >> Visit other Kitware open-source projects at
>> >> http://www.kitware.com/opensource/opensource.html
>> >>
>> >> Search the list archives at:
>> http://markmail.org/search/?q=vtk-developers
>> >>
>> >> Follow this link to subscribe/unsubscribe:
>> >> http://public.kitware.com/mailman/listinfo/vtk-developers
>> >>
>> >
>>
>>
>>
>> --
>> Unpaid intern in BillsBasement at noware dot com
>>
>
>
>
>
-- 
___________________________________________
Andrew J. P. Maclean

___________________________________________
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/vtk-developers/attachments/20150617/d56aaa6c/attachment-0001.html>