Sharif Digital Repository

In this paper we investigate a three-dimensional finite-strain phenomenological constitutive model and propose an efficient solution algorithm by properly defining the variables and by avoiding extensively-used regularization schemes which increase the solution time. We define a nucleation-completion criterion and modify the regularized solution algorithm. Implementation of the proposed integration algorithm within a user-defined subroutine UMAT in the commercial nonlinear finite element software ABAQUS has made possible the solution of boundary value problems. The obtained results show the efficiency of the proposed solution algorithm and confirm the improved efficiency (in terms of... 

Theoretical and numerical analyses of rotating disks with non-uniform thickness and material properties subjected to thermo-mechanical loadings have been carried out by variable material properties (VMP), Runge-Kutta's (RK) and finite element (FE) methods. The material is assumed to be elastic-linear hardening. A power form function is used to describe the temperature gradient with the higher temperature at outer surface. Von-Mises theory has been used as failure criterion. The effects of geometry, material and thermal loading parameters as well as boundary conditions on radial, hoop and equivalent stress distributions which have not been studied in much detail in previous works have been... 

This paper investigates implementation of upwind compact implicit and explicit high-order finite difference schemes for solution of the level set equation. The upwind compact implicit and explicit high-order finite difference schemes are well-known techniques to descritize spatial derivatives for convection term in hyperbolic equations. Applying of upwind high-order schemes on the level set equation leads to less error and CPU time reduction compared to essential non-oscillatory (ENO), weighted essential non-oscillatory schemes (WENO), and even different particle level set methods. The results indicate the error based on area loss decreases drastically with applying high-order upwind,... 

The Weibull distribution plays an important role in failure distribution modeling of reliability research. While there are three parameters in the general form of this distribution, for simplicity, one of its parameters is usually omitted and as a result, the others are estimated easily. However, due to its more flexibility, when the general form of the Weibull distribution is of interest, the estimation procedure is not an easy task anymore. For example, in the maximum likelihood estimation method, the likelihood function that is formed for a three-parameter Weibull distribution is very hard to maximize. In this paper, a new hybrid methodology based on a variable neighborhood search and a... 

Numerical approximation based on different forms of the governing partial differential equation can lead to significantly different results for two-phase flow in porous media. Selecting the proper primary variables is a critical step in efficiently modeling the highly nonlinear problem of multiphase subsurface flow. A comparison of various forms of numerical approximations for two-phase flow equations is performed in this work. Three forms of equations including the pressure-based, mixed pressure-saturation and modified pressure-saturation are examined. Each of these three highly nonlinear formulations is approximated using finite difference method and is linearized using both Picard and... 

In this paper, a technique for the study of nonlinear performance of structures in different levels of earthquakes is developed. In this method, the Incremental Dynamic Analysis (IDA) curves are not achieved from nonlinear dynamic analysis of multi-degree-of-freedom (MDF) structure. However, the procedure of constructing these curves is based on modelling of the entire structure with several single-degrees-of-freedom (SDF) structures and evaluating them through the modal pushover analysis method. An innovative idea for approximating pushover curves that is based on error distribution is introduced in this investigation. Furthermore, the total input energy applied towards the SDF oscillator,... 

In this paper, a multi-product multi-chance constraint joint single-vendor multi-buyers inventory problem is considered in which the demand follows a uniform distribution, the lead-time is assumed to vary linearly with respect to the lot size, and the shortage in combination of backorder and lost-sale is assumed. Furthermore, the orders are placed in multiple of packets, there is a limited space available for the vendor, there are chance constraints on the vendor service rate to supply the products, and there is a limited budget for each buyer to purchase the products. While the elements of the buyers' cost function are holding, shortage, order and transportation costs, the set up and... 

Overall turbine analysis requires large CPU time and computer memory, even in the present days. As a result, choosing an appropriate computational domain accompanied by a suitable boundary condition can dramatically reduce the time cost of computations. This work compares different geometries for numerical investigation of the 3D flow in the runner of a Francis turbine, and presents an optimum geometry with least computational effort and desirable numerical accuracy. The numerical results are validated with a GAMM Francis Turbine runner, which was used as a test case (GAMM workshop on 3D computation of incompressible internal flows, 1989) in which the geometry and detailed best efficiency...