Civil Engineering Reference
In-Depth Information
L
E=10
5
kN/m
2
u
= 0.3
1 kN/m
2
1
41
0 m
21
26
31
36
6
2
22
7
3
4
5
-2 m
45
0 m
1 m
6 m
type_2d
'plane'
element nod dir
'triangle' 15 'y'
nxe nye nip np_types
2 1 12 1
prop(e,v)
1.0e5 0.2
etype(not needed)
x_coords, y_coords
0.0 1.0 6.0
0.0 -2.0
nr,(k,nf(:,k),i=1,nr)
17
1 0 1 2 0 1 3 0 1 4 0 1 5 0 0
10 0 0 15 0 0 20 0 0 25 0 0 30 0 0
35 0 0 40 0 0 41 0 1 42 0 1 43 0 1
44 0 1 45 0 0
loaded_nodes,(k,loads(nf(:,k)),i=1,loaded_nodes)
5
1 0.0 -0.0778 6 0.0 -0.3556 11 0.0 -0.1333
16 0.0 -0.3556 21 0.0 -0.0778
fixed_freedoms
0
Figure 5.7 Mesh and data for second Program 5.1 example
a roller boundary at six times the load width from the centreline. The data indicate that a
15-noded triangle is to be used in a plane strain analysis, with node and element numbering
in the
y
-direction. The relatively high order of the interpolating polynomials associated with
this element, suggests that fewer elements would be required for a typical boundary value
problem than if working with a lower order element. The mesh shown in Figure 5.7 consists
of two columns of elements (
nxe=2
) and one row of elements (
nye=1
). The recommended
