• International Journal of Naval Architecture and Ocean Engineering
• /
• 제12권1호
• /
• pp.376-386
• /
• 2020

The aim of this study was to develop a new efficient strategy that uses the Vector form Intrinsic Finite-element (VFIFE) method to conduct the static and dynamic analyses of marine pipes. Nonlinear problems, such as large displacement, small strain, and contact and collision, can be analyzed using a unified calculation process in the VFIFE method according to the fundamental theories of point value description, path element, and reverse motion. This method enables analysis without the need to integrate the stiffness matrix of the structure, because only motion equations of particles established according to Newton's second law are required. These characteristics of the VFIFE facilitate the modeling and computation efficiencies in analyzing the nonlinear dynamic problem of flexible pipe with large deflections. In this study, a three-dimensional (3-D) dynamical model based on 3-D beam element was established according to the VFIFE method. The deep-sea flexible pipe was described by a set of spatial mass particles linked by 3-D beam element. The motion and configuration of the pipe are determined by these spatial particles. Based on this model, a simulation procedure to predict the 3-D dynamical behavior of flexible pipe was developed and verified. It was found that the spatial configuration and static internal force of the mining pipe can be obtained by calculating the stationary state of pipe motion. Using this simulation procedure, an analysis was conducted on the static and dynamic behaviors of the flexible mining pipe based on a 1000-m sea trial system. The results of the analysis proved that the VFIFE method can be efficiently applied to the static and dynamic analyses of marine pipes.

• 대한기계학회논문집A
• /
• 제30권1호
• /
• pp.66-75
• /
• 2006

Nonlinear response structural optimization using equivalent static loads (NROESL) has been proposed. Nonlinear response optimization is solved by sequential linear response optimization with equivalent static loads which are generated from the nonlinear responses and linear stiffness matrix. The linear stiffness matrix should be obtained in NROESL, and this process can be fairly difficult for some applications. Proportional transformation of loads (PTL) is proposed to overcome the difficulties. Equivalent static loads are obtained by PTL. It is the same as NROESL except for the process of calculating equivalent static loads. PTL is developed for large-scale probems. First, linear and nonlinear responses are evaluated from linear and nonlinear analyses, respectively. At a DOF of the finite element method, the ratio of the two responses is calculated and an equivalent static load is made by multiplying the ratio and the loads for linear analysis. Therefore, the mumber of the equivalent static loads is as many as that of DOF's and an equivalent static load is used with the reponse for the corresponding DOF in the optimization process. All the equivalent static loads are used as multiple loading conditions during linear response optimization. The process iterates until it converges. Examples are solved by using the proposed method and the results are compared with conventional methods.

• 해양환경안전학회지
• /
• 제27권1호
• /
• pp.193-200
• /
• 2021

본 연구에서는 국부 비선형 정적 해석을 효율적으로 수행하기 위하여 강성응축(Static condensation)을 활용한 해석기법을 제시하였다. 강성응축기법은 자유도 기반의 유한요소 모델 축소기법이며, 해석 모델을 관심 대상(Target) 부분과 응축되어 생략될(Omitted) 부분으로 구분한다. 본 연구에서는, 관심 대상 부분에는 비선형 영역, 생략될 부분에는 선형 영역으로 지정하였고, 선형 영역에 대응되는 강성 행렬 및 하중 벡터를 비선형 영역, 즉 관심 대상 부분으로 모두 응축하였다. 모델 응축 후에는 비선형 영역에 대한 강성 행렬 및 하중 벡터만으로 이루어진 축소 모델을 구성하였으며, 이 축소 모델만을 뉴턴-랩슨 반복(Newton-Raphson iteration)을 통해 갱신하여 효율적으로 비선형 해석을 수행하였다. 끝으로, 제안된 기법을 다양한 수치 예제에 적용하여 해석기법의 효율성과 신뢰성을 제시하였다.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2008년도 추계학술대회논문집
• /
• pp.71-81
• /
• 2008

This paper proposes a spectral element which can represent dynamic responses in high frequency domain such as Lamb waves on a thin plate. A two layer beam model under 2-D plane strain condition is introduced to simulate high-frequency dynamic responses induced by piezoelectric layer (PZT layer) bonded on a base plate. In the two layer beam model, a PZT layer is assumed to be rigidly bonded on a base beam. Mindlin-Herrmann and Timoshenko beam theories are employed to represent the first symmetric and anti-symmetric Lamb wave modes on a base plate, respectively. The Bernoulli beam theory and 1-D linear piezoelectricity are used to model the electro-mechanical behavior of a PZT layer. The equations of motions of a two layer beam model are derived through Hamilton's principle. The necessary boundary conditions associated with electro mechanical properties of a PZT layer are formulated in the context of dual functions of a PZT layer as an actuator and a sensor. General spectral shape functions of response field and the associated boundary conditions are formulated through equations of motions converted into frequency domain. A detailed spectrum element formulation for composing the dynamic stiffness matrix of a two layer beam model is presented as well. The validity of the proposed spectral element is demonstrated through comparison results with the conventional 2-D FEM and the previously developed spectral elements.

• 전산구조공학
• /
• 제9권1호
• /
• pp.33-45
• /
• 1996

전단변형을 고려한 보강재요소를 p-version 유한요소법을 사용하여 정식화 하였다. 적분형 르장드르 다항식으로부터 유도된 계층적 C/sup 0/-형상함수를 5자유도를 갖는 보강재와 평판요소의 조립강성도 행렬을 정의하는데 사용하였다. 보강재와 평판의 접속부에서 변위의 적합성을 만족시키기 위해 적절한 좌표변환행렬을 사용하여 국부좌표계에서 정의된 보강재의 강성도 행렬을 기준좌표계인 평판의 좌표계로 변환시켰다. 평판의 기준좌표계에 대한 보강재의 방향과 편심효과를 설명할 수 있는 변환행렬이 평판과 보강재의 접속부에서의 국부적인 거동과 합성구조로 된 보강판에서 평판과 보강재가 감당하는 상대적인 강도 분담을 파악하기 위해 사용되었다. p-version 유한요소법에 의한 결과를 기존의 연구결과와 비교하였으며, 특히 h-version유한요소해석 프로그램인 MICROFEAP-II의 결과를 비교하였다.

• 한국전산구조공학회논문집
• /
• 제30권2호
• /
• pp.127-135
• /
• 2017

멀티코어 CPU와 BLAS, LAPACK을 구현한 최적 수치라이브러리, 직접 희소 솔버의 대중화 등 PC나 워크스테이션 수준에서도 대규모 유한요소 모델을 해석할 수 있도록 컴퓨팅 환경이 급속도로 변화되었다. 이 논문에서는 멀티코어 CPU를 갖는 공유 메모리 구조에 대한 병렬 유한요소 프로그램 설계시 고려사항으로 (1) 최적화된 수치라이브러리의 사용, (2) 최신 직접 희소 솔버의 사용, (3) OpenMP를 이용한 병렬 요소 강성 행렬의 계산, (4) 희소행렬 저장방식의 일종인 triplet을 이용한 어셈블 기법 등을 제시하였다. 또한 대규모 수치모델을 통해 많은 시간이 소요되는 작업을 기준으로 병렬화 효과를 검토하였다.

• 한국소음진동공학회논문집
• /
• 제18권11호
• /
• pp.1157-1169
• /
• 2008

This paper presents spectral element formulation which approximates Lamb wave propagation by PZT transducers bonded on a thin plate. A two layer beam model under 2-D plane strain condition is introduced to simulate high-frequency dynamic responses induced by a piezoelectric (PZT) layer rigidly bonded on a base plate. Mindlin-Herrmann and Timoshenko beam theories are employed to represent the first symmetric and anti-symmetric Lamb wave modes on a base plate, respectively. The Euler-Bernoulli beam theory and 1-D linear piezoelectricity are used to model the electro-mechanical behavior of a PZT layer. The equations of motions of a two layer beam model are derived through Hamilton's principle. The necessary boundary conditions associated with the electro-mechanical properties of a PZT layer are formulated in the context of dual functions of a PZT layer as an actuator and a sensor. General spectral shape functions of response field and the associated boundary conditions are obtained through equations of motions converted into frequency domain. Detailed spectrum element formulation for composing the dynamic stiffness matrix of a two layer beam model is presented as well. The validity of the proposed spectral element is demonstrated through numerical examples.

• 한국지반공학회논문집
• /
• 제38권9호
• /
• pp.45-52
• /
• 2022

암반내 불연속면이나 지반-구조물의 접촉면은 열-수리-역학적으로 연계된 거동 특성을 보이므로, 온전한 지배방정식에 근거한 경계면 요소의 개발이 필요하다. 본 논문은 경계면 요소에 대한 힘평형 방정식, 유체의 연속방정식 그리고 에너지 평형 방정식을 유도하였다. 그리고 경계면 요소에 대한 탄소성 역학 모델의 강성행렬을 제시하였다. 개발된 유한요소는 2차원 조건에서 변위는 6절점, 수압과 온도는 4절점을 사용한다. 단층내 유체 주입에 대한 완전연계된 THM 해석은 단층내의 유효응력 감소와 주위 암반의 온도 수축에 의한 주입압의 복합적인 변화을 모델링 할수 있었다. 하지만, 열적 현상을 무시한 HM해석은 수리-역학적 변수를 과다하게 산정하였다.

• Journal of Mechanical Science and Technology
• /
• 제20권1호
• /
• pp.110-124
• /
• 2006

The multi scale wavelet-Galerkin method implemented in an adaptive manner has an advantage of obtaining accurate solutions with a substantially reduced number of interpolation points. The method is becoming popular, but its numerical efficiency still needs improvement. The objectives of this investigation are to present a new numerical scheme to improve the performance of the multi scale adaptive wavelet-Galerkin method and to give detailed implementation procedure. Specifically, the subdomain technique suitable for multiscale methods is developed and implemented. When the standard wavelet-Galerkin method is implemented without domain subdivision, the interaction between very long scale wavelets and very short scale wavelets leads to a poorly-sparse system matrix, which considerably worsens numerical efficiency for large-sized problems. The performance of the developed strategy is checked in terms of numerical costs such as the CPU time and memory size. Since the detailed implementation procedure including preprocessing and stiffness matrix construction is given, researchers having experiences in standard finite element implementation may be able to extend the multi scale method further or utilize some features of the multiscale method in their own applications.

• International Journal of CAD/CAM
• /
• 제12권1호
• /
• pp.20-28
• /
• 2012

Finite Element Analysis (FEA) is often used to identify local stress/strain concentrations where a component is likely to fail. In order to reduce the degree of strain concentration, component thickness can be increased in those regions, or a stronger material can be used. In short fiber reinforced composite materials, strength and stiffness can be increased through proper fiber alignment. The field-aided microtailoring (FAiMTa) process is one promising method for doing this. FAiMTa uses principles of dielectrophoresis to preferentially align particles or fibers within a matrix. To achieve the preferred fiber orientation, an interdigitated electrode network must be integrated into the mold halves which can be fabricated by additive manufacturing (AM) processes. However, the process of determining the preferred fiber arrangements and electrode locations can be very challenging. This paper presents algorithms to semi-automate the interdigitated electrode design process. The algorithm has been implemented in the Solidworks CAD system and is demonstrated in this paper.