Element.quad4i

From Nemesis
Revision as of 16:34, 9 January 2008 by Fotis (Talk | contribs)

(diff) ←Older revision | Current revision (diff) | Newer revision→ (diff)
Jump to: navigation, search

The following article describes the use of the element.quad4i command. This command defines a four-node non-conforming quad based on the work of R.L. Taylor, P.J Beresford & E.L. Wilson<ref>R.L. Taylor, P.J Beresford and E.L. Wilson, "A non-conforming element for stress analysis", International Journal for Numerical Methods in Engineering, Vol. 10, pp.1211-1219 (1976)</ref>.

Contents

Syntax

element.quad4i(id,node1,node2,node3,node4,mat)

where

Parameter description type default
id The elemental id integer -
node1-4 The elemental nodes integer -
mat The element material integer -

Examples

Mesh distortion test

Mesh distortion test.

The example code can be found here.

===========================
|alpha|STD      |QM6      |
+-----+---------+---------+
| 0.00|   0.2424|   1.0000|
| 1.00|   0.1274|   0.6174|
| 2.00|   0.0846|   0.5337|
| 3.00|   0.0703|   0.5271|
| 4.00|   0.0597|   0.5044|
===========================

Patch Test

Patch test.

A typical patch test. The outer nodes are ordered anti-clockwise as 1(0.0, 0.0), 2(9.0, 0.0), 3(8.0, 8.0) and 4(0.0, 6.0) while the inner nodes are arranged similarly as 5(2.5, 2.0), 6(6.5, 1.5), 7(5.5, 5.0) and 8(2.5, 5.0). A displacement field is given in the outer nodes defined by ux=0.0001+0.0061x+0.0049y and uy=-0.0005+0.0042x+0.0038y, where x and y are the nodal coordinates. The calculated displacements in the inner nodes must coincide with the previous defined displacement field.

The displacement field (displacement vector magnitude) is plotted in the picture on the left, while the results are shown in the next table.

The example code can be found here.

=================================================
|   |     Analytical      |     Numerical       |
+---+----------+----------+----------+----------+
|id |ux        |uy        |ux        |uy        |
+---+----------+----------+----------+----------+
|  5|  0.025150|  0.017600|  0.025150|  0.017600|
|  6|  0.047100|  0.032500|  0.047100|  0.032500|
|  7|  0.058150|  0.041600|  0.058150|  0.041600|
|  8|  0.039850|  0.029000|  0.039850|  0.029000|
=================================================

Comments

Nodes, material and section must be already defined. This element can be only used in plane strain/stress analysis.

Theory


Implementation

The implementation of this element follows closely the work of T.J.R. Hughes<ref>T.J.R. Hughes, "The Finite Element Method: Linear Static and Dynamic Finite Element Analysis"', Prentice-Hall Inc., Englewood Cliffs, New Jersey, 1997, pp.243-248.</ref>.

Notes

<references/>


See also

Personal tools