• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.3
• /
• pp.615-622
• /
• 1995

It may be inefficient to conduct the aeroelastic analysis by using full-scale conventional finite-element analyses or experiments, from the initial design phase, for an aircraft wing which can be considered as the discontinuum complex structure with composite laminated skins. In this paper, therefore more efficient aeroelastic analysis has been conducted for a box-beam typed aircraft wing by using the equivalent continuum beam-rod model which is derived from the concept of energy equivalence. Equivalent structural properties of the continuum beam-rod model are obtained from the direct comparison of the finite-element matrices of continuum beam-rod model with those of box-beam typed aircraft wing. Numerical results by the continuum beam-rod model approach are compared with those by the conventional finite-element analysis approach to show that the continuum beam-rod model proposed herein is quite satisfactory as a simplified model of aircraft wing structure for aeroelastic analyses.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.11
• /
• pp.2704-2710
• /
• 1993

A simple and straightforward method is introduced for developing continuum beam-rod model of a box-beam typed aircraft wing with composite layered skin based on "energy equivalence." The equivalent continuum structral properties are obtained from the direct comparison of the reduced stiffness and mass matrices for box-beam typed wing with those for continuum beam-rod model. The stiffness and mass matrices are all represented in terms of the continuum degrees-of freedom defined in this paper. The finite-element method. The advantage of the present continuum method is to give every continuum structural properties including all possible coupling terms which represent the couplings between different deformations. To evaluate the continuum method developed in this paper, free vibration analyses for both continuum beam-rod and box-beam are conducted. Numerical tests show that the present continuum method gives very reliable structural and dynamic properties compared to the results by the conventional finite-element analysis. analysis.

• Structural Engineering and Mechanics
• /
• v.65 no.1
• /
• pp.33-41
• /
• 2018

In this paper, using modified couple stress theory in place of classical continuum theory, and using shell model in place of beam model, vibrational behavior of nanotubes is investigated via the finite element method. Accordingly classical continuum theory is unable to correctly compute stiffness and account for size effects in micro/nanostructures, higher order continuum theories such as modified couple stress theory have taken on great appeal. In the present work the mass-stiffness matrix for cylindrical shell element is developed, and by means of size-dependent finite element formulation is extended to more precisely account for nanotube vibration. In addition to modified couple stress cylindrical shell element, the classical cylindrical shell element can also be defined by setting length scale parameter to zero in the equations. The boundary condition were assumed simply supported at both ends and it is shown that the natural frequency of nano-scale shell using the modified coupled stress theory is larger than that using the classical shell theory and the results of Ansys. The results have indicated using the modified couple stress cylindrical shell element, the rigidity of the nano-shell is greater than that in the classical continuum theory, which results in increase in natural frequencies. Besides, in addition to reducing the number of elements required, the use of this type of element also increases convergence speed and accuracy.

• Journal of the korean Society of Automotive Engineers
• /
• v.13 no.3
• /
• pp.58-67
• /
• 1991

A theory of continuum damage mechanics based on the theory of materials of type N was developed and its nonlinear finite element approximation and numerical simulation was carried out. To solve the finite elastoplasticity problems, reasonable kinematics of large deformed solids was introduced and constitutive relations based on the theory of materials of type-N were derived. These highly nonlinear equations were reduced to the incremental weak formulation and approximated by the theory of nonlinear finite element method. Two types of problems, compression moulding problem and pure bending problem, were solved for aluminum 2024.

• Structural Engineering and Mechanics
• /
• v.35 no.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.

• Journal of Theoretical and Applied Mechanics
• /
• v.1 no.1
• /
• pp.97-109
• /
• 1995

In this paper, a finite viscoelastic continuum model for rubber and its finite element analysis are presented. This finite viscoelatic model based on continuum mechanics is an extended model of Johnson and Wuigley's 1-D model. In this extended model, continuum based kinematic measures are rigorously defied and by using this kinematic measures, elastic stage law and flow rule are introduced. In kinematics, three configuration are introduced. In kinematics, three configuration are introduced. They are reference, current and virtual visco configurations. In elastic state law, it is assumed that at a certain time, there exists an elastic potential which describes the recoverable elastic energy. From this elastic potential, elastic state law is derived. The proposed flow rule is based on phenomenological observation. The flow rule gives precise relaxation response. In finite element approximation, mixed Lagrangian description is used, where total and similar method of updated Lagrangian descriptions are used together. This approach reduces the numerical job and gives simple nonlinear syatem of equations. To satisfy the incompressible condition, penalty-type modified Mooney-Rivlin energy function is adopted. By this method nearly incompressible condition is obtain the virtual visco configuration. For verification, uniaxial stretch tests are simulated for various stretch rates. The simulated results show good agreement with experiments. As a practical experiments. As a preactical example, pressurized rubber plate is simulated. The result shows finite viscoelastic effects clearly.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1997.10a
• /
• pp.24-27
• /
• 1997

In the present work, rigid-plastic continuum elements employing the shape change and anisotropic effects are derived for the purpose of applying more realistic blankholding force condition in three-dimensional finite element analysis of sheet metal forming process. In order to incorporate the effect of shape change effectively in the derivation of finite element equation using continuum element for sheet metal forming, the convected coordinate system is introduced, rendering the analysis more rigorous and accurate. The formulation is extended to cover the orthotropic material using Hill's quadratic yield function. For the purpose of applying more realistic blankholding force condition, distributed normal and associated frictional tangent forces are employed in the blankholder, which is pressed normal and associated frictional tangent forces are employed in the blankholder, which is pressed against the flange until the resultant contact force with the blank reaches the prescribed value. As an example of sheet metal forming process coupling the effect of planar anisotropy and that of blankholding boundary condition, circular cup deep drawing has been analyzed considering both effects together.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.6
• /
• pp.1021-1029
• /
• 1992

A rational and straightforward method is introduced for developing continuum models of large platelike periodic lattice structures based on energy equivalence. The procedure for developing continuum plate models involves the use of existing well-defined finite element matrices for the easy calculation of strain and kinetic energies of a repeating cell, from which the reduced stiffness and mass matrices are obtained in terms of continuum degrees- of-freedom defined in this paper. The equivalent continuum plate properties are obtained from the direct comparison of the reduced matrices for continuum plate with those for lattice plate. The advantages of the present continuum method are that it may be applied to arbitrary lattice configurations and may give most diverse equivalent continuum plate properties including all kinds of coupling, while other methods may give only limited structural properties. To evaluate the continuum method developed in this paper, free vibration analyses for both of continuum and lattice plates are conducted. Numerical results show that the present continuum method gives very reliable structural and dynamic properties compared to other well-recognized methods.

• Interaction and multiscale mechanics
• /
• v.2 no.4
• /
• pp.353-374
• /
• 2009

A new seamless multiscale simulation was developed for coupling the continuum model with its molecular dynamics. Kriging-based Finite Element Method (K-FEM) is employed to model the continuum base of the entire domain, while the molecular dynamics (MD) is confined in a localized domain of interest. In the coupling zone, where the MD domain overlaps the continuum model, the overall Hamiltonian is postulated by contributions from the continuum and the molecular overlays, based on a quartic spline scaling parameter. The displacement compatibility in this coupling zone is then enforced by the Lagrange multiplier technique. A multiple-time-step velocity Verlet algorithm is adopted for its time integration. The validation of the present method is reported through numerical tests of one dimensional atomic lattice. The results reveal that at the continuum/MD interface, the commonly reported spurious waves in the literature are effectively eliminated in this study. In addition, the smoothness of the transition from MD to the continuum can be significantly improved by either increasing the size of the coupling zone or expanding the nodal domain of influence associated with K-FEM.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.22 no.2
• /
• pp.83-100
• /
• 1980

This study was intended to computerize the structural analysis of two-dimensional continuum by finite element method, and to provide a preparatory basis for more sophisticated and more generalized computer programs of this kind. A computer program, applicable to any shape of two-dimensional continuum, was formulated on the basis of 16-degree-of- freedom rectangular element. Various computational aspects pertaining to the implementation of finite element method were reviewed and settled in the course of programming. The validity of the program was checked through several case studies. To assess the accuracy and the convergence characteristics of the method, the results computed by the program were compared with solutions by other methods, namely the analytical Navier's method and the framework method. Through actual programming and analysis of the computed results, the following facts were recognized; 1) The stiffness matrix should necessarily be assembled in a condensed form in order to make it possible to discretize the continuum into practically adequate number of elements without using back-up storage. 2) For minimization of solution time, in-core solution of the equilibrium equation is essential. LDLT decomposition is recommended for stiffness matrices condensed by the compacted column storage scheme. 3) As for rectangular plates, the finite element method shows better performances both in the accuracy and in the rate of convergence than the framework method. As the number of elements increases, the error of the finite element method approaches around 1%. 4) Regardless of the structural shape, there is a uniform tendency in convergence characteristics dependent on the shape of element. Square elements show the best performance. 5) The accuracy of computation is independent of the interpolation function selected.