• Journal of Power System Engineering
• /
• v.3 no.1
• /
• pp.74-79
• /
• 1999

This paper describes formulation for algorithm of time historical response analysis of vibration for straight-line structure. This method is derived from a combination of the transfer stiffness coefficient method and the Newmark method. And this present method improves the computational accuracy of the transient vibration response analysis remarkably owing to several advantages of the transfer stiffness coefficient method. We regarded the structure as a lumped mass system here. The analysis algorithm for the time historical response was formulated for the straight-line structure containing crooked, tree type system. The validity of the present method compared with the transfer matrix method and the Finite Element Method for transient vibration analysis is demonstrated through the numerical computations.

• Journal of Mechanical Science and Technology
• /
• v.19 no.12
• /
• pp.2187-2196
• /
• 2005

In this paper, the stiffness and the mass matrices for the in-plane motion of a thin circular beam element are derived respectively from the strain energy and the kinetic energy by using the natural shape functions of the exact in-plane displacements which are obtained from an integration of the differential equations of a thin circular beam element in static equilibrium. The matrices are formulated in the local polar coordinate system and in the global Cartesian coordinate system with the effects of shear deformation and rotary inertia. Some numerical examples are performed to verify the element formulation and its analysis capability. The comparison of the FEM results with the theoretical ones shows that the element can describe quite efficiently and accurately the in-plane motion of thin circular beams. The stiffness and the mass matrices with respect to the coefficient vector of shape functions are presented in appendix to be utilized directly in applications without any numerical integration for their formulation.

• Proceeding of KASS Symposium
• /
• 2006.05a
• /
• pp.256-264
• /
• 2006

This paper describes the formulation of rectangular flat shell element that is modeled with the six degree of freedoms including a rotational degree of freedom. The rectangular finite element matrix with a rotational degree of freedom is developed using a beam stiffness matrix and compared with other methods. The outputs of the quantity of vertical deflection of cantilever beam show us the improving evidence of the Frame-Shell finite element matrix in a calculation of vertical deflections of cantilever beam.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.11 no.7
• /
• pp.2661-2669
• /
• 2010

In this study, a finite element formulation based first-order shear deformation theory is developed for non-linear behaviors of laminated composite plates containing matrix cracking. The multi-directional stiffness degradation is developed for adopting the stiffness variation induced from matrix cracking, which is proposed by Duan and Yao. The matrix cracking can be expressed in terms of the variation of material properties, such as Young's modulus, shear modulus and Possion ratio of plates, and sequently it is possible to predict the variation of the local stiffness. Using the assumed natural strain method, the present shell element generates neither membrane nor shear locking behavior. Numerical examples demonstrate that the present element behaves quite satisfactorily either for the linear or geometrical nonlinear analysis of laminated composite plates. The results of laminated composite plates with matrix cracking may be the benchmark test for the non-linear analysis of damaged laminated composite plates.

• Proceedings of the KSME Conference
• /
• 2001.06b
• /
• pp.355-360
• /
• 2001

By using Frenet formulation and Timoshenko beam theory, the partial differential equations of motion are derived for a helical spring having a doubly symmetrical cross section subjected to the pre-load axially. These equations of motion are solved to give the dispersion relationship and dynamic stiffness matrix is assembled. Natural frequencies are obtained from the receptance of the system. The results of the dynamic stiffness method are compared with those of the transfer matrix method from published examples and finite element method.

• Architectural research
• /
• v.13 no.1
• /
• pp.17-24
• /
• 2011

Isogeometric solution of Poisson equation is provided. NURBS (NonUniform B-spline Surface) is introduced to express both geometry of structure and unknown field of governing equation. The terms of stiffness matrix and load vector are consistently derived with very accurate geometric definition. The validity of the isogeometric formulation is demonstrated by using two numerical examples such as square plate and L-shape plate. From numerical results, the present solutions have a good agreement with analytical and finite element (FE) solutions with the use of a few cells in isogeometric analysis.

• Steel and Composite Structures
• /
• v.39 no.5
• /
• pp.543-564
• /
• 2021

In this study free vibration analysis of a cracked Goland composite wing is investigated. The wing is modelled as a cantilevered beam based on Euler- Bernoulli equations. Also, composite material is modelled based on lamina fiber-reinforced. Edge crack is modelled by additional boundary conditions and local flexibility matrix in crack location, Castigliano's theorem and energy release rate formulation. Governing differential equations are extracted by Hamilton's principle. Using the separation of variables method, general solution in the normalized form for bending and torsion deflection is achieved then expressions for the cross-sectional rotation, the bending moment, the shear force and the torsional moment for the cantilevered beam are obtained. The cracked beam is modelled by separation of beam into two interconnected intact beams. Free vibration analysis of the beam is performed by applying boundary conditions at the fixed end, the free end, continuity conditions in the crack location of the beam and dynamic stiffness matrix determinant. Also, the effects of various parameters such as length and location of crack and fiber angle on natural frequencies and mode shapes are studied. Modal analysis results illustrate that natural frequencies and mode shapes are affected by depth and location of edge crack and coupling parameter.

• Tunnel and Underground Space
• /
• v.10 no.2
• /
• pp.173-179
• /
• 2000

The influence of the plastic flow rules on the elasto-plastic behaviour of a discrete joint element was investigated by performing the numerical direct shear tests under both constant normal displacement and normal displacement conditions. The finite interface elements obeying Plesha’s joint constitutive law was used to allow the relative motion of the rock blocks on the joint surface. Realistic results were obtained in the tests adopting the non-associated flow rule, while the associated flow rule overestimated the joint dilation. To overcome the computational drawbacks coming from the non-symmetric element stiffness matrix in the conventional non-associated plasticity, the symmetric formulation of the tangential stiffness matrix for a non-associated joint element was proposed. The symmetric elasto-plastic matrix it derived by assuming an imaginary equivalent joint with associated flow rule which shows the same plastic response as that of original Joint with non-associated flow rule. The validity of the formulation was confirmed through the numerical direct shear tests under constant normal stress condition.

• Structural Engineering and Mechanics
• /
• v.17 no.1
• /
• pp.113-140
• /
• 2004

A mixed eight-node hexahedral element formulated via the Hu-Washizu principle as well as the field extrapolation technique is presented. The mixed element with only three translational degrees of freedom at each node can provide extremely accurate and reliable performance for popular benchmark problems such as spacial beams, plates, shells as well as general three-dimensional elasticity problems. Numerical calculations also show that when extremely skewed and coarse meshes and nearly incompressible materials are used, the proposed mixed element can still possess excellent behaviour. The mixed formulation starts with introduction of a parallelepiped domain associated with the given general eight-node hexahedral element. Then, the assumed strain field at the nodal level is constructed via the Hu-Washizu variational principle for that associated parallelepiped domain. Finally, the assumed strain field at the nodal level of the given hexahedral element is established by using the field extrapolation technique, and then by using the trilinear shape functions the assumed strain field of the whole element domain is obtained. All matrices involved in establishing the element stiffness matrix can be evaluated analytically and expressed explicitly; however, a 24 by 24 matrix has to be inverted to construct the displacement extrapolation matrix. The proposed hexahedral element satisfies the patch test as long as the element with a shape of parallelepiped.

• Interaction and multiscale mechanics
• /
• v.1 no.2
• /
• pp.231-250
• /
• 2008

The statistical two-order and two-scale method is developed for predicting the mechanics parameters, such as stiffness and strength of core-shell particle-filled polymer composites. The representation and simulation on meso-configuration of random particle-filled polymers are stated. And the major statistical two-order and two-scale analysis formulation is briefly given. The two-order and two-scale expressions for the strains and stresses of conventionally strength experimental components, including the tensional or compressive column, the twist bar and the bending beam, are developed by means of their classical solutions with orthogonal-anisotropic coefficients. Then a new effective mesh generation algorithm is presented. The mechanics parameters of core-shell particle-filled polymer composites, including the expected stiffness parameters, minimum stiffness parameters, and the expected elasticity limit strength and the minimum elasticity limit strength, are defined by means of the stiffness coefficients and elasticity strength criterions for core, shell and matrix. Finally, the numerical results for predicting both stiffness and elasticity limit strength parameters are compared with the experimental data.