• 대한조선학회논문집
• /
• 제51권5호
• /
• pp.409-418
• /
• 2014

Simulation of flow past a complex marine structure requires a fine resolution in the vicinity of the structure, whereas a coarse resolution is enough far away from it. Therefore, a lot of grid cells may be wasted, when a simple Cartesian grid system is used for an Immersed Boundary Method (IBM). To alleviate this problems while maintaining the Cartesian frame work, we adopted an Adaptive Mesh Refinement (AMR) scheme where the grid system dynamically and locally refines as needed. In this study, We implemented a moving IBM and an AMR technique in our basic 3D incompressible Navier-Stokes solver. A Volume Of Fluid (VOF) method was used to effectively treat the free surface, and a recently developed Lagrangian Dynamic Subgrid-scale Model (LDSM) was incorporated in the code for accurate turbulence modeling. To capture vortex induced vibration accurately, the equation for the structure movement and the governing equations for fluid flow were solved at the same time implicitly. Also, We have developed an interface by using AutoLISP, which can properly distribute marker particles for IBM, compute the geometrical information of the object, and transfer it to the solver for the main simulation. To verify our numerical methodology, our results were compared with other authors' numerical and experimental results for the benchmark problems, revealing excellent agreement. Using the verified code, we investigated the following cases. (1) simulating flow around a floating sphere. (2) simulating flow past a marine structure.

• Structural Engineering and Mechanics
• /
• 제24권6호
• /
• pp.647-664
• /
• 2006

In this study, an improved finite element formulation with a scheme of solution for the large deflection analysis of inextensible prismatic and nonprismatic slender beams is developed. For this purpose, a three-noded Lagrangian beam-element with two dependent degrees of freedom per node (i.e., the vertical displacement, y, and the actual slope, $dy/ds=sin{\theta}$, where s is the curved coordinate along the deflected beam) is used to derive the element stiffness matrix. The element stiffness matrix in the global xy-coordinate system is achieved by means of coordinate transformation of a highly nonlinear ($6{\times}6$) element matrix in the local sy-coordinate. Because of bending with large curvature, highly nonlinear expressions are developed within the global stiffness matrix. To achieve the solution after specifying the proper loading and boundary conditions, an iterative quasi-linearization technique with successive corrections are employed considering these nonlinear expressions to remain constant during all iterations of the solution. In order to verify the validity and the accuracy of this study, the vertical and the horizontal displacements of prismatic and nonprismatic beams subjected to various cases of loading and boundary conditions are evaluated and compared with analytic solutions and numerical results by available references and the results by ADINA, and excellent agreements were achieved. The main advantage of the present technique is that the solution is directly obtained, i.e., non-incremental approach, using few iterations (3 to 6 iterations) and without the need to split the stiffness matrix into elastic and geometric matrices.

• Structural Engineering and Mechanics
• /
• 제15권1호
• /
• pp.83-114
• /
• 2003

A new simple 3-node triangular flat-shell element with standard nodal DOF (6 DOF per node) is proposed for the linear and geometrically nonlinear analysis of very thin to thick plate and shell structures. The formulation of element GT9 (Long and Xu 1994), a generalized conforming membrane element with rigid rotational freedoms, is employed as the membrane component of the new shell element. Both one-point reduced integration scheme and a corresponding stabilization matrix are adopted for avoiding membrane locking and hourglass phenomenon. The bending component of the new element comes from a new generalized conforming Kirchhoff-Mindlin plate element TSL-T9, which is derived in this paper based on semiLoof constrains and rational shear interpolation. Thus the convergence can be guaranteed and no shear locking will happen. Furthermore, a simple hybrid procedure is suggested to improve the stress solutions, and the Updated Lagrangian formulae are also established for the geometrically nonlinear problems. Numerical results with solutions, which are solved by some other recent element models and the models in the commercial finite element software ABAQUS, are presented. They show that the proposed element, denoted as GMST18, exhibits excellent and better performance for the analysis of thin-think plates and shells in both linear and geometrically nonlinear problems.

• KSII Transactions on Internet and Information Systems (TIIS)
• /
• 제12권5호
• /
• pp.2233-2252
• /
• 2018

For the problem of cross-layer joint resource allocation (JRA) in the Long-Term Evolution (LTE)-Advanced standard using carrier aggregation (CA) technology, it is difficult to obtain the optimal resource allocation scheme. This paper proposes a joint resource allocation algorithm based on the weights of user's average quality of experience (JRA-WQOE). In contrast to prevalent algorithms, the proposed method can satisfy the carrier aggregation abilities of different users and consider user fairness. An optimization model is established by considering the user quality of experience (QoE) with the aim of maximizing the total user rate. In this model, user QoE is quantified by the mean opinion score (MOS) model, where the average MOS value of users is defined as the weight factor of the optimization model. The JRA-WQOE algorithm consists of the iteration of two algorithms, a component carrier (CC) and resource block (RB) allocation algorithm called DABC-CCRBA and a subgradient power allocation algorithm called SPA. The former is used to dynamically allocate CC and RB for users with different carrier aggregation capacities, and the latter, which is based on the Lagrangian dual method, is used to optimize the power allocation process. Simulation results showed that the proposed JRA-WQOE algorithm has low computational complexity and fast convergence. Compared with existing algorithms, it affords obvious advantages such as improving the average throughput and fairness to users. With varying numbers of users and signal-to-noise ratios (SNRs), the proposed algorithm achieved higher average QoE values than prevalent algorithms.

• 한국전산구조공학회논문집
• /
• 제27권1호
• /
• pp.9-16
• /
• 2014

레벨셋 기법과 위상민감도를 이용하여 선형 탄성 구조물에 대하여, 초기 설계형상에 의존성이 없는 위상 및 형상 최적설계 기법을 개발하였다. 레벨셋 기법에서는 복잡한 위상 형상변화를 쉽게 다루기 위해 초기 영역은 고정한 채 레벨셋 함수로 표현되는 암시적 이동경계로 경계를 표현한다. 해밀턴-자코비(H-J) 방정식과 수치적으로 강건한 기법인 'up-wind scheme'은 컴플라이언스 목적함수를 최소화시키고 허용체적 제약조건을 만족시키면서, 초기 암시적 경계를 법선 속도장에 따라 최적의 형상으로 이끌어 낸다. 점근적인 정규화 개념에 근거하여, 구멍의 반지름을 0으로 접근시켜 형상 미분의 극한을 취한 위상민감도를 고려하였다. 최적조건으로부터 유도된 라그란지안의 감소 방향을 이용하여 H-J 방정식을 갱신하기 위한 속도장을 결정하였다. 개발한 방법에서는 위상민감도로부터 얻어지는 지표를 이용하여 구멍을 언제든지 어디에서나 생성가능하기 때문에 초기 구멍이 최적 형상을 얻기 위해 요구되지 않는다는 사실을 확인하였다. 또한 효율적인 최적화 과정을 위해서는 구멍 생성을 위한 조정변수의 적절한 선택이 중요함을 확인하였다.

• 한국농공학회지
• /
• 제39권4호
• /
• pp.106-113
• /
• 1997

Dynamic loading of structures often causes excursions of stresses well into the inelastic range, and the influence of the geometric changes on the dynamic response is also significant in many cases. Therefore, both material and geometric nonlinearity effects should be considered in case that a dynamic load acts on the structure. A structure in a nuclear power plant is a structure of importance which puts emphasis on safety. A nuclear container is a pressure vessel subject to internal pressure and this structure is constructed by a reinforced concrete or a pre-stressed concrete. In this study, the material nonlinearity effect on the dynamic response is formulated by the elasto-viscoplastic model highly corresponding to the real behavior of the material. Also, the geometrically nonlinear behavior is taken into account using a total Lagrangian coordinate system, and the equilibrium equation of motion is numerically solved by a central difference scheme. The constitutive relation of concrete is modeled according to a Drucker-Prager yield criterion in compression. The reinforcing bars are modeled by a smeared layer at the location of reinforcements, and the steel layer model under Von Mises yield criteria is adopted to represent an elastic-plastic behavior. To investigate the dynamic response of a nuclear reinforced concrete containment structure, the steel-ratios of 0, 3, 5 and 10 percent, are considered. The results obtained from the analysis of an example were summarized as follows 1. As the steel-ratio increases, the amplitude and the period of the vertical displacements in apex of dome decreased. The Dynamic Magnification Factor(DMF) was some larger than that of the structure without steel. However, the regular trend was not found in the values of DMF. 2. The dynamic response of the vertical displacement and the radial displacement in the dome-wall junction were shown that the period of displacement in initial step decreased with the steel-ratio increases. Especially, the effect of the steel on the dynamic response of radial displacement disapeared almost. The values of DMF were 1.94, 2.5, 2.62 and 2.66, and the values increased with the steel-ratio. 3. The characteristics of the dynamic response of radial displacement in the mid-wall were similar to that of dome-wall junction. The values of DMF were 1.91, 2.11, 2.13 and 2.18, and the values increased with the steel-ratio. 4. The amplitude and the period of the hoop-stresses in the dome, the dome-wall junction, and the mid-wall were shown the decreased trend with the steel-ratio. The values of DMF were some larger than those of the structure without steel. However, the regular trend was not found in the values of DMF.

• 지구물리와물리탐사
• /
• 제11권1호
• /
• pp.68-77
• /
• 2008

2.5차원 전자탐사 적분방정식의 확장된 Born 근사해 또는 국소 비선형 근사에 기초하여 루프-루프 전자탐사 역산 알고리듬이 개발되었다 송수신 배열은 수평 동일면(HCP) 또는 수직 동일면(VCP) 방식이고, 다중 주파수 및 다중 송수신 간격을 포함할 수 있으며 PC에서 작동된다. 안정적이고 고해상도를 유지하는 역산이 가능하도록 변수분해 행렬과 Backus-Gilbert 분산 함수 분석을 통해 감도 분포의 함수로서의 공간적으로 변화하는 최적 Lagrange 곱수 결정 알고리듬을 포함하였다. HCP와 VCP 배열 자료가 지하 전기비저항 구조에 따라 서로 다른 감도를 가짐에 따라 동시 역산에서 안정성과 해상도에 영향을 미치게 되므로, 계산값과 측정값 차의 분산에 따라 가중치를 적용하는 방식을 도입하였다. 모델링 코드의 정확성은 통상적으로 루프-루프 전자탐사에서 사용하는 주파수 및 송수신 간격 범위에서 유한차분법에 의해 계산된 결과와의 비교를 통하여 증명되었다. 개발된 역산 알고리듬은 먼저 반무한 공간내 전도체 및 저항체가 포함된 모델에 대한 계산자료에 적용되어 성능이 입증되었다. 현장자료에 적용하고 그 결과 영상을 전기비저항 탐사자료에 대한 역산 결과와 비교하여, 의미있는 지하구조의 영상을 얻을 수 있음을 확인하였다.

• 지구물리와물리탐사
• /
• 제6권4호
• /
• pp.199-206
• /
• 2003

MT 탐사자료의 역산에 있어서 지하의 전기비저항과 함께 정적효과를 파라미터로 설정하여 동시 역산을 수행하는 알고리듬을 하나의 지하구조 모델에 각기 다른 양의 정적효과를 포함시킨 4개의 자료에 대하여 적용시키고 이를 정적효과가 전혀 고려되지 않은 경우와 비교하여, 3차원 역산에서 정적효과가 미치는 영향 및 그 특성에 대하여 분석하였다. 일반적으로 현장자료에 정적효과가 어느 정도 포함되어 있는지에 대한 사전 정보가 전혀 없으므로 역산과정에서 이를 조절하는 trade-off 파라미터의 적절한 선택이 매우 중요하며, 본 연구에서는 모델의 smoothness와 static shift의 양을 조절하는 각각의 파라미터의 크기를 매 반복마다 구하는 알고리듬을 동시역산에 적용하였으며 4개의 이론자료에 적용한 결과 만족할 만한 결과를 얻었다. 정적효과가 포함된 자료에 대하여 정적효과를 고려하지 않은 역산(기존의 MT 역산)에서는 지표 block의 전기비저항을 바꿔 역산 스스로가 정적효과를 유발하려는 경향을 보였으며 이의 결과로 저주파수에서는 상당한 정적효과를 발생시켜 정적효과가 그리 크지 않은 경우 심부구조를 어느 정도 규명해 내는 것으로 나타났다. 그러나 고주파수에서는 이들 지표 block의 영향이 주파수에 무관하지 않게 되어 정적효과를 포함하는 자료의 겉보기 전기비저항과 위상을 동시에 만족시키지 못하게 된다. 그러나 정적효과를 파라미터로 하는 동시역산의 경우, 매우 심한 정적효과를 포함하는 자료에 대해서도 지하구조를 매우 정확히 영상화 하는 것이 가능하였다.