Civil Engineering Reference
In-Depth Information
L
0.25 kN
0.5 kN
0.25 kN
0,1
2,3
4,5
1
1
2
2
3
3
E=10
6
kN/m
2
υ
= 0.3
0,6
7,8
9,10
1 m
4
5
6
7
8
9
0,0
11,0
12,0
type_2d
'plane'
element nod dir
'triangle' 3 'x'
nxe nye nip np_types
2 2 1 1
prop(e,v)
1.0e6 0.3
etype(not needed)
x_coords, y_coords
0.0 0.5 1.0
0.0 -0.5 -1.0
nr,(k,nf(:,k),i=1,nr)
1 0 1 4 0 1 7 0 0 8 1 0 9 1 0
loaded_nodes,(k,loads(nf(:,k)),i=1,loaded_nodes)
1 0.0 -0.25 2 0.0 -0.50 3 0.0 -0.25
fixed_freedoms
0
1m
Figure 5.2 Mesh and data for first Program 5.1 example
nye=2
, which sets the mesh to consist of two columns and two rows of elements respec-
tively, with the diagonals forming triangles as referred to above,
nip=1
which sets the
numerical integration to use one integrating point per element and
np types
which indi-
cates that there is only one property group in this homogeneous example. As usual, since
np types
equals 1, the
etype
data is not required. The next line reads the properties
which in an elastic analysis consist of Young's modulus and Poisson's ratio, set respec-
tively to 1
10
6
kN/m
2
and 0
.
3. The next two lines read the
x
-coordinates (
x coords
)
and
y
-coordinates (
y coords
) of the vertical and horizontal lines that make up the mesh.
Nodal freedom data concerning boundary restraints is read next, consisting of the number
of restrained nodes,
nr=5
, followed by the restrained node number and a binary “on-off”
switch corresponding to the
x
-and
y
-displacement components. For example
101
means
×
