Revision of Problem with fortran code for solving large strain 2D elastic problem by FEM using Updated Lagrangian Framework from Sun, 2009-12-13 10:05

The revisions let you track differences between multiple versions of a post.

Suman.Guha's picture

Normal
0
false
false
false
EN-US
X-NONE
X-NONE
MicrosoftInternetExplorer4

/* Style Definitions */
table.MsoNormalTable
{mso-style-name:"Table Normal";
mso-tstyle-rowband-size:0;
mso-tstyle-colband-size:0;
mso-style-noshow:yes;
mso-style-priority:99;
mso-style-qformat:yes;
mso-style-parent:"";
mso-padding-alt:0in 5.4pt 0in 5.4pt;
mso-para-margin-top:0in;
mso-para-margin-right:0in;
mso-para-margin-bottom:10.0pt;
mso-para-margin-left:0in;
line-height:115%;
mso-pagination:widow-orphan;
font-size:11.0pt;
font-family:"Calibri","sans-serif";
mso-ascii-font-family:Calibri;
mso-ascii-theme-font:minor-latin;
mso-fareast-font-family:"Times New Roman";
mso-fareast-theme-font:minor-fareast;
mso-hansi-font-family:Calibri;
mso-hansi-theme-font:minor-latin;}

I am using my own fortran code to solve a large strain 2D elastic
problem using FEM in Updated Lagrangian framework and I am using incremental
formulation. It works well for 2D pull of a square block with 2 elements and I've
compared the results with a commercial FEA package. But while solving bending
of a cantilever beam the results starts to blow up and the elements distort
abruptly after certain time step. This happens only when the load or the
displacement given at the free end of the beam is large otherwise it works
fine.

Basically, I am supposed to solve a gradient visco-plasticity
problem using the formulation given by Niordson and Fleck (see the paper
attached herewith). I am facing the same problem when using the viscoplastic
formulation. The results start to blow up after few steps only.

Does anyone have any idea what could be the possible
problem?

Thanks in advance.