• Structural Engineering and Mechanics
• /
• v.61 no.1
• /
• pp.15-29
• /
• 2017

The recently introduced index Ratio Of Torsion (ROT) quantifies the base shear amplification due to torsional effects on shear cantilever types of building structures. In this work, a theoretical proof based on the theory of elasticity is provided, depicting that the ratio of torsion (ROT) is independent of the forces acting on the structure, although its definition stems from the shear forces. This is a particular attribute of other design and evaluation criteria against torsion such as center of rigidity and center of strength. In the case of ROT, this evidence could be considered as inconsistent, as ROT is a function solely of the forces acting on structural members, nevertheless it is proven to be independent of them. As ROT is the amplification of the shear forces due to in-plan irregularities, this work depicts that this increase of internal shear forces rely only on the structural topology. Moreover, a numerical verification of this theoretical finding was accomplished, using linear statistics interpretation and nonlinear neural networks simulation for an adequate database of structures.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.6
• /
• pp.615-626
• /
• 2010

The discrete-layer model is proposed to analyze the edge-effect problem of laminates under extension and flexure. Based on three-dimensional elasticity theory, the displacement fields of each layer in a laminate have been treated discretely in terms of three displacement components across the thickness. The displacement fields at bottom and top surfaces within a layer are approximated by two-dimensional shape functions. Then two surfaces are connected by one-dimensional high order shape functions. Thus the p-convergent refinement on approximated one- and two-dimensional shape functions can be implemented independently of each other. The quality of present model is mostly determined by polynomial degrees of shape functions for given displacement fields. For nodal modes with physical meaning, the linear Lagrangian polynomials are considered. Additional modes without physical meaning, which are created by increasing nodeless degrees of shape functions, are derived from integrals of Legendre polynomials which have an orthogonality property. Also, it is assumed that mapping functions are linear in the light of shape of laminated plates. The results obtained by this proposed model are compared with those available in literatures. Especially, three-dimensional out-of-plane stresses in the interior and near the free edges are evaluated and convergence performance of the present model is established with the stress results.

• Structural Engineering and Mechanics
• /
• v.40 no.2
• /
• pp.257-277
• /
• 2011

In this paper the geometrically nonlinear continuum plate finite element model, hitherto not reported in the literature, is developed using the total Lagrange formulation. With the layerwise displacement field of Reddy, nonlinear Green-Lagrange small strain large displacements relations (in the von Karman sense) and linear elastic orthotropic material properties for each lamina, the 3D elasticity equations are reduced to 2D problem and the nonlinear equilibrium integral form is obtained. By performing the linearization on nonlinear integral form and then the discretization on linearized integral form, tangent stiffness matrix is obtained with less manipulation and in more consistent form, compared to the one obtained using laminated element approach. Symmetric tangent stiffness matrixes, together with internal force vector are then utilized in Newton Raphson's method for the numerical solution of nonlinear incremental finite element equilibrium equations. Despite of its complex layer dependent numerical nature, the present model has no shear locking problems, compared to ESL (Equivalent Single Layer) models, or aspect ratio problems, as the 3D finite element may have when analyzing thin plate behavior. The originally coded MATLAB computer program for the finite element solution is used to verify the accuracy of the numerical model, by calculating nonlinear response of plates with different mechanical properties, which are isotropic, orthotropic and anisotropic (cross ply and angle ply), different plate thickness, different boundary conditions and different load direction (unloading/loading). The obtained results are compared with available results from the literature and the linear solutions from the author's previous papers.

• Steel and Composite Structures
• /
• v.44 no.3
• /
• pp.371-388
• /
• 2022

The present study tackles the problem of forced vibration of imperfect axially functionally graded shell structure with truncated conical geometry. The linear and nonlinear large-deflection of the structure are considered in the mathematical formulation using von-Kármán models. Modified coupled stress method and principle of minimum virtual work are employed in the modeling to obtain the final governing equations. In addition, formulations of classical elasticity theory are also presented. Different functions, including the linear, convex, and exponential cross-section shapes, are considered in the grading material modeling along the thickness direction. The grading properties of the material are a direct result of the porosity change in the thickness direction. Vibration responses of the structure are calculated using the semi-analytical method of a couple of homotopy perturbation methods (HPM) and the generalized differential quadrature method (GDQM). Contradicting effects of small-scale, porosity, and volume fraction parameters on the nonlinear amplitude, frequency ratio, dynamic deflection, resonance frequency, and natural frequency are observed for shell structure under various boundary conditions.

• Advances in Computational Design
• /
• v.9 no.1
• /
• pp.39-52
• /
• 2024

If the governing differential equation arising from engineering problems is treated as an analytic, continuous and derivable function, it can be expanded by one point as a series of finite numbers. For the function to be zero for each value of its domain, the coefficients of each term of the same power must be zero. This results in a recursive relationship which, after applying the natural conditions or the boundary conditions, makes it possible to obtain the values of the derivatives of the function with acceptable accuracy. The elastoplastic analysis of an inhomogeneous thick sphere of metallic materials with linear variation of the modulus of elasticity, yield stress and Poisson's ratio as a function of radius subjected to internal pressure is presented. The Beltrami-Michell equation is established by combining equilibrium, compatibility and constitutive equations. Assuming axisymmetric conditions, the spherical coordinate parameters can be used as principal stress axes. Since there is no analytical solution, the natural boundary conditions are applied and the governing equations are solved using a proposed new method. The maximum effective stress of the von Mises yield criterion occurs at the inner surface; therefore, the negative sign of the linear yield stress gradation parameter should be considered to calculate the optimal yield pressure. The numerical examples are performed and the plots of the numerical results are presented. The validation of the numerical results is observed by modeling the elastoplastic heterogeneous thick sphere as a pressurized multilayer composite reservoir in Abaqus software. The subroutine USDFLD was additionally written to model the continuous gradation of the material.

• International Journal of Concrete Structures and Materials
• /
• v.4 no.1
• /
• pp.37-43
• /
• 2010

Theoretical models based on modern interpretations of the morphology and interactions of cement hydration products are developed for prediction of the mechanical properties of hydrated cement paste (hcp). The models are based on the emerging nanostructural vision of calcium silicate hydrate (C-S-H) morphology, and account for the intermolecular interactions between nano-scale calcium C-S-H particles. The models also incorporate the effects of capillary porosity and microcracking within hydrated cement paste. The intrinsic modulus of elasticity and tensile strength of hydrated cement paste are determined based on intermolecular interactions between C-S-H nano-particles. Modeling of fracture toughness indicates that frictional pull-out of the micro-scale calcium hydroxide (CH) platelets makes major contributions to the fracture energy of hcp. A tensile strength model was developed for hcp based on the linear elastic fracture mechanics theories. The predicted theoretical models are in reasonable agreements with empirical models developed based on the experimental performance of hcp.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.7 no.1
• /
• pp.52-63
• /
• 1983

The purpose of the study is to find development of contact area, contact pressure and friction forces occurring at joints or connection areas inbetween structural members or mechanical parts. The problem has a pair of difficulties intrinsically; a constraint of displacement due to contact, and presence of work term by nonconservative friction force in the variational principle of the problem. Because of these difficulties, the variational principle remains in the form of inequality. It is resolved by penalty method and perturbation method making the inequality to an equality which is proper for computational purposes. A contact problem without friction is solved to find contact area and contact pressure, which are to be used as data for the analysis of the friction problem using perturbed variational principle. For numerical experiments, a Hertz problem, a rigid punch problem, and the latter one with friction effects are solved using $Q_2$-finite elements.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.33 no.2
• /
• pp.135-144
• /
• 2009

In this present work, spline FEA for the trimmed NURBS surface of the 2D linear elasticity problem is presented. The main benefit of the proposed method is that no additional efforts for modeling of trimmed NURBS surfaces are needed and the information of the trimming curves and trimmed surfaces exported from the CAD system can be directly used for analysis. For this, trimmed elements are searched by using NURBS projection scheme. The integration of the trimmed elements is performed by using the NURBS-enhanced integration scheme. The formulation of constructing stiffness matrix of trimmed elements is presented. In this formulation, the information of the trimming curve is used for calculating the Jacobian as well as for obtaining integration points. The robustness and effectiveness of the proposed method are investigated through various numerical examples.

• Journal of the Korean Society for Precision Engineering
• /
• v.23 no.12 s.189
• /
• pp.80-87
• /
• 2006

Seals are usually made from elastomer, a kind of rubber, and it has the non-linearity and hyper-elasticity. U-type seals are used to prevent the leakage of internal fluid sealed in hydraulic actuator because they have more excellent performance than O-rings or rectangular seals. As a core part of hydraulic actuator, U-type seal gives much effect on performance and reliability of actuator. This study considers an NBR U-type seal under high pressure of a hydraulic actuator, and provides its deformation, stress-strain characteristic and contact force using the non-linear finite element analysis. Analysis results are compared with the experimental ones performed by the self-developed testing equipment. Verification result shows that this study presents a good application process for the effective design of U-type seals under high operation pressure.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2002.11b
• /
• pp.190-197
• /
• 2002

A new design methodology is presented for the multi-degree-of-freedom vibration absorber for an elastically suspended rigid body with planes of symmetry in general motion. Unlike the common single degree-of-freedom vibration absorber, the presented methodology makes use of both linear and rotational properties of the absorber. It is suggested that an absorber is designed separately for the in-plane and out-of-plane vibration modes and thereby combined the two cases for a six-degree-of-freedom absorber. The nine possible design methods are suggested for the six-degree-of-freedom absorber when an elastically suspended rigid body has one, two, or three planes of symmetry.