• Structural Engineering and Mechanics
• /
• 제17권1호
• /
• pp.1-14
• /
• 2004

Numerical solution to linear bending analysis of circular plates is obtained by the method of harmonic differential quadrature (HDQ). In the method of differential quadrature (DQ), partial space derivatives of a function appearing in a differential equation are approximated by means of a polynomial expressed as the weighted linear sum of the function values at a preselected grid of discrete points. The method of HDQ that was used in the paper proposes a very simple algebraic formula to determine the weighting coefficients required by differential quadrature approximation without restricting the choice of mesh grids. Applying this concept to the governing differential equation of circular plate gives a set of linear simultaneous equations. Bending moments, stresses values in radial and tangential directions and vertical deflections are found for two different types of load. In the present study, the axisymmetric bending behavior is considered. Both the clamped and the simply supported edges are considered as boundary conditions. The obtained results are compared with existing solutions available from analytical and other numerical results such as finite elements and finite differences methods. A comparison between the HDQ results and the finite difference solutions for one example plate problem is also made. The method presented gives accurate results and is computationally efficient.

• Nuclear Engineering and Technology
• /
• 제54권7호
• /
• pp.2625-2634
• /
• 2022

Discretization errors are extremely challenging conundrums of discrete ordinates calculations for radiation transport problems with void regions. In previous work, we have presented a multi-collision source method (MCS) to overcome discretization errors, but the efficiency needs to be improved. This paper proposes a goal-oriented algorithm for the MCS method to adaptively determine the partitioning of the geometry and dynamically change the angular quadrature in remaining iterations. The importance factor based on the adjoint transport calculation obtains the response function to get a problem-dependent, goal-oriented spatial decomposition. The difference in the scalar fluxes from one high-order quadrature set to a lower one provides the error estimation as a driving force behind the dynamic quadrature. The goal-oriented algorithm allows optimizing by using ray-tracing technology or high-order quadrature sets in the first few iterations and arranging the integration order of the remaining iterations from high to low. The algorithm has been implemented in the 3D transport code ARES and was tested on the Kobayashi benchmarks. The numerical results show a reduction in computation time on these problems for the same desired level of accuracy as compared to the standard ARES code, and it has clear advantages over the traditional MCS method in solving radiation transport problems with reflective boundary conditions.

• Geomechanics and Engineering
• /
• 제29권6호
• /
• pp.667-681
• /
• 2022

This paper presents a thermoelastic analysis of variable thickness plates made of functionally graded materials (FGM) subjected to mechanical and thermal loads. The thermal load is applied to the plate as a temperature difference between the top and bottom surfaces. Temperature distribution in the plate is obtained using the steady-state heat equation. Except for Poisson's ratio, all mechanical properties of the plate are assumed to vary linearly along the thickness direction based on the volume fractions of ceramic and metal. The plate is resting on an elastic foundation modeled based on the Winkler foundation model. The governing equations are derived based on the third-order shear deformation theory (TSDT) and are solved numerically for various boundary conditions using the differential quadrature method (DQM). The effects of various parameters on the stress distribution and deflection of the plate are investigated such as the value of thermal and mechanical loads, volume fractions of ceramic and metal, and the stiffness coefficients of the foundation.

• Wind and Structures
• /
• 제24권5호
• /
• pp.431-446
• /
• 2017

As concrete is most usable material in construction industry it's been required to improve its quality. Nowadays, nanotechnology offers the possibility of great advances in construction. In this study, buckling of horizontal concrete columns reinforced with Zinc Oxide (ZnO) nanoparticles is analyzed. Due to the presence of ZnO nanoparticles which have piezoelectric properties, the structure is subjected to electric field for intelligent control. The Column is located in foundation with vertical springs and shear modulus constants. Sinusoidal shear deformation beam theory (SSDBT) is applied to model the structure mathematically. Micro-electro-mechanic model is utilized for obtaining the equivalent properties of system. Using the nonlinear stress-strain relation, energy method and Hamilton's principal, the motion equations are derived. The buckling load of the column is calculated by Difference quadrature method (DQM). The aim of this study is presenting a mathematical model to obtain the buckling load of structure as well as investigating the effect of nanotechnology and electric filed on the buckling behavior of structure. The results indicate that the negative external voltage applied to the structure, increases the stiffness and the buckling load of column. In addition, reinforcing the structure by ZnO nanoparticles, the buckling load of column is increased.

• Journal of applied mathematics & informatics
• /
• 제29권3_4호
• /
• pp.813-829
• /
• 2011

In this article, we present a numerical scheme for solving singularly perturbed (i.e. highest -order derivative term multiplied by small parameter) Burgers-Huxley equation with appropriate initial and boundary conditions. Most of the traditional methods fail to capture the effect of layer behavior when small parameter tends to zero. The presence of perturbation parameter and nonlinearity in the problem leads to severe difficulties in the solution approximation. To overcome such difficulties the present numerical scheme is constructed. In construction of the numerical scheme, the first step is the dicretization of the time variable using forward difference formula with constant step length. Then, the resulting non linear singularly perturbed semidiscrete problem is linearized using quasi-linearization process. Finally, differential quadrature method is used for space discretization. The error estimate and convergence of the numerical scheme is discussed. A set of numerical experiment is carried out in support of the developed scheme.

• 전산구조공학
• /
• 제10권2호
• /
• pp.139-148
• /
• 1997

A finite element method is programmed to analyse the nonlinear behavior of axisymmetric structures. The lst order Mindlin shell theory which takes into account the transversal shear deformation is used to formulate a conical two node element with six degrees of freedom. To evade the shear locking phenomenon which arises in Mindlin type element when the effect of shear deformation tends to zero, the reduced integration of one point Gauss Quadrature at the center of element is employed. This method is the Updated Lagrangian formulation which refers the variables to the state of the most recent iteration. The solution is searched by Newton-Raphson iteration method. The tangent matrix of this method is obtained by a finite difference method by perturbating the degrees of freedom with small values. For the moment this program is limited to the analyses of non-linear elastic problems. For structures which could have elastic stability problem, the calculation is controled by displacement.

• 한국전산구조공학회논문집
• /
• 제25권2호
• /
• pp.139-148
• /
• 2012

MLS(Moving Least Squares) 차분법은 무요소법의 이동최소제곱법과 Taylor 전개를 이용하여 요소망의 제약 및 수치 적분이 없이 절점만을 이용하여 미분방정식을 수치해석할 수 있는 방법이다. 본 연구에서는 고체역학 문제의 동적해석을 위하여 MLS 차분법의 시간이력해석 알고리즘을 제시한다. 개발된 알고리즘은 Newmark 방법으로 시간적분을 하였으며, 강형식을 그대로 이산화하여 해석을 수행했다. 이동최소제곱법을 이용해 Taylor 전개식을 근사하여 실제 미분계산없이 미분근사식을 얻기 때문에 고차까지 Taylor 다항식의 차수를 증가하는 것이 용이하다. 1차원과 2차원 수치예제들을 통하여 동적해석을 위한 MLS 차분법의 정확성과 효율성을 검증하였다. 수치결과들이 정확해에 잘 수렴하였으며, 유한요소법(FEM)의 해석결과와 비교하여 떨림현상(oscillation) 및 주기성(periodicity) 오차에 대해 보다 안정적인 모습을 보였다.

• 대한수학회논문집
• /
• 제37권3호
• /
• pp.939-955
• /
• 2022

In this paper, we study a first-order non-linear singularly perturbed Volterra integro-differential equation (SPVIDE). We discretize the problem by a uniform difference scheme on a Bakhvalov-Shishkin mesh. The scheme is constructed by the method of integral identities with exponential basis functions and integral terms are handled with interpolating quadrature rules with remainder terms. An effective quasi-linearization technique is employed for the algorithm. We establish the error estimates and demonstrate that the scheme on Bakhvalov-Shishkin mesh is O(N^-1) uniformly convergent, where N is the mesh parameter. The numerical results on a couple of examples are also provided to confirm the theoretical analysis.

• Computers and Concrete
• /
• 제25권6호
• /
• pp.537-549
• /
• 2020

Geometrically nonlinear axisymmetric bending analysis of shear deformable circular plates on a nonlinear three-parameter elastic foundation was made. Plates ranging from "thin" to "moderately thick" were investigated for three types of material: isotropic, transversely isotropic, and orthotropic. The differential equations were discretized by means of the finite difference method (FDM) and the differential quadrature method (DQM). The Newton-Raphson method was applied to find the solution. A parametric investigation using seven unknowns per node was presented. The novelty of the paper is that detailed numerical simulations were made to highlight the combined effects of the material properties and the boundary conditions on (i) the deflection, (ii) the stress resultants, and (iii) the external load. The formulation was verified through comparison studies. It was observed that the results are highly influenced from the boundary conditions, and from the material properties.

• 대한토목학회논문집
• /
• 제29권5A호
• /
• pp.457-465
• /
• 2009

본 연구는 균열선단에서 응력특이성을 갖는 탄성균열문제를 해석하기 위한 이동최소제곱 유한차분법을 제시한다. 응력특이성을 유발하는 균열선단 주변장을 모형화하기 위해 근사식에 선단주변함수를 내재적으로 도입하여 이동최소제곱 근사의 틀을 그대로 유지하면서 실제 미분계산을 거의 하지 않고 미분근사를 할 수 있는 이동최소제곱 Taylor 다항식 근사의 장점을 살렸다. 균열문제 정식화시 시간소모적인 적분과정이 필요한 약정식화 대신 해석영역에 배치된 절점에서 지배 미분방정식에 대한 차분식을 직접 구성하는 강정식화를 적용하여 계산 효율성을 향상시켰다. 균열문제 해석을 통해 내적확장된 이동최소제곱 유한차분법이 응력 특이성을 내포한 선단주변 변위장을 정확히 묘사할 수 있을 뿐만 아니라 응력확대계수를 정확히 계산 할 수 있음을 보였다.