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

In this paper, three dimensional linear conforming variable-finite elements are presented with the aid of a smoothed integration (a class of stabilized conforming nodal integration), for mnltiscale mechanics problems. These elements meet the desirable properties of an interpolation such as the Kronecker delta condition, the partition of unity condition and the positiveness of interpolation function. The necessary condition of linear exactness is fully relaxed by employing the smoothed integration, which renders us to meet the linear exactness in a straightforward manner. This novel element description extend the category of the conventional finite elements space to ration type function space and give the flexibility on the number of nodes of element which are fixed in the conventional finite elements. Several examples are provided to show the convergence and the accuracy of the proposed elements, and to demonstrate their potential with emphasis on the multiscale mechanics problems such as global/local analysis, nonmatching contact problems, and modeling of composite material with defects.

• Structural Engineering and Mechanics
• /
• 제10권5호
• /
• pp.491-504
• /
• 2000

This paper is a presentation of a physical model for the elastic-partly plastic behavior of rectangular hollow section pinned struts subjected to static cyclic axial loading and the evaluation of the compressive strength of retrofitted damaged struts. Retrofitting is achieved by welding stiffening plates along the webs of damaged struts. The shape of the elastic and permanent deformations of the strut axis satisfy the conditions at the ends and midspan. Continuous functions of the geometric variables of stress distributions in the yielded zone are evaluated by interpolation between three points along each partly plastic zone. Permanent deformations of the partly plastic region are computed and used to update the shape of the unloaded strut. The necessity of considering geometric nonlinearity is discussed. The sensitivity of the results to the location of interpolation points, the shape of the permanent deformation and material hysteretic properties is investigated.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2010년 춘계학술대회논문집
• /
• pp.551-555
• /
• 2010

Electrowetting is a versatile tool to handle tiny droplets and forms a backbone of digital microfluidics. Numerical analysis is necessary to fully understand the dynamics of electrowetting, especially in designing electrowetting-based devices, such as liquid lenses and reflective displays. We developed a numerical method to analyze the general contact-line problems, incorporating dynamic contact angle models. The method is based on the conservative level set method to capture the interface of two fluids without loss of mass. We applied the method to the analysis of spreading process of a sessile droplet for step input voltages and oscillation of the droplet for alternating input voltages in electrowetting. The result was compared with experimental data. It is shown that contact line friction significantly affects the contact line motion and the oscillation amplitude. The pinning process of contact line was well represented by including the hysteresis effect in the contact angle models. In level set method, in the mean time, material properties are made to change smoothly across an interface of two materials with different properties by introducing an interpolation or smoothing scheme. So far, the weighted arithmetic mean (WAM) method has been exclusively adopted in level set method, without complete assessment for its validity. We viscosity, thermal conductivity, electrical conductivity, and permittivity, can be an alternative. I.e., the WHM gives more accurate results than the WAM method in certain circumstances. The interpolation scheme should be selected considering various characteristics including type of property, ratio of property of two fluids, geometry of interface, and so on.

• Structural Engineering and Mechanics
• /
• 제12권5호
• /
• pp.507-526
• /
• 2001

In standard finite element algorithms, the local stability conditions are not accounted for in the formulation of the tangent stiffness matrix. As a result, the loss of the local stability is not adequately related to the onset of the global instability. The phenomenon typically arises with material-type localizations, such as shear bands and plastic hinges. This paper addresses the problem in the context of the planar, finite-strain, rate-independent, materially non-linear beam theory, although the proposed technology is in principle not limited to beam structures. A weak formulation of Reissner's finite-strain beam theory is first presented, where the pseudocurvature of the deformed axis is the only unknown function. We further derive the local stability conditions for the large deformation case, and suggest various possible combinations of the interpolation and numerical integration schemes that trigger the simultaneous loss of the local and global instabilities of a statically determined beam. For practical applications, we advice on a procedure that uses a special numerical integration rule, where interpolation nodes and integration points are equal in number, but not in locations, except for the point of the local instability, where the interpolation node and the integration point coalesce. Provided that the point of instability is an end-point of the beam-a condition often met in engineering practice-the procedure simplifies substantially; one of such algorithms uses the combination of the Lagrangian interpolation and Lobatto's integration. The present paper uses the Galerkin finite element discretization, but a conceptually similar technology could be extended to other discretization methods.

• Steel and Composite Structures
• /
• 제48권5호
• /
• pp.583-597
• /
• 2023

This paper introduces a novel approach to multi-material topology optimization (MTO) targeting in-plane bi-directional functionally graded (IBFG) non-uniform thickness Reissner-Mindlin plates, employing an alternative active phase approach. The mathematical formulation integrates a first shear deformation theory (FSDT) to address compliance minimization as the objective function. Through an alternating active-phase algorithm in conjunction with the block Gauss-Seidel method, the study transforms a multi-phase topology optimization challenge with multi-volume fraction constraints into multiple binary phase sub-problems, each with a single volume fraction constraint. The investigation focuses on IBFG materials that incorporate adequate local bulk and shear moduli to enhance the precision of material interactions. Furthermore, the well-established mixed interpolation of tensorial components 4-node elements (MITC4) is harnessed to tackle shear-locking issues inherent in thin plate models. The study meticulously presents detailed mathematical formulations for IBFG plates in the MTO framework, underscored by numerous numerical examples demonstrating the method's efficiency and reliability.

• Steel and Composite Structures
• /
• 제46권1호
• /
• pp.33-51
• /
• 2023

This research presents a multi-material topology optimization for functionally graded material (FGM) and nonFGM with elastic buckling criteria. The elastic buckling based multi-material topology optimization of functionally graded steels (FGSs) uses a Jacobi scheme and a Method of Moving Asymptotes (MMA) as an expansion to revise the design variables shown first. Moreover, mathematical expressions for modified interpolation materials in the buckling framework are also described in detail. A Solid Isotropic Material with Penalization (SIMP) as well as a modified penalizing material model is utilized. Based on this investigation on the buckling constraint with homogenization material properties, this method for determining optimal shape is presented under buckling constraint parameters with non-homogenization material properties. For optimal problems, minimizing structural compliance like as an objective function is related to a given material volume and a buckling load factor. In this study, conflicts between structural stiffness and stability which cause an unfavorable effect on the performance of existing optimization procedures are reduced. A few structural design features illustrate the effectiveness and adjustability of an approach and provide some ideas for further expansions.

• 한국전산구조공학회논문집
• /
• 제14권4호
• /
• pp.495-503
• /
• 2001

기능 경사 복합재에서는 열적 탄성 거동이 성분의 구성분포에 의해 명백하게 구분되어지도록 두가지의 성분입자들을 특정한 체적분율 분포에 따라서 혼합한다. 따라서, 설계자는 주어진 제약조건에 대해 목적하는 성능에 적합한 기능 경사 복합재를 설계하기 위해서 최적의 체적분율 분포를 결정해야만 한다. 본 연구에서는 금속과 세라믹으로 구성된 내열 기능 경사 복합재의 2차원 체적분율을 최적화하기 위하여 내부벌칙함수법과 유한차분법을 사용한 수치 최적화기법을 제안하였다. 최적화 효율을 위해 단일 설계변수의 유한개의 균질 사각형 격자로 기능 경사 복합재의 영역을 나누었다. 그렇지만, 연속적인 체적분율을 구현하기 위하여 최적설계 후에 전체적으로 연속적인 이차원 선형함수로써 불연속적인 체적분율을 보간하였다.

• International Journal of Fuzzy Logic and Intelligent Systems
• /
• 제8권4호
• /
• pp.264-269
• /
• 2008

In this paper, two neural networks as a material model, which are based on the state-space method, have been proposed. One outputs the rates of inelastic strain and material internal variables whereas the outputs of the other are the next state of the inelastic strain and material internal variables. Both the neural networks were trained using input-output data generated from Chaboche's model and successfully converged. The former neural network could reproduce the original stress-strain curve. The neural network also demonstrated its ability of interpolation by generating untrained curve. It was also found that the neural network can extrapolate in close proximity to the training data.

• Structural Engineering and Mechanics
• /
• 제45권4호
• /
• pp.569-594
• /
• 2013

The free vibration of axially functionally graded tapered arches including shear deformation and rotatory inertia are studied through solving the governing differential equation of motion. Numerical results are presented for circular, parabolic, catenary, elliptic and sinusoidal arches with hinged-hinged, hinged-clamped and clamped-clamped end restraints. In this study Differential Quadrature element of lowest order (DQEL) or Lagrangian Interpolation technique is applied to solve the problems. Three general taper types for rectangular section are considered. The lowest four natural frequencies are calculated and compared with the published results.

• 한국콘크리트학회:학술대회논문집
• /
• 한국콘크리트학회 1997년도 봄 학술발표회 논문집
• /
• pp.398-405
• /
• 1997

The finite element program is presented for the analysis of the moment distribution in continuous slabs of simply supported girders. The program includes the material nonlinearity of the continuous and steel reinforcements of the RC slabs, but assumes that the PC beam and cross beam behave linearly. In modeling the PC slabs and girders, the four node degenerated shell element formulated based on the assumed strain interpolation and the 3D beam element are used, respectively. Using the program, the influence of the existence of the cross bean, the filling of open joints on the continuous at supports, and perfect continuation of precast girder elements are examined.