You do not have permission to edit this page, for the following reason:

You are not allowed to execute the action you have requested.


You can view and copy the source of this page.

x
 
1
2
3
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
4
<big>'''ANALYSIS OF FLUID-SOIL-STRUCTURE INTERACTION PROBLEMS WITH THE PARTICLE FINITE ELEMENT METHOD (PFEM)'''</big></div>
5
6
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
7
'''E. Oñate, '''<sup>1</sup>;''' M.A. Celigueta, '''<sup>1</sup>;''' S.R. Idelsohn, '''<sup>2</sup>;''' F. Salazar,'''<sup>1</sup>;''' A. Larese, '''<sup>1</sup>;''' R. Rossi, '''<sup>1</sup>; '''B. Suárez, '''<sup>1</sup>; '''R. Morán, '''<sup>3</sup></div>
8
9
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
10
<sup>1</sup> Centre Internacional de Metodes Numerics a l'Enginyeria (CIMNE)</div>
11
12
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
13
Gran Capitán s/n, 08034 Barcelona, Spain, [http://www.cimne.com www.cimne.com]</div>
14
15
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
16
<sup>2</sup> ICREA Research Professor at CIMNE</div>
17
18
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
19
<sup>3</sup> School of Civil Engineering. Universidad Politécnica de Madrid, Spain</div>
20
21
<span id='_GoBack'></span>
22
23
==Summary ==
24
25
We present some developments in the Particle Finite Element Method (PFEM) for analysis of coupled problems in mechanics involving fluid-soil-structure interaction (FSSI). The PFEM uses an updated Lagrangian description to model the motion of material points in both the fluid and the solid domains (the later including soil/rocks and structures). A mesh connects the particles (nodes) defining the discretized domain where the governing equations for each of the constituent materials are solved as in the standard FEM. The procedure to model frictional contact conditions and material erosion at fluid-solid and solid-solid interfaces is described. We present several examples of application of the PFEM to solve FSSI problems such as the motion of rocks by water streams, the erosion of river beds, the stability of breakwaters and constructions under sea waves, the falling of landslides on houses and into a reservoir and the failure of rockfill dams in overspill situations.
26
27
<big>'''<br/>1. Introduction'''</big>
28
29
The analysis of problems involving the interaction of fluids, soil/rocks and structures is of relevance in many areas of engineering. Examples are common in the study of landslides and their effect on reservoirs and adjacent structures, off-shore and harbour structures under large waves, constructions hit by floods and tsunamis, soil erosion and stability of rock-fill dams in overspill situations, etc.
30
31
The authors have developed in previous works a particular class of Lagrangian formulation for solving problems involving the interaction between free surface fluids and solids. The so-called particle finite element method (PFEM, [http://www.cimne.com/pfem www.cimne.com/pfem]), treats the mesh nodes in the fluid and solid domains as particles which can freely move and even separate from the main fluid domain representing, for instance, the effect of water drops. A mesh connects the nodes discretizing the domain where the governing equations are solved using a stabilized FEM.
32
33
An advantage of the PFEM is that, as a Lagrangian formulation, the convective terms disappear from the fluid equations [4,25]. The difficulty is however transferred to the problem of adequately (and efficiently) moving the mesh nodes. We use a mesh regeneration procedure blending elements of different shapes using an extended Delaunay tesselation with special shape functions [8]. The theory and applications of the PFEM are reported in [2,3,9-22].
34
35
In this paper we describe recent advances of the PFEM for fluid-soil-structure interaction (FSSI) problems and present several examples of application to the study of the motion of rocks by water streams, the stability of breakwaters and constructions under sea waves, the falling of landslides on houses and into a reservoir and the failure of a rockfill dam in an overspill situation.
36
37
<big>'''2. The Basis of the Particle Finite Element method'''</big>
38
39
In the PFEM both the fluid and the solid domains are modelled using an '' updated'' '' Lagrangian formulation'' [9,11,16,25]. That is, all variables are assumed to be known in the '' current configuration'' at time [[Image:draft_Samper_231355424-picture-x0000_i1025.svg|9px]] . The new set of variables in both domains are sought for in the '' next or updated configuration'' at time [[Image:draft_Samper_231355424-picture-x0000_i1026.svg|39px]] . The finite element method (FEM) is used to solve the equations of continuum mechanics for each of the subdomains. Hence a mesh discretizing these domains must be generated in order to solve the governing equations for each subdomain in the standard FEM fashion. The quality of the numerical solution depends on the discretization chosen as in the standard FEM. Adaptive mesh refinement techniques can be used to improve the solution in zones where large motions of the fluid or the structure occur.
40
41
A typical solution with the PFEM involves the following steps.
42
43
1.  The starting point at each time step is the cloud of points (particles) in the fluid and solid domains. For instance [[Image:draft_Samper_231355424-picture-x0000_i1027.svg|22px]] denotes the cloud at time [[Image:draft_Samper_231355424-picture-x0000_i1028.svg|36px]] (Fig. 1).
44
45
2.  Identify the boundaries for both the fluid and solid domains defining the analysis domain [[Image:draft_Samper_231355424-picture-x0000_i1029.svg|22px]] in the fluid and the solid. This is an essential step as some boundaries (such as the free surface in fluids) may be severely distorted during the solution, including separation and re-entering of nodes. The Alpha Shape method [6,9,17] is used for the boundary definition.
46
47
3.  Discretize the fluid and solid domains with a finite element mesh [[Image:draft_Samper_231355424-picture-x0000_i1030.svg|28px]] . We use an innovative mesh generation scheme based on the extended Delaunay tesselation [8].
48
49
4.  Solve the coupled Lagrangian equations of motion for the fluid and the solid domains. Compute the state variables at the next (updated) configuration for [[Image:draft_Samper_231355424-picture-x0000_i1031.svg|39px]] : velocities, pressure and viscous stresses in the fluid and displacements, stresses and strains in the solid.
50
51
5.  Move the mesh nodes to a new position [[Image:draft_Samper_231355424-picture-x0000_i1032.svg|31px]] where [[Image:draft_Samper_231355424-picture-x0000_i1033.svg|31px]] denotes the time [[Image:draft_Samper_231355424-picture-x0000_i1034.svg|45px]] , in terms of the time increment size.
52
53
6.  Go back to step 1 and repeat the process for the next time step to obtain [[Image:draft_Samper_231355424-picture-x0000_i1035.svg|33px]] (Fig. 1).
54
55
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
56
'''
57
[[Image:draft_Samper_231355424-image12.png|center|480px]]
58
'''</div>
59
60
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
61
Figure 1: Sequence of steps to update a ``cloud´´ of nodes representing a domain containing fluid and a solid subdomains from time [[Image:draft_Samper_231355424-picture-x0000_i1036.svg|13px]] ([[Image:draft_Samper_231355424-picture-x0000_i1037.svg|36px]] ) to time [[Image:draft_Samper_231355424-picture-x0000_i1038.svg|36px]] ([[Image:draft_Samper_231355424-picture-x0000_i1039.svg|74px]] )</div>
62
63
<big>'''3. FIC/FEM formulation for a Lagrangian continuum'''</big>
64
65
==3.1. Governing equations==
66
67
The equations to be solved are the standard ones in continuum mechanics, written in the Lagrangian frame of reference:
68
69
{|
70
|-
71
| style="width: 60.0075mm;vertical-align: top;" | '''Momentum: '''
72
| style="" | [[Image:draft_Samper_231355424-picture-x0000_i1040.svg|169px]] 
73
| style="text-align: right; margin-left: 1em;" |  (1)
74
|}
75
76
{|
77
|-
78
| style="width: 60.0075mm;vertical-align: top;" | '''Pressure-velocity relationship:     '''[[Image:draft_Samper_231355424-picture-x0000_i1041.svg|161px]] 
79
| style="" | (2)
80
| style="text-align: right; margin-left: 1em;" | 
81
| style="" | 
82
| style="" | 
83
| style="" | 
84
|}
85
86
In above equations [[Image:draft_Samper_231355424-picture-x0000_i1042.svg|14px]] is the velocity along the ''i''th global (cartesian) axis, [[Image:draft_Samper_231355424-picture-x0000_i1043.svg|16px]] is the pressure (assumed to be positive in compression) [[Image:draft_Samper_231355424-picture-x0000_i1044.svg|16px]] and [[Image:draft_Samper_231355424-picture-x0000_i1045.svg|17px]] are the density and bulk modulus of the material, respectively, [[Image:draft_Samper_231355424-picture-x0000_i1046.svg|14px]] and [[Image:draft_Samper_231355424-picture-x0000_i1047.svg|21px]] are the body forces and the (Cauchy) stresses. Equations 1 and 2 are completed with the constitutive relationships. For details see [11,24,25].
87
88
==3.2 Discretization of the equations==
89
90
A key problem in the solution of Equations 1--2 is the satisfaction of the mass balance condition for the incompressible case (i.e. [[Image:draft_Samper_231355424-picture-x0000_i1048.svg|45px]] ). A number of procedures to solve his problem exist in the finite element literature [4,25]. In our approach we use a stabilized formulation based in the so-called finite calculus procedure [15]. The essence of this method is the solution of a ''modified mass balance'' equation which is written as
91
92
{|
93
|-
94
| style="width: 77.5229mm;vertical-align: top;" | 
95
| style="margin: 1em auto 1em auto;" | [[Image:draft_Samper_231355424-picture-x0000_i1049.svg|3px]] 
96
| style="text-align: right; margin-left: 1em;" | (5)
97
|}
98
99
where [[Image:draft_Samper_231355424-picture-x0000_i1050.svg|13px]] are weighting functions, [[Image:draft_Samper_231355424-picture-x0000_i1051.svg|13px]] is a stabilization parameter. 6 [[Image:draft_Samper_231355424-picture-x0000_i1052.svg|17px]] are auxiliary pressure projection variables chosen so as to ensure that the second term in Equation 6 can be interpreted as weighted sum of the residuals of the momentum equations and therefore it vanishes for the exact solution. The set of governing equations is completed by adding the following constraint equation
100
101
{|
102
|-
103
| style="width: 72.5135mm;vertical-align: top;" | 
104
| style="margin: 1em auto 1em auto;" | [[Image:draft_Samper_231355424-picture-x0000_i1053.svg|600px]] 
105
| style="text-align: right; margin-left: 1em;" | (6)
106
|}
107
108
where [[Image:draft_Samper_231355424-picture-x0000_i1054.svg|19px]] are arbitrary weighting functions.
109
110
The rest of the integral equations are obtained by applying the standard weighted residual technique to the governing equations (1), (2), (3) and (5) and the corresponding boundary conditions [11,18,21]. We interpolate next in the standard finite element fashion (for 3D problems) the three velocities [[Image:draft_Samper_231355424-picture-x0000_i1055.svg|14px]] , the pressure [[Image:draft_Samper_231355424-picture-x0000_i1056.svg|19px]] and the three pressure gradient projections [[Image:draft_Samper_231355424-picture-x0000_i1057.svg|17px]] . In our work we use equal order ''linear interpolation'' for all variables over meshes of 3-noded triangles (in 2D) and 4-noded tetrahedra (in 3D). The resulting set of discretized equations using the standard Galerkin technique has the following form
111
112
{|
113
|-
114
| style="width: 84.6667mm;vertical-align: top;" | '''Momentum:'''
115
| style="margin: 1em auto 1em auto;" | [[Image:draft_Samper_231355424-picture-x0000_i1058.svg|122px]] 
116
| style="text-align: right; margin-left: 1em;" | (7)
117
|}
118
119
{|
120
|-
121
| style="width: 84.6667mm;vertical-align: top;" | '''Pressure-velocity relationship:'''
122
| style="margin: 1em auto 1em auto;" | [[Image:draft_Samper_231355424-picture-x0000_i1059.svg|154px]] 
123
| style="text-align: right; margin-left: 1em;" | (8)
124
|}
125
126
{|
127
|-
128
| style="width: 77.5229mm;vertical-align: top;" | '''Pressure gradient projection:'''
129
| style="margin: 1em auto 1em auto;" |   [[Image:draft_Samper_231355424-picture-x0000_i1060.svg|96px]] 
130
| style="text-align: right; margin-left: 1em;" | (9)
131
|}
132
133
In Eqs.(7)-(9)  [[Image:draft_Samper_231355424-picture-x0000_i1061.svg|22px]] denotes nodal variables, [[Image:draft_Samper_231355424-picture-x0000_i1062.svg|73px]] .
134
135
The solution of Eqs.(7)-(9) can be performed using any time integration scheme typical of the updated Lagrangian FEM. Details of the element matrices and the solution procedure can be found in [11,18,21].
136
137
<big>'''4. Generation of a new mesh'''</big>
138
139
One of the key points for the success of the PFEM is the fast regeneration of a mesh at every time step on the basis of the position of the nodes in the space domain. Indeed, any fast meshing algorithm can be used for this purpose. In our work the mesh is generated at each time step using the so called extended Delaunay tesselation (EDT) presented in [8].
140
141
As a general rule for large 3D problems meshing consumes around 15% of the total CPU time for each time step, while the solution of the equations (with typically 3 iterations to reach convergence within a time step) and the assembling of the system consume approximately 70% and 15% of the CPU time for each time step, respectively. These figures refer to solutions obtained in a standard single processor Pentium IV PC for all the computations and prove that the generation of the mesh has an acceptable cost in the PFEM [21].
142
143
<big>'''5. Identification of boundary surfaces'''</big>
144
145
One of the main tasks in the PFEM is the correct definition of the boundary domain. Boundary nodes are sometimes explicitly identified. In other cases, the total set of nodes is the only information available and the algorithm must recognize the boundary nodes.
146
147
In our work we use an extended Delaunay partition for recognizing boundary nodes [8,16]. Considering that the nodes follow a variable distribution, where is typically the minimum distance between two nodes. All nodes on an empty sphere with a radius greater than are considered as boundary nodes. In practice is a parameter close to, but greater than one. Values of ranging between 1.3 and 1.5 have been found to be optimal in all examples analyzed. This criterion is coincident with the Alpha Shape concept [6].
148
149
The boundary recognition method is also useful for detecting isolated fluid particle outside the main fluid domain and contact conditions between the fluid domain and a fixed boundary, as well as between different interacting solids.
150
151
<big>'''6. Treatment of contact conditions in the PFEM'''</big>
152
153
==6.1 Contact between the fluid and a fixed boundary==
154
155
The condition of prescribed velocities at the fixed boundaries in the PFEM are applied in strong form to the boundary nodes. Nodes might belong to fixed external boundaries or to moving boundaries linked to the interacting solids. Contact between fluid particles and fixed boundaries is accounted for by the incompressibility condition which naturally prevents the fluid nodes to penetrate into the solid boundaries [3,17,18,19,21].
156
157
==6.2 Contact between solid-solid interfaces==
158
159
The contact between two solid interfaces is treated by introducing a layer of '' contact elements'' between the two interacting solid interfaces. This layer is ''created during the mesh generation step'' by prescribing a minimum distance ([[Image:draft_Samper_231355424-picture-x0000_i1063.svg|17px]] ) between two solid boundaries. If the distance exceeds the minimum value ([[Image:draft_Samper_231355424-picture-x0000_i1064.svg|17px]] ) then the generated elements are treated as fluid elements. Otherwise the elements are treated as contact elements where a relationship between the tangential and normal forces and the displacements is introduced (Fig. 2).
160
161
This algorithm can also be used effectively to model frictional contact conditions between rigid or elastic solids in structural mechanics applications [2].
162
163
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
164
165
[[Image:draft_Samper_231355424-image42.png|center|384px]]
166
</div>
167
168
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
169
Figure 2: Modelling of contact conditions at a solid-solid interface with the PFEM</div>
170
171
<big>'''7. Modeling of bed erosion'''</big>
172
173
Prediction of bed erosion and sediment transport in open channel flows are important tasks in many areas of river and environmental engineering. Modeling of bed erosion is also relevant for predicting the evolution of surface material dragged in earth dams in overspill situations. Bed erosion is one of the main causes of environmental damage in floods.
174
175
The algorithm for modeling the erosion of soil/rock particles at the fluid bed is the following:
176
177
1.  Compute at every point of the bed surface the resultant tangential stress [[Image:draft_Samper_231355424-picture-x0000_i1065.svg|13px]] induced by the fluid motion. In 3D problems [[Image:draft_Samper_231355424-picture-x0000_i1066.svg|85px]] where [[Image:draft_Samper_231355424-picture-x0000_i1067.svg|16px]] and [[Image:draft_Samper_231355424-picture-x0000_i1068.svg|14px]] are the tangential stresses in the plane defined by the normal direction [[Image:draft_Samper_231355424-picture-x0000_i1069.svg|13px]] at the bed node. The value of [[Image:draft_Samper_231355424-picture-x0000_i1070.svg|13px]] for 2D problems can be estimated as shown in Fig.3.
178
179
2.  Compute the frictional work [[Image:draft_Samper_231355424-picture-x0000_i1071.svg|22px]] due to the tangential stresses at the bed surface (Fig. 3).
180
181
3.  The onset of erosion at a bed point occurs when the frictional work [[Image:draft_Samper_231355424-picture-x0000_i1072.svg|22px]] exceeds a critical threshold value [[Image:draft_Samper_231355424-picture-x0000_i1073.svg|21px]] . Then the node is detached from the bed region and it is allowed to move with the fluid flow. As a consequence, the mass of the patch of bed elements surrounding the bed node vanishes in the bed domain and it is transferred to the new fluid node. This mass is transported with the fluid.
182
183
4.  Sediment deposition can be modeled by an inverse process to that described above. Hence, a suspended node adjacent to the bed surface with a velocity below a threshold value is attached to the bed surface.
184
185
Details of the bed erosion algorithm described can be found in [17,18].
186
187
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
188
189
[[Image:draft_Samper_231355424-image52.png|center|414px]]
190
</div>
191
192
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
193
Figure 3: Modeling of bed erosion with the PFEM by dragging of bed material</div>
194
195
<big>'''8. Examples'''</big>
196
197
==8.1 Dragging of rocks by a water stream==
198
199
Predicting the speed at which a rock will be dragged by a water stream is of importance in many problems in hydraulic, harbour, civil and environmental engineering.
200
201
Figure 4 shows the study of the motion of a 1Tn quasi-spherical rock due to a water stream. Frictional conditions between the analyzed rock and the rest of the rocks (assumed to be rigid) are assumed. Figure 4a shows that a water stream of 1m/s can not displace the individual rock. An increase of the water speed to 2m/s induces the motion of the rock (Figure 4b).
202
203
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
204
<big>'''
205
{|
206
|-
207
| [[Image:draft_Samper_231355424-image53.png|center|252px]]
208
| [[Image:draft_Samper_231355424-image54.png|center|252px]]
209
|}
210
'''</big></div>
211
212
{|
213
|-
214
| style="width: 35.0132mm;vertical-align: top;" | 
215
| style="margin: 1em auto 1em auto;" | (a)                        
216
| style="margin: 1em auto 1em auto;" | (b)
217
|}
218
219
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
220
Figure 4: Study of the drag of an individual rock of 1Tn under a water stream at speeds of a) 1m/s and b) 2m/s</div>
221
222
==8.2 Impact of sea waves on piers and breakwaters==
223
224
Figure 5 shows the analysis of the effect of breaking waves on two different sites of a breakwater containing reinforced concrete blocks (each one of [[Image:draft_Samper_231355424-picture-x0000_i1074.svg|33px]] m). The figures correspond to the study of Langosteira harbour in A Coruña, Spain using PFEM.
225
226
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
227
<big>'''
228
{|
229
|-
230
| [[Image:draft_Samper_231355424-image56-c.png|center|234px]]
231
| [[Image:draft_Samper_231355424-image57-c.png|center|294px]]
232
|}
233
'''</big></div>
234
235
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
236
Figure 5: Breaking waves on breakwater slope containing reinforced concrete blocks</div>
237
238
==8.3 Soil erosion problems==
239
240
Figure 6 shows the capacity of the PFEM to modelling soil erosion, sediment transport and material deposition in a river bed. The soil particles are first detached from the bed surface under the action of the jet stream. Then they are transported by the flow and eventually fall down due to gravity forces into the bed surface at a downstream point.
241
242
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
243
244
{|
245
|-
246
| [[Image:draft_Samper_231355424-image58-c.png|center|222px]]
247
| [[Image:draft_Samper_231355424-image59-c.png|center|228px]]
248
|}
249
</div>
250
251
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
252
Figure 6: Erosion, transport and deposition of soil particles at a river bed due to a jet stream</div>
253
254
Figure 7 shows the progressive erosion of the unprotected part of a breakwater slope in the Langosteira harbour in A Coruña, Spain. The non protected upper shoulder zone is progressively eroded under the sea waves.
255
256
Other applications of the PFEM to bed erosion problems can be found in [17,18,21].
257
258
==8.4 Falling of a lorry into the sea by sea wave erosion of the road slope==
259
260
Figure 8 shows a representative example of the progressive erosion of a soil mass adjacent to the shore due to sea waves and the subsequent falling into the sea of a 2D object representing the section of a lorry. The object has been modeled as a rigid solid.
261
262
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
263
264
[[Image:draft_Samper_231355424-image60-c.png|center|336px]]
265
</div>
266
267
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
268
269
[[Image:draft_Samper_231355424-image61-c.png|center|336px]]
270
</div>
271
272
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
273
Figure 7: Erosion of an unprotected shoulder of a breakwater due to sea waves</div>
274
275
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
276
277
{|
278
|-
279
| [[Image:draft_Samper_231355424-image62-c.png|center|276px]]
280
| [[Image:draft_Samper_231355424-image63-c.png|center|276px]]
281
|}
282
</div>
283
284
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
285
286
[[Image:draft_Samper_231355424-image64-c.png|center|276px]]
287
</div>
288
289
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
290
Figure 8: Erosion of a soil mass due to sea waves and the subsequent falling into the sea of an adjacent lorry</div>
291
292
==8.5 Simulation of landslides==
293
294
Figure 9 shows a 3D simulation with the PFEM of a landslide falling on adjacent constructions. The landslide material has been modelled as a viscous incompressible fluid.
295
296
==8.6 The landslide in Lituya Bay==
297
298
We present some results of the 3D analysis of the landslide produced in Lituya Bay (Alaska) on July 9th 1958 (Fig. 10). The landslide was originated by an earthquake and movilized 90 millions tons of rocks that fell on the bay originating a large wave that reached a hight on the opposed slope of 524 m.
299
300
Figures 11 show images of the simulation of the landslide with PFEM. The sliding mass has been modelled as a continuum with a prescribed shear modulus. No frictional effect between the sliding mass and the underneath soil has been considered. The analysis has not taken into account the erosion and dragging of soil material induced by the landslide mass.
301
302
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
303
304
{|
305
|-
306
| [[Image:draft_Samper_231355424-image65-c.png|center|144px]]
307
| [[Image:draft_Samper_231355424-image66-c.png|center|132px]]
308
|}
309
310
{|
311
|-
312
| [[Image:draft_Samper_231355424-image67-c.png|center|144px]]
313
| [[Image:draft_Samper_231355424-image68-c.png|center|162px]]
314
|}
315
</div>
316
317
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
318
Figure 9: Landslide falling on constructions. 3D simulation using PFEM</div>
319
320
321
[[Image:draft_Samper_231355424-image69-c.jpeg|center|546px]]
322
323
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
324
Figure 10: Lituya Bay landslide. Landslide direction and maximum wave level [7]</div>
325
326
PFEM results have been compared with observed values of the maximum water level in the north hill adjacent to the reservoir. The maximum water level in this hill obtained with PFEM was 551 m. This is 5% higher than the observed value of 524 m. [7]. The maximum height location differs in 300 m from the observed value. In the south slope the maximum water height observed was 208 m, while the PFEM result (not shown here) was 195 m (6% error).
327
328
More information on the PFEM solutions of this example can be found in [20,22].
329
330
==8.7 Failure of a rockfill dam due to an overspill==
331
332
The PFEM has been used in combination with the classical (Eulerian) FEM for the analysis of the seapage and eventual failure of rockfill dams under overspill situations.
333
334
The evolution of the free surface of the fluid within the dam and over the rockfill slope is calculated using an Eulerian fixed mesh approach. A level set technique is used to track the evolution of the free surface. The Navier-Stokes equations for the fluid are modified to consider the presence of the porous media via a nonlinear (quadratic) Darcy law using Ergun coefficients that automatically go to zero when the flow exits the porous material. The balance equations are solved in term of Darcy velocity.
335
336
The structural response is evaluated using a continuum visco-rigid model and a Mohr-Coulomb failure criterion (with no cohesion) via the PFEM. Its specific features make it appropriate to treat the rockfill material undergoing large deformations and shape changes.
337
338
An alternative algorithm for projecting variables between non-matching meshes has been developed for passing information between the fixed mesh for the fluid and the moving mesh for the dams. Fig. 12 shows a scheme of the algorithm. For more details see [13,14].
339
340
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
341
'''
342
[[Image:draft_Samper_231355424-image70-c.jpeg|center|504px]]
343
'''</div>
344
345
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
346
Figure 11: Lituya Bay landslide. Evolution of the landslide into the reservoir obtained with the PFEM. Maximum level of generated wave (551 m) in the north slope</div>
347
348
Prototype dams like the one shown in Figure 13 have been analyzed and monitored during their failing process. Results are compared with experiments carried at the Technical University of Madrid (UPM) and the Center for Hydrographic studies of CEDEX. The bottom pressure distribution is registered when the failure achieves its stationary value for a given incoming discharge. The length of failure is also measured. Both these parameters are compared to numerical results (Figures 14 and 15).
349
350
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
351
352
[[Image:draft_Samper_231355424-image71-c.jpeg|center|414px]]
353
</div>
354
355
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
356
Figure 12: Failure of rockfill dam. Alternative solution algorithms.</div>
357
358
{| style="width: 100%;"
359
|-
360
| style="text-align: center;"|
361
[[Image:draft_Samper_231355424-image72.jpeg|center|252px]]
362
363
| style="text-align: center;"|
364
[[Image:draft_Samper_231355424-image73.jpeg|center|258px]]
365
366
|}
367
368
{| style="margin: 1em auto 1em auto;"
369
|-
370
| style="vertical-align: top;"|
371
[[Image:draft_Samper_231355424-image74.png|center|600px]]
372
373
[[Image:draft_Samper_231355424-image75-c.png|center|180px]]
374
375
|}
376
377
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
378
Figure 13. Experimental setting of rockfill dam failure carried out at UPM and CEDEX.</div>
379
380
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
381
(a) Failure length (B)</div>
382
383
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
384
385
[[Image:draft_Samper_231355424-image76.jpeg|center|402px]]
386
</div>
387
388
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
389
(b) Bottom pressure distribution</div>
390
391
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
392
393
[[Image:draft_Samper_231355424-image77-c.png|center|468px]]
394
</div>
395
396
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
397
Figure 14. Comparison between: (a) Experimental length of failure B<sub>exp</sub> = 1.50m and numerical one B<sub>num</sub> = 1.58m. (b) Experimental and numerical bottom pressure distributions.</div>
398
399
<big>'''Acknowledgements'''</big>
400
401
This research was supported by projects SEDUREC of the Consolider Programme of the Ministerio de Educación y Ciencia of Spain, SAFECON of the European Research Council of European Commission (EC) and the X-PRES and E-DAMS projects of the R+D National Plan of the Spanish Ministry for Science and Innovation [5,23]. Thanks are given to the Spanish construction company Dragados for financial support for the study of harbour engineering problems with PFEM.
402
403
<big>'''References'''</big>
404
405
:[1]  Aubry, R., Idelsohn, S.R. & Oñate, E. (2005). Particle finite element method in fluid mechanics including thermal convection-diffusion, Comput. & Struct. 83(17-18), 1459-1475.
406
407
:[2]  Carbonell, J.M., Oñate, E. & Suárez, B. (2010). Modeling of ground excavation with the Particle Finite Element Method. Journal of Engineering Mechanics (ASCE) 136(4), 455-463.
408
409
:[3]  Del Pin F., Idelsohn, S.R., Oñate, E. & Aubry, R. (2007). The ALE/Lagrangian particle finite element method: A new approach to computation of free-surface flows and fluid-object interactions. Computers & Fluids 36, 27-38
410
411
:[4]  Donea, J. & Huerta, A. (2003). Finite element method for flow problems. J. Wiley, 350 pp.
412
413
:[5]  E-DAMS. Numerical and experimental techniques for safety assessment and protection of embankment dams in overtopping scenarios. Project of National Plan for TRD of the Spanish Ministry for Science and Innovation BIA2010-21350-C03-00, 2010-2013.
414
415
:[6]  Edelsbrunner, H. & Mucke, E.P. (1999). Three dimensional alpha shapes. ACM Trans. Graphics 13, 43.
416
417
:[7]  Fritz, H.M., Hager, W.H. & Minor, H.E. (2001). Lituya Bay Case: Rockslide impact and wave run-up. Science of Tsunami Hazards 19(1), 3-22.
418
419
:[8]  Idelsohn, S.R., Calvo, N. & Oñate, E.  (2003c). Polyhedrization of an arbitrary point set. Comput. Method Appl. Mech. Engng. 192(22-24), 2649-2668.
420
421
:[9]  Idelsohn, S.R., Oñate, E. & Del Pin, F. (2004). The particle finite element method: a powerful tool to solve incompressible flows with free-surfaces and breaking waves. Int. J. Num. Meth. Engng. 61, 964-989.
422
423
:[10]  Idelsohn, S.R., Oñate, E., Del Pin, F. & Calvo, N. (2006). Fluid-structure interaction using the particle finite element method. Comput. Meth. Appl. Mech. Engng. 195, 2100-2113.
424
425
:[11]  Idelsohn, S.R., Marti, J., Limache, A., Oñate, E. (2008). Unified Lagrangian formulation for elastic solids and incompressible fluids: Application to fluid-structure interaction problems via the PFEM. Comput Methods Appl Mech Engrg. 197, 1762-1776.
426
427
:[12]  Idelsohn, S.R., Mier-Torrecilla, M. & Oñate, E. (2009). Multi-fluid flows with the Particle Finite Element Method. Comput Methods Appl Mech Engrg. 198, 2750-2767.
428
429
:[13]  Larese, A., Rossi, R., Oñate, E. & Idelsohn, S.R. (2008). Validation of the Particle Finite Element Method (PFEM) for simulation of free surface flows. Engng. Comput. 25(4), 385-425.
430
431
:[14]  Larese, A., Rossi, R., Oñate, E. & Toledo, M.A. (2010). Physical and Numerical Modelization of the Behavior of rockfill dams during overtopping scenarios. Dam Maintenance and Rehabilitation II. CRCpress Balkema.
432
433
:[15]  Oñate, E. (2004). Possibilities of finite calculus in computational mechanics. Int. J. Num. Meth. Engng. 60(1), 255-281.
434
435
:[16]  Oñate, E., Idelsohn, S.R., Del Pin, F. & Aubry, R. (2004b). The particle finite element method. An overview. Int. J. Comput. Methods 1(2), 267-307.
436
437
:[17]  Oñate, E., Celigueta, M.A. & Idelsohn, S.R. (2006c). Modeling bed erosion in free surface flows by the Particle Finite Element Method, Acta Geotechnia 1(4), 237-252.
438
439
:[18]  Oñate, E., Idelsohn, S.R., Celigueta, M.A. & Rossi, R. (2008). Advances in the particle finite element method for the analysis of fluid-multibody interaction and bed erosion in free surface flows. Comput. Meth. Appl. Mech. Engng. 197(19-20), 1777-1800.
440
441
:[19]  Oñate, E., Rossi, R., Idelsohn, S.R. & Butler, K. (2010). Melting and spread of polymers in fire with the Particle Finite Element Method. Int. J. Numerical Methods Engng., 81(8), 1046-1072.
442
443
:[20]  Oñate, E., Salazar, F. & Morán, R. (2011). Modeling of landslides into reservoirs with the Particle Finite Element Method. Research Report CIMNE No. PI355. Submitted to Int. J. Numerical Methods in Geomechanics.
444
445
:[21]  Oñate, E., Celigueta, M.A., Idelsohn, S.R., Salazar, F. & Suárez, B. (2011). Possibilities of the Particle Finite Element Method for fuid-soil-structure interaction problems. Computational Mechanics. Accepted for publication.
446
447
:[22]  Salazar, F., Oñate, E. & Morán, R. (2011). Modelación numérica de deslizamientos de ladera en embalses mediante el método de partículas y elementos finitos (PFEM). Rev. Int. Mét. Num. Cálc. Dis. Ing. Accepted for publication
448
449
:[23]  XPRES. Development of a method for study of the failure process of rockfill dams in overtopping situations combining particle based method and finite element method. National Plan for RTD of Spanish Ministry for Science and Innovation I+D BIA2007-68120-C03-01, 2007-2010.
450
451
:[24]  Zienkiewicz, O.C. & Taylor, R.L. (2005). The finite element method for solid and structural mechanics, Elsevier.
452
453
:[25]  Zienkiewicz, O.C., Taylor, R.L. & Nithiarasu, P. (2006). The finite element method for fluid dynamics, Elsevier.
454

Return to Onate et al 2012c.

Back to Top