• Steel and Composite Structures
• /
• v.51 no.3
• /
• pp.261-270
• /
• 2024

This paper proposes an effective method for optimizing the structure of functionally graded isotropic and incompressible linear elastic materials. The main emphasis is on utilizing a specialized polytopal composite finite element (PCE) technique capable of handling a broad range of materials, addressing common volumetric locking issues found in nearly incompressible substances. Additionally, it employs a continuum model for bi-directional functionally graded (BFG) material properties, amalgamating these aspects into a unified property function. This study thus provides an innovative approach that tackles diverse material challenges, accommodating various elemental shapes like triangles, quadrilaterals, and polygons across compressible and nearly incompressible material properties. The paper thoroughly details the mathematical formulations for optimizing the topology of BFG structures with various materials. Finally, it showcases the effectiveness and efficiency of the proposed method through numerous numerical examples.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.12 no.1
• /
• pp.13-23
• /
• 2008

In [7, 8], they proposed a new singular function(NSF) method to compute singular solutions of Poisson equations on a polygonal domain with re-entrant angles. Singularities are eliminated and only the regular part of the solution that is in $H^2$ is computed. The stress intensity factor and the solution can be computed as a post processing step. This method was extended to the interface problem and Poisson equations with the mixed boundary condition. In this paper, we give NSF method for the Helmholtz equations ${\Delta}u+Ku=f$ with homogeneous Dirichlet boundary condition. Examples with a singular point are given with numerical results.

• Smart Structures and Systems
• /
• v.18 no.1
• /
• pp.155-181
• /
• 2016

In the present paper we discuss the stability and the post-buckling behaviour of thin piezoelastic plates. The first part of the paper is concerned with the modelling of such plates. We discuss the constitutive modelling, starting with the three-dimensional constitutive relations within Voigt's linearized theory of piezoelasticity. Assuming a plane state of stress and a linear distribution of the strains with respect to the thickness of the thin plate, two-dimensional constitutive relations are obtained. The specific form of the linear thickness distribution of the strain is first derived within a fully geometrically nonlinear formulation, for which a Finite Element implementation is introduced. Then, a simplified theory based on the von Karman and Tsien kinematic assumption and the Berger approximation is introduced for simply supported plates with polygonal planform. The governing equations of this theory are solved using a Galerkin procedure and cast into a non-dimensional formulation. In the second part of the paper we discuss the stability and the post-buckling behaviour for single term and multi term solutions of the non-dimensional equations. Finally, numerical results are presented using the Finite Element implementation for the fully geometrically nonlinear theory. The results from the simplified von Karman and Tsien theory are then verified by a comparison with the numerical solutions.

• Proceedings of the Korean Society of Computational and Applied Mathematics Conference
• /
• 2003.09a
• /
• pp.9-9
• /
• 2003

Let $\Omega$ be a bounded, open, and polygonal domain in $R^2$ with re-entrant corners. We consider the following Partial Differential Equations: $$(I-\nabla\nabla\cdot+\nabla^{\bot}\nabla\times)u\;=\;f\;in\;\Omega$$, $$n\cdotu\;0\;0\;on\;{\Gamma}_{N}$$, $${\nabla}{\times}u\;=\;0\;on\;{\Gamma}_{N}$$, $$\tau{\cdot}u\;=\;0\;on\;{\Gamma}_{D}$$, $$\nabla{\cdot}u\;=\;0\;on\;{\Gamma}_{D}$$ where the symbol $\nabla\cdot$ and $\nabla$ stand for the divergence and gradient operators, respectively; $f{\in}L^2(\Omega)^2$ is a given vector function, $\partial\Omega=\Gamma_{D}\cup\Gamma_{N}$ is the partition of the boundary of $\Omega$; nis the outward unit vector normal to the boundary and $\tau$represents the unit vector tangent to the boundary oriented counterclockwise. For simplicity, assume that both $\Gamma_{D}$ and $\Gamma_{N}$ are nonempty. Denote the curl operator in $R^2$ by $$\nabla\times\;=\;(-{\partial}_2,{\partial}_1$$ and its formal adjoint by $${\nabla}^{\bot}\;=\;({-{\partial}_1}^{{\partial}_2}$$ Consider a weak formulation(WF): Find $u\;\in\;V$ such that $$a(u,v):=(u,v)+(\nabla{\cdot}u,\nabla{\cdot}v)+(\nabla{\times}u,\nabla{\times}V)=(f,v),\;A\;v{\in}V$$. (2) We assume there is only one singular corner. There are many methods to deal with the domain singularities. We introduce them shortly and we suggest a new Finite Element Methods by using Singular representation for the solution.

• Structural Monitoring and Maintenance
• /
• v.4 no.3
• /
• pp.195-220
• /
• 2017

Nowadays, railway system plays a significant role in transportation, conveying cargo, passengers, minerals, grains, and so forth. Railway ballasted track is a conventional railway track as can be seen all over the world. Ballast, located underneath the sleepers, is the most important elements on ballasted track, which has many functions and requires routine maintenance. Ballast needs to be maintained frequently to prevent rail buckling, settlement, misalignment so that ballast has to be modelled accurately. Continuum model was introduced to model granular material and was extended in ballast. However, ballast is a heterogeneous material with highly nonlinear behaviour. Hence, ballast could not be modelled accurately in continuum model due to the discontinuities nature and material degradation of ballast. Discrete element modelling (DEM) is proposed as an alternative approach that provides insight into constitutive model, realistic particle, and contact algorithm between each particle. DEM has been studied in many recent decades. However, there are limitations due to the high computational time and memory consumption, which cause the lack of using in high range. This paper presents a review of recent ballast modelling with benefits and drawbacks. Ballast particles are illustrated either circular, circular crump, spherical, spherical crump, super-quadric, polygonal and polyhedral. Moreover, the gaps and limitations of previous studies are also summarized. The outcome of this study will help the understanding into different ballast modelling and particle. The insight information can be used to improve ballast modelling and monitoring for condition-based track maintenance.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1998.03a
• /
• pp.205-209
• /
• 1998

Unlike the drawing of round section from round bar, the shaped drawing like polygonal section is known to have influence not only drawing stress but also comer filling. Therefore, this study analyze the drawing process of suqare rod from round bar using nonsteady state rigid-plastic FEM. To investigate effects of process variables of the drawing process of square rod from round bar, FE-simulations with variety of reduction in area and semi-die angle for a given frictional condition have been conduction. By this results, it has to suggest optimal process condition on the drawing stress and the comer filling. In addition, it has determined forming limit considering necking and bulging.

• East Asian mathematical journal
• /
• v.14 no.1
• /
• pp.63-76
• /
• 1998

we consider the hp-version to solve non-constant coefficients elliptic equations $-div(a{\nabla}u)=f$ with Dirichlet boundary conditions on a bounded polygonal domain $\Omega$ in $R^2$. In [6], M. Suri obtained an optimal error-estimate for the hp-version: ${\parallel}u-u^h_p{\parallel}_{1,\Omega}{\leq}Cp^{(\sigma-1)}h^{min(p,\sigma-1)}{\parallel}u{\parallel}_{\sigma,\Omega}$. This optimal result follows under the assumption that all integrations are performed exactly. In practice, the integrals are seldom computed exactly. The numerical quadrature rule scheme is needed to compute the integrals in the variational formulation of the discrete problem. In this paper we consider a family $G_p=\{I_m\}$ of numerical quadrature rules satisfying certain properties, which can be used for calculating the integrals. Under the numerical quadrature rules we will give the variational form of our non-constant coefficients elliptic problem and derive an error estimate of ${\parallel}u-\tilde{u}^h_p{\parallel}_{1,\Omega}$.

• Journal of the Korean Society for Precision Engineering
• /
• v.15 no.11
• /
• pp.145-151
• /
• 1998

Unlike the drawing of round section from round bar, the shaped drawing like polygonal section is known to have influence not only drawing stress but also corner filling. Therefore. this study analyze the drawing process of suqare rod from round bar using nonsteady state rigid-plastic FEM. To investigate effects of process variables of the drawing process of square rod from round bar, FE-simulations with variety of reduction in area and semi-die angle for a given frictional condition have been conducted. By this results, it has to suggest optimal process condition on the drawing stress and the corner filling. In addition, it has determined forming limit considering necking and bulging.

• International Journal of Ocean System Engineering
• /
• v.1 no.2
• /
• pp.95-101
• /
• 2011

Perforated plates with cutouts (or holes) are widely used in structural members. These cutouts provide stress concentration in plates. Extensive studies have been carried out on stress concentration in perforated plates, which consider cutout shapes, boundary conditions, bluntness of cutouts, and more. This study presents stress concentration analyses of perforated plates with not only various cutouts and bluntness but also different cutout orientations. Especially, the effect of cutout orientation on stress concentration is emphasized since structural members have become more complicated recently. To obtain stress concentration patterns, a finite element program, ANSYS, is used. For the designated goal, three parameters are considered as follows: the shapes of polygonal cutouts (circle, triangle, and square), bluntness (a counter measure of radius ratio, r/R), and rotation of cutouts (${\theta}$). From the analyses, it is shown that, in general, as bluntness increases, the stress concentration increases, regardless of the shape and rotation. A more important finding is that the stress concentration increases as the cutouts become more oriented from the baseline, which is the positive horizontal axis (+x). This fact demonstrates that the orientation is also a relatively significant design factor to reduce stress concentration. In detail, in the case of the triangle cutout, orienting one side of the triangle cutout to be perpendicular to the applied tensile forces is preferable. Similarly, in the case of the square cutout, it is more advantageous to orient two sides of square cutout to be perpendicular to the applied tensile force. Therefore, at the design stage, determining the direction of a major tensile force is required. Then, by aligning those polygon cutouts properly, we can reduce stress concentration.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.9 no.1
• /
• pp.56-63
• /
• 1985

Methods for evaluating the added masses of square and hexagonal structures with fluid gap are presented. For a sufficiently small fluid gap, an analytical expression for the added mass is found using the method of matched asymptotic expansion. Experimental data and numerical results using finite element method are also obtained for various sizes of fluid gap. It is shown that added masses increase in inverse proportion to the fluid gap as it becomes smaller. Experimental data, theoretical and numerical results are in good agreement.