Entrar

É necessária uma assinatura da JoVE para visualizar este conteúdo. Faça login ou comece sua avaliação gratuita.

Neste Artigo

  • Resumo
  • Resumo
  • Introdução
  • Protocolo
  • Resultados
  • Discussão
  • Divulgações
  • Agradecimentos
  • Materiais
  • Referências
  • Reimpressões e Permissões

Resumo

Um método de controle semi-ativo ótimo estocástico com compensação de atraso de tempo (SOSC-PSO) é proposto neste artigo, projetado para manter a confiabilidade de estruturas controladas por amortecedores de RM.

Resumo

O uso de amortecedores magnetoreológicos (MR) em sistemas de controle semi-ativos enfrenta um desafio importante: atraso de tempo causado por processos de feedback, o que reduz a confiabilidade das estruturas de engenharia civil sob excitações estocásticas. Este artigo propõe um método de Controle Semi-ativo Ótimo Estocástico com compensação de atraso de tempo (SOSC-PSO), aproveitando a teoria de controle Ótimo Estocástico Físico (PSO) para resolver esse problema e manter a confiabilidade estrutural. O método proposto deriva a força de controle semi-ativa em função dos estados atual e anterior, compensando os atrasos no processo de controle. Para otimizar a eficácia do controle, os principais parâmetros são ajustados com base em um critério de confiabilidade para o sistema. Análises de validação em estruturas de grau único de liberdade e multigrau de liberdade sob excitações sísmicas estocásticas mostram que atrasos de tempo prejudicam significativamente o desempenho dos amortecedores de RM. No entanto, o método SOSC-PSO com compensação de atraso de tempo melhora significativamente a eficácia do controle e, com parâmetros otimizados, aumenta a confiabilidade do sistema de controle estrutural além dos métodos sem otimização de parâmetros.

Introdução

Melhorar o desempenho das estruturas de engenharia diante de eventos catastróficos, como terremotos e ventos extremos, continua sendo uma preocupação primordial na comunidade de engenharia civil. O controle estrutural, uma técnica eficaz para reduzir vibrações, provou aumentar a segurança e a funcionalidade de tais estruturas 1,2,3. Nas últimas décadas, vários métodos e tecnologias avançadas foram desenvolvidos para esse fim. Esses métodos podem ser amplamente classificados em quatro categorias com base no tipo de energia usada para acionar os dispositivos de controle: sistemas de controle ativos, semiativos, passivos e híbridos 4,5,6,7.

No controle ativo, a força de controle necessária é aplicada diretamente através dos dispositivos de controle, o que requer uma quantidade significativa de energia 8,9,10. O controle semiativo, por outro lado, envolve o ajuste das propriedades dos dispositivos de controle (como amortecimento ou rigidez) com base nos sinais de controle, exigindo muito menos energia em comparação com os sistemas ativos11. O controle passivo, por outro lado, depende da dissipação de energia sem qualquer entrada de energia externa no sistema12 , 13 , 14 . Os sistemas híbridos combinam as características das estratégias de controle ativo/semiativo e passivo para alcançar um desempenho mais eficaz15. Dentre essas abordagens, o controle semiativo é considerado particularmente promissor devido ao seu equilíbrio entre baixo consumo de energia e alta eficiência 16,17,18. O amortecedor magnetoreológico (MR), com suas características superiores de amortecimento dinâmico, é considerado um dos dispositivos de controle semi-ativo mais eficazes 19,20,21,22.

No entanto, surge um desafio em sistemas de controle semi-ativos que utilizam amortecedores de RM, pois a lógica de feedback inevitavelmente introduz atrasos de tempo. Esses atrasos são normalmente causados por vários fatores 23,24,25, incluindo: (i) a aquisição e processamento de dados do sensor, (ii) o cálculo da força de controle necessária pelo controlador 26, (iii) transmissão de sinal através de transdutores para os amortecedores MR27 e (iv) a geração de força real pelos próprios amortecedores MR28. Tais atrasos podem levar a discrepâncias entre a força de controle produzida e a força esperada pela estrutura, comprometendo significativamente a eficácia do controle29. Infelizmente, a maioria dos algoritmos de controle existentes para amortecedores de RM não leva em conta esses atrasos.

Além disso, devido à aleatoriedade inerente de eventos catastróficos, qualquer algoritmo de controle semi-ativo eficaz deve ser capaz de manter o desempenho sob excitações estocásticas. O método de controle Gaussiano Quadrático Linear (LQG), uma técnica clássica de otimização estocástica, tem sido explorado para mitigar vibrações em estruturas equipadas com amortecedores MR. Por exemplo, Dyke et al. propuseram uma estratégia de controle ideal de corte LQG para melhorar a resposta sísmica de estruturas equipadas com amortecedores MR, usando feedback de aceleração30. Ying et al. introduziram uma estratégia de controle ótimo estocástico semi-ativo não cortado para estruturas não lineares com amortecedores MR, incorporando técnicas de média estocástica e programação dinâmica31. Outros estudos aplicaram o controle LQG baseado em modal para melhorar o desempenho sísmico de edifícios isolados na base19 e analisaram sua eficácia para mitigação de vibração em torres de turbinas eólicas32. No entanto, os métodos clássicos de LQG, que assumem ruído gaussiano branco, não são adequados para modelar a natureza não estacionária e não gaussiana de excitações do mundo real, como terremotos ou ventos fortes. Para resolver isso, o conceito de controle estocástico ótimo (PSO) baseado em física foi desenvolvido33 , 34 , que supera as limitações dos métodos LQG tradicionais e oferece uma estrutura mais precisa para lidar com diversas excitações estocásticas35 . Estudos mostraram que o controle ótimo estocástico semi-ativo baseado em PSO melhora significativamente a mitigação de vibrações para estruturas lineares e não lineares sujeitas a forças sísmicas estocásticas36.

Assim, há uma necessidade urgente de um método de controle semi-ativo robusto que compense os atrasos, aumentando a confiabilidade das estruturas controladas por amortecedores de RM25,37. Além disso, para garantir o desempenho ideal do controle sob excitações estocásticas, é essencial otimizar os parâmetros críticos do método de controle proposto usando uma abordagem baseada em confiabilidade. Portanto, um método de controle semi-ativo ótimo estocástico com compensação de atraso de tempo (SOSC-PSO) é proposto neste artigo para melhorar a confiabilidade de estruturas com amortecedores de RM.

Protocolo

1. Método de controle semi-ativo ótimo estocástico

Como o efeito de controle de um amortecedor MR é significativamente influenciado pelo inevitável atraso de tempo, um método de controle semi-ativo com um algoritmo de compensação de atraso de tempo é desenvolvido para aumentar o desempenho da estrutura controlada. Além disso, a aleatoriedade inerente às excitações externas causa as incertezas óbvias das respostas dinâmicas. O controle PSO é introduzido para otimizar os parâmetros críticos do método proposto para garantir a confiabilidade do sistema de controle da estrutura.

1.1 Algoritmo de compensação de controle semi-ativo

Sem perda de generalidade, considere uma estrutura n-Degree-Of-Freedoms (n-DOF) controlada por amortecedores MR com atraso de tempo, da qual a equação de movimento submetida a excitações estocásticas é dada por:

figure-protocol-1031 Características (1)

Onde figure-protocol-1192 representam as matrizes de massa, amortecimento e rigidez n x n da estrutura, respectivamente; n representa o número de graus de liberdade da estrutura controlada. Os vetores de deslocamento, velocidade e aceleração da estrutura são representados por figure-protocol-1548, respectivamente. O pequeno e os dois pontos acima dos símbolos representam a primeira e a segunda derivadas. Us (t - tD) representa o vetor r-dimensional da força de controle retardada no tempo dos amortecedores MR e r representa o número de amortecedores; t representa o tempo; tD denota atraso de tempo, assumido como uniforme em todos os amortecedores neste estudo, onde figure-protocol-2064, com aqui o número figure-protocol-2188integral , e figure-protocol-2305 denotando o período de amostragem. figure-protocol-2431 representa o vetor p-dimensional de excitações estocásticas e figure-protocol-2584 representa o vetor de parâmetro estocástico que caracteriza a aleatoriedade associada a excitações externas. A dimensão de figure-protocol-2798 depende do modelo usado para as excitações externas, mas não está relacionada aos graus de liberdade mecânicos da estrutura. Bs e D são matrizes n x r e n x p que especificam as localizações dos amortecedores e excitações externas, respectivamente. Na representação do espaço de estados, a Eq. (1) é escrita como:

figure-protocol-3293 (2)

Onde figure-protocol-3438 representa o vetor de estado 2n-dimensional; figure-protocol-3574 representa a matriz do sistema 2n x 2n; figure-protocol-3719 representa a matriz de localização 2n x r dos amortecedores MR; figure-protocol-3888 Representa a matriz de localização 2n x p de excitações externas. Esses parâmetros são expressos como:

figure-protocol-4121 , figure-protocol-4214 , figure-protocol-4307 figure-protocol-4398 , (3)

Onde figure-protocol-4545 denota uma matriz de identidade com a mesma ordem de figure-protocol-4689. Para facilitar o cálculo, a equação de espaço de estados contínuos Eq. (2) pode ser expressa de forma discreta como:

figure-protocol-4943 Eletrônicos (4)

Onde o ponto de figure-protocol-5111 tempo é simplificado como figure-protocol-5228. E figure-protocol-5322, figure-protocol-5414, figure-protocol-5506 denotam matrizes 2n x 2n, 2n x r e 2n x p, que são expressas como:

figure-protocol-5706, figure-protocol-5798, figure-protocol-5890 (5)

Onde figure-protocol-6035 denota o período de amostragem.

Para alcançar uma eficácia semelhante à do controle ativo, um método de controle simples e eficiente baseado no algoritmo Hrovat 38 é proposto para o controle baseado em amortecedor MR com atraso de tempo:

figure-protocol-6420 Eletrônicos (6)

onde figure-protocol-6579 representa o sinal de força de controle semi-ativo no ponto de tempo\figure-protocol-6764 para o amortecedor MR; figure-protocol-6892 representa a força de comando activa de referência com compensação de atraso de tempo; figure-protocol-7081 representa a força máxima de amortecimento do amortecedor MR, que é mutável; figure-protocol-7249 representa o símbolo de operação de valor absoluto; figure-protocol-7392 representam as forças coulombianas máximas e mínimas do amortecedor MR; figure-protocol-7580 representa o coeficiente de amortecimento viscoso; figure-protocol-7747 representa a entrada de velocidade no amortecedor MR, que é a velocidade de movimento do pistão em relação ao cilindro do amortecedor. Na Eq. (6), figure-protocol-8002 representam os parâmetros dimensionados do amortecedor MR.

A Eq. (6) mostra o cálculo da força de controle semi-ativa do amortecedor MR com atraso de tempo. Vê-se que a força figure-protocol-8313de controle semi-ativa, no intervalo figure-protocol-8440 de tempo na Eq. (4), é calculada com base na força figure-protocol-8599 de controle ativa no intervalo figure-protocol-8721 de tempo e na variável figure-protocol-8835 de estado no intervalo figure-protocol-8949 de tempo desde o atraso dos figure-protocol-9068 intervalos de tempo. Para obter a força figure-protocol-9199de controle ativa, a função de custo convencional é estabelecida por

figure-protocol-9386 (7)

Onde figure-protocol-9531 representa a matriz de ponderação semidefinida positiva simétrica 2n x 2n do estado do sistema; figure-protocol-9718 representa a matriz de ponderação definida positiva simétrica r x r da força de controle. Vê-se que a excitação externa é insignificante. Na verdade, com base na Eq. (7), o sistema de controle da estrutura pode alcançar o efeito de controle ideal, seja qual for o tipo de excitação externa. Como forma discreta, a Eq. (7) é expressa como39:

figure-protocol-10202 (8)

A minimização da função figure-protocol-10383 de custo leva a um problema de valor extremo condicional, e a força figure-protocol-10542 de controle ativa é calculada por25:

figure-protocol-10736 Telefone (9)

Onde figure-protocol-10890 denotam os ganhos de controle para a variável figure-protocol-11027 de estado no intervalo figure-protocol-11141 de tempo e força figure-protocol-11249 de controle ativa no intervalo figure-protocol-11371de tempo , que é determinado pelas matrizes figure-protocol-11505 de ponderação e figure-protocol-1161237. Desde o inevitável atraso de tempo, a força figure-protocol-11774 de controle obtida é atuada na estrutura em pontos figure-protocol-11916de tempo . Portanto, a força de controle ativo na Eq. (9) é calculada por:

figure-protocol-12140 Eletrônicos (10)

Para maior clareza de concepção, a força de controle semi-ativa na Eq. (10) é expressa da seguinte forma com base na Eq. (6):

figure-protocol-12447 Artigos (11)

Onde figure-protocol-12601 denota a força de controle ativa no ponto de figure-protocol-12737tempo . Substituindo a Eq. (11) na Eq. (10),

figure-protocol-12948 Acessórios (12)

NOTA: Todos os valores de estado figure-protocol-13135 do ponto figure-protocol-13235 de tempo até são figure-protocol-13360 necessários para calcular figure-protocol-13477. Portanto, o seguinte método de previsão de estado é introduzido40.

figure-protocol-13702 Acessórios (13)

Onde figure-protocol-13859 denota a variável prevista. Pela iteração da Eq. (13), o estado figure-protocol-14031 em figure-protocol-14142 é calculado. A força ativa de comando Ua(k) é expressa como:

figure-protocol-14347Telefone (14)

A força de controle semi-ativa pode ser obtida combinando as Eqs. (6), (12) e (13).

1.2 Análise de confiabilidade do sistema de controle da estrutura

De acordo com o princípio da preservação da probabilidade, o sistema figure-protocol-14824 aumentado e figure-protocol-14927 a estrutura com amortecedores de RM são conservados em probabilidade e são governados pelas seguintes equações de evolução de densidade de probabilidade generalizada (GDEEs)38:

figure-protocol-15247 Eletrônicos (15)

figure-protocol-15400 Acessórios (16)

Onde figure-protocol-15557, figure-protocol-15666 denotam os componentes correspondentes do estado do sistema e da força de controle semi-ativa com atraso de tempo, respectivamente; figure-protocol-15900 e figure-protocol-15993 denotam as funções de densidade de probabilidade conjunta de sistemas figure-protocol-16154 de amostra aumentada e figure-protocol-16268, respectivamente; figure-protocol-16377 denota a amostra dentro do espaço amostral de figure-protocol-16514; figure-protocol-16606 e figure-protocol-16699) denotam as derivadas parciais de tempo de primeira ordem do estado do sistema e componentes de força de controle semi-ativos, respectivamente. A densidade de probabilidade conjunta funciona figure-protocol-16981 e figure-protocol-17074 pode ser alcançada resolvendo as Eqs. (15) e (16) com as seguintes condições iniciais:

figure-protocol-17314 Acessórios (17)

figure-protocol-17466 Acessórios (18)

Onde figure-protocol-17623 denota a função delta de Dirac; figure-protocol-17746 e figure-protocol-17850 denotam os valores iniciais determinísticos de figure-protocol-17988 e figure-protocol-18081, respectivamente; figure-protocol-18190 denota a função de densidade de probabilidade da amostra figure-protocol-18338.

As funções de densidade de probabilidade instantânea de figure-protocol-18514 e figure-protocol-18607 podem ser obtidas pela integração unidimensional de figure-protocol-18750 e figure-protocol-18843 sobre o domínio do espaço amostral:

figure-protocol-18998 (19)

figure-protocol-19139 (20)

Onde figure-protocol-19285 denota o domínio de distribuição de figure-protocol-19412.

Com base nas Eqs. (19) e (20), vê-se que as informações probabilísticas completas das grandezas físicas em questão são prontamente derivadas, se suas relações associadas aos parâmetros do algoritmo forem definidas. A confiabilidade das grandezas físicas em questão pode ser calculada pela seguinte função de desempenho pertinente:

figure-protocol-19925 Acessórios (21)

figure-protocol-20077 Acessórios (22)

Onde figure-protocol-20234 e figure-protocol-20327 denota a confiabilidade calculada da quantidade de estado em questão e da força de controle semi-ativa, respectivamente; figure-protocol-20539 denota o vetor de valor extremo equivalente da i-ésima quantidade física em questão da estrutura controlada; figure-protocol-20739 denota o vetor de valor extremo equivalente da j-ésima força de controle; o chapéu '-' nos símbolos indica o vetor de valor extremo equivalente41; figure-protocol-21001 denota a i-ésima quantidade física em questão do k-ésimo DOF da estrutura controlada; figure-protocol-21203 denota a j-ésima força de controle; figure-protocol-21341, figure-protocol-21447, figure-protocol-21539,figure-protocol-21630, figure-protocol-21722 e figure-protocol-21826 denota o número de grandezas físicas em questão e DOF da estrutura e amortecedores MR instalados na estrutura, respectivamente. figure-protocol-22056 e figure-protocol-22149 denota os limiares de figure-protocol-22262 e figure-protocol-22355; figure-protocol-22447 denota a probabilidade do evento aleatório.

1.3 Otimização de parâmetros

A razão das amplitudes e figure-protocol-22692 figure-protocol-22783 influencia fortemente a eficácia do controle38 na modalidade de controle de feedback. Portanto, para alcançar a melhor eficácia de controle,figure-protocol-23038 e figure-protocol-23131 precisa ser otimizado. Além disso, como um método de controle baseado em amortecedor de RM, a eficácia do controle também é afetada pelos parâmetros do dispositivo, figure-protocol-23404 na Eq. (6). Para alcançar a confiabilidade ideal do sistema de controle da estrutura, o critério baseado em confiabilidade é proposto na Eq. (25).

figure-protocol-23775
Figura 1: Fluxograma de otimização de parâmetros para controle de estrutura semi-ativa ótima estocástica com compensação de atraso de tempo. Abreviaturas: GDEE = Equação Diferencial Generalizada do Evento; PDF = Função de densidade de probabilidade. Clique aqui para ver uma versão maior desta figura.

Em resumo, são necessárias duas etapas para realizar o controle estocástico da estrutura semiativa ótima com base em amortecedores MR com compensação de atraso de tempo:

1.3.1 Minimizando a função figure-protocol-24525 de custo mostrada na Eq. (7), a lógica de feedback de estado no sentido de amostras é alcançada, ou seja, Eq. 12)

1.3.2 Minimizando a função figure-protocol-24819de desempenho baseada em confiabilidade, os parâmetros ideais no sentido de estatísticas são alcançados. A Figura 1 apresenta o fluxograma de otimização dos parâmetros do controlador na etapa 1.3.2.

As etapas a seguir na etapa 1.3.2, que envolvem duas camadas de loops, são as seguintes:

1.3.2.1 O particionamento do espaço atribuído à probabilidade de excitações estocásticas, que é caracterizado pelo vetor figure-protocol-25374de parâmetro estocástico , permite a identificação de um conjunto de pontos figure-protocol-25540 representativos junto com suas probabilidades atribuídas associadasPq's. Essa abordagem permite a geração eficiente de processos de amostra para a excitação estocástica, denotada como figure-protocol-25826, a ser prontamente obtida.

1.3.2.2 Para a otimização de parâmetros do sistema de controle de estrutura semi-ativa, inicialize ou atualize os valores dos pesos figure-protocol-26104 da função de custo e figure-protocol-26233. A força de controle de realimentação de estado associada com atraso de tempo, consulte a Eq. (12), é então calculada.

1.3.2.3 Obtenção das funções de densidade de probabilidade (PDFs) de respostas estruturais e força de controle resolvendo os GDEEs mostrados nas Eqs. (17)–(20):

A otimização é a conclusão offline, mas não online. Na aplicação do método proposto em sistemas de controle de estrutura real, o ótimo figure-protocol-26863 e figure-protocol-26973 foram alcançados, e a otimização não é necessária no processo real aplicável.

Análise dinâmica determinística do sistema de controle de estrutura semi-ativa submetido à excitação da amostra, pela qual as grandezas figure-protocol-27323 físicas em questão e seus derivados figure-protocol-27461 são calculados.

Empregando um método de diferenças finitas, como o esquema de Lax-Wendroff modificado com propriedades de Variação Total Diminuindo (TVD), as Equações Diferenciais Generalizadas (GDEEs) podem ser resolvidas e soluções numéricas para as funções figure-protocol-27840de densidade de probabilidade conjunta , podem ser derivadas.

Repetindo as duas etapas acima da análise dinâmica determinística e do método das diferenças finitas, e percorrendo todos os pontos figure-protocol-28152representativos, as funções de densidade de probabilidade podem ser obtidas por soma:

figure-protocol-28356 Acessórios (23)

figure-protocol-28508 Eletrônicos (24)

Onde figure-protocol-28666 representa a medida de área do subdomínio associado ao ponto figure-protocol-28832amostral . O método utilizado para resolver a Equação Diferencial Generalizada do Evento (GDEE) é conhecido como método de evolução da densidade de probabilidade (PDEM)42.

1.3.2.4 Os PDFs das grandezas físicas relevantes são usados para avaliar sua confiabilidade, que são então incorporados à função figure-protocol-29275 de desempenho do critério probabilístico.

1.3.2.5 Avaliar se os critérios de terminação para o processo de otimização de parâmetros foram atendidos. Se as condições não forem satisfeitas, retorne à Etapa 1.3.2.2; Se forem, os parâmetros ideais para o sistema de controle retardado podem ser determinados. Neste trabalho, a caixa de ferramentas do algoritmo genético (GA) no MATLAB é utilizada para otimização de parâmetros, o que facilita as atualizações de parâmetros e define as condições de terminação. Por algoritmo genético, os valores ótimos podem ser alcançados em dez iterações, e a convergência é estável sem problemas mínimos locais. Embora o custo computacional do GA seja maior do que a otimização por enxame de partículas ou métodos baseados em gradiente, o GA é bom para lidar com problemas complexos e não diferenciáveis, como o deste manuscrito. Portanto, o GA é aplicado. E como a propriedade de convergência rápida, o custo computacional do GA é aceitável.

1.3.2.6 Para obter um controle de confiabilidade ideal, a seguinte função figure-protocol-30471 de desempenho é formulada. Minimize figure-protocol-30598 para determinar os valores ideais dos parâmetros figure-protocol-30738 e figure-protocol-30831:

figure-protocol-30951} (25)

figure-protocol-31093 Acessórios (26)

Como mencionado anteriormente, as matrizes de ponderação do controlador são simétricas, onde os elementos na matriz figure-protocol-31361 correspondem aos pesos atribuídos ao deslocamento, velocidade e sua interação, enquanto os elementos na matriz figure-protocol-31577 representam os pesos relacionados à força de controle. Com base nos achados de estudos anteriores38, observou-se que a variação das configurações dos pesos figure-protocol-31848 da função de custo e figure-protocol-31991 tem impacto mínimo na eficácia do controle. Consequentemente, adota-se uma configuração simplificada neste estudo, conforme descrito a seguir:

figure-protocol-32253, figure-protocol-32345 (27)

Onde figure-protocol-32505 denota o coeficiente da matriz de ponderação de estado a ser definida; figure-protocol-32681 denota o coeficiente da matriz de ponderação de controle a ser definida; figure-protocol-32845 denota a matriz de identidade. O efeito de controle é influenciado pela proporção de figure-protocol-33021 para figure-protocol-33131, mas não pelos valores em si. Além disso, figure-protocol-33264 na Eq. (6) denota a força Coulombica mínima do amortecedor MR, o que significa que figure-protocol-33455 com a entrada de corrente no amortecedor sendo zero. Então, a Eq. (26) pode ser expressa como:T

figure-protocol-33701 Acessórios (28)

Onde figure-protocol-33858 estão os coeficientes ideais de preocupação. Para torná-los razoáveis em engenharia, as faixas ideais deles são baseadas na capacidade de produção das fábricas de amortecedores MR.

Obviamente, os parâmetros ótimos obtidos pela Eq. (28) mantêm a eficácia de controle ideal do algoritmo de controle semi-ativo proposto com atraso de tempo, o que permite que o sistema de controle da estrutura atinja a confiabilidade ideal equilibrada sujeita a excitações estocásticas.

2. Estudo de caso

Para analisar a eficácia do método de controle semi-ativo proposto com compensação de atraso de tempo para um amortecedor de RM, uma estrutura plana de cisalhamento de um andar presa a um amortecedor de RM foi submetida ao movimento sísmico estocástico horizontal do solo, conforme mostrado na Figura 2. Os parâmetros do sistema de estrutura semi-ativamente controlada são os seguintes: massa figure-protocol-34944da estrutura, frequência figure-protocol-35073circular natural, taxa figure-protocol-35186de amortecimento. De acordo com a experiência passada e as escalas de amortecedores existentes no mercado, os valores limite de deslocamento estrutural, velocidade, aceleração e força de controle são 10 mm, 100 mm/s, 1.500 mm/s2 e 150 kN, respectivamente. Para representar o movimento sísmico estocástico do solo, o modelo de movimento do solo estocástico fisicamente motivado foi empregado43:

figure-protocol-35747 Acessórios (29)

Onde figure-protocol-35904 representa o movimento do solo no domínio do tempo no local de engenharia e figure-protocol-36071 denota o movimento do solo no domínio da frequência no leito rochoso. O vetor figure-protocol-36240 caracteriza a natureza estocástica do movimento do solo na superfície do local de engenharia. Os parâmetros figure-protocol-36439 são variáveis estocásticas que descrevem as características do solo do local, incluindo a frequência figure-protocol-36631 predominante e o amortecimento figure-protocol-36753equivalente. O vetor figure-protocol-36864 representa a natureza estocástica do movimento do solo no leito rochoso, que é influenciado pelas propriedades da fonte e pelo caminho de propagação, indicando figure-protocol-37115 o número de variáveis estocásticas envolvidas neste estágio. figure-protocol-37267 refere-se à frequência circular e i é a unidade imaginária.

A frequência figure-protocol-37459 predominante e a taxa figure-protocol-37572 de amortecimento equivalente do canteiro de obras são parâmetros-chave que caracterizam as propriedades dinâmicas do solo do sítio. A distribuição probabilística e os parâmetros estatísticos dessas grandezas podem ser determinados com base em registros de aceleração sísmica coletados de uma classe específica de locais de engenharia. Para fins ilustrativos, uma classe de local com uma faixa de velocidade de onda de cisalhamento de [150, 250] m / s foi considerada, e o período característico de projeto do movimento do solo foi definido como 0,45 s.

Os parâmetros figure-protocol-38258críticos , figure-protocol-38376 foram tratados como variáveis estocásticas mutuamente independentes, ambos seguindo uma distribuição log-normal para identificação dos parâmetros. Os parâmetros estatísticos foram os seguintes: as médias de , figure-protocol-38693 foram de figure-protocol-3881012 rad/s e 0,1, respectivamente. Os coeficientes de variação de figure-protocol-38981 e figure-protocol-39091 foram de 0,42 e 0,35, respectivamente. O movimento do solo no leito rochoso foi modelado como um processo de ruído branco gaussiano com uma amplitude de Fourier de 0,20 m/s2, correspondendo a um pico de aceleração do solo de 0,11 g. O ângulo de fase usado para gerar o movimento do leito rochoso foi denotado como figure-protocol-39524. Portanto, Sb = 1 e figure-protocol-39646. Essa abordagem pode ser considerada como modelagem de movimentos de solo estocásticos condicionais com uma determinada probabilidade de excedência44. A suposição de distribuição log-normal é amplamente utilizada na modelagem do movimento do solo devido à sua capacidade de capturar a natureza distorcida dos parâmetros de movimento do solo observados.

Para avaliar a sensibilidade dessa suposição, as análises adicionais usando distribuições probabilísticas alternativas, incluindo as distribuições normal e gama, foram realizadas em outras pesquisas45,46. Os resultados indicam que, embora as tendências gerais permaneçam consistentes, a distribuição log-normal fornece o melhor ajuste aos dados observados, particularmente para capturar o comportamento da cauda das medidas de intensidade do movimento do solo

Em virtude do método das esferas tangentes47, uma coleção de 221 pontos representativos com probabilidades atribuídas associadas foi selecionada e as acelerações representativas do solo foram sintetizadas. A frequência de amostragem foi de 50 Hz e a duração dos movimentos do solo foi de 20,48 s. Para atribuir uma intensidade não estacionária ao movimento simulado do solo, uma função de modulação uniforme foi utilizada e sua formulação foi a seguinte42:

figure-protocol-41222 (30)

Onde figure-protocol-41368 tomar 2 e 16 s, respectivamente.

A média e o desvio padrão do movimento sísmico estocástico do solo e um histórico de tempo do movimento sísmico representativo do solo são mostrados na Figura 3. A amplitude da média (0,06 m/s2) foi ~8% da amplitude do desvio padrão (0,8 m/s2), indicando que o modelo de movimento do solo estocástico fisicamente motivado exibiu a propriedade de média zero. Enquanto isso, o movimento sísmico do solo exibiu comportamentos não estacionários notáveis nos domínios temporal e de frequência.

figure-protocol-42105
Figura 2: Esboço de uma estrutura de cisalhamento de um andar com um amortecedor magnetoreológico. figure-protocol-42340 denotam massa da estrutura, taxa de amortecimento e rigidez; figure-protocol-42506 denota a excitação sísmica estocástica. Clique aqui para ver uma versão maior desta figura.

figure-protocol-42853
Figura 3: Estatísticas e amostras representativas de movimentos sísmicos selecionados do solo. (A) Média e desvio padrão do movimento sísmico estocástico do solo; (B) Histórico de tempo do movimento sísmico representativo do solo. Clique aqui para ver uma versão maior desta figura.

2.1 Análise de influência do atraso de tempo

Para analisar a influência do atraso de tempo na eficácia de controle semi-ativo de um amortecedor de RM, a Figura 4 mostra o deslocamento, a velocidade e a aceleração RMS (resposta quadrática média da raiz) da estrutura controlada. A razão da matriz de ponderação e os parâmetros do amortecedor foram figure-protocol-43798, figure-protocol-43890, e figure-protocol-43998, que são os resultados otimizados sem atrasode tempo 38. As respostas da estrutura controlada com atraso de tempo foram maiores do que a estrutura controlada sem atraso de tempo, e a influência do atraso de tempo no efeito de controle representou periodicidade junto com o aumento do atraso de tempo. Como o método de controle semi-ativo com amortecedor MR é um tipo de método de feedback, a periodicidade da influência foi considerada relacionada ao período natural da estrutura controlada (T = 0,56 s).

Em comparação com as respostas da estrutura descontrolada, das quais o deslocamento máximo de MRS, velocidade e aceleração foram de 24,6 mm, 270,0 mm/s e 3111,3 mm/s2, a estrutura controlada por amortecedor MR alcançou um efeito notável. Ao contrário da instabilidade da estrutura controlada ativa com atraso de tempo25, o método de controle semi-ativo baseado em amortecedor MR ainda obteve algum efeito mesmo com atraso de tempo.

figure-protocol-45198
Figura 4: Os quadrados médios máximos das respostas da estrutura controlada com o tempo. (A) Deslocamento; (B) Velocidade; (C) Aceleração. Clique aqui para ver uma versão maior desta figura.

Para analisar a influência do atraso de tempo na confiabilidade de uma estrutura controlada semiativamente, a Figura 5 mostra a confiabilidade do deslocamento, velocidade e aceleração com o aumento do atraso de tempo. A confiabilidade das respostas da estrutura controlada com qualquer atraso de tempo foi menor do que os valores sem atraso de tempo, o que significa que a confiabilidade da estrutura controlada foi diminuída pelo atraso de tempo. Enquanto isso, semelhante ao RMS máximo, o efeito do controle de confiabilidade apresentou periodicidade com o aumento do atraso de tempo. Além disso, a confiabilidade da estrutura controlada, com ou sem atraso de tempo, foi maior do que os valores de estruturas não controladas com confiabilidade de deslocamento, velocidade e aceleração 0,0954, 0,1058 e 0,1111.

figure-protocol-46583
Figura 5: Confiabilidade das respostas da estrutura controlada com atraso de tempo crescente. Clique aqui para ver uma versão maior desta figura.

Em resumo, as respostas da estrutura foram mitigadas pelo método de controle semi-ativo. No entanto, a eficácia do controle, independentemente do RMS ou confiabilidade, do método de controle semi-ativo foi diminuída pelo atraso no tempo. Portanto, é necessário compensar o atraso.

2.1.1 Análise do método de compensação de atraso de tempo

Para analisar a eficácia do método de compensação de atraso de tempo proposto, os históricos de tempo RMS de deslocamento, velocidade e aceleração para sistema não controlado (Unc), de atraso de tempo com não compensação controlada (TDN-SAC-PSO) e de tempo de atraso com compensação controlada (TDC-SAC-PSO) são comparados na Figura 6, onde o atraso de tempo foi de 0,1 s. Os parâmetros da matriz de ponderação e do método de controlo semi-activo são os mesmos que na etapa 2.1.

figure-protocol-47934
Figura 6: Comparação do histórico de tempo quadrático médio de estruturas controladas por diferentes métodos. (A) Deslocamento; (B) Velocidade; (c) aceleração; (D) Força de controle. Abreviaturas: Unc = não controlado; TDC-SAC-PSO = sistema de retardo de tempo com compensação controlada; TDN -SAC-PSO = sistema de retardo de tempo com não compensação controlada. Clique aqui para ver uma versão maior desta figura.

Ambos os métodos de controle diminuíram significativamente as respostas da estrutura em comparação com a estrutura não controlada, revelando as vantagens do método de controle baseado em amortecedores MR. Em comparação com a estrutura não controlada, o RMS máximo de deslocamento, velocidade e aceleração é reduzido em 75,79%, 73,75% e 61,22% com o método TDN-SAC-PSO. Além disso, com a compensação de atraso de tempo, as respostas da estrutura controlada são diminuídas em 82,59%, 80,40% e 73,04% com o método TDC-SAC-PSO. O RMS máximo da força de controle com o método TDC-SAC-PSO diminuiu 8,43% em comparação com o método TDN-SAC-PSO, embora as respostas de vibração do primeiro método tenham sido menores do que o segundo. A análise acima indica que a influência do atraso de tempo no método de controle baseado em amortecedor MR foi efetivamente diminuída pelo método de compensação, que mostra a necessidade de compensação de atraso de tempo para o método de controle semi-ativo do amortecedor MR.

figure-protocol-49716
Figura 7: Comparação em PDF em momentos típicos de respostas de estrutura. (A) Deslocamento; (B) Velocidade; (C) Aceleração. Abreviaturas: PDF = função de densidade de probabilidade; Unc = não controlado; TDC-SAC-PSO = sistema de retardo de tempo com compensação controlada; TDN -SAC-PSO = sistema de retardo de tempo com não compensação controlada. Clique aqui para ver uma versão maior desta figura.

Para revelar de forma abrangente o efeito do método de compensação proposto na incerteza das respostas da estrutura controlada, as comparações da função de densidade de probabilidade (PDF) de deslocamento, velocidade e aceleração em tempos típicos, 3, 7 e 11s, são mostradas na Figura 7. Os PDFs para o método TDN-SAC-PSO foram mais estreitos do que os da estrutura não controlada, o que significa que a incerteza das respostas da estrutura foi diminuída pelo método de controle semi-ativo, mesmo com atraso de tempo. Quando o atraso foi compensado, os PDFs das respostas foram reduzidos ainda mais. Portanto, o método de compensação proposto é necessário para melhorar o efeito de controle semi-ativo.

Método de controleDeslocamentoVelocidadeAceleraçãoForça de controleValor mínimo
Unc0.09540.10580.1111--0.1111
TDC-SAC-PSO0.95650.81070.76540.98450.7654
TDN-SAC-PSO0.77470.5310.20540.8840.2054

Tabela 1: Respostas de vibração e confiabilidade da força de controle de estruturas controladas por diferentes métodos. Abreviaturas: Unc = não controlado; TDC-SAC-PSO = sistema de retardo de tempo com compensação controlada; TDN -SAC-PSO = sistema de retardo de tempo com não compensação controlada.

As respostas de vibração e a confiabilidade da força de controle de estruturas controladas por métodos não controladas, TDC-SAC-PSO e TDN-SAC-PSO são mostradas na Tabela 1. Com o controle do método TDN-SAC-PSO, a confiabilidade da estrutura foi drasticamente aumentada em comparação com a estrutura não controlada, que valida a eficácia do método de controle baseado em amortecedores MR, mesmo com o tempo. No entanto, com a compensação de atraso de tempo, a confiabilidade foi significativamente aumentada em comparação com sem compensação, especialmente para a aceleração. Destaca-se também que a confiabilidade das diferentes respostas e a força de controle apresentaram diferenças significativas para o método TDC-SAC-PSO. Isso indica que os valores ideais de parâmetros críticos para o controle sem atraso de tempo não são ideais para o controle com atraso de tempo. Para obter o efeito de controle ideal das respostas de vibração e da força de controle com influência de atraso de tempo, os parâmetros precisam ser otimizados ainda mais. Com base na análise acima, conclui-se que o método de controle semi-ativo com um amortecedor MR pode efetivamente diminuir as respostas de vibração da estrutura, enquanto os parâmetros precisam ser otimizados devido à influência do atraso de tempo.

2.2 Otimização de parâmetros

Para atingir os valores ótimos dos parâmetros, foi aplicado o método de otimização integrado baseado no critério de confiabilidade38. A otimização dos parâmetros críticos foi realizada com figure-protocol-53786, figure-protocol-53878, figure-protocol-53970, figure-protocol-54062, e os tempos sintonizáveis da força do amortecedor MR s = 8. A caixa de ferramentas GA dentro do MATLAB foi usada para implementar a otimização.

Os resultados de otimização dos parâmetros são figure-protocol-54387, figure-protocol-54493, e figure-protocol-54601 para o atraso figure-protocol-54706de tempo . A Figura 8 mostra os históricos de tempo RMS de deslocamento, velocidade, aceleração e força de controle para as estruturas controladas pelo método Uncontrolled (Unc), TDC-SAC-PSO e SOSC-PSO. O método SOSC-PSO denota o método de controle semi-ativo com compensação de atraso de tempo e valores de parâmetros críticos otimizados.

figure-protocol-55237
Figura 8: Históricos de tempo RMS de deslocamento, velocidade, aceleração e força de controle para estruturas controladas pelo método não controlado e TDC-SAC-PSO e SOSC-PSO. (A) Deslocamento; (B) Velocidade; (c) aceleração; (D) Força de controle. Abreviaturas: Unc = não controlado; TDC-SAC-PSO = sistema de retardo de tempo com compensação controlada; SOSC-PSO = Método de controle semi-ativo ótimo estocástico com compensação de atraso de tempo. Clique aqui para ver uma versão maior desta figura.

A partir da Figura 8, observa-se que os efeitos de controle dos métodos TDC-SAC-PSO e SOSC-PSO apresentaram pouca diferença. O RMS máximo da velocidade de deslocamento e aceleração para o primeiro método foi de 81,60%, 81,21%, 73,62% da estrutura não controlada e 82,59%, 80,40%, 73,04% para o segundo método. Para analisar de forma abrangente as características de probabilidade das respostas para ambos os métodos de controle, a Figura 9 mostra os PDFs de deslocamento, velocidade e aceleração para as estruturas controladas pelo método não controlado, TDC-SAC-PSO e SOSC-PSO. Ambos os métodos de controle tinham quase os mesmos PDFs de deslocamento e velocidade em momentos típicos. Os PDFs de aceleração do método SOSC-PSO foram mais estreitos do que os do método TDC-SAC-PSO, o que significa que a incerteza da aceleração alcançou um melhor efeito de controle pelo método anterior.

figure-protocol-57047
Figura 9: Comparação de PDFs de deslocamento, velocidade e aceleração para diferentes estruturas controladas. (A) Deslocamento; (B) Velocidade; (C) Aceleração. Abreviaturas: PDFs = funções de densidade de probabilidade; Unc = não controlado; TDC-SAC-PSO = sistema de retardo de tempo com compensação controlada; SOSC-PSO = Método de controle semi-ativo ótimo estocástico com compensação de atraso de tempo. Clique aqui para ver uma versão maior desta figura.

A confiabilidade do deslocamento, velocidade, aceleração e força de controle para as estruturas controladas pelos métodos TDC-SAC-PSO e SOSC-PSO não controladas são mostradas na Tabela 2. Ambos os métodos de controle alcançaram aumento significativo na confiabilidade em comparação com a estrutura não controlada. A confiabilidade do deslocamento mostrou pouca diferença entre os métodos TDC-SAC-PSO e SOSC-PSO. Enquanto a confiabilidade da aceleração, a menor confiabilidade da estrutura controlada, alcançou um aumento notável. Isso significa que o método SOSC-PSO atinge o efeito de controle ideal do equilíbrio.

Método de controleDeslocamentoVelocidadeAceleraçãoForça de controleValor mínimo
Unc0.09540.10580.1111--0.1111
TDC-SAC-PSO0.95650.81070.76540.98450.7654
SOSC-PSO0.9530.87040.89340.95860.8934

Tabela 2: Confiabilidade para estruturas controladas por diferentes métodos. Abreviaturas: Unc = não controlado; TDC-SAC-PSO = sistema de retardo de tempo com compensação controlada; SOSC-PSO = Método de controle semi-ativo ótimo estocástico com compensação de atraso de tempo.

A eficácia do método SOSC-PSO proposto com atraso figure-protocol-59491 de tempo foi analisada. Para verificar ainda mais o método proposto para diferentes atrasos de tempo, a Tabela 3 mostra os valores ideais dos parâmetros e a confiabilidade correspondente para o intervalo de figure-protocol-59820atraso.

tD(s)lg(r*)CD (kN · s / mm)U dc,max (kN)DeslocamentoVelocidadeAceleraçãoForça de controle
0-2.5241.44116.1630.99970.99980.97521
0.02-4.4130.414144.0660.99890.99880.8720.9011
0.04-4.4160.25141.7590.99980.98460.88810.9135
0.06-3.2260.203145.6130.9930.93180.87560.9556
0.08-5.2070.345133.6920.96460.89980.89980.953
0.1-4.6690.702130.7390.95310.87040.89340.9586
0.12-4.8950.83369.720.95260.8630.8880.9686
0.14-4.2310.92175.2210.95310.86810.86970.9671
0.16-5.5940.662138.9220.95290.84610.82620.931
0.18-4.2620.88276.3350.95160.78730.76140.9618
0.2-3.9260.78367.9510.9040.72770.71160.9653
0.22-5.9651.07152.5540.83770.64110.62050.965
0.24-4.4811.216128.1480.75670.4740.4740.9687
0.26-3.530.67253.7060.530.40960.41580.9921
0.28-4.4330.63864.2250.35470.37450.34640.9928
0.3-3.5360.45957.7050.27350.3160.27360.9929
Unc------0.09540.10580.1111--

Tabela 3: Os valores ideais dos parâmetros e os valores de confiabilidade correspondentes para diferentes atrasos de tempo. figure-protocol-63210 denota atraso de tempo; figure-protocol-63356 denotar o valor logarítmico ótimo de r,r denota o coeficiente da matriz de ponderação de controle; figure-protocol-63577 denota o coeficiente de amortecimento viscoso do amortecedor MR; figure-protocol-63775 denotam as forças Coulombic máximas e mínimas do amortecedor MR.

Pode-se observar que a confiabilidade do deslocamento, velocidade e aceleração foi melhorada em comparação com a estrutura não controlada, mesmo na presença de atraso de tempo, demonstrando a robustez do método proposto em relação a tais atrasos. No entanto, a confiabilidade das respostas estruturais com o método de controle proposto diminuiu à medida que o atraso aumentava, indicando que, embora o impacto do atraso na eficácia do controle possa ser mitigado, ele não pode ser totalmente eliminado. Notavelmente, a confiabilidade da força de controle permaneceu acima de 90% em diferentes atrasos de tempo.

3. Exemplo numérico

Para verificar o método SOSC-PSO proposto na estrutura MDOF, foi analisada uma estrutura de seis pavimentos com dois amortecedores MR instalados no primeiro e terceiro pavimentos, conforme mostrado na Figura 10. As amostras sísmicas produzidas pelo modelo sísmico estocástico na Seção 2 foram aplicadas, e os limiares de confiabilidade de deslocamento, velocidade, aceleração e força de controle foram de 20 mm, 200 mm/s, 3.000 mm/s2 e 200 kN. Eqs. (27) foram utilizados como matriz de ponderação com figure-protocol-65220. A caixa de ferramentas GA do MATLAB foi usada para otimizar figure-protocol-65372,figure-protocol-65463, e figure-protocol-65557, com intervalos figure-protocol-65678ótimos , figure-protocol-65791, e figure-protocol-65885.

figure-protocol-66048
Figura 10: Estrutura de seis andares com dois amortecedores MR. Abreviatura: MR = magnetorreológico. Clique aqui para ver uma versão maior desta figura.

3.1 Análise de influência do atraso de tempo

Para analisar a influência do atraso de tempo na eficácia de controle do amortecedor MR na estrutura MDOF, a Figura 11 mostra o RMS máximo de deslocamento (deslocamento entre andares), velocidade (velocidade entre andares) e aceleração (aceleração de andar) junto com o atraso de tempo para o método TDC-SAC-PSO. Os parâmetros foram figure-protocol-66863,figure-protocol-66968, e figure-protocol-67076 para o amortecedor MR no 1º andar; figure-protocol-67202 e figure-protocol-67295 para o amortecedor MR no 3º andar, que são os resultados otimizados sem atraso de tempo. O RMS máximo das respostas da estrutura controlada exibiu flutuação do período junto com o atraso de tempo. O RMS máximo de deslocamento e velocidade estava no 1º andar, a aceleração máxima no 6º andar. Além disso, todas as respostas máximas de MRS com qualquer atraso de tempo foram maiores do que os valores sem atraso de tempo, o que significa que o atraso de tempo diminui o efeito de controle do amortecedor de RM.

figure-protocol-67980
Figura 11: RMS máximo de respostas junto com atraso de tempo para a estrutura controlada pelo método TDC-SAC-PSO. (A) Deslocamento; (B) Velocidade; (C) Aceleração. Abreviaturas: RMS = raiz quadrada média; TDC-SAC-PSO = sistema de retardo de tempo com compensação controlada. Clique aqui para ver uma versão maior desta figura.

A confiabilidade do deslocamento, velocidade e aceleração, juntamente com o atraso de tempo, é mostrada na Figura 12. A confiabilidade da velocidade e da aceleração mostra flutuação do período junto com o atraso de tempo, enquanto a confiabilidade do deslocamento diminuiu drasticamente junto com o atraso de tempo; Todos os valores de confiabilidade em qualquer atraso de tempo foram menores do que aqueles sem atraso de tempo. Portanto, o método de compensação de atraso de tempo é necessário para reduzir a influência do atraso de tempo na confiabilidade da estrutura controlada pelo amortecedor MR.

figure-protocol-69288
Figura 12: Confiabilidade das respostas da estrutura controlada junto com o atraso de tempo. Clique aqui para ver uma versão maior desta figura.

3.2 Efeito de controle da compensação de atraso de tempo

Com base no método TDC-SAC-PSO proposto, foi analisada a estrutura de 6 graus de profundidade com dois amortecedores de RM. A Figura 13 mostra os históricos de tempo RMS de deslocamento e aceleração no 1º e 3º andares para as estruturas controladas pelo método TDC-SAC-PSO (Unc), controladas pelo método TDC-SAC-PSO e controladas pelo método TDN-SAC-PSO. O atraso de tempo foi figure-protocol-70151, e a razão figure-protocol-70253da matriz de ponderação; os parâmetros do amortecedor de RM foram figure-protocol-70423, figure-protocol-70529 para o 1º andar; figure-protocol-70637 e figure-protocol-70730 para o 3º andar.

figure-protocol-70923
Figura 13: Históricos de tempo RMS de deslocamento e aceleração no 1º e 3º andares para estruturas controladas por diferentes métodos. (A) Deslocamento no 1º andar; (B) Deslocamento no 3º andar; (C) Aceleração no 1º andar; (D) Aceleração no 3º andar. Abreviaturas: RMS = raiz quadrada média; Unc = não controlado; TDC-SAC-PSO = sistema de retardo de tempo com compensação controlada; TDN -SAC-PSO = sistema de retardo de tempo com não compensação controlada. Clique aqui para ver uma versão maior desta figura.

As RMS máximas de deslocamento no 1º e 3º andares foram reduzidas em 35,42% e 30,44% pelo método TDN-SAC-PSO, onde o atraso de tempo não foi compensado. Já com a compensação, o RMS máximo de deslocamento no 1º e 3º andares foi reduzido em 49,33% e 53,39%, respectivamente. Sem compensação, o RMS máximo de aceleração no 1º e 3º andares aumentou 16,22% e 2,88%, mas diminuiu 25,77% e 36,00%, respectivamente, com compensação. Portanto, a compensação do atraso de tempo no método de controle do amortecedor MR é necessária para reduzir as respostas da estrutura, especialmente para a aceleração.

figure-protocol-72393
Figura 14: PDFs de deslocamento e aceleração no 1º e 3º andares para diferentes métodos de controle. (A) Deslocamento no 1º andar; (B) Deslocamento no 3º andar; (C) Aceleração no 1º andar; (D) Aceleração no 3º andar. Abreviaturas: PDFs = funções de densidade de probabilidade; Unc = não controlado; TDC-SAC-PSO = sistema de retardo de tempo com compensação controlada; TDN -SAC-PSO = sistema de retardo de tempo com não compensação controlada. Clique aqui para ver uma versão maior desta figura.

Os PDFs de deslocamento e aceleração no 1º e 3º pavimentos são mostrados na Figura 14, revelando os diferentes efeitos de controle dos métodos TDC-SAC-PSO e TDN-SAC-PSO em comparação com a estrutura não controlada. Sem compensação, os PDFs de deslocamento no 1º e 3º andares ainda eram mais estreitos do que os da estrutura não controlada. Além disso, com a compensação do atraso de tempo, os PDFs de deslocamento tornaram-se mais estreitos do que os da estrutura controlada pelo método TDN-SAC-PSO. Ao contrário do deslocamento, sem compensação, os PDFs de aceleração no 1º pavimento para o método TDN-SAC-PSO foram mais largos do que os da estrutura não controlada. No entanto, com compensação, os PDFs de aceleração no 1º e 3º andares foram mais estreitos do que os da estrutura não controlada. Assim, verifica-se a necessidade de compensação de atraso de tempo para aumentar a certeza das respostas da estrutura, especialmente para a aceleração.

Método de controleDeslocamentoVelocidadeAceleraçãoForça de controle no andarForça de controle no andar
Unc0.01140.58220.6372----
TDN-SAC-PSO0.58510.52580.39811
TDC-SAC-PSO0.8060.86740.803711

Tabela 4: Confiabilidade para estruturas controladas por diferentes métodos. Abreviaturas: Unc = não controlado; TDC-SAC-PSO = sistema de retardo de tempo com compensação controlada; TDN -SAC-PSO = sistema de retardo de tempo com não compensação controlada.

Os valores de confiabilidade de deslocamento, velocidade, aceleração e força de controle com os sistemas de estrutura controlados pelos diferentes métodos são mostrados na Tabela 4. Os valores de confiabilidade das respostas da estrutura foram calculados usando o método dos valores extremos equivalentes25. Pelo controle do método TDN-SAC-PSO, a confiabilidade do deslocamento aumentou em comparação com a estrutura não controlada, enquanto a confiabilidade da velocidade e aceleração diminuiu. Com a compensação de atraso de tempo no método TDC-SAC-PSO, a confiabilidade do deslocamento, velocidade e aceleração aumentou significativamente. Além disso, os valores de confiabilidade da força de controle para os métodos TDN-SAC-PSO e TDN-SAC-PSO são suficientemente mantidos. Assim, o método de controle de compensação de atraso de tempo proposto alcança um efeito de controle favorável para a estrutura controlada por amortecedor MR.

3.3 Otimização de parâmetros

Para alcançar o efeito de controle ideal, os parâmetros do método de controle proposto foram otimizados com base no critério de confiabilidade. Com o atraso figure-protocol-76417de tempo, o resultado ótimo foi figure-protocol-76553; figure-protocol-76645, e figure-protocol-76753 para o amortecedor MR do 1º andar; figure-protocol-76893 e figure-protocol-77000 para o amortecedor MR do 3º andar.

Os históricos de tempo RMS de deslocamento e aceleração no 1º e 3º andares são mostrados na Figura 15 , onde os efeitos de controle dos métodos TDC-SAC-PSO e SOSC-PSO e a estrutura não controlada são comparados. Com o controle do amortecedor MR, com ou sem compensação de atraso de tempo, o RMS de deslocamento e aceleração diminuiu significativamente. Com o método SOSC-PSO, o RMS máximo de deslocamento no 1º e 3º andares diminuiu 65,15% e 63,16% em relação à estrutura não controlada, respectivamente. Em comparação com o método TDC-SAC-PSO, a resposta de deslocamento foi ainda mais reduzida. O RMS máximo de aceleração no 1º e 3º andares é reduzido em 23,39% e 35,60% para o método SOSC-PSO. Houve pouca diferença no efeito de controle da aceleração para os métodos SOSC-PSO e TDC-SAC-PSO, indicando que o efeito de controle da compensação do atraso de tempo foi diferente para deslocamento e aceleração.

figure-protocol-78165
Figura 15: Históricos de tempo RMS de deslocamento e aceleração no 1º e 3º andares para diferentes métodos de controle. (A) Deslocamento no 1º andar; (B) Deslocamento no 3º andar; (C) Aceleração no 1º andar; (D) Aceleração no 3º andar. Abreviaturas: RMS = raiz quadrada média; Unc = não controlado; TDC-SAC-PSO = sistema de retardo de tempo com compensação controlada; SOSC-PSO = Método de controle semi-ativo ótimo estocástico com compensação de atraso de tempo. Clique aqui para ver uma versão maior desta figura.

Os PDFs de deslocamento e aceleração no 1º e 3º andares são apresentados na Figura 16; os PDFs de deslocamento e aceleração foram reduzidos pelos métodos TDC-SAC-PSO e SOSC-PSO em comparação com a estrutura não controlada. Além disso, os PDFs de deslocamento e aceleração para o método SOSC-PSO foram ainda mais reduzidos em comparação com os da estrutura controlada pelo método TDC-SAC-PSO. Portanto, com a otimização de parâmetros, o método de controle de compensação de atraso de tempo proposto alcançou um efeito de controle melhor do que sem otimização.

figure-protocol-79678
Figura 16: PDFs de deslocamento e aceleração no 1º e 3º andares para diferentes métodos de controle. (A) Deslocamento no 1º andar; (B) Deslocamento no 3º andar; (C) Aceleração no 1º andar; (D) Aceleração no 3º andar. Abreviaturas: PDFs = funções de densidade de probabilidade; Unc = não controlado; TDC-SAC-PSO = sistema de retardo de tempo com compensação controlada; SOSC-PSO = Método de controle semi-ativo ótimo estocástico com compensação de atraso de tempo Clique aqui para ver uma versão maior desta figura.

Os valores de confiabilidade de deslocamento, velocidade, aceleração e força de controle são mostrados na Tabela 5. Os valores de confiabilidade para o método SOSC-PSO foram maiores do que para o método TDC-SAC-PSO. Enquanto isso, a força de controle ainda mantinha confiabilidade suficiente. Assim, a otimização dos parâmetros é necessária para obter o melhor efeito de controle para o método de controle baseado em amortecedor MR.

Método de controleDeslocamentoVelocidadeAceleraçãoForça de controle no andarForça de controle no andar
Unc0.01140.58220.6372----
TDC-SAC-PSO0.8060.86740.803711
SOSC-PSO0.95440.93140.878811

Tabela 5: Confiabilidade para estruturas controladas por diferentes métodos. Abreviaturas: Unc = não controlado; TDC-SAC-PSO = sistema de retardo de tempo com compensação controlada; SOSC-PSO = Método de controle semi-ativo ótimo estocástico com compensação de atraso de tempo.

Resultados

Visando a influência do atraso de tempo na eficácia de controle do método baseado em amortecedores de RM, um método de controle semi-ativo com compensação de atraso de tempo é proposto neste artigo. No método proposto, os parâmetros críticos são otimizados com base no critério de confiabilidade. Ao comparar a eficácia do controle, as seguintes conclusões são abordadas:

(1) O método de controle semi-ativo baseado em um amortecedor MR é mais rob...

Discussão

Com a introdução da teoria de controle ótimo estocástico físico (PSO), um método de controle semi-ativo ótimo estocástico com compensação de atraso de tempo (SOSC-PSO) é proposto neste artigo, projetado para manter a confiabilidade das estruturas controladas por amortecedores de RM. Para compensar o atraso de tempo no método proposto, a força de controle semi-ativa é derivada como a função não apenas dos estados presentes, mas também dos estados e forças de controle em...

Divulgações

Todos os autores não têm conflitos de interesse a declarar.

Agradecimentos

Os autores agradecem o apoio da Fundação de Ciências Naturais da Província de Hebei (Concessão nº E2023210007).

Materiais

NameCompanyCatalog NumberComments
MATLABMathWorks2016Using for the calculation of the proposed method

Referências

  1. Yao, J. T. P. Concept of structural control. J Struct Div. 98 (7), 1567-1574 (1972).
  2. Liu, J., Silva, C. E., Dyke, S. J., Wu, Y., Liu, H. Using real-time hybrid simulation for active mass damper experimentation and validation. Mech Mach Theory. 191, 105474 (2024).
  3. Pan, H., Li, C., Cao, L. Bidirectional pendulum-type tuned tendem mass dampers-inerters for integrated control of alongwind and crosswind responses in super-tall buildings. J Build Eng. 100, 111676 (2025).
  4. Housner, G. W. et al. Structural control: past, present, and future. J Eng Mech. 123 (9), 897-971 (1997).
  5. Rahman, M., Ong, Z. C., Chong, W. T., Julai, S. Smart Semi-active PID-ACO control strategy for tower vibration reduction in Wind Turbines with MR damper. Earthq Eng Eng Vib. 18 (4), 887-902 (2019).
  6. Morales-Beltran, M., Paul, J. Active and semi-active strategies to control building structures under large earthquake motion. J Earthquake Eng. 19 (7), 1086-1111 (2015).
  7. Shu, Z., You, R., Xie, Y. Viscoelastic dampers for vibration control of building structures: A state-of-art review. J Earthq Eng. 28 (12), 3558-3585 (2024).
  8. Soong, T. T. Active structural control: Theory and practice. Longman Scientific & Technical. New York (1990).
  9. Wan, H. P., Ma, Q., Dong, G. S., Luo, Y., Ni, Y. Q. Data-driven model reduction approach for active vibration control of cable-strut structures. Eng Struct. 302, 117434 (2024).
  10. Cao, L., Li, X., Huang, Y., Li, C., Pan, H. High robust eddy current tuned tandem mass dampers-inerters for structures under the ground acceleration. Soil Dyn Earthq Eng. 188, 109040 (2025).
  11. Chha, H., Peng, Y. Adaptive semiactive control of structure with magnetorheological dampers using wavelet packet transform. Adv Struct Eng. 27 (9), 1509-1527 (2024).
  12. Zeng, X., Peng, Y., Chen, J. Serviceability-based damping optimization of randomly wind-excited high-rise buildings. Struct Des Tall Spec. 26 (11), e1371 (2017).
  13. Li, C., Chang, K., Cao, L., Huang, Y. Performance of a nonlinear hybrid base isolation system under the ground motions. Soil Dyn Earthq Eng. 143, 106589 (2021).
  14. Cao, L., Li, C. A high performance hybrid passive base-isolated system. Struct Control Health Monitoring. 29 (3), e2887 (2022).
  15. Zareie, S., Issa, A. S., Seethaler, R., Zabihollah, A., Ahmad, R. A novel SMA-magnetorheological hybrid bracing system for seismic control. Eng Struct. 244, 112709 (2021).
  16. Chu, S. Y., Soong, T. T., Reinhorn, A. M. Active, hybrid and semi-active structural control. John Wiley & Sons, New York (2005).
  17. Danielian, S. A., Jabbari, F., Zareian, F. Effects and distribution of semi-active resettable springs on the seismic response of asymmetric structures. Earthq Eng Struct D. 52 (10). 2966-2982 (2023).
  18. Shuang, Z., Wenliuhan, H. S., Liu, Y. H., Inoue, N., Zhai, Z. P. Seismic response control for bridge piers with semi-active MR damper based on displacement feedback. J Earthq Eng. 28 (3). 849-865 (2024).
  19. Wang, Y., Dyke, S. Modal-based LQG for smart base isolation system design in seismic response control. Struct Control Health Monitoring. 20 (5), 753-768 (2013).
  20. Desai, R. M. et al. Performance evaluation of a single sensor control scheme using a twin-tube MR damper based semi-active suspension. J Vib Eng Technol. 9, 1193-1210 (2021).
  21. Abdi, M. S., Nekooei, M. Jafari, M. A. Seismic control of multi-degrees-of-freedom structures by vertical mass isolation method using MR dampers. Earthq Eng Eng Vib. 23 (2), 503-510 (2024).
  22. Bhowmik, K., Debnath, N. Semi-active vibration control of soft-storey building with magnetorheological damper under seismic excitation. J Vib Eng Technol. 12 (4), 6943-6961 (2024).
  23. Yang, J. N., Akbarpour, A., Askar, G. Effect of time delay on control of seismic-excited buildings. J Struct Eng. 116 (10), 2801-2814 (1990).
  24. Temimi, H., Ben-Romdhane, M., El-Borgi, S., Cha, Y. Time-delay effects on controlled seismically excited linear and nonlinear structures. Int J Struct Stab Dyn. 16 (7), 1550031 (2016).
  25. Peng, J., Wang, L., Zhao, Y., Lenci, S. Time-delay dynamics of the MR damper-cable system with one-to-one internal resonances. Nonlinear Dynamics. 105 (2), 1343-1356 (2021).
  26. Cha, Y. J., Agrawal, A. K., Dyke, S. J. Time delay effects on large-scale MR damper based semi-active control strategies. Smart Mater Struct. 22 (1), 015011 (2013).
  27. Yu, Y. et al. Experimental study of wireless structural vibration control considering different time delays. Smart Mater Struct. 24, 045005 (2015).
  28. Ahmadizadeh, M., Mosqueda, G., Reinhorn, A. M. Compensation of actuator delay and dynamics for real-time hybrid structural simulation. Earthq Eng Struct Dyn. 37 (1), 21-42 (2008).
  29. Bathaei, A., Zahrai, S. M. Compensating time delay in semi-active control of a SDOF structure with MR damper using predictive control. Struc Eng Mech. 82 (4), 445-458 (2022).
  30. Dyke, S. J., Spencer, Jr, B. F., Sain, M. K., Carlson, J. D. Modeling and control of magnetorheological dampers for seismic response reduction. Smart Mater Struct. 5 (5), 565-575 (1996).
  31. Ying, Z. G., Ni, Y. Q., Ko, J. M. A semi-active stochastic optimal control strategy for nonlinear structural systems with MR dampers. Smart Struct Syst. 5 (1), 69-79 (2009).
  32. Rosol, M., Martynowicz, P. Implementation of the LQG controller for a wind turbine tower-nacelle model with a MR tuned vibration absorber. J Theor App Mech-Pol. 54 (4), 1109-1123 (2016).
  33. Lyu, M. Z., Feng, D. C., Chen, J. B., Li, J. A decoupled approach for determination of the joint probability density function of a high-dimensional nonlinear stochastic dynamical system via the probability density evolution method. Comput Method Appl M. 418 (A), 116443 (2024).
  34. Chha, H., Peng, Y. Multiscale stochastic optimal control of hysteretic structures based on wavelet transform and probability density evolution method. Eng Computation. 41 (4), 865-896 (2024).
  35. Li, J., Peng, Y. B., Chen, J. B. A physical approach to structural stochastic optimal controls. Probabilist Eng Mech. 25 (1), 127-141 (2010).
  36. Peng, Y. B., Yang, J. G., Li, J. Seismic risk-based stochastic optimal control of structures using magnetorheological dampers. Nat Hazards Rev. 18 (1), B4016001 (2017).
  37. Dong, X. M., Yu, M., Li, Z., Liao, C., Chen, W. Neural network compensation of semi-active Control for magneto-rheological suspension with time delay uncertainty. Smart Mater Struct. 18 (1), 015014 (2008).
  38. Peng, Y. B., Zhang, Z. Optimal MR damper-based semiactive control scheme for strengthening seismic capacity and structural reliability. J Eng Mech. 146 (6), 04020045 (2020).
  39. Luo, M. X., Gao, M. Cai, G. Delayed full-state feedback control of airfoil flutter using sliding mode control method. J Fluid Struct. 61, 262-273 (2016).
  40. Xu, J., Chung, K. W., Chan, C. L. An efficient method for studying weak resonant double Hopf bifurcation in nonlinear systems with delayed feedbacks. Siam J Appl Dyn Syst. 6 (1), 29-60 (2007).
  41. Yang, J. S., Jensen, H., Chen, J. B. Structural optimization under dynamic reliability constraints utilizing probability density evolution method and metamodels in augmented input space. Struct Multidiscip Optimization. 65, 107 (2022).
  42. Li, J., Chen, J. B. Stochastic dynamics of structures. John Wiley & Sons, Singapore (2009).
  43. Peng, Y. B., Li, J. Stochastic optimal control of structures. Springer (2019).
  44. Dowrick, C. Beyond depression: A new approach to understanding and management. Oxford University Press (2009).
  45. Wang, D., Li, J. Physical random function model of ground motions for engineering purposes. Sci China Technol Sci. 54 (1), 175-182 (2011).
  46. Wang, D., Li, J. A random physical model of seismic ground motion field on local engineering site. Sci China Technol Sci. 55 (7), 2057-2065 (2012).
  47. Chen, J. B., Li, J. Strategy for selecting representative points via tangent spheres in the probability density evolution method. Int J Numer Meth Eng. 74 (13), 1988-2014 (2008).

Reimpressões e Permissões

Solicitar permissão para reutilizar o texto ou figuras deste artigo JoVE

Solicitar Permissão

Explore Mais Artigos

Este m s em JoVEEdi o 218controle semi ativoamortecedor MRatraso de tempocrit rio de confiabilidadeexcita es estoc sticas

This article has been published

Video Coming Soon

JoVE Logo

Privacidade

Termos de uso

Políticas

Pesquisa

Educação

SOBRE A JoVE

Copyright © 2025 MyJoVE Corporation. Todos os direitos reservados