• Structural Engineering and Mechanics
• /
• v.62 no.6
• /
• pp.759-769
• /
• 2017

The effect of central axial load on natural frequencies of various thin-walled beams, are investigated by some researchers using different methods such as finite element, transfer matrix and dynamic stiffness matrix methods. However, there are situations that the load will be off centre. This type of loading is called eccentric load. The effect of the eccentricity of axial load on the natural frequencies of asymmetric thin-walled beams is a subject that has not been investigated so far. In this paper, the mentioned effect is studied using exact dynamic stiffness matrix method. Flexure and torsion of the aforesaid thin-walled beam is based on the Bernoulli-Euler and Vlasov theories, respectively. Therefore, the intended thin-walled beam has flexural rigidity, saint-venant torsional rigidity and warping rigidity. In this paper, the Hamilton‟s principle is used for deriving governing partial differential equations of motion and force boundary conditions. Throughout the process, the uniform distribution of mass in the member is accounted for exactly and thus necessitates the solution of a transcendental eigenvalue problem. This is accomplished using the Wittrick-Williams algorithm. Finally, in order to verify the accuracy of the presented theory, the numerical solutions are given and compared with the results that are available in the literature and finite element solutions using ABAQUS software.

• Coupled systems mechanics
• /
• v.7 no.1
• /
• pp.95-109
• /
• 2018

The relevance of turbulent mixing in estuarine numerical models for stratified two-layer shallow water flows is analysed in this paper. A one-dimensional numerical model was developed for this purpose by extending an immiscible two-layer model with an additional source term, which accounts for turbulent mixing effects, namely the entrainment of fluid from the lower to the upper layer. The entrainment rate is quantified by an empirical equation as a function of the bulk Richardson number. A finite volume method based on an approximated Roe solver was used to solve the governing coupled system of partial differential equations. A comparison of numerical results with and without entrainment is presented to illustrate the influence of entrainment on both the salt-water intrusion length and lower layer dynamics. Furthermore, one example is given to demonstrate how entrainment terms may help to stabilize the numerical scheme and prevent a possible loss of hyperbolicity. Finally, the model with entrainment is validated by comparing the numerical results to field measurements.

• Coupled systems mechanics
• /
• v.8 no.2
• /
• pp.129-146
• /
• 2019

The main goal of this work is to develop a numerical simulator to study an isothermal single-phase two-component flow in a naturally fractured oil reservoir, taking into account advection and diffusion effects. We use the Peng-Robinson equation of state with a volume translation to evaluate the properties of the components, and the discretization of the governing partial differential equations is carried out using the Finite Difference Method, along with implicit and first-order upwind schemes. This process leads to a coupled non-linear algebraic system for the unknowns pressure and molar fractions. After a linearization and the use of an operator splitting, the Conjugate Gradient and Bi-conjugated Gradient Stabilized methods are then used to solve two algebraic subsystems, one for the pressure and another for the molar fraction. We studied the effects of fractures in both the flow field and mass transport, as well as in computing time, and the results show that the fractures affect, as expected, the flow creating a thin preferential path for the mass transport.

• Journal of the Korean Physical Society
• /
• v.73 no.11
• /
• pp.1774-1786
• /
• 2018

Phenotypic variation among clones (individuals with identical genes, i.e. isogenic individuals) has been recognized both theoretically and experimentally. We investigate the effects of phenotypic variation on evolutionary dynamics of a population. In a population, the individuals are assumed to be haploid with two genotypes : one genotype shows phenotypic variation and the other does not. We use an individual-based Moran model in which the individuals reproduce according to their fitness values and die at random. The evolutionary dynamics of an individual-based model is formulated in terms of a master equation and is approximated as the Fokker-Planck equation (FPE) and the coupled non-linear stochastic differential equations (SDEs) with multiplicative noise. We first analyze the deterministic part of the SDEs to obtain the fixed points and determine the stability of each fixed point. We find that there is a discrete phase transition in the population distribution when the probability of reproducing the fitter individual is equal to the critical value determined by the stability of the fixed points. Next, we take demographic stochasticity into account and analyze the FPE by eliminating the fast variable to reduce the coupled two-variable FPE to the single-variable FPE. We derive a quasi-stationary distribution of the reduced FPE and predict the fixation probabilities and the mean fixation times to absorbing states. We also carry out numerical simulations in the form of the Gillespie algorithm and find that the results of simulations are consistent with the analytic predictions.

• Coupled systems mechanics
• /
• v.11 no.4
• /
• pp.357-371
• /
• 2022

This paper presents a parametric study on the free vibration analysis of a functionally graded material (FGM) circular plate with non-uniform thickness resting on a variable Pasternak elastic foundation. The mechanical properties of the material vary in the transverse direction through the thickness of the plate according to the power-law distribution to represent the constituent components. The equation of motion of the circular plate has been carried out based on the classical plate theory (CPT), and the differential quadrature method (DQM) is employed to solve the governing equations as a semi-analytical method. The grid points are chosen based on Chebyshev-Gauss-Lobatto distribution to achieve acceptable convergence and better accuracy. The influence of geometric parameters, variable elastic foundation, and functionally graded variation for clamped and simply supported boundary conditions on the first three natural frequencies are investigated. Comparisons of results with similar studies in the literature have been presented and two-dimensional mode shapes for particular plates have been plotted to illustrate the effect of variable thickness profile.

• Steel and Composite Structures
• /
• v.47 no.5
• /
• pp.569-585
• /
• 2023

This paper proposes an efficient approach for the structural topology optimization of bi-directional functionally graded structures by incorporating popular radial basis functions (RBFs) into an implicit level set (ILS) method. Compared to traditional element density-based methods, a level set (LS) description of material boundaries produces a smoother boundary description of the design. The paper develops RBF implicit modeling with multiquadric (MQ) splines, thin-plate spline (TPS), exponential spline (ES), and Gaussians (GS) to define the ILS function with high accuracy and smoothness. The optimization problem is formulated by considering RBF-based nodal densities as design variables and minimizing the compliance objective function. A LS-RBF optimization method is proposed to transform a Hamilton-Jacobi partial differential equation (PDE) into a system of coupled non-linear ordinary differential equations (ODEs) over the entire design domain using a collocation formulation of the method of lines design variables. The paper presents detailed mathematical expressions for BiDFG beams topology optimization with two different material models: continuum functionally graded (CFG) and mechanical functionally graded (MFG). Several numerical examples are presented to verify the method's efficiency, reliability, and success in accuracy, convergence speed, and insensitivity to initial designs in the topology optimization of two-dimensional (2D) structures. Overall, the paper presents a novel and efficient approach to topology optimization that can handle bi-directional functionally graded structures with complex geometries.

• Steel and Composite Structures
• /
• v.45 no.3
• /
• pp.305-330
• /
• 2022

This article aims to investigate the static deflection and stress analysis of bi-directional functionally graded porous plate (BDFGPP) modeled by unified higher order kinematic theories to include the shear stress effects, which not be considered before. Different shear functions are described according to higher order models that satisfy the zero-shear influence at the top and bottom surfaces, and hence refrain from the need of shear correction factor. The material properties are graded through two spatial directions (i.e., thickness and length directions) according to the power law distribution. The porosities and voids inside the material constituent are described by different cosine functions. Hamilton's principle is implemented to derive the governing equilibrium equation of bi-directional FG porous plate structures. An efficient numerical differential integral quadrature method (DIQM) is exploited to solve the coupled variable coefficients partial differential equations of equilibrium. Problem validation and verification have been proven with previous prestigious work. Numerical results are illustrated to present the significant impacts of kinematic shear relations, gradation indices through thickness and length, porosity type, and boundary conditions on the static deflection and stress distribution of BDFGP plate. The proposed model is efficient in design and analysis of many applications used in nuclear, mechanical, aerospace, naval, dental, and medical fields.

• Steel and Composite Structures
• /
• v.50 no.4
• /
• pp.429-441
• /
• 2024

This paper studies the free vibration behavior of trapezoidal shaped coupled double-layered graphene sheets (DLGS) system using first-order shear deformation theory (FSDT) and incorporating nonlocal elasticity theory. Two nanoplates are assumed to be bonded by an interlayer van der walls force and surrounded by an external kelvin-voight viscoelastic medium. The governing equations together with related boundary condition are discretized using a mapping-differential quadrature method (DQM) in the spatial domain. Then the natural frequency of the system is obtained by solving the eigen value matrix equation. The validity of the current study is evaluated by comparing its numerical results with those available in the literature and then a parametric study is thoroughly performed, concentrating on the series effects of angles and aspect ratio of GS, viscoelastic medium, and nonlocal parameter. The model is used to study the vibration of DLGS for two typical deformation modes, the in-phase and out-of-phase vibrations, which are investigated. Numerical results indicate that due to Increasing the damping parameter of the viscoelastic medium has reduced the frequency of both modes and this medium has been able to overdamped the oscillations and by increasing stiffness parameters both in-phase and out-of-phase vibration frequencies increased.

• The Bulletin of The Korean Astronomical Society
• /
• v.40 no.1
• /
• pp.49.2-49.2
• /
• 2015

IntraCluster Medium (ICM) located at the galaxy cluster is in the state of very hot, tenuous, magnetized, and highly ionized X-ray emitting plasmas. High temperature and low density make ICM very viscous and conductive. In addition to the high conductivity, fluctuating random plasma motions in ICM, occurring at all evolution stages, generate and amplify the magnetic fields in such viscous ionized gas. The amplified magnetic fields in reverse drive and constrain the plasma motions beyond the viscous scale through the magnetic tension. Moreover, without the influence of resistivity viscous damping effect gets balanced only with the magnetic tension in the extended viscous scale leading to peculiar ICM energy spectra. This overall collisionless magnetohydrodynamic (MHD) turbulence in ICM was simulated using a hyper diffusivity method. The results show the plasma motions and frozen magnetic fields have power law of $E_V^k{\sim}k^{-3}$, $E_M^k{\sim}k^{-1}$. To explain these abnormal power spectra we set up two simultaneous differential equations for the kinetic and magnetic energy using an Eddy Damped Quasi Normal Markovianized (EDQNM) approximation. The solutions and dimensions of leading terms in the coupled equations derive the power spectra and tell us how the spectra are formed. We also derived the same results with a more intuitive balance relation and stationary energy transport rate.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.20 no.3
• /
• pp.243-259
• /
• 2016

The Richards equation for water movement in unsaturated soil is highly nonlinear partial differential equations which are not solvable analytically unless unrealistic and oversimplifying assumptions are made regarding the attributes, dynamics, and properties of the physical systems. Therefore, conventionally, numerical solutions are the only feasible procedures to model flow in partially saturated porous media. The standard Finite element numerical technique is usually coupled with an Euler time discretizations scheme. Except for the fully explicit forward method, any other Euler time-marching algorithm generates nonlinear algebraic equations which should be solved using iterative procedures such as Newton and Picard iterations. In this study, lumped mass and distributed mass in the frame of Picard and Newton iterative techniques were evaluated to determine the most efficient method to solve the Richards equation with finite element model. The accuracy and computational efficiency of the scheme and of the Picard and Newton models are assessed for three test problems simulating one-dimensional flow processes in unsaturated porous media. Results demonstrated that, the conventional mass distributed finite element method suffers from numerical oscillations at the wetting front, especially for very dry initial conditions. Even though small mesh sizes are applied for all the test problems, it is shown that the traditional mass-distributed scheme can still generate an incorrect response due to the highly nonlinear properties of water flow in unsaturated soil and cause numerical oscillation. On the other hand, non oscillatory solutions are obtained and non-physics solutions for these problems are evaded by using the mass-lumped finite element method.