Application of Hagstrom-Warburton high-order truncation boundary condition on time harmonic analysis of concrete arch dam-reservoir systems

2021 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Vahid Lotfi ◽  
Ali Lotfi

Purpose Subsequently, the response of idealized Morrow Point arch dam is studied due to stream, vertical and cross-stream ground motions for reservoir bottom/sidewalls conditions of both fully reflective and absorptive. For each combination, different orders of Hagstrom–Warburton (HW) condition are evaluated from accuracy point of view by comparing them against exact solutions. It should be emphasized that normalized length of reservoir near-field region is taken as a very low value of L/H = 0.2 during this process which makes it a very challenging test for any kind of truncation boundary condition. Design/methodology/approach In present study, dynamic analysis of concrete arch dam-reservoir systems is formulated by FE-(FE-TE) approach [i.e. finite element-(finite element-truncation element)]. In this technique, dam and reservoir are discretized by solid and fluid finite elements. Moreover, the HW high-order condition imposed at the reservoir truncation boundary. This task is formulated by employing a truncation element at that boundary. It is emphasized that reservoir far-field is excluded from the discretized model. The formulation is initially explained in details. Findings The trend in gaining accuracy with increase in order of HW condition were basically the same for all three types of excitations under both full reflective and absorptive reservoir bottom/sidewalls conditions. The only exception was for cross-stream excitation response which was showing less accurate results near the first major peak for moderate orders of HW (e.g. O3-2) in comparison to what was observed for responses due to symmetric excitations (stream and vertical). This is mainly attributed to the selection of evanescent-type parameters of HW condition which is based on the first symmetric mode of reservoir. However, it is noted that error diminishes even for cross-stream excitation as order increases. High orders of HW condition, such as O5-5 considered herein, generate highly accurate responses for all three possible excitations under both types of full reflective and absorptive reservoir bottom/sidewalls conditions. It is such that responses are hardly distinguishable from corresponding exact responses. This reveals that proposed FE-(FE-TE) analysis technique based on HW condition is quite successful, and one may fully rely on that for accurate and efficient analysis of concrete arch dam-reservoir systems. Originality/value Dynamic analysis of concrete arch dam-reservoir system is formulated by new method. HW high-order condition is applied for a very low and challenging reservoir length. Different orders are evaluated against exact solution with excellent agreement. Generalized matrices of truncation element are derived for FEM programmers. The method is discussed for all types of excitation and reservoir base conditions.

2021 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Ouadie Koubaiti ◽  
Said EL Fakkoussi ◽  
Jaouad El-Mekkaoui ◽  
Hassan Moustachir ◽  
Ahmed Elkhalfi ◽  
...  

Purpose This paper aims to propose a new boundary condition and a web-spline basis of finite element space approximation to remedy the problems of constraints due to homogeneous and non-homogeneous; Dirichlet boundary conditions. This paper considered the two-dimensional linear elasticity equation of Navier–Lamé with the condition CAB. The latter allows to have a total insertion of the essential boundary condition in the linear system obtained; without using a numerical method as Lagrange multiplier. This study have developed mixed finite element; method using the B-splines Web-spline space. These provide an exact implementation of the homogeneous; Dirichlet boundary conditions, which removes the constraints caused by the standard; conditions. This paper showed the existence and the uniqueness of the weak solution, as well as the convergence of the numerical solution for the quadratic case are proved. The weighted extended B-spline; approach have become a much more workmanlike solution. Design/methodology/approach In this paper, this study used the implementation of weighted finite element methods to solve the Navier–Lamé system with a new boundary condition CA, B (Koubaiti et al., 2020), that generalises the well-known basis, especially the Dirichlet and the Neumann conditions. The novel proposed boundary condition permits to use a single Matlab code, which summarises all kind of boundary conditions encountered in the system. By using this model is possible to save time and programming recourses while reap several programs in a single directory. Findings The results have shown that the Web-spline-based quadratic-linear finite elements satisfy the inf–sup condition, which is necessary for existence and uniqueness of the solution. It was demonstrated by the existence of the discrete solution. A full convergence was established using the numerical solution for the quadratic case. Due to limited regularity of the Navier–Lamé problem, it will not change by increasing the degree of the Web-spline. The computed relative errors and their rates indicate that they are of order 1/H. Thus, it was provided their theoretical validity for the numerical solution stability. The advantage of this problem that uses the CA, B boundary condition is associated to reduce Matlab programming complexity. Originality/value The mixed finite element method is a robust technique to solve difficult challenges from engineering and physical sciences using the partial differential equations. Some of the important applications include structural mechanics, fluid flow, thermodynamics and electromagnetic fields (Zienkiewicz and Taylor, 2000) that are mainly based on the approximation of Lagrange. However, this type of approximation has experienced a great restriction in the level of domain modelling, especially in the case of complicated boundaries such as that in the form of curvilinear graphs. Recently, the research community tried to develop a new way of approximation based on the so-called B-spline that seems to have superior results in solving the engineering problems.


Sign in / Sign up

Export Citation Format

Share Document