For the more homogenous specimen, the models with local size dependent parameters and models with global uniform parameters give similar response. The models with uniform parameters evaluated according to the averages of the geometric parameters used in the evaluation of local parameters ensure better agreement with the respective models with size-dependent parameters than the models with uniform parameters evaluated according to the particle radii. Simulations using the more heterogenous specimen reveal differences between the considered models. There are significant differences in stress–strain curves as well as in the failure pattern. The models with local size-dependent parameters are more sensitive to the change of heterogeneity than the model with global uniform parameters.

The idea is to use the information existing at time t=tn in the velocity streamlines as well as in the acceleration streamlines to update the particle position as well as the velocity in an updated Lagrangian frame. The method may be used with moving or fixed meshes.

In this report some validation cases are described, where the impact of sliding blocks against a mass of water is simulated, as well as the generation and propagation of the subsequent impulse wave. They are based on the experiments performed by Sælevik et al. [21].

- To develop a simulation program of the behaviour of upwind sails and
rigging, to help the crew to optimize the performance of the sailing yacht in
real time. For this purpose, it will be necessary to formulate a fluid-structure
interaction algorithm to compute the performance of a particular sail/rigging
configuration. Since the crew dynamically trims the rigging and sails, in order
to evaluate the performance of the actual configuration, a tool to monitor
the rigging and sails will be necessary, too.
- To adjust a monitoring element to quantify in physic values the manoeuvre of
the crew. This will be our monitoring tool.

- To reproduce the crew manoeuvre in the simulation program with the data
obtained with the monitoring tool. Once the sail/rigging configuration has
been adapted ‘in real time’ to the actual one, the performance of this new
configuration can be computed. For this purpose the simulation program and
the monitoring tool must communicate among them.

Por otro lado, se ha simulado el fenómeno estimando principalmente la iniciación y la
propagación del daño interlaminar e igualmente los efectos inducidos en la respuesta
dinámica. En el enfoque numérico se trata el material compuesto como un sistema
microestructural, en el cual los fallos surgen en la interacción entre los materiales
constituyentes. Para reproducir la degradación se emplea una estrategia de reducción
localizada de la rigidez elástica del material que no demanda una intervención durante
el preproceso.
Los resultados empíricos aportan conclusiones relevantes en relación al grado de
sensibilidad y a la adecuación de los diferentes criterios de correlación modal para la
identificación del daño inducido por un impacto. Los nuevos criterios de correlación
definidos y el análisis conjunto de las propiedades estáticas y dinámicas residuales, han
permitido acotar el intervalo de incertidumbre y reducir la limitación actual en relación
a la deformación máxima admisible a compresión de los laminados. La herramienta
de cálculo desarrollada permite simular el comportamiento vibratorio de laminados
compuestos definiendo el material y su estado de degradación en la microescala. Los
resultados corroboran la viabilidad del enfoque microestructural para la simulación
del fenómeno.

Todavía el trabajo aporta la conclusión más importante en cuanto a la distribución de
las tensiones tangenciales y con ésta el cálculo del ffacttor de cortte, que desarrollamos
en la PARTE SEGUNDA denominándolo como TEOREMA GENERAL. Afirma que una
vez hallada la distribución tangencial en algún tipo de barra gruesa –recta, curva o
ficticia– un cálculo directo, con apropiados intercambios de parámetros físico–geométricos
que dependen del tipo de barra y de la distribución de las no homogeneidades y forma de la
sección transversal de la barra, permite conocer la distribución de las tensiones
tangenciales en los otros dos tipos de barra con la misma sección.
Se incluyen otros dos resultados originales que se denominan métodos I y II de
superposición –para los tres tipos de barras gruesas que el trabajo aborda– que, cuando
las secciones de las barras modifican “a saltos” sus propiedades elásticas y de densidad,
puede hallarse el régimen tensional combinando linealmente regímenes conocidos de
secciones homogéneas. Esto simplifica notablemente el trabajo de búsqueda para estos
tipos de sección no homogénea bajo estudio. Se encuentra también un acople flexo–axial
del movimiento y se trabaja y se generaliza el concepto del corrimiento del eje neutro,
imposición clásica en barras homogéneas de gran curvatura que permite no sólo
simplificar el proceso algebraico sino que extiende naturalmente la definición de eje
neutro que se utiliza tradicionalmente en barras rectas homogéneas.
Al desarrollar el ítem de vibraciones naturales de las barras gruesas partiendo de las
ecuaciones de movimiento, tanto rectas como curvas, se encuentran las condiciones de
ortogonalidad entre formas modales (extensión ad-hoc del tradicional modelo de Sturm–
Liouville [6]) para ser utilizadas en una eventual superposición modal clásica en problemas
lineales y separables, para hallar la respuesta dinámica del sistema (Vibración Forzada).
Algunos Apéndices y varios Ejemplos resueltos analítica y numéricamente completan el
trabajo que permite inferir que los resultados encontrados coinciden con muy buena
precisión con los hallados con otras metodologías aproximadas como son las de elementos
finitos en 2D y 3D pero por medio de un encuadre mucho más directo, sencillo y abarcativo.
Cabe todavía agregar que la propuesta para llegar al régimen de tensiones
tangenciales de las secciones de formas arbitrarias constitutivamente no homogéneas,
y aún múltiplemente conexas, reemplaza a las utilizadas comúnmente para barras
rectas homogéneas y que son conocidas como de Collignon o de Jourawsky. Por otro
lado estas metodologías tradicionales serían prácticamente inadecuadas de utilizar
para ciertos casos dentro del espectro de aplicaciones que presentamos.
Fundamentalmente, lo dicho permite hallar el factor de corte de las secciones de
barras gruesas y entonces completar los coeficientes del sistema diferencial que
gobierna el movimiento de estos tipos estructurales.
Por último entendemos que el Resullttado Fundamenttall, el TEOREMA GENERAL, los
Métodos de Superposición, las Ecuaciones de Movimiento y las Condiciones de
Ortogonalidad presentadas, son resultados originales dentro de la bibliografía afín.

It was found that the use of some standard practices (e. g.M-Matrices theory) for the design of essentially non-oscillatory numerical methods is not appropriate when consistency recovery methods are employed. Hence, with respect to convective stabilization, such recovery methods are not preferred. Next, we present the design of a high-resolution Petrov–Galerkin (HRPG) method for the CDR problem. The structure
of the method in 1d is identical to the consistent approximate upwind (CAU) Petrov–Galerkin method [68] except for the definitions of the stabilization parameters.

Such a structure may also be attained via the Finite Calculus (FIC) procedure [141] by an appropriate definition of the characteristic length. The prefix high-resolution is used here in the sense popularized by Harten, i. e.second order accuracy for smooth/regular regimes and good shock-capturing in non-regular regimes. The design procedure in 1d embarks on the problem of circumventing the Gibbs phenomenon observed in L2 projections. Next, we study the conditions on the stabilization parameters to circumvent the global oscillations due to the convective term. A conjuncture of the two results is made to deal with the problem at hand that is usually plagued by Gibbs, global and dispersive oscillations in the numerical solution. A multi dimensional extension of the HRPG method using multi-linear block finite elements is also presented.
Next, we propose a higher-order compact scheme (involving two parameters) on structured meshes for the Helmholtz equation. Making the parameters equal, we recover the alpha-interpolation of the Galerkin finite element method (FEM) and the classical central finite difference method. In 1d this scheme is identical to the alphainterpolation method [140] and in 2d choosing the value 0.5 for both the parameters, we recover the generalized fourth-order compact Padé approximation [81, 168] (therein using the parameter
gamma = 2). We follow [10] for the analysis of this scheme and its performance on square meshes is compared with that of the quasi-stabilized FEM [10]. Generic expressions for the parameters are given that guarantees a dispersion accuracy of sixth-order should the parameters be distinct and fourth-order should they be equal. In the later case, an expression for the parameter is given that minimizes the maximum relative phase error in 2d. A Petrov–Galerkin formulation that yields the aforesaid scheme on structured meshes is also presented. Convergence studies of the
error in the L2 norm, the H1 semi-norm and the l1 Euclidean norm is done and the pollution effect is found to be small.

Finally, we present a collection of stabilized FE methods derived via first-order and second-order FIC procedures for the Stokes problem. It is shown that several well known existing stabilized FE methods such as the penalty technique, the Galerkin Least Square (GLS) method [93], the Pressure Gradient Projection (PGP) method [35] and the orthogonal sub-scales (OSS) method [34] are recovered from the general
iii residual-based FIC stabilized form. A new family of Pressure Laplacian Stabilization (PLS) FE methods with consistent nonlinear forms of the stabilization parameters are derived. The distinct feature of the family of PLS methods is that they are nonlinear and residual-based, i. e.the stabilization terms depend on the discrete residuals of the momentum and/or the incompressibility equations. The advantages and disadvantages of these stabilization techniques are discussed and several examples of application are presented.

(2) Estado MCI ocasionado por
AD; y (3) Etapa de demencia ocasionada por AD. Por otro lado, los modelos SVM
estarán enfocados a dos tareas principales: (1) Diagnóstico de AD mediante la discriminación entre sujetos sanos y sujetos con AD; y (2) Predicción de AD, orientada
a la discriminación de sujetos MCI con riesgo de convertirse a AD y sujetos MCI
sin riesgo de conversión. Asimismo, estos modelos deberán garantizar resultados
aceptables, respecto a la sensibilidad y especificidad de las tareas de clasificación.

Los resultados obtenidos en esta investigación son prometedores. Por un lado, el
subconjunto de variables seleccionadas como relevantes para el diagnóstico de AD,
tienen correlación con los resultados de investigaciones previas. Asimismo, en la
etapa de testeo, los resultados demostraron que los modelos SVM son de gran utilidad
para el soporte diagnóstico clínico de esta enfermedad, siendo capaces de discriminar
sujetos con AD de sujetos sanos (diagnóstico) con una exactitud mayor al 99% y
distinguir a los sujetos MCI con riesgo de conversión a AD de los sujetos MCI sin
riego de conversión (predicción) con una exactitud superior al 94%.

The principal aim of the present work is the development of a computational method
to simulate the overtopping and the beginning of failure of the downstream shoulder of
a rockfill dam. The whole phenomenon is treated in a continuous framework.

The fluid free surface problem outside and inside the rockfill slop e is treated using a
unique Eulerian fixed mesh formulation. A level set technique is employed to track the
evolution of the free surface. The traditional Navier-Stokes equations are modified in
order to automatically detect the presence of the porous media. The non-linear seepage
is evaluated using a quadratic form of the resistance law for which the Ergun's coefficients
have been chosen.

This work presents new developments of the discrete element method improving efficiency and accuracy of modelling of rock-like materials, especially in excavation processes. All the numerical algorithms has been implemented in an in-house software, which was then used to run numerical examples.

The basic formulation of DEM with linear elastic-perfectly brittle contact model is presented. The main difference with other models found in the literature is the consideration of global stiffness and strength parameters that are constants in the whole model.
The result of a simulations is strongly related with the configuration of the particle assembly used. Particle assemblies should be suficiently compact and ensure the isotropy to reproduce the physical properties of the modelled material. This work presents a novel technique for the generation of highly dense particle assemblies in arbitrary geometries, satisfying all the requirements for accurate discrete element simulations.

One of the key issues in the use of the DEM is the estimation of the contact model parameters. A methodology is proposed for the estimation of the contact model parameters yielding required macroscopic properties of the material. The relationships between the contact model parameters and the mechanical properties of brittle materials, as well as the influence of the particles assembly configuration on the macroscopic properties, are analysed.

A major dificulty in the application of the DEM to real engineering problems is the high computational cost in simulation involving a large number of particles. The most common way to solve this is the use of parallel computing techniques, where multiple processors are used. As an alternative, a coupling scheme between DEM and the finite element method (FEM) is proposed in the work. Within the hybrid DEM/FEM model, DEM is only used in the region of the domain where it provides an advantage over a continuum-based approach, as the FEM. The coupling is dynamically adapted, starting with the whole domain discretized with FEM. During the simulation, in the regions where a high stress level are found, a change of modelling method from continuum FEM to the discrete DEM is employed.

Finally, all the developments are applied to the simulation of a real excavation process. An analysis of the rock cutting process with TBM disc cutters is performed, where DEM and the DEM/FEM coupling technique presents an important advantage over other simulation techniques.