Aula4_2 Adição Pto Ct Em 2k_Parte II

Embed Size (px)

Citation preview

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    1/28

    Universidade Federal de Campina Grande - Prof. André Luiz Fiquenede Brito

    Aula 4.2 – Metodologia de Superfície de

    Resposta (MSR) Parte II: Adição de Ptt e! 2"

    (Método de Máxima Inclinação Ascendente)

    1

    Universidade Federal de Campina Grande –UFCG

    Centro de Cin!ias e "e!nolo#ia – CC"

    Unidade A!admi!a de $n#en%aria &u'mi!a –UA$&

    Planejamento Experimental e Otimização de ProcessosANDRÉ LUIZ I!UENE DE "RI#O$ Dr %

    UFCG/UAEQ/CCT [email protected]  ou [email protected]

    Campina Grande PB –

    mailto:[email protected]:[email protected]:[email protected]:[email protected]

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    2/28

    Exemplo: Otimização de um Modelo de 1a Ordem

    ConversãoPercentual?

    Fator ATempo

    Fator B

    Temperatura

    Exemplo: Indústria Química

    m !rupo de pes"uisadores

    estão estudando a conversãopercentual#$% de um processo&

    '( duas vari(veis de interesse:

     A: tempo de rea)ão

    B: Temperatura de rea)ão

    &%  Inicialmente

    não '( certeza daLinearidade da

    re)ião de exploração****

    +% Decide,se cond-zir-m planejamento

    .atorial com + n/0eis$com -ma 1nica

    r2plica de corrida emcada ponto .atorial e

    a-mentando 3pontos centrais%

    Universidade Federal de Campina Grande - Prof. André Luiz Fiquene de Brito

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    3/28

    3

    ii4 .azer ossinais

    al)25ricosiii4 Representar 

    )eom2tricamente++63P7

    I4Ac'ar o n%

    de exp%para o planj%atorial com

    + n/0eis%

    i04 Ac'ar a8!erro

    04 Ac'ara

     9!c-r0at-ra

    0i4 Interpretar ores-ltado:

    comparar9!c-r0at-racom 8!erro;razão4< e 0ii4Usar o 8INI#A"para ANO=A%

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    4/28

    Solução:

    i% *úmero deexperimentopara o ++ , -

    PtCt

    i

    % *úmero deexperimentopara o ++ , -

    PtCt

    #$p % Ptos fatoriais & Ptt

    #$p % 2" & nc

    #$p % 22 & '

    #$p % 4 & '.:#$p%e$peri!etos(o*s)

    Universidade Federal de Campina Grande - Prof. André Luiz Fiquene de Brito

    &% #empo: >?;,4 e @?;64 min%%%%%%%%%%%%%%%%%%%%%%%%%%%Pt7t:

    >3 min-tos ;?4

    +% #emperat-ra: &3?;,4 e &?;64%%%%%%%%%%%%%%%%%%%%% Pt7t:&33 )ra-s ;?4

    +

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    5/28

    Sinais algébricos para clculos dos efeitos

    Tratamento

    (1)

    a

     b

    ab

    Efeito Fatorial

    I A B AxB+ −   −   +

    + + −   −

    + −  + −

    + + + +

    Y (%)

    y1=39,3

    y2=4,9

    y3=4,

    y4=41,!

    y!= 4,3y"=4,!

    y#=4,#

    y$=4,2

    !

    "

    #

    $

    9

    n1

    n2

    n3

    n4

    +

    +

    +

    +

    +

    En&aio

    1

    2

    3

    4

    n!

    Universidade Federal de Campina Grande - Prof. André Luiz Fiquene de Brito

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    6/28

    40.3

    40.5

    40.7

    40.240.6

    41.540.0

    40.939.3150 -1

    155 0

    160 1

    -1

    30

    0

    35

    1

    40

    tempo

         t    e    m    p    e    r    a     t    u    r    a

    "epresentação geométrica para 2 2  # !$t%t 

    #.% a

    / a/

    0

    Universidade Federal de Campina Grande - Prof. André Luiz Fiquene de Brito

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    7/28

     'char '()*' com 'u+,lio do S)- /('1

     

    1alorP

    ou

    23AFIC45

    1alor F

    4u

    Teste t

     A*41A   A    3    4

         6   I    *   I    T   A    B

     +A,#-A A A/01A22 & ' Rep. Ptt

     +A,#-A A A/01A22 & ' Rep. Ptt

    Crit7rio de 8ecisãoCrit7rio de 8ecisão

     +estes para o#feito pricipalcur3atura(efeito

    uadr5tico) eiteração

     +estes para o#feito pricipalcur3atura(efeito

    uadr5tico) eiteração

    Con9orme o α esta/elecido

    5A8A 6I*ITAB5A8A 6I*ITAB.& E9eitos e

    Coe9icientesestimados#si!%

    +& 5Q curvatura

    ;& Q6 erro

    sar 6inita/sar 6inita/5TAT84EFAT43IA< C3EATE4=

    +>

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    8/28

    SAÍDA DO MINITAB PAA ANO!A

    @

    Universidade Federal de Campina Grande - Prof. André Luiz Fiquene de Brito

    O"#1.$ P%a&e'ame&to ()p. *om Po&to *e&tra% te#tar+ #e

    o mo,e%o ua,r+t/o ou %/&ear

    O"#2.$ Se o te#te mo#trar a %/&ear/,a,e ,ee-#e e&errar

    a a&+%/#e e aar o po&to t/mo po/# po##/e%me&te omo,e%o a,ota,o #er/r+ para /&terpretar e#t/mar e o"tero&%u#6e# #/&/8/at/a# #o"re o e)per/me&to

    O"#3.$ Se o te#te mo#trar ue oueuratura9ua,r+t/o: ,ee-#e o&t/&uar o e)per/me&toamp%/a&,o o# &;e/#.

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    9/28

    B E.eito !-adr(tico C E.eito na 7-r0at-ra

    B 7-r0at-re: Representa o e.eito -adr(tico Indica -ede0e,se ampliar n/0eis%

    9e P!alor   F ?$?3  N9 para o teste -adr(tico ;não '(c-r0a4<

    Lo)o$ ?$G&@ F ?$?3% Então de0e,se adotar o modelo linear%

     Analsis o9 1ariance 9or #coded units%

    5ource 8F 5e" 55 AdD 55 AdD 65 F P6ain E99ects + +@+-00 +@+-00 ..+-0 ;+@- 000; si!  A . +0+-0 +0+-0 +0+-0 --@G 000+ si!  B . 0++-0 0++-0 0++-0 @; 00;- si!+>Ha Interactions . 000+-0 000+-0 000+-0 00 0@+. ns  AJB . 000+-0 000+-0 000+-0 00 0@+. ns  Curvature . 000+G+ 000+G+ 000+G+ 00 0@. *ão K( curva3esidual Error 0.G+00 0.G+00 00;00

      Pure Error 0.G+00 0.G+00 00;00Total @ ;00+++

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    10/28

    Termo Efeito 'oefiiente  p valor

    'on&tante 4,42 , * ,1 i Entra nomo-elo)

    Fator A 1,!! ,##! ,! * ,1 i Entra no mo-elo)

    Fator B ,"! ,32! ,#4 * ,1 i Entra no mo-elo)

    AxB ,! ,2! ,$"2 . ,1  /0o i   /0o Entra no mo-elo)

    't ,$! ,""4 . ,1 /0o i  /0o Entra no mo-elo)

     y=β 0+β 

    1 A+β 

    2 B+β 

    3 AB

    4 modelo !en7rico 7:

    O modelo codi.icado est( apresentado a5aixo:

    .0

    ,32!B,##!A4,42   ++

    2132!,##!,42,4   x x  y   ++=

    =  y

    Universidade Federal de Campina Grande - Prof. André Luiz Fiquene de Brito

      p Valor 

    MOD(

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    11/28

    (-, - ire/o da 02ima 3n!lina/o As!endente -03A

    11

    Continua/o daParte 3

    Universidade Federal de Campina Grande - Prof. André Luiz Fiquene de Brito

     Exemplo Conversão Percentual .$ 22 = 5 Pt*t

      O# ,a,o# ,o e)per/me&to mo#traram ue o mo,e%o ,e 1>or,em po,e #er u#a,o. Para e&o&trar a MáximaInclinação Ascendente – MIA ,ee-#e o/,erar ue$

    40.3

    40.5

    40.7

    40.2

    40.6

    41.540.0

    40.939.3150 -1

    155 0

    160 1

    -1

    30

    0

    35

    1

    40

    tempo

         t    e    m    p    e    r    a     t    u    r    a

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    12/28

    1 Etaa5 'om lane6amento 22  = 5 Pt*t, o

    mo-elo obti-o foi5 y = 40,42 +0,775 t + 0,325 T 

    7b&tit7in-o o& maiore& n8ei& o-ifia-o& (t=+1) e (T = +1) o7 t = 4 min e T = 1", tem&e5

     y ? 41,50%. (apenas codificado 1 e !1"

    1

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    13/28

    2  Etaa) /a& on-i:;e& -etemo, t = 3! min() eTemerat7ra = 1!!(), oren-imento < exatamente5

     y # 40,4$%. (apenas codificado1 e !1"

    1

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    14/28

    3 Etaa)o-e&e onl7ir 7e moen-o&e -o entro -o

    exerimento, om t= e T=, o ren-imentoa7menta>? o&&8el 7e oorra a m@xima inlina:0o

    a&e-ente ara obten:0o -o timo

    1

    Exlia:0o5

    'om t = e T = 5 y  4,42%'om t = +1 e T = +1 5  y = 41,!%(Co ontoentral mai& 7m a&&o D -ireita o alor -e y

    a7mento7 -e 4,42 ara 41,!%)

    Universidade Federal de Campina Grande Prof André Luiz Fiquene de Brito

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    15/28

    4  Etaa)  A&&im, n0o exi&te nen7ma ra0o ara 7e&tionar o mo-elo -e rimeira or-em>

    ! Etaa) ara obtermo& a Máxima

    Inclinação Ascendente  de3e9seaalisar o !odelo do plaea!eto22 = 5 Pt*t e o mo,e%o$

    1*

    4042 =0775 t = 0325 T.$ o&,e $ t ? x 1 e T ? x 2= y2

    Universidade Federal de Campina Grande - Prof. André Luiz Fiquene de Brito

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    16/28

    I&terpreta@o ,o Mo,e%o$

    i) ara moer&e -o entro -o -elineamento ( x1= e x2 = ), no amino a&e-ente -eer8amo& moer,##! 7ni-a-e& na -ire:0o -e  x1  (temo) ara

    a-a ,32! 7ni-a-e& -e x2 (temerat7ra)

      4042 =0775 t = 0325 T.$

    ii) A -ire:0o -a inlina:0o a&e-ente a&&a elo onto entral ( x1= e  x2  = ) e tem inlina:0o 32!G ##! = 42

    = y2

    1

    Universidade Federal de Campina Grande - Prof André Luiz Fiquene de Brito

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    17/281,

    ///:

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    18/28

     - reião experimental ser/ (0 e 40"min para o tempo de rea+ão e (150 e10")* para temperatura

      A" vari.vei" i$!epe$!e$te" /atore"

     oram #o!ii#a!a" para /(1-+1  para "imp)ii#ar o" #.)#')o"

    16

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    19/28

     1"Como a2ustar o 3odelo de 1rdem, para encontrar a reião

    do timo

     e"po"ta ea)iar $ovo"experime$to" a partir !o 6  

    rea)ia!o

    67E89:E8

    1

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    20/28

      $" Como ;eali

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    21/28

    "Como -c=ar o timo Comodefinir os tempos E as

    temperaturas do experimento 

    e#po#ta$ aar a re%a@o e&trear/+e/# o,/8/a,a# e &atura/#9ou e)per/me&ta/#:.

    (1

    67E89:E8

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    22/28

    ((

    /: Se 1  e 2  repre#e&tam a# ar/+e/#

    &atura/# para tempo e temperaturare#pet/ame&te o# a%ore# #oo,/8/a,o#$

    !1!!2

    2!3!1

    1−

    =−

    =

    ξ ξ  x x

    CTama&o,o Pa##o

    M,/a ,oI&tera%o

    )()(   #ó!igo "em ;at'ra) 

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    23/28

    (4

    ! : D/re@o ,a m+)/ma /&%/&a@oa#e&,e&te$

      ara an-ar (moer&e) -o entro -o-elineamento ( x1=0 e x2=0  no amino -ainlina:0o a&en-ente, -eer8amo& moer

    ,##! 7ni-a-e& na -ire:0o x1 ara a-a ,32!7ni-a-e& na -ire:0o -e x2-

      A##/m a ,/re@o ,a /&%/&a@oa#e&,e&te pa##a pe%o po&to e&tra% (x 1#0

    e x $#0" e tem /&%/&a@o 0325E0775?042

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    24/28

    H re&on&@el elo exerimento -ei-e 7&ar 7m temo-e rea:0o -e ! min7to& omo tamano -o a&&o iniial>

    F#a&,o a re%a@o e&tre 1 e x 1 5 m/&uto# &o tempo ,erea@o orre#po&,e a um /&tera%o 9pa##o: &a

    ar/+e% o,/8/a,a x 1 ,e  x 1#1. 

    O# pa##o# &o am/&o ,a /&%/&a@o a#e&,e&te #o$

     ∆ x1 = 1 /1? pa""o para t

     ∆ x2 = 0,325@0,775  x ∆ x1  ∆ x2 = 0,42 x 1

    ∆ x2 = 0,42 /1? pa""o para T (+

    Universidade Federal de Campina Grande - Prof. André Luiz Fiquene de Brito

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    25/28

      Para o&erter o tama&o ,o# pa##o# o,/8/a,o#9  x 1#1 e x $# 0,$5>0,??5 # 0,4$: para a# u&/,a,e#

    or//&a/# ,e tempo e temperatura re#pet/ame&te a#er u#a,o &o e)per/me&to ,ee-#e u#ar$

      O# po&to# e)per/me&ta/# #o o"t/,o# e a pro,u@opara e#te# po&to# o"#era,o# at ue #e pere"a um

    ,er#/mo &a pro,u@o.

    (*

     ∆ x1 = 1 /1? pa""o para t  ∆ x2 = 0,42 /1? pa""o para T

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    26/28

    !2

    !1

    1

    ξ ξ    ∆=∆

    ∆=∆   x x

    (5

      O# re#u%ta,o# #o mo#tra,o# &a Ta"e%a a

    #eu/r$

    !!(11

    !11

    ==∆

    ∆=∆

     x

     x x

    ξ 

    ξ 

    2!42,2

    !22

    ==∆

    ∆=∆

     x

     x x

    ξ 

    ξ 

    7on0ersão

    #empo #emperat-ra

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    27/28

    (,$stes valores s/o o8tidos e2perimentalmenteem la8orat9rio a artir das variveis ori inais

    ( 6Universidade Federal de Campina Grande - Prof. André Luiz Fiquene de Brito

  • 8/17/2019 Aula4_2 Adição Pto Ct Em 2k_Parte II

    28/28

    (6

    Com o planejamento foi observado umacréscimo na resposta até o passo 1!

    $ntretanto: depois deste passo tm-se umde!rés!imo na resposta;11 e 1"ustado na vizin%ana do ponto

    ; 1 ?6* e ( ?1,*