• 한국지구과학회지
• /
• 제35권5호
• /
• pp.333-341
• /
• 2014

A new set of basis functions was constructed using the Rossby-Haurwitz waves, which are the eigenfunctions of nondivergent barotropic vorticity equations on the sphere. The basis functions were designed to be non-separable, that is, not factored into functions of either the longitude or the latitude. Due to this property, the nodal lines of the functions are aligned neither along with the meridian nor the parallel. The basis functions can be categorized into groups of which members have the same degree or the total wavenumber-like index on the sphere. The orthonormality of the basis functions were found to be close to the machine roundoffs, giving the error of $O(10^{-15})$ or $O(10^{-16})$ for double-precision computation (64 bit arithmetic). It was demonstrated through time-stepping procedure that the basis functions were also the eigenfunctions of the non-divergent barotropic vorticity equations. The projection of the basis functions was carried out onto the low-resolution geopotential field of Gaussian bell, and compared with the theory. The same projections were performed for the observed atmospheric-geopotential height field of 500 hPa surface to demonstrate decomposition into the fields that contain disturbance of certain range of horizontal scales. The usefulness of the new basis functions was thus addressed for application to the eigenmode analysis of the atmospheric motions on the global domain.

• 한국전산구조공학회논문집
• /
• 제23권6호
• /
• pp.615-626
• /
• 2010

The discrete-layer model is proposed to analyze the edge-effect problem of laminates under extension and flexure. Based on three-dimensional elasticity theory, the displacement fields of each layer in a laminate have been treated discretely in terms of three displacement components across the thickness. The displacement fields at bottom and top surfaces within a layer are approximated by two-dimensional shape functions. Then two surfaces are connected by one-dimensional high order shape functions. Thus the p-convergent refinement on approximated one- and two-dimensional shape functions can be implemented independently of each other. The quality of present model is mostly determined by polynomial degrees of shape functions for given displacement fields. For nodal modes with physical meaning, the linear Lagrangian polynomials are considered. Additional modes without physical meaning, which are created by increasing nodeless degrees of shape functions, are derived from integrals of Legendre polynomials which have an orthogonality property. Also, it is assumed that mapping functions are linear in the light of shape of laminated plates. The results obtained by this proposed model are compared with those available in literatures. Especially, three-dimensional out-of-plane stresses in the interior and near the free edges are evaluated and convergence performance of the present model is established with the stress results.

• Structural Engineering and Mechanics
• /
• 제4권1호
• /
• pp.25-35
• /
• 1996

This paper describes a continuum variational formulation for design optimization of nonlinear structures in the elastic-plastic domain, where unloading and reloading of the structures are allowed to occur. The Total Lagrangian procedure is used for the description of the structural deformation. The direct differentiation approach is used to derive the sensitivities of the various structural response measures with respect to the design parameters. Since the material goes into the inelastic range and unloading and reloading of the structure are allowed to occur, the structural response is path dependent and an additional step is needed to integrate the constitutive equations. It can be shown, consequently, that design sensitivity analysis is also path-dependent. The theory has been discretized by the finite element technique and implemented in a structural analysis code. Mathematical programming approach is used for the optimization process. Numerical applications on trusses are performed, where cross-sectional areas and nodal point coordinates are treated as design variables. Optimal designs have been obtained and compared by using two different strategies: a two level strategy where the levels are defined accordingly the type of design variables, cross sectional areas or node coordinates, and optimizing simultaneously with respect to both types of design variables.

• Nuclear Engineering and Technology
• /
• 제25권1호
• /
• pp.20-27
• /
• 1993

이 논문은 가압경수로의 노심내 핵연료 관리용 탐색코드를 개발하기 위한 것이다. 이 목적으로 점반응도모형을 사용하여 핵연료주기 결정을 위한 FCYPRM코드를 제작하였고, 수정형 Borresen의 소격확산모형과 노달전개법에 의한 중성자 공간 해석용 CMSNAP코드를 개발하였다. 또한 수치 실험을 통하여 일련의 경험칙을 수립하고 이들을 이용하여 재장전노심 핵연료집합체 배치코드로서 ALPS코드를 개발하였다. 수치계산결과를 예시함으로서 개개 코드들의 유용성과 응용성을 입증하였으며, 이들 코드들을 가압경수로의 재장전노심 설계문제를 해결하기 위한 코드로 합성, 응용함으로서 상기 코드들이 효과적인 탐색코드가 될 수 있음을 보였다.

• 한국공간구조학회논문집
• /
• 제11권1호
• /
• pp.121-130
• /
• 2011

스페이스 프레임 구조물은 연속체 쉘 구조물의 원리를 이용하여 매우 넓은 공간을 효과적 으로 덮을 수 있는 구조물이지만 뜀좌굴 및 분기좌굴 등과 같은 불안정거동은 돔형 구조물에서는 더욱 복잡하게 나타난다. 또한 붕괴메커니즘의 이론적 연구와 실험적 연구결과들 사이에서도 많은 차이를 보인다. 본 논문에서는 미적 효과가 크며 단층의 대공간을 확보하기에 적합한 돔형 공간 구조물의 구조 불안정 특성을 접선강성방정식을 이용하여 비선형 증분해석을 수행하고, Rise-span(${\mu}$)비 및 하중모드($R_L$)에 따른 임계점과 분기점의 특성을 돔형 공간구조물의 예제를 통해 고찰하였다. 여기서 불안정점은 증분해석과정을 통해서 예측할 수 있었으며, 예제에서 낮은 ${\mu}$에서는 전체좌굴이, 높은 ${\mu}$의 경우는 절점좌굴이 지배적이며, 낮은 $R_L$에서 정점좌굴이, 높은 $R_L$에서는 전체좌굴이 지배적이고, 전체좌굴이 나타나는 경우, 분기좌굴하중은 완전형상의 극한점좌굴하중의 약 50%에서 70%의 분포를 보였다.

• 한국전산구조공학회논문집
• /
• 제17권4호
• /
• pp.375-387
• /
• 2004

팻취 보강된 철근콘크리트 구조물 해석을 위한 p-version 비선형 유한요소 모델이 제시되었다. 이방성 적층평판이론에 기초를 둔 제안된 모델은 Total Lagrangian기법에 기초한 von Karman의 대변형-소변형률 이론과 증분소성이론(incremental theory of plasticity)을 적용하였다. 콘크리트의 경화법칙(hardening rule)과 그에 따른 파괴기준을 고려하고, 단부 계면 층분리 모델(plate-end interfacial debonding model) 즉, 보강판 끝 부분에서의 콘크리트 탈락에 대한 기준으로서 Oehlers Model과 Raoof and Zhang Model을 사용하였다. 콘크리트는 두께 방향으로 층상화기법(layered model)이 이용되며, 철근과 보강판은 환산층(smeared reinforcing layer)으로 계산되도록 하였다 적분형 르장드르 다항식이 형상함수로 사용되며, 절점에서의 응력값 산출을 위해 Gauss Lobatto 수치적분법을 사용하였다. 본 연구의 목적은 p-version 유한요소법을 사용하여 RC구조물에 대한 수피해의 정확도 및 모델의 단순성을 높인 수 있도록 하였다. 따라서, 철근과 콘크리트모델에 대한 이론적 근거는 기존의 연구문헌에 근거를 두었으며, 수치해석의 적정성은 팻취 보강된 RC보와 슬래브에 대한 문헌의 실험치 및 해석치와 비교 분석되었다.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1991년도 가을 학술발표회 논문집
• /
• pp.133-138
• /
• 1991

The aim of the present paper is to develop a computer program predicting ultimate fracture strength of initially cracked structure under monotonically increasing external loads. For this purpose, two kinds of 3-D isoparametric solid elements, one 6-node wedge element and another 8-node brick element are formulated along the small deformation theory. Plasticity in the element is checked using von Mises' yield criterion. Elasto-plastic stiffness matrix of the element is calculated taking account of strain hardening effect. If the principal strain at crack tip which is one nodal point exceeds the critical strain dependin on the material property, crack tip is supposed to be opened and the crack tip node which was previously constrained in the direction perpendicular to the crack line is released. After that, the crack lay be propagated to the adjacent node. Once a crack tip node is fractured, the energy of the newly fractured node should be released which is to be absorbed by the remaining part. The accumulated reaction force which was carried by the newly fractured node so far is then applied in the opposite direction. During the action of crack tip relief force, since unloading may be occured in the plastic element, unloading check should be made. If a plastic element unloads, elastic stress-strain equation is used in the calculation of the stiffness matrix of the element, while for a loading element, elasto-plastic stress-strain equation is continuously used. Verification of the computer program is made comparing with the experimental results for center cracked panel subjected to uniform tensile load. Also some factors affecting ultimate fracture strength of initially cracked plate are investigated. It is concluded that the computer program developed here gives an accurate solution and becomes useful tool for predicting ultimate fracture load of initially cracked structural system under monotonically increasing external loads.

• 대한전기학회:학술대회논문집
• /
• 대한전기학회 1991년도 하계학술대회 논문집
• /
• pp.753-756
• /
• 1991

This paper aimed to analyse dose sensitivity to the controllable parameters of in-door radon $(^{222}Rn)$ and its decay products(Rn-D) by applying the input-output linear system theory. Physical behaviors of $^{222}Rn$ & Rn-D were analyzed in terms of $^{222}Rn$ gas generation, -migation and - infiltration to indoor environments, and the performance output-function(i.e. mean dose equivalent to Tracho-Bronchial(TB) lung region was assessed to the following ranges of the controllable parameters; a) the ventilation rate constant $({\lambda}_v)$ : $0{\sun}500[h^{-1}]$. b) the attachment rate constant$({\lambda}_a)$ : 0-500 $[h^{-1}]$. c) deposition rate constant $({\lambda}{_{d}^{u}})$: 0-50$[h^{-1}]$. A linear input-output model was reconstructed from the original models in literatures, as follows, which was modified into the matrices consisting of 111 nodal equations. a) indoor ${222}Rn$ & Rn-D Behaviour: jacobi- Porstendorfer- Bruno model. b) lung dosimerty : Jacobi-Eisfeld model. Some of the major findings, which identify the effectiveness of this model, were as follows. a) ${\lambda}_v$ is most effective, dominant controllable parameters in dose reduction, if mechanical ventilation is applied. b) ${\lambda}_v$, depending on the air particle-concentration, reduces the dose somewhat within ${\lambda}_v$<1 $h^{-1}R range. However, the dose increases conversely, ${\lambda}_v$>1 $h^{-1}R range range. c) ${\lambda}{_{d}^{4}}$ reduces the dose linearly as ${\lambda}_v$ dose. Such dose(z-axis) sentivities are shown with three-dimensional plots whoes x,y-axes are combined 2out the 3 parameter${\lambda}_v{\lambda}_s,\;{\lambda}_d^s$.

• 대한토목학회논문집
• /
• 제13권2호
• /
• pp.41-49
• /
• 1993

2차원 곡선보요소를 이용하여 대변형효과 및 탄소성거동을 고려한 유한요소법을 제시한다. Total Lagrangian 방법을 적용하며, 보요소는 평면응력요소로부터 유도된다. 또한, 임의의 단면을 갖는 평면뼈대 구조의 탄소성 해석을 수행하기 위하여 layered approach를 사용한다. 비선형 유한요소해석은 하중증분 또는 변위증분법을 적용한다. 본 연구에서 제안한 보요소의 거동과 정확도를 검증하기 위하여, 다른 문헌의 결과와 본 연구의 결과를 비교 분석한 해석예제를 제시한다.

• Structural Engineering and Mechanics
• /
• 제58권1호
• /
• pp.59-91
• /
• 2016

Laminated composite shells are commonly used in various engineering applications including aerospace and marine structures. In this paper, using semi-analytical finite strip method, the buckling behavior of laminated composite deep as well as thick shells of revolution under follower forces which remain normal to the shell is investigated. The stiffness caused by pressure is calculated for the follower forces subjected to external fibers in thick shells. The shell is divided into several closed strips with alignment of their nodal lines in the circumferential direction. The governing equations are derived based on first-order shear deformation theory which accounts for through thickness-shear flexibility. Displacements and rotations in the middle surface of shell are approximated by combining polynomial functions in the meridional direction as well as truncated Fourier series with an appropriate number of harmonic terms in the circumferential direction. The load stiffness matrix which accounts for variation of loads direction will be derived for each strip of the shell. Assembling of these matrices results in global load stiffness matrix which may be un-symmetric. Upon forming linear elastic stiffness matrix called constitutive stiffness matrix, geometric stiffness matrix and load stiffness matrix, the required elements for the second step analysis which is an eigenvalue problem are provided. In this study, different parameter effects are investigated including shell geometry, material properties, and different boundary conditions. Afterwards, the outcomes are compared with other researches. By considering the results of this article, it can be concluded that the deformation-dependent pressure assumption can entail to decrease the calculated buckling load in shells. This characteristic is studied for different examples.