• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.11
• /
• pp.1844-1851
• /
• 1997

Displacement fields and interface stresses are obtained by modifying the potential energy functional with a penalty function which enforces the continuity of stresses at the interface of two-materials. Based on the displacement field and the interface stresses, a new methodology to generate a continuous stress field over the entire domain including the interface of the dissimilar materials has been proposed by combining the L$^{2}$ projection method of stress-smoothing and the Loubignac's iterative method of improving the displacement field. Stress analysis was carried out on two examples which are made of highly dissimilar materials. As a result of the analysis, it is found that the proposed method provides improved continuity of the stress field over the entire domain as well as predicting accurate nodal stresses at the interface. In contrast, the conventional displacement-based finite element method provides significant stress discontinuties at the interfaces. In addition, it was found that the total strain energy evaluated from the improved continuous stress field converge to the exact value as increasing the number of iterations in the proposed method.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.10
• /
• pp.1598-1608
• /
• 1997

Based on the application of te theory of conjugate approximations and the Loubignac's iterative method in a local region, a method to improve the stress filed in a displacement-formulated finite element solution has been proposed. The validity of the proposed method has been tested through two examples : a thick cylinder under internal pressure loading and an infinite plate with a central circular hole subjected to uniaxial tension. As a result of analysis of the examples, it was found that the stress field obtained for the local region model by the proposed method approximates well for the whole domain model. In addition, it was found that because of a significant decrease in the computing time to obtain the improved stress field, the proposed method is efficient and useful for the detailed stress analysis in local regions.

• Structural Engineering and Mechanics
• /
• v.22 no.2
• /
• pp.151-167
• /
• 2006

A p-version of the finite element method in conjunction with the modeling dynamic method using the arc-length stretch deformation is considered to determine the bending natural frequencies of a cantilever flexible plate mounted on the periphery of a rotating hub. The plate Fourier p-element is used to set up the linear equations of motion. The transverse displacements are formulated in terms of cubic polynomials functions used generally in FEM plus a variable number of trigonometric shapes functions representing the internals DOF for the plate element. Trigonometric enriched stiffness, mass and centrifugal stiffness matrices are derived using symbolic computation. The convergence properties of the rotating plate Fourier p-element proposed and the results are in good agreement with the work of other investigators. From the results of the computation, the influences of rotating speed, aspect ratio, Poisson's ratio and the hub radius on the natural frequencies are investigated.

• Journal of Welding and Joining
• /
• v.18 no.4
• /
• pp.104-110
• /
• 2000

An integral equation representation of cracks was presented, which differs from well-known "dislocation layer" representation. In this new representation, an integral equation representation of cracks was developed and coupled to the direct boundary-element method for treatment of cracks in plane finite bodies. The method was developed for in-plane (modes I and II) loadings only. In this paper, the method is formulated and applied to mode III problems involving smooth or kinked cracks in finite region. The results are compared to exact solutions where available and the method is shown to be very accurate despite of its simplicity.implicity.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1995.10a
• /
• pp.145-152
• /
• 1995

A general stiffener element which includes transverse shear deformation(TSD) is formulated using the p-version of finite element method. Hierarchic C"-shape functions, derived from Integrals of Legendre polynomials, are used to define the assembled stiffness matrix of the stiffener and plate on the basis of 5 D.0.F displacement fields. The stiffness matrix for the stiffener with respect to the local reference frame is transformed to the plate reference system by applying the appropriate transformation matrices in order to insure compatibility of displacements at the junction of the stiffener and plate. The transformation matrices which account for the orientation and the eccentricity effects of the stiffener with respect to the plate reference axes are used to find local behavior at the junction of the stiffener and the relative contributions of the plate and stiffener to the strength of the composite system. The results obtained by the p-version of the finite element method are compared with the results in literatures, especially those by the h-version software, MICROFEAP-II.P-II.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.1
• /
• pp.163-171
• /
• 1996

For the effective analysis of two dimensional plane problems, Treffiz finite elements and cavity elements have been proposed. These element matrix equaitons were formulated on the basis of hybrid variational principle and Treffiz function sets derived consitstently from the complex theoy of plane elasticity. In order to suggest the accuracy chatacteristics of the proposed Treffiz elements typical plane problems were analyzed and these results were compared with ones obtained by using the conveintional displacement type elements. The accuracy of the proposed elements is less sensitive to the element size and shape than the conventional displacement type elements. These elements, being able to be formed with multi-nodes, give the convenient modeling of an analytic domain. The cavity elements give the comparatively exact values of stress concentration factors of stress intensity factors and can be effectively used for the analysis of mechanical stuctures containing various cavities.

• Coupled systems mechanics
• /
• v.2 no.4
• /
• pp.389-410
• /
• 2013

Three-dimensional Lagrangian fluid finite element is applied to seismic response analysis of an oil storage tank with a floating roof. The fluid element utilized in the present analysis is formulated based on the displacement finite element method considering only volumetric elasticity and its element stiffness matrix is derived by using one-point integration method in order to avoid volumetric locking. The method usually adds a rotational penalty stiffness to satisfy the irrotational condition for fluid motion and modifies element mass matrices through the projected mass method to suppress spurious hourglass-mode appeared in compensation for one-point integration. In the fluid element utilized in the present paper, a small hourglass stiffness is employed. The fluid and structure domains for the objective oil storage tank are modeled by eight-node solid elements and four-node shell elements, respectively, and the transient response of the floating roof structure or the free surface are evaluated by implicit direct time integration method. The results of seismic response analyses are compared with those by other method and the validation of the present analysis using three-dimensional Lagrangian fluid finite elements is shown.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.17 no.4
• /
• pp.63-71
• /
• 2009

The elastic effects of the valve train components are analyzed by using the finite element models of the rocker arm and valve. The whole equations of motion of the valvetrain of an internal combustion engine formulated by finite element techniques are solved by imposing the contact conditions with the augmented Lagrange multiplier method. The velocity and acceleration constraints as well as the displacement constraints are imposed on the contact points. The numerical simulations show that, even if the magnitude of the elastic deformation of the components is very small, it may have large effects on the valvetrain dynamics of a high-speed engine.

• Computers and Concrete
• /
• v.15 no.3
• /
• pp.373-389
• /
• 2015

The development of a finite element model for the geometric and material nonlinear analysis of bonded prestressed concrete continuous beams is presented. The nonlinear geometric effect is introduced by the coupling of axial and flexural fields. A layered approach is applied so as to consider different material properties across the depth of a cross section. The proposed method of analysis is formulated based on the Euler-Bernoulli beam theory. According to the total Lagrangian description, the constructed stiffness matrix consists of three components, namely, the material stiffness matrix reflecting the nonlinear material effect, the geometric stiffness matrix reflecting the nonlinear geometric effect and the large displacement stiffness matrix reflecting the large displacement effect. The analysis is capable of predicting the nonlinear behaviour of bonded prestressed concrete continuous beams over the entire loading stage up to failure. Some numerical examples are presented to demonstrate the validity and applicability of the proposed model.

• Journal of Power System Engineering
• /
• v.9 no.1
• /
• pp.42-49
• /
• 2005

The vibration analysis of a rectangular plate with stiffeners is formulated by using the transfer stiffness coefficient method (TSCM). This method is based on the concept of the successive transmission of stiffness coefficients which are defined as the relationship between the force vector and the displacement vector at an arbitrary nodal line. In order to confirm the validity of the present method, bending vibration analysis for a rectangular plate with stiffener is carried out on a personal computer by using the present method and the finite element method (FEM). Through comparing computational results of the TSCM and the FEM, the effectivness of the TSCM from the viewpoint of computational cost, that is, computational time and storage is demonstrated.