• 제목/요약/키워드: natural element method (NEM)

검색결과 17건 처리시간 0.019초

효율적이고 신뢰성있는 자연요소 균열해석을 위한 균열선단 그리드 세분화기법 (A Near-tip Grid Refinement for the Effective and Reliable Crack Analysis by Natural Element Method)

  • 조진래
    • 한국전산구조공학회논문집
    • /
    • 제32권3호
    • /
    • pp.183-190
    • /
    • 2019
  • 본 논문은 균열선단 그리드 세분화기법을 소개하고 자연요소법을 이용한 균열해석에 이 기법을 적용함으로서 그 유효성을 고찰하였다. 유한요소법에 있어서의 국부적 h-세분화와 같이 높은 응력 특이성을 보이는 균열선단 주위를 따라 자연요소법 그리드를 국부적으로 세분화하였다. 본 논문에서 소개되는 그리드 세분화기법은 2단계로 구성되며, 1단계에서는 그리드 점들이 추가되고 2단계에서는 균열선단 절점을 공유하는 델라우니 삼각형들이 나뉘게 된다. 제안하는 그리드 세분화기법의 타당성과 균열해석에서의 유효성을 입증하기 위해 대칭 엣지 균열을 갖는 평면 변형률 상태의 사각 평판을 해석하였다. 수치해석 결과의 상대비교를 위해 균일한 자연요소 그리드를 이용한 균열해석도 수행하였으며, 균열선단이 세분화된 그리드는 균일한 그리드와는 달리 이론해와 조밀한 그리드와 유사한 균열선단 응력분포를 나타내었다. 또한, 총 그리드 절점수에 대한 해석결과의 전역 상대오차에서도 세분화된 그리드가 균일한 그리드에 비해 높은 수렴율 나타내었다.

Nonlinear bending analysis of functionally graded CNT-reinforced composite plates

  • Cho, Jin-Rae
    • Steel and Composite Structures
    • /
    • 제42권1호
    • /
    • pp.23-32
    • /
    • 2022
  • In this paper, a nonlinear numerical method to solve the large deflection problem is introduced. And the non-dimensional load-deflection behavior of functionally graded carbon nanotube-reinforced composite (FG-CNTRC) plates is parametrically investigated. The large deflection problem is formulated according to the von Kármán nonlinear theory and the (1,1,0)* hierarchical model, and it is approximated by 2-D natural element method (NEM). The shear locking phenomenon is suppressed by the selectively reduced integration method. The nonlinear matrix equations are solved by combining the incremental loading scheme and the Newton-Raphson iteration method. The proposed method is validated from the benchmark experiments, where the propose method shows an excellent agreement with the reference methods. The nonlinear behavior of FG-CNTRC plates is evaluated in terms of the non-dimensional load-deflection curve, and it is parametrically investigated with respect to the existence/non-existence and gradient pattern of CNTs, the width-to-thickness and aspect ratios of plates and the type of boundary conditions. The non-dimensional central deflection is significantly reduced when CNTs and added, and it decreases with the volume fraction of CNTs. But, it shows a uniform increase in proportion to the width-to-thickness and aspect ratios. Both the gradient pattern of CNTs and the type of boundary conditions do also show the remarkable effects.

Free Vibration and Dynamic Response Analysis by Petrov-Galerkin Natural Element Method

  • Cho, Jin-Rae;Lee, Hong-Woo
    • Journal of Mechanical Science and Technology
    • /
    • 제20권11호
    • /
    • pp.1881-1890
    • /
    • 2006
  • In this paper, a Petrov-Galerkin natural element method (PG-NEM) based upon the natural neighbor concept is presented for the free vibration and dynamic response analyses of two-dimensional linear elastic structures. A problem domain is discretized with a finite number of nodes and the trial basis functions are defined with the help of the Voronoi diagram. Meanwhile, the test basis functions are supported by Delaunay triangles for the accurate and easy numerical integration with the conventional Gauss quadrature rule. The numerical accuracy and stability of the proposed method are verified through illustrative numerical tests.

Neutral surface-based static and free vibration analysis of functionally graded porous plates

  • J.R. Cho
    • Steel and Composite Structures
    • /
    • 제49권4호
    • /
    • pp.431-440
    • /
    • 2023
  • The functionally graded (FG) porous plates are usually characterized by the non-symmetric elastic modulus distribution through the thickness so that the plate neutral surface does not coincide with the mid-surface. Nevertheless, the conventional analysis models were mostly based on the plate mid-surface so that the accuracy of resulting numerical results is questionable. In this context, this paper presents the neutral surface-based static and free vibration analysis of FG porous plates and investigates the differences between the mid- and neutral surface-based analysis models. The neutral surface-based numerical method is formulated using the (3,3,2) hierarchical model and approximated by the last introduced natural element method (NEM). The volume fractions of metal and ceramic are expressed by the power-law function and the cosine-type porosity distributions are considered. The proposed numerical method is demonstrated through the benchmark experiment, and the differences between two analysis models are parametrically investigated with respect to the thickness-wise material and porosity distributions. It is found from the numerical results that the difference cannot be negligible when the material and porosity distributions are remarkably biased in the thickness direction.

Strain recovery-based equilibrated transverse shear stresses in functionally graded shell-like structures

  • Jin-Rae Cho
    • Structural Engineering and Mechanics
    • /
    • 제91권5호
    • /
    • pp.527-538
    • /
    • 2024
  • The standard numerical approximation of structural displacement field leads to the thickness-wise transverse shear stress distributions which are quite different from the exact ones. To overcome this inherent problem, an effective and reliable post-processing method is presented based on the strain recovery and the stress equilibrium, particularly for functionally graded cylindrical and conical elastic panels. The present method is developed in the framework of locking-free 2-D natural element method. Through the recovery of displacement component-wise derivatives, the element-wise discontinuous in-plane strain distributions are enhanced to be globally continuous and smoothened. And, using the continuous in-plane strains, the troublesome poor transverse shear stress distributions are enhanced through the thickness-wise integration of static equilibrium equations. The validity of present post-processing method is verified through the comparison with the reference solutions. In addition, the comparative experiments are also performed to investigate the difference between the present method and other available post-processing methods. The numerical results confirm that the present method provides the accurate transverse shear stress distributions which are consistent with the reference solutions and much better than other available methods.

Investigation of nonlinear free vibration of FG-CNTRC cylindrical panels resting on elastic foundation

  • J.R. Cho
    • Structural Engineering and Mechanics
    • /
    • 제88권5호
    • /
    • pp.439-449
    • /
    • 2023
  • Non-linear vibration characteristics of functionally graded CNT-reinforced composite (FG-CNTRC) cylindrical shell panel on elastic foundation have not been sufficiently examined. In this situation, this study aims at the profound numerical investigation of the non-linear vibration response of FG-CNTRC cylindrical panels on Winkler-Pasternak foundation by introducing an accurate and effective 2-D meshfree-based non-linear numerical method. The large-amplitude free vibration problem is formulated according to the first-order shear deformation theory (FSDT) with the von Karman non-linearity, and it is approximated by Laplace interpolation functions in 2-D natural element method (NEM) and a non-linear partial derivative operator HNL. The complex and painstaking numerical derivation on the curved surface and the crucial shear locking are overcome by adopting the geometry transformation and the MITC3+ shell elements. The derived nonlinear modal equations are iteratively solved by introducing a three-step iterative solving technique which is combined with Lanczos transformation and Jacobi iteration. The developed non-linear numerical method is estimated through the benchmark test, and the effects of foundation stiffness, CNT volume fraction and functionally graded pattern, panel dimensions and boundary condition on the non-linear vibration of FG-CNTRC cylindrical panels on elastic foundation are parametrically investigated.

Memory Organization for a Fuzzy Controller.

  • Jee, K.D.S.;Poluzzi, R.;Russo, B.
    • 한국지능시스템학회:학술대회논문집
    • /
    • 한국퍼지및지능시스템학회 1993년도 Fifth International Fuzzy Systems Association World Congress 93
    • /
    • pp.1041-1043
    • /
    • 1993
  • Fuzzy logic based Control Theory has gained much interest in the industrial world, thanks to its ability to formalize and solve in a very natural way many problems that are very difficult to quantify at an analytical level. This paper shows a solution for treating membership function inside hardware circuits. The proposed hardware structure optimizes the memoried size by using particular form of the vectorial representation. The process of memorizing fuzzy sets, i.e. their membership function, has always been one of the more problematic issues for the hardware implementation, due to the quite large memory space that is needed. To simplify such an implementation, it is commonly [1,2,8,9,10,11] used to limit the membership functions either to those having triangular or trapezoidal shape, or pre-definite shape. These kinds of functions are able to cover a large spectrum of applications with a limited usage of memory, since they can be memorized by specifying very few parameters ( ight, base, critical points, etc.). This however results in a loss of computational power due to computation on the medium points. A solution to this problem is obtained by discretizing the universe of discourse U, i.e. by fixing a finite number of points and memorizing the value of the membership functions on such points [3,10,14,15]. Such a solution provides a satisfying computational speed, a very high precision of definitions and gives the users the opportunity to choose membership functions of any shape. However, a significant memory waste can as well be registered. It is indeed possible that for each of the given fuzzy sets many elements of the universe of discourse have a membership value equal to zero. It has also been noticed that almost in all cases common points among fuzzy sets, i.e. points with non null membership values are very few. More specifically, in many applications, for each element u of U, there exists at most three fuzzy sets for which the membership value is ot null [3,5,6,7,12,13]. Our proposal is based on such hypotheses. Moreover, we use a technique that even though it does not restrict the shapes of membership functions, it reduces strongly the computational time for the membership values and optimizes the function memorization. In figure 1 it is represented a term set whose characteristics are common for fuzzy controllers and to which we will refer in the following. The above term set has a universe of discourse with 128 elements (so to have a good resolution), 8 fuzzy sets that describe the term set, 32 levels of discretization for the membership values. Clearly, the number of bits necessary for the given specifications are 5 for 32 truth levels, 3 for 8 membership functions and 7 for 128 levels of resolution. The memory depth is given by the dimension of the universe of the discourse (128 in our case) and it will be represented by the memory rows. The length of a world of memory is defined by: Length = nem (dm(m)+dm(fm) Where: fm is the maximum number of non null values in every element of the universe of the discourse, dm(m) is the dimension of the values of the membership function m, dm(fm) is the dimension of the word to represent the index of the highest membership function. In our case then Length=24. The memory dimension is therefore 128*24 bits. If we had chosen to memorize all values of the membership functions we would have needed to memorize on each memory row the membership value of each element. Fuzzy sets word dimension is 8*5 bits. Therefore, the dimension of the memory would have been 128*40 bits. Coherently with our hypothesis, in fig. 1 each element of universe of the discourse has a non null membership value on at most three fuzzy sets. Focusing on the elements 32,64,96 of the universe of discourse, they will be memorized as follows: The computation of the rule weights is done by comparing those bits that represent the index of the membership function, with the word of the program memor . The output bus of the Program Memory (μCOD), is given as input a comparator (Combinatory Net). If the index is equal to the bus value then one of the non null weight derives from the rule and it is produced as output, otherwise the output is zero (fig. 2). It is clear, that the memory dimension of the antecedent is in this way reduced since only non null values are memorized. Moreover, the time performance of the system is equivalent to the performance of a system using vectorial memorization of all weights. The dimensioning of the word is influenced by some parameters of the input variable. The most important parameter is the maximum number membership functions (nfm) having a non null value in each element of the universe of discourse. From our study in the field of fuzzy system, we see that typically nfm 3 and there are at most 16 membership function. At any rate, such a value can be increased up to the physical dimensional limit of the antecedent memory. A less important role n the optimization process of the word dimension is played by the number of membership functions defined for each linguistic term. The table below shows the request word dimension as a function of such parameters and compares our proposed method with the method of vectorial memorization[10]. Summing up, the characteristics of our method are: Users are not restricted to membership functions with specific shapes. The number of the fuzzy sets and the resolution of the vertical axis have a very small influence in increasing memory space. Weight computations are done by combinatorial network and therefore the time performance of the system is equivalent to the one of the vectorial method. The number of non null membership values on any element of the universe of discourse is limited. Such a constraint is usually non very restrictive since many controllers obtain a good precision with only three non null weights. The method here briefly described has been adopted by our group in the design of an optimized version of the coprocessor described in [10].

  • PDF