Probabilistic nonlinear analysis of CFR dams by MCS using Response Surface Method

This paper presents the probabilistic analysis of concrete-faced rockfill (CFR) dams according to the Monte Carlo Simulation (MCS) results which are obtained through the Response Surface Method (RSM). ANSYS finite element program is used to get displacement and principal stress components. First of all, some parametric studies are performed according to the simple and representative finite element model of dam body to obtain the optimum approximate model. Secondly, a sensitivity analysis is performed to get the most effective parameters on dam response. Then, RSM is used to obtain the approximate function through the selected parameters. After the performed analyses, star experimental design with quadratic function without mixed terms according to the k = 1 is determined as the most appropriate model. Finally, dam-foundation-reservoir interaction finite element model is constituted and probabilistic analyses are performed with MCS using the selected parameters, sampling method, function and arbitrary factor under gravity load for empty and full reservoir conditions. Geometrically and materially nonlinearity are considered in the analysis of dam-foundation-reservoir interaction system. Reservoir water is modeled by fluid finite elements based on the Lagrangian approach. Structural connections are modeled as welded contact and friction contact based on Coulomb’s friction law. Probabilistic displacements and stresses are presented and compared with deterministic results.     

Shape optimal design of arch dams including dam-water–foundation rock interaction using a grading strategy and approximation concepts

Optimal design of arch dams including dam-water–foundation rock interaction is achieved using the soft computing techniques. For this, linear dynamic behavior of arch dam-water–foundation rock system subjected to earthquake ground motion is simulated using the finite element method at first and then, to reduce the computational cost of optimization process, a wavelet back propagation neural network (WBPNN) is designed to predict the arch dam response instead of directly evaluating it by a time-consuming finite-element analysis (FEA). In order to enhance the performance generality of the neural network, a dam grading technique (DGT) is also introduced. To assess the computational efficiency of the proposed methodology for arch dam optimization, an actual arch dam is considered. The optimization is implemented via the simultaneous perturbation stochastic approximation (SPSA) algorithm for the various conditions of the interaction problem. Numerical results show the merits of the suggested techniques for arch dam optimization. It is also found that considering the dam-water–foundation rock interaction has an important role for safely designing an arch dam.     

Advances in concrete arch dams shape optimization

This paper presents an efficient methodology to find the optimum shape of arch dams. In order to create the geometry of arch dams a new algorithm based on Hermit Splines is proposed. A finite element based shape sensitivity analysis for design-dependent loadings involving body force, hydrostatic pressure and earthquake loadings is implemented. The sensitivity analysis is performed using the concept of mesh design velocity. In order to consider the practical requirements in the optimization model such as construction stages, many geometrical and behavioral constrains are included in the model in comparison with previous researches. The optimization problem is solved via the sequential quadratic programming (SQP) method. The proposed methods are applied successfully to an Iranian arch dam, and good results are achieved. By using such methodology, efficient software for shape optimization of concrete arch dams for practical and reliable design now is available.     

Shape optimal design of arch dams using an adaptive neuro-fuzzy inference system and improved particle swarm optimization

An efficient methodology is proposed to find the optimal shape of arch dams including fluid–structure interaction subject to earthquake ground motion. In order to reduce the computational cost of optimization process, an adaptive neuro-fuzzy inference system (ANFIS) is built to predict the dam effective response instead of directly evaluating it by a time-consuming finite element analysis (FEA). The presented ANFIS is compared with a widespread neural network termed back propagation neural network (BPNN) and it appears a better performance generality for estimating the dam response. The optimization task is implemented using an improved version of particle swarm optimization (PSO) named here as IPSO. In order to assess the effectiveness of the proposed methodology, the optimization of a real world arch dam is performed via both IPSO–ANFIS and PSO–BPNN approaches. The numerical results demonstrate the computational advantages of the proposed IPSO–ANFIS for optimal design of arch dams when compared with the PSO–BPNN approach.     

Non-linear seismic response of concrete gravity dams to near-fault ground motions including dam-water-sediment-foundation interaction

This study focuses on non-linear seismic response of a concrete gravity dam subjected to near-fault and far-fault ground motions including dam-water-sediment-foundation rock interaction. The elasto-plastic behavior of the dam concrete is idealized using Drucker–Prager yield criterion based on associated flow rule assumption. Water in the reservoir is represented by 9-noded isoparametric quadrilateral fluid finite elements while the dam, the foundation rock and the sediment layer are modeled by using 8-noded isoparametric quadrilateral solid finite elements. The program NONSAP modified for elasto-plastic analysis of fluid-structure systems using the Lagrangian fluid finite element is employed in the response calculations. The fluid element includes the effects of surface waves and sloshing behavior of fluids. Non-linear seismic analyses of the selected concrete dam subjected to both near-fault and far-fault ground motions are performed. The results obtained from linear and non-linear analyses are compared with each other.     

Comparison of coupled and decoupled modal approaches in seismic analysis of concrete gravity dams in time domain

Different methods are available for dynamic analysis of concrete dams. Among these, modal approach is highly popular due to the efficiency of the method. This becomes more significant if the response is to be calculated for several earthquake ground motion records as required in most practical cases. In this study, two different modal approaches have been considered for dynamic analysis of concrete gravity dam–reservoir systems. These are coupled modal method and decoupled modal technique. The former approach utilizes the coupled modes of the system. It is well known that calculation of these modes involves some complications due to its corresponding unsymmetrical eigenproblem. However, the response at each step can be obtained very efficiently in this method. The latter technique, relies on the decoupled modes of the system, which are easily obtained by standard eigenvalue solvers. The equation of motion is also solved with reasonable efficiency in this approach. In the latter part of this paper, analysis of a typical dam–reservoir system is performed by both methods mentioned above. These analyses are compared from accuracy and efficiency point of view.     

Damage analysis of arch concrete dams subjected to underwater explosion

The nonlinear behavior of arch concrete dams subjected to underwater explosions was investigated. Pressure waves were spherically modeled using a 3-D finite element method. The nonlinear behavior of concrete was modeled using the concrete damage plasticity model. In addition, absolute wave formulation was used to develop a fully-coupled model incorporating the propagation of the shock wave and water–dam interaction. Analysis of an arch concrete dam subject to underwater explosion was performed upon model verification. The dynamic response of the dam subject to the underwater explosion was studied for different sizes of explosions at different depths. The results showed that the closer the point to the explosive source, the sooner the impression and damage was caused by the blast. It was clear that the middle part of the dam facing the explosion was the first location to incur damage and then cracks extended to the downstream face and abutments. Finally, cracks extended to the bottom of the dam. The results of the dynamic analysis and dissipated energy by the loss of elements’ resistance indicated that the time slice of 150 ms was appropriate to analyze arch concrete dam structure subject to an underwater explosion. In addition, for a given amount of explosives, an increase in the depth of explosion corresponded to an increase in the maximum displacement of the upstream face of the dam. The results also confirmed that the damage distribution to the arch concrete dam subject to an underwater explosion depended upon the amount and depth of explosive materials.     

Axisymmetric contact of a punch of polynomial profile with an elastic half-space when there is friction and adhesion

The axisymmetric problem of the contact interaction of a punch of polynomial profile and an elastic half-space when there is friction and partial adhesion in the contact area is considered. Using the Wiener–Hopf method the problem is reduced to an infinite system of algebraic Poincare–Koch equations, the solution of which is obtained in series. The radii of the contact area and of the adhesion zone, the distribution of the contact pressures and the indentation of the punch are obtained.     

Modelación del movimiento del agua en medios porosos no lineales mediante un esquema de diferencias finitas. Aplicación al sobrevertido en presas de escollera

Flow through coarse porous media, as rockfill, gravel or sand of big size, exhibit a non linear relation between hydraulic gradient and seepage velocity. This behavior is observed when water flows over the crest of a rockfill dam, seeps into the downstream shoulder and moves through the rockfill voids. A finite difference scheme is used here to solve this problem through numerical modeling. It is implemented in the code MNLEE.RDS for solving dam overtopping problems. In addition, a procedure is provided to solve the coupled problem of flow over and through the rockfill dam. For that purpose, the concept of «law of flow exchange»is introduced.     

一个第二类变分不等式的有限元逼近

本短文讨论下述第二类变分不等式（见 ［２, ４］）的有限元逼近及其误差分析：其中是平面凸多边形区域的的边界, 且而 ． 诸如热量控制问题,流体通过半可透性壁的扩散问题以及简化库仑摩擦接触问题的正则化方法等均可归为上述变分不等式（１）（见［２,３］）．在文［２］中给出了上述变分不等式的有限元逼近格式,作出了收敛性分析及误差估计．本文的目的是进一步用数值积分简化上述有限元逼近格式并改进原有的估计误差． 设Ｔｈ是的拟一致三角形部分,Ｖｈ是对应的线性元空间,且使得ｖｈ＝０在上．［２］中用数值积分代替其中 Ｍｉ…     

Conjugate gradient techniques for the optimal control evolution dam problem

This paper considers the numerical simulation of optimal control evolution dam problem by using conjugate gradient method.The paper considers the free boundary value problem related to time dependent fluid flow in a homogeneous earth rectangular dam.The dam is taken to be sufficiently long that the flow is considered to be two dimensional.On the left and right walls of the dam there is a reservoir of fluid at a level dependent on time.This problem can be transformed into a variational inequality on a fixed domain.The numerical techniques we use are based on a linear finite element method to approximate the state equations and a conjugate gradient algorithm to solve the discrete optimal control problem.This algorithm is based on Armijo's rule in the unconstrained optimization theory.The convergence of the discrete optimal solutions to the continuous optimal solutions,and the convergence of the conjugate gradient algorithm are proved.A numerical example is given to determine the location of the minimum surface     

溃坝问题的间断有限元方法

本文研究90年代初提出的Runge-Kutta间断Galerkin有限元方法，给出该方法的精度分析，通过经典算例验证该方法处理间断问题、捕捉锐利波形的能力，并将其推广到求解浅水问题．针对坝底无摩擦，无坡度的理想情形进行讨论，给出方溃坝和圆溃坝问题的数值模拟结果．     

Modeling of large deformation frictional contact in powder compaction processes

In this paper, the computational aspects of large deformation frictional contact are presented in powder forming processes. The influence of powder–tool friction on the mechanical properties of the final product is investigated in pressing metal powders. A general formulation of continuum model is developed for frictional contact and the computational algorithm is presented for analyzing the phenomena. It is particularly concerned with the numerical modeling of frictional contact between a rigid tool and a deformable material. The finite element approach adopted is characterized by the use of penalty approach in which a plasticity theory of friction is incorporated to simulate sliding resistance at the powder–tool interface. The constitutive relations for friction are derived from a Coulomb friction law. The frictional contact formulation is performed within the framework of large FE deformation in order to predict the non-uniform relative density distribution during large deformation of powder die pressing. A double-surface cap plasticity model is employed together with the nonlinear contact friction behavior in numerical simulation of powder material. Finally, the numerical schemes are examined for efficiency and accuracy in modeling of several powder compaction processes.     

Efficient Reduction Method based on Trial Vector Derivatives for Nonlinear Transient Analysis of Multilayer Sheet Structures

For the dynamic response of multilayer sheet structures it is important to take nonlinear contact and friction forces between the layers into account. The nonlinearities inside the joints can significantly influence the dynamic behavior of the entire structure. For the dynamic analysis of such structures conventional computational approaches have an unbalanced ratio between computational time and accuracy. In order to overcome this problem, so-called “Joint Interface Modes” (JIMs) are used. These JIMs, which are problem-oriented trial vectors, extend the reduction base of classical trial vector based reduction methods. In this contribution a new calculation approach for JIMs based on trial vector derivatives is introduced. The result quality of the presented method is comparable to the direct finite element method but the computational effort is extremely low due to model order reduction. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The element‐free Galerkin method for a quasistatic contact problem with the Tresca friction in elastic materials

This paper is proposed for the error estimates of the element‐free Galerkin method for a quasistatic contact problem with the Tresca friction. The penalty method is used to impose the clamped boundary conditions. The duality algorithm is also given to deal with the non‐differentiable term in the quasistatic contact problem with the Tresca friction. The error estimates indicate that the convergence order is dependent on the nodal spacing, the time step, the largest degree of basis functions in the moving least‐squares approximation, and the penalty factor. Numerical examples demonstrate the effectiveness of the element‐free Galerkin method and verify the theoretical analysis.     

求解摩擦接触问题的一个非内点光滑化算法

给出了一个求解三维弹性有摩擦接触问题的新算法,即基于NCP函数的非内点光滑化算法．首先通过参变量变分原理和参数二次规划法,将三维弹性有摩擦接触问题的分析归结为线性互补问题的求解;然后利用NCP函数,将互补问题的求解转换为非光滑方程组的求解;再用凝聚函数对其进行光滑化,最后用NEWTON法解所得到的光滑非线性方程组．方法具有易于理解及实现方便等特点．通过线性互补问题的数值算例及接触问题实例证实了该算法的可靠性与有效性．     

Numerical Analysis of a Dynamic Contact Problem with History-Dependent Operators

In this paper, we study a dynamic contact model with long memory which allows both the convex potential and nonconvex superpotentials to depend on history-dependent operators. The deformable body consists of a viscoelastic material with long memory and the process is assumed to be dynamic. The contact involves a nonmonotone Clarke subdifferential boundary condition and the friction is modeled by a version of the Coulomb's law of dry friction with the friction bound depending on the total slip. We introduce and study a fully discrete scheme of the problem, and derive error estimates for numerical solutions. Under appropriate solution regularity assumptions, an optimal order error estimate is derived for the linear finite element method. This theoretical result is illustrated numerically.     

弹性接触问题及线性有限元逼近

1．引言及预备知识最早研究简化了的弹性接触的数学问题（称为单边问题）应属A．Signorini（1933年）．到70年代初，G．Fichera（1972）[4]及G．Duvaut和J．L．Lions（1972）问对此类问题有进一步的研究，建立了问题的解的存在性，正则性以及唯一性．R．Glowinski等的专著[5]中讨论了问题的数值逼近（也可见问中所列文献）．特别应提到捷克数学家J．Ne6as，1．HlavA6ek等所进行的工作（见［7］及其中的文献）以及Kikuchi和Oden［9］，从数学问题的形成到与之等价的变分问题解的存在性，以及有限元逼近的收敛性和误差估计都有研究．本…     

Contact characteristics of viscoelastic bonded layers

A viscoelastic layered contact model has successfully been developed and solved analytically. The single layered linear viscoelastic material is assumed to be perfectly bonded to a rigid substrate in contact with a rigid indenter without friction under a step load. Two cases are considered: (a) a compressible layered material with a typical Poisson's ratio of 0.4 and (b) an incompressible layer with a Poisson's ratio of 0.5. Two viscoelastic models: Maxwell and three element standard linear solid are investigated. This paper highlights the methodology employed and the results obtained under various conditions.