• Structural Engineering and Mechanics
• /
• v.25 no.5
• /
• pp.501-518
• /
• 2007

In this work, a new approach is developed for dynamic analysis of a composite beam with an interply crack, based on finite element solution of partial differential equations with the use of the COMSOL Multiphysics package, allowing for fast and simple change of geometric characteristics of the delaminated area. The use of COMSOL Multiphysics package facilitates automatic mesh generation, which is needed if the problem has to be solved many times with different crack lengths. In the model, a physically impossible interpenetration of the crack faces is prevented by imposing a special constraint, leading to taking account of a force of contact interaction of the crack faces and to nonlinearity of the formulated boundary value problem. The model is based on the first order shear deformation theory, i.e., the longitudinal displacement is assumed to vary linearly through the beam's thickness. The shear deformation and rotary inertia terms are included into the formulation, to achieve better accuracy. Nonlinear partial differential equations of motion with boundary conditions are developed and written in the format acceptable by the COMSOL Multiphysics package. An example problem of a clamped-free beam with a piezoelectric actuator is considered, and its finite element solution is obtained. A noticeable difference of forced vibrations of the delaminated and undelaminated beams due to the contact interaction of the crack's faces is predicted by the developed model.

• Journal of the Korean Mathematical Society
• /
• v.50 no.5
• /
• pp.1129-1163
• /
• 2013

Based on finite element discretization, two linearization approaches to the defect-correction method for the steady incompressible Navier-Stokes equations are discussed and investigated. By applying $m$ times of Newton and Picard iterations to solve an artificial viscosity stabilized nonlinear Navier-Stokes problem, respectively, and then correcting the solution by solving a linear problem, two linearized defect-correction algorithms are proposed and analyzed. Error estimates with respect to the mesh size $h$, the kinematic viscosity ${\nu}$, the stability factor ${\alpha}$ and the number of nonlinear iterations $m$ for the discrete solution are derived for the linearized one-step defect-correction algorithms. Efficient stopping criteria for the nonlinear iterations are derived. The influence of the linearizations on the accuracy of the approximate solutions are also investigated. Finally, numerical experiments on a problem with known analytical solution, the lid-driven cavity flow, and the flow over a backward-facing step are performed to verify the theoretical results and demonstrate the effectiveness of the proposed defect-correction algorithms.

• Advances in nano research
• /
• v.9 no.3
• /
• pp.133-145
• /
• 2020

In this study, nonlinear vibrations and dynamic instabilities of a smart embedded micro shell conveying varied fluid flow and subjected to the combined electro-thermo-mechanical loadings are investigated. With the aim of designing new hydraulic sensors and actuators, the piezoelectric materials are employed for the body and the effects of applying electric field on the stability of the system as well as the induced voltage due to the dynamic behavior of the system are studied. The nonlocal piezoelasticity theory and the nonlinear cylindrical shell model in conjunction with the energy approach are utilized to mathematically modeling of the structure. The fluid flow is assumed to be isentropic, incompressible and fully develop, and for more generality of the problem both steady and time dependent flow regimes are considered. The mathematical modeling of fluid flow is also carried out based on a scalar potential function, time mean Navier-Stokes equations and the theory of slip boundary condition. Employing the modified Lagrange equations for open systems, the nonlinear coupled governing equations of motion are achieved and solved via the state space problem; forth order numerical integration and Bolotin's method. In the numerical results, a comprehensive discussion is made on the dynamical instabilities of the system (such as divergence, flutter and parametric resonance). We found that applying positive electric potential field will improve the stability of the system as an actuator or vibration amplitude controller in the micro electro mechanical systems.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.4
• /
• pp.427-436
• /
• 2000

This study presents the nonlinear responses of a hinged-clamped beam under broadband random excitation. By using Galerkin's method the governing equation is reduced to a system or nonautonomous nonlinear ordinary differential equations. The Fokker-Planck equation is used to generate a general first-order differential equation in the joint moments of response coordinates. Gaussian and non-Gaussian closure schemes are used to close the infinite coupled moment equations. The closed equations are then solved for response statistics in terms of system and excitation parameters. The case of two mode interaction is considered in order to compare it with the case of three mode interaction. Monte Carlo simulation is used for numerical verification.

• Journal of Advanced Research in Ocean Engineering
• /
• v.1 no.3
• /
• pp.157-167
• /
• 2015

In order to provide simple and accurate wave theory in design of offshore structure, an analytical approximation is introduced in this paper. The solution is limited to flat bottom having a constant water depth. Water is considered as inviscid, incompressible and irrotational. The solution satisfies the continuity equation, bottom boundary condition and non-linear kinematic free surface boundary condition exactly. Error for dynamic condition is quite small. The solution is suitable in description of breaking waves. The solution is presented with closed form and dispersion relation is also presented with closed form. In the last century, there have been two main approaches to the nonlinear problems. One of these is perturbation method. Stokes wave and Cnoidal wave are based on the method. The other is numerical method. Dean's stream function theory is based on the method. In this paper, power series method was considered. The power series method can be applied to certain nonlinear differential equations (initial value problems). The series coefficients are specified by a nonlinear recurrence inherited from the differential equation. Because the non-linear wave problem is a boundary value problem, the power series method cannot be applied to the problem in general. But finite number of coefficients is necessary to describe the wave profile, truncated power series is enough. Therefore the power series method can be applied to the problem. In this case, the series coefficients are specified by a set of equations instead of recurrence. By using the set of equations, the nonlinear wave problem has been solved in this paper.

• Structural Engineering and Mechanics
• /
• v.88 no.2
• /
• pp.159-168
• /
• 2023

Cost-effective high precision hybrid elements are presented in a hierarchical form for dynamic analysis of plates. The costs associated with controlling the vibrations of ferromagnetic plates can be minimized by adequate determination of the amount of electric current and magnetic field. In the present study, the effect of magnetic field and electric current on nonlinear vibrations of ferromagnetic plates is investigated. The general form of Lorentz forces and Maxwell's equations have been considered for the first time to present new relationships for electromagnetic interaction forces with ferromagnetic plates. In order to derive the governing nonlinear differential equations, the theory of third-order shear deformations of three-dimensional plates has been applied along with the von Kármán large deformation strain-displacement relations. Afterward, the nonlinear equations are discretized using the Galerkin method, and the effect of various parameters is investigated. According to the results, electric current and magnetic field have different effects on the equivalent stiffness of ferromagnetic plates. As the electric current increases and the magnetic field decreases, the equivalent stiffness of the plate decreases. This is a phenomenon reported here for the first time. Furthermore, the magnetic field has a more significant effect on the steady-state deflection of the plate compared to the electric current. Increasing the magnetic field and electric current by 10-times results in a reduction of about 350% and an increase of 3.8% in the maximum steady-state deflection, respectively. Furthermore, the nonlinear frequency decreases as time passes, and these changes become more intense as the magnetic field increases.

• Journal of applied mathematics & informatics
• /
• v.31 no.1_2
• /
• pp.221-228
• /
• 2013

The method of upper and lower solutions and the generalized quasilinearization technique is developed for the existence and approximation of solutions to boundary value problems for higher order fractional differential equations of the type $^c\mathcal{D}^qu(t)+f(t,u(t))=0$, $t{\in}(0,1),q{\in}(n-1,n],n{\geq}2$ $u^{\prime}(0)=0,u^{\prime\prime}(0)=0,{\ldots},u^{n-1}(0)=0,u(1)={\xi}u({\eta})$, where ${\xi},{\eta}{\in}(0,1)$, the nonlinear function f is assumed to be continuous and $^c\mathcal{D}^q$ is the fractional derivative in the sense of Caputo. Existence of solution is established via the upper and lower solutions method and approximation of solutions uses the generalized quasilinearization technique.

• Journal of the Chungcheong Mathematical Society
• /
• v.21 no.3
• /
• pp.339-345
• /
• 2008

We prove the existence of the positive solution for the nonlinear system of suspension bridge equations with Dirichlet boundary condition and periodic condition $$\{u_{tt}+u_{xxxx}+av^+=1+{\epsilon}_1h_1(x,t)\text{ in }(-\frac{\pi}{2},\frac{\pi}{2}){\times}R,\\v_{tt}+v_{xxxx}+bu^+=1+{\epsilon}_2h_2(x,t)\text{ in }(-\frac{\pi}{2},\frac{\pi}{2}){\times}R,$$ where $u^+={\max}\{u,0\},\;{\epsilon}_1,\;{\epsilon}_2$ are small numbers and $h_1(x,t)$, $h_2(x,t)$ are bounded, ${\pi}$-periodic in t and even in x and t and ${\parallel}h_1{\parallel}={\parallel}h_2{\parallel}=1$.

• International Journal of Automotive Technology
• /
• v.8 no.3
• /
• pp.289-298
• /
• 2007

A nonlinear dynamic simulation model from cranking to idle speed is developed to optimize the cold start process of a diesel engine. Physically-based first order nonlinear differential equations and some algebraic equations describing engine dynamics and starter motor dynamics are used to model the performance of cold starting process which is very complex and involves many components including the cold start aiding method. These equations are solved using numerical schemes to describe the starting process of a diesel engine and to study the effects of cold starting parameters. The validity of this model is examined by a cold start test at $-20^{\circ}C$. Using the developed model the effects of the important starting variables on the cold starting processes were investigated. This model can be served as a tool for designing computer aided control systems that improve cold start performance.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.18 no.3
• /
• pp.245-256
• /
• 2014

This work is focused on the boundary layer and heat transfer characteristics of hydromagnetic flow over a stretching sheet with variable thickness. Steady, two dimensional, nonlinear, laminar flow of an incompressible, viscous and electrically conducting fluid over a stretching sheet with variable thickness and power law velocity in the presence of variable magnetic field and variable temperature is considered. Governing equations of the problem are converted into ordinary differential equations utilizing similarity transformations. The resulting non-linear differential equations are solved numerically by utilizing Nachtsheim-Swigert shooting iterative scheme for satisfaction of asymptotic boundary conditions along with fourth order Runge-Kutta integration method. Numerical computations are carried out for various values of the physical parameters and the effects over the velocity and temperature are analyzed. Numerical values of dimensionless skin friction coefficient and non-dimensional rate of heat transfer are also obtained.