• 한국전산구조공학회논문집
• /
• 제30권4호
• /
• pp.335-341
• /
• 2017

본 논문에서는 고차미분 연속성을 가지는 형상함수에 기초하여 오일러-베르누이 보 유한요소모델을 정식화하였으며, 다양한 경계조건들에 대하여 그 성능을 평가하였다. 이러한 유한요소 모델들은 새로이 개발되는 고차 보 이론들과 논로컬 탄성이론에 기초한 보 이론들의 유한요소해석에 필요하다. 그러나 고차 연속성을 가지는 유한요소에 대한 성능평가는 문헌에서 찾아보기 어렵다. 따라서 본 연구에서는 $C^2$ 및 $C^3$ 두 종류의 고차 유한요소들을 정식화하여 외팔보, 단순지지, 고정-힌지 등의 경계조건들을 적용하고 정적해석을 수행하였다. 고전적인 경계조건들 이외에도 고차 경계조건들이 보의 거동에 미치는 영향을 비교분석하였다. 경계조건에 따라서는 처짐의 미분 값들이 경계주변에서 진동하는 현상이 관찰되었으며, 이는 기하학적 경계조건들에 대하여 뚜렷이 나타난다. 특히 고정단과 같은 경계에서의 변위의 고차미분 조건은 이러한 불안정한 현상을 유발한다. 본 연구에서 얻어진 결과들은 고차 미분 연속성을 가지는 유한요소 이용에 가이드라인으로서 역할을 할 수 있을 것으로 기대된다.

• Steel and Composite Structures
• /
• 제23권1호
• /
• pp.41-51
• /
• 2017

A $C^0$ FE model developed based on an efficient higher order zigzag theory is used for hygrothermal analysis of laminated composite plates. The $C^0$ FE model satisfies the inter-laminar shear stress continuity at the interfaces and zero transverse shear stress conditions at plate top and bottom. In this model the first derivatives of transverse displacement have been treated as independent variables to circumvent the problem of $C^1$ continuity associated with the above plate theory. In the present theory the above mentioned $C^0$ continuity of the present element is compensated in the stiffness matrix formulation by using penalty parameter approach. In order to avoid stress oscillations observed in the displacement based finite element, the stress field derived from temperature/moisture fields (initial strains) must be consistent with total strain field. Special steps are introduced by field consistent approach (e.g., sampling at gauss points) to compensate this problem. A nine noded $C^0$ continuous isoparametric element is used in the proposed FE model. Comparison of present numerical results with other existing solutions shows that the proposed FE model is efficient, accurate and free of locking.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2006년도 정기 학술대회 논문집
• /
• pp.849-856
• /
• 2006

This paper is to examine the applicability of ordinary Kriging interpolation(OK) to the p-adaptivity of the finite element analysis that is based on variogram. In the p-refinement, the analytical domain has to be refined automatically to obtain an acceptable level of accuracy by increasing the p-level non-uniformly or selectively. In case of non-uniform p-distribution, the continuity between elements with different polynomial orders is achieved by assigning zero higher-order derivatives associated with the edge in common with the lower-order derivatives. It is demonstrated that the validity of the proposed approach by analyzing results for stress singularity problem.

• Steel and Composite Structures
• /
• 제38권3호
• /
• pp.281-291
• /
• 2021

An equivalent single-layer theory (EST) is put forward for analyzing free vibrations of steel-concrete composite beams (SCCB) based on a higher-order beam theory. In the EST, the effect of partial interaction between sub-beams and the transverse shear deformation are taken into account. After using the interlaminar shear force continuity condition and the shear stress free conditions at the top and bottom surface, the displacement function of the EST does not contain the first derivatives of transverse displacement. Therefore, the C0 interpolation functions are just demanded during its finite element implementation. Finally, the EST is validated by comparing the results of two simply-supported steel-concrete composite beams which are tested in laboratory and calculated by ANSYS software. Then, the influencing factors for free vibrations of SCCB are analyzed, such as, different boundary conditions, depth to span ratio, high-order shear terms, and interfacial shear connector stiffness.

• Structural Engineering and Mechanics
• /
• 제85권2호
• /
• pp.147-161
• /
• 2023

Based on the ideas of variational differential quadrature (VDQ) and finite element method (FEM), a numerical approach named as VDQFEM is applied herein to study the large deformations of plate-type structures under static loading with arbitrary shape hole made of functionally graded graphene platelet-reinforced composite (FG-GPLRC) in the context of higher-order shear deformation theory (HSDT). The material properties of composite are approximated based upon the modified Halpin-Tsai model and rule of mixture. Furthermore, various FG distribution patterns are considered along the thickness direction of plate for GPLs. Using novel vector/matrix relations, the governing equations are derived through a variational approach. The matricized formulation can be efficiently employed in the coding process of numerical methods. In VDQFEM, the space domain of structure is first transformed into a number of finite elements. Then, the VDQ discretization technique is implemented within each element. As the last step, the assemblage procedure is performed to derive the set of governing equations which is solved via the pseudo arc-length continuation algorithm. Also, since HSDT is used herein, the mixed formulation approach is proposed to accommodate the continuity of first-order derivatives on the common boundaries of elements. Rectangular and circular plates under various boundary conditions with circular/rectangular/elliptical cutout are selected to generate the numerical results. In the numerical examples, the effects of geometrical properties and reinforcement with GPL on the nonlinear maximum deflection-transverse load amplitude curve are studied.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2008년도 정기 학술대회
• /
• pp.216-221
• /
• 2008

Shape design optimization for linear elasticity problem is performed using isogeometric analysis method. In many design optimization problems for real engineering models, initial raw data usually comes from CAD modeler. Then designer should convert this CAD data into finite element mesh data because conventional design optimization tools are generally based on finite element analysis. During this conversion there is some numerical error due to a geometry approximation, which causes accuracy problems in not only response analysis but also design sensitivity analysis. As a remedy of this phenomenon, the isogeometric analysis method is one of the promising approaches of shape design optimization. The main idea of isogeometric analysis is that the basis functions used in analysis is exactly same as ones which represent the geometry, and this geometrically exact model can be used shape sensitivity analysis and design optimization as well. In shape design sensitivity point of view, precise shape sensitivity is very essential for gradient-based optimization. In conventional finite element based optimization, higher order information such as normal vector and curvature term is inaccurate or even missing due to the use of linear interpolation functions. On the other hands, B-spline basis functions have sufficient continuity and their derivatives are smooth enough. Therefore normal vector and curvature terms can be exactly evaluated, which eventually yields precise optimal shapes. In this article, isogeometric analysis method is utilized for the shape design optimization. By virtue of B-spline basis function, an exact geometry can be handled without finite element meshes. Moreover, initial CAD data are used throughout the optimization process, including response analysis, shape sensitivity analysis, design parameterization and shape optimization, without subsequent communication with CAD description.