• Structural Engineering and Mechanics
• /
• v.75 no.6
• /
• pp.747-769
• /
• 2020

This paper provides a semi-analytical approach to investigate the variations of 3D displacement components, electric potential, stresses, electric displacements and transverse vibration frequencies in laminated piezoelectric composite plates based on the scaled boundary finite element method (SBFEM) and the precise integration algorithm (PIA). The proposed approach can analyze the static and dynamic responses of multilayered piezoelectric plates with any number of laminae, various geometrical shapes, boundary conditions, thickness-to-length ratios and stacking sequences. Only a longitudinal surface of the plate is discretized into 2D elements, which helps to improve the computational efficiency. Comparing with plate theories and other numerical methods, only three displacement components and the electric potential are set as the basic unknown variables and can be represented analytically through the transverse direction. The whole derivation is built upon the three dimensional key equations of elasticity for the piezoelectric materials and no assumptions on the plate kinematics have been taken. By virtue of the equilibrium equations, the constitutive relations and the introduced set of scaled boundary coordinates, three-dimensional governing partial differential equations are converted into the second order ordinary differential matrix equation. Furthermore, aided by the introduced internal nodal force, a first order ordinary differential equation is obtained with its general solution in the form of a matrix exponent. To further improve the accuracy of the matrix exponent in the SBFEM, the PIA is employed to make sure any desired accuracy of the mechanical and electric variables. By virtue of the kinetic energy technique, the global mass matrix of the composite plates constituted by piezoelectric laminae is constructed for the first time based on the SBFEM. Finally, comparisons with the exact solutions and available results are made to confirm the accuracy and effectiveness of the developed methodology. What's more, the effect of boundary conditions, thickness-to-length ratios and stacking sequences of laminae on the distributions of natural frequencies, mechanical and electric fields in laminated piezoelectric composite plates is evaluated.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.38 no.5
• /
• pp.74-84
• /
• 1996

Dynamic loading of structures often causes excursions of stresses well into the inelastic range and the influence of geometric changes on the dynamic response is also significant in many cases. Therefore, both material and geometric nonlinearity effects should be considered in case that a dynamic load acts on the structure. For developing a program to analyze the dynamic response of an axisymmetric shell in this study, the material nonlinearity effect on the dynamic response was formulated by the elasto-viscoplastic model highly corresponding to the real behavior of the material. Also, the geometrically nonlinear behavior is taken into account using a total Lagrangian coordinate system, and the equilibrium equation of motion was numerically solved by a central difference scheme. A complete finite element program has been developed and the results obtained by it are compared with those in the references 1 and 2. The results are in good agreement with each other. As a case study of its application, the developed program was applied to a dynamic response analysis of a nuclear reinforced concrete containment structure. The results obtained from the' numerical examples are summarized as follows : 1. The dynamic magnification factor of the displacement and the stress were unrelated with the concrete strength. 2. As shown by the results that the displacement dynamic magnification factor were form 1.7 to 2.3 and the stress dynamic magnification factor from 1.8 to 2.5, the dynamic magnification factor of stress were larger than that of displacement. 3. The dynamic magnification factor of stress on the exterior surface was larger than that on the interior surface of the structure.

• Structural Engineering and Mechanics
• /
• v.31 no.2
• /
• pp.181-210
• /
• 2009

This paper presents the theoretical developments of an exact finite strip for the buckling and initial post-buckling analyses of isotropic flat plates. The so-called exact finite strip is assumed to be simply supported out-of-plane at the loaded ends. The strip is developed based on the concept that it is effectively a plate. The present method, which is designated by the name Full-analytical Finite Strip Method in this paper, provides an efficient and extremely accurate buckling solution. In the development process, the Von-Karman's equilibrium equation is solved exactly to obtain the buckling loads and the corresponding form of out-of-plane buckling deflection modes. The investigation of thin flat plate buckling behavior is then extended to an initial post-buckling study with the assumption that the deflected form immediately after the buckling is the same as that obtained for the buckling. It is noted that in the present method, only one of the calculated out-of-plane buckling deflection modes, corresponding to the lowest buckling load, i.e., the first mode is used for the initial post-buckling study. Thus, the postbuckling study is effectively a single-term analysis, which is attempted by utilizing the so-called semi-energy method. In this method, the Von-Karman's compatibility equation governing the behavior of isotropic flat plates is used together with a consideration of the total strain energy of the plate. Through the solution of the compatibility equation, the in-plane displacement functions which are themselves related to the Airy stress function are developed in terms of the unknown coefficient in the assumed out-of-plane deflection function. These in-plane and out-of-plane deflected functions are then substituted in the total strain energy expressions and the theorem of minimum total potential energy is applied to solve for the unknown coefficient. The developed method is subsequently applied to analyze the initial postbuckling behavior of some representative thin flat plates for which the results are also obtained through the application of a semi-analytical finite strip method. Through the comparison of the results and the appropriate discussion, the knowledge of the level of capability of the developed method is significantly promoted.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.39 no.4
• /
• pp.106-113
• /
• 1997

Dynamic loading of structures often causes excursions of stresses well into the inelastic range, and the influence of the geometric changes on the dynamic response is also significant in many cases. Therefore, both material and geometric nonlinearity effects should be considered in case that a dynamic load acts on the structure. A structure in a nuclear power plant is a structure of importance which puts emphasis on safety. A nuclear container is a pressure vessel subject to internal pressure and this structure is constructed by a reinforced concrete or a pre-stressed concrete. In this study, the material nonlinearity effect on the dynamic response is formulated by the elasto-viscoplastic model highly corresponding to the real behavior of the material. Also, the geometrically nonlinear behavior is taken into account using a total Lagrangian coordinate system, and the equilibrium equation of motion is numerically solved by a central difference scheme. The constitutive relation of concrete is modeled according to a Drucker-Prager yield criterion in compression. The reinforcing bars are modeled by a smeared layer at the location of reinforcements, and the steel layer model under Von Mises yield criteria is adopted to represent an elastic-plastic behavior. To investigate the dynamic response of a nuclear reinforced concrete containment structure, the steel-ratios of 0, 3, 5 and 10 percent, are considered. The results obtained from the analysis of an example were summarized as follows 1. As the steel-ratio increases, the amplitude and the period of the vertical displacements in apex of dome decreased. The Dynamic Magnification Factor(DMF) was some larger than that of the structure without steel. However, the regular trend was not found in the values of DMF. 2. The dynamic response of the vertical displacement and the radial displacement in the dome-wall junction were shown that the period of displacement in initial step decreased with the steel-ratio increases. Especially, the effect of the steel on the dynamic response of radial displacement disapeared almost. The values of DMF were 1.94, 2.5, 2.62 and 2.66, and the values increased with the steel-ratio. 3. The characteristics of the dynamic response of radial displacement in the mid-wall were similar to that of dome-wall junction. The values of DMF were 1.91, 2.11, 2.13 and 2.18, and the values increased with the steel-ratio. 4. The amplitude and the period of the hoop-stresses in the dome, the dome-wall junction, and the mid-wall were shown the decreased trend with the steel-ratio. The values of DMF were some larger than those of the structure without steel. However, the regular trend was not found in the values of DMF.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.3
• /
• pp.305-319
• /
• 2000

Numerical solution lot frictional contact problems of nonlinearly deformable bodies having large deformation is presented. The contact conditions on the possible contact points are expressed by using the contact error vector, and the iterative scheme is used to reduce the contact error vector monotonically toward zero. At each iteration the solution consists of two steps : The first step is to revise the contact force by using the contact error vector given by the previous geometry, and the second step is to compute the displacement and the contact error vector by solving the equilibrium equation with the contact force given at the first step. Convergence of the iterative scheme to the correct solution is analyzed, and the numerical simulations we performed with a rigid-plastic membrane and a nonlinear elastic beam.

• Steel and Composite Structures
• /
• v.17 no.3
• /
• pp.287-303
• /
• 2014

8 steel reinforced concrete (SRC) columns with the encased steel ratio of 13.12% and 15.04% respectively were tested under the test axial load ratio of 0.33-0.80 and the low-frequency cyclic lateral loading. The cross sectional area of composite columns was $500mm{\times}500mm$. The mechanical properties, failure modes and deformabilities were studied. All the specimens produced flexure failure subject to combined axial force, bending moment and shear. Force-displacement hysteretic curves, strain curves of encased steels and rebars were obtained. The interaction behavior of encased steel and concrete were verified. The hysteretic curves of columns were plump in shapes. Hysteresis loops were almost coincident under the same levels of lateral loading, and bearing capacities did not change much, which indicated that the columns had good energy-dissipation performance and seismic capacity. Based on the equilibrium equation, the suggested practical calculation method could accurately predict the flexural strength of SRC columns with cross-shaped section encased steel. The obtained M-N curves of SRC columns can be used as references for further studies.

• Journal of Korean Society of Steel Construction
• /
• v.13 no.1
• /
• pp.41-52
• /
• 2001

In this paper the analysis of laminate composite folded plates with arbitrary angle connection like box girder is studied by finite strip method Total stiffness of laminated plate is obtained by integration of the stiffness in each layer or lamina through laminate thickness and total stiffness in each layer or lamina through laminate thickness and total tiffness matrix is obtained by substitutionto equilibrium equation derived from the minimum total potential energy theorem. The assumed displacement functions for a finite strip method in plate or box girder analysis are combinations of one-way polynomial functions in the transverse direction and harmonic functions in the span-wise direction. Finite strip method with the merits of the simplification in modeling and the reduction of analytical time is accurate in the analysis of laminate composite folded plates shaped like box firders.

• Structural Engineering and Mechanics
• /
• v.5 no.1
• /
• pp.105-113
• /
• 1997

This paper proposes a model for simulating concrete shrinkage taking into account aggregate restraint. In the model, concrete is regarded as a two-phase material based on shrinkage property. One is paste phase which undergoes shrinkage. Another is aggregate phase which is much more volumetrically stable. In the concrete, the aggregate phase is considered to restrain the paste shrinkage by particle interaction. Strain compatibility was derived under the assumption that there is no relative macroscopic displacement between both phases. Stresses on both phases were derived based on the shrinking stress of the paste phase and the resisting stress of the aggregate phase. Constitutive relation of paste phase was adopted from the study of Yomeyama, K. et al., and that of the aggregate phase was adopted from the author's particle contact density model. The equation for calculating concrete shrinkage considering aggregate restraint was derived from the equilibrium of the two phases. The concrete shrinkage was found to be affected by the free shrinkage of the paste phase, aggregate content and the stiffness of both phases. The model was then verified to be effective for simulating concrete shrinkage by comparing the predicted results with the autogeneous and drying shrinkage test results on mortar and concrete specimens.

• Smart Structures and Systems
• /
• v.20 no.5
• /
• pp.595-605
• /
• 2017

The present article reported the thermal buckling strength of the sandwich shell panel structure and subsequent improvement of the same by embedding shape memory alloy (SMA) fibre via a general higher-order mathematical model in conjunction with finite element method. The geometrical distortion of the panel structure due to the temperature is included using Green-Lagrange strain-displacement relations. In addition, the material nonlinearity of SMA fibre due to the elevated thermal environment also incorporated in the current analysis through the marching technique. The final form of the equilibrium equation is obtained by minimising the total potential energy functional and solved computationally with the help of an original MATLAB code. The convergence and the accuracy of the developed model are demonstrated by solving similar kind of published numerical examples including the necessary input parameter. After the necessary establishment of the newly developed numerical solution, the model is extended further to examine the effect of the different structural parameters (side-to-thickness ratios, curvature ratios, core-to-face thickness ratios, volume fractions of SMA fibre and end conditions) on the buckling strength of the SMA embedded sandwich composite shell panel including the different geometrical configurations.

• Structural Engineering and Mechanics
• /
• v.51 no.4
• /
• pp.685-705
• /
• 2014

This paper presents a new formulation for forward scalar wave simulations in semi-infinite media. Perfectly-Matched-Layers (PMLs) are used as a wave absorbing boundary layer to surround a finite computational domain truncated from the semi-infinite domain. In this work, a hybrid formulation was developed for the simulation of scalar wave motion in two-dimensional PML-truncated domains. In this formulation, displacements and stresses are considered as unknowns in the PML domain, while only displacements are considered to be unknowns in the interior domain. This formulation reduces computational cost compared to fully-mixed formulations. To obtain governing wave equations in the PML region, complex coordinate stretching transformation was introduced to equilibrium, constitutive, and compatibility equations in the frequency domain. Then, equations were converted back to the time-domain using the inverse Fourier transform. The resulting equations are mixed (contain both displacements and stresses), and are coupled with the displacement-only equation in the regular domain. The Newmark method was used for the time integration of the semi-discrete equations.