Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

Mindlin Plate Theory and Abaqus UEL

Implementation

Page 1 of 97
 Nikhil Padhye

[email protected]

Subodh Kalia

[email protected]

2
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

Contents

1 Mindlin Plate Theory 7

1.1 Stress and strain components . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8

1.2 Shape functions and isoparametric element formulation . . . . . . . . . . . . . 13

1.3 Computation of Jacobian . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15

1.4 Calculation of stress and strain at integration points . . . . . . . . . . . . . . . . 18

1.5 Principle of virtual work and derivation of element stiffness matrix . . . . . . . 24

2 Sample Hand Calculations 29

2.1 Bending Stiffness . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30

2.2 Shear Stiffness . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34

3 Some Consequences of Mindlin Plate Theory 37

4 ABAQUS UEL Implementation 40

4.1 UEL input variables . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41

4.2 UEL output variables . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42

4.3 Pseudo code for the UEL subroutine . . . . . . . . . . . . . . . . . . . . . . . . . . 42

5 Benchmarking of UEL Implementation for Mindlin Plate Element 50

5.1 Single element test . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50

5.2 Testing multiple elements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56

6 Acknowledgments 59

Page 3 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

List of Figures

1 W , L and t represent the width, length and thickness of the plate, respectively. 7

2 Rotation of a material line (normal to the neutral plane in the undeformed

configuration), about Y-axis in anti-clockwise direction when viewed from the

positive Y-axis, and denoted by θ y . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10

3 Isoparametric formulation for 4-noded element using 2x2 quad integration. . . 14

4 Four-noded element with integration points in isoparametric formulation. . . . 30

5 UEL subroutine header from Abaqus 2016 documentation, Abaqus User Sub-

routines Reference Guide, Section 1.1.28. . . . . . . . . . . . . . . . . . . . . . . . 40

6 Input-Output block diagram for the UEL subroutine. . . . . . . . . . . . . . . . . 41

7 Boundary conditions for a single element test. . . . . . . . . . . . . . . . . . . . . 50

8 Reaction forces in Z-direction with Abaqus S4 element. . . . . . . . . . . . . . . 51

9 Reaction moments in X-direction with Abaqus S4 element. . . . . . . . . . . . . 52

10 Reaction moments in Y-direction with Abaqus S4 element. . . . . . . . . . . . . 52

11 Plots for reaction forces in Z-direction at the nodes with respect to nodal dis-

placement. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54

12 Plots for reaction moments in X-direction at the nodes with respect to nodal

displacements. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54

13 Plots of reaction moments in Y-direction at the nodes with respect to nodal

displacements. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55

14 Plate deflection from Matlab code. . . . . . . . . . . . . . . . . . . . . . . . . . . . 57

15 Plate deflection from Abaqus S4 element analysis with a scale factor of 106 . . . 57

Page 4 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

16 Plate deflection from UEL with a scale factor of 106 . . . . . . . . . . . . . . . . . . 58

17 Comparison of deflections for a laterally loaded plate clamped at the edges. . . 58

Page 5 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

A BSTRACT

In this report, we present a model for plate bending based on the Mindlin plate theory for

small elastic deformations. The related field and constitutive equations, and finite element

discretization (FEM) at an element level are also presented. The implementation of the

Mindlin plate element is carried out for the UEL module of ABAQUS. The overall goal of

this technical report is to facilitate the understanding and implementation of finite element

discretization for the UEL module of Abaqus, so that users can define their own elements. The

overall procedure described in this technical report is expected to aid the implementation

of custom defined user elements for ABAQUS, a topic that is currently not discussed well

in the opinion of the authors. For example, other than official Abaqus documentation [1],

only one comprehensive documentation for UEL in Abaqus is available online [9]. The

current demonstration is carried out for small deformations, for the sake of simplicity, and

applications involving other plate/shell theories (along with large deformation or rotation

effects) will be done in the subsequent reports and publications 1 . Although benchmarking

on simple test problems is performed in this study, detailed analyses on the performance of

current plate formulation such as patch test performance, shear locking issues, etc. are not

discussed. Reader is referred to [3], [7], [4], [2] and [6] for detailed discussions on these topics,

including more sophisticated finite element methods.

1
The readers are encouraged to contact the authors in case they need any further clarification. Use of this
document or any content/code presented here for any commercial purposes is strictly prohibited. Readers are
welcome to cite this report if it is found useful to them in anyway.

Page 6 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

1 M INDLIN P LATE T HEORY

Top surface

z L
t

y
Mid surface

Figure 1: W , L and t represent the width, length and thickness of the plate, respectively.

A plate is commonly identified as a thin flat structure as shown in Figure 1. In order

to qualify as a plate, we require that the thickness (t ) of the structure is small compared

to its length (L) and width (W ), i.e. t << L,W . This assumption simplifies the underlying

governing equilibrium equations and also allows for several simplifications on the kinematics

of deformation due to the “thin" nature of the structure at hand. The Mindlin plate theory

can account for homogeneous through thickness shear deformations (as opposed to another

elementary theory of Kirchhoff-Love plates).

The two key features of this theory are: (a) we assume by construction that σz = 0, and

(b) since this is a first order shear deformation theory; therefore, allowing σxz and σ y z to be

non-zero. We estimate σxz , σ y z from γxz , γ y z which themselves are taken to be non-zero and

constant through the thickness of the plate. Membrane effects (i.e. in plane stretching of the

plate mid-surface) are not included in this theory. Some other consequences of such ad-hoc

Page 7 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

assumptions made in this theory are discussed later in Section 3.

1.1 Stress and strain components

The deformation in the Mindlin plate can be characterized by the displacement of the mid-

surface of the plate in the Z-direction (denoted by w), and rotation of the normal to the

mid-surface (denoted by θx and θ y ) about X-axis and Y-axis, respectively. Thus we require

three degrees of freedom w, θx , θ y at any point on the mid-surface of the plate to describe the

complete deformation.

Let us consider a material line L 1 L 2 , across thickness of the plate, perpendicular to the

neutral plane in the undeformed configuration as shown in Figure 2. Upon deformation

we assume that points on the material line L 1 L 2 still remain on a straight line l 1 l 2 in the

deformed configuration. The material line L 1 L 2 is rotated by an angle θ y in the anti-clockwise

direction with respect to the Y-axis. For small deformations, the displacement of any material

point lying on the line L 1 L 2 , relative to the point ‘O’ on the mid-surface, in the X-direction,

can be approximated as 2 :

u x = zθ y , (1)

The extensional strain in the X-direction is given as:

∂u x ∂θ y
²x = =z . (2)
∂x ∂x
2
We have assumed no membrane action on the plate, i.e., no in-plane strains on the mid-surface of the
plate, therefore the mid-surface displacement does not vary spatially, nor does it plays any role in strains across
thickness as it shall merely enter as a constant sum to u x in equation 2.

Page 8 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

In the above equation, z indicates the coordinate of a material point along Z-axis in

the undeformed configuration, and for the sake of clarity we repeat that θ y is measured in

anti-clockwise direction about Y-axis when observing from positive direction to the negative

direction.

It is worth highlighting that a rather tacit assumption has been made in this formulation

(while locating a material point along the thickness direction) - regarding no change in

thickness (implying ²z = 0), which seems inconsistent with the prior assumption of σz = 0. In

opinion of the authors, it is quite important to understand the limitations of such classical

plate/shell theories (even with their more elegant extensions as made in the literature), since

they can lead to rather spurious predictions in scenarios where the underlying assumptions

are quite an overreach.

Following a similar procedure one can estimate u y = −zθx . Please note the minus sign in

the estimation of u y , which was not present in the case of u x . Thus, ² y is given as

∂u y ∂θx
²y = = −z (3)
∂y ∂y

and accordingly the shear strains in the XY-plane are given as

∂θx ∂θ y
· ¸
γx y = z − + . (4)
∂x ∂y

The shear strains across the thickness (γxz and γ y z ) can be calculated as

∂u x ∂u z
γxz = + , (5)
∂z ∂x

Page 9 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

N
θy
L1

‘O’
L2
z
Top surface
N
z L1
y ux = z θy
x ‘O’
L2

Bottom surface Center line

Figure 2: Rotation of a material line (normal to the neutral plane in the undeformed configuration),
about Y-axis in anti-clockwise direction when viewed from the positive Y-axis, and denoted by θ y .

Page 10 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

substituting u x = z θ y and u z = w(x, y) in equation 5, we get

∂w(x, y)
γxz = θ y (x, y) + . (6)
∂x

Similarly, we can obtain

∂w(x, y)
γ y z = −θx (x, y) + . (7)
∂y

Now we define a bending strain vector ²b and shear strain vector ²s as follows:

 
 ²x 
 
 
²b =  ² y , (8)
 
 
 
 
γx y

 
 γxz 
²s =  . (9)
 

γy z
 

From linear elasticity we have

1£
²z = σz − ν(σ y + σx ) ,
¤
(10)
E

where E is Young’s modulus and ν is Poisson’s ratio. Under the assumption of plane stress,

i.e. σz = 0, we get

−ν
²z = (σx + σ y ). (11)
E

Page 11 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

Similarly, we get

1£
²x = σx − νσ y ,
¤
(12)
E

1£
²y = σ y − νσx .
¤
(13)
E

The shear strain γx y is given as

γx y = Gσx y , (14)

E
where G = 2(1+ν) is the shear modulus. Solving for σx and σ y from equations 12 and 13,

we get

E £
σx = ² ν²
¤
x + y , (15)
1 − ν2

E £
σy = ν²x + ² y .
¤
(16)
1−ν2

Now we assemble all the (x and y) stress and strain components (and refer to them as

bending components):

    
 σx   1 ν 0   ²x 
    
  E   
σb =  σ y = (17)
 1−ν  ν 1   ²y
    
 2
 0   

    
    
1−ν
σx y 0 0 2
γx y

Page 12 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

The through thickness shear components of are given as

σxz = α G γxz (18)

σy z = α G γy z (19)

where α is some shear correction factor, and we choose α = 5/6 in this study. Shear stress

components in terms of shear strain components now can be assembled and written as

    
 σxz   1 0   γxz 
σs =  =G . (20)
    
 
σy z γy z
    
0 1

Thus far we have formulated stress and strain components for both the bending and

shear separately, and this was possible due to several kinematic assumptions and small

strain setting. Next, we present the finite element discretization at an element level such

that Mindlin theory can be simulated. We once again wish to highlight that the current finite

element plate formulation is rather naive, and more effective and consistent proposals have

been made in the literature. Our goal here is to focus on the simplicity so that implementation

aspects for UEL in Abaqus can be demonstrated easily.

1.2 Shape functions and isoparametric element formulation

We will use a four-noded bilinear isoparametric element where each node has three degrees

of freedom (w, θx and θ y ). Four integration points (as per 2x2 quad integration rule) are used

Page 13 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

1 1 1 1
- , ,
3 3
Node 4 3 3
η
Node 3
(x4, y4) y (-1, 1) (1, 1)
(x3, y3)
z x x
ζ
x x
x
(-1, -1) (1, -1)

Node 1 Node 2 -
1
3
,-
1
3
1
3
,-
1
3
(x1, y1) (x2, y2)

Figure 3: Isoparametric formulation for 4-noded element using 2x2 quad integration.

to capture bending response, and single integration point (reduced integration) is used to

account for shear response. The shape functions for this isoparametric element (with ξ and η

as isoparametric coordinates) are given as:

1
N1 = (1 − ξ)(1 − η) , (21)
4

1
N2 = (1 + ξ)(1 − η) , (22)
4

1
N3 = (1 + ξ)(1 + η) , (23)
4

1
N4 = (1 − ξ)(1 + η) . (24)
4

Page 14 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

We denote the shape function vector as N (which has a size of 4x1 and suppress ξ and η

for the sake of convenience in notation):

 
 N1 
 
 
 
 N2 
N = . (25)
 
 
 N 
 3 
 
 
N4

The derivatives of shape functions, denoted by d N , with respect to isoparametric coordi-

nates are written in a matrix form as:

   
∂N1 ∂N1
 ∂ξ ∂η   −(1 − η) −(1 − ξ) 
   
   
∂N2 ∂N2
 1  (1 − η) −(1 + ξ)
   
∂ξ ∂η
 
dN =  =  . (26)
   
 ∂N3 ∂N3
 4 
   (1 + η) (1 + ξ) 
 ∂ξ ∂η   
   
∂N4 ∂N4
   
∂ξ ∂η
−(1 + η) (1 − ξ)

1.3 Computation of Jacobian

The derivatives of the shape functions (Ni ) with respect to the isoparametric coordinates can

be written as

∂Ni ∂Ni ∂x ∂Ni ∂y

= + , (27)
∂ξ ∂x ∂ξ ∂y ∂ξ

∂Ni ∂Ni ∂x ∂Ni ∂y

= + . (28)
∂η ∂x ∂η ∂y ∂η

Page 15 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

By combining equation 27 and 28 in a matrix form we obtain

    
∂Ni ∂x ∂y ∂Ni

 ∂ξ  
  ∂ξ ∂ξ 
 ∂x 

    
= . (29)
    
 
    
    
∂Ni ∂x ∂y ∂Ni
    
∂η ∂η ∂η ∂y

Now we define the Jacobian as

 
∂x ∂y

 ∂ξ ∂ξ 

 
J = (30)
 

 
 
∂x ∂y
 
∂η ∂η

where x and y are spatial coordinates of any point within the element. Using equation

30, the equation 29 can be re-written as

   
∂Ni ∂Ni

 ∂x 


 ∂ξ 

   
 = J −1  (31)
   
 
   
   
∂Ni ∂Ni
   
∂y ∂η

Now we can assemble a matrix which contains the spatial derivatives of all the shape

functions in terms of the Jacobian and isoparametric shape function derivatives:

   
∂N1 ∂N2 ∂N3 ∂N4 ∂N1 ∂N2 ∂N3 ∂N4
∂x ∂x ∂x ∂x ∂ξ ∂ξ ∂ξ ∂ξ
   
−1 
= J   = J −1 d N > . (32)
  

∂N1 ∂N2 ∂N3 ∂N4 ∂N1 ∂N2 ∂N3 ∂N4
   
∂y ∂y ∂y ∂y ∂η ∂η ∂η ∂η

We see from equation 31, that the derivatives of the shape functions with the x, y spatial

coordinates can be obtained from the derivates of the shape functions with respect to the

Page 16 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

isoparametric coordinates (ξ, η) if the Jacobian is known. Next we derive more explicit

expression for the Jacobian in terms of the nodal point coordinates.

Given the spatial nodal coordinates (x i , y i ), we can use the shape functions (Ni ) to map

any spatial point (x, y) in the element by using following equations

4
X
x= Ni x i , (33)
i =1

4
X
y= Ni y i . (34)
i =1

Using equation 33 we can write

∂x ∂N1 ∂N2 ∂N3 ∂N4

= x1 + x2 + x3 + x4 , (35)
∂ξ ∂ξ ∂ξ ∂ξ ∂ξ

∂x ∂N1 ∂N2 ∂N3 ∂N4

= x1 + x2 + x3 + x4 . (36)
∂η ∂η ∂η ∂η ∂η

Similarly, by using equation 34 we have

∂y ∂N1 ∂N2 ∂N3 ∂N4

= y1 + y2 + y3 + y4 , (37)
∂ξ ∂ξ ∂ξ ∂ξ ∂ξ

∂y ∂N1 ∂N2 ∂N3 ∂N4

= y1 + y2 + y3 + y4 . (38)
∂η ∂η ∂η ∂η ∂η

Page 17 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

Writing equations 35 - 38 in a matrix form we get

 
 x1 y1 
   
∂x ∂y ∂N1 ∂N2 ∂N3 ∂N4  

 ∂ξ ∂ξ  
  ∂ξ ∂ξ ∂ξ ∂ξ 



     x2 y2 
J = = . (39)
    

    
 x y3 
 3
  
∂x ∂y ∂N1 ∂N2 ∂N3 ∂N4
   
∂η ∂η ∂η ∂η ∂η ∂η
. 


x4 y4

Using equation 39 we can calculate J from the known quantities on the right hand side.

The above equation can be written more compactly as

J = dNT X , (40)

where d N T and X are the first and second matrices on the right hand side of equation

39.

1.4 Calculation of stress and strain at integration points

If the nodal point displacements are given then the corresponding displacement field values

within the element can be obtained using the shape functions. The displacement values at

the integration points of the elements allow us to calculate strain and stress values at those

integration points and thereby leading to calculation of stiffness matrix through numerical

integration.

Page 18 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

The displacement fields are given as

n
X
w= Ni w i , (41)
i =1

n
θx = Ni θxi ,
X
(42)
i =1

n
θy = Ni θ yi .
X
(43)
i =1

Using equations 41 to 43 in equation 8, the bending strain vector ²b is computed as

 
∂ Pn
∂x i =1 Ni (ξ, η)θ yi
 
 
 
 
 
 
 
²b = z 
 
∂ Pn . (44)


− ∂y i =1 Ni (ξ, η)θxi 

 
 
 
 
 
∂ ∂
 Pn Pn 
− ∂x i =1 Ni (ξ, η)θxi + ∂y i =1 Ni (ξ, η)θ yi

Page 19 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

By expanding the terms in equation 44 we obtain

 
∂N1 θ y1 ∂N2 θ y2 ∂N3 θ y3 ∂N4 θ y4
 ∂x
+ ∂x
+ ∂x
+ ∂x

 
 
 
 
 
 
²b = z  − ∂N∂y
1 θx1
− ∂N∂y
2 θx2
− ∂N∂y
3 θx3
− ∂N∂y
4 θx4
 
. (45)
 
 
 
 
 
 
 
∂N1 θ y1 ∂N2 θ y2 ∂N3 θ y3 ∂N4 θ y4
− ∂N∂x
1 θx1
− ∂N∂x
2 θx2
− ∂N∂x
3 θx3
− ∂N∂x
4 θx4
 
+ ∂y
+ ∂y
+ ∂y
+ ∂y

We can split the matrix on the right hand side of equation 45 as product of a matrix B b

Page 20 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

and nodal point displacement vector d as follows

 
 w1 
 
 
θx1 
 

 
 
 

 θ y1 

 
 
w2 
  

∂N1 ∂N2 ∂N3 ∂N4
 0 0 0 0 0 0 0 0
 
∂x ∂x ∂x ∂x  
  
θx2 
 
 
  
  
  
  θ y2 
²b = z  ∂N
0 − ∂N 0 − ∂N 0 − ∂N
  
3 , (46)
 0 − ∂y1 0 2
0 0 4
0 

∂y ∂y ∂y 
w3 
  
 
  
  
θx3 
  
 
0 − ∂N ∂N1
0 − ∂N ∂N2
0 − ∂N ∂N3
0 − ∂N ∂N4
  
1 2 3 4  
∂x ∂y ∂x ∂y ∂x ∂y ∂x ∂y  

 θ y3 

 
 
 

 w4  
 
 

 θx4 
 
 
θ y4

i.e.,

²b = z B b d . (47)

∂Ni ∂Ni
The spatial derivatives of the shape functions ∂x and ∂y , in equation 46, can be obtained

from equation 31. A similar procedure can be adopted to compute the components of the

shear strain vector (²s ).

Page 21 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

By substituting w, θx , θ y from equations 41 - 43 in equation 9 we obtain

   
 γxz
Pn ∂ Pn
i =1 Ni θ yi + ∂x i =1 Ni w i
  
²s =  = . (48)
   
Pn ∂ Pn
γy z i =1 Ni θxi
   
− + ∂y i =1 Ni w i

Upon expanding the terms in equation 48

 
 ∂N1 w 1
∂x
+ N1 θ y1 + ∂N∂x2 w2
+ N2 θ y2 + ∂N∂x3 w3
+ N3 θ y3 + ∂N∂x4 w4
+ N4 θ y4 
²s =  . (49)
 
∂N1 w 1 ∂N2 w 2 ∂N3 w 3 ∂N4 w 4
θ θ θ θ
 
∂y − N 1 x1 + ∂y − N 2 x2 + ∂y − N 3 x3 + ∂y − N 4 x4

The matrix on the right hand side in equation 49 can again be written as product of a

Page 22 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

matrix B s and nodal displacement vectors d as

 
 w1 
 
 
θx1 
 

 
 
 

 θ y1 

 
 
w2 
 

 
 
 θx2 
 

∂N1 ∂N2 ∂N3 ∂N4
 
 ∂x 0 N1 ∂x 0 N2 ∂x 0 N3 ∂x 0 N4   
 
θ y2 

 
²s = 
  
 , (50)

  
 w3 
 

∂N1 ∂N2 ∂N3 ∂N4
  
∂y
−N1 0 ∂y
−N2 0 ∂y
−N3 0 ∂y
−N4 0  
θx3 
 

 
 
 

 θ y3 

 
 
 

 w4  
 
 

 θx4 
 
 
θ y4

i.e.,

²s = B s d . (51)

∂Ni ∂Ni
As before ∂x and ∂y in equation 50 can be obtained from equation 31.

Page 23 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

1.5 Principle of virtual work and derivation of element stiffness matrix

Since we have considered small elastic deformations only, therefore both K b and K s (stiffness

matrices corresponding to bending and shearing response, respectively) are constant. Non-

linear material behavior or finite deformation may require updating of the element stiffness

matrix during iterative solution procedure of the global equilibrium equations.

At each node we have considered a total three degrees of freedom, one translational

(w in z direction) and two rotational (θx and θ y about X and Y directions, respectively).

Corresponding to these degrees of freedom, at every node we have a work conjugate force

F z and moments M x and M y . Thus, the nodal force vector of an element consists of 12

Page 24 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

components and we represent it as

 
 F z1 
 
 
M x1 
 

 
 
 

 M y1 

 
 
F z2 
 

 
 
 

 M x2 

 
 

 M y2 

F=

.
 (52)
F z3 
 

 
 
 

 M x3 

 
 

 M y3 

 
 
 

 F z4 

 
 

 M x4 

 
 
M y4

To derive the expressions for bending stiffness and shear stiffness we decompose the

internal work done into bending and shear components. Let us consider an element which

has undergone a certain level of bending (²b ) and shear deformation (²s ). In this state, the

associated stress vectors are denoted by σb and σs , the nodal point displacement vector is

represented by U and the external (equilibrating) nodal force vector by F . If we imagine a

virtual displacement of δU , then according to the principle of virtual work, external virtual

Page 25 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

work being equal to the internal virtual work, we can write

Z Z
>
F δU = σ>
b δ²b d v + σ>
s δ²s d v. (53)
V V

Since ²b = zB b U (using equation 47) and ²s = zB s U (using equation 51), therefore

δ²b = zB b δU (54)

and,

δ²s = B s δU . (55)

Using σb (= D b ²b ) = zD b B b U , σs (= D s ²s ) = D s B s U and equations 54 and 55 on the right

hand side of equation 53 we obtain

Z Z
>
F δUb = U >
B b> D b> B b δU z 2 d v + U > B s> D s> B s δU d v. (56)
V V

Since above equation holds good for all virtual displacements δU , thus we have

Z Z
>
F = U >
B b> D b> B b z 2 d v + U > B s> D s> B s δU d v. (57)
V V

By taking transpose of both sides in the above equation, and using F = K b U + K s U we can

Page 26 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

arrive at:

Z
Kb = B b> D b B b z 2 d v (58)
V

and

Z
Ks = B s> D s B s d v. (59)
V

Now we shall analyze the computation of both K b and K s . In equation 58 the quantities

within the volume integral are not a function of thickness, except the ‘z term’ itself, and hence

it can be integrated across the thickness from −h/2 to h/2 to yield a result in terms of area

integral over the element as

h3
Z
Kb = B b> D b B b d Ω. (60)
12 Ω

In finite element analysis the integration of the stiffness matrix is done numerically, and

since we are using isoparametric elements, we can re-write the integral on the right hand

side of the equation 60 in terms of isoparametric coordinates as follows

h3 1 1
Z Z
Kb = B b (ξ, η)> D b B b (ξ, η) |J (ξ, η)| d ξd η. (61)
12 −1 −1

Here, B b and J depend on the integration point coordinates (ξ, η). We shall use 2x2

−1 −1
quad integration rule with four integration points with coordinates of (ξ1 = p ,η = p ),
3 1 3

(ξ2 = p1 , η 2 = −1
p ), (ξ3 = p1 , η 3 = p1 ) and (ξ4 = −1
p ,η = p1 ). The right hand side of the
3 3 3 3 3 4 3

Page 27 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

equation 61 can be numerically integrated as follows

h3 X 4
Kb = B b (ξi , η i )> D b B b (ξi , η i ) |J (ξi , η i )| ∆ξ∆η, (62)
12 i =1

where ∆ξ = 1 and ∆η = 1. For the sake of simplicity we express the above equation as

4
IPi
X
Kb = (63)
i =1

where,

h3
IPi = B b (ξi , η i )> D b B b (ξi , η i ) |J (ξi , η i )| ∆ξ∆η. (64)
12

Next, the integral on the right hand side of equation 59 for shear stiffness K s can be

written as

Z Z h/2
Ks = B s> D s B s d z d Ω. (65)
Ω −h/2

As all the quantities inside the integral on the right hand side of the above equation are

independent of the thickness coordinate, therefore we can integrate along that direction to

arrive at

Z
Ks = h B s> D s B s d Ω. (66)
Ω

Page 28 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

We can re-write right hand side of equation 66 in terms of isoparametric coordinates as

Z 1 Z 1
Ks = h B s (ξ, η)> D s B s (ξ, η) |J (ξ, η)| d ξd η. (67)
−1 −1

In order avoid shear locking, we perform reduced numerical integration for the shear

stiffness by choosing a single integration point whose coordinates are given by ξ = 0, η = 0.

Thus, the expression for K s can be written as

K s = h B s (ξ = 0, η = 0)> D s B s (ξ = 0, η = 0) |J (ξ = 0, η = 0)| ∆ξ∆η (68)

where ∆ξ = 2 and ∆η = 2. Finally, overall stiffness K of the element in terms of bending

stiffness K b and shear stiffness K s can be written as

K = Kb + K s . (69)

2 S AMPLE H AND C ALCULATIONS

In this section we demonstrate through hand calculations how K b and K s are calculated for

an element. Since we are working under the assumption of small strain elasticity, the stiffness

matrix for each element will be the same and constant. E , ν and h are chosen as 2 × 108 , 0.3

and 0.1, respectively throughout this study. The shear correction factor is taken to be 56 .

Page 29 of 97
 Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

2.1 Bending Stiffness

The expression for K b was derived in equation 62 based on four integration points. The

isoparametric co-ordinates (ξi ,η i ) for these four integration points are ( p , −1 ), ( p1 , p

−1 p −1
),
3 3 3 3

( p1 , p1 ) and ( p
−1 p1
, ). We will show the calculation of I P 1 (i = 1 in equation 64) for the first
3 3 3 3

−1 −1
integration point with ξ = p and η = p . Substituting these values in equations 25 and 26 we
3 3

obtain

1 1 1 1
- , ,
3 3
3 3
η
Node 3
(-1, 1) (1, 1)
(x3, y3)
x x
ζ
x x
x
(-1, -1) (1, -1)

Node 2 -
1
3
,-
1
3
1
3
,-
1
3
(x2, y2)
Figure 4: Four-noded element with integration points in isoparametric formulation.

 
 0.622008 
 
 
 
 0.166666 
N = , (70)
 
 
 0.044658 
 
 
 
0.166666

Page 30 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

and

 
 −0.39433 −0.39433 
 
 
 
 0.39433 −0.10566 
dN =  . (71)
 
 
 0.10566 0.10566 
 
 
 
−0.10566 0.39433

If we choose the spatial nodal coordinates of the element nodes in X-Y coordinates as

(0, 0), (1, 0), (1, 1) and (0, 1) then we can construct X as

 
 0 0 
 
 
 
 1 0 
X = . (72)
 
 
 1 1 
 
 
 
0 1

The Jacobian J according to equation 39 is computed as

 
 0 0 
  
 
 
 −0.39433 0.39433 0.10566 −0.10566   1 0 
J = (73)
  
 
  
−0.39433 −0.10566 0.10566 0.39433 
 1 1 

 
 
0 1

Page 31 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

i.e.,

 
 0.5 0 
J = . (74)
 
 
0 0.5

J −1 is computed by inverting J and turns out to be

 
 2 0 
J −1 =  . (75)
 
 
0 2

Spatial derivatives of the all shape functions are denoted in a matrix X Y 0 and with aid of

equation 31 as

 
∂N1 ∂N1
 ∂x ∂y 
 
 
 
 
 
 
∂N2 ∂N2
 
 
 ∂x ∂y 
 
0  = d N J −> .
 
XY =
  (76)
 
 
 ∂N3 ∂N3 
∂x ∂y
 
 
 
 
 
 
 
 
∂N4 ∂N4
∂x ∂y

The matrix of derivatives with respect to spatial coordinates (X Y 0 ) can be calculated

Page 32 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

based on equation 76 as

 
 −0.78867 −0.78867 
 
 
 
 0.78867 −0.211325 
XY 0 = . (77)
 
 
 0.211325 0.211325 
 
 
 
−0.211325 0.78867

Now we can assemble the elements of B b matrix (as defined by equations 46 and 47) using

the elements of X Y 0 matrix obtained in equation 77,

 
 0 0 −.78867 0 0 .78867 0 0 .211325 0 0 −.211325 
 
 
 
 
 
 
 
Bb = 
 0 .78867 0 0 .211325 0 0 −.21135 0 0 −.78867 0
.

 
 
 
 
 
 
 
0 .78867 −.78867 0 −.78867 −.211325 0 −.211325 .21135 0 .211325 .78867
(78)

Now I P 1 (as defined by equation 64) can be calculated with all the known quantities.

Similar procedure can be carried to compute I P 2 , I P 3 and I P 4 , and finally K b is calculated

according to equation 62. It is worth mentioning that B b and J will vary as a function of (ξ, η)

Page 33 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

and will be different for each integration point. We get K b as

 
 0 0 0 0 0 0 0 0 0 0 0 0 
 
 

 0 0.824 −0.298 0 0.091 −0.022 0 −0.412 0.298 0 −0.504 0.022 

 
 
0 −0.298 0.824 0 0.022 −0.504 0 0.298 −0.412 0 −0.022 0.091
 
 
 
 
 

 0 0 0 0 0 0 0 0 0 0 0 0 

 
 

 0 0.091 0.022 0 0.824 0.298 0 −0.504 −0.022 0 −0.412 −0.298 

 
 
 0 −0.022 −0.504 0 0.298 0.824 0 0.022 0.091 0 −0.298 −0.412 
 × 104
 
Kb = 
  (79)
0 0 0 0 0 0 0 0 0 0 0 0
 
 
 
 
 

 0 −0.412 0.298 0 −0.504 0.022 0 0.824 −0.298 0 0.091 −0.022 

 
 

 0 0.298 −0.412 0 −0.022 0.091 0 −0.298 0.824 0 0.022 −0.504 

 
 

 0 0 0 0 0 0 0 0 0 0 0 0 

 
 
0 −0.504 −0.022 0 −0.412 −0.298 0 0.091 0.022 0 0.824 0.298
 
 
 
 
0 0.022 0.091 0 −0.298 −0.412 0 −0.022 −0.504 0 0.298 0.824

2.2 Shear Stiffness

For calculation of stiffness associated with shear deformation, we use the reduced integration

approach, i.e., only one integration point with ξ = 0 and η = 0 is chosen. According to

equations 25 and 26 N and d N are computed as:

 
 0.25 
 
 
 
 0.25 
N = , (80)
 
 
 0.25 
 
 
 
0.25

 
 −0.25 −0.25 
 
 
 
 0.25 −0.25 
dN =  . (81)
 
 
 0.25 0.25 
 
 
 
−0.25 0.25

Matrix of nodal coordinates X is chosen same as in equation 72, and the Jacobian J is

Page 34 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

calculated using the equation 39 as

 
 0 0 
  
 
 
 −0.25 0.25 0.25 −0.25   1 0 
J = , (82)
  

  
−0.25 −0.25 0.25 0.25  1 1 
 
 
 
0 1

 
 0.5 0 
J = . (83)
 
 
0 0.5

J −1 turns out to be

 
 2 0 
J −1 =  . (84)
 
 
0 2

The matrix of X and Y derivatives of shape functions (X Y 0 ) is again calculated according

to equation 76 as

 
 −0.5 −0.5 
 
 
 
 0.5 −0.5 
XY 0 = . (85)
 
 
 0.5 0.5 
 
 
 
−0.5 0.5

Using the components of X Y 0 from equation 85, we can now construct the B s matrix

Page 35 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

(based on equations 50 and 51) as

 
 −0.5 0 0.25 0.5 0 0.25 0.5 0 0.25 −0.5 0 0.25 
Bs =  .
 
 
−0.5 −0.25 0 −0.5 −0.25 0 0.5 −0.25 0 0.5 −0.25 0

Finally, the shear stiffness matrix K s can be calculated from equation 68 and we get

 
 0.385 0.096 −0.096 0 0.096 −0.096 −0.385 0.096 −0.096 0 0.096 −0.096 
 
 

 0.096 0.048 0 0.096 0.048 0 −0.096 0.048 0 −0.096 0.048 0 

 
 
−0.096 0 0.048 0.096 0 0.048 0.096 0 0.048 −0.096 0 0.048 
 

 
 
 

 0 0.096 0.096 0.385 0.096 0.096 0 0.096 0.096 −0.385 0.096 0.096 
 
 

 0.096 0.048 0 0.096 0.048 0 −0.096 0.048 0 −0.096 0.048 0 

 
 
 −0.096 0 0.048 0.096 0 0.048 0.096 0 0.048 −0.096 0 0.048 
 × 107 .
 
Ks = 
  (86)
−0.385 −0.096 0.096 0 −0.096 0.096 0.385 −0.096 0.096 0 −0.096 0.096 
 

 
 
 

 0.096 0.048 0 0.096 0.048 0 −0.096 0.048 0 −0.096 0.048 0 

 
 

 −0.096 0 0.048 0.096 0 0.048 0.096 0 0.048 −0.096 0 0.048 
 
 

 0 −0.096 −0.096 −0.385 −0.096 −0.096 0 −0.096 −0.096 0.385 −0.096 −0.096 

 
 
0.096 0.048 0 0.096 0.048 0 −0.096 0.048 0 −0.096 0.048 0
 
 
 
 
−0.096 0 0.048 0.096 0 0.048 0.096 0 0.048 −0.096 0 0.048

The total stiffness matrix for the element is then obtained by summing the bending and

shear stiffness response as

 
 0.385 0.096 −0.096 0 0.096 −0.096 −0.385 0.096 −0.096 0 0.096 −0.096 
 
 

 0.096 0.048 −0.298 0.096 0.048 −0.022 −0.096 0.047 0.298 −0.096 0.047 0.022 

 
 
−0.096 −0.298 0.048 0.096 0.022 0.047 0.096 0.298 0.047 −0.096 −0.022 0.048 
 

 
 
 

 0 0.096 0.096 0.385 0.096 0.096 0 0.096 0.096 −0.385 0.096 0.096 

 
 

 0.096 0.048 0.022 0.096 0.048 0.298 −0.096 0.047 −0.022 −0.096 0.047 −0.298 

 
 
 −0.096 −0.022 0.047 0.096 0.298 0.048 0.096 0.022 0.048 −0.096 −0.298 0.047 
 × 107 .
 
K =
  (87)
−0.385 −0.096 0.096 0 −0.096 0.096 0.385 −0.096 0.096 0 −0.096 0.096 
 

 
 
 

 0.096 0.047 0.298 0.096 0.047 0.022 −0.096 0.048 −0.298 −0.096 0.048 −0.022 

 
 

 −0.096 0.298 0.047 0.096 −0.022 0.048 0.096 −0.298 0.048 −0.096 0.022 0.047 
 
 

 0 −0.096 −0.096 −0.385 −0.096 −0.096 0 −0.096 −0.096 0.385 −0.096 −0.096 

 
 
0.096 0.047 −0.022 0.096 0.047 −0.298 −0.096 0.048 0.022 −0.096 0.048 0.298 
 

 
 
−0.096 0.022 0.048 0.096 −0.298 0.047 0.096 −0.022 0.047 −0.096 0.298 0.048

Page 36 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

3 S OME C ONSEQUENCES OF M INDLIN P LATE T HEORY

Here we discuss consequences of the assumptions made in the Mindlin theory and associated

inconsistencies that arise in the formulation.

First we assumed that transverse normal stress σz = 0, by construction, through the plate

thickness. On the other hand we allowed for non-zero transverse shear stresses through

thickness. To investigate this let us study the equilibrium equation in the Z direction:

∂σzx ∂σz y ∂σz

+ + = 0, (88)
∂x ∂y ∂z

substituting σz = 0 in the above equation we obtain

∂σzx ∂σz y
+ = 0. (89)
∂x ∂y

Substitute σzx = Gγzx and σz y = Gγz y in equation 89, where G is the shear modulus, to

arrive at

∂γzx ∂γz y
+ = 0. (90)
∂x ∂y

From the previous kinematic assumptions made we have

∂w
γzx = + θy ,
∂x

and

Page 37 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

∂w
γz y = − θx .
∂y

Substituting γzx and γz y in equation 90 we obtain

∂ ∂w ∂ ∂w
· ¸ · ¸
+ θy + − θx = 0 (91)
∂x ∂x ∂y ∂y

∂2 w ∂θ y ∂2 w ∂θx
+ + − = 0. (92)
∂x 2 ∂x ∂y 2 ∂y

∂2 w
For small deformation, we can approximate ∂x 2
= K x (curvature in X-direction) and

∂2 w
∂y 2
= K y (curvature in Y-direction), and therefore the equation 92 becomes

∂θx ∂θ y
Kx + K y = − . (93)
∂y ∂x

In the Mindlin theory we have (from equations 2 and 3)

∂θx ²y ∂θ y ²x
=− , = . (94)
∂y z ∂x z

Thus using above expressions in equation 93 we obtain

z(K x + K y ) = −(²x + ² y ). (95)

The above condition may not be exactly satisfied in general throughout the thickness and

is a consequence of ad-hoc assumptions made on the kinematics of deformation and stress

Page 38 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

state in the plate.

As mentioned in the beginning, the Mindlin theory does not take into account any mem-

brane strains, i.e. if we assume a case of homogeneous biaxial straining, with ²x = ² y = δ > 0,

then at z = 0 equation 95 is not satisfied. Lastly, we had assumed σz = 0, but this in turn

implies ²z 6= 0 which again has not been accounted exactly.

Page 39 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

4 ABAQUS UEL I MPLEMENTATION

UEL is a programming interface provided with Abaqus Standard using which one can define

customized finite elements. In this section, we discuss several implementation details in

regards to the UEL. In order to work with UEL, one must start by using the default UEL

header from the Abaqus User Subroutines Reference Guide, Section 1.1.28 in Abaqus 2016

documentation, as shown in Figure 5. The user specified program is contained under this

header and the end of the UEL is marked by the ‘END’ syntax.

Figure 5: UEL subroutine header from Abaqus 2016 documentation, Abaqus User Subroutines Refer-
ence Guide, Section 1.1.28.

Conceptually it is important to understand the overall operation of the ‘UEL’ block. Figure

6 shows a schematic describing the ‘input’ and ‘output’ flow from the UEL subroutine, that is

Page 40 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

COORDS RHS
Unew
dUold UEL AMATRX
State variables Updated state variables

Figure 6: Input-Output block diagram for the UEL subroutine.

called for every element. For a current time increment, Abaqus provides the incremental and

final nodal point displacements. State variables at the start of that time increment are also

provided. UEL is then required to return the updated element stiffness matrix, internal force

and state variables all at the end of the current time step. Next we discuss in greater detail,

the incoming quantifies provided by Abaqus into the UEL and associated outgoing quantities

provided back to Abaqus.

4.1 UEL input variables

U : is a vector of size number-of-element-nodes×degrees-of-freedom-per-node, and contains

the degree of freedom values at each node at the end of the current time increment. This is

shown as Unew in Figure 6.

DU : is a vector of size number-of-nodes×degrees-of-freedom-per-node and contains the

incremental degree of freedom values for all the nodes in the previous increment. This is

dUol d in the Figure 6.

COORDS: is a matrix of size-of-spatial-dimensions×number-of-nodes and contains the

Page 41 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

coordinates of the nodes of an element, expressed in (x, y, z) format. A related parameter,

COORD I N AT E S needs to be defined in the UEL section of the Abaqus input file which sets

the value of ‘size-of-spatial-dimensions’.

4.2 UEL output variables

SV ARS: is a vector which stores the state variables of an element and needs to be passed

back to Abaqus after being updated. Typically stress and strains at all integration points

within an element are stored.

AM AT R X : is the stiffness matrix for an element. For our case, small deformation and linear

elasticity, the stiffness matrix is constant, however, if there are non-linear effects then the

stiffness matrix needs to be updated before returning to Abaqus.

R H S: is the internal force vector for the element which needs to be returned at the end of

the current increment. Since, we are solving a linear elasticity problem, for simplicity, the

internal force vector can be calculated as

R H S = −AM AT R X ×U (96)

4.3 Pseudo code for the UEL subroutine

The overall implementation of the UEL is summarized in form of a pseudo code as shown in

Algorithm 1. The pseudo code is based on the theory presented earlier and is self-explanatory.

Page 42 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

We also provide the actual UEL written in Fortran corresponding to this pseudo code subse-

quently.

Page 43 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

Algorithm 1 Pseudo code for UEL.

1: procedure UEL
2: U ← Abaqus . D.O.F. values at current increment
3: DU ← Abaqus . Previous incremental D.O.F. values
4: SV ARS ← Abaqus . Previous state variables for the element
5: COORDS ← Abaqus . Current nodal coordinates of the element
6: X ← COORDS . Converted 2D nodal coordinates of the element
7: Isoparametric formulation:
8: Initialize I P _COORDS . Initialize integration point coordinates
9: Db ← E , ν
10: i =1 . First integration point
11: Kb = 0 . Initialize bending stiffness matrix to zero
12: loop : . Loop over 4 integration point
13: (ξi , η i ) ← I P _COORDS . Get (ξ, η) for the current integration point
14: N ← (ξi , η i ) . Calculate shape functions
15: d N ← (ξi , η i ) . Calculate derivatives of shape functions
16: J ← d N >X . Calculate Jacobian
17: J −1 . Calculate inverse of Jacobian
18: X Y 0 ← d N J −> . Calculate XY derivatives of shape function
19: Bb ← X Y 0
. Assemble B-matrix for bending
i h3
20: >
I P ← 12 × B b D b B b × |J | . Stiffness contribution for one integration point
21: if i == 4 then
22: break . Exit the loop once all four integration points are completed
23: else
24: Kb ← Kb + I P i . Form bending stiffness matrix
25: i = i +1
26: goto loop . Iterate for next integration point
27: close;
28: Ds ← E , ν
29: Ks = 0 . Initialize shear stiffness matrix to zero
30: (ξ, η) ← (0, 0) . Reduced integration for shear part
31: N ← (ξ, η) . Shape function values for reduced integration at origin
32: d N ← (ξ, η) . Derivatives of shape function
33: J ← d N >X . Calculate Jacobian
34: J −1
. Calculate inverse of Jacobian
35: X Y 0 ← d N J −> . Calculate XY derivates of shape function
36: Bs ← X Y 0 . Assemble B-matrix for shear
37: >
K s ← h × B s D s B s × |J | × 4 . Calculate shear stiffness matrix
38: AM AT R X ← K b + K s . Total stiffness
39: R H S ← −AM AT R X U . Calculate internal force vector
40:
41: SV ARS → To Abaqus . Updated state variables for the element
42: AM AT R X → To Abaqus . Stiffness matrix for the element
43: R H S → To Abaqus . Internal force vector
Page 44 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

1 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
2 CCCCCCC THIS ROUTINE IS WRITTEN BY SUBODH KALIA, NIKHIL PADHYE CCCCCCCCC
3 CCCCCCC FOR ABAQUS 2016. THE DISTRIBUTION OF THIS CODE WITHOUT PRIOR CCC
4 CCCCCCC PERMISSION OR ANY KIND OF COMMERICAL USE IS STRICTLY CCCCCCCCCCC
5 CCCCCCC PROHIBITED. PLEASE CITE THE WORK IF YOU WANT TO USE IT FOR CCCCC
6 CCCCCCC YOUR RESEARCH. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
7 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
8 SUBROUTINE UEL(RHS,AMATRX,SVARS,ENERGY,NDOFEL,NRHS,NSVARS,
9 1 PROPS,NPROPS,COORDS,MCRD,NNODE,U,DU,V,A,JTYPE,TIME,DTIME,
10 2 KSTEP,KINC,JELEM,PARAMS,NDLOAD,JDLTYP,ADLMAG,PREDEF,NPREDF,
11 3 LFLAGS,MLVARX,DDLMAG,MDLOAD,PNEWDT,JPROPS,NJPROP,PERIOD)
12 C
13 INCLUDE 'ABA_PARAM.INC'
14 C
15 DIMENSION RHS(MLVARX,*),AMATRX(NDOFEL,NDOFEL),PROPS(*),
16 1 SVARS(*),ENERGY(8),COORDS(MCRD,NNODE),U(NDOFEL),
17 2 DU(MLVARX,*),V(NDOFEL),A(NDOFEL),TIME(2),PARAMS(*),
18 3 JDLTYP(MDLOAD,*),ADLMAG(MDLOAD,*),DDLMAG(MDLOAD,*),
19 4 PREDEF(2,NPREDF,NNODE),LFLAGS(*),JPROPS(*)
20
21 Double precision, Dimension(12, 1) :: FORCE, V_TEMP
22 Double precision, Dimension(12, 1) :: U_CONV, DU_CONV
23 Double precision, Dimension(4, 1) :: N
24 Double precision, Dimension(4, 2) :: DN, IP_COORDS, COORDS_T
25 Double precision, Dimension(4, 2) :: XY_DERIVATIVES
26 Double precision, Dimension(2, 4) :: XY_DERIVATIVES_T
27 Double precision, Dimension(2, 4) :: DN_TRANSPOSE
28 Double precision, Dimension(2, 2) :: J, INV_J
29 Double precision, Dimension(3, 12) :: M_B_MAT_BENDING, M_TEMP
30 Double precision, Dimension(12, 12) :: M_STIFFNESS, M_TEMP_STIFF
31 Double precision, Dimension(12, 12) :: M_STIFF
32 Double precision, Dimension(12, 3) :: M_B_MAT_BENDING_TRANSPOSE
33 Double precision, Dimension(3, 1) :: V_STRAIN_INCREMENT_BENDING
34 Double precision, Dimension(3, 1) :: V_STRESS_INCREMENT_BENDING
35 Double precision, Dimension(3, 3) :: M_D_BENDING
36 Double precision, Dimension(2, 12) :: M_B_MAT_SHEAR
37 Double precision, Dimension(2, 1) :: V_STRAIN_INCREMENT_SHEAR
38 Double precision, Dimension(2, 1) :: V_STRESS_INCREMENT_SHEAR
39 Double precision, Dimension(12, 2) :: M_B_MAT_SHEAR_TRANSPOSE
40 Double precision, Dimension(2, 2) :: M_D_SHEAR
41
42 Double precision :: CHI, ETA, S_E, S_NU, S_INERTIA, S_THICKNESS
43 Double precision :: S_J_DETERMINANT, S_KAPPA
44
45 INTEGER :: NDOFEL, NSVARS, NPROPS, MCRD, NNODE, K1, K2, K3, NUM_IP
46 INTEGER :: NDOF, NUM_SVAR_PER_IP, POSITION_POINTER
47
48 PARAMETER(NUM_IP=4, NUM_SVAR_PER_IP=10)
49
50 COORDS_T(1,1) = COORDS(1,1)
51 COORDS_T(1,2) = COORDS(2,1)
52
53 COORDS_T(2,1) = COORDS(1,2)
54 COORDS_T(2,2) = COORDS(2,2)
55
56 COORDS_T(3,1) = COORDS(1,3)
57 COORDS_T(3,2) = COORDS(2,3)
58
59 COORDS_T(4,1) = COORDS(1,4)
60 COORDS_T(4,2) = COORDS(2,4)
61
62 DO K1 = 1, 12
63 U_CONV(K1, 1) = U(K1)
64 DU_CONV(K1, 1) = DU(K1,1)
65 END DO
66
67 IF (LFLAGS(3).EQ.4) THEN
68 DO K1=1, 12
69 DO K2=1, 12

Page 45 of 97
 Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

70 AMATRX(K1, K2) = 0.D0

71 END DO
72 AMATRX(K1, K1) = 1.D0
73 END DO
74 GOTO 999
75 END IF
76
77 CCCCCCCCCCCCCCCCCCCCCCCCCC INITIALIZE PARAMETERS CCCCCCCCCCCCCCCCCCCCCCC
78 DO K1 = 1, 12
79 DO K2 = 1, 12
80 M_STIFFNESS(K1, K2) = 0.D0
81 M_STIFF(K1, K2) = 0.D0
82 END DO
83 END DO
84
85 DO K1 = 1, 12
86 FORCE(K1, 1) = 0.D0
87 V_TEMP(K1, 1) = 0.D0
88 END DO
89
90 S_E = 200000000
91 S_NU = 0.3D0
92 S_G = S_E/(2.D0*(1.D0 + S_NU))
93 S_KAPPA = 6.D0/6.D0
94 S_THICKNESS = 0.1D0
95 S_INERTIA = (S_THICKNESS**3.D0)/12.D0
96
97 M_D_BENDING(1, 1) = 1.D0
98 M_D_BENDING(1, 2) = S_NU
99 M_D_BENDING(1, 3) = 0.D0
100 M_D_BENDING(2, 1) = S_NU
101 M_D_BENDING(2, 2) = 1.D0
102 M_D_BENDING(2, 3) = 0.D0
103 M_D_BENDING(3, 1) = 0.D0
104 M_D_BENDING(3, 2) = 0.D0
105 M_D_BENDING(3, 3) = (1.D0 - S_NU)/2.D0
106
107 M_D_BENDING = (S_E/(1.D0 - (S_NU**2.D0)))*M_D_BENDING
108
109 M_D_SHEAR(1, 1) = S_G
110 M_D_SHEAR(1, 2) = 0.D0
111 M_D_SHEAR(2, 1) = 0.D0
112 M_D_SHEAR(2, 2) = S_G
113
114 CALL TWO_BY_TWO_INTEGRATION_POINT_COORDINATES(IP_COORDS, NUM_IP)
115
116 DO K1 = 1, NUM_IP
117 CHI = IP_COORDS(K1, 1)
118 ETA = IP_COORDS(K1, 2)
119 CALL FOUR_NODED_BILINEAR_SHAPE_FUNCTIONS(CHI, ETA, N, DN)
120
121 CALL GENERIC_TRANSPOSE_OF_MATRIX(DN, 4, 2, DN_TRANSPOSE)
122
123 CALL GENERIC_MATRIX_MULTIPLY(DN_TRANSPOSE,2,4,COORDS_T,4,2,J)
124
125 CALL MAT_INVERSE_TWO_BY_TWO(J, INV_J, S_J_DETERMINANT)
126
127 CALL GENERIC_MATRIX_MULTIPLY(INV_J,2,2,DN_TRANSPOSE,
128 1 2,4,XY_DERIVATIVES_T)
129
130 CALL GENERIC_TRANSPOSE_OF_MATRIX(XY_DERIVATIVES_T,2,4,
131 1 XY_DERIVATIVES)
132
133 CCCCCCCCCCCCCCCCCCCCCC ASSEMBLE B_BENDING (3X12) CCCCCCCCCCCCCCCCCCCCCCC
134 DO K2 = 1, 3
135 DO K3 = 1, 12
136 M_B_MAT_BENDING(K2, K3) = 0.D0
137 END DO
138 END DO

Page 46 of 97
 Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

139 M_B_MAT_BENDING(1, 3) = XY_DERIVATIVES(1, 1)

140 M_B_MAT_BENDING(1, 6) = XY_DERIVATIVES(2, 1)
141 M_B_MAT_BENDING(1, 9) = XY_DERIVATIVES(3, 1)
142 M_B_MAT_BENDING(1, 12) = XY_DERIVATIVES(4, 1)
143
144 M_B_MAT_BENDING(2, 2) = -XY_DERIVATIVES(1, 2)
145 M_B_MAT_BENDING(2, 5) = -XY_DERIVATIVES(2, 2)
146 M_B_MAT_BENDING(2, 8) = -XY_DERIVATIVES(3, 2)
147 M_B_MAT_BENDING(2, 11) = -XY_DERIVATIVES(4, 2)
148
149 M_B_MAT_BENDING(3, 2) = -XY_DERIVATIVES(1, 1)
150 M_B_MAT_BENDING(3, 5) = -XY_DERIVATIVES(2, 1)
151 M_B_MAT_BENDING(3, 8) = -XY_DERIVATIVES(3, 1)
152 M_B_MAT_BENDING(3, 11) = -XY_DERIVATIVES(4, 1)
153
154 M_B_MAT_BENDING(3, 3) = XY_DERIVATIVES(1, 2)
155 M_B_MAT_BENDING(3, 6) = XY_DERIVATIVES(2, 2)
156 M_B_MAT_BENDING(3, 9) = XY_DERIVATIVES(3, 2)
157 M_B_MAT_BENDING(3, 12) = XY_DERIVATIVES(4, 2)
158 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
159 CALL GENERIC_MATRIX_MULTIPLY(M_B_MAT_BENDING, 3, 12,
160 1 DU_CONV, 12, 1,
161 2 V_STRAIN_INCREMENT_BENDING)
162
163 CALL GENERIC_MATRIX_MULTIPLY(M_D_BENDING, 3, 3,
164 1 V_STRAIN_INCREMENT_BENDING, 3, 1,
165 2 V_STRESS_INCREMENT_BENDING)
166
167 CALL GENERIC_MATRIX_MULTIPLY(S_INERTIA*S_J_DETERMINANT*
168 1 M_D_BENDING, 3, 3, M_B_MAT_BENDING,
169 2 3, 12, M_TEMP)
170
171 CALL GENERIC_TRANSPOSE_OF_MATRIX(M_B_MAT_BENDING, 3, 12,
172 1 M_B_MAT_BENDING_TRANSPOSE)
173
174 CALL GENERIC_MATRIX_MULTIPLY(M_B_MAT_BENDING_TRANSPOSE, 12, 3,
175 1 M_TEMP, 3, 12, M_TEMP_STIFF)
176
177 CALL GENERIC_ADD_MATRIX(M_STIFFNESS, M_TEMP_STIFF, 12, 12,
178 1 M_STIFF)
179
180 M_STIFFNESS = M_STIFF
181
182 END DO
183 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
184 CHI = 0.D0
185 ETA = 0.D0
186 CALL FOUR_NODED_BILINEAR_SHAPE_FUNCTIONS(CHI, ETA, N, DN)
187 CALL GENERIC_TRANSPOSE_OF_MATRIX(DN, 4, 2, DN_TRANSPOSE)
188 CALL GENERIC_MATRIX_MULTIPLY(DN_TRANSPOSE,2,4,COORDS_T,4,2,J)
189 CALL MAT_INVERSE_TWO_BY_TWO(J, INV_J, S_J_DETERMINANT)
190 CALL GENERIC_MATRIX_MULTIPLY(INV_J,2,2,DN_TRANSPOSE,
191 1 2,4,XY_DERIVATIVES_T)
192
193 CALL GENERIC_TRANSPOSE_OF_MATRIX(XY_DERIVATIVES_T,2,4,
194 1 XY_DERIVATIVES)
195
196 CCCCCCCCCCCCCCCCCCCCCC ASSEMBLE B_SHEAR (2X12) CCCCCCCCCCCCCCCCCCCCCCC
197 DO K2 = 1, 2
198 DO K3 = 1, 12
199 M_B_MAT_SHEAR(K2, K3) = 0.D0
200 END DO
201 END DO
202 M_B_MAT_SHEAR(1, 1) = XY_DERIVATIVES(1, 1)
203 M_B_MAT_SHEAR(1, 4) = XY_DERIVATIVES(2, 1)
204 M_B_MAT_SHEAR(1, 7) = XY_DERIVATIVES(3, 1)
205 M_B_MAT_SHEAR(1, 10) = XY_DERIVATIVES(4, 1)
206
207 M_B_MAT_SHEAR(1, 3) = N(1, 1)

Page 47 of 97
 Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

208 M_B_MAT_SHEAR(1, 6) = N(2, 1)

209 M_B_MAT_SHEAR(1, 9) = N(3, 1)
210 M_B_MAT_SHEAR(1, 12) = N(4, 1)
211
212 M_B_MAT_SHEAR(2, 1) = XY_DERIVATIVES(1, 2)
213 M_B_MAT_SHEAR(2, 4) = XY_DERIVATIVES(2, 2)
214 M_B_MAT_SHEAR(2, 7) = XY_DERIVATIVES(3, 2)
215 M_B_MAT_SHEAR(2, 10) = XY_DERIVATIVES(4, 2)
216
217 M_B_MAT_SHEAR(2, 2) = -N(1, 1)
218 M_B_MAT_SHEAR(2, 5) = -N(2, 1)
219 M_B_MAT_SHEAR(2, 8) = -N(3, 1)
220 M_B_MAT_SHEAR(2, 11) = -N(4, 1)
221 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
222 CALL GENERIC_MATRIX_MULTIPLY(M_B_MAT_SHEAR, 2, 12,
223 1 DU_CONV, 12, 1,
224 2 V_STRAIN_INCREMENT_SHEAR)
225
226 CALL GENERIC_MATRIX_MULTIPLY(M_D_SHEAR, 2, 2,
227 1 V_STRAIN_INCREMENT_SHEAR, 2, 1,
228 2 V_STRESS_INCREMENT_SHEAR)
229
230 CALL GENERIC_TRANSPOSE_OF_MATRIX(M_B_MAT_SHEAR, 2, 12,
231 1 M_B_MAT_SHEAR_TRANSPOSE)
232
233
234 CALL GENERIC_MATRIX_MULTIPLY(S_KAPPA*S_J_DETERMINANT*
235 1 S_THICKNESS*4.D0*M_D_SHEAR, 2, 2, M_B_MAT_SHEAR,
236 2 2, 12, M_TEMP)
237
238 CALL GENERIC_MATRIX_MULTIPLY(M_B_MAT_SHEAR_TRANSPOSE, 12, 2,
239 1 M_TEMP, 2, 12, M_TEMP_STIFF)
240
241 CALL GENERIC_ADD_MATRIX(M_STIFFNESS, M_TEMP_STIFF, 12, 12,
242 1 M_STIFF)
243
244 M_STIFFNESS = M_STIFF
245
246 CCCCCCCCCCCCCCCCCCCCCC UPDATE STRESS, STRAIN STATE VARIABLES CCCCCCCCC
247 DO K1 = 1, NUM_IP
248
249 POSITION_POINTER = (K1 - 1)*NUM_SVAR_PER_IP + 1
250 C SIGMA_X
251 SVARS(POSITION_POINTER) = SVARS(POSITION_POINTER)+
252 1 V_STRESS_INCREMENT_BENDING(1, 1)
253 C SIGMA_Y
254 SVARS(POSITION_POINTER+1) = SVARS(POSITION_POINTER+1)+
255 1 V_STRESS_INCREMENT_BENDING(2, 1)
256 C SIGMA_XY
257 SVARS(POSITION_POINTER+2) = SVARS(POSITION_POINTER+2)+
258 1 V_STRESS_INCREMENT_BENDING(3, 1)
259 C SIGMA_XZ
260 SVARS(POSITION_POINTER+3) = SVARS(POSITION_POINTER+3)+
261 1 V_STRESS_INCREMENT_SHEAR(1, 1)
262 C SIGMA_YZ
263 SVARS(POSITION_POINTER+4) = SVARS(POSITION_POINTER+4)+
264 1 V_STRESS_INCREMENT_SHEAR(2, 1)
265 C EPSILON_X
266 SVARS(POSITION_POINTER+5) = SVARS(POSITION_POINTER+5)+
267 1 V_STRAIN_INCREMENT_BENDING(1, 1)
268 C EPSILON_Y
269 SVARS(POSITION_POINTER+6) = SVARS(POSITION_POINTER+6)+
270 1 V_STRAIN_INCREMENT_BENDING(2, 1)
271 C EPSILON_XY
272 SVARS(POSITION_POINTER+7) = SVARS(POSITION_POINTER+7)+
273 1 V_STRAIN_INCREMENT_BENDING(3, 1)
274 C EPSILON_XZ
275 SVARS(POSITION_POINTER+8) = SVARS(POSITION_POINTER+8)+
276 1 V_STRAIN_INCREMENT_SHEAR(1, 1)

Page 48 of 97
 Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

277 C EPSILON_YZ
278 SVARS(POSITION_POINTER+9) = SVARS(POSITION_POINTER+9)+
279 1 V_STRAIN_INCREMENT_SHEAR(2, 1)
280
281 END DO
282
283 CCCCCCCCCCCCCCCCCCCCCCCC CALCULATE FORCE CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
284 CALL GENERIC_MATRIX_MULTIPLY(M_STIFFNESS, 12, 12,
285 1 U_CONV, 12, 1, FORCE)
286
287 CCCCCCCCCCCCCCCCCCCCCCCCCCC RETURN TO ABAQUS CCCCCCCCCCCCCCCCCCCCCCCCCCC
288 DO K1 = 1, 12
289 RHS(K1, 1) = -FORCE(K1, 1)
290 END DO
291
292 DO K1 = 1, 12
293 DO K2 = 1, 12
294 AMATRX(K1, K2) = M_STIFFNESS(K1, K2)
295 END DO
296 END DO
297
298 999 CONTINUE
299
300
301 RETURN
302 END

Page 49 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

5 B ENCHMARKING OF UEL I MPLEMENTATION FOR M INDLIN P LATE

E LEMENT

In this section, we compare the results of plate bending from the Abaqus with the default S4

element, Matlab implementation of Mindlin plate from [5], and our developed 4-noded UEL

subroutine.

5.1 Single element test

Figure 7: Boundary conditions for a single element test.

A single element of length and width of 2 units each and a thickness of 0.1 units is

considered. Boundary conditions for this element are shown in Figure 7, where we have

restricted w, θx , θ y at nodes 1, 2, and 3. At node 4, we restrict θx and θ y and apply a maximum

displacement of 0.1 in the positive Z-direction. Thus, this is a well defined displacement

Page 50 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

driven boundary value problem for which we seek to find the solution. It is important to note

that Abaqus S4 element has a total of 6 degrees of freedom at each node, therefore we simply

constrain the degrees of freedom other than w, θx and θ y at all nodes.

Figure 8, 9 and 10 shows the reaction forces (RF3) in Z-direction, and reaction moments

RM1 and RM2 in X- and Y-directions, respectively from Abaqus S4 element analysis.

Figure 8: Reaction forces in Z-direction with Abaqus S4 element.

Results from the referenced Matlab code, ABAQUS UEL implementation of Mindlin plate,

and Abaqus S4 element analysis were gathered and compared for the nodal point reaction

forces and moments and shown in Figures 11, 12 and 13. Our UEL implementation yields

results exactly similar to those predicted by the referenced Matlab code. In general there

is a good agreement with Abaqus S4 element analysis, however, there is deviation in the

reaction forces at nodes 2 and 3. A detailed investigation reveals that reaction moments for

the UEL implementation and Matlab code were same in magnitude at four nodes, whereas

Page 51 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

Figure 9: Reaction moments in X-direction with Abaqus S4 element.

Figure 10: Reaction moments in Y-direction with Abaqus S4 element.

Page 52 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

in case of Abaqus S4 element the reaction moments in both directions varied within the

element. A simple explanation for this observation stems from the fact that in the Mindlin

plate element (UEL or Matlab implementation) the transverse shear strains are assumed to

be constant through the thickness, and interpolated based on the nodal point rotations and

displacements. The bending strains are interpolated based on nodal point rotations only.

In the current case, we have set the rotations and displacements at all the nodes equal to

zero except for the imposed displacement in the vertical direction at node 4. This causes the

vanishing of the bending strains all together and linearly varying (constant through thickness)

transverse shear strains. Thus the reaction moments at the nodes in the Mindlin element

appear to simply maintain the global moment-equilibrium and have the same magnitude in

both X- and Y-directions, at all nodes, due to symmetric shape of the element. On the other

hand the S4 element of Abaqus is more sophisticated and picks up the variation of bending

and shear more accurately.

Exact matching of the results between our UEL implementation and referenced Matlab

code, and overall consistency with respect to the Abaqus S4 element (foregoing certain

deviations) indicates that our procedure is working well.

Page 53 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

0 10,000
Abaqus S4 element
Reaction force in z-direction at node

Reaction force in z-direction at node

-50,000 UEL 4-noded element 0
MATLAB 4-noded element
-100,000 -10,000 Abaqus S4 element
UEL 4-noded element
-150,000 -20,000 MATLAB 4-noded element
1

2
-200,000 -30,000

-250,000 -40,000

-300,000 -50,000

-350,000 -60,000
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1
Displacement Displacement

(a) (b)
10,000 400,000
Abaqus S4 element
Reaction force in z-direction at node

Reaction force in z-direction at node

0 350,000 UEL 4-noded element
MATLAB 4-noded element
Abaqus S4 element 300,000
-10,000
UEL 4-noded element 250,000
-20,000 MATLAB 4-noded element
200,000
3

4
-30,000
150,000
-40,000
100,000

-50,000 50,000

-60,000 0
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1
Displacement Displacement

(c) (d)

Figure 11: Plots for reaction forces in Z-direction at the nodes with respect to nodal displacement.

0 0
Abaqus S4 element Abaqus S4 element
Reaction moment in x-direction at
Reaction moment in x-direction at

-20,000 UEL 4-noded element -20,000 UEL 4-noded element

-40,000 MATLAB 4-noded element -40,000 MATLAB 4-noded element
-60,000
-60,000
-80,000
node 2
node 1

-80,000
-100,000
-100,000
-120,000
-120,000
-140,000
-140,000 -160,000
-160,000 -180,000
-180,000 -200,000
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1
Displacement Displacement

(a) (b)
0 0
Abaqus S4 element Abaqus S4 element
Reaction moment in x-direction at
Reaction moment in x-direction at

-20,000 UEL 4-noded element -20,000 UEL 4-noded element

-40,000 MATLAB 4-noded element -40,000 MATLAB 4-noded element
-60,000
-60,000
-80,000
node 3

node 4

-80,000
-100,000
-100,000
-120,000
-120,000
-140,000
-140,000 -160,000
-160,000 -180,000
-180,000 -200,000
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1
Displacement Displacement

(c) (d)

Figure 12: Plots for reaction moments in X-direction at the nodes with respect to nodal displacements.

Page 54 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

180,000 180,000
Abaqus S4 element Abaqus S4 element

Reaction moment in y-direction at

Reaction moment in y-direction at

160,000 UEL 4-noded element 160,000 UEL 4-noded element

140,000 MATLAB 4-noded element 140,000 MATLAB 4-noded element

120,000 120,000
node 2
node 1

100,000 100,000

80,000 80,000

60,000 60,000

40,000 40,000

20,000 20,000

0 0
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1
Displacement Displacement

(a) (b)
200,000 200,000
Abaqus S4 element Abaqus S4 element
Reaction moment in y-direction at

Reaction moment in y-direction at

180,000 180,000
UEL 4-noded element UEL 4-noded element
160,000 MATLAB 4-noded element 160,000 MATLAB 4-noded element
140,000 140,000
120,000 120,000
node 3

node 4

100,000 100,000
80,000 80,000
60,000 60,000
40,000 40,000
20,000 20,000
0 0
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1
Displacement Displacement

(c) (d)

Figure 13: Plots of reaction moments in Y-direction at the nodes with respect to nodal displacements.

Page 55 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

5.2 Testing multiple elements

Now we shall compare the results from Abaqus, our UEL implementation and the Matlab

code on a laterally loaded plate, which is clamped at the edges, and made up of multiple

elements. 400 elements in a (20 × 20) grid along X and Y directions are employed. The overall

dimensions of the plate in X-, Y- and Z-directions are again taken as 2 units, 2 units and 0.1

units, respectively. Maximum lateral pressure of 1 unit per unit area is applied3 .

Matlab code yields a maximum displacement of −1.15765 × 10−6 in Z-direction at the

central node. This is shown in Figure 14. With our implemented UEL subroutine,a similar

maximum displacement of −1.15765 × 10−6 in Z-direction is noted, shown in Figure 16.

Using S4 element analysis in Abaqus, we get a maximum displacement of −1.15751 × 10−6

in Z-direction, Figure 15. For Abaqus and UEL results, we use a scale factor of 106 to give a

consistent visual comparison with Matlab generated displacements. These results reveal an

excellent agreement.

Figure 17 shows the variation of maximum displacement in Z-direction as we gradually

vary the lateral pressure to 1 unit. The obtained results confirm that our implementation is

accurate and working as desired.

3
This is achieved by applying an equivalent level of concentrated force of 0.01 units over all the element
nodes.

Page 56 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

#10 -6
0

-0.2

-0.4
Z-axis

-0.6

-0.8

-1

-1.2
0
0.5 2
1 1.5
1
1.5
0.5
X-axis 2 0 Y-axis

Figure 14: Plate deflection from Matlab code.

Figure 15: Plate deflection from Abaqus S4 element analysis with a scale factor of 106 .

Page 57 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

Figure 16: Plate deflection from UEL with a scale factor of 106 .

Abaqus S4 element
-1E-07
UEL 4-noded element
Maximum displacement

-3E-07 MATLAB 4-noded element

-5E-07

-7E-07

-9E-07

-1.1E-06

-1.3E-06
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Lateral pressure
Figure 17: Comparison of deflections for a laterally loaded plate clamped at the edges.

Page 58 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

6 A CKNOWLEDGMENTS

Nikhil Padhye appreciates past technical discussions with his advisor Professor David M.

Parks for motivating the usage of Abaqus in finite element studies during the course of his

doctoral thesis [8].

Page 59 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

R EFERENCES

[1] Ver ABAQUS. 6.14 documentation. Dassault Systemes Simulia Corporation, 2014.

[2] F Auricchio and RL Taylor. A shear deformable plate element with an exact thin limit.

Computer Methods in Applied Mechanics and Engineering, 118(3-4):393–412, 1994.

[3] Klaus-Jürgen Bathe and Eduardo N Dvorkin. A four-node plate bending element based

on mindlin/reissner plate theory and a mixed interpolation. International Journal for

Numerical Methods in Engineering, 21(2):367–383, 1985.

[4] Jean-Louis Batoz and Mabrouk Ben Tahar. Evaluation of a new quadrilateral thin

plate bending element. International Journal for Numerical Methods in Engineering,

18(11):1655–1677, 1982.

[5] Ferreira. MATLAB Codes for Finite Element Analysis - Solids and Structures. Solid mechan-

ics & it applications 157. Springer, 2009.

[6] Bo Häggblad and Klaus-Jürgen Bathe. Specifications of boundary conditions for reiss-

ner/mindlin plate bending finite elements. International Journal for Numerical Methods

in Engineering, 30(5):981–1011, 1990.

[7] Thomas JR Hughes and TEi Tezduyar. Finite elements based upon mindlin plate theory

with particular reference to the four-node bilinear isoparametric element. Journal of

applied mechanics, 48(3):587–596, 1981.

[8] Nikhil Padhye. Sub-Tg, solid-state, plasticity-induced bonding of polymeric films and

continuous forming. PhD thesis, Massachusetts Institute of Technology, 2015.

Page 60 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

[9] Erik Saether and Alexander Tessler. User-defined subroutine for implementation of

higher-order shell element in abaqus. Technical report, ARMY RESEARCH LAB ADELPHI

MD, 1993.

Page 61 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

Abaqus job input file for UEL with single element

1 *Heading
2 ** Job name: abaqus_4_elements_lift_up Model name: Model-1
3 ** Generated by: Abaqus/CAE 2016.HF2
4 *Preprint, echo=NO, model=NO, history=NO, contact=NO
5 **
6 ** PARTS
7 **
8 *Part, name=Plate
9 *Node
10 1, 0., 0., 0.
11 2, 2., 0., 0.
12 3, 0., 2., 0.
13 4, 2., 2., 0.
14 ******************************************************
15 ******************************************************
16 ******************************************************
17 *USER ELEMENT, TYPE=U1, NODES=4, COORDINATES=3, VAR=40,
18 INTEGRATION=4
19 3,4,5
20 *ELEMENT,TYPE=U1,ELSET=SOLID
21 1, 1,2,4,3
22 *UEL PROPERTY, ELSET=SOLID
23 **************************************
24 *End Part
25 **
26 **
27 ** ASSEMBLY
28 **
29 *Assembly, name=Assembly
30 **
31 *Instance, name=Plate-1, part=Plate
32 *End Instance
33 **
34 *Nset, nset=all_elements, instance=Plate-1, generate
35 1, 4, 1
36 *Elset, elset=all_elements, instance=Plate-1
37 1,
38 *Nset, nset=_PickedSet13, internal, instance=Plate-1
39 4,
40 *Nset, nset=_PickedSet14, internal, instance=Plate-1, generate
41 1, 3, 1
42 *End Assembly
43 **
44 **
45 ** BOUNDARY CONDITIONS
46 **
47 ** Name: fix_nodes Type: Symmetry/Antisymmetry/Encastre
48 *Boundary
49 _PickedSet14, ENCASTRE
50 ** ----------------------------------------------------------------
51 **
52 ** STEP: loading
53 **
54 *Step, name=loading, nlgeom=NO, inc=100000
55 *Static, direct
56 0.1, 1.,
57 **
58 ** BOUNDARY CONDITIONS
59 **
60 ** Name: lift_up Type: Displacement/Rotation
61 *Boundary
62 _PickedSet13, 3, 3, 0.1
63 _PickedSet13, 4, 4
64 _PickedSet13, 5, 5
65 **
66 ** OUTPUT REQUESTS
67 **
68 *Restart, write, frequency=0
69 **

Page 62 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

70 ** FIELD OUTPUT: F-Output-2

71 **
72 *Output, field
73 *Node Output
74 CF, RF, RM, RT, TF, VF
75 *Element Output, directions=YES
76 BF, CENTMAG, CENTRIFMAG, CORIOMAG, ESF1, GRAV, HP, NFORC, NFORCSO, P, ROTAMAG, SF,
TRNOR, TRSHR
77 **
78 ** FIELD OUTPUT: F-Output-1
79 **
80 *Output, field, variable=PRESELECT
81 **
82 ** HISTORY OUTPUT: H-Output-2
83 **
84 *Output, history
85 *Element Output
86 IRF1, IRF2, IRF3, IRM1, IRM2, IRM3
87 **
88 ** HISTORY OUTPUT: H-Output-1
89 **
90 *Output, history, variable=PRESELECT
91 *End Step
92

Page 63 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

Modified Matlab main routine for single element

1 %................................................................
2
3 % MATLAB codes for Finite Element Analysis
4 % problem19.m
5 % Mindlin plate in bending
6 % antonio ferreira 2008
7
8 % clear memory
9 clc
10 clear *;
11 clear all;
12 colordef white;
13 format long e;
14 close all;
15
16 % materials
17 E = 200000000; poisson = 0.30; kapa=5/6;
18 thickness=0.1;
19 I=thickness^3/12;
20
21 % matrix C
22 % bending part
23 C_bending=I*E/(1-poisson^2)*...
24 [1 poisson 0;poisson 1 0;0 0 (1-poisson)/2];
25 % shear part
26 C_shear=kapa*thickness*E/2/(1+poisson)*eye(2);
27
28 % load
29 P = -1;
30
31 %Mesh generation
32 L = 2;
33 numberElementsX=1;
34 numberElementsY=1;
35 numberElements=numberElementsX*numberElementsY;
36 %
37 [nodeCoordinates, elementNodes] = ...
38 rectangularMesh(L,L,numberElementsX,numberElementsY);
39 xx=nodeCoordinates(:,1);
40 yy=nodeCoordinates(:,2);
41
42 %drawingMesh(nodeCoordinates,elementNodes,'Q4','k-');
43 %axis off
44 numberNodes=size(xx,1);
45
46 % GDof: global number of degrees of freedom
47 GDof=3*numberNodes;
48
49 % computation of the system stiffness matrix and force vector
50 [stiffness]=...
51 formStiffnessMatrixMindlinQ4(GDof,numberElements,...
52 elementNodes,numberNodes,nodeCoordinates,C_shear,...
53 C_bending,thickness,I);
54
55 timestep = 1;
56
57 U=[
58 0
59 0
60 0
61 0.1*timestep
62 0
63 0
64 0
65 0
66 0
67 0
68 0
69 0];

Page 64 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

70
71
72 stiffness*U

Page 65 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

Abaqus job input file for UEL with multiple elements

1 *Heading
2 ** Job name: push_center_Down Model name: Model-1
3 ** Generated by: Abaqus/CAE 2016.HF2
4 *Preprint, echo=NO, model=NO, history=NO, contact=NO
5 **
6 ** PARTS
7 **
8 *Part, name=Plate
9 *Node
10 1, 0., 0., 0.
11 2, 0.100000001, 0., 0.
12 3, 0.200000003, 0., 0.
13 4, 0.300000012, 0., 0.
14 5, 0.400000006, 0., 0.
15 6, 0.5, 0., 0.
16 7, 0.600000024, 0., 0.
17 8, 0.699999988, 0., 0.
18 9, 0.800000012, 0., 0.
19 10, 0.899999976, 0., 0.
20 11, 1., 0., 0.
21 12, 1.10000002, 0., 0.
22 13, 1.20000005, 0., 0.
23 14, 1.29999995, 0., 0.
24 15, 1.39999998, 0., 0.
25 16, 1.5, 0., 0.
26 17, 1.60000002, 0., 0.
27 18, 1.70000005, 0., 0.
28 19, 1.79999995, 0., 0.
29 20, 1.89999998, 0., 0.
30 21, 2., 0., 0.
31 22, 0., 0.100000001, 0.
32 23, 0.100000001, 0.100000001, 0.
33 24, 0.200000003, 0.100000001, 0.
34 25, 0.300000012, 0.100000001, 0.
35 26, 0.400000006, 0.100000001, 0.
36 27, 0.5, 0.100000001, 0.
37 28, 0.600000024, 0.100000001, 0.
38 29, 0.699999988, 0.100000001, 0.
39 30, 0.800000012, 0.100000001, 0.
40 31, 0.899999976, 0.100000001, 0.
41 32, 1., 0.100000001, 0.
42 33, 1.10000002, 0.100000001, 0.
43 34, 1.20000005, 0.100000001, 0.
44 35, 1.29999995, 0.100000001, 0.
45 36, 1.39999998, 0.100000001, 0.
46 37, 1.5, 0.100000001, 0.
47 38, 1.60000002, 0.100000001, 0.
48 39, 1.70000005, 0.100000001, 0.
49 40, 1.79999995, 0.100000001, 0.
50 41, 1.89999998, 0.100000001, 0.
51 42, 2., 0.100000001, 0.
52 43, 0., 0.200000003, 0.
53 44, 0.100000001, 0.200000003, 0.
54 45, 0.200000003, 0.200000003, 0.
55 46, 0.300000012, 0.200000003, 0.
56 47, 0.400000006, 0.200000003, 0.
57 48, 0.5, 0.200000003, 0.
58 49, 0.600000024, 0.200000003, 0.
59 50, 0.699999988, 0.200000003, 0.
60 51, 0.800000012, 0.200000003, 0.
61 52, 0.899999976, 0.200000003, 0.
62 53, 1., 0.200000003, 0.
63 54, 1.10000002, 0.200000003, 0.
64 55, 1.20000005, 0.200000003, 0.
65 56, 1.29999995, 0.200000003, 0.
66 57, 1.39999998, 0.200000003, 0.
67 58, 1.5, 0.200000003, 0.
68 59, 1.60000002, 0.200000003, 0.
69 60, 1.70000005, 0.200000003, 0.

Page 66 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

70 61, 1.79999995, 0.200000003, 0.

71 62, 1.89999998, 0.200000003, 0.
72 63, 2., 0.200000003, 0.
73 64, 0., 0.300000012, 0.
74 65, 0.100000001, 0.300000012, 0.
75 66, 0.200000003, 0.300000012, 0.
76 67, 0.300000012, 0.300000012, 0.
77 68, 0.400000006, 0.300000012, 0.
78 69, 0.5, 0.300000012, 0.
79 70, 0.600000024, 0.300000012, 0.
80 71, 0.699999988, 0.300000012, 0.
81 72, 0.800000012, 0.300000012, 0.
82 73, 0.899999976, 0.300000012, 0.
83 74, 1., 0.300000012, 0.
84 75, 1.10000002, 0.300000012, 0.
85 76, 1.20000005, 0.300000012, 0.
86 77, 1.29999995, 0.300000012, 0.
87 78, 1.39999998, 0.300000012, 0.
88 79, 1.5, 0.300000012, 0.
89 80, 1.60000002, 0.300000012, 0.
90 81, 1.70000005, 0.300000012, 0.
91 82, 1.79999995, 0.300000012, 0.
92 83, 1.89999998, 0.300000012, 0.
93 84, 2., 0.300000012, 0.
94 85, 0., 0.400000006, 0.
95 86, 0.100000001, 0.400000006, 0.
96 87, 0.200000003, 0.400000006, 0.
97 88, 0.300000012, 0.400000006, 0.
98 89, 0.400000006, 0.400000006, 0.
99 90, 0.5, 0.400000006, 0.
100 91, 0.600000024, 0.400000006, 0.
101 92, 0.699999988, 0.400000006, 0.
102 93, 0.800000012, 0.400000006, 0.
103 94, 0.899999976, 0.400000006, 0.
104 95, 1., 0.400000006, 0.
105 96, 1.10000002, 0.400000006, 0.
106 97, 1.20000005, 0.400000006, 0.
107 98, 1.29999995, 0.400000006, 0.
108 99, 1.39999998, 0.400000006, 0.
109 100, 1.5, 0.400000006, 0.
110 101, 1.60000002, 0.400000006, 0.
111 102, 1.70000005, 0.400000006, 0.
112 103, 1.79999995, 0.400000006, 0.
113 104, 1.89999998, 0.400000006, 0.
114 105, 2., 0.400000006, 0.
115 106, 0., 0.5, 0.
116 107, 0.100000001, 0.5, 0.
117 108, 0.200000003, 0.5, 0.
118 109, 0.300000012, 0.5, 0.
119 110, 0.400000006, 0.5, 0.
120 111, 0.5, 0.5, 0.
121 112, 0.600000024, 0.5, 0.
122 113, 0.699999988, 0.5, 0.
123 114, 0.800000012, 0.5, 0.
124 115, 0.899999976, 0.5, 0.
125 116, 1., 0.5, 0.
126 117, 1.10000002, 0.5, 0.
127 118, 1.20000005, 0.5, 0.
128 119, 1.29999995, 0.5, 0.
129 120, 1.39999998, 0.5, 0.
130 121, 1.5, 0.5, 0.
131 122, 1.60000002, 0.5, 0.
132 123, 1.70000005, 0.5, 0.
133 124, 1.79999995, 0.5, 0.
134 125, 1.89999998, 0.5, 0.
135 126, 2., 0.5, 0.
136 127, 0., 0.600000024, 0.
137 128, 0.100000001, 0.600000024, 0.
138 129, 0.200000003, 0.600000024, 0.

Page 67 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

139 130, 0.300000012, 0.600000024, 0.

140 131, 0.400000006, 0.600000024, 0.
141 132, 0.5, 0.600000024, 0.
142 133, 0.600000024, 0.600000024, 0.
143 134, 0.699999988, 0.600000024, 0.
144 135, 0.800000012, 0.600000024, 0.
145 136, 0.899999976, 0.600000024, 0.
146 137, 1., 0.600000024, 0.
147 138, 1.10000002, 0.600000024, 0.
148 139, 1.20000005, 0.600000024, 0.
149 140, 1.29999995, 0.600000024, 0.
150 141, 1.39999998, 0.600000024, 0.
151 142, 1.5, 0.600000024, 0.
152 143, 1.60000002, 0.600000024, 0.
153 144, 1.70000005, 0.600000024, 0.
154 145, 1.79999995, 0.600000024, 0.
155 146, 1.89999998, 0.600000024, 0.
156 147, 2., 0.600000024, 0.
157 148, 0., 0.699999988, 0.
158 149, 0.100000001, 0.699999988, 0.
159 150, 0.200000003, 0.699999988, 0.
160 151, 0.300000012, 0.699999988, 0.
161 152, 0.400000006, 0.699999988, 0.
162 153, 0.5, 0.699999988, 0.
163 154, 0.600000024, 0.699999988, 0.
164 155, 0.699999988, 0.699999988, 0.
165 156, 0.800000012, 0.699999988, 0.
166 157, 0.899999976, 0.699999988, 0.
167 158, 1., 0.699999988, 0.
168 159, 1.10000002, 0.699999988, 0.
169 160, 1.20000005, 0.699999988, 0.
170 161, 1.29999995, 0.699999988, 0.
171 162, 1.39999998, 0.699999988, 0.
172 163, 1.5, 0.699999988, 0.
173 164, 1.60000002, 0.699999988, 0.
174 165, 1.70000005, 0.699999988, 0.
175 166, 1.79999995, 0.699999988, 0.
176 167, 1.89999998, 0.699999988, 0.
177 168, 2., 0.699999988, 0.
178 169, 0., 0.800000012, 0.
179 170, 0.100000001, 0.800000012, 0.
180 171, 0.200000003, 0.800000012, 0.
181 172, 0.300000012, 0.800000012, 0.
182 173, 0.400000006, 0.800000012, 0.
183 174, 0.5, 0.800000012, 0.
184 175, 0.600000024, 0.800000012, 0.
185 176, 0.699999988, 0.800000012, 0.
186 177, 0.800000012, 0.800000012, 0.
187 178, 0.899999976, 0.800000012, 0.
188 179, 1., 0.800000012, 0.
189 180, 1.10000002, 0.800000012, 0.
190 181, 1.20000005, 0.800000012, 0.
191 182, 1.29999995, 0.800000012, 0.
192 183, 1.39999998, 0.800000012, 0.
193 184, 1.5, 0.800000012, 0.
194 185, 1.60000002, 0.800000012, 0.
195 186, 1.70000005, 0.800000012, 0.
196 187, 1.79999995, 0.800000012, 0.
197 188, 1.89999998, 0.800000012, 0.
198 189, 2., 0.800000012, 0.
199 190, 0., 0.899999976, 0.
200 191, 0.100000001, 0.899999976, 0.
201 192, 0.200000003, 0.899999976, 0.
202 193, 0.300000012, 0.899999976, 0.
203 194, 0.400000006, 0.899999976, 0.
204 195, 0.5, 0.899999976, 0.
205 196, 0.600000024, 0.899999976, 0.
206 197, 0.699999988, 0.899999976, 0.
207 198, 0.800000012, 0.899999976, 0.

Page 68 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

208 199, 0.899999976, 0.899999976, 0.

209 200, 1., 0.899999976, 0.
210 201, 1.10000002, 0.899999976, 0.
211 202, 1.20000005, 0.899999976, 0.
212 203, 1.29999995, 0.899999976, 0.
213 204, 1.39999998, 0.899999976, 0.
214 205, 1.5, 0.899999976, 0.
215 206, 1.60000002, 0.899999976, 0.
216 207, 1.70000005, 0.899999976, 0.
217 208, 1.79999995, 0.899999976, 0.
218 209, 1.89999998, 0.899999976, 0.
219 210, 2., 0.899999976, 0.
220 211, 0., 1., 0.
221 212, 0.100000001, 1., 0.
222 213, 0.200000003, 1., 0.
223 214, 0.300000012, 1., 0.
224 215, 0.400000006, 1., 0.
225 216, 0.5, 1., 0.
226 217, 0.600000024, 1., 0.
227 218, 0.699999988, 1., 0.
228 219, 0.800000012, 1., 0.
229 220, 0.899999976, 1., 0.
230 221, 1., 1., 0.
231 222, 1.10000002, 1., 0.
232 223, 1.20000005, 1., 0.
233 224, 1.29999995, 1., 0.
234 225, 1.39999998, 1., 0.
235 226, 1.5, 1., 0.
236 227, 1.60000002, 1., 0.
237 228, 1.70000005, 1., 0.
238 229, 1.79999995, 1., 0.
239 230, 1.89999998, 1., 0.
240 231, 2., 1., 0.
241 232, 0., 1.10000002, 0.
242 233, 0.100000001, 1.10000002, 0.
243 234, 0.200000003, 1.10000002, 0.
244 235, 0.300000012, 1.10000002, 0.
245 236, 0.400000006, 1.10000002, 0.
246 237, 0.5, 1.10000002, 0.
247 238, 0.600000024, 1.10000002, 0.
248 239, 0.699999988, 1.10000002, 0.
249 240, 0.800000012, 1.10000002, 0.
250 241, 0.899999976, 1.10000002, 0.
251 242, 1., 1.10000002, 0.
252 243, 1.10000002, 1.10000002, 0.
253 244, 1.20000005, 1.10000002, 0.
254 245, 1.29999995, 1.10000002, 0.
255 246, 1.39999998, 1.10000002, 0.
256 247, 1.5, 1.10000002, 0.
257 248, 1.60000002, 1.10000002, 0.
258 249, 1.70000005, 1.10000002, 0.
259 250, 1.79999995, 1.10000002, 0.
260 251, 1.89999998, 1.10000002, 0.
261 252, 2., 1.10000002, 0.
262 253, 0., 1.20000005, 0.
263 254, 0.100000001, 1.20000005, 0.
264 255, 0.200000003, 1.20000005, 0.
265 256, 0.300000012, 1.20000005, 0.
266 257, 0.400000006, 1.20000005, 0.
267 258, 0.5, 1.20000005, 0.
268 259, 0.600000024, 1.20000005, 0.
269 260, 0.699999988, 1.20000005, 0.
270 261, 0.800000012, 1.20000005, 0.
271 262, 0.899999976, 1.20000005, 0.
272 263, 1., 1.20000005, 0.
273 264, 1.10000002, 1.20000005, 0.
274 265, 1.20000005, 1.20000005, 0.
275 266, 1.29999995, 1.20000005, 0.
276 267, 1.39999998, 1.20000005, 0.

Page 69 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

277 268, 1.5, 1.20000005, 0.

278 269, 1.60000002, 1.20000005, 0.
279 270, 1.70000005, 1.20000005, 0.
280 271, 1.79999995, 1.20000005, 0.
281 272, 1.89999998, 1.20000005, 0.
282 273, 2., 1.20000005, 0.
283 274, 0., 1.29999995, 0.
284 275, 0.100000001, 1.29999995, 0.
285 276, 0.200000003, 1.29999995, 0.
286 277, 0.300000012, 1.29999995, 0.
287 278, 0.400000006, 1.29999995, 0.
288 279, 0.5, 1.29999995, 0.
289 280, 0.600000024, 1.29999995, 0.
290 281, 0.699999988, 1.29999995, 0.
291 282, 0.800000012, 1.29999995, 0.
292 283, 0.899999976, 1.29999995, 0.
293 284, 1., 1.29999995, 0.
294 285, 1.10000002, 1.29999995, 0.
295 286, 1.20000005, 1.29999995, 0.
296 287, 1.29999995, 1.29999995, 0.
297 288, 1.39999998, 1.29999995, 0.
298 289, 1.5, 1.29999995, 0.
299 290, 1.60000002, 1.29999995, 0.
300 291, 1.70000005, 1.29999995, 0.
301 292, 1.79999995, 1.29999995, 0.
302 293, 1.89999998, 1.29999995, 0.
303 294, 2., 1.29999995, 0.
304 295, 0., 1.39999998, 0.
305 296, 0.100000001, 1.39999998, 0.
306 297, 0.200000003, 1.39999998, 0.
307 298, 0.300000012, 1.39999998, 0.
308 299, 0.400000006, 1.39999998, 0.
309 300, 0.5, 1.39999998, 0.
310 301, 0.600000024, 1.39999998, 0.
311 302, 0.699999988, 1.39999998, 0.
312 303, 0.800000012, 1.39999998, 0.
313 304, 0.899999976, 1.39999998, 0.
314 305, 1., 1.39999998, 0.
315 306, 1.10000002, 1.39999998, 0.
316 307, 1.20000005, 1.39999998, 0.
317 308, 1.29999995, 1.39999998, 0.
318 309, 1.39999998, 1.39999998, 0.
319 310, 1.5, 1.39999998, 0.
320 311, 1.60000002, 1.39999998, 0.
321 312, 1.70000005, 1.39999998, 0.
322 313, 1.79999995, 1.39999998, 0.
323 314, 1.89999998, 1.39999998, 0.
324 315, 2., 1.39999998, 0.
325 316, 0., 1.5, 0.
326 317, 0.100000001, 1.5, 0.
327 318, 0.200000003, 1.5, 0.
328 319, 0.300000012, 1.5, 0.
329 320, 0.400000006, 1.5, 0.
330 321, 0.5, 1.5, 0.
331 322, 0.600000024, 1.5, 0.
332 323, 0.699999988, 1.5, 0.
333 324, 0.800000012, 1.5, 0.
334 325, 0.899999976, 1.5, 0.
335 326, 1., 1.5, 0.
336 327, 1.10000002, 1.5, 0.
337 328, 1.20000005, 1.5, 0.
338 329, 1.29999995, 1.5, 0.
339 330, 1.39999998, 1.5, 0.
340 331, 1.5, 1.5, 0.
341 332, 1.60000002, 1.5, 0.
342 333, 1.70000005, 1.5, 0.
343 334, 1.79999995, 1.5, 0.
344 335, 1.89999998, 1.5, 0.
345 336, 2., 1.5, 0.

Page 70 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

346 337, 0., 1.60000002, 0.

347 338, 0.100000001, 1.60000002, 0.
348 339, 0.200000003, 1.60000002, 0.
349 340, 0.300000012, 1.60000002, 0.
350 341, 0.400000006, 1.60000002, 0.
351 342, 0.5, 1.60000002, 0.
352 343, 0.600000024, 1.60000002, 0.
353 344, 0.699999988, 1.60000002, 0.
354 345, 0.800000012, 1.60000002, 0.
355 346, 0.899999976, 1.60000002, 0.
356 347, 1., 1.60000002, 0.
357 348, 1.10000002, 1.60000002, 0.
358 349, 1.20000005, 1.60000002, 0.
359 350, 1.29999995, 1.60000002, 0.
360 351, 1.39999998, 1.60000002, 0.
361 352, 1.5, 1.60000002, 0.
362 353, 1.60000002, 1.60000002, 0.
363 354, 1.70000005, 1.60000002, 0.
364 355, 1.79999995, 1.60000002, 0.
365 356, 1.89999998, 1.60000002, 0.
366 357, 2., 1.60000002, 0.
367 358, 0., 1.70000005, 0.
368 359, 0.100000001, 1.70000005, 0.
369 360, 0.200000003, 1.70000005, 0.
370 361, 0.300000012, 1.70000005, 0.
371 362, 0.400000006, 1.70000005, 0.
372 363, 0.5, 1.70000005, 0.
373 364, 0.600000024, 1.70000005, 0.
374 365, 0.699999988, 1.70000005, 0.
375 366, 0.800000012, 1.70000005, 0.
376 367, 0.899999976, 1.70000005, 0.
377 368, 1., 1.70000005, 0.
378 369, 1.10000002, 1.70000005, 0.
379 370, 1.20000005, 1.70000005, 0.
380 371, 1.29999995, 1.70000005, 0.
381 372, 1.39999998, 1.70000005, 0.
382 373, 1.5, 1.70000005, 0.
383 374, 1.60000002, 1.70000005, 0.
384 375, 1.70000005, 1.70000005, 0.
385 376, 1.79999995, 1.70000005, 0.
386 377, 1.89999998, 1.70000005, 0.
387 378, 2., 1.70000005, 0.
388 379, 0., 1.79999995, 0.
389 380, 0.100000001, 1.79999995, 0.
390 381, 0.200000003, 1.79999995, 0.
391 382, 0.300000012, 1.79999995, 0.
392 383, 0.400000006, 1.79999995, 0.
393 384, 0.5, 1.79999995, 0.
394 385, 0.600000024, 1.79999995, 0.
395 386, 0.699999988, 1.79999995, 0.
396 387, 0.800000012, 1.79999995, 0.
397 388, 0.899999976, 1.79999995, 0.
398 389, 1., 1.79999995, 0.
399 390, 1.10000002, 1.79999995, 0.
400 391, 1.20000005, 1.79999995, 0.
401 392, 1.29999995, 1.79999995, 0.
402 393, 1.39999998, 1.79999995, 0.
403 394, 1.5, 1.79999995, 0.
404 395, 1.60000002, 1.79999995, 0.
405 396, 1.70000005, 1.79999995, 0.
406 397, 1.79999995, 1.79999995, 0.
407 398, 1.89999998, 1.79999995, 0.
408 399, 2., 1.79999995, 0.
409 400, 0., 1.89999998, 0.
410 401, 0.100000001, 1.89999998, 0.
411 402, 0.200000003, 1.89999998, 0.
412 403, 0.300000012, 1.89999998, 0.
413 404, 0.400000006, 1.89999998, 0.
414 405, 0.5, 1.89999998, 0.

Page 71 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

415 406, 0.600000024, 1.89999998, 0.

416 407, 0.699999988, 1.89999998, 0.
417 408, 0.800000012, 1.89999998, 0.
418 409, 0.899999976, 1.89999998, 0.
419 410, 1., 1.89999998, 0.
420 411, 1.10000002, 1.89999998, 0.
421 412, 1.20000005, 1.89999998, 0.
422 413, 1.29999995, 1.89999998, 0.
423 414, 1.39999998, 1.89999998, 0.
424 415, 1.5, 1.89999998, 0.
425 416, 1.60000002, 1.89999998, 0.
426 417, 1.70000005, 1.89999998, 0.
427 418, 1.79999995, 1.89999998, 0.
428 419, 1.89999998, 1.89999998, 0.
429 420, 2., 1.89999998, 0.
430 421, 0., 2., 0.
431 422, 0.100000001, 2., 0.
432 423, 0.200000003, 2., 0.
433 424, 0.300000012, 2., 0.
434 425, 0.400000006, 2., 0.
435 426, 0.5, 2., 0.
436 427, 0.600000024, 2., 0.
437 428, 0.699999988, 2., 0.
438 429, 0.800000012, 2., 0.
439 430, 0.899999976, 2., 0.
440 431, 1., 2., 0.
441 432, 1.10000002, 2., 0.
442 433, 1.20000005, 2., 0.
443 434, 1.29999995, 2., 0.
444 435, 1.39999998, 2., 0.
445 436, 1.5, 2., 0.
446 437, 1.60000002, 2., 0.
447 438, 1.70000005, 2., 0.
448 439, 1.79999995, 2., 0.
449 440, 1.89999998, 2., 0.
450 441, 2., 2., 0.
451 ******************************************************
452 ******************************************************
453 ******************************************************
454 *USER ELEMENT, TYPE=U1, NODES=4, COORDINATES=3, VAR=40,
455 INTEGRATION=4
456 3,4,5
457 *ELEMENT,TYPE=U1,ELSET=SOLID
458 1, 1, 2, 23, 22
459 2, 2, 3, 24, 23
460 3, 3, 4, 25, 24
461 4, 4, 5, 26, 25
462 5, 5, 6, 27, 26
463 6, 6, 7, 28, 27
464 7, 7, 8, 29, 28
465 8, 8, 9, 30, 29
466 9, 9, 10, 31, 30
467 10, 10, 11, 32, 31
468 11, 11, 12, 33, 32
469 12, 12, 13, 34, 33
470 13, 13, 14, 35, 34
471 14, 14, 15, 36, 35
472 15, 15, 16, 37, 36
473 16, 16, 17, 38, 37
474 17, 17, 18, 39, 38
475 18, 18, 19, 40, 39
476 19, 19, 20, 41, 40
477 20, 20, 21, 42, 41
478 21, 22, 23, 44, 43
479 22, 23, 24, 45, 44
480 23, 24, 25, 46, 45
481 24, 25, 26, 47, 46
482 25, 26, 27, 48, 47
483 26, 27, 28, 49, 48

Page 72 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

484 27, 28, 29, 50, 49

485 28, 29, 30, 51, 50
486 29, 30, 31, 52, 51
487 30, 31, 32, 53, 52
488 31, 32, 33, 54, 53
489 32, 33, 34, 55, 54
490 33, 34, 35, 56, 55
491 34, 35, 36, 57, 56
492 35, 36, 37, 58, 57
493 36, 37, 38, 59, 58
494 37, 38, 39, 60, 59
495 38, 39, 40, 61, 60
496 39, 40, 41, 62, 61
497 40, 41, 42, 63, 62
498 41, 43, 44, 65, 64
499 42, 44, 45, 66, 65
500 43, 45, 46, 67, 66
501 44, 46, 47, 68, 67
502 45, 47, 48, 69, 68
503 46, 48, 49, 70, 69
504 47, 49, 50, 71, 70
505 48, 50, 51, 72, 71
506 49, 51, 52, 73, 72
507 50, 52, 53, 74, 73
508 51, 53, 54, 75, 74
509 52, 54, 55, 76, 75
510 53, 55, 56, 77, 76
511 54, 56, 57, 78, 77
512 55, 57, 58, 79, 78
513 56, 58, 59, 80, 79
514 57, 59, 60, 81, 80
515 58, 60, 61, 82, 81
516 59, 61, 62, 83, 82
517 60, 62, 63, 84, 83
518 61, 64, 65, 86, 85
519 62, 65, 66, 87, 86
520 63, 66, 67, 88, 87
521 64, 67, 68, 89, 88
522 65, 68, 69, 90, 89
523 66, 69, 70, 91, 90
524 67, 70, 71, 92, 91
525 68, 71, 72, 93, 92
526 69, 72, 73, 94, 93
527 70, 73, 74, 95, 94
528 71, 74, 75, 96, 95
529 72, 75, 76, 97, 96
530 73, 76, 77, 98, 97
531 74, 77, 78, 99, 98
532 75, 78, 79, 100, 99
533 76, 79, 80, 101, 100
534 77, 80, 81, 102, 101
535 78, 81, 82, 103, 102
536 79, 82, 83, 104, 103
537 80, 83, 84, 105, 104
538 81, 85, 86, 107, 106
539 82, 86, 87, 108, 107
540 83, 87, 88, 109, 108
541 84, 88, 89, 110, 109
542 85, 89, 90, 111, 110
543 86, 90, 91, 112, 111
544 87, 91, 92, 113, 112
545 88, 92, 93, 114, 113
546 89, 93, 94, 115, 114
547 90, 94, 95, 116, 115
548 91, 95, 96, 117, 116
549 92, 96, 97, 118, 117
550 93, 97, 98, 119, 118
551 94, 98, 99, 120, 119
552 95, 99, 100, 121, 120

Page 73 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

553 96, 100, 101, 122, 121

554 97, 101, 102, 123, 122
555 98, 102, 103, 124, 123
556 99, 103, 104, 125, 124
557 100, 104, 105, 126, 125
558 101, 106, 107, 128, 127
559 102, 107, 108, 129, 128
560 103, 108, 109, 130, 129
561 104, 109, 110, 131, 130
562 105, 110, 111, 132, 131
563 106, 111, 112, 133, 132
564 107, 112, 113, 134, 133
565 108, 113, 114, 135, 134
566 109, 114, 115, 136, 135
567 110, 115, 116, 137, 136
568 111, 116, 117, 138, 137
569 112, 117, 118, 139, 138
570 113, 118, 119, 140, 139
571 114, 119, 120, 141, 140
572 115, 120, 121, 142, 141
573 116, 121, 122, 143, 142
574 117, 122, 123, 144, 143
575 118, 123, 124, 145, 144
576 119, 124, 125, 146, 145
577 120, 125, 126, 147, 146
578 121, 127, 128, 149, 148
579 122, 128, 129, 150, 149
580 123, 129, 130, 151, 150
581 124, 130, 131, 152, 151
582 125, 131, 132, 153, 152
583 126, 132, 133, 154, 153
584 127, 133, 134, 155, 154
585 128, 134, 135, 156, 155
586 129, 135, 136, 157, 156
587 130, 136, 137, 158, 157
588 131, 137, 138, 159, 158
589 132, 138, 139, 160, 159
590 133, 139, 140, 161, 160
591 134, 140, 141, 162, 161
592 135, 141, 142, 163, 162
593 136, 142, 143, 164, 163
594 137, 143, 144, 165, 164
595 138, 144, 145, 166, 165
596 139, 145, 146, 167, 166
597 140, 146, 147, 168, 167
598 141, 148, 149, 170, 169
599 142, 149, 150, 171, 170
600 143, 150, 151, 172, 171
601 144, 151, 152, 173, 172
602 145, 152, 153, 174, 173
603 146, 153, 154, 175, 174
604 147, 154, 155, 176, 175
605 148, 155, 156, 177, 176
606 149, 156, 157, 178, 177
607 150, 157, 158, 179, 178
608 151, 158, 159, 180, 179
609 152, 159, 160, 181, 180
610 153, 160, 161, 182, 181
611 154, 161, 162, 183, 182
612 155, 162, 163, 184, 183
613 156, 163, 164, 185, 184
614 157, 164, 165, 186, 185
615 158, 165, 166, 187, 186
616 159, 166, 167, 188, 187
617 160, 167, 168, 189, 188
618 161, 169, 170, 191, 190
619 162, 170, 171, 192, 191
620 163, 171, 172, 193, 192
621 164, 172, 173, 194, 193

Page 74 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

622 165, 173, 174, 195, 194

623 166, 174, 175, 196, 195
624 167, 175, 176, 197, 196
625 168, 176, 177, 198, 197
626 169, 177, 178, 199, 198
627 170, 178, 179, 200, 199
628 171, 179, 180, 201, 200
629 172, 180, 181, 202, 201
630 173, 181, 182, 203, 202
631 174, 182, 183, 204, 203
632 175, 183, 184, 205, 204
633 176, 184, 185, 206, 205
634 177, 185, 186, 207, 206
635 178, 186, 187, 208, 207
636 179, 187, 188, 209, 208
637 180, 188, 189, 210, 209
638 181, 190, 191, 212, 211
639 182, 191, 192, 213, 212
640 183, 192, 193, 214, 213
641 184, 193, 194, 215, 214
642 185, 194, 195, 216, 215
643 186, 195, 196, 217, 216
644 187, 196, 197, 218, 217
645 188, 197, 198, 219, 218
646 189, 198, 199, 220, 219
647 190, 199, 200, 221, 220
648 191, 200, 201, 222, 221
649 192, 201, 202, 223, 222
650 193, 202, 203, 224, 223
651 194, 203, 204, 225, 224
652 195, 204, 205, 226, 225
653 196, 205, 206, 227, 226
654 197, 206, 207, 228, 227
655 198, 207, 208, 229, 228
656 199, 208, 209, 230, 229
657 200, 209, 210, 231, 230
658 201, 211, 212, 233, 232
659 202, 212, 213, 234, 233
660 203, 213, 214, 235, 234
661 204, 214, 215, 236, 235
662 205, 215, 216, 237, 236
663 206, 216, 217, 238, 237
664 207, 217, 218, 239, 238
665 208, 218, 219, 240, 239
666 209, 219, 220, 241, 240
667 210, 220, 221, 242, 241
668 211, 221, 222, 243, 242
669 212, 222, 223, 244, 243
670 213, 223, 224, 245, 244
671 214, 224, 225, 246, 245
672 215, 225, 226, 247, 246
673 216, 226, 227, 248, 247
674 217, 227, 228, 249, 248
675 218, 228, 229, 250, 249
676 219, 229, 230, 251, 250
677 220, 230, 231, 252, 251
678 221, 232, 233, 254, 253
679 222, 233, 234, 255, 254
680 223, 234, 235, 256, 255
681 224, 235, 236, 257, 256
682 225, 236, 237, 258, 257
683 226, 237, 238, 259, 258
684 227, 238, 239, 260, 259
685 228, 239, 240, 261, 260
686 229, 240, 241, 262, 261
687 230, 241, 242, 263, 262
688 231, 242, 243, 264, 263
689 232, 243, 244, 265, 264
690 233, 244, 245, 266, 265

Page 75 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

691 234, 245, 246, 267, 266

692 235, 246, 247, 268, 267
693 236, 247, 248, 269, 268
694 237, 248, 249, 270, 269
695 238, 249, 250, 271, 270
696 239, 250, 251, 272, 271
697 240, 251, 252, 273, 272
698 241, 253, 254, 275, 274
699 242, 254, 255, 276, 275
700 243, 255, 256, 277, 276
701 244, 256, 257, 278, 277
702 245, 257, 258, 279, 278
703 246, 258, 259, 280, 279
704 247, 259, 260, 281, 280
705 248, 260, 261, 282, 281
706 249, 261, 262, 283, 282
707 250, 262, 263, 284, 283
708 251, 263, 264, 285, 284
709 252, 264, 265, 286, 285
710 253, 265, 266, 287, 286
711 254, 266, 267, 288, 287
712 255, 267, 268, 289, 288
713 256, 268, 269, 290, 289
714 257, 269, 270, 291, 290
715 258, 270, 271, 292, 291
716 259, 271, 272, 293, 292
717 260, 272, 273, 294, 293
718 261, 274, 275, 296, 295
719 262, 275, 276, 297, 296
720 263, 276, 277, 298, 297
721 264, 277, 278, 299, 298
722 265, 278, 279, 300, 299
723 266, 279, 280, 301, 300
724 267, 280, 281, 302, 301
725 268, 281, 282, 303, 302
726 269, 282, 283, 304, 303
727 270, 283, 284, 305, 304
728 271, 284, 285, 306, 305
729 272, 285, 286, 307, 306
730 273, 286, 287, 308, 307
731 274, 287, 288, 309, 308
732 275, 288, 289, 310, 309
733 276, 289, 290, 311, 310
734 277, 290, 291, 312, 311
735 278, 291, 292, 313, 312
736 279, 292, 293, 314, 313
737 280, 293, 294, 315, 314
738 281, 295, 296, 317, 316
739 282, 296, 297, 318, 317
740 283, 297, 298, 319, 318
741 284, 298, 299, 320, 319
742 285, 299, 300, 321, 320
743 286, 300, 301, 322, 321
744 287, 301, 302, 323, 322
745 288, 302, 303, 324, 323
746 289, 303, 304, 325, 324
747 290, 304, 305, 326, 325
748 291, 305, 306, 327, 326
749 292, 306, 307, 328, 327
750 293, 307, 308, 329, 328
751 294, 308, 309, 330, 329
752 295, 309, 310, 331, 330
753 296, 310, 311, 332, 331
754 297, 311, 312, 333, 332
755 298, 312, 313, 334, 333
756 299, 313, 314, 335, 334
757 300, 314, 315, 336, 335
758 301, 316, 317, 338, 337
759 302, 317, 318, 339, 338

Page 76 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

760 303, 318, 319, 340, 339

761 304, 319, 320, 341, 340
762 305, 320, 321, 342, 341
763 306, 321, 322, 343, 342
764 307, 322, 323, 344, 343
765 308, 323, 324, 345, 344
766 309, 324, 325, 346, 345
767 310, 325, 326, 347, 346
768 311, 326, 327, 348, 347
769 312, 327, 328, 349, 348
770 313, 328, 329, 350, 349
771 314, 329, 330, 351, 350
772 315, 330, 331, 352, 351
773 316, 331, 332, 353, 352
774 317, 332, 333, 354, 353
775 318, 333, 334, 355, 354
776 319, 334, 335, 356, 355
777 320, 335, 336, 357, 356
778 321, 337, 338, 359, 358
779 322, 338, 339, 360, 359
780 323, 339, 340, 361, 360
781 324, 340, 341, 362, 361
782 325, 341, 342, 363, 362
783 326, 342, 343, 364, 363
784 327, 343, 344, 365, 364
785 328, 344, 345, 366, 365
786 329, 345, 346, 367, 366
787 330, 346, 347, 368, 367
788 331, 347, 348, 369, 368
789 332, 348, 349, 370, 369
790 333, 349, 350, 371, 370
791 334, 350, 351, 372, 371
792 335, 351, 352, 373, 372
793 336, 352, 353, 374, 373
794 337, 353, 354, 375, 374
795 338, 354, 355, 376, 375
796 339, 355, 356, 377, 376
797 340, 356, 357, 378, 377
798 341, 358, 359, 380, 379
799 342, 359, 360, 381, 380
800 343, 360, 361, 382, 381
801 344, 361, 362, 383, 382
802 345, 362, 363, 384, 383
803 346, 363, 364, 385, 384
804 347, 364, 365, 386, 385
805 348, 365, 366, 387, 386
806 349, 366, 367, 388, 387
807 350, 367, 368, 389, 388
808 351, 368, 369, 390, 389
809 352, 369, 370, 391, 390
810 353, 370, 371, 392, 391
811 354, 371, 372, 393, 392
812 355, 372, 373, 394, 393
813 356, 373, 374, 395, 394
814 357, 374, 375, 396, 395
815 358, 375, 376, 397, 396
816 359, 376, 377, 398, 397
817 360, 377, 378, 399, 398
818 361, 379, 380, 401, 400
819 362, 380, 381, 402, 401
820 363, 381, 382, 403, 402
821 364, 382, 383, 404, 403
822 365, 383, 384, 405, 404
823 366, 384, 385, 406, 405
824 367, 385, 386, 407, 406
825 368, 386, 387, 408, 407
826 369, 387, 388, 409, 408
827 370, 388, 389, 410, 409
828 371, 389, 390, 411, 410

Page 77 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

829 372, 390, 391, 412, 411

830 373, 391, 392, 413, 412
831 374, 392, 393, 414, 413
832 375, 393, 394, 415, 414
833 376, 394, 395, 416, 415
834 377, 395, 396, 417, 416
835 378, 396, 397, 418, 417
836 379, 397, 398, 419, 418
837 380, 398, 399, 420, 419
838 381, 400, 401, 422, 421
839 382, 401, 402, 423, 422
840 383, 402, 403, 424, 423
841 384, 403, 404, 425, 424
842 385, 404, 405, 426, 425
843 386, 405, 406, 427, 426
844 387, 406, 407, 428, 427
845 388, 407, 408, 429, 428
846 389, 408, 409, 430, 429
847 390, 409, 410, 431, 430
848 391, 410, 411, 432, 431
849 392, 411, 412, 433, 432
850 393, 412, 413, 434, 433
851 394, 413, 414, 435, 434
852 395, 414, 415, 436, 435
853 396, 415, 416, 437, 436
854 397, 416, 417, 438, 437
855 398, 417, 418, 439, 438
856 399, 418, 419, 440, 439
857 400, 419, 420, 441, 440
858 *UEL PROPERTY, ELSET=SOLID
859 **************************************
860 *End Part
861 **
862 **
863 ** ASSEMBLY
864 **
865 *Assembly, name=Assembly
866 **
867 *Instance, name=Plate-1, part=Plate
868 *End Instance
869 **
870 *Nset, nset=all_elements, instance=Plate-1, generate
871 1, 441, 1
872 *Elset, elset=all_elements, instance=Plate-1, generate
873 1, 400, 1
874 *Nset, nset=all_nodes, instance=Plate-1, generate
875 1, 441, 1
876 *Nset, nset=encastre_nodes, instance=Plate-1
877 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16
878 17, 18, 19, 20, 21, 22, 42, 43, 63, 64, 84, 85, 105, 106, 126, 127
879 147, 148, 168, 169, 189, 190, 210, 211, 231, 232, 252, 253, 273, 274, 294, 295
880 315, 316, 336, 337, 357, 358, 378, 379, 399, 400, 420, 421, 422, 423, 424, 425
881 426, 427, 428, 429, 430, 431, 432, 433, 434, 435, 436, 437, 438, 439, 440, 441
882 *End Assembly
883 **
884 **
885 ** BOUNDARY CONDITIONS
886 **
887 ** Name: fix_nodes Type: Symmetry/Antisymmetry/Encastre
888 *Boundary
889 encastre_nodes, ENCASTRE
890 ** ----------------------------------------------------------------
891 **
892 ** STEP: loading
893 **
894 *Step, name=loading, nlgeom=NO, inc=100000
895 *Static, direct
896 0.1, 1.,
897 **

Page 78 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

898 ** LOADS
899 **
900 ** Name: Load-1 Type: Concentrated force
901 *Cload
902 all_nodes, 3, -0.01
903 **
904 ** OUTPUT REQUESTS
905 **
906 *Restart, write, frequency=0
907 **
908 ** FIELD OUTPUT: F-Output-2
909 **
910 *Output, field
911 *Node Output
912 CF, RF, RM, RT, TF, VF
913 *Element Output, directions=YES
914 BF, CENTMAG, CENTRIFMAG, CORIOMAG, ESF1, GRAV, HP, NFORC, NFORCSO, P, ROTAMAG, SF,
TRNOR, TRSHR
915 **
916 ** FIELD OUTPUT: F-Output-1
917 **
918 *Output, field, variable=PRESELECT
919 **
920 ** HISTORY OUTPUT: H-Output-2
921 **
922 *Output, history
923 *Element Output
924 IRF1, IRF2, IRF3, IRM1, IRM2, IRM3
925 **
926 ** HISTORY OUTPUT: H-Output-1
927 **
928 *Output, history, variable=PRESELECT
929 *End Step
930

Page 79 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

Abaqus job input file for default S4 element

1 *Heading
2 ** Job name: abaqus_4_elements_lift_up Model name: Model-1
3 ** Generated by: Abaqus/CAE 2016.HF2
4 *Preprint, echo=NO, model=NO, history=NO, contact=NO
5 **
6 ** PARTS
7 **
8 *Part, name=Plate
9 *Node
10 1, 0., 0., 0.
11 2, 2., 0., 0.
12 3, 0., 2., 0.
13 4, 2., 2., 0.
14 *Element, type=S4
15 1, 1, 2, 4, 3
16 *Nset, nset=_PickedSet2, internal, generate
17 1, 4, 1
18 *Elset, elset=_PickedSet2, internal
19 1,
20 ** Section: Section
21 *Shell Section, elset=_PickedSet2, material=Steel
22 0.1, 5
23 *End Part
24 **
25 **
26 ** ASSEMBLY
27 **
28 *Assembly, name=Assembly
29 **
30 *Instance, name=Plate-1, part=Plate
31 *End Instance
32 **
33 *Nset, nset=all_elements, instance=Plate-1, generate
34 1, 4, 1
35 *Elset, elset=all_elements, instance=Plate-1
36 1,
37 *Nset, nset=_PickedSet13, internal, instance=Plate-1
38 4,
39 *Nset, nset=_PickedSet14, internal, instance=Plate-1, generate
40 1, 3, 1
41 *End Assembly
42 **
43 ** MATERIALS
44 **
45 *Material, name=Steel
46 *Elastic
47 2e+08, 0.3
48 **
49 ** BOUNDARY CONDITIONS
50 **
51 ** Name: fix_nodes Type: Symmetry/Antisymmetry/Encastre
52 *Boundary
53 _PickedSet14, ENCASTRE
54 ** ----------------------------------------------------------------
55 **
56 ** STEP: loading
57 **
58 *Step, name=loading, nlgeom=NO, inc=100000
59 *Static, direct
60 0.1, 1.,
61 *****************************************************
62 *ELEMENT MATRIX
OUTPUT,ELSET=all_elements,STIFF=YES,MASS=YES,DLOAD=YES,OUTPUTFILE=USER,FILENAME=Dump
63 *****************************************************
64 ** BOUNDARY CONDITIONS
65 **
66 ** Name: lift_up Type: Displacement/Rotation
67 *Boundary
68 _PickedSet13, 3, 3, 0.1

Page 80 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

69 _PickedSet13, 4, 4
70 _PickedSet13, 5, 5
71 **
72 ** OUTPUT REQUESTS
73 **
74 *Restart, write, frequency=0
75 **
76 ** FIELD OUTPUT: F-Output-2
77 **
78 *Output, field
79 *Node Output
80 CF, RF, RM, RT, TF, VF
81 *Element Output, directions=YES
82 BF, CENTMAG, CENTRIFMAG, CORIOMAG, ESF1, GRAV, HP, NFORC, NFORCSO, P, ROTAMAG, SF,
TRNOR, TRSHR
83 **
84 ** FIELD OUTPUT: F-Output-1
85 **
86 *Output, field, variable=PRESELECT
87 **
88 ** HISTORY OUTPUT: H-Output-2
89 **
90 *Output, history
91 *Element Output
92 IRF1, IRF2, IRF3, IRM1, IRM2, IRM3
93 **
94 ** HISTORY OUTPUT: H-Output-1
95 **
96 *Output, history, variable=PRESELECT
97 *End Step
98 *STEP
99 *MATRIX GENERATE, STIFFNESS
100 *MATRIX OUTPUT, STIFFNESS, FORMAT=COORDINATE
101 *END STEP
102

Page 81 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

Modified Matlab main routine for multiple elements

1 %................................................................
2
3 % MATLAB codes for Finite Element Analysis
4 % problem19.m
5 % Mindlin plate in bending
6 % antonio ferreira 2008
7
8 % clear memory
9 clc
10 clear *;
11 clear all;
12 colordef white;
13 format long;
14 close all;
15
16 % materials
17 E = 200000000; poisson = 0.30; kapa=5/6;
18 thickness=0.1;
19 I=thickness^3/12;
20
21 % matrix C
22 % bending part
23 C_bending=I*E/(1-poisson^2)*...
24 [1 poisson 0;poisson 1 0;0 0 (1-poisson)/2];
25 % shear part
26 C_shear=kapa*thickness*E/2/(1+poisson)*eye(2);
27
28 % load
29 P = -1;
30
31 %Mesh generation
32 L = 2;
33 numberElementsX=20;
34 numberElementsY=20;
35 numberElements=numberElementsX*numberElementsY;
36 %
37 [nodeCoordinates, elementNodes] = ...
38 rectangularMesh(L,L,numberElementsX,numberElementsY);
39 xx=nodeCoordinates(:,1);
40 yy=nodeCoordinates(:,2);
41
42 %drawingMesh(nodeCoordinates,elementNodes,'Q4','k-');
43 %axis off
44 numberNodes=size(xx,1);
45
46 % GDof: global number of degrees of freedom
47 GDof=3*numberNodes;
48
49 % computation of the system stiffness matrix and force vector
50 [stiffness]=...
51 formStiffnessMatrixMindlinQ4(GDof,numberElements,...
52 elementNodes,numberNodes,nodeCoordinates,C_shear,...
53 C_bending,thickness,I);
54
55
56 [force]=...
57 formForceVectorMindlinQ4(GDof,numberElements,...
58 elementNodes,numberNodes,nodeCoordinates,P);
59
60 % % boundary conditions
61 [prescribedDof,activeDof]=...
62 EssentialBC('cccc',GDof,xx,yy,nodeCoordinates,numberNodes);
63
64 % solution
65 displacements=solution(GDof,prescribedDof,stiffness,force);
66
67 % displacements
68 disp('Displacements')
69 jj=1:GDof; format

Page 82 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

70 f=[jj; displacements'];
71 fprintf('node U\n')
72 fprintf('%3d %12.8f\n',f)
73
74 % deformed shape
75 figure
76 plot3(xx,yy,displacements(1:numberNodes),'.')
77 xlabel('X-axis') % x-axis label
78 ylabel('Y-axis') % y-axis label
79 zlabel('Z-axis') % y-axis label
80 format long
81 D1=E*thickness^3/12/(1-poisson^2);
82 min(displacements(1:numberNodes))

Page 83 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

Abaqus job input file for multiple default S4 elements

1 *Heading
2 ** Job name: push_center_Down Model name: Model-1
3 ** Generated by: Abaqus/CAE 2016.HF2
4 *Preprint, echo=NO, model=NO, history=NO, contact=NO
5 **
6 ** PARTS
7 **
8 *Part, name=Plate
9 *Node
10 1, 0., 0., 0.
11 2, 0.100000001, 0., 0.
12 3, 0.200000003, 0., 0.
13 4, 0.300000012, 0., 0.
14 5, 0.400000006, 0., 0.
15 6, 0.5, 0., 0.
16 7, 0.600000024, 0., 0.
17 8, 0.699999988, 0., 0.
18 9, 0.800000012, 0., 0.
19 10, 0.899999976, 0., 0.
20 11, 1., 0., 0.
21 12, 1.10000002, 0., 0.
22 13, 1.20000005, 0., 0.
23 14, 1.29999995, 0., 0.
24 15, 1.39999998, 0., 0.
25 16, 1.5, 0., 0.
26 17, 1.60000002, 0., 0.
27 18, 1.70000005, 0., 0.
28 19, 1.79999995, 0., 0.
29 20, 1.89999998, 0., 0.
30 21, 2., 0., 0.
31 22, 0., 0.100000001, 0.
32 23, 0.100000001, 0.100000001, 0.
33 24, 0.200000003, 0.100000001, 0.
34 25, 0.300000012, 0.100000001, 0.
35 26, 0.400000006, 0.100000001, 0.
36 27, 0.5, 0.100000001, 0.
37 28, 0.600000024, 0.100000001, 0.
38 29, 0.699999988, 0.100000001, 0.
39 30, 0.800000012, 0.100000001, 0.
40 31, 0.899999976, 0.100000001, 0.
41 32, 1., 0.100000001, 0.
42 33, 1.10000002, 0.100000001, 0.
43 34, 1.20000005, 0.100000001, 0.
44 35, 1.29999995, 0.100000001, 0.
45 36, 1.39999998, 0.100000001, 0.
46 37, 1.5, 0.100000001, 0.
47 38, 1.60000002, 0.100000001, 0.
48 39, 1.70000005, 0.100000001, 0.
49 40, 1.79999995, 0.100000001, 0.
50 41, 1.89999998, 0.100000001, 0.
51 42, 2., 0.100000001, 0.
52 43, 0., 0.200000003, 0.
53 44, 0.100000001, 0.200000003, 0.
54 45, 0.200000003, 0.200000003, 0.
55 46, 0.300000012, 0.200000003, 0.
56 47, 0.400000006, 0.200000003, 0.
57 48, 0.5, 0.200000003, 0.
58 49, 0.600000024, 0.200000003, 0.
59 50, 0.699999988, 0.200000003, 0.
60 51, 0.800000012, 0.200000003, 0.
61 52, 0.899999976, 0.200000003, 0.
62 53, 1., 0.200000003, 0.
63 54, 1.10000002, 0.200000003, 0.
64 55, 1.20000005, 0.200000003, 0.
65 56, 1.29999995, 0.200000003, 0.
66 57, 1.39999998, 0.200000003, 0.
67 58, 1.5, 0.200000003, 0.
68 59, 1.60000002, 0.200000003, 0.
69 60, 1.70000005, 0.200000003, 0.

Page 84 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

70 61, 1.79999995, 0.200000003, 0.

71 62, 1.89999998, 0.200000003, 0.
72 63, 2., 0.200000003, 0.
73 64, 0., 0.300000012, 0.
74 65, 0.100000001, 0.300000012, 0.
75 66, 0.200000003, 0.300000012, 0.
76 67, 0.300000012, 0.300000012, 0.
77 68, 0.400000006, 0.300000012, 0.
78 69, 0.5, 0.300000012, 0.
79 70, 0.600000024, 0.300000012, 0.
80 71, 0.699999988, 0.300000012, 0.
81 72, 0.800000012, 0.300000012, 0.
82 73, 0.899999976, 0.300000012, 0.
83 74, 1., 0.300000012, 0.
84 75, 1.10000002, 0.300000012, 0.
85 76, 1.20000005, 0.300000012, 0.
86 77, 1.29999995, 0.300000012, 0.
87 78, 1.39999998, 0.300000012, 0.
88 79, 1.5, 0.300000012, 0.
89 80, 1.60000002, 0.300000012, 0.
90 81, 1.70000005, 0.300000012, 0.
91 82, 1.79999995, 0.300000012, 0.
92 83, 1.89999998, 0.300000012, 0.
93 84, 2., 0.300000012, 0.
94 85, 0., 0.400000006, 0.
95 86, 0.100000001, 0.400000006, 0.
96 87, 0.200000003, 0.400000006, 0.
97 88, 0.300000012, 0.400000006, 0.
98 89, 0.400000006, 0.400000006, 0.
99 90, 0.5, 0.400000006, 0.
100 91, 0.600000024, 0.400000006, 0.
101 92, 0.699999988, 0.400000006, 0.
102 93, 0.800000012, 0.400000006, 0.
103 94, 0.899999976, 0.400000006, 0.
104 95, 1., 0.400000006, 0.
105 96, 1.10000002, 0.400000006, 0.
106 97, 1.20000005, 0.400000006, 0.
107 98, 1.29999995, 0.400000006, 0.
108 99, 1.39999998, 0.400000006, 0.
109 100, 1.5, 0.400000006, 0.
110 101, 1.60000002, 0.400000006, 0.
111 102, 1.70000005, 0.400000006, 0.
112 103, 1.79999995, 0.400000006, 0.
113 104, 1.89999998, 0.400000006, 0.
114 105, 2., 0.400000006, 0.
115 106, 0., 0.5, 0.
116 107, 0.100000001, 0.5, 0.
117 108, 0.200000003, 0.5, 0.
118 109, 0.300000012, 0.5, 0.
119 110, 0.400000006, 0.5, 0.
120 111, 0.5, 0.5, 0.
121 112, 0.600000024, 0.5, 0.
122 113, 0.699999988, 0.5, 0.
123 114, 0.800000012, 0.5, 0.
124 115, 0.899999976, 0.5, 0.
125 116, 1., 0.5, 0.
126 117, 1.10000002, 0.5, 0.
127 118, 1.20000005, 0.5, 0.
128 119, 1.29999995, 0.5, 0.
129 120, 1.39999998, 0.5, 0.
130 121, 1.5, 0.5, 0.
131 122, 1.60000002, 0.5, 0.
132 123, 1.70000005, 0.5, 0.
133 124, 1.79999995, 0.5, 0.
134 125, 1.89999998, 0.5, 0.
135 126, 2., 0.5, 0.
136 127, 0., 0.600000024, 0.
137 128, 0.100000001, 0.600000024, 0.
138 129, 0.200000003, 0.600000024, 0.

Page 85 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

139 130, 0.300000012, 0.600000024, 0.

140 131, 0.400000006, 0.600000024, 0.
141 132, 0.5, 0.600000024, 0.
142 133, 0.600000024, 0.600000024, 0.
143 134, 0.699999988, 0.600000024, 0.
144 135, 0.800000012, 0.600000024, 0.
145 136, 0.899999976, 0.600000024, 0.
146 137, 1., 0.600000024, 0.
147 138, 1.10000002, 0.600000024, 0.
148 139, 1.20000005, 0.600000024, 0.
149 140, 1.29999995, 0.600000024, 0.
150 141, 1.39999998, 0.600000024, 0.
151 142, 1.5, 0.600000024, 0.
152 143, 1.60000002, 0.600000024, 0.
153 144, 1.70000005, 0.600000024, 0.
154 145, 1.79999995, 0.600000024, 0.
155 146, 1.89999998, 0.600000024, 0.
156 147, 2., 0.600000024, 0.
157 148, 0., 0.699999988, 0.
158 149, 0.100000001, 0.699999988, 0.
159 150, 0.200000003, 0.699999988, 0.
160 151, 0.300000012, 0.699999988, 0.
161 152, 0.400000006, 0.699999988, 0.
162 153, 0.5, 0.699999988, 0.
163 154, 0.600000024, 0.699999988, 0.
164 155, 0.699999988, 0.699999988, 0.
165 156, 0.800000012, 0.699999988, 0.
166 157, 0.899999976, 0.699999988, 0.
167 158, 1., 0.699999988, 0.
168 159, 1.10000002, 0.699999988, 0.
169 160, 1.20000005, 0.699999988, 0.
170 161, 1.29999995, 0.699999988, 0.
171 162, 1.39999998, 0.699999988, 0.
172 163, 1.5, 0.699999988, 0.
173 164, 1.60000002, 0.699999988, 0.
174 165, 1.70000005, 0.699999988, 0.
175 166, 1.79999995, 0.699999988, 0.
176 167, 1.89999998, 0.699999988, 0.
177 168, 2., 0.699999988, 0.
178 169, 0., 0.800000012, 0.
179 170, 0.100000001, 0.800000012, 0.
180 171, 0.200000003, 0.800000012, 0.
181 172, 0.300000012, 0.800000012, 0.
182 173, 0.400000006, 0.800000012, 0.
183 174, 0.5, 0.800000012, 0.
184 175, 0.600000024, 0.800000012, 0.
185 176, 0.699999988, 0.800000012, 0.
186 177, 0.800000012, 0.800000012, 0.
187 178, 0.899999976, 0.800000012, 0.
188 179, 1., 0.800000012, 0.
189 180, 1.10000002, 0.800000012, 0.
190 181, 1.20000005, 0.800000012, 0.
191 182, 1.29999995, 0.800000012, 0.
192 183, 1.39999998, 0.800000012, 0.
193 184, 1.5, 0.800000012, 0.
194 185, 1.60000002, 0.800000012, 0.
195 186, 1.70000005, 0.800000012, 0.
196 187, 1.79999995, 0.800000012, 0.
197 188, 1.89999998, 0.800000012, 0.
198 189, 2., 0.800000012, 0.
199 190, 0., 0.899999976, 0.
200 191, 0.100000001, 0.899999976, 0.
201 192, 0.200000003, 0.899999976, 0.
202 193, 0.300000012, 0.899999976, 0.
203 194, 0.400000006, 0.899999976, 0.
204 195, 0.5, 0.899999976, 0.
205 196, 0.600000024, 0.899999976, 0.
206 197, 0.699999988, 0.899999976, 0.
207 198, 0.800000012, 0.899999976, 0.

Page 86 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

208 199, 0.899999976, 0.899999976, 0.

209 200, 1., 0.899999976, 0.
210 201, 1.10000002, 0.899999976, 0.
211 202, 1.20000005, 0.899999976, 0.
212 203, 1.29999995, 0.899999976, 0.
213 204, 1.39999998, 0.899999976, 0.
214 205, 1.5, 0.899999976, 0.
215 206, 1.60000002, 0.899999976, 0.
216 207, 1.70000005, 0.899999976, 0.
217 208, 1.79999995, 0.899999976, 0.
218 209, 1.89999998, 0.899999976, 0.
219 210, 2., 0.899999976, 0.
220 211, 0., 1., 0.
221 212, 0.100000001, 1., 0.
222 213, 0.200000003, 1., 0.
223 214, 0.300000012, 1., 0.
224 215, 0.400000006, 1., 0.
225 216, 0.5, 1., 0.
226 217, 0.600000024, 1., 0.
227 218, 0.699999988, 1., 0.
228 219, 0.800000012, 1., 0.
229 220, 0.899999976, 1., 0.
230 221, 1., 1., 0.
231 222, 1.10000002, 1., 0.
232 223, 1.20000005, 1., 0.
233 224, 1.29999995, 1., 0.
234 225, 1.39999998, 1., 0.
235 226, 1.5, 1., 0.
236 227, 1.60000002, 1., 0.
237 228, 1.70000005, 1., 0.
238 229, 1.79999995, 1., 0.
239 230, 1.89999998, 1., 0.
240 231, 2., 1., 0.
241 232, 0., 1.10000002, 0.
242 233, 0.100000001, 1.10000002, 0.
243 234, 0.200000003, 1.10000002, 0.
244 235, 0.300000012, 1.10000002, 0.
245 236, 0.400000006, 1.10000002, 0.
246 237, 0.5, 1.10000002, 0.
247 238, 0.600000024, 1.10000002, 0.
248 239, 0.699999988, 1.10000002, 0.
249 240, 0.800000012, 1.10000002, 0.
250 241, 0.899999976, 1.10000002, 0.
251 242, 1., 1.10000002, 0.
252 243, 1.10000002, 1.10000002, 0.
253 244, 1.20000005, 1.10000002, 0.
254 245, 1.29999995, 1.10000002, 0.
255 246, 1.39999998, 1.10000002, 0.
256 247, 1.5, 1.10000002, 0.
257 248, 1.60000002, 1.10000002, 0.
258 249, 1.70000005, 1.10000002, 0.
259 250, 1.79999995, 1.10000002, 0.
260 251, 1.89999998, 1.10000002, 0.
261 252, 2., 1.10000002, 0.
262 253, 0., 1.20000005, 0.
263 254, 0.100000001, 1.20000005, 0.
264 255, 0.200000003, 1.20000005, 0.
265 256, 0.300000012, 1.20000005, 0.
266 257, 0.400000006, 1.20000005, 0.
267 258, 0.5, 1.20000005, 0.
268 259, 0.600000024, 1.20000005, 0.
269 260, 0.699999988, 1.20000005, 0.
270 261, 0.800000012, 1.20000005, 0.
271 262, 0.899999976, 1.20000005, 0.
272 263, 1., 1.20000005, 0.
273 264, 1.10000002, 1.20000005, 0.
274 265, 1.20000005, 1.20000005, 0.
275 266, 1.29999995, 1.20000005, 0.
276 267, 1.39999998, 1.20000005, 0.

Page 87 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

277 268, 1.5, 1.20000005, 0.

278 269, 1.60000002, 1.20000005, 0.
279 270, 1.70000005, 1.20000005, 0.
280 271, 1.79999995, 1.20000005, 0.
281 272, 1.89999998, 1.20000005, 0.
282 273, 2., 1.20000005, 0.
283 274, 0., 1.29999995, 0.
284 275, 0.100000001, 1.29999995, 0.
285 276, 0.200000003, 1.29999995, 0.
286 277, 0.300000012, 1.29999995, 0.
287 278, 0.400000006, 1.29999995, 0.
288 279, 0.5, 1.29999995, 0.
289 280, 0.600000024, 1.29999995, 0.
290 281, 0.699999988, 1.29999995, 0.
291 282, 0.800000012, 1.29999995, 0.
292 283, 0.899999976, 1.29999995, 0.
293 284, 1., 1.29999995, 0.
294 285, 1.10000002, 1.29999995, 0.
295 286, 1.20000005, 1.29999995, 0.
296 287, 1.29999995, 1.29999995, 0.
297 288, 1.39999998, 1.29999995, 0.
298 289, 1.5, 1.29999995, 0.
299 290, 1.60000002, 1.29999995, 0.
300 291, 1.70000005, 1.29999995, 0.
301 292, 1.79999995, 1.29999995, 0.
302 293, 1.89999998, 1.29999995, 0.
303 294, 2., 1.29999995, 0.
304 295, 0., 1.39999998, 0.
305 296, 0.100000001, 1.39999998, 0.
306 297, 0.200000003, 1.39999998, 0.
307 298, 0.300000012, 1.39999998, 0.
308 299, 0.400000006, 1.39999998, 0.
309 300, 0.5, 1.39999998, 0.
310 301, 0.600000024, 1.39999998, 0.
311 302, 0.699999988, 1.39999998, 0.
312 303, 0.800000012, 1.39999998, 0.
313 304, 0.899999976, 1.39999998, 0.
314 305, 1., 1.39999998, 0.
315 306, 1.10000002, 1.39999998, 0.
316 307, 1.20000005, 1.39999998, 0.
317 308, 1.29999995, 1.39999998, 0.
318 309, 1.39999998, 1.39999998, 0.
319 310, 1.5, 1.39999998, 0.
320 311, 1.60000002, 1.39999998, 0.
321 312, 1.70000005, 1.39999998, 0.
322 313, 1.79999995, 1.39999998, 0.
323 314, 1.89999998, 1.39999998, 0.
324 315, 2., 1.39999998, 0.
325 316, 0., 1.5, 0.
326 317, 0.100000001, 1.5, 0.
327 318, 0.200000003, 1.5, 0.
328 319, 0.300000012, 1.5, 0.
329 320, 0.400000006, 1.5, 0.
330 321, 0.5, 1.5, 0.
331 322, 0.600000024, 1.5, 0.
332 323, 0.699999988, 1.5, 0.
333 324, 0.800000012, 1.5, 0.
334 325, 0.899999976, 1.5, 0.
335 326, 1., 1.5, 0.
336 327, 1.10000002, 1.5, 0.
337 328, 1.20000005, 1.5, 0.
338 329, 1.29999995, 1.5, 0.
339 330, 1.39999998, 1.5, 0.
340 331, 1.5, 1.5, 0.
341 332, 1.60000002, 1.5, 0.
342 333, 1.70000005, 1.5, 0.
343 334, 1.79999995, 1.5, 0.
344 335, 1.89999998, 1.5, 0.
345 336, 2., 1.5, 0.

Page 88 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

346 337, 0., 1.60000002, 0.

347 338, 0.100000001, 1.60000002, 0.
348 339, 0.200000003, 1.60000002, 0.
349 340, 0.300000012, 1.60000002, 0.
350 341, 0.400000006, 1.60000002, 0.
351 342, 0.5, 1.60000002, 0.
352 343, 0.600000024, 1.60000002, 0.
353 344, 0.699999988, 1.60000002, 0.
354 345, 0.800000012, 1.60000002, 0.
355 346, 0.899999976, 1.60000002, 0.
356 347, 1., 1.60000002, 0.
357 348, 1.10000002, 1.60000002, 0.
358 349, 1.20000005, 1.60000002, 0.
359 350, 1.29999995, 1.60000002, 0.
360 351, 1.39999998, 1.60000002, 0.
361 352, 1.5, 1.60000002, 0.
362 353, 1.60000002, 1.60000002, 0.
363 354, 1.70000005, 1.60000002, 0.
364 355, 1.79999995, 1.60000002, 0.
365 356, 1.89999998, 1.60000002, 0.
366 357, 2., 1.60000002, 0.
367 358, 0., 1.70000005, 0.
368 359, 0.100000001, 1.70000005, 0.
369 360, 0.200000003, 1.70000005, 0.
370 361, 0.300000012, 1.70000005, 0.
371 362, 0.400000006, 1.70000005, 0.
372 363, 0.5, 1.70000005, 0.
373 364, 0.600000024, 1.70000005, 0.
374 365, 0.699999988, 1.70000005, 0.
375 366, 0.800000012, 1.70000005, 0.
376 367, 0.899999976, 1.70000005, 0.
377 368, 1., 1.70000005, 0.
378 369, 1.10000002, 1.70000005, 0.
379 370, 1.20000005, 1.70000005, 0.
380 371, 1.29999995, 1.70000005, 0.
381 372, 1.39999998, 1.70000005, 0.
382 373, 1.5, 1.70000005, 0.
383 374, 1.60000002, 1.70000005, 0.
384 375, 1.70000005, 1.70000005, 0.
385 376, 1.79999995, 1.70000005, 0.
386 377, 1.89999998, 1.70000005, 0.
387 378, 2., 1.70000005, 0.
388 379, 0., 1.79999995, 0.
389 380, 0.100000001, 1.79999995, 0.
390 381, 0.200000003, 1.79999995, 0.
391 382, 0.300000012, 1.79999995, 0.
392 383, 0.400000006, 1.79999995, 0.
393 384, 0.5, 1.79999995, 0.
394 385, 0.600000024, 1.79999995, 0.
395 386, 0.699999988, 1.79999995, 0.
396 387, 0.800000012, 1.79999995, 0.
397 388, 0.899999976, 1.79999995, 0.
398 389, 1., 1.79999995, 0.
399 390, 1.10000002, 1.79999995, 0.
400 391, 1.20000005, 1.79999995, 0.
401 392, 1.29999995, 1.79999995, 0.
402 393, 1.39999998, 1.79999995, 0.
403 394, 1.5, 1.79999995, 0.
404 395, 1.60000002, 1.79999995, 0.
405 396, 1.70000005, 1.79999995, 0.
406 397, 1.79999995, 1.79999995, 0.
407 398, 1.89999998, 1.79999995, 0.
408 399, 2., 1.79999995, 0.
409 400, 0., 1.89999998, 0.
410 401, 0.100000001, 1.89999998, 0.
411 402, 0.200000003, 1.89999998, 0.
412 403, 0.300000012, 1.89999998, 0.
413 404, 0.400000006, 1.89999998, 0.
414 405, 0.5, 1.89999998, 0.

Page 89 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

415 406, 0.600000024, 1.89999998, 0.

416 407, 0.699999988, 1.89999998, 0.
417 408, 0.800000012, 1.89999998, 0.
418 409, 0.899999976, 1.89999998, 0.
419 410, 1., 1.89999998, 0.
420 411, 1.10000002, 1.89999998, 0.
421 412, 1.20000005, 1.89999998, 0.
422 413, 1.29999995, 1.89999998, 0.
423 414, 1.39999998, 1.89999998, 0.
424 415, 1.5, 1.89999998, 0.
425 416, 1.60000002, 1.89999998, 0.
426 417, 1.70000005, 1.89999998, 0.
427 418, 1.79999995, 1.89999998, 0.
428 419, 1.89999998, 1.89999998, 0.
429 420, 2., 1.89999998, 0.
430 421, 0., 2., 0.
431 422, 0.100000001, 2., 0.
432 423, 0.200000003, 2., 0.
433 424, 0.300000012, 2., 0.
434 425, 0.400000006, 2., 0.
435 426, 0.5, 2., 0.
436 427, 0.600000024, 2., 0.
437 428, 0.699999988, 2., 0.
438 429, 0.800000012, 2., 0.
439 430, 0.899999976, 2., 0.
440 431, 1., 2., 0.
441 432, 1.10000002, 2., 0.
442 433, 1.20000005, 2., 0.
443 434, 1.29999995, 2., 0.
444 435, 1.39999998, 2., 0.
445 436, 1.5, 2., 0.
446 437, 1.60000002, 2., 0.
447 438, 1.70000005, 2., 0.
448 439, 1.79999995, 2., 0.
449 440, 1.89999998, 2., 0.
450 441, 2., 2., 0.
451 *Element, type=S4R
452 1, 1, 2, 23, 22
453 2, 2, 3, 24, 23
454 3, 3, 4, 25, 24
455 4, 4, 5, 26, 25
456 5, 5, 6, 27, 26
457 6, 6, 7, 28, 27
458 7, 7, 8, 29, 28
459 8, 8, 9, 30, 29
460 9, 9, 10, 31, 30
461 10, 10, 11, 32, 31
462 11, 11, 12, 33, 32
463 12, 12, 13, 34, 33
464 13, 13, 14, 35, 34
465 14, 14, 15, 36, 35
466 15, 15, 16, 37, 36
467 16, 16, 17, 38, 37
468 17, 17, 18, 39, 38
469 18, 18, 19, 40, 39
470 19, 19, 20, 41, 40
471 20, 20, 21, 42, 41
472 21, 22, 23, 44, 43
473 22, 23, 24, 45, 44
474 23, 24, 25, 46, 45
475 24, 25, 26, 47, 46
476 25, 26, 27, 48, 47
477 26, 27, 28, 49, 48
478 27, 28, 29, 50, 49
479 28, 29, 30, 51, 50
480 29, 30, 31, 52, 51
481 30, 31, 32, 53, 52
482 31, 32, 33, 54, 53
483 32, 33, 34, 55, 54

Page 90 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

484 33, 34, 35, 56, 55

485 34, 35, 36, 57, 56
486 35, 36, 37, 58, 57
487 36, 37, 38, 59, 58
488 37, 38, 39, 60, 59
489 38, 39, 40, 61, 60
490 39, 40, 41, 62, 61
491 40, 41, 42, 63, 62
492 41, 43, 44, 65, 64
493 42, 44, 45, 66, 65
494 43, 45, 46, 67, 66
495 44, 46, 47, 68, 67
496 45, 47, 48, 69, 68
497 46, 48, 49, 70, 69
498 47, 49, 50, 71, 70
499 48, 50, 51, 72, 71
500 49, 51, 52, 73, 72
501 50, 52, 53, 74, 73
502 51, 53, 54, 75, 74
503 52, 54, 55, 76, 75
504 53, 55, 56, 77, 76
505 54, 56, 57, 78, 77
506 55, 57, 58, 79, 78
507 56, 58, 59, 80, 79
508 57, 59, 60, 81, 80
509 58, 60, 61, 82, 81
510 59, 61, 62, 83, 82
511 60, 62, 63, 84, 83
512 61, 64, 65, 86, 85
513 62, 65, 66, 87, 86
514 63, 66, 67, 88, 87
515 64, 67, 68, 89, 88
516 65, 68, 69, 90, 89
517 66, 69, 70, 91, 90
518 67, 70, 71, 92, 91
519 68, 71, 72, 93, 92
520 69, 72, 73, 94, 93
521 70, 73, 74, 95, 94
522 71, 74, 75, 96, 95
523 72, 75, 76, 97, 96
524 73, 76, 77, 98, 97
525 74, 77, 78, 99, 98
526 75, 78, 79, 100, 99
527 76, 79, 80, 101, 100
528 77, 80, 81, 102, 101
529 78, 81, 82, 103, 102
530 79, 82, 83, 104, 103
531 80, 83, 84, 105, 104
532 81, 85, 86, 107, 106
533 82, 86, 87, 108, 107
534 83, 87, 88, 109, 108
535 84, 88, 89, 110, 109
536 85, 89, 90, 111, 110
537 86, 90, 91, 112, 111
538 87, 91, 92, 113, 112
539 88, 92, 93, 114, 113
540 89, 93, 94, 115, 114
541 90, 94, 95, 116, 115
542 91, 95, 96, 117, 116
543 92, 96, 97, 118, 117
544 93, 97, 98, 119, 118
545 94, 98, 99, 120, 119
546 95, 99, 100, 121, 120
547 96, 100, 101, 122, 121
548 97, 101, 102, 123, 122
549 98, 102, 103, 124, 123
550 99, 103, 104, 125, 124
551 100, 104, 105, 126, 125
552 101, 106, 107, 128, 127

Page 91 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

553 102, 107, 108, 129, 128

554 103, 108, 109, 130, 129
555 104, 109, 110, 131, 130
556 105, 110, 111, 132, 131
557 106, 111, 112, 133, 132
558 107, 112, 113, 134, 133
559 108, 113, 114, 135, 134
560 109, 114, 115, 136, 135
561 110, 115, 116, 137, 136
562 111, 116, 117, 138, 137
563 112, 117, 118, 139, 138
564 113, 118, 119, 140, 139
565 114, 119, 120, 141, 140
566 115, 120, 121, 142, 141
567 116, 121, 122, 143, 142
568 117, 122, 123, 144, 143
569 118, 123, 124, 145, 144
570 119, 124, 125, 146, 145
571 120, 125, 126, 147, 146
572 121, 127, 128, 149, 148
573 122, 128, 129, 150, 149
574 123, 129, 130, 151, 150
575 124, 130, 131, 152, 151
576 125, 131, 132, 153, 152
577 126, 132, 133, 154, 153
578 127, 133, 134, 155, 154
579 128, 134, 135, 156, 155
580 129, 135, 136, 157, 156
581 130, 136, 137, 158, 157
582 131, 137, 138, 159, 158
583 132, 138, 139, 160, 159
584 133, 139, 140, 161, 160
585 134, 140, 141, 162, 161
586 135, 141, 142, 163, 162
587 136, 142, 143, 164, 163
588 137, 143, 144, 165, 164
589 138, 144, 145, 166, 165
590 139, 145, 146, 167, 166
591 140, 146, 147, 168, 167
592 141, 148, 149, 170, 169
593 142, 149, 150, 171, 170
594 143, 150, 151, 172, 171
595 144, 151, 152, 173, 172
596 145, 152, 153, 174, 173
597 146, 153, 154, 175, 174
598 147, 154, 155, 176, 175
599 148, 155, 156, 177, 176
600 149, 156, 157, 178, 177
601 150, 157, 158, 179, 178
602 151, 158, 159, 180, 179
603 152, 159, 160, 181, 180
604 153, 160, 161, 182, 181
605 154, 161, 162, 183, 182
606 155, 162, 163, 184, 183
607 156, 163, 164, 185, 184
608 157, 164, 165, 186, 185
609 158, 165, 166, 187, 186
610 159, 166, 167, 188, 187
611 160, 167, 168, 189, 188
612 161, 169, 170, 191, 190
613 162, 170, 171, 192, 191
614 163, 171, 172, 193, 192
615 164, 172, 173, 194, 193
616 165, 173, 174, 195, 194
617 166, 174, 175, 196, 195
618 167, 175, 176, 197, 196
619 168, 176, 177, 198, 197
620 169, 177, 178, 199, 198
621 170, 178, 179, 200, 199

Page 92 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

622 171, 179, 180, 201, 200

623 172, 180, 181, 202, 201
624 173, 181, 182, 203, 202
625 174, 182, 183, 204, 203
626 175, 183, 184, 205, 204
627 176, 184, 185, 206, 205
628 177, 185, 186, 207, 206
629 178, 186, 187, 208, 207
630 179, 187, 188, 209, 208
631 180, 188, 189, 210, 209
632 181, 190, 191, 212, 211
633 182, 191, 192, 213, 212
634 183, 192, 193, 214, 213
635 184, 193, 194, 215, 214
636 185, 194, 195, 216, 215
637 186, 195, 196, 217, 216
638 187, 196, 197, 218, 217
639 188, 197, 198, 219, 218
640 189, 198, 199, 220, 219
641 190, 199, 200, 221, 220
642 191, 200, 201, 222, 221
643 192, 201, 202, 223, 222
644 193, 202, 203, 224, 223
645 194, 203, 204, 225, 224
646 195, 204, 205, 226, 225
647 196, 205, 206, 227, 226
648 197, 206, 207, 228, 227
649 198, 207, 208, 229, 228
650 199, 208, 209, 230, 229
651 200, 209, 210, 231, 230
652 201, 211, 212, 233, 232
653 202, 212, 213, 234, 233
654 203, 213, 214, 235, 234
655 204, 214, 215, 236, 235
656 205, 215, 216, 237, 236
657 206, 216, 217, 238, 237
658 207, 217, 218, 239, 238
659 208, 218, 219, 240, 239
660 209, 219, 220, 241, 240
661 210, 220, 221, 242, 241
662 211, 221, 222, 243, 242
663 212, 222, 223, 244, 243
664 213, 223, 224, 245, 244
665 214, 224, 225, 246, 245
666 215, 225, 226, 247, 246
667 216, 226, 227, 248, 247
668 217, 227, 228, 249, 248
669 218, 228, 229, 250, 249
670 219, 229, 230, 251, 250
671 220, 230, 231, 252, 251
672 221, 232, 233, 254, 253
673 222, 233, 234, 255, 254
674 223, 234, 235, 256, 255
675 224, 235, 236, 257, 256
676 225, 236, 237, 258, 257
677 226, 237, 238, 259, 258
678 227, 238, 239, 260, 259
679 228, 239, 240, 261, 260
680 229, 240, 241, 262, 261
681 230, 241, 242, 263, 262
682 231, 242, 243, 264, 263
683 232, 243, 244, 265, 264
684 233, 244, 245, 266, 265
685 234, 245, 246, 267, 266
686 235, 246, 247, 268, 267
687 236, 247, 248, 269, 268
688 237, 248, 249, 270, 269
689 238, 249, 250, 271, 270
690 239, 250, 251, 272, 271

Page 93 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

691 240, 251, 252, 273, 272

692 241, 253, 254, 275, 274
693 242, 254, 255, 276, 275
694 243, 255, 256, 277, 276
695 244, 256, 257, 278, 277
696 245, 257, 258, 279, 278
697 246, 258, 259, 280, 279
698 247, 259, 260, 281, 280
699 248, 260, 261, 282, 281
700 249, 261, 262, 283, 282
701 250, 262, 263, 284, 283
702 251, 263, 264, 285, 284
703 252, 264, 265, 286, 285
704 253, 265, 266, 287, 286
705 254, 266, 267, 288, 287
706 255, 267, 268, 289, 288
707 256, 268, 269, 290, 289
708 257, 269, 270, 291, 290
709 258, 270, 271, 292, 291
710 259, 271, 272, 293, 292
711 260, 272, 273, 294, 293
712 261, 274, 275, 296, 295
713 262, 275, 276, 297, 296
714 263, 276, 277, 298, 297
715 264, 277, 278, 299, 298
716 265, 278, 279, 300, 299
717 266, 279, 280, 301, 300
718 267, 280, 281, 302, 301
719 268, 281, 282, 303, 302
720 269, 282, 283, 304, 303
721 270, 283, 284, 305, 304
722 271, 284, 285, 306, 305
723 272, 285, 286, 307, 306
724 273, 286, 287, 308, 307
725 274, 287, 288, 309, 308
726 275, 288, 289, 310, 309
727 276, 289, 290, 311, 310
728 277, 290, 291, 312, 311
729 278, 291, 292, 313, 312
730 279, 292, 293, 314, 313
731 280, 293, 294, 315, 314
732 281, 295, 296, 317, 316
733 282, 296, 297, 318, 317
734 283, 297, 298, 319, 318
735 284, 298, 299, 320, 319
736 285, 299, 300, 321, 320
737 286, 300, 301, 322, 321
738 287, 301, 302, 323, 322
739 288, 302, 303, 324, 323
740 289, 303, 304, 325, 324
741 290, 304, 305, 326, 325
742 291, 305, 306, 327, 326
743 292, 306, 307, 328, 327
744 293, 307, 308, 329, 328
745 294, 308, 309, 330, 329
746 295, 309, 310, 331, 330
747 296, 310, 311, 332, 331
748 297, 311, 312, 333, 332
749 298, 312, 313, 334, 333
750 299, 313, 314, 335, 334
751 300, 314, 315, 336, 335
752 301, 316, 317, 338, 337
753 302, 317, 318, 339, 338
754 303, 318, 319, 340, 339
755 304, 319, 320, 341, 340
756 305, 320, 321, 342, 341
757 306, 321, 322, 343, 342
758 307, 322, 323, 344, 343
759 308, 323, 324, 345, 344

Page 94 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

760 309, 324, 325, 346, 345

761 310, 325, 326, 347, 346
762 311, 326, 327, 348, 347
763 312, 327, 328, 349, 348
764 313, 328, 329, 350, 349
765 314, 329, 330, 351, 350
766 315, 330, 331, 352, 351
767 316, 331, 332, 353, 352
768 317, 332, 333, 354, 353
769 318, 333, 334, 355, 354
770 319, 334, 335, 356, 355
771 320, 335, 336, 357, 356
772 321, 337, 338, 359, 358
773 322, 338, 339, 360, 359
774 323, 339, 340, 361, 360
775 324, 340, 341, 362, 361
776 325, 341, 342, 363, 362
777 326, 342, 343, 364, 363
778 327, 343, 344, 365, 364
779 328, 344, 345, 366, 365
780 329, 345, 346, 367, 366
781 330, 346, 347, 368, 367
782 331, 347, 348, 369, 368
783 332, 348, 349, 370, 369
784 333, 349, 350, 371, 370
785 334, 350, 351, 372, 371
786 335, 351, 352, 373, 372
787 336, 352, 353, 374, 373
788 337, 353, 354, 375, 374
789 338, 354, 355, 376, 375
790 339, 355, 356, 377, 376
791 340, 356, 357, 378, 377
792 341, 358, 359, 380, 379
793 342, 359, 360, 381, 380
794 343, 360, 361, 382, 381
795 344, 361, 362, 383, 382
796 345, 362, 363, 384, 383
797 346, 363, 364, 385, 384
798 347, 364, 365, 386, 385
799 348, 365, 366, 387, 386
800 349, 366, 367, 388, 387
801 350, 367, 368, 389, 388
802 351, 368, 369, 390, 389
803 352, 369, 370, 391, 390
804 353, 370, 371, 392, 391
805 354, 371, 372, 393, 392
806 355, 372, 373, 394, 393
807 356, 373, 374, 395, 394
808 357, 374, 375, 396, 395
809 358, 375, 376, 397, 396
810 359, 376, 377, 398, 397
811 360, 377, 378, 399, 398
812 361, 379, 380, 401, 400
813 362, 380, 381, 402, 401
814 363, 381, 382, 403, 402
815 364, 382, 383, 404, 403
816 365, 383, 384, 405, 404
817 366, 384, 385, 406, 405
818 367, 385, 386, 407, 406
819 368, 386, 387, 408, 407
820 369, 387, 388, 409, 408
821 370, 388, 389, 410, 409
822 371, 389, 390, 411, 410
823 372, 390, 391, 412, 411
824 373, 391, 392, 413, 412
825 374, 392, 393, 414, 413
826 375, 393, 394, 415, 414
827 376, 394, 395, 416, 415
828 377, 395, 396, 417, 416

Page 95 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

829 378, 396, 397, 418, 417

830 379, 397, 398, 419, 418
831 380, 398, 399, 420, 419
832 381, 400, 401, 422, 421
833 382, 401, 402, 423, 422
834 383, 402, 403, 424, 423
835 384, 403, 404, 425, 424
836 385, 404, 405, 426, 425
837 386, 405, 406, 427, 426
838 387, 406, 407, 428, 427
839 388, 407, 408, 429, 428
840 389, 408, 409, 430, 429
841 390, 409, 410, 431, 430
842 391, 410, 411, 432, 431
843 392, 411, 412, 433, 432
844 393, 412, 413, 434, 433
845 394, 413, 414, 435, 434
846 395, 414, 415, 436, 435
847 396, 415, 416, 437, 436
848 397, 416, 417, 438, 437
849 398, 417, 418, 439, 438
850 399, 418, 419, 440, 439
851 400, 419, 420, 441, 440
852 *Nset, nset=_PickedSet2, internal, generate
853 1, 441, 1
854 *Elset, elset=_PickedSet2, internal, generate
855 1, 400, 1
856 ** Section: Section
857 *Shell Section, elset=_PickedSet2, material=Steel
858 0.1, 5
859 *End Part
860 **
861 **
862 ** ASSEMBLY
863 **
864 *Assembly, name=Assembly
865 **
866 *Instance, name=Plate-1, part=Plate
867 *End Instance
868 **
869 *Nset, nset=all_elements, instance=Plate-1, generate
870 1, 441, 1
871 *Elset, elset=all_elements, instance=Plate-1, generate
872 1, 400, 1
873 *Nset, nset=encastre_nodes, instance=Plate-1
874 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16
875 17, 18, 19, 20, 21, 22, 42, 43, 63, 64, 84, 85, 105, 106, 126, 127
876 147, 148, 168, 169, 189, 190, 210, 211, 231, 232, 252, 253, 273, 274, 294, 295
877 315, 316, 336, 337, 357, 358, 378, 379, 399, 400, 420, 421, 422, 423, 424, 425
878 426, 427, 428, 429, 430, 431, 432, 433, 434, 435, 436, 437, 438, 439, 440, 441
879 *Elset, elset=__PickedSurf13_SPOS, internal, instance=Plate-1, generate
880 1, 400, 1
881 *Surface, type=ELEMENT, name=_PickedSurf13, internal
882 __PickedSurf13_SPOS, SPOS
883 *End Assembly
884 **
885 ** MATERIALS
886 **
887 *Material, name=Steel
888 *Elastic
889 2e+08, 0.3
890 **
891 ** BOUNDARY CONDITIONS
892 **
893 ** Name: fix_nodes Type: Symmetry/Antisymmetry/Encastre
894 *Boundary
895 encastre_nodes, ENCASTRE
896 ** ----------------------------------------------------------------
897 **

Page 96 of 97
Nikhil Padhye, Subodh Kalia Implementation of Mindlin plate element

898 ** STEP: loading

899 **
900 *Step, name=loading, nlgeom=NO, inc=100000
901 *Static
902 0.0001, 1., 1e-08, 0.01
903 **
904 ** LOADS
905 **
906 ** Name: Load-1 Type: Pressure
907 *Dsload
908 _PickedSurf13, P, 1.
909 **
910 ** OUTPUT REQUESTS
911 **
912 *Restart, write, frequency=0
913 **
914 ** FIELD OUTPUT: F-Output-2
915 **
916 *Output, field
917 *Node Output
918 CF, RF, RM, RT, TF, VF
919 *Element Output, directions=YES
920 BF, CENTMAG, CENTRIFMAG, CORIOMAG, ESF1, GRAV, HP, NFORC, NFORCSO, P, ROTAMAG, SF,
TRNOR, TRSHR
921 **
922 ** FIELD OUTPUT: F-Output-1
923 **
924 *Output, field, variable=PRESELECT
925 **
926 ** HISTORY OUTPUT: H-Output-2
927 **
928 *Output, history
929 *Element Output
930 IRF1, IRF2, IRF3, IRM1, IRM2, IRM3
931 **
932 ** HISTORY OUTPUT: H-Output-1
933 **
934 *Output, history, variable=PRESELECT
935 *End Step
936

Page 97 of 97