• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.319-326
• /
• 2002

A p-version finite element model based on degenerate shell element is proposed for the analysis of orthotropic laminated plates. In the nonlinear formulation of the model, the total Lagrangian formulation is adopted with large deflection and moderate rotation being accounted for in the sense of von Karman hypothesis. The material model Is based on the Huber-Mises yield criterion and Prandtl-Reuss flow rule in accordance with the theory of strain hardening yield function, which is generalized for anisotropic materials by introducing the parameters of anisotropy. The model is also based on extension of equivalent-single layer laminate theory(ESL theory) with shear deformation, leading to continuous shear strain at the interface of two layers. The Integrals of Legendre Polynomials we used for shape functions with p-level varying from 1 to 10. Gauss-Lobatto numerical quadrature is used to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed p-version finite element model is demonstrated through several comparative points of view in terms of ultimate load, convergence characteristics, nonlinear effect, and shape of plastic zone

• International Journal of Aeronautical and Space Sciences
• /
• v.2 no.2
• /
• pp.65-72
• /
• 2001

This paper demonstrates an explicit-implicit, finite element analysis for linear as well as nonlinear hygrothermal stress problems. Additional features, such as moisture diffusion equation, crack element and virtual crack extension(VCE) method for evaluating J-integral are implemented in this program. The Linear Elastic Fracture Mechanics(LEFM) Theory is employed to estimate the crack driving force under the transient condition for an existing crack. Pores in materials are assumed to be saturated with moisture in the liquid form at the room temperature, which may vaporize as the temperature increases. The vaporization effects on the crack driving force are also studied. The ideal gas equation is employed to estimate the thermodynamic pressure due to vaporization at each time step after solving basic nodal values. A set of field equations governing the time dependent response of porous media are derived from balance laws based on the mixture theory. Darcy's law is assumed for the fluid flow through the porous media. Perzyna's viscoplastic model incorporating the Von-Mises yield criterion are implemented. The Green-Naghdi stress rate is used for the invariant of stress tensor under superposed rigid body motion. Isotropic elements are used for the spatial discretization and an iterative scheme based on the full Newton-Raphson method is used for solving the nonlinear governing equations.

• International Journal of Automotive Technology
• /
• v.7 no.3
• /
• pp.353-360
• /
• 2006

The nominal structural stress calculation method proposed by Radaj has included some problems as follows: (a) How the value of the diameter D is decided in the method; (b) It is not possible to estimate nominal structural stress of the spot welded joints with the balanced sheet in-plane load that no general loads are obtained by FE shell analysis. In this paper, the new method for calculating nominal structural stress was proposed to solve above-mentioned problems. The proposed method calculates the nominal structural stress through the circular plate theory in theory of elasticity. This theoretical analysis uses not only general loads but also nodal displacements around spot welding provided by FE shell analysis as boundary condition. Fatigue test data of various spot-welded joints could be arranged in a narrow bandwidth on S-N chart using the nominal structural stresses calculated by proposed method. The fatigue life prediction method using the proposed method for calculating nominal structural stress is useful for the prior evaluation technique that can predict the fatigue life of spot welding by CAE.

• Journal of Korean Association for Spatial Structures
• /
• v.3 no.2 s.8
• /
• pp.85-90
• /
• 2003

If we use a third order approximation for the displacement function of beam element in finite element methods, finite element solutions of beams yield nodal displacement values matching to beam theory results to have no connection with the number increasing of elements of beams. It is assumed that, as the member displacement value at beam nodes are correct, the calculation procedure of beam element stiffness matrix have no numerical errors. A the member forces are calculated by the equations of $\frac{-M}{EI}=\frac{{d^2}{\omega}}{dx^2}\;and\;\frac{dM}{dx}=V$, the member forces at nodes of beams have errors in a moment and a shear magnitudes in the case of smaller number of element. The nodal displacement value of plate subject to the lateral load converge to the exact values according to the increase of the number of the element. So it is assumed that the procedures of plate element stiffness matrix calculations has a error in the fundamental assumptions. The beam methods for the high accuracy ratio solution Is also applied to the plate analysis. The method of reducing a error ratio of member forces and element stiffness matrix in the finite element methods is studied. Results of study were as follows. 1. The matrixes of EI[B] and [K] in the equations of M(x)=EI[B]{q} and M(x) = [K]{q}+{Q} of beams are same. 2. The equations of $\frac{-M}{EI}=\frac{{d^2}{\omega}}{dx^2}\;and\;\frac{dM}{dx}=V$ for the member forces have a error ratio in a finite element method of uniformly loaded structures, so equilibrium node loads {Q} must be substituted in the equation of member forces as the numerical examples of this paper revealed.

• Advances in aircraft and spacecraft science
• /
• v.3 no.3
• /
• pp.351-366
• /
• 2016

Macroperforations improve the sound absorption performance of porous materials in acoustic cavities and in waveguides. In an acoustic cavity, enhanced noise reduction is achieved using porous materials having macroperforations. Double porosity materials are obtained by filling these macroperforations with different poroelastic materials having distinct physical properties. The locations of macroperforations in porous layers can be chosen based on cavity mode shapes. In this paper, the effect of variation of macroporosity and double porosity in porous materials on noise reduction in an acoustic cavity is presented. This analysis is done keeping each perforation size constant. Macroporosity of a porous material is the fraction of area covered by macro holes over the entire porous layer. The number of macroperforations decides macroporosity value. The system under investigation is an acoustic cavity having a layer of poroelastic material rigidly attached on one side and excited by an internal point source. The overall sound pressure level (SPL) inside the cavity coupled with porous layer is calculated using mixed displacement-pressure finite element formulation based on Biot-Allard theory. A 32 node, cubic polynomial brick element is used for discretization of both the cavity and the porous layer. The overall SPL in the cavity lined with porous layer is calculated for various macroporosities ranging from 0.05 to 0.4. The results show that variation in macroporosity of the porous layer affects the overall SPL inside the cavity. This variation in macroporosity is based on the cavity mode shapes. The optimum range of macroporosities in poroelastic layer is determined from this analysis. Next, SPL is calculated considering periodic and nodal line based optimum macroporosity. The corresponding results show that locations of macroperforations based on mode shapes of the acoustic cavity yield better noise reduction compared to those based on nodal lines or periodic macroperforations in poroelastic material layer. Finally, the effectiveness of double porosity materials in terms of overall sound pressure level, compared to equivolume double layer poroelastic materials is investigated; for this the double porosity material is obtained by filling the macroperforations based on mode shapes of the acoustic cavity.

• Nuclear Engineering and Technology
• /
• v.25 no.2
• /
• pp.193-203
• /
• 1993

A new heuristic method for loading pattern search has been developed to overcome short-comings of the algorithmic approach. To reduce the size of vast solution space, general shuffling rules, a regionwise shuffling method, and a pattern grouping method were introduced. The entropy theory was applied to classify possible loading patterns into groups with similarity between them. The pattern search program was implemented with use of the PROLOG language. A two-group nodal code MEDIUM-2D was used for analysis of power distribution in the core. The above mentioned methodology has been tested to show effectiveness in reducing of solution space down to a few hundred pattern groups. Burnable poison rods were then arranged in each pattern group in accordance with burnable poison distribution rules, which led to further reduction of the solution space to several scores of acceptable pattern groups. The method of maximizing cycle length(MCL) and minimizing power-peaking factor(MPF) were applied to search for specific useful loading patterns from the acceptable pattern groups. Thus, several specific loading patterns that have low power-peaking factor and large cycle length were successfully searched from the selected pattern groups.

• Structural Engineering and Mechanics
• /
• v.48 no.1
• /
• pp.77-91
• /
• 2013

The paper presents the formulation of 3-nodal line semi-analytical Mindlin-Reissner finite strip in the buckling analysis of thin-walled members, which are subjected to arbitrary loads. The finite strip is simply supported in two opposite edges. The general loading and in-plane rotation techniques are used to develop this finite strip. The linear stiffness matrix and the geometric stiffness matrix of the finite strip are given in explicit forms. To validate the proposed model and study its performance, numerical examples of some thin-walled sections have been performed and the results obtained have been compared with finite element models and the published ones.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.28 no.12
• /
• pp.2019-2031
• /
• 2004

The least-squares formulation for rigid-plasticity based on J$_2$-flow rule and infinitesimal theory and its meshfree implementation using moving least-squares approximation are proposed. In the least-squares formulation the squared residuals of the constitutive and equilibrium equations are minimized. Those residuals are represented in a form of first-order differential system using the velocity and stress components as independent variables. For the enforcement of the boundary and frictional contact conditions, penalty scheme is employed. Also the reshaping of nodal supports is introduced to avoid the difficulties due to the severe local deformation near the contact interface. The proposed least-squares meshfree method does not require any structure of extrinsic cells during the whole process of analysis. Through some numerical examples of metal forming processes, the validity and effectiveness of the method are investigated.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.31 no.3
• /
• pp.233-245
• /
• 2011

A new BWIM(bridge weigh-in-motion) algorithm using both measured strain and acceleration data is proposed. To consider the effects of bridge vibration on the estimation of moving loads, the dynamic governing equation is applied with the known stiffness and mass properties but damping is ignored. Dynamic displacements are computed indirectly from the measured strains using the beam theory and accelerations are measured directly by accelerometers. To convert a unit moving load to its equivalent nodal force, a transformation matrix is determined. The incompleteness in the measured responses is considered in developing the algorithm. To examine the proposed BWIM algorithm, simulation studies, laboratory experiments and field tests were carried. In the simulation study, effects of measurement noise and estimation error in the vehicle speed on the results were investigated.

• Bulletin of the Korean Space Science Society
• /
• 1992.10a
• /
• pp.11-11
• /
• 1992

The most celebrated problem in artificial satellite theory is undoubtedly the critical inclination problem. The KITSAT-1 satellite launched by Arian 42P from Guiana in August 11, 1992 has orbital inclination close to the critical value cos-1(1/rs). In that case, there is a singularity in some perturbation terms and therefore perigee will be fixed because do/dt is equals to zero. But actually the long periodic behaviour in argument of perigee, u is affected by luni-solar gravity andrelativistic effect, etc. Luni-solar gravity induces periodic perturbations in allorbital elements except the semi-major axis, and secular variations in 0, u and M.We have obtained nodal rate and inclination variations in case of the KITSAT-1.In this paper, we will also show the geopotential perturbations in 0, M of the Satellite in CIUding J2, J3, J4, J22 terms.