Vés al contingut

Mètode de Newmark

De la Viquipèdia, l'enciclopèdia lliure

El Mètode beta de Newmark és un mètode d'integració numèrica emprat per resoldre equacions diferencials. S'utilitza en anàlisi d'element finits per modelar sistemes dinàmics.

Fent memòria de l'equació del moviment d'un cos en temps continu i en una dimensió,

Emprant el teorema del valor mitjà, el mètode de Newmark diu que la primera derivada temporal (la velocitat a l' equació del moviment) es pot calcular de la següent manera,

on i són dos instants de temps consecutius i és la proporció entre ambdós instants de temps a la qual el valor de l'acceleració satisfà el teorema del valor mitjà. Aquest valor de l'acceleració es pot estimar com

per tant

Tanmateix, com que l'acceleració varia en el temps, per tal de no perdre precisió en l'aproximació numèrica i obtenir el desplaçament correcte, cal calcular també la posició emprant el teorema del valor mitja ha d'ésser estès també a la segona derivada temporal (és a dir, les acceleracions):

i l'acceleració es pot estimar com

per tant

Aquest mètode rep el nom de Nathan M. Newmark, que el va introduir 1959.[1]

Aplicació

[modifica]

Les regles d'actualització de l'algorisme per a uns valors de i donats queda


Les característiques de l'algorisme depenen dels dos paràmetres i , en particular l'estabilitat i el seu caràcter implícit o explícit. Newmark va mostrar que el valor raonable de és 0.5, i normalment, encara que no sempre, s'utilitza aquest valor. Valors diferents de β entre 0 i 1 poden donar un ampli ventall de resultats. Normalment es pren β = 1/4, que resulta en el mètode d'acceleració mitjana constant.

Domini Estabilitat
inestable
i condicionalment estable
i incondicionalment estable

El mètode ha estat aplicat a casos particulars emprant valors adequats de i també :[2]

Nom del mètode Propietats
Explícit pur 0 0 explícit. Límit d'estabilitat
Diferències centrades 0 explícit i condicionalment estable
Fox Goodwin condicionalment estable
Acceleració lineal condicionalment estable
Acceleració mitjana incondicionalment estable

Aplicació a l'equació dinàmica

[modifica]

El mètode s'aplica a la simulació de sistemes mecànics que respon al sistema d'equacions diferencials

que, per a un pas de temps donat, es pot expressar com

on

Cal doncs recalcular les forces per cada coordenada (o grau de llibertat) i velocitat segons l'algoritme següent emprant els paràmetres i

INICI
      ; {Els valors del pas anterior es converteixen en els valors antics}
      ; {Predicció de les noves acceleracions en funció dels valors antics}
      ; {Predicció de les noves velocitats en funció de les noves acceleracions}
      ; {Predicció de les noves posicions en funció de les noves acceleracions}
      REPETIR
            ; {Els valors nous esdevenen valors antics}
            ; {Càlcul de les noves acceleracions en funció dels valors antics}
            ; {Càlcul de les noves velocitats en funció de les noves acceleracions}
            ; {Càlcul de les noves posicions en funció de les noves acceleracions}
      FINS QUE  {La nova aproximació és prou igual a l'antiga}
FINAL

El mètode de Newmark és una eina pràctica aplicable a problemes no lineals que poden arribar a ésser extremadament complexes. Convergeix raonablement quan el pas de temps és prou petit respecte les freqüències del sistema modelat. Tanmateix, a diferència del seu ús en problemes lineals, no es pot assegurar pas la seva convergència incondicional.

Referències i notes

[modifica]
  1. Newmark N. M., A Method of Computation for Structural Dynamics, Proc. ASCE 85(EM3), 67-94, 1959
  2. Géradin M., Rixen D. Mechanicl Vibrations. Theory and Application to Structural Dynamics, 2a Edició, Wiley, 1997

Vegeu també

[modifica]