A Complete Parametric Study of Pull-In Voltage by Nonlinear Differential Equation

2014 ◽  
Author(s):  
M. Amin Changizi ◽  
Ion Stiharu
Keyword(s):  
Differential Equation ◽  
Electric Field ◽  
Cantilever Beam ◽  
Parametric Study ◽  
Nonlinear Ordinary Differential Equation ◽  
Cantilever Beams ◽  
Time Step ◽  
Robust Adaptive ◽  
Linear Differential ◽  
Micro Cantilever

Micro-cantilever beams are interested structures in MEMS because of their fabrication is very easy and its versatility. The importance of micro-cantilevers beam in MEMS has driven various investigations like static and dynamic performances under different loading such as potential fields. In this research the non-linear differential equation which models dynamics of a micro-cantilever beams vibration subjected to electrostatic field has been studied. The model which has one degree of freedom is used to calculate the pull-in voltage. This model adopted based on different method of calculating stiffness of micro-cantilever beam. The nonlinear ordinary differential equation which used to model the dynamics of the cantilever subjected to electric field close to snap on is highly stiff. Investigation on solving of nonlinear stiff ordinary equation showed that only Lsode algorithm yield to correct solution to the problem. Lsode is equipped with a robust adaptive time step selection mechanism that enables solutions to very stiff problems, as the one under discussion. The best match in the resonant frequency for equivalent stiffness based on four different models was considered. The stiffness model suitable for the best match in deflection is proved to be different from the model that yields. Pull-in voltage under electric field was studied. Pull-in voltage has been investigated from the analytical and numerical perspective. A complete parametric study of structural damping effect on large deflection of micro-cantilever beam was studied was done numerically in this work. Different kind of impulse voltages were considered and effect of them on pulling voltage numerically was studied. A cumbersome mathematical method, Lie symmetry, was used to drive a closed from of time response to step voltage for undamped system and pull in voltage of such system was calculated. Finally, a closed form driven from the nonlinear ODE for calculating pulling voltage was presented.

High Precision Numerical Analysis of Nonlinear Beam Bending Problems under Large Deflection

2014 ◽  
Vol 638-640 ◽  
pp. 1705-1709
Author(s):  
Zhao Qing Wang ◽  
Jian Jiang ◽  
Bing Tao Tang ◽  
Wei Zheng
Keyword(s):  
Differential Equation ◽  
High Precision ◽  
Large Deflection ◽  
Nonlinear Ordinary Differential Equation ◽  
Initial Guess ◽  
Computational Accuracy ◽  
Nonlinear Bending ◽  
Barycentric Interpolation ◽  
Bending Problems ◽  
Linear Differential

The high precision numerical method for solving nonlinear bending problems of large deflection beam is presented. The governing equation of large deflection beam is a strongly nonlinear ordinary differential equation. Using the solution of linear bending beam as an initial guess function, the nonlinear bending equation of beam can be transferred into a linear differential equation. The improvement solution of nonlinear bending beam is obtained by solving the linearized bending equation using barycentric interpolation collocation method. Then, the solution of nonlinear bending beam can be given by iterative method. Some examples demonstrate the validity and computational accuracy of proposed method.

MEMS Wind Speed Sensor: Large Deflection of Curved Micro-Cantilever Beam Under Uniform Horizontal Force

2015 ◽  
Author(s):  
M. Amin Changizi ◽  
Ali Abolfathi ◽  
Ion Stiharu
Keyword(s):  
Wind Speed ◽  
Numerical Solution ◽  
Cantilever Beam ◽  
Large Deflection ◽  
Total Force ◽  
Curved Beam ◽  
Cantilever Beams ◽  
Canonical Coordinates ◽  
Nonlinear Ode ◽  
Micro Cantilever

Micro-cantilever beams are currently employed as sensor in various fields. Of main applications, is using such beams in wind speed sensors. For this purpose, curved out of plane micro-cantilever beams are used. Uniform pressure on such beams causes a large deflection of beam. General mechanics of material theory deals with small deflection and thus cannot be used for explaining this deflection. Although there are a body of works on analysing of large deflection [1], nonlinear deflection, of curved beams [2], yet there is no research on large deflection of curved beam under horizontal uniform distributed force. Theoretically, the wind force is applying horizontally on curved micro-cantilever beam. Here, we neglect the effect of moving weather from beam sides. We first aim how to drive the governed equation. A curved beam does not have a calculable centroid. Also large deflection of beam changes its curvature which would change the centroid of beam consciously. The variation of centroid makes very though calculating the bending moment of each cross section in the beam. To address this issue, an integral equation will be used. The total force will be considered as a single force applied at the centroid. The second challenge is solving the governed nonlinear ordinary differential equation (ODE). Although there are several methods to solve analytically nonlinear ODE, Lie symmetry method, with all its complication, is a general method for this kind of equations. This approach covers all current methods in analytical solving nonlinear ODEs. In this method, an infinitesimal transformation should be calculated. All transformations under one parameter creates a group that called Lie group. A value of parameter which transfers the equation onto itself is called invariant of ODE. One can calculate canonical coordinates ODEs by the invariant. Solving the canonical coordinates ODEs yields to calculating the canonical coordinates. Canonical coordinate are used to reduce the order of nonlinear ODE [3]. By repeating this method one can solve high order ODEs. Our last question is how to do numerical solution of ODE. The possible answer will help to explain the phenomena of deflection clearly and compare the analytical solution with numerical results. Small dimensions of beam, small values of applied force from one side and Young modules value from the other side, will create a stiff ODE. Authors experience in this area shows that the best method to sole these kind of equations is LSODE. This method can be used in Maple. Here, primary calculations show that the governed equation is second order nonlinear ODE and we propose two possible invariants to solve ODE. Overall, the primary numerical solution has shown perfect match with the exact solution.

scholarly journals On the Numerical Simulation of HPDEs Using θ-Weighted Scheme and the Galerkin Method

Mathematics ◽  
2020 ◽  
Vol 9 (1) ◽  
pp. 78
Author(s):  
Haifa Bin Jebreen ◽  
Fairouz Tchier
Keyword(s):  
Differential Equation ◽  
Galerkin Method ◽  
Linear Ordinary Differential Equation ◽  
Time Interval ◽  
Hyperbolic Partial Differential Equation ◽  
Time Step ◽  
Numerical Examples ◽  
One Dimensional ◽  
The Stability ◽  
The Galerkin Method

Herein, an efficient algorithm is proposed to solve a one-dimensional hyperbolic partial differential equation. To reach an approximate solution, we employ the θ-weighted scheme to discretize the time interval into a finite number of time steps. In each step, we have a linear ordinary differential equation. Applying the Galerkin method based on interpolating scaling functions, we can solve this ODE. Therefore, in each time step, the solution can be found as a continuous function. Stability, consistency, and convergence of the proposed method are investigated. Several numerical examples are devoted to show the accuracy and efficiency of the method and guarantee the validity of the stability, consistency, and convergence analysis.

scholarly journals Correctness conditions for high-order differential equations with unbounded coefficients

2021 ◽  
Vol 2021 (1) ◽  
Author(s):  
Kordan N. Ospanov
Keyword(s):  
Differential Equation ◽  
Hilbert Space ◽  
Existence And Uniqueness ◽  
Sufficient Conditions ◽  
Integral Operators ◽  
Order Linear Differential Equation ◽  
Unbounded Coefficients ◽  
Weighted Norms ◽  
Uniqueness Of The Solution ◽  
Linear Differential

AbstractWe give some sufficient conditions for the existence and uniqueness of the solution of a higher-order linear differential equation with unbounded coefficients in the Hilbert space. We obtain some estimates for the weighted norms of the solution and its derivatives. Using these estimates, we show the conditions for the compactness of some integral operators associated with the resolvent.

Asymptotic behavior of solutions of half-linear differential equations and generalized Karamata functions

2020 ◽  
Vol 0 (0) ◽  
Author(s):  
Kusano Takaŝi ◽  
Jelena V. Manojlović
Keyword(s):  
Differential Equation ◽  
Asymptotic Behavior ◽  
Continuous Function ◽  
Differential Equations ◽  
Linear Differential Equation ◽  
Regular Variation ◽  
Asymptotic Formulas ◽  
Asymptotic Behavior Of Solutions ◽  
Positive Continuous Function ◽  
Linear Differential

AbstractWe study the asymptotic behavior of eventually positive solutions of the second-order half-linear differential equation(p(t)\lvert x^{\prime}\rvert^{\alpha}\operatorname{sgn}x^{\prime})^{\prime}+q(% t)\lvert x\rvert^{\alpha}\operatorname{sgn}x=0,where q is a continuous function which may take both positive and negative values in any neighborhood of infinity and p is a positive continuous function satisfying one of the conditions\int_{a}^{\infty}\frac{ds}{p(s)^{1/\alpha}}=\infty\quad\text{or}\quad\int_{a}^% {\infty}\frac{ds}{p(s)^{1/\alpha}}<\infty.The asymptotic formulas for generalized regularly varying solutions are established using the Karamata theory of regular variation.

scholarly journals Hyperbolastic Models from a Stochastic Differential Equation Point of View

Mathematics ◽  
2021 ◽  
Vol 9 (16) ◽  
pp. 1835
Author(s):  
Antonio Barrera ◽  
Patricia Román-Román ◽  
Francisco Torres-Ruiz
Keyword(s):  
Differential Equation ◽  
Simulated Data ◽  
Real Data ◽  
Point Of View ◽  
Likelihood Method ◽  
Numerical Resolution ◽  
The Family ◽  
The Common ◽  
Linear Differential ◽  
Mean Function

A joint and unified vision of stochastic diffusion models associated with the family of hyperbolastic curves is presented. The motivation behind this approach stems from the fact that all hyperbolastic curves verify a linear differential equation of the Malthusian type. By virtue of this, and by adding a multiplicative noise to said ordinary differential equation, a diffusion process may be associated with each curve whose mean function is said curve. The inference in the resulting processes is presented jointly, as well as the strategies developed to obtain the initial solutions necessary for the numerical resolution of the system of equations resulting from the application of the maximum likelihood method. The common perspective presented is especially useful for the implementation of the necessary procedures for fitting the models to real data. Some examples based on simulated data support the suitability of the development described in the present paper.

scholarly journals Frequency-Dependent FDTD Algorithm Using Newmark’s Method

2014 ◽  
Vol 2014 ◽  
pp. 1-6 ◽  
Author(s):  
Bing Wei ◽  
Le Cao ◽  
Fei Wang ◽  
Qian Yang
Keyword(s):  
Differential Equation ◽  
Electric Field ◽  
Field Intensity ◽  
Time Domain ◽  
Difference Method ◽  
Electric Field Intensity ◽  
Polarization Vector ◽  
Solution Stability ◽  
Central Difference ◽  
Central Difference Method

According to the characteristics of the polarizability in frequency domain of three common models of dispersive media, the relation between the polarization vector and electric field intensity is converted into a time domain differential equation of second order with the polarization vector by using the conversion from frequency to time domain. Newmarkβγdifference method is employed to solve this equation. The electric field intensity to polarizability recursion is derived, and the electric flux to electric field intensity recursion is obtained by constitutive relation. Then FDTD iterative computation in time domain of electric and magnetic field components in dispersive medium is completed. By analyzing the solution stability of the above differential equation using central difference method, it is proved that this method has more advantages in the selection of time step. Theoretical analyses and numerical results demonstrate that this method is a general algorithm and it has advantages of higher accuracy and stability over the algorithms based on central difference method.

scholarly journals Steady Motion of Ice Sheets

1980 ◽  
Vol 25 (92) ◽  
pp. 229-246 ◽  
Author(s):  
L. W. Morland ◽  
I. R. Johnson
Keyword(s):  
Differential Equation ◽  
Linear Differential Equation ◽  
Power Law ◽  
Perturbation Expansion ◽  
Ice Sheet ◽  
Leading Order ◽  
Plane Flow ◽  
General Law ◽  
Non Linear ◽  
Linear Differential

AbstractSteady plane flow under gravity of a symmetric ice sheet resting on a horizontal rigid bed, subject to surface accumulation and ablation, basal drainage, and basal sliding according to a shear-traction-velocity power law, is treated. The surface accumulation is taken to depend on height, and the drainage and sliding coefficient also depend on the height of overlying ice. The ice is described as a general non-linearly viscous incompressible fluid, with illustrations presented for Glen’s power law, the polynomial law of Colbeck and Evans, and a Newtonian fluid. Uniform temperature is assumed so that effects of a realistic temperature distribution on the ice response are not taken into account. In dimensionless variables a small paramter ν occurs, but the ν = 0 solution corresponds to an unbounded sheet of uniform depth. To obtain a bounded sheet, a horizontal coordinate scaling by a small factor ε(ν) is required, so that the aspect ratio ε of a steady ice sheet is determined by the ice properties, accumulation magnitude, and the magnitude of the central thickness. A perturbation expansion in ε gives simple leading-order terms for the stress and velocity components, and generates a first order non-linear differential equation for the free-surface slope, which is then integrated to determine the profile. The non-linear differential equation can be solved explicitly for a linear sliding law in the Newtonian case. For the general law it is shown that the leading-order approximation is valid both at the margin and in the central zone provided that the power and coefficient in the sliding law satisfy certain restrictions.

Sign in / Sign up

Export Citation Format

Share Document