• Journal of Korean Society of Steel Construction
• /
• v.16 no.3 s.70
• /
• pp.295-303
• /
• 2004

In this paper, an improved four-node Reissner-Mindlin plate-bending element with enhanced assumed strain field is presented for the analysis of isotropic and laminated composite plates. To avoid the shear locking and spurious zero energy modes, the transverse shear behavior is improved by the addition of a new enhanced shear strain based on the incompatible displacement mode approach and bubble function. The "standard" enhanced strain fields (Andelfinger and Ramm, 1993) are also employed to improve the in-plane behaviors of the plate elements. The four-node quadrilateral element derived using the first-order shear deformation theory is designated as "14EASP". Several applications are investigated to assess the features and the performances of the proposed element. The results are compared with other finite element solutions and analytical solutions. Numerical examples show that the element is stable, invariant, passes the patch test, and yields good results especially in highly distorted regimes.

• Steel and Composite Structures
• /
• v.46 no.4
• /
• pp.471-483
• /
• 2023

This paper presents an analytical hyperbolic theory based on the refined shear deformation theory for mechanical stability analysis of the simply supported advanced composites plates (exponentially, sigmoidal and power-law graded) under triangular, trapezoidal and uniform uniaxial and biaxial loading. The developed model ensures the boundary condition of the zero transverse stresses at the top and bottom surfaces without using the correction factor as first order shear deformation theory. The mathematical formulation of displacement contains only four unknowns in which the transverse deflection is divided to shear and bending components. The current study includes the effect of the geometric imperfection of the material. The modeling of the micro-void presence in the structure is based on the both true and apparent density formulas in which the porosity will be dense in the mid-plane and zero in the upper and lower surfaces (free surface) according to a logarithmic function. The analytical solutions of the uniaxial and biaxial critical buckling load are determined by solving the differential equilibrium equations of the system with the help of the Navier's method. The correctness and the effectiveness of the proposed HyRPT is confirmed by comparing the results with those found in the open literature which shows the high performance of this model to predict the stability characteristics of the FG structures employed in various fields. Several parametric analyses are performed to extract the most influenced parameters on the mechanical stability of this type of advanced composites plates.

• Atmosphere
• /
• v.25 no.2
• /
• pp.339-351
• /
• 2015

The purpose of this study is for the fundamental understandings about building morphological parameters and aerodynamic roughness parameters of Seoul, Korea using the detailed urban geographic information datasets. Applied roughness parameter calculations are based on a digital map of buildings with lot area polygons. The quality of the developed roughness length ($z_0$) of Seoul was evaluated with densely installed 107 automatic weather stations. The correlation coefficient results between averaged wind speeds of AWS data and averaged $z_0$ is -0.303 in night and -0.398 in day (200 m radii circles case). Further $z_0$ enhancement should follow by considering other surface features such as high tree and orography of Seoul. However, this study would meet the needs to for local- or meso-scale meteorological modeling applications of Seoul. However, further studies would require for enhancing the $z_0$ applications of Seoul.

• Structural Engineering and Mechanics
• /
• v.8 no.3
• /
• pp.233-242
• /
• 1999

In the present study, the elasto-plastic analysis of prismatic plate structures subjected to pure bending is carried out using the finite strip method. The end cross-sections of the structure are assumed to remain plane during deformation, and the compatibility along corner lines is ensured by choosing proper displacement functions. The effects of both the initial geometrical imperfections and residual stresses due to fabrication are included in the combined geometrically and materially nonlinear simulation. The von-Mises yield criterion and the Prandtl-Reuss flow theory of plasticity are applied in modelling the elasto-plastic behavior of material. Newton-Raphson iterations are carried out as the rotation of the end cross sections of the structure is increased step by step. The parameter representing the overall axial strain of structure is adjusted constantly during the iteration process in order to eliminate the resulting overall axial force on any cross-section of the structure in correspondence with the assumption of zero axial force in pure bending. Several numerical examples are presented to validate the present method and to investigate the effects of some material and geometrical parameters.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.16 no.4
• /
• pp.37-52
• /
• 2012

In this study, a new 4-node general shell element with 6 DOFs per node is presented. Drilling rotational degrees of freedom are introduced by the variational principle with an independent rotation field. In formulation of the element, substitute transverse shear strain fields are used to avoid shear locking, while four nonconforming modes are applied in the in-plane displacement fields as a remedy for membrane locking. In addition, a direct modification method for nonconforming modes is employed in the numerical implementation of nonconforming modes to represent constant strain states. A 9-points integration rule is adopted for volume integration in the computation of the element stiffness matrix. With the combined use of these techniques, the developed shell element has no spurious zero energy modes, and can represent a constant strain state. Several numerical tests are carried out to evaluate the performance of the new element developed. The test results show that the behavior of the elements is satisfactory.

• The Journal of Engineering Geology
• /
• v.18 no.4
• /
• pp.447-457
• /
• 2008

Laboratory tests for single plane sliding were conducted using the model rock slope to investigate the cut slope deformability and failure mechanism due to combined effect of engineering characteristics such as angle of sliding plane, water force, joint roughness and infillings. Also the possibility of prediction of slope failure through displacement monitoring was explored. The joint roughness was prepared in forms of saw-tooth type having different roughness specifications. The infillings was maintained between upper and lower roughness plane from zero to 1.2 times of the amplitude of the surface projections. Water force was expressed as the percent filling of tension crack from dry (0%) to full (100%), and constantly increased from 0% at the rate of 0.5%/min and 1%/min upto failure. Total of 50 tests were performed at sliding angles of $30^{\circ}$ and $35^{\circ}$ based on different combinations of joint roughness, infilling thickness and water force increment conditions. For smooth sliding plane, it was found that the linear type of deformability exhibited irrespective of the infilling thickness and water force conditions. For sliding planes having roughness, stepping or exponential types of deformability were predominant under condition that the infilling thickness is lower or higher than asperity height, respectively. These arise from the fact that, once the infilling thickness exceeds asperities, strength and deformability of the sliding plane is controlled by the engineering characteristics of the infilling materials. The results obtained in this study clearly show that the water force at failure was found to increase with increasing joint roughness, and to decrease with increasing filling thickness. It seems possible to estimate failure time using the inverse velocity method for sliding plane having exponential type of deformability. However, it is necessary to estimate failure time by trial and error basis to predict failure of the slope accurately.

• Wind and Structures
• /
• v.11 no.4
• /
• pp.275-289
• /
• 2008

Large eddy simulations have been performed within and over different types of urban building arrays. This paper adopted three dimensionless parameters, building frontal area density (${\lambda}_f$) the variation degree of building height (${\sigma}_h$), and the staggered degree of building range ($r_s$), to study the systematic influence of building spacing, height and layout on wind and turbulent characteristics. The following results have been achieved: (1) As ${\lambda}_f$ decrease from 0.25 to 0.18, the mean flow patterns transfer from "skimming" flow to "wake interference" flow, and as ${\lambda}_f$ decrease from 0.06 to 0.04, the mean flow patterns transfer from "wake interference" flow to "isolated roughness" flow. With increasing ${\lambda}_f$, wind velocity within arrays increases, and the vortexes in front of low buildings would break, even disappear, whereas the vortexes in front of tall buildings would strengthen and expand. Tall buildings have greater disturbance on wind than low buildings do. (2) All the wind velocity profiles and the upstream profile converge at the height of 2.5H approximately. The decay of wind velocity within the building canopy was in positive correlation with ${\lambda}_f$ and $r_s$. If the height of building arrays is variable, Macdonald's wind velocity model should be modified through introducing ${\sigma}_h$, because wind velocity decreases at the upper layers of the canopy and increases at the lower layers of the canopy. (3) The maximum of turbulence kinetic energy (TKE) always locates at 1.2 times as high as the buildings. TKE within the canopy decreases with increasing ${\lambda}_f$ and $r_s$ but the maximum of TKE are very close though ${\sigma}_h$ varies. (4) Wind velocity profile follows the logarithmic law approximately above the building canopy. The Zero-plane displacement $z_d$ heighten with increasing ${\lambda}_f$, whereas the maximum of and Roughness length $z_0$ occurs when ${\lambda}_f$ is about 0.14. $z_d$ and $z_0$ heighten linearly with ${\sigma}_h$ and $r_s$, If ${\sigma}_h$ is large enough, $z_d$ may become higher than the average height of buildings.