• Smart Structures and Systems
• /
• v.20 no.4
• /
• pp.509-524
• /
• 2017

In this article, a free vibration analysis of functionally graded (FG) plates resting on elastic foundations is presented using a quasi-3D hybrid-type higher order shear deformation theory. Undetermined integral terms are employed in the proposed displacement field and modeled based on a hybrid-type (sinusoidal and parabolic) quasi-3D HSDT with five unknowns in which the stretching effect is taken into account. Thus, it can be said that the significant feature of this theory is that it deals with only 5 unknowns as the first order shear deformation theory (FSDT). The elastic foundation parameters are introduced in the present formulation by following the Pasternak (two-parameter) mathematical model. Equations of motion are obtained via the Hamilton's principles and solved using Navier's method. Accuracy of the proposed theory is confirmed by comparing the results of numerical examples with the ones available in literature.

• Steel and Composite Structures
• /
• v.27 no.1
• /
• pp.13-25
• /
• 2018

The main objective of this study is to develop a dual approach for geometrically nonlinear finite element analysis of plane truss structures. The geometric nonlinearity is considered using the Total Lagrangian formulation. The nonlinear solution is obtained by introducing and minimizing an objective function subjected to displacement-type constraints. The proposed method can fully trace the whole equilibrium path of geometrically nonlinear plane truss structures not only before the limit point but also after it. No stiffness matrix is used in the main approach and the solution is acquired only based on the direct classical stress-strain formulations. As a result, produced errors caused by linearization and approximation of the main equilibrium equation will be eliminated. The suggested algorithm can predict both pre- and post-buckling behavior of the steel plane truss structures as well as any arbitrary point of equilibrium path. In addition, an equilibrium path with multiple limit points and snap-back phenomenon can be followed in this approach. To demonstrate the accuracy, efficiency and robustness of the proposed procedure, numerical results of the suggested approach are compared with theoretical solution, modified arc-length method, and those of reported in the literature.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.10 no.5
• /
• pp.583-595
• /
• 2018

Stiffened panels are widely used in naval architecture and ocean engineering, and knowledge about their dynamic behaviour represents important issue in the design procedure. Ordinary vibration analysis consists of natural frequencies and mode shapes determination and can be extended to forced response assessment, while the Structural Intensity (SI) analysis, assessing magnitude and direction of vibrational energy flow provides information on dominant transmission paths and energy distribution including sink positions. In this paper, vibrational energy flow in stiffened panels under harmonic loading is analyzed by the SI technique employing the finite element method. Structural intensity formulation for plate and beam element is outlined, and developed system combining in-house code and general finite element tool is described. As confirmed within numerical examples, the developed tool enables separation of SI components, enabling generation of novel SI patterns and providing deeper insight in the vibrational energy flow in stiffened panels, comparing to existing works.

• Advances in aircraft and spacecraft science
• /
• v.3 no.4
• /
• pp.471-502
• /
• 2016

This paper presents the dynamic instability analysis of un-damped elastically supported piezoelectric functionally graded (FG) beams subjected to in-plane static and dynamic periodic thermomechanical loadings with uncertain system properties. The elastic foundation model is assumed as one parameter Pasternak foundation with Winkler cubic nonlinearity. The piezoelectric FG beam is subjected to non-uniform temperature distribution with temperature dependent material properties. The Young's modulus and Poison's ratio of ceramic, metal and piezoelectric, density of respective ceramic and metal, volume fraction exponent and foundation parameters are taken as uncertain system properties. The basic nonlinear formulation of the beam is based on higher order shear deformation theory (HSDT) with von-Karman strain kinematics. The governing deterministic static and dynamic random instability equation and regions is solved by Bolotin's approach with Newmark's time integration method combined with first order perturbation technique (FOPT). Typical numerical results in terms of the mean and standard deviation of dynamic instability analysis are presented to examine the effect of slenderness ratios, volume fraction exponents, foundation parameters, amplitude ratios, temperature increments and position of piezoelectric layers by changing the random system properties. The correctness of the present stochastic model is examined by comparing the results with direct Monte Caro simulation (MCS).

• Journal of the Institute of Electronics Engineers of Korea TC
• /
• v.39 no.4
• /
• pp.201-208
• /
• 2002

In this paper, we present a new formulation using time-domain electric field integral equation (TD-EFIE) to obtain transient scattering response from arbitrarily shaped three-dimensional conducting bodies. The time derivative of the magnetic vector potential is approximated with a central finite difference and the scalar potential is time averaged by dividing it into two terms. This approach with an implicit method using central difference results in accurate and more stable transient scattering responses from conducting objects. Detailed mathematical steps are included and several numerical results are presented and compared with the inverse discrete Fourier transform (IDFT) of the frequency-domain solution.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.14 no.5
• /
• pp.458-465
• /
• 2003

In this paper, we analyze the transient electromagnetic response from three-dimensional(3-D) dielectric bodies using a time-domain electric field integral equation formulation. The solution method in this paper is based on the Galerkin＇s method that involves separate spatial and temporal testing procedures. Triangular patch basis functions are used for spatial expansion and testing functions for arbitrarily shaped 3-D dielectric structures. The time-domain unknown coefficients of the equivalent electric and magnetic currents are approximated as an orthonormal basis function set that is derived from the Laguerre functions. These basis functions are also used as the temporal testing. Numerical results involving equivalent currents and far fields computed by the proposed method are presented.

• Proceedings of the Korean Society of Computational and Applied Mathematics Conference
• /
• 2003.09a
• /
• pp.6-6
• /
• 2003

In this work we consider the mathematical formulation and numerical resolution of the linear feedback control problem for Boussinesq equations. The controlled Boussinesq equations is given by $$\frac{{\partial}u}{{\partial}t}-{\nu}{\Delta}u+(u{\cdot}{\nabla}u+{\nabla}p={\beta}{\theta}g+f+F\;\;in\;(0,\;T){\times}\;{\Omega}$$, $${\nabla}{\cdot}u=0\;\;in\;(0,\;T){\times}{\Omega}$$, $$u|_{{\partial}{\Omega}=0,\;u(0,x)=\;u_0(x)$$ $$\frac{{\partial}{\theta}}{{\partial}t}-k{\Delta}{\theta}+(u{\cdot}){\theta}={\tau}+T,\;\;in(0,\;T){\times}{\Omega}$$ $${\theta}|_{{\partial}{\Omega}=0,\;\;{\theta}(0,X)={\theta}_0(X)$$, where $\Omega$ is a bounded open set in $R^{n}$, n=2 or 3 with a $C^{\infty}$ boundary ${\partial}{\Omega}$. The control is achieved by means of a linear feedback law relating the body forces to the velocity and temperature field, i.e., $$f=-{\gamma}_1(u-U),\;\;{\tau}=-{\gamma}_2({\theta}-{\Theta}}$$ where (U,$\Theta$) are target velocity and temperature. We show that the unsteady solutions to Boussinesq equations are stabilizable by internal controllers with exponential decaying property. In order to compute (approximations to) solution, semi discrete-in-time and full space-time discrete approximations are also studied. We prove that the difference between the solution of the discrete problem and the target solution decay to zero exponentially for sufficiently small time step.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.6
• /
• pp.987-992
• /
• 2002

In an attempt to solve multiobjective optimization problems, many traditional methods scalarize the objective vector into a single objective. In those cases, the obtained solution is highly sensitive to the weight vector used in the scalarization process and demands the user to have knowledge about the underlying problem. Moreover, in solving multiobjective problems, designers may be interested in a set of Pareto-optimal points, instead of a single point. In this paper, pareto-based Continuous Evolutionary Algorithms for Multiobjective Optimization problems having continuous search space are introduced. This algorithm is based on Continuous Evolutionary Algorithms to solve single objective optimization problems with a continuous function and continuous search space efficiently. For multiobjective optimization, a progressive reproduction operator and a niche-formation method fur fitness sharing and a storing process for elitism are implemented in the algorithm. The operator and the niche formulation allow the solution set to be distributed widely over the Pareto-optimal tradeoff surface. Finally, the validity of this method has been demonstrated through a numerical example.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.10
• /
• pp.1971-1978
• /
• 1992

A turbulent flow in a conical diffuser with total divergence angle of 8.deg. was numerically studied. The low Reynolds number k-.epsilon. model(Launder-Sharma model) was adopted to simulate the turbulence. The continuity and time averaged Navier-Stokes equations in a nonorthogonal coordinate system were solved by a finite volume method based on the fully elliptic formulation. The low Reynolds number k-.epsilon. model reasonably simulates the pressure recovery and the mean velocity components. However, there are also considerable discrepancies between predicted and measured shear stress distribution on the wall and turbulent kinetic energy distributions. It is necessary to investigate the flow structure at the entry of the diffuser, numerically as well as experimentally.

• Journal of the Korean Society of Food Culture
• /
• v.27 no.6
• /
• pp.710-718
• /
• 2012

The purpose of this study was to investigate the antioxidative effects of Opuntia humifusa and determine the optimal mixing ratio of Opuntia humifusa powder and eggs for the preparation of pasta. In regard to its antioxidant effects, Opuntia humifusa powder had a total phenolic content and DPPH free radical scavenging activity of 6.64 mg/g and 90.12%, respectively. We used response surface methodology to obtain ten experimental points (including two replicates for Opuntia humifusa powder and egg) and Opuntia humifusa pasta formulation was optimized using rheology. Yellowness (p<0.001) displayed a linear model pattern, whereas lightness (p<0.01) and redness (p<0.01) were described by a quadratic model pattern. Mechanical properties showed significant values in hardness (p<0.5), chewiness (p<0.5), gummiess (p<0.5), and cohesiveness (p<0.5). The sensory evaluation parameters showed significant differences in color (p<0.5), appearance (p<0.5), flavor (p<0.5), texture (p<0.5) and overall quality (p<0.5). The optimum formulations processed by numerical and graphical optimization were found to be 7.30 g of Opuntia humifusa powder and 31.55 g of egg.