• Journal of applied mathematics & informatics
• /
• v.30 no.5_6
• /
• pp.951-969
• /
• 2012

In this paper, we propose a class of meshless collocation approaches for the solution of time dependent partial differential algebraic equations (PDAEs) in terms of a radial basis function interpolation numerical scheme. Kansa's method and the Hermite collocation method (HCM) for PDAEs are given. A sensitivity analysis of the solutions from different shape parameter c is obtained by numerical experiments. With use of the random collocation points, we have obtain the more accurate solution by the methods than those by the finite difference method for the PDAEs with index-2, i.e, we avoid the influence from an index jump of PDAEs in some degree. Several numerical experiments show that the methods are efficient.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.12 s.255
• /
• pp.1626-1633
• /
• 2006

In the present paper, in the context of the meshless interpolation of a moving least squares (MLS) type, a novel method which uses primary and secondary nodes in the domain and on the global boundary is introduced, in order to improve the accuracy of solution. The secondary nodes can be placed at any location where one needs to obtain a better resolution. The support domains for the shape functions in the MLS approximation are defined from the primary nodes, and the secondary nodes use the same support domains. The shape functions based on the MLS approximation, in an integration domain, have a single type of a rational function, which reduces the difficulty of numerical integration to evaluate the weak form. The present method is very useful in an adaptive calculation, because the secondary nodes can be easily added and moved without an additional mesh. Several numerical examples are presented to illustrate the effectiveness of the present method.

• Structural Engineering and Mechanics
• /
• v.11 no.2
• /
• pp.221-236
• /
• 2001

A local point interpolation method (LPIM) is presented for the stress analysis of two-dimensional solids. A local weak form is developed using the weighted residual method locally in two-dimensional solids. The polynomial interpolation, which is based only on a group of arbitrarily distributed nodes, is used to obtain shape functions. The LPIM equations are derived, based on the local weak form and point interpolation. Since the shape functions possess the Kronecker delta function property, the essential boundary condition can be implemented with ease as in the conventional finite element method (FEM). The presented LPIM method is a truly meshless method, as it does not need any element or mesh for both field interpolation and background integration. The implementation procedure is as simple as strong form formulation methods. The LPIM has been coded in FORTRAN. The validity and efficiency of the present LPIM formulation are demonstrated through example problems. It is found that the present LPIM is very easy to implement, and very robust for obtaining displacements and stresses of desired accuracy in solids.

• Structural Engineering and Mechanics
• /
• v.62 no.1
• /
• pp.43-54
• /
• 2017

In this study, the failure behavior of composite material in the biaxial and off-axis loading is studied based on a computational micromechanical model. The model is developed so that the combination of mechanical and thermal loading conditions can be considered in the analysis. The modified generalized plane strain assumption of the theory of elasticity is used for formulation of the micromechanical modeling of the problem. A truly meshless method is employed to solve the governing equation and predict the distribution of micro-stresses in the selected RVE of composite. The fiber matrix interface is assumed to be perfect until the interface failure occurs. The biaxial and off-axis loading of the SiC/Ti and Kevlar/Epoxy composite is studied. The failure envelopes of SiC/Ti and Kevlar/Epoxy composite in off-axis loading, biaxial transverse-transverse and axial-transverse loading are predicted based on the micromechanical approach. Various failure criteria are considered for fiber, matrix and fiber-matrix interface. Comparison of results with the available results in the litreture shows excellent agreement with experimental studies.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2001.04a
• /
• pp.77-80
• /
• 2001

The element-free Galerkin (EFG) method is one of meshless methods, which is an efficient method of modeling problems of fluid or solid mechanics with complex boundary shapes and large changes in boundary conditions. This paper discusses the theory of the EFG method and its applications to modeling of groundwater flow. In the EFG method, shape functions are constructed based on the moving least square (MLS) approximation, which requires only set of nodes. The EFG method can eliminate time-consuming mesh generation procedure with irregular shaped boundaries because it does not require any elements. The coupled EFG-FEM technique was introduced to treat Dirichlet boundary conditions. A computer code EFGG was developed and tested for the problems of steady-state and transient groundwater flow in homogeneous or heterogeneous aquifers. The accuracy of solutions by the EFG method was similar to that by the FEM. The EFG method has the advantages in convenient node generation and flexible boundary condition implementation.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.35-42
• /
• 2002

In this paper, effective analysis of beam is studied using the RKPM in meshless methods. So, RKPM is extended for solving moderately thick and thin beam. General Timoshenko beam theory is used for formulation. Shear locking is the main difficulty in analysis of beam structures. The shear relaxation factor and corrected shear rigidity are introduced to overcome shear locking. Analysis results obtained reveal that RKPM using introduced methods Is free of locking and very effectively applicable to deeply as well as shallowly beams.

• Structural Engineering and Mechanics
• /
• v.11 no.2
• /
• pp.123-132
• /
• 2001

An improved version of the Element-free Galerkin method (EFGM) is presented here for addressing the problem of transverse shear locking in shear-deformable beams with a high length over thickness ratio. Based upon Timoshenko's theory of thick beams, it has been recognized that shear locking will be completely eliminated if the rotation field is constructed to match the field of slope, given by the first derivative of displacement. This criterion is applied directly to the most commonly implemented version of EFGM. However in the numerical process to integrate strain energy, the second derivative of the standard Moving Least Square (MLS) shape functions must be evaluated, thus requiring at least a $C^1$ continuity of MLS shape functions instead of $C^0$ continuity in the conventional EFGM. Yet this hindrance is overcome effortlessly by only using at least a $C^1$ weight function. One-dimensional quartic spline weight function with $C^2$ continuity is therefore adopted for this purpose. Various numerical results in this work indicate that the modified version of the EFGM does not exhibit transverse shear locking, reduces stress oscillations, produces fast convergence, and provides a surprisingly high degree of accuracy even with coarse domain discretizations.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.4
• /
• pp.26-34
• /
• 2003

In this paper, a novel technique is proposed to arbitrarily activate the nodal points in finite element model through the meshless approximation methods such as MLS(moving least squares method), and theoretical investigations are carried out including the consistency and boundeness of numerical solution to prove the validity of the proposed method. By using the proposed node activation technique, one can activate and handle only the concerned nodes as unknown variables among the large number of nodal points in the finite element model. Therefore, the proposed technique has a great potential in design and reanalysis procedure.

• International Journal of CAD/CAM
• /
• v.8 no.1
• /
• pp.1-10
• /
• 2009

This paper presents a topology optimization approach using element-free Galerkin method (EFGM) for the optimal design of compliant mechanisms with geometrically non-linearity. Meshless method has an advantage over the finite element method(FEM) because it is more capable of handling large deformation resulted from geometrical nonlinearity. Therefore, in this paper, EFGM is employed to discretize the governing equations and the bulk density field. The sensitivity analysis of the optimization problem is performed by incorporating the adjoint approach with the meshless method. The Lagrange multipliers method adjusted for imposition of both the concentrated and continuous essential boundary conditions in the EFGM is proposed in details. The optimization mathematical formulation is developed to convert the multi-criteria problem to an equivalent single-objective problem. The popularly applied interpolation scheme, solid isotropic material with penalization (SIMP), is used to indicate the dependence of material property upon on pseudo densities discretized to the integration points. A well studied numerical example has been applied to demonstrate the proposed approach works very well and the non-linear EFGM can obtain the better topologies than the linear EFGM to design large-displacement compliant mechanisms.

• Structural Engineering and Mechanics
• /
• v.88 no.3
• /
• pp.239-249
• /
• 2023

The Radial Point Interpolation Method (RPIM) has been proposed to overcome the difficulties associated with the use of the Radial Basis Functions (RBFs). The RPIM has the following properties: Simple implementation in terms of boundary conditions as in the Finite Element Method (FEM). A less expensive CPU time compared to other collocation meshless methods such as the Moving Least Square (MLS) collocation method. In this work, we propose an adaptive high-order numerical algorithm based on RPIM to simulate the thermoviscoplastic behavior of a material mixing observed in the Friction Stir Welding (FSW) process. The proposed adaptive meshfree RPIM algorithm adapts well to the geometric and physical data by choosing a good shape parameter with a good precision. Our numerical approach combines the RPIM and the Asymptotic Numerical Method (ANM). A numerical procedure is also proposed in this work to automatically determine an improved shape parameter for the RBFs. The efficiency of the proposed algorithm is analyzed in comparison with an iterative algorithm.