• 한국항공우주학회지
• /
• 제36권11호
• /
• pp.1049-1055
• /
• 2008

최근 다분야 전산유체 역학에서는 설계 최적화, 공탄성, 강제 경계 운동 등에서 요구되어지는 이동경계문제를 다루게 된다. 이동경계의 변위가 클 경우 강건하고 효율적인 격자 변형 알고리즘의 개발이 필요하다. 본 연구에서는 유한 대형요소와 초월유한보간에 근거한 체적격자 변형 코드를 개발하였고, 정렬격자 다중 블록 Navier-Stokes 코드와 연계하였다. 개발된 코드의 검증을 위해 주기적으로 진동 운동을 하는 에어포일 문제에 대해 계산을 수행 하였고 양력, 항력, 모멘트 계수의 이력 계산 결과가 실험 결과와 잘 일치하였다.

• International Journal of CAD/CAM
• /
• 제8권1호
• /
• pp.1-10
• /
• 2009

This paper presents a topology optimization approach using element-free Galerkin method (EFGM) for the optimal design of compliant mechanisms with geometrically non-linearity. Meshless method has an advantage over the finite element method(FEM) because it is more capable of handling large deformation resulted from geometrical nonlinearity. Therefore, in this paper, EFGM is employed to discretize the governing equations and the bulk density field. The sensitivity analysis of the optimization problem is performed by incorporating the adjoint approach with the meshless method. The Lagrange multipliers method adjusted for imposition of both the concentrated and continuous essential boundary conditions in the EFGM is proposed in details. The optimization mathematical formulation is developed to convert the multi-criteria problem to an equivalent single-objective problem. The popularly applied interpolation scheme, solid isotropic material with penalization (SIMP), is used to indicate the dependence of material property upon on pseudo densities discretized to the integration points. A well studied numerical example has been applied to demonstrate the proposed approach works very well and the non-linear EFGM can obtain the better topologies than the linear EFGM to design large-displacement compliant mechanisms.

• 한국정밀공학회지
• /
• 제20권2호
• /
• pp.192-198
• /
• 2003

A finite element analysis of end-to-end artery/PTFE anastomosis has been presented in this study to evaluate the distribution of compliance and stresses in the vicinity of the anastomosis due to any mismatch in compliance characteristics. The artery wall was assumed to be made of linear isotropic material in this simplified model and a nonlinear analysis and convergency study with respect to increasing meshed element numbers were performed with a mean artery pressure loading of the artery-PTFE model. Also, sub-modeling method was Introduced to progress the accuracy of the finite element analysis. The results are as follow : 1. A hypercompliant zone on the artery side was observed around 4.Omm from the anastomosis and a high hoop stresses in the wall of artery and PTFE was dominant. 2. An artery displays large deformation so that nonlinear analysis and sub-modeling method was used. 3. An anastomosis with the thinner thickness and larger diameter PTFE (B type) could reduce the compliance disagreement.

• 한국CDE학회논문집
• /
• 제14권3호
• /
• pp.197-206
• /
• 2009

In a finite element analysis of the metal forming processes having large plastic deformation, largely distorted elements are unstable and hence they influence upon the result toward negative way so that adaptive remeshing is required to avoid a failure in the numerical computation. Therefore automatic mesh generation and regeneration is very important to avoid a numerical failure in a finite element analysis. In case of generating quadrilateral mesh, the automation is more difficult than that of triangular mesh because of its geometric complexity. However its demand is very high due to the precision of analysis. Thus, in this study, an automatic quadrilateral mesh generation and regeneration method using grid-based approach is developed. The developed method contains decision of grid size to generate initial mesh inside a two dimensional domain, classification of boundary angles and inner boundary nodes to improve element qualities in case of concave domains, and boundary projection to construct the final mesh.

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

The finite element analysis of plate and shell structures has been one of the major research interests for many years because of the technological importance of such structures. Quite often these structures are constructed by laminated composites. This is due to the high specific stiffness and strength of composite structures. The main objective of this paper is to extend the use of an improved degenerated shell element to the large displacement analysis of plates and shells with laminated composites. The total Lagrangian approach has been chosen for the definition of the deformation and the solution to the nonlinear equilibrium equations is obtained by the Newton-Raphson method.

• 대한기계학회논문집A
• /
• 제29권4호
• /
• pp.534-539
• /
• 2005

Nodal displacements are referred to the initial configuration in the total Lagrangian formulation and to the last converged configuration in the updated Lagrangian furmulation. This research proposes a relative nodal displacement method to represent the position and orientation for a node in truss structures. Since the proposed method measures the relative nodal displacements relative to its adjacent nodal reference frame, they are still small for a truss structure undergoing large deformations for the small size elements. As a consequence, element formulations developed under the small deformation assumption are still valid for structures undergoing large deformations, which significantly simplifies the equations of equilibrium. A structural system is represented by a graph to systematically develop the governing equations of equilibrium for general systems. A node and an element are represented by a node and an edge in graph representation, respectively. Closed loops are opened to form a spanning tree by cutting edges. Two computational sequences are defined in the graph representation. One is the forward path sequence that is used to recover the Cartesian nodal displacements from relative nodal displacement sand traverses a graph from the base node towards the terminal nodes. The other is the backward path sequence that is used to recover the nodal forces in the relative coordinate system from the known nodal forces in the absolute coordinate system and traverses from the terminal nodes towards the base node. One open loop and one closed loop structure undergoing large deformations are analyzed to demonstrate the efficiency and validity of the proposed method.

• Structural Engineering and Mechanics
• /
• 제35권6호
• /
• pp.677-697
• /
• 2010

This paper focuses on geometrically non-linear static analysis of a simply supported beam made of hyperelastic material subjected to a non-follower transversal uniformly distributed load. As it is known, the line of action of follower forces is affected by the deformation of the elastic system on which they act and therefore such forces are non-conservative. The material of the beam is assumed as isotropic and hyperelastic. Two types of simply supported beams are considered which have the following boundary conditions: 1) There is a pin at left end and a roller at right end of the beam (pinned-rolled beam). 2) Both ends of the beam are supported by pins (pinned-pinned beam). In this study, finite element model of the beam is constructed by using total Lagrangian finite element model of two dimensional continuum for a twelve-node quadratic element. The considered highly non-linear problem is solved by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. In order to use the solution procedures of Newton-Raphson type, there is need to linearized equilibrium equations, which can be achieved through the linearization of the principle of virtual work in its continuum form. In the study, the effect of the large deflections and rotations on the displacements and the normal stress and the shear stress distributions through the thickness of the beam is investigated in detail. It is known that in the failure analysis, the most important quantities are the principal normal stresses and the maximum shear stress. Therefore these stresses are investigated in detail. The convergence studies are performed for various numbers of finite elements. The effects of the geometric non-linearity and pinned-pinned and pinned-rolled support conditions on the displacements and on the stresses are investigated. By using a twelve-node quadratic element, the free boundary conditions are satisfied and very good stress diagrams are obtained. Also, some of the results of the total Lagrangian finite element model of two dimensional continuum for a twelve-node quadratic element are compared with the results of SAP2000 packet program. Numerical results show that geometrical nonlinearity plays very important role in the static responses of the beam.

• 한국정밀공학회:학술대회논문집
• /
• 한국정밀공학회 2003년도 춘계학술대회 논문집
• /
• pp.1569-1572
• /
• 2003

Meshfree approximations exhibit significant potential to solve partial differential equations. Meshfree methods have been successfully applied to various problems which the traditional finite element methods have difficulties to handle, including the quasi-static and dynamic fracture. large deformation problems, contact problems, and strain localization problems. A meshfree method based on the reproducing kernel particle approximation(RKPM) is applied to sheet metal forming analysis in this research. Metal forming examples, such as stretch forming and flanging operation, are analyzed to demonstrate the performance of the proposed meshfree method for largely deformed elasto-plastic material.

• Geomechanics and Engineering
• /
• 제11권1호
• /
• pp.1-39
• /
• 2016

The finite element method (FEM), discrete element method (DEM), and Discontinuous deformation analysis (DDA) are among the standard numerical techniques applied in computational geo-mechanics. However, in some cases there no possibility for modelling by traditional finite analytical techniques or other mesh-based techniques. The solution presented in the current study as a completely Lagrangian and mesh-free technique is smoothed particle hydrodynamics (SPH). This method was basically applied for simulation of fluid flow by dividing the fluid into several particles. However, several researchers attempted to simulate soil-water interaction, landslides, and failure of soil by SPH method. In fact, this method is able to deal with behavior and interaction of different states of materials (liquid and solid) and multiphase soil models and their large deformations. Soil indicates different behaviors when interacting with water, structure, instrumentations, or different layers. Thus, study into these interactions using the mesh based grids has been facilitated by mesh-less SPH technique in this work. It has been revealed that the fast development, computational sophistication, and emerge of mesh-less particle modeling techniques offer solutions for problems which are not modeled by the traditional mesh-based techniques. Also it has been found that the smoothed particle hydrodynamic provides advanced techniques for simulation of soil materials as compared to the current traditional numerical methods. Besides, findings indicate that the advantages of applying this method are its high power, simplicity of concept, relative simplicity in combination of modern physics, and particularly its potential in study of large deformations and failures.

• 소성∙가공
• /
• 제4권2호
• /
• pp.131-140
• /
• 1995

In this paper damage propagation of a material during forming is investigated with the concept of continuum damage mechanics. An isotropic damage model based on the theory of materials of type N is adopted to describe the damage process of a ductile material with large elasto-viscoplastic deformation. The stiffness degradation of the loaded material is chosen as a damage measure. The highly nonlinear equilibrium equations are reduced to the incremental weak form and approximated by the total Lagrangian finite element method. To simulate contact condition, extended interior penalty method with modified coulomb friction law is adopted. The displacement control method along with the modified Riks' continuation technique is used to solve the incremental iterative equations. As numerical examples, upsetting problem and backward extrusion problem are simulated and the results of damage propagation and $J_2$ stress contours with and without friction are presented.