• Structural Engineering and Mechanics
• /
• v.85 no.5
• /
• pp.621-633
• /
• 2023

The major objective of this paper is to study the receding contact problem between two functional graded layers under a flat indenter. The gravity is assumed negligible, and the shear moduli of both layers are assumed to vary exponentially along the thickness direction. In the absence of body forces, the problem is reduced to a system of Fredholm singular integral equations with the contact pressure and contact size as unknowns via Fourier integral transform, which is transformed into an algebraic one by the Gauss-Chebyshev quadratures and polynomials of both the first and second kinds. Then, an iterative speediest descending algorithm is proposed to numerically solve the system of algebraic equations. Both semi-analytical and finite element method, FEM solutions for the presented problem validate each other. To improve the accuracy of the numerical result of FEM, a graded FEM solution is performed to simulate the FGM mechanical characteristics. The results reveal the potential links between the contact stress/size and the indenter size, the thickness, as well as some other material properties of FGM.

• Wind and Structures
• /
• v.37 no.3
• /
• pp.179-189
• /
• 2023

Multiple studies conducted in the past evaluated the conductor response under one tornado wind field, while the performance of transmission lines under different tornado wind fields still remains unknown. Thus, the objective of this paper is to estimate the variation in the conductor's critical longitudinal and transverse reactions under different tornado wind fields, as well as providing the corresponding critical tornado configurations. The considered full-scale tornadoes are the Spencer, South Dakota, 1998, the Stockton, Kansas, 2005 and the Goshen County, Wyoming, 2009. Computational Fluid Dynamics (CFD) simulations were previously conducted to develop these wind fields. All tornadoes have been rescaled to have a common velocity matching the upper limit of the F2 Fujita scale. Eight conductor systems, each including six spans, are considered in this paper. For each conductor, parametric studies are conducted by varying the location of the three tornado wind fields relative to the tower of interest, therefore the peak reactions associated with each tornado are determined. A semi-analytical closed-form solution, previously developed and validated, is used to calculate the reactions. The study conducted in this paper can be divided into two parts: In the first part, a parametric study considering a wide range of tornado locations is conducted. In the second part, the parametric study focuses on the tornado location leading to the critical tangential velocity on the tower. Based on this extensive parametric study, a critical tornado defined as the Design Tornado and its critical locations, tornado distance R = 125 m, tornado angle 𝜃 = 15° and 30°, are recommended for design purposes.

• Advances in aircraft and spacecraft science
• /
• v.10 no.3
• /
• pp.281-302
• /
• 2023

In this paper, bending analysis of exponentially varying functionally graded (FG) plate is presented using trigonometric shear deformation theory (TSDT) considering both transverse shear and normal deformation effects. The in-plane displacement field consists of sinusoidal functions in thickness direction to include transverse shear strains and transverse displacement include the effect of transverse normal strain using the cosine function in thickness coordinate. The governing equations and boundary conditions of the theory are derived using the virtual work principle. System of governing equations, for simply supported conditions, Navier's solution technique is used to obtain results. Plate material properties vary across thickness direction according to exponential distribution law. In the current theory, transverse shear stresses are distributed accurately through the plate thickness, hence obviates the need for a shear correction factor. TSDT results are compared with those from other theories to ensure the accuracy and effectiveness of the present theory. The current theory is in excellent agreement with the semi-analytical theory.

• Computers and Concrete
• /
• v.33 no.1
• /
• pp.13-25
• /
• 2024

In this study, contact analysis in a functionally graded (FG) layer loaded with two circular punches is solved using the finite element method (FEM). The problem is consisted of a functionally graded layer that resting on an elastic semi-infinite plane and is loaded with two rigid punches of circular geometry. External loads P and Q are transferred to the layer via two rigid punches. The finite element model of the functionally graded layer is created using the ANSYS package program and a 2-dimensional analysis of the problem is analyzed. The contact lengths, obtained as a result of the analysis are compared with the analytical solution in the literature. In the study, the effects of parameters such as distances between punches, loads, inhomogenity parameter on contact zones, initial separation loads and distances, normal stresses, stresses across depth and contact stresses are investigated. As a result, in this study, it can be said that the magnitude of the stresses occurring in the FG layer is less than the homogeneous layer, therefore the life of FG materials will be longer than the homogeneous layer. When the distance between the punches is 2.25, the initial separation distance is 6.98, and when the distance between the punches is 4, the initial separation distance decreases to 6.10. In addition, when the load increased in the second punch, the initial separation load decreased from 55 to 18. The obtained results are presented in the form of graphs and tables.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.1
• /
• pp.117-124
• /
• 2009

The p-convergent boundary element method has been proposed to analyze two-dimensional potential problem on the basis of high order Legendre shape functions that have different property comparing with the shape functions in conventional boundary element method. The location of nodes corresponding to high order shape function are not defined along the boundary, called by nodeless node, similar to the p-convergent finite element method. As the order of shape function increases, the collocation point method is used to solve linear simultaneous equations. The collocation patterns of p-convergent boundary element method consist of non-symmetric hierarchial or symmetric non-hierarchical. As the order of shape function increases, the number of collocation point increases. The singular integral that appears in p-convergent boundary element has been calculated by special numeric quadrature technique and semi-analytical integration technique. The L-shape domain problem including singularity in the vicinity of reentrant comer is analyzed and the numerical results show that the relative error is smaller than $10^{-2}%$ range as compared with other results in literatures. In case of same condition, the symmetric p-collocation point pattern shows high accuracy of solution.

• Journal of the Institute of Electronics and Information Engineers
• /
• v.49 no.9
• /
• pp.3-10
• /
• 2012

In this paper, CDD(Cyclic Delay Diversity)-DSF(decode-and-selectively forward) relay schemes are proposed for MC-CDMA(multicarrier-code division multiple access) systems over multipath Rayleigh fading channels. On the contrary to general DSF schemes, each relay equipped a CDD method can transmit on single channel. Therefore, all link channels between destination and relays can be considered as a single channel which is widely delay spread. The cooperative diversity of relay networks can be obtained as the frequency diversity for MC-CDMA systems with proposed CDD schemes. Furthermore, a semi-analytical approach for the averaged bit error rare(BER) performance evaluation is proposed for the proposed CDD-DSF-Relay schemes. By simulation and numerical results, it is confirmed that the proposed one can be a possible solution to achieve cooperative diversity gain without a reduction of spectral efficiency.

• Smart Structures and Systems
• /
• v.14 no.6
• /
• pp.1173-1196
• /
• 2014

Real-time dynamic substructuring tests have been conducted on a cable-deck system. The cable is representative of a full scale cable for a cable-stayed bridge and it interacts with a deck, numerically modelled as a single-degree-of-freedom system. The purpose of exciting the inclined cable at the bottom is to identify its nonlinear dynamics and to mark the stability boundary of the semi-trivial solution. The latter physically corresponds to the point at which the cable starts to have an out-of-plane response when both input and previous response were in-plane. The numerical and the physical parts of the system interact through a transfer system, which is an actuator, and the input signal generated by the numerical model is assumed to interact instantaneously with the system. However, only an ideal system manifests a perfect correspondence between the desired signal and the applied signal. In fact, the transfer system introduces into the desired input signal a delay, which considerably affects the feedback force that, in turn, is processed to generate a new input. The effectiveness of the control algorithm is measured by using the synchronization technique, while the online adaptive forward prediction algorithm is used to compensate for the delay error, which is present in the performed tests. The response of the cable interacting with the deck has been experimentally observed, both in the presence of delay and when delay is compensated for, and it has been compared with the analytical model. The effects of the interface delay in real-time dynamic substructuring tests conducted on the cable-deck system are extensively discussed.

• 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.

• Advances in aircraft and spacecraft science
• /
• v.8 no.4
• /
• pp.345-372
• /
• 2021

The analysis of nonlinear vibrations, buckling, post-buckling, flutter boundary determination and post-flutter behavior of a homogeneous curved plate assuming cylindrical bending is conducted in this article. Other assumptions include simply-supported boundary conditions, supersonic aerodynamic flow at the top of the plate, constant pressure conditions below the plate, non-viscous flow model (using first- and third-order piston theory), nonlinear structural model with large deformations, and application of mechanical and thermal loads on the curved plate. The analysis is performed with constant environmental indicators (flow density, heat, Reynolds number and Mach number). The material properties (i.e., coefficient of thermal expansion and modulus of elasticity) are temperature-dependent. The equations are derived using the principle of virtual displacement. Furthermore, based on the definitions of virtual work, the potential and kinetic energy of the final relations in the integral form, and the governing nonlinear differential equations are obtained after fractional integration. This problem is solved using two approaches. The frequency analysis and flutter are studied in the first approach by transferring the handle of ordinary differential equations to the state space, calculating the system Jacobin matrix and analyzing the eigenvalue to determine the instability conditions. The second approach discusses the nonlinear frequency analysis and nonlinear flutter using the semi-analytical solution of governing differential equations based on the weighted residual method. The partial differential equations are converted to ordinary differential equations, after which they are solved based on the Runge-Kutta fourth- and fifth-order methods. The comparison between the results of frequency and flutter analysis of curved plate is linearly and nonlinearly performed for the first time. The results show that the plate curvature has a profound impact on the instability boundary of the plate under supersonic aerodynamic loading. The flutter boundary decreases with growing thermal load and increases with growing curvature.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.21 no.3
• /
• pp.503-509
• /
• 2020

In offshore floating structures, the viscous mean drift force due to drag is considered a design part that has not been considered until recently. In this paper, an analytical solution for the viscous mean drift forces on a floating vertical cylinder considering the waves and currents was obtained. The area was considered by dividing it into a splash zone above the free surface and a submerged zone below the free surface. In the case of waves, only the splash zone was considered; in the case of waves and currents, equations were obtained in both the splash zone and the submerged zone. The RAO results of previous studies were used to compare the calculated results with the drift forces acting on the fixed cylinder. Except for the case in only waves in the splash zone, the viscous mean drift force acting on the floating cylinder was larger than the drift force acting on the relatively fixed cylinder in most frequencies. In particular, the increase was greater when the currents were considered to be more important. Therefore, these results provide the inference for the viscous drift force due to drag in the design of floating offshore structures.