H. Arzani, E. Khoshbavar Rad,
Volume 37, Issue 2 (3-2019)
Abstract
In this paper, a method is proposed to improve the results of the standard finite element method. L2 norm is used to determine the nodal error. In the next step, the appropriate order of the interpolation cover is seclected to be proportional to the nodal error and the results are corrected. The error computation procedure and the use of covering enrichment functions will continue until the error reaches the specified value. Cover enrichment interpolation functions will consider the effects of the adjacent elements of each node, in addition to the values obtained from the standard interpolation for each element. Computation rules are programmed in the matlab program and considered for the same examples. Comparison of the results of the proposed method with the exact solutions and the results of the methods proposed by the other researchers in the field of linear elasticity indicates the efficiency and accuracy of the proposed method.
S. Saravani, B. Keshtegar,
Volume 37, Issue 2 (3-2019)
Abstract
The computational burdens and more accurate approximations for the estimation of the failure probability are the main concerns in the structural reliability analyses. The Monte Carlo simulation (MCS) method can simply provide an accurate estimation for the failure probability, but it is a time-consuming method for complex reliability engineering problems with a low failure probability and may efficiently approximate the failure probability. In this paper, the efficiency of MCS for the computations of the performance function is improved using a random-weighted method known as the random-weighted MCS (RWMC) method. By using the weighted exponential function, the weights of random data points generated by MCS are adjusted by selecting the random point in the design space. The convergence performances including the computational burdens for evaluating the limit sate function and the accuracy of failure probabilities of RWMC are compared with MCS by using several nonlinear and complex mathematical and structural problems with normal and no-normal random variables. The results indicate that the proposed RWMC method can estimate the accurate results with the less computational burdens, about 100 to 1000 times faster than MCS
M. H. Sadeghi, S. Lotfan,
Volume 38, Issue 1 (8-2019)
Abstract
In this paper, nonlinear modal interactions caused by one-to-three internal resonance in a beam-mass-spring-damper system are investigated based on nonlinear system identification. For this purpose, the equations governing the transverse vibrations of the beam and mass are analyzed via the multiple scale method and the vibration response of the system under primary resonance is extracted. Then, the frequency behavior of the vibration response is studied by Fourier and Morlet wavelet transforms. In order to perform the nonparametric identification of the time response, mono-frequency intrinsic mode functions are derived by the advanced empirical mode decomposition. In this approach, masking signals are utilized in order to avoid mode mixing caused by modal interaction. After analyzing the frequency behavior of each mode function, slow flow dynamics of the system is established and intrinsic modal oscillators for reconstructing the corresponding intrinsic mode are extracted. Finally, by analyzing the beating phenomenon in a simple one-degree-of-freedom system, it is shown that the internal resonance causes beating only under the circumstance that the slope of the logarithmic amplitude of oscillator force is nonzero. The results, therefore, show that depending on the periodic, pseudo-periodic, and chaotic behavior of the response, modal interactions might be stationary or non-stationary. Moreover, the chaotic behavior occurs mostly in the vibration mode excited by the internal resonance mechanism
M. Khashei, F. Chahkoutahi,
Volume 38, Issue 1 (8-2019)
Abstract
Nowadays, electricity load forecasting, as one of the most important areas, plays a crucial role in the economic process. What separates electricity from other commodities is the impossibility of storing it on a large scale and cost-effective construction of new power generation and distribution plants. Also, the existence of seasonality, nonlinear complexity, and ambiguity pattern in electricity data set makes it more difficult to forecast by using the traditional methods. Therefore, new models, computational intelligence and soft computing tools and combining models are the most accurate and widely used methods for modeling the complexity and uncertainty in the data set. In this paper, a parallel optimal hybrid model using computational intelligence tools and soft computations is proposed to forecast the electricity load forecasting. The main idea of this model is the use of the advantages of the individual models in the modeling of complex systems in a structure and elimination of the limitations of them, simultaneously. The experimental results indicate that the proposed hybrid model has a higher performance accuracy in comparison to iterative suboptimal hybrid models and its computational cost is lower than the other hybrid models; also, the proposed model can achieve more accurate results, as compared with its component and some other seasonal hybrid models.
H. Bazai, A. Azari, M. Moshtagh,
Volume 38, Issue 1 (8-2019)
Abstract
The purpose of this article is the numerical study of flow and heat transfer characteristics of Nanofluids inside a cylindrical microchannel with rectangular, triangular, and circular cross-sections. The size and shape of these sections have a significant impact on the thermal and hydraulic performance of the microchannel heat exchanger. The Nanofluids used in this work include water and De-Ethylene Glycol (DEG) as the base fluids and Al2O3, Cu, SiO2 and CuO as the nanoparticles. To solve the problem and extract the required data, a 3-D simulation was performed for the microchannel using ANSYS FLUENT 15.0 software and the effect of the cross-sectional shape of the fluid flow and the type of nanoparticles on the thermal transfer and fluid flow parameters was studied. From the obtained results, it can be observed that the addition of nanoparticles to the base fluid increases the heat transfer and pressure drop. The results also show that rectangular channels have the best performance among the three geometries examined as its heat transfer coefficient was 19.26% higher than the triangular cross section which had the worst performance.
F. Kalateh, F. Hosseinejad,
Volume 38, Issue 1 (8-2019)
Abstract
Biot equations that consider fluid and soil interaction at the same time are the most applicable relationships in the soil dynamic analysis. However, in dynamic analysis, due to the sudden increase in the excess pore pressure caused by seismic excitation and the occurrence of high hydraulic gradients, the assumption of the Darcy flow used in these equations is questionable. In the present study, in the u-p form of Biot equations, non-Darcy flow is considered. Also, the nonlinear behavior of soil is modeled by the Pastor-Zienkiewicz -Chan model. For validation, the VELACS No.1 experiment is modeled and the effect of the nonlinear fluid flow assumption on the results is examined. The results indicate that in the low permeability coefficients, the obtained results of the non-Darcy and Darcy flow are in agreement; however, in high permeability coefficients, these two methods differ by time and depth.
I. Ahmadi, D. Kouhbor, R. Taghiloo,
Volume 38, Issue 1 (8-2019)
Abstract
In this paper, a finite element model is presented for the transient analysis of low velocity impact, and the impact induced damage in the composite plate subjected to low velocity impact is studied. The failure criteria suggested by Choi and Chang and the Tsai-Hill failure criteria are used for the prediction of the damage in the composite plate; then the effect of various parameters on the impact induced damage is investigated. The first order shear deformation plate theory and the Ritz finite element method are employed for modeling the behavior of plate, and the modified Hertz contact low is used for the prediction of the contact force through the impact. In the numerical results, the time history of indentation, contact force and stress during the impact and the impact induced damage is investigated. The matrix cracking and delamination in the plies of the laminated composite plate subjected to low velocity impact are studied and the effects of various parameters are investigated.
S. A. Ghazi Mirsaeed, V. Kalatjari,
Volume 38, Issue 1 (8-2019)
Abstract
In this paper, finite element analysis of thin viscoelastic plates is performed by proposing new plate elements using complex Fourier shape functions. New discrete Kirchhoff Fourier Theory (DKFT) plate elements are constructed by the enrichment of quadratic function fields in a six-noded triangular plate element with complex Fourier radial basis functions. In order to illustrate the validity and accuracy of the presented approach and robustness of the proposed elements in viscoelasticity, finite element analysis of square and elliptical viscoelastic thin plates is performed and the results are compared to those of analytical solutions and with those obtained by discrete Kirchhoff Theory (DKT) elements and the commercial software ABAQUS. The results show that FE solutions using DKFT elements have an excellent agreement with the analytical solutions and ABAQUS solutions, even though noticeably fewer elements, in comparison to DKT and classic plate elements, are employed, which means that the computational costs are reduced effectively.
A. Zamani Nouri, P. Ebrahimi,
Volume 38, Issue 2 (2-2020)
Abstract
With respect to the great application of pipes conveying fluid in civil engineering, presenting a mathematical model for their stability analysis is essential. For this purpose, a concrete pipe, reinforced by iron oxide (Fe2O3) nanoparticles, conveying fluid is considered. The goal of this study is to investigate the structural stability to show the effects of the inside fluid and the nanoparticles. The structure was modeled by a cylindrical shell and using Reddy theory. To obtain the force induced by the inside fluid, the Navier-Stokes equation was used. To assume the effect of the nanoparticles in the pipe, the Mori-Tanaka model was utilized so that the effects of agglomeration of nanoparticles could be considered. Finally, by applying energy method and the Hamilton's principle, the governing equations were derived. For the stability analysis of the structure, differential quadrature method (DQM) was proposed and the effects of different parameters such as volume fraction of the nanoparticles and agglomeration of the nanoparticles inside fluid and geometrical parameters were investigated. The results showed that the existence of the nanoparticles as the reinforcement for the pipe led to the delay in the pipe instability.
S. M. Navabi, M. Reisi-Nafchi, Gh. Moslehi,
Volume 38, Issue 2 (2-2020)
Abstract
Nowadays, outpatient providers are struggling to reduce the current costs and improve the service quality. A part of the outpatient service provider is a hemodialysis department with expensive supplies and equipment. Therefore, in the present paper, the scheduling of hemodialysis patients with their preferences has been studied. The aim of scheduling hemodialysis patients in this study is to minimize the normalized weighted sum of deviations from the patients' preferences and the total completion time. It should be noted that the patient's preferences include beds, treatment combination of days and their turn. To solve the problem, two mathematical models have been presented. Performence of the models in solving the real data of the hemodyalisis department of Imam Khomeini Hospital, in Kermanshah, was investigated. The results showed the efficiency of the proposed models in considering the preferences of patients; however, these preferences in the hospital schedule were considered in few cases, as far as it was possible. So, these preferences has no priority in the hospital schedule. In addition to considering the patients’ preferences, the solution of models reduced the total completion time of the pationts treatment. Also, one of the proposed models in this papercould optimally solve the instances three times larger than the hospital cases
Z. Shafiei, S. Sarrami-Foroushani, M. Azhari,
Volume 38, Issue 2 (2-2020)
Abstract
Graphene is one of the nanostructured materials that has recently attracted the attention of many researchers. This is due to the increasing expansion of nanotechnology and the application of this nanostructure in technology and industry owing to its mechanical, electrical and thermal properties. Thermal buckling behavior of single-layered graphene sheets is studied in this paper. Given the failure of classical theories to consider the scale effects and the limitations of the nano-sized experimental investigations of nano-materials, the small-scale effect is taken into account in this study, by employing the modified couple stress theory which has only one scale parameter. On the other hand, the two-variable refined plate theory, which considers the shear deformations in addition to bending deformations, is used to define the displacement field and to formulate the problem. The developed finite strip method formulation is used to evaluate the critical buckling temperature of the nanoplates. The validity of the proposed method is confirmed by comparing the results of this study with the those in the literature. The effects of different boundary conditions, temperature changing patterns, aspect ratio, and the ratio of length parameter to thickness on the critical buckling temperature are considered and the results are presented in the form of Tables and Figures
Z. Arefinia,
Volume 38, Issue 2 (2-2020)
Abstract
As thermalisation loss is the dominant loss process in the quantum dot intermediate band solar cells (QD-IBSCs), it has been investigated and calculated for a QD-IBSC, where IB is created by embedding a stack of InAs(1-x) Nx QDs with a square pyramid shape in the intrinsic layer of the AlPySb(1-y) p-i-n structure. IB, which is an optically coupled but electrically isolated mini-band, divides the total band gap of AlPySb(1-y) into two sub-band gaps. To obtain the thermalisation loss of AlPySb(1-y)/InAs(1-x)Nx QD-IBSCs, the position and width of IB in the band gap of AlPySb(1-y) should be calculated. The position of IB, which is equal to the first eigen-energy of a unit cell of QD, is obtained by solving the 3D Schrödinger equation with a finite-element method and the width of IB is obtained by the absorption characteristics. Then, with the investigation of the effect of nitrogen and phosphorous molar fraction, QDs size and the distance between the QDs on the thermalisation loss, the minimized loss for the optimized structure of AlPySb(1-y)/InAs(1-x)Nx QD-IBSCs is obtained
A. M. Salehizadeh, A. Shafiei,
Volume 38, Issue 2 (2-2020)
Abstract
This paper presents a numerical analysis of granular column collapse phenomenon using a two-dimensional smoothed particle hydrodynamics model and a local constitutive law proposed by Jop et al. This constitutive law, which is based on the viscoplastic behaviour of dense granular material flows, is characterized by an apparent viscosity depending both on the local strain rate and the local pressure. The rheological parameters are directly derived from the experiments. A simple proposed regularization method used in the viscosity relation to reproduce the stopping condition and the free surface of a granular flow where the pressure is disappeared. Pressure oscillation, as the main disadvantage of the weakly compressible SPH method, leads to an inaccurate pressure distribution. In this research, a new algorithm is proposed to remove the nonphysical oscillations by relating the divergence of velocity to the Laplacian of pressure. The simulations based on the proposed SPH algorithm satisfactorily capture the dynamics of gravity-driven granular flows observed in the experiments. The maximum thickness of a granular flowing on a rough inclined plane is obtained based on the local rheology model and compared with the experimental results. The run-out distances and the slopes of the deposits in the simulations showed a good agreement with the values found in the experiments. The results of the simulation proved that the initial column ratio played an important role in spreading the granular mass
F. Moradpouri,
Volume 38, Issue 2 (2-2020)
Abstract
Wave-field extrapolation based on solving the wave equation is an important step in seismic modeling and needs a high level of accuracy. It has been implemented through a various numerical methods such as finite difference method as the most popular and conventional one. Moreover, the main drawbacks of the finite difference method are the low level of accuracy and the numerical dispersion for large time intervals (∆t). On the other hand, the symplectic integrators due to their structure can cope with this problem and act more accurately in comparison to the finite difference method. They reduce the computation cost and do not face numerical dispersion when time interval is increased. Therefore, the aim of the current paper is to present a symplectic integrator for wave-field extrapolation using the Euler method. Then, the extrapolation is implemented for rather large time intervals using a simple geological model. The extrapolation employed for both symplectic Euler and finite difference methods showed a better quality image for the proposed method. Finally the accuracy was compared to the finite difference method
M. Jamei, H. R. Ghafouri,
Volume 38, Issue 2 (2-2020)
Abstract
In this study, we present a numerical solution for the two-phase incompressible flow in the porous media under isothermal condition using a hybrid of the linear lower-order nonconforming finite element and the interior penalty discontinuous Galerkin (DG) method. This hybridization is developed for the first time in the two-phase modeling and considered as the main novelty of this research.The pressure equation and convection dominant saturation equation are discretized using the nonconforming Crouziex-Raviart finite element (CR FEM) and the weighed interior penalty discontinuous Galerkin (SWIP) method, respectively. Utilizing the nonconforming finite element method for solving the flow equation made the pressure and velocity values be consistent with respect to the degrees of freedom arrangement at the midpoint of the neighboring element edges. The boundary condition governing the simulation is the Robin type at entrance boundaries, and the time marching discretization for the governing equations is the sequential solution scheme. An H (div) projection using Raviart-Thomas element is implemented to improve the results’ resolution and preserve the continuity of the normal component of the velocity field. At the end of each time step, the non-physical oscillation is omitted using a slope limiter, namely, modified Chavent-Jaffre limiter, in each element. Also, in this study, the developed algorithm is verified using some benchmark problems and the test cases are considered to demonstrate the efficiency of the developed model in capturing the shock front at the interface of fluid phases and discontinuities.
R. Ghiasi , M. R. Ghasemi ,
Volume 39, Issue 1 (8-2020)
Abstract
This paper focuses on the processing of structural health monitoring (SHM) big data. Extracted features of a structure are reduced using an optimization algorithm to find a minimal subset of salient features by removing noisy, irrelevant and redundant data. The PSO-Harmony algorithm is introduced for feature selection to enhance the capability of the proposed method for processing the measured big data, which have been collected from sensors of the structure and uncertainties associated with this process. Structural response signals under ambient vibration are preprocessed according to wavelet packet decomposition (WPD) and statistical characteristics for feature extraction. It optimizes feature vectors to be used as inputs to surrogate models based on the wavelet weighted support vector machine (WWLS-SVM) and radial basis function neural network (RBFNN). Two illustrative test examples are considered, the benchmark dataset from IASC-ASCE SHM group and a 120-bar dome truss. The results indicate that the features acquired by WPT from vibrational signal have higher sensitivity to the damage of the structure. Furthermore, the proposed PSO-Harmony is compared with four well-known metaheuristic optimization algorithms. The obtaind results show that the proposed method has a better performance and convergence rate. Finally, the proposed feature subset selection method has the capability of 90% data reduction
H. Mohammadiun, M. Mohammadiun, M. H. Dibaee Bonab, M. Darabi, S. R. Hejazi, V. Janipour Bidsardareh,
Volume 39, Issue 1 (8-2020)
Abstract
: In this research, dimensionless temperature and entropy generation for the steady state flow in the stagnation point of incompressible nanofluid impinging on an infinite cylinder have been investigated. The impinging free stream is steady with a constant strain rate k. Similarity solution of the Navier-Stokes equations and energy equation is derived in this problem. A reduction of these equations is obtained using appropriate transformations introduced in this research. The general self similar solution is obtained when the heat flux on the cylinder wall is constant. All solutions brought above are presented for Reynolds numbers Re=ka^2/2vf that range from 0.1 to 1000 and the selected values of particle fractions, where a is the radius of the cylinder and υf is the kinematic viscosity of the base fluid. Results show that for Reynolds numbers examined, as the particle fraction increases, the depth of diffusion of the fluid velocity field in axial direction decreases, whereas Nusselt number is raised. Also, the maximum value of entropy generation has been calculated.
M. Rezaee, M. M. Ettefagh, R. Fathi ,
Volume 39, Issue 1 (8-2020)
Abstract
Although the traditional automatic ball balancer (ABB) has numerous advantages, it has two major deficiencies, i.e., it has a limited balance stable region and it increases the vibration amplitude of the rotor at transient state. These deficiencies limit the applicability of ABBs. In this regard, a new type of ABB called “the Ball-spring autobalancer” has been proposed to resolve the mentioned deficiencies of the traditional ABBs. In order to investigate the capability of the Ball-spring AB in balancing rotors, it is necessary to study its dynamics accurately. The dynamics of a rotor with linear bearing equipped with a Ball-spring AB has been studied previously; however, in real situations, the bearings have nonlinear characteristics. Here, the dynamics of a rotor with nonlinear bearings equipped with a Ball-spring AB is investigated by the multiple scales method for the first time. The results show that the nonlinearity at the rotor bearings does not impair the advantages of the Ball-spring AB.
N. Cheraghi, M. Miri, M. Rashki,
Volume 39, Issue 1 (8-2020)
Abstract
This paper presents a probabilistic assessment on the free vibration analysis of functionally graded material plates, including layers with magneto-electro-elastic properties, using the 3D solution and surrogate models. The plate is located on an elastic foundation and the intra-layer slipping effect is also considered in the analysis by employing the generalized intra-layer spring model. Due to the high computational cost of the 3D solution in calculating the free vibration frequency of the plate, surrogate models are used. The meta models including kriging method, radial fundamental function method and polynomial response surface method are used to construct the surrogate model. For surrogate models training, the results of the three-dimensional solving method are used. The elastic foundation hardness, the intra-layer slipping effect, the material properties index, and the layer densities are considered as the variables with uncertainty. The three-dimensional solution method is validated through a comparison with other available reference. The results obtained through the surrogate models have been compared to those of the 3D solution formulation, showing a good agreement. The effects of some parameters including the elastic foundation hardness, the intra-layer slipping effect, the density of each layer, and the material properties index on the fundamental frequency of functionally graded material plates are investigated. By using three-dimensional solution method and Kriging Surrogate Model, it is shown that the shear and transverse components of elastic foundation hardness and the density of each layer have the greatest effect on the fundamental frequency of the functionally graded material plates.
S. A. Ahmadi, M. H. Pashaei, R. A. Jafari-Talookolaeilokoolaei,
Volume 39, Issue 1 (8-2020)
Abstract
In this paper, three-dimensional displacement response of a cylindrical sandwich panel with compressible core under the action of dynamic pulse loading is addressed using the extended high order sandwich panel theory. Also, local dynamic pulse buckling of facesheets is studied by considering the Budiansky-Roth buckling criterion. It is assumed that the sandwich panels consist of orthotropic face sheets and an isotropic viscoelastic foam core layer. The effects of various parameters including the panel span, core and facing thickness, pulse duration and maximum pressure on the non-linear dynamic response and buckling strength of the sandwich cylindrical panel are studied. The results obtained from the present method are compared with finite element solutions using the commercial software ANSYS and those reported in the literature, showing a good agreement. It is revealed that applied core non-linear theory could be satisfactory for the dynamic pulse response of sandwich viscoelastic panels. It is also shown that the pulse buckling strength of panel increases with a decrease of the panel radius or an increase of the panel thickness.
