Hysteresis model for magnetic materials using the Jiles-Atherton model

Size: px
Start display at page:

Download "Hysteresis model for magnetic materials using the Jiles-Atherton model"

Transcription

1 Hysteresis odel for agnetic aterials using the Jiles-Atherton odel Predrag Petrovic Technical faculty Svetog Save Cacak, Nebojsa itrovic Technical faculty Svetog Save Cacak, ilorad Stevovic Technical faculty Svetog Save Cacak, Predrag Pejovic Electrotechnical faculty Bulevar Revolucije Belgrade, Abstract This paper presents a hysteresis odel of agnetic aterials, which is based on the theory of ferroagnetic hysteresis. Results obtained by nuerical siulation of hysteresis curves of the ferroagnetic, hard agnetic core d isotropic aterials have been given. The atheatical odel consists of only nine paraeters. The hysteresis curve obtained is in good agreeent with the curve defined for these aterials in the PSpice progra packet. The proposed ethod for paraeter deterination enables siulation of y agnetic aterial, which has a characteristic R, Z or F-shaped hysteresis curve. 1. Introduction The PSpice progra packet for the siulation of electronic circuits is widely used in the design of electronic devices. The Jiles-Atherton odel [1] for defining the hysteresis curves of ferroagnetic aterials is built into this packet. However, ost users of this progra packet are not failiar with odeling hysteresis curves d they have probles when defining paraeters of the odel. A description of possibilities for the alysis of coponents with a ferroagnetic core is included in the instruction ual [2], but it only gives recoendations to the user on how to deterine the five paraeters related to agnetic properties of the aterial: The PSpice progra packet includes coplete data on only one agnetic soft aterial: the n-zn ferrite-3c8 (with the so-called R-round shaped hysteresis curve d for a series of coon core shapes. Stdard catalogues of agnetic aterials contain data on only 5 paraeters (we need 9. The reaining four paraeters have to be additionally deterined. Three paraeters available in catalogues c be used: r - residual agnetization, H C - coercive field d χ i - starting agnetic susceptibility. It c be concluded that the procedure for deterining these paraeters is not singular, as a syste of three equations with four unknowns need to be solved. The result obtained exclusively depends on the user skill in attaining satisfactory agreeent with the real hysteresis curve using ultiple paraeter fitting. The error that c be ade using this ethod c seriously influence validity of the PSpice alysis, especially when alyzing hard agnetic aterials. A satisfactory solution of this proble requires knowledge of the Jiles-Atherton odel, which itself deds isotropy of the alyzed aterial d R- shaped hysteresis curve. However, agnetic aterials with Z or F-shaped hysteresis curves are often used, which requires certain odification of the odel presented here [3]. Input data for the n-zn ferrite-3c8 has been used for evaluation of the nuerical siulation developed here [4], while agnetic characteristics obtained have been copared with corresponding characteristics read fro the hysteresis curve of this aterial obtained using paraeters given in the PSpice progra packet. 2. Theoretical odel The basis of the theoretical alysis of ferroagnetic aterials is founded on the fact that the effective agnetic field (H C in each individual doain is obtained as the su of the external agnetic field (H d the influence of surrounding agnetic doains: He = H + α (1 Where: - agnetization of the aterial d α ALPHA. The hysteresis curve represents the dependence of agnetization on the agnetic field in the case when the pinning effect of the doain wall is disregarded; i.e. it represents the iniu doain energy. It c be deterined using Lengevins function

2 for copletely isotropic aterials L(H C [5] in the for of: He a ( H = s L( He = s coth a He (2 Where a A. However, in a real aterial, the oveent of doain walls is prevented (grain boundaries, grain inhoogeneneities, dislocations, non-agnetic inclusions, regions of inhoogeneous strain, which is defined as the pinning effect of a doain wall. This effect is active until the agnetic potential becoes high enough to cause wall oveent, until the next pinning center is fored. oveent of the doain wall c be reversible d irreversible. Irreversible oveent of the doain wall is described by the following equation: dirr irr = (3 dh δ k α ( irr μo Where: irr - irreversible agnetization coponent, k - a coefficient which depends on the energy required for doain wall oveent, μ - agnetic pereability of vacuu, δ= ±1 (+ denotes the direction of the increase of the agnetic field d - denotes the decrease of the agnetic field d finally, paraeter K is defined as: K k/μ. If bending of the doain wall is assued, the reversible coponent of agnetization rev c be described by the equation: ( = c (4 rev irr Where c C. The total agnetization is obtained as the su of: = + (5 irr rev So, using eqs. (3, (4 d (5 the expression for differential agnetic susceptibility χ is obtained as: d ( 1 c χ = = + c d irr (6 dh K δ α ( irr dh Variables which c be read fro the hysteresis curve d starting agnetization curve ( r, H C,, H, χ i... do not appear in the given equations, so further odel developent is required for their solution. The atheatical odel for hard agnetic devices akes use of the doain coupling paraeter α of the theory of ferroagnetic hysteresis. The strong coupling of the doain walls is considered as the ain reason of the hysteresis of hard agnetic aterials. Keeping the odel as siple as possible eq. (4 is neglected. Therefore, (1 is replaced by other hysteresis using power functions. Furtherore, the feedback loop in the differential equation is drooped too. The atheatical odel consist only of one nonlinear algebraic d one linear ordinary differential equation. H / s + α = + f ( / s (7 Kc s Rμ where g ( / s, f ( / s = g( / s < (8 g ( H β H β1 / $ $ s = + 1 s s d if sign dh, = 1d dt d = if sign dh (9, = 1d μ dh dt sign dh, elsewere K dt The paraeters of the odel are K c, α, s, R μ, H, H 1, β, β 1, d K. H, β are obtained by a curve fit of the second d H 1, β 1, of the first quadrt. The additional paraeter K c is redundt. It was included in the odel to ake the paraeter deterination easier. The existing isotropic odel equations of hysteresis in agnetic aterials have been extended to include the effects of isotropy d texture. This has proved particularly iportt as the odel is been used to describe increasing rge of agnetic aterials in which isotropy plays a significt role. Anisotropy d texture in polycrystalline agnetic aterials c be adequately described by odifying the equation for the hysteretic curve to account for these effects. To include isotropy effects into the odel, the isotropy energy ust be incorporated into the total energy of the oents. Following the developent of the generalized hysteretic function described previously [6]: iso = e E/ kbt all oents s E/ kbt e all oents cosθ (1 Where θ is the gle between the direction of the agnetic oent d the direction of the applied field, d E = μ ( H+ α + E iso (11 And E iso is the isotropy energy, which depends on the isotropy structure of aterial. In the case of cubic isotropy,

3 E iso = K1 cos i cos j i j θ θ (12 With the noral convention on sybols. In the present description, we only used the first isotropy coefficient, because this approxiation is in ost case sufficient to provide accurate description of the different agnetization curves along different directions. To incorporate the texture effect into the odel, we introduce texture coefficient t, which is a statistical evaluation of the fraction of the textured portion of the aterial. Th the hysteretic agnetization c be given as: = t + ( t iso iso 1 (13 Where iso is given by (1 d is the isotropic hysteretic agnetization. The general equation of hysteresis that has been described previously c now be solved with the incorporation of the isotropy d textured hysteretic agnetization given above to obtain the agnetization curves along particular directions. 3. Nuerical odel First, the equation for deterining the paraeter of reversible oveent of the doain wall - c, is defined. This paraeter c be deterined using the starting differential noral susceptibility χ in : χ in d ( c = 1 = dh H= = K δ α + c d dh (14 Where the fact that agnetization is copletely reversible for the starting agnetization curve = rev (i.e. d irr /dh is used. If eq. (2 d its derivative (for = are introduced in eq. (14, d the developent of the coth(h/a function for H is used, the following is obtained: c χ s in = a in, E.g. c = 3 χ. (15 3 a s The starting differential hysteresis susceptibility χ i defines a connection between paraeter a d α: d s χ i = lih = dh 3 (16 a α s So that: s a = 1 + α (17 3 χ i This equation is used as the connection between paraeters a d α. Differential susceptibility in the coercive point χ Hc, representing axial differential susceptibility on the whole agnetization curve, i.e. χ Hc =χ ax is used for deterining paraeter K. If eq. (6 is used to define the coercive point H = +H C d =, d it is trsfored, paraeter K is obtained in the following for: ( HC irr ( HC K = d ( HC d irr HC c dh ( χ ax dh (18 + α [ ( H ( H ] C irr C Eqs. (4 d (5 are used to deterine variables irr d d irr /dh in the coercive point as: c irr ( HC = ( HC, (19 1 c dirr ( HC 1 c d ( HC = χ ax dh 1 c 1 c dh If these expressions are introduced into eq.(18, the final expression used for the calculation of the pinning paraeter of the doain wall, K is obtained: ( HC 1 K = α + (2 1 c χ ax c d HC c ( 1 1 c dh A siilar procedure is repeated for the residual point (=+ r d H = using eqs. in the following order: (6 to define the so-called differential residual susceptibility χ r, then (4 d (5, then expressions for irr d d irr /dh are defined in this point, d the following is obtained: K r = ( r + (21 α c c d ( r χ r dh The axial excitation point ( = + d H = H c also be used. If this point is close to the saturation point, the noral d hysteresis differential susceptibility are alost the sae (χ n χ =χ d agnetization is practically equal to the irreversible value ( = irr. Differential susceptibility in the point of axial excitation c be defined using eq.(6 as: χ = K α d ( H + c dh ( H ( ( H d irr ( H dh (22

4 The approxiations ade in this point enable: dirr ( H d( H d ( H = = (23 dh dh dh And after exchging into (22 d trsforations: K χ = ( H (24. α χ The procedure used to solve the nuerical odel Having in ind that the equations used for defining paraeters of the Jiles-Atherton odel c only be iplicitly expressed as a function of the paraeter being calculated d the reaining odel paraeters, the procedure used for deterining these paraeters ust be iterative. Coparisons of values obtained in successive iterations are ade, d if there is no significt chge (for exaple: the deviation is saller th ξ=.1 the procedure is copleted, d the solutions obtained are treated as final values. As input values we use 1 paraeters (table 1. The progra was written using the atheatica environent, as it requires the nuerical solution of paraeter differential equations, which are supported by this progra packet. The paraeter of reversible oveent of the doain wall - c c be directly calculated using the starting noral susceptibility (eq. (15. The starting value of the averaging paraeter of the agnetic field is defined as α=.1, as this value is ost often found for all isotropic aterials, though other value (of the sae order of agnitude c be defined. The shape paraeter of the hysteresis curve - a c be initially calculated using eq. (17, while K is defined as. Calculation of the pinning paraeter of the doain wall K using eq. (2, followed by deterination of paraeter α using eq. (21 is perfored in a progra loop (which is exited when the difference between calculated paraeter values in two successive loops is less th ξ. Paraeter a is obtained in the sae loop using eq. (24. In this case it is necessary to ake a correct approxiation of hyperbolic functions. In this concrete case developing this function using Bernoulli nubers approxiates the coth function. This was later proved to be a correct ove, as this function easily diverges, thus disrupting the coplete calculation. By noralizing its arguent, we have retained the arguent value close to zero. The cubic equation obtained in this way was later solved using the atheatica 3.1 progra packet. At the sae tie, investigation of the nature of solutions obtained d selection of real ones were essential. The procedure for calculating paraeters K, a d α is repeated the necessary nuber of ties until the defined accuracy is attained, though the developed progra showed that results obtained converge to the final result in a sall nuber of steps (2-5. Negative values were also obtained as final paraeter values, which is noral, having in ind the nature of the hysteresis curve d the oveent direction on the curve, when the sign of coefficients used in the calculation is chged. This is why the absolute value of the calculated paraeter was taken in soe cases. This did not have y influence of further accuracy of the procedure. Fig. 1 shows the hysteresis curve for the n-zn ferrite - 3C8 obtained by introducing odel paraeters (calculated using our progra packet into the PSpice progra packet. Fig. 2 shows the hysteresis curve obtained using odel paraeters defined in the PSpice progra packet for the sae aterial. 5.K -5.K B(k1 H(k1 Fig. 1 Starting agnetization curve d hysteresis curve for the n-zn ferrite - 3C8 obtained using odel paraeter values (calculated by our progra in the PSpice packet. 5.K -5.K B(k1 H(k1 Fig. 2 Starting agnetization curve d hysteresis curve for the n-zn ferrite 3C8 obtained using odel paraeter values defined in the PSpice progra packet Good agreeent of both curves c be seen. agnetic characteristics read fro these curves are presented d copared with input values in table 1. We should p p

5 ention that the PSpice progra also gives the dependence of agnetic induction B on the agnetic field strength H, which is the case when experiental hysteresis curve is obtained. As these values do not figure in presented odel equations (1-24, their connection with odel variables is as follows: B( H ( H = H (25 μ o χ = μr 1 ; χ μr (26 Where μ r - the relative agnetic pereability (for soft agnetic aterials χ μ r. The unit for agnetic induction in the PSpice progra packet is Gauss - G (1G = 1-4 T, while Oersted - Oe (1Oe A/ is used as the agnetic field strength unit. Table 1 Input values of agnetic characteristics for the n-zn ferrite - 3C8, their values read fro the curve obtained using calculated odel paraeter (Fig. 1 d values obtained using odel paraeters defined in the PSpice progra packet input values wrote fro wrote fro Fig.2 Fig.2 B S (T χ i χ in H (A/ B (T χ B r (T χ r H c (A/ χ ax For hard agnetic core, due to the fact that one paraeter is redundt, we chose for the feedback paraeter a value α=1 without y loss of generality. R μ is obtained by a curve fit of the initial agnetization curve. The paraeters H 1, β 1 d H, β are deterined such that they approxiate the first d second quadrt, respectively. Fitting the second quadrt, initial guess for Kc, K, d H1 starts the coputation. Typical start values are K c =1...1, K=1, d H 1 =. For a atheatical treatent is hence replaced by in (7. We ded that the second ter of g(. is neglectible copared to the first in the second quadrt. Collects the linear part of (7, the curve fit is effected by selecting two points on deagnetization curve of the second quadrt, e.g. (H 1 =, 1 = r, where r is the reence agnetization d (H 2 =-H c /2, 2 =(H 2. The paraeters β, H are approxiately given by β ( c( H1 1 ( c( H 2 2 ln( 1 ln( 2 (, ln, ln, $ c H H 1 1 β 1 s H / s (, = + α $ c H H K c s R μ s If the slope of (7 becoes singular chges abruptly fro the upper to the lower brch d vice versa. The agnetization for which this singularity occurs is denoted as Hc. The field strength at Hc is therefore approxiately H c. The coercivity is nearly given by α 1 / R μ ( β H / s 1 c β H$ (28 H K H β H c H c $ c c ( α 1 / R μ s s After choosing a better guess for K c fro (28, the procedure is repeated aking use of (27. K is deterined by approxiating the slope of the agnetization curve near the corcivety aking use of (9. Now β 1 d H 1 are coputed by fitting the agnetization curve in the first quadrt. The paraeters ust be chosen such that they do not disturb the curve fit of the second quadrt, e.g., the second ter of g( / s in (7 ust be negligible for < r. We require that for (H 1 =, 1 = r β β (27 H$ 1 r 1 s H$ = k << 1 (29 β r s e.g., k=1-2. The paraeters β 1, H 1 are calculated by the sae reasoning as before by (27 replacing β by β 1 d H by H 1. The extended odel has been used to fit easured curves of Nd 2 Fe 14 B aterials, which have uniaxial isotropy. The odeled d easured curves are shown in Fig.3 d 4. The odeled curves show good agreeent with the experiental results. For the isotropic saple, the field was applied along its hard axis. The odeled texture level is t=.32, which es statistically 32 per cent of the saple was oriented along a specific direction d the rest was rdoly oriented.

6 Figure 3: odeled d easured curves of isotropic Nd 2 Fe 14 B aterial (s=61ka/; k=75ka/; a=42ka/; alpha=1.5; c=.6 Figure 5: Siulation curve of Alnico 2 core (β =1.2; H =12; K=8; K c =45; R μ =1; α=1 3 Figure 4: odeled d easured curves of textured Nd 2 Fe 14 B aterial (s=61ka/; k=75ka/; a=42ka/; alpha=2.3; c=.6; Ka=4266J/cu; t=.32; easy axis =1; field axis=1 Fig.5 d 6 copares easured d siulated hysteresis curve of Alnico 2 d Alnico 5 alloys. The paraeters were coputed using (27 d (28. The deviation of the easured to the siulated values is not significt for all two core aterials. Figure 6: Siulation curve of Alnico 5 core (β =16; H =216; K=2; K c =6; R μ =1; α=1 3 It c be noted that this algorith required a significt nuber of input data, not given by y producers of agnetic aterials, which necessitated the developent of a procedure for their calculation. We are able to experientally record the hysteresis curve d starting agnetization curve for y agnetic aterial. These curves are input into the atlab progra packet, which is used to deterine values of input data required by the iterative procedure. The required slope coefficients are obtained using slope values for two near points. The value of saturation induction B s c be calculated (predicted using the axial induction B value obtained fro the experiental curve, as the converging point. 5. Conclusion The proble of deterining paraeters of the Jiles- Atherton odel for R, Z, F-shaped hysteresis curves of agnetic aterials was alyzed in this paper. The algorith was tested at the siulation level, using the PSpice progra packet. Results obtained showed good agreeent with expected values. The presented algorith is of adaptive type in accordce with type of agnetic aterials. Expsion of this progra to classes

7 of isotropy (i.e. for Z d F-shaped hysteresis curves d hard agnetic aterials is ade, as corresponding physical-atheatical odels have been developed [3,6]. This enable the calculation of paraeters required in the alysis of coplex electronic circuits using the PSpice progra packet for a wide class of agnetic aterials. The proposed algorith suitable for on-line easureents. Finally, we obtained results which uch better then known results. References [1] D.C.Jiles d D.L.Atherton, Theory of Ferroagnetic Hysteresis, Journ. agn. agn. at. vol. 61, pp.48-6, [2]PSPICE User s Guide, icrosi Corp. Version 4.5., 1991 [3] A.Raesh, D.C.Jiles d J..Roderick, A odel of Anisotropic Anhysteretic agnetization, IEEE Trs. on agnetics AG vol. 32(5, pp , [4] D.C.Jiles, J.B.Thoelke, d.k.devine, Nuerical Deterination of Hysteresis Paraeters Using the Theory of Ferroagnetic Hysteresis, IEEE Trs. on agnetics AG vol. 28(1, pp.27-35, [5].Lengevin, Ann. de Che. et Phys. vol. 5, pp. 7, 195. [6] D.C.Jiles, A.Raesh, Y.Shi d X.Fg, Application of the Anisotropic Extension of the Theory of Hysteresis to the agnetization Curves of Crystalline d Textured agnetic aterials, IEEE Trs. on agnetics AG vol. 33(5, pp , 1997.

8.1 Force Laws Hooke s Law

8.1 Force Laws Hooke s Law 8.1 Force Laws There are forces that don't change appreciably fro one instant to another, which we refer to as constant in tie, and forces that don't change appreciably fro one point to another, which

More information

Data-Driven Imaging in Anisotropic Media

Data-Driven Imaging in Anisotropic Media 18 th World Conference on Non destructive Testing, 16- April 1, Durban, South Africa Data-Driven Iaging in Anisotropic Media Arno VOLKER 1 and Alan HUNTER 1 TNO Stieltjesweg 1, 6 AD, Delft, The Netherlands

More information

Experimental Design For Model Discrimination And Precise Parameter Estimation In WDS Analysis

Experimental Design For Model Discrimination And Precise Parameter Estimation In WDS Analysis City University of New York (CUNY) CUNY Acadeic Works International Conference on Hydroinforatics 8-1-2014 Experiental Design For Model Discriination And Precise Paraeter Estiation In WDS Analysis Giovanna

More information

Some Perspective. Forces and Newton s Laws

Some Perspective. Forces and Newton s Laws Soe Perspective The language of Kineatics provides us with an efficient ethod for describing the otion of aterial objects, and we ll continue to ake refineents to it as we introduce additional types of

More information

Kernel Methods and Support Vector Machines

Kernel Methods and Support Vector Machines Intelligent Systes: Reasoning and Recognition Jaes L. Crowley ENSIAG 2 / osig 1 Second Seester 2012/2013 Lesson 20 2 ay 2013 Kernel ethods and Support Vector achines Contents Kernel Functions...2 Quadratic

More information

Spine Fin Efficiency A Three Sided Pyramidal Fin of Equilateral Triangular Cross-Sectional Area

Spine Fin Efficiency A Three Sided Pyramidal Fin of Equilateral Triangular Cross-Sectional Area Proceedings of the 006 WSEAS/IASME International Conference on Heat and Mass Transfer, Miai, Florida, USA, January 18-0, 006 (pp13-18) Spine Fin Efficiency A Three Sided Pyraidal Fin of Equilateral Triangular

More information

Model Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon

Model Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon Model Fitting CURM Background Material, Fall 014 Dr. Doreen De Leon 1 Introduction Given a set of data points, we often want to fit a selected odel or type to the data (e.g., we suspect an exponential

More information

Ph 20.3 Numerical Solution of Ordinary Differential Equations

Ph 20.3 Numerical Solution of Ordinary Differential Equations Ph 20.3 Nuerical Solution of Ordinary Differential Equations Due: Week 5 -v20170314- This Assignent So far, your assignents have tried to failiarize you with the hardware and software in the Physics Coputing

More information

Optical Properties of Plasmas of High-Z Elements

Optical Properties of Plasmas of High-Z Elements Forschungszentru Karlsruhe Techni und Uwelt Wissenschaftlishe Berichte FZK Optical Properties of Plasas of High-Z Eleents V.Tolach 1, G.Miloshevsy 1, H.Würz Project Kernfusion 1 Heat and Mass Transfer

More information

Ch 12: Variations on Backpropagation

Ch 12: Variations on Backpropagation Ch 2: Variations on Backpropagation The basic backpropagation algorith is too slow for ost practical applications. It ay take days or weeks of coputer tie. We deonstrate why the backpropagation algorith

More information

Proc. of the IEEE/OES Seventh Working Conference on Current Measurement Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES

Proc. of the IEEE/OES Seventh Working Conference on Current Measurement Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES Proc. of the IEEE/OES Seventh Working Conference on Current Measureent Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES Belinda Lipa Codar Ocean Sensors 15 La Sandra Way, Portola Valley, CA 98 blipa@pogo.co

More information

This is a repository copy of Analytical optimisation of electromagnetic design of a linear (tubular) switched reluctance motor.

This is a repository copy of Analytical optimisation of electromagnetic design of a linear (tubular) switched reluctance motor. This is a repository copy of Analytical optiisation of electroagnetic design of a linear (tubular) switched reluctance otor. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk/907/

More information

DESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS *

DESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS * IJST, Transactions of Mechanical Engineering, Vol. 39, No. M1, pp 89-100 Printed in The Islaic Republic of Iran, 2015 Shira University DESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS

More information

1 Bounding the Margin

1 Bounding the Margin COS 511: Theoretical Machine Learning Lecturer: Rob Schapire Lecture #12 Scribe: Jian Min Si March 14, 2013 1 Bounding the Margin We are continuing the proof of a bound on the generalization error of AdaBoost

More information

ma x = -bv x + F rod.

ma x = -bv x + F rod. Notes on Dynaical Systes Dynaics is the study of change. The priary ingredients of a dynaical syste are its state and its rule of change (also soeties called the dynaic). Dynaical systes can be continuous

More information

N-Point. DFTs of Two Length-N Real Sequences

N-Point. DFTs of Two Length-N Real Sequences Coputation of the DFT of In ost practical applications, sequences of interest are real In such cases, the syetry properties of the DFT given in Table 5. can be exploited to ake the DFT coputations ore

More information

EE5900 Spring Lecture 4 IC interconnect modeling methods Zhuo Feng

EE5900 Spring Lecture 4 IC interconnect modeling methods Zhuo Feng EE59 Spring Parallel LSI AD Algoriths Lecture I interconnect odeling ethods Zhuo Feng. Z. Feng MTU EE59 So far we ve considered only tie doain analyses We ll soon see that it is soeties preferable to odel

More information

NUMERICAL MODELLING OF THE TYRE/ROAD CONTACT

NUMERICAL MODELLING OF THE TYRE/ROAD CONTACT NUMERICAL MODELLING OF THE TYRE/ROAD CONTACT PACS REFERENCE: 43.5.LJ Krister Larsson Departent of Applied Acoustics Chalers University of Technology SE-412 96 Sweden Tel: +46 ()31 772 22 Fax: +46 ()31

More information

Keywords: Estimator, Bias, Mean-squared error, normality, generalized Pareto distribution

Keywords: Estimator, Bias, Mean-squared error, normality, generalized Pareto distribution Testing approxiate norality of an estiator using the estiated MSE and bias with an application to the shape paraeter of the generalized Pareto distribution J. Martin van Zyl Abstract In this work the norality

More information

Extension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels

Extension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels Extension of CSRSM for the Paraetric Study of the Face Stability of Pressurized Tunnels Guilhe Mollon 1, Daniel Dias 2, and Abdul-Haid Soubra 3, M.ASCE 1 LGCIE, INSA Lyon, Université de Lyon, Doaine scientifique

More information

e-companion ONLY AVAILABLE IN ELECTRONIC FORM

e-companion ONLY AVAILABLE IN ELECTRONIC FORM OPERATIONS RESEARCH doi 10.1287/opre.1070.0427ec pp. ec1 ec5 e-copanion ONLY AVAILABLE IN ELECTRONIC FORM infors 07 INFORMS Electronic Copanion A Learning Approach for Interactive Marketing to a Custoer

More information

Analytical Study of AC Magnetic Susceptibility of (Bi, Pb) Sr-Ca-Cu-O Superconducting Systems Using Bean Critical State Model

Analytical Study of AC Magnetic Susceptibility of (Bi, Pb) Sr-Ca-Cu-O Superconducting Systems Using Bean Critical State Model Kasetsart J. (Nat. Sci.) 43 : 353-359 (9) Analytical Study of AC Magnetic Susceptibility of (Bi, Pb) Sr-Ca-Cu-O Superconducting Systes Using Bean Critical State Model Supreya Trivijitkase *, Adisak Boonchun,

More information

lecture 36: Linear Multistep Mehods: Zero Stability

lecture 36: Linear Multistep Mehods: Zero Stability 95 lecture 36: Linear Multistep Mehods: Zero Stability 5.6 Linear ultistep ethods: zero stability Does consistency iply convergence for linear ultistep ethods? This is always the case for one-step ethods,

More information

A Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair

A Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair Proceedings of the 6th SEAS International Conference on Siulation, Modelling and Optiization, Lisbon, Portugal, Septeber -4, 006 0 A Siplified Analytical Approach for Efficiency Evaluation of the eaving

More information

An Approximate Model for the Theoretical Prediction of the Velocity Increase in the Intermediate Ballistics Period

An Approximate Model for the Theoretical Prediction of the Velocity Increase in the Intermediate Ballistics Period An Approxiate Model for the Theoretical Prediction of the Velocity... 77 Central European Journal of Energetic Materials, 205, 2(), 77-88 ISSN 2353-843 An Approxiate Model for the Theoretical Prediction

More information

Kinetic Theory of Gases: Elementary Ideas

Kinetic Theory of Gases: Elementary Ideas Kinetic Theory of Gases: Eleentary Ideas 17th February 2010 1 Kinetic Theory: A Discussion Based on a Siplified iew of the Motion of Gases 1.1 Pressure: Consul Engel and Reid Ch. 33.1) for a discussion

More information

COS 424: Interacting with Data. Written Exercises

COS 424: Interacting with Data. Written Exercises COS 424: Interacting with Data Hoework #4 Spring 2007 Regression Due: Wednesday, April 18 Written Exercises See the course website for iportant inforation about collaboration and late policies, as well

More information

Reading from Young & Freedman: For this topic, read the introduction to chapter 25 and sections 25.1 to 25.3 & 25.6.

Reading from Young & Freedman: For this topic, read the introduction to chapter 25 and sections 25.1 to 25.3 & 25.6. PHY10 Electricity Topic 6 (Lectures 9 & 10) Electric Current and Resistance n this topic, we will cover: 1) Current in a conductor ) Resistivity 3) Resistance 4) Oh s Law 5) The Drude Model of conduction

More information

Multi Degrees of Freedom Maglev System with Permanent Magnet Motion Control

Multi Degrees of Freedom Maglev System with Permanent Magnet Motion Control Multi Degrees of Freedo Maglev Syste with Peranent Magnet Motion Control Cui Tianshi A dissertation subitted to Kochi University of Technology in partial fulfillent of the requireents for the degree of

More information

In this chapter, we consider several graph-theoretic and probabilistic models

In this chapter, we consider several graph-theoretic and probabilistic models THREE ONE GRAPH-THEORETIC AND STATISTICAL MODELS 3.1 INTRODUCTION In this chapter, we consider several graph-theoretic and probabilistic odels for a social network, which we do under different assuptions

More information

The Wilson Model of Cortical Neurons Richard B. Wells

The Wilson Model of Cortical Neurons Richard B. Wells The Wilson Model of Cortical Neurons Richard B. Wells I. Refineents on the odgkin-uxley Model The years since odgkin s and uxley s pioneering work have produced a nuber of derivative odgkin-uxley-like

More information

Figure 1: Equivalent electric (RC) circuit of a neurons membrane

Figure 1: Equivalent electric (RC) circuit of a neurons membrane Exercise: Leaky integrate and fire odel of neural spike generation This exercise investigates a siplified odel of how neurons spike in response to current inputs, one of the ost fundaental properties of

More information

Interactive Markov Models of Evolutionary Algorithms

Interactive Markov Models of Evolutionary Algorithms Cleveland State University EngagedScholarship@CSU Electrical Engineering & Coputer Science Faculty Publications Electrical Engineering & Coputer Science Departent 2015 Interactive Markov Models of Evolutionary

More information

Kinetic Theory of Gases: Elementary Ideas

Kinetic Theory of Gases: Elementary Ideas Kinetic Theory of Gases: Eleentary Ideas 9th February 011 1 Kinetic Theory: A Discussion Based on a Siplified iew of the Motion of Gases 1.1 Pressure: Consul Engel and Reid Ch. 33.1) for a discussion of

More information

On the characterization of non-linear diffusion equations. An application in soil mechanics

On the characterization of non-linear diffusion equations. An application in soil mechanics On the characterization of non-linear diffusion equations. An application in soil echanics GARCÍA-ROS, G., ALHAMA, I., CÁNOVAS, M *. Civil Engineering Departent Universidad Politécnica de Cartagena Paseo

More information

Ufuk Demirci* and Feza Kerestecioglu**

Ufuk Demirci* and Feza Kerestecioglu** 1 INDIRECT ADAPTIVE CONTROL OF MISSILES Ufuk Deirci* and Feza Kerestecioglu** *Turkish Navy Guided Missile Test Station, Beykoz, Istanbul, TURKEY **Departent of Electrical and Electronics Engineering,

More information

Chapter 4: Hypothesis of Diffusion-Limited Growth

Chapter 4: Hypothesis of Diffusion-Limited Growth Suary This section derives a useful equation to predict quantu dot size evolution under typical organoetallic synthesis conditions that are used to achieve narrow size distributions. Assuing diffusion-controlled

More information

2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all

2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all Lecture 6 Introduction to kinetic theory of plasa waves Introduction to kinetic theory So far we have been odeling plasa dynaics using fluid equations. The assuption has been that the pressure can be either

More information

Solving initial value problems by residual power series method

Solving initial value problems by residual power series method Theoretical Matheatics & Applications, vol.3, no.1, 13, 199-1 ISSN: 179-9687 (print), 179-979 (online) Scienpress Ltd, 13 Solving initial value probles by residual power series ethod Mohaed H. Al-Sadi

More information

Intelligent Systems: Reasoning and Recognition. Perceptrons and Support Vector Machines

Intelligent Systems: Reasoning and Recognition. Perceptrons and Support Vector Machines Intelligent Systes: Reasoning and Recognition Jaes L. Crowley osig 1 Winter Seester 2018 Lesson 6 27 February 2018 Outline Perceptrons and Support Vector achines Notation...2 Linear odels...3 Lines, Planes

More information

a a a a a a a m a b a b

a a a a a a a m a b a b Algebra / Trig Final Exa Study Guide (Fall Seester) Moncada/Dunphy Inforation About the Final Exa The final exa is cuulative, covering Appendix A (A.1-A.5) and Chapter 1. All probles will be ultiple choice

More information

Analysis of Impulsive Natural Phenomena through Finite Difference Methods A MATLAB Computational Project-Based Learning

Analysis of Impulsive Natural Phenomena through Finite Difference Methods A MATLAB Computational Project-Based Learning Analysis of Ipulsive Natural Phenoena through Finite Difference Methods A MATLAB Coputational Project-Based Learning Nicholas Kuia, Christopher Chariah, Mechatronics Engineering, Vaughn College of Aeronautics

More information

Numerical Solution of the MRLW Equation Using Finite Difference Method. 1 Introduction

Numerical Solution of the MRLW Equation Using Finite Difference Method. 1 Introduction ISSN 1749-3889 print, 1749-3897 online International Journal of Nonlinear Science Vol.1401 No.3,pp.355-361 Nuerical Solution of the MRLW Equation Using Finite Difference Method Pınar Keskin, Dursun Irk

More information

ESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS. A Thesis. Presented to. The Faculty of the Department of Mathematics

ESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS. A Thesis. Presented to. The Faculty of the Department of Mathematics ESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS A Thesis Presented to The Faculty of the Departent of Matheatics San Jose State University In Partial Fulfillent of the Requireents

More information

Sharp Time Data Tradeoffs for Linear Inverse Problems

Sharp Time Data Tradeoffs for Linear Inverse Problems Sharp Tie Data Tradeoffs for Linear Inverse Probles Saet Oyak Benjain Recht Mahdi Soltanolkotabi January 016 Abstract In this paper we characterize sharp tie-data tradeoffs for optiization probles used

More information

Accuracy of the Scaling Law for Experimental Natural Frequencies of Rectangular Thin Plates

Accuracy of the Scaling Law for Experimental Natural Frequencies of Rectangular Thin Plates The 9th Conference of Mechanical Engineering Network of Thailand 9- October 005, Phuket, Thailand Accuracy of the caling Law for Experiental Natural Frequencies of Rectangular Thin Plates Anawat Na songkhla

More information

Using a De-Convolution Window for Operating Modal Analysis

Using a De-Convolution Window for Operating Modal Analysis Using a De-Convolution Window for Operating Modal Analysis Brian Schwarz Vibrant Technology, Inc. Scotts Valley, CA Mark Richardson Vibrant Technology, Inc. Scotts Valley, CA Abstract Operating Modal Analysis

More information

EFFECT OF SURFACE ASPERITY TRUNCATION ON THERMAL CONTACT CONDUCTANCE

EFFECT OF SURFACE ASPERITY TRUNCATION ON THERMAL CONTACT CONDUCTANCE EFFECT OF SURFACE ASPERITY TRUNCATION ON THERMAL CONTACT CONDUCTANCE Fernando H. Milanez *, M. M. Yovanovich, J. R. Culha Microelectronics Heat Transfer Laboratory Departent of Mechanical Engineering University

More information

ACTIVE VIBRATION CONTROL FOR STRUCTURE HAVING NON- LINEAR BEHAVIOR UNDER EARTHQUAKE EXCITATION

ACTIVE VIBRATION CONTROL FOR STRUCTURE HAVING NON- LINEAR BEHAVIOR UNDER EARTHQUAKE EXCITATION International onference on Earthquae Engineering and Disaster itigation, Jaarta, April 14-15, 8 ATIVE VIBRATION ONTROL FOR TRUTURE HAVING NON- LINEAR BEHAVIOR UNDER EARTHQUAE EXITATION Herlien D. etio

More information

A note on the multiplication of sparse matrices

A note on the multiplication of sparse matrices Cent. Eur. J. Cop. Sci. 41) 2014 1-11 DOI: 10.2478/s13537-014-0201-x Central European Journal of Coputer Science A note on the ultiplication of sparse atrices Research Article Keivan Borna 12, Sohrab Aboozarkhani

More information

( ') ( ) 3. Magnetostatic Field Introduction

( ') ( ) 3. Magnetostatic Field Introduction 3. Magnetostatic Field 3.. Introduction A agnetostatic field is a agnetic field produced by electric charge in peranent unifor oveent, i.e. in a peranent oveent with constant velocity. Any directed oveent

More information

Analyzing Simulation Results

Analyzing Simulation Results Analyzing Siulation Results Dr. John Mellor-Cruey Departent of Coputer Science Rice University johnc@cs.rice.edu COMP 528 Lecture 20 31 March 2005 Topics for Today Model verification Model validation Transient

More information

Research and Experiments on Electromagnetic Field Induced by Two. Coaxial Solenoid Coils of Axially Mag-lev Driving Needle

Research and Experiments on Electromagnetic Field Induced by Two. Coaxial Solenoid Coils of Axially Mag-lev Driving Needle 3rd International Conference on Mechatronics and Inforation Technology (ICMIT 16) Research and Experients on Electroagnetic Field Induced by Two Coaxial Solenoid Coils of Axially Mag-lev Driving Needle

More information

Computer Model For Sieves Vibrations Analysis, Using an Algorithm Based on the False-Position Method

Computer Model For Sieves Vibrations Analysis, Using an Algorithm Based on the False-Position Method Aerican Journal of Applied Sciences 6 (): 48-56, 9 ISSN 546-939 9 Science Publications Coputer Model For Sieves Vibrations Analysis, Using an Algorith Based on the False-Position Method Dinu I. Stoicovici,

More information

Self-Erecting Inverted Pendulum: Swing up and Stabilization Control

Self-Erecting Inverted Pendulum: Swing up and Stabilization Control Self-Erecting Inverted Pendulu: Swing up and Stabilization Control S. McGilvray (Winner of the IEEE Life Meber ward for best paper fro the Central Canada Council for the IEEE Student Paper Contest, ) Contents

More information

MODIFIED SERIES RESISTANCE MODEL - DETERMINATION OF MEAN CONCENTRATION BY INTEGRAL TRANSFORMATION

MODIFIED SERIES RESISTANCE MODEL - DETERMINATION OF MEAN CONCENTRATION BY INTEGRAL TRANSFORMATION MODIFIED SERIES RESISTANCE MODEL - DETERMINATION OF MEAN CONCENTRATION BY INTEGRAL TRANSFORMATION A. L. Venezuela a, R. F. Cantão a, R. S. Ongaratto b, and R. N. Haneda c a Universidade Federal de São

More information

Efficient Filter Banks And Interpolators

Efficient Filter Banks And Interpolators Efficient Filter Banks And Interpolators A. G. DEMPSTER AND N. P. MURPHY Departent of Electronic Systes University of Westinster 115 New Cavendish St, London W1M 8JS United Kingdo Abstract: - Graphical

More information

Topic 5a Introduction to Curve Fitting & Linear Regression

Topic 5a Introduction to Curve Fitting & Linear Regression /7/08 Course Instructor Dr. Rayond C. Rup Oice: A 337 Phone: (95) 747 6958 E ail: rcrup@utep.edu opic 5a Introduction to Curve Fitting & Linear Regression EE 4386/530 Coputational ethods in EE Outline

More information

Symbolic Analysis as Universal Tool for Deriving Properties of Non-linear Algorithms Case study of EM Algorithm

Symbolic Analysis as Universal Tool for Deriving Properties of Non-linear Algorithms Case study of EM Algorithm Acta Polytechnica Hungarica Vol., No., 04 Sybolic Analysis as Universal Tool for Deriving Properties of Non-linear Algoriths Case study of EM Algorith Vladiir Mladenović, Miroslav Lutovac, Dana Porrat

More information

Statistical Logic Cell Delay Analysis Using a Current-based Model

Statistical Logic Cell Delay Analysis Using a Current-based Model Statistical Logic Cell Delay Analysis Using a Current-based Model Hanif Fatei Shahin Nazarian Massoud Pedra Dept. of EE-Systes, University of Southern California, Los Angeles, CA 90089 {fatei, shahin,

More information

Easy Evaluation Method of Self-Compactability of Self-Compacting Concrete

Easy Evaluation Method of Self-Compactability of Self-Compacting Concrete Easy Evaluation Method of Self-Copactability of Self-Copacting Concrete Masanori Maruoka 1 Hiroi Fujiwara 2 Erika Ogura 3 Nobu Watanabe 4 T 11 ABSTRACT The use of self-copacting concrete (SCC) in construction

More information

Support Vector Machine Classification of Uncertain and Imbalanced data using Robust Optimization

Support Vector Machine Classification of Uncertain and Imbalanced data using Robust Optimization Recent Researches in Coputer Science Support Vector Machine Classification of Uncertain and Ibalanced data using Robust Optiization RAGHAV PAT, THEODORE B. TRAFALIS, KASH BARKER School of Industrial Engineering

More information

Supplementary Information for Design of Bending Multi-Layer Electroactive Polymer Actuators

Supplementary Information for Design of Bending Multi-Layer Electroactive Polymer Actuators Suppleentary Inforation for Design of Bending Multi-Layer Electroactive Polyer Actuators Bavani Balakrisnan, Alek Nacev, and Elisabeth Sela University of Maryland, College Park, Maryland 074 1 Analytical

More information

ABSTRACT INTRODUCTION

ABSTRACT INTRODUCTION Wave Resistance Prediction of a Cataaran by Linearised Theory M.INSEL Faculty of Naval Architecture and Ocean Engineering, Istanbul Technical University, TURKEY A.F.MOLLAND, J.F.WELLICOME Departent of

More information

Effects of an Inhomogeneous Magnetic Field (E =0)

Effects of an Inhomogeneous Magnetic Field (E =0) Effects of an Inhoogeneous Magnetic Field (E =0 For soe purposes the otion of the guiding centers can be taken as a good approxiation of that of the particles. ut it ust be recognized that during the particle

More information

Optimization of ripple filter for pencil beam scanning

Optimization of ripple filter for pencil beam scanning Nuclear Science and Techniques 24 (2013) 060404 Optiization of ripple filter for pencil bea scanning YANG Zhaoxia ZHANG Manzhou LI Deing * Shanghai Institute of Applied Physics, Chinese Acadey of Sciences,

More information

IN modern society that various systems have become more

IN modern society that various systems have become more Developent of Reliability Function in -Coponent Standby Redundant Syste with Priority Based on Maxiu Entropy Principle Ryosuke Hirata, Ikuo Arizono, Ryosuke Toohiro, Satoshi Oigawa, and Yasuhiko Takeoto

More information

A Better Algorithm For an Ancient Scheduling Problem. David R. Karger Steven J. Phillips Eric Torng. Department of Computer Science

A Better Algorithm For an Ancient Scheduling Problem. David R. Karger Steven J. Phillips Eric Torng. Department of Computer Science A Better Algorith For an Ancient Scheduling Proble David R. Karger Steven J. Phillips Eric Torng Departent of Coputer Science Stanford University Stanford, CA 9435-4 Abstract One of the oldest and siplest

More information

THE KALMAN FILTER: A LOOK BEHIND THE SCENE

THE KALMAN FILTER: A LOOK BEHIND THE SCENE HE KALMA FILER: A LOOK BEHID HE SCEE R.E. Deain School of Matheatical and Geospatial Sciences, RMI University eail: rod.deain@rit.edu.au Presented at the Victorian Regional Survey Conference, Mildura,

More information

Chapter 6 1-D Continuous Groups

Chapter 6 1-D Continuous Groups Chapter 6 1-D Continuous Groups Continuous groups consist of group eleents labelled by one or ore continuous variables, say a 1, a 2,, a r, where each variable has a well- defined range. This chapter explores:

More information

12 Towards hydrodynamic equations J Nonlinear Dynamics II: Continuum Systems Lecture 12 Spring 2015

12 Towards hydrodynamic equations J Nonlinear Dynamics II: Continuum Systems Lecture 12 Spring 2015 18.354J Nonlinear Dynaics II: Continuu Systes Lecture 12 Spring 2015 12 Towards hydrodynaic equations The previous classes focussed on the continuu description of static (tie-independent) elastic systes.

More information

Measuring Temperature with a Silicon Diode

Measuring Temperature with a Silicon Diode Measuring Teperature with a Silicon Diode Due to the high sensitivity, nearly linear response, and easy availability, we will use a 1N4148 diode for the teperature transducer in our easureents 10 Analysis

More information

Use of PSO in Parameter Estimation of Robot Dynamics; Part One: No Need for Parameterization

Use of PSO in Parameter Estimation of Robot Dynamics; Part One: No Need for Parameterization Use of PSO in Paraeter Estiation of Robot Dynaics; Part One: No Need for Paraeterization Hossein Jahandideh, Mehrzad Navar Abstract Offline procedures for estiating paraeters of robot dynaics are practically

More information

Block designs and statistics

Block designs and statistics Bloc designs and statistics Notes for Math 447 May 3, 2011 The ain paraeters of a bloc design are nuber of varieties v, bloc size, nuber of blocs b. A design is built on a set of v eleents. Each eleent

More information

Anisotropic reference media and the possible linearized approximations for phase velocities of qs waves in weakly anisotropic media

Anisotropic reference media and the possible linearized approximations for phase velocities of qs waves in weakly anisotropic media INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS D: APPLIED PHYSICS J. Phys. D: Appl. Phys. 5 00 007 04 PII: S00-770867-6 Anisotropic reference edia and the possible linearized approxiations for phase

More information

2.003 Engineering Dynamics Problem Set 2 Solutions

2.003 Engineering Dynamics Problem Set 2 Solutions .003 Engineering Dynaics Proble Set Solutions This proble set is priarily eant to give the student practice in describing otion. This is the subject of kineatics. It is strongly recoended that you study

More information

Kinematics and dynamics, a computational approach

Kinematics and dynamics, a computational approach Kineatics and dynaics, a coputational approach We begin the discussion of nuerical approaches to echanics with the definition for the velocity r r ( t t) r ( t) v( t) li li or r( t t) r( t) v( t) t for

More information

Ştefan ŞTEFĂNESCU * is the minimum global value for the function h (x)

Ştefan ŞTEFĂNESCU * is the minimum global value for the function h (x) 7Applying Nelder Mead s Optiization Algorith APPLYING NELDER MEAD S OPTIMIZATION ALGORITHM FOR MULTIPLE GLOBAL MINIMA Abstract Ştefan ŞTEFĂNESCU * The iterative deterinistic optiization ethod could not

More information

IAENG International Journal of Computer Science, 42:2, IJCS_42_2_06. Approximation Capabilities of Interpretable Fuzzy Inference Systems

IAENG International Journal of Computer Science, 42:2, IJCS_42_2_06. Approximation Capabilities of Interpretable Fuzzy Inference Systems IAENG International Journal of Coputer Science, 4:, IJCS_4 6 Approxiation Capabilities of Interpretable Fuzzy Inference Systes Hirofui Miyajia, Noritaka Shigei, and Hiroi Miyajia 3 Abstract Many studies

More information

TABLE FOR UPPER PERCENTAGE POINTS OF THE LARGEST ROOT OF A DETERMINANTAL EQUATION WITH FIVE ROOTS. By William W. Chen

TABLE FOR UPPER PERCENTAGE POINTS OF THE LARGEST ROOT OF A DETERMINANTAL EQUATION WITH FIVE ROOTS. By William W. Chen TABLE FOR UPPER PERCENTAGE POINTS OF THE LARGEST ROOT OF A DETERMINANTAL EQUATION WITH FIVE ROOTS By Willia W. Chen The distribution of the non-null characteristic roots of a atri derived fro saple observations

More information

Compression and Predictive Distributions for Large Alphabet i.i.d and Markov models

Compression and Predictive Distributions for Large Alphabet i.i.d and Markov models 2014 IEEE International Syposiu on Inforation Theory Copression and Predictive Distributions for Large Alphabet i.i.d and Markov odels Xiao Yang Departent of Statistics Yale University New Haven, CT, 06511

More information

Magnetorheological Fluid for Levitation Migration Technology

Magnetorheological Fluid for Levitation Migration Technology 66 The Open Materials Science Journal, 2011, 5, 66-71 Magnetorheological Fluid for Levitation Migration Technology Guang Yang * and ZhuPing Chen * Open Access Faculty of Mechanical Engineering, Jiei University,

More information

Density and structure of undercooled liquid titanium

Density and structure of undercooled liquid titanium Article Condensed Matter Physics March 2012 Vol.57 No.7: 719 723 doi: 10.1007/s11434-011-4945-6 Density and structure of undercooled liquid titaniu WANG HaiPeng, YANG ShangJing & WEI BingBo * Departent

More information

IN A SENSE, every material is a composite, even if the

IN A SENSE, every material is a composite, even if the IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES, VOL. 47, NO. 11, NOVEMBER 1999 2075 Magnetis fro Conductors and Enhanced Nonlinear Phenoena J. B. Pendry, A. J. Holden, D. J. Robbins, and W. J. Stewart,

More information

Generalized Rayleigh Wave Dispersion in a Covered Half-space Made of Viscoelastic Materials

Generalized Rayleigh Wave Dispersion in a Covered Half-space Made of Viscoelastic Materials Copyright 7 Tech Science Press CMC vol.53 no.4 pp.37-34 7 Generalized Rayleigh Wave Dispersion in a Covered Half-space Made of Viscoelastic Materials S.D. Akbarov and M. Negin 3 Abstract: Dispersion of

More information

Engineering of Composite Media for Shields at Microwave Frequencies

Engineering of Composite Media for Shields at Microwave Frequencies Engineering of Coposite Media for Shields at Microwave Frequencies Marina Koledintseva, Poorna Chander avva, ichard DuBroff, Jaes Drewniak EMC Laboratory, ECE Departent, University of Missouri-olla, olla,

More information

Nonmonotonic Networks. a. IRST, I Povo (Trento) Italy, b. Univ. of Trento, Physics Dept., I Povo (Trento) Italy

Nonmonotonic Networks. a. IRST, I Povo (Trento) Italy, b. Univ. of Trento, Physics Dept., I Povo (Trento) Italy Storage Capacity and Dynaics of Nononotonic Networks Bruno Crespi a and Ignazio Lazzizzera b a. IRST, I-38050 Povo (Trento) Italy, b. Univ. of Trento, Physics Dept., I-38050 Povo (Trento) Italy INFN Gruppo

More information

Mutual capacitor and its applications

Mutual capacitor and its applications Mutual capacitor and its applications Chun Li, Jason Li, Jieing Li CALSON Technologies, Toronto, Canada E-ail: calandli@yahoo.ca Published in The Journal of Engineering; Received on 27th October 2013;

More information

Scattering and bound states

Scattering and bound states Chapter Scattering and bound states In this chapter we give a review of quantu-echanical scattering theory. We focus on the relation between the scattering aplitude of a potential and its bound states

More information

3D acoustic wave modeling with a time-space domain dispersion-relation-based Finite-difference scheme

3D acoustic wave modeling with a time-space domain dispersion-relation-based Finite-difference scheme P-8 3D acoustic wave odeling with a tie-space doain dispersion-relation-based Finite-difference schee Yang Liu * and rinal K. Sen State Key Laboratory of Petroleu Resource and Prospecting (China University

More information

Inspection; structural health monitoring; reliability; Bayesian analysis; updating; decision analysis; value of information

Inspection; structural health monitoring; reliability; Bayesian analysis; updating; decision analysis; value of information Cite as: Straub D. (2014). Value of inforation analysis with structural reliability ethods. Structural Safety, 49: 75-86. Value of Inforation Analysis with Structural Reliability Methods Daniel Straub

More information

A New Algorithm for Reactive Electric Power Measurement

A New Algorithm for Reactive Electric Power Measurement A. Abiyev, GAU J. Soc. & Appl. Sci., 2(4), 7-25, 27 A ew Algorith for Reactive Electric Power Measureent Adalet Abiyev Girne Aerican University, Departernt of Electrical Electronics Engineering, Mersin,

More information

A model reduction approach to numerical inversion for a parabolic partial differential equation

A model reduction approach to numerical inversion for a parabolic partial differential equation Inverse Probles Inverse Probles 30 (204) 250 (33pp) doi:0.088/0266-56/30/2/250 A odel reduction approach to nuerical inversion for a parabolic partial differential equation Liliana Borcea, Vladiir Drusin

More information

Detection and Estimation Theory

Detection and Estimation Theory ESE 54 Detection and Estiation Theory Joseph A. O Sullivan Sauel C. Sachs Professor Electronic Systes and Signals Research Laboratory Electrical and Systes Engineering Washington University 11 Urbauer

More information

Chapter 2: Introduction to Damping in Free and Forced Vibrations

Chapter 2: Introduction to Damping in Free and Forced Vibrations Chapter 2: Introduction to Daping in Free and Forced Vibrations This chapter ainly deals with the effect of daping in two conditions like free and forced excitation of echanical systes. Daping plays an

More information

2.141 Modeling and Simulation of Dynamic Systems Assignment #2

2.141 Modeling and Simulation of Dynamic Systems Assignment #2 2.141 Modeling and Siulation of Dynaic Systes Assignent #2 Out: Wednesday Septeber 20, 2006 Due: Wednesday October 4, 2006 Proble 1 The sketch shows a highly siplified diagra of a dry-dock used in ship

More information

Chapter 1: Basics of Vibrations for Simple Mechanical Systems

Chapter 1: Basics of Vibrations for Simple Mechanical Systems Chapter 1: Basics of Vibrations for Siple Mechanical Systes Introduction: The fundaentals of Sound and Vibrations are part of the broader field of echanics, with strong connections to classical echanics,

More information

A remark on a success rate model for DPA and CPA

A remark on a success rate model for DPA and CPA A reark on a success rate odel for DPA and CPA A. Wieers, BSI Version 0.5 andreas.wieers@bsi.bund.de Septeber 5, 2018 Abstract The success rate is the ost coon evaluation etric for easuring the perforance

More information

Smith Predictor Based-Sliding Mode Controller for Integrating Process with Elevated Deadtime

Smith Predictor Based-Sliding Mode Controller for Integrating Process with Elevated Deadtime Sith Predictor Based-Sliding Mode Controller for Integrating Process with Elevated Deadtie Oscar Caacho, a, * Francisco De la Cruz b a Postgrado en Autoatización e Instruentación. Grupo en Nuevas Estrategias

More information

COULD A VARIABLE MASS OSCILLATOR EXHIBIT THE LATERAL INSTABILITY?

COULD A VARIABLE MASS OSCILLATOR EXHIBIT THE LATERAL INSTABILITY? Kragujevac J. Sci. 3 (8) 3-44. UDC 53.35 3 COULD A VARIABLE MASS OSCILLATOR EXHIBIT THE LATERAL INSTABILITY? Nebojša Danilović, Milan Kovačević and Vukota Babović Institute of Physics, Faculty of Science,

More information