Numeric reconstruction of cytoskeleton with finite element method and topology optimization method

The importance of mechanical aspects related to cell activity and its environment is becoming more evident due to their influence in stem cell differentiation and in the development of diseases such as atherosclerosis. The mechanical tension homeostasis is related to normal tissue behavior and its lack may be related to the formation of cancer, which shows a higher mechanical tension. Due to the complexity of cellular activity, the application of simplified models may elucidate which factors are really essential and which have a marginal effect. The development of a systematic method to reconstruct the elements involved in the perception of mechanical aspects by the cell may accelerate substantially the validation of these models. This work proposes the development of a routine capable of reconstructing the topology of focal adhesions and the actomyosin portion of the cytoskeleton from the displacement field generated by the cell on a flexible substrate. Another way to think of this problem is to develop an algorithm to reconstruct the forces applied by the cell from the measurements of the substrate displacement, which would be characterized as an inverse problem. For these kind of problems, the Topology Optimization Method (TOM) is suitable to find a solution. TOM is consisted of an iterative application of an optimization method and an analysis method to obtain an optimal distribution of material in a fixed domain. One way to experimentally obtain the substrate displacement is through Traction Force Microscopy (TFM), which also provides the forces applied by the cell. Along with systematically generating the distributions of focal adhesion and actin-myosin for the validation of simplified models, the algorithm also represents a complementary and more phenomenological approach to TFM. As a first approximation, actin fibers and flexible substrate are represented through two-dimensional linear Finite Element Method. Actin contraction is modeled as an initial stress of the FEM elements. Focal adhesions connecting actin and substrate are represented by springs. The algorithm was applied to data obtained from experiments regarding cytoskeletal prestress and micropatterning, comparing the numerical results to the experimental onesCNPqFAPES

Design of graded and homogeneous piezoelectric multi-actuators using the topology optimization method.

Microdispositivos piezelétricos tem uma vasta aplicação em mecânica de precisão, como, por exemplo, manipulação de células, microcirurgias, equipamentos de nanotecnologia e principalmente em microeletromecanismos (MEMS). Os microdispositivos piezelétricos considerados nesta tese essencialmente consistem de uma estrutura multi-flexível atuada por duas ou mais piezocerâmicas, que geram deslocamentos e forças em direções e regiões pré-determinadas do domínio, ou seja, a estrutura multi-flexível atua como um transformador mecânico amplificando e alterando os deslocamentos gerados pelas piezocerâmicas nos movimentos de atuação. O desenvolvimento destes microdispositivos piezelétricos em sua grande maioria não utiliza ferramentas sistemáticas e genéricas. A complexidade dos movimentos de atuação torna o desenvolvimento dos microdispositivos piezelétricos complexo, principalmente devido ao surgimento de movimentos indesejados ou acoplados durante a sua atuação. Portanto, é necessário um método sistemático e eficiente como o método de otimização topológica (MOT), que incorpore na sua formulação as principais exigências de projeto dos microdispositivos, como apresentado nesse trabalho. O MOT implementado é baseado na abordagem CAMD (Distribuição Contínua da Distribuição de Material), onde as pseudo-densidades são interpoladas nos nós de cada elemento finito, resultando numa distribuição contínua de material no domínio. Um método adjunto foi implementado para o cálculo das sensibilidades. São consideradas três formulações. A primeira denominada de MAPs (Multi-Atuadores Piezelétricos) considera as regiões piezocerâmicas fixas, otimizando apenas a estrutura multi-flexível no domínio de projeto. Nesta formulação materiais não-piezelétricos (como, por exemplo, Alumínio) e vazio são distribuídos no domínio de projeto, mantendo as regiões piezocerâmicas fixas e homogêneas. Para validar os resultados obtidos com essa formulação foram fabricados protótipos de nanoposicionadores $XY$, que foram caracterizados experimentalmente utilizando técnicas de interferometria laser, considerando excitação quasi-estática. No entanto, essa primeira formulação impõe restrições no problema, limitando a optimalidade da solução obtida pela otimização topológica. Assim, surgiu a necessidade de desenvolver uma segunda formulação, que permite distribuir simultaneamente material não-piezelétrico, piezelétrico e vazio no domínio de projeto, denominada de LOMPs (Localização Ótima do Material Piezelétrico). A formulação dos LOMPs obtém simultaneamente a localização do material piezelétrico na estrutura flexível otimizada pela OT, e inclui também uma variável de projeto para determinar o ângulo ótimo entre as direções de polarização e do campo elétrico. Nesta formulação como as posições dos eletrodos não são conhecidas, ``a priori\'\', é utilizado como abordagem aplicar um campo elétrico constante para determinar a localização do material piezelétrico e conseqüentemente dos eletrodos. Finalmente, foi explorado o conceito de materiais com gradação funcional (MGFs) no projeto dos MAPs. Os MGFs apresentam uma distribuição contínua de materiais na sua microestrutura, não possuindo interface entre os materiais distribuídos, o que possibilita aumentar a vida útil do dispositivo piezelétrico. Assim, foi implementado uma terceira formulação denominada de MAPs MGFs, que permite obter a gradação ótima de materiais piezelétricos e não-piezelétricos no domínio piezocerâmico dos MAPs, conjuntamente com a topologia da estrutura multi-flexível. Essa formulação foi estendida para projetar atuadores bilaminares MGFs. Todas as formulações desenvolvidas utilizam uma função multi-objetivo, que permite controlar a rigidez e a flexibilidade minimizando o movimento acoplado, de cada movimento de atuação. Os exemplos numéricos são limitados a modelos bi-dimensionais, utilizando o estado plano de tensões e deformações mecânicas e elétricas, uma vez que a grande maioria das aplicações dos microdispositivos piezelétricos são bi-dimensionais.Microtools offer significant promise in a wide range of applications such as cell manipulation, microsurgery, nanotechnology processes, and many other fields. The microtools considered in this doctoral thesis essentially consist of a multi-flexible structure actuated by two or more piezoceramic devices that when each piezoceramic is actuated, it generates an output displacement and force at a specified point of the domain and direction. The multi-flexible structure acts as a mechanical transformer by amplifying and changing the direction of the piezoceramic output displacements. Thus, the development of microtools requires the design of actuated flexible structures that can perform complex movements. The development of these microtools is still in the beginning and it can be strongly enhanced by using design tools. In addition, when multiple piezoceramic devices are involved, coupling effects in their movements become critical, especially the appearance of undesired movements, which makes the design task very complex. One way to avoid such undesirable effects is the use of a systematic design method, such as topology optimization, with appropriate formulation of the optimization problem. The topology optimization method implemented is based on the CAMD (Continuous Approximation of Material Distribution) approach where fictitious densities are interpolated at each finite element, providing a continuum material distribution in the domain. The corresponding sensitivity analysis is presented using the adjoint method. Three formulations are considered. The first formulation, called Piezoelectric Multi-Actuators (PMAs), keeps fixed piezoceramic positions in the design domain and only the flexible structure is designed by distributing some non-piezoelectric material (Aluminum, for example). $XY$ Piezoelectric Nanopositioner are manufactured and experimentally analyzed to validate the results of the topology optimization obtained using this formulation. Experimental analyses are conducted using laser interferometry to measure displacement, while considering a quasi-static excitation. However, this first formulation imposes a constraint to the position of piezoelectric material in the optimization problem limiting the optimality of the solution. Thus, the second formulation presented, called LOMPs, allows the simultaneous distribution of non-piezoelectric and piezoelectric material in the design domain, to achieve certain specified actuation movements. The optimization problem is posed as the simultaneous search for an optimal topology of a flexible structure as well as the optimal position of piezoceramics in the design domain and optimal rotation angle of piezoceramic material axes that maximize output displacements or output forces at a specified point of the domain and direction. When the distribution of a non-piezoelectric conductor material and a piezoceramic material is considered in the design domain, the electrode positions are not known ``a priori\'\'. To circumvent this problem, an electric field is applied as electrical excitation. Finally, the concept of functionally graded materials (FGM) is applied to PMAs design. FGMs are special materials that possess continuously graded properties without interfaces which can increase lifetime of piezoelectric devices. Thus, a third formulation is implemented to find the optimum gradation and polarization sign variation of piezoceramic FGMs, while simultaneously optimizing the multi-flexible structural configuration. This formulation is extended to design bimorph type FGM actuators. For all developed formulations, a multi-objective function is defined that controls the stiffness and flexibility, minimizing the coupling movement of each actuated movement. The present examples are limited to two-dimensional models because most part of the applications for such micro-tools are planar devices

Design of flextensional piezoelectric actuator using the topology optimization method.

Atuadores Piezelétricos Flextensionais consistem de uma estrutura flexível atuada por cerâmicas piezelétricas (ou “pilhas” de cerâmicas). A estrutura flexível conectada a piezocerâmica deve gerar deslocamentos e forças em diferentes pontos específicos do domínio, para uma direção especificada. Estes atuadores são usados em aplicações de mecânica de precisão, tal como, sistemas microeletromecânicos (MEMS), manipulador de células, interferometria laser, equipamentos de nanotecnologia, equipamentos de microcirurgias, nanoposicionadores, sonda de varredura microscópica, e etc. Porém, devido ao fato destes atuadores consistirem principalmente de um mecanismo flexível, seu projeto é complexo. A estrutura flexível comporta-se como um transformador mecânico pela amplificação para converter, direcionar e amplificar os pequenos deslocamentos gerados pela piezocerâmica (ordem de nanômetros). A estrutura flexível é projetada distribuindo-se flexibilidade e rigidez no domínio de projeto, o que pode ser obtido usando a otimização topológica. Portanto, o objetivo deste trabalho é implementar um método sistemático baseado no método de otimização topológica para projetar atuadores piezelétricos flextensionais. Essencialmente, o método de otimização topológica consiste em encontrar a distribuição ótima de material perfurando o domínio de projeto com infinitos microfuros. O material em cada ponto pode alterar de vazio a total presença de material, também assumindo material intermediário (ou compósito). A implementação do método de otimização topológica é baseado no modelo de material SIMP (Simple Isotropic Material with Penalization). O problema de otimização é posto como a maximização dos deslocamentos gerados (ou força de blocagem) em diferentes pontos e direções especificadas do domínio. Considerando o comportamento linear da piezocerâmica. Alterando a flexibilidade e a rigidez da estrutura flexível conectada a piezocerâmica obtém-se diferentes tipos de atuadores piezelétricos flextensionais, que podem ser projetados para determinadas aplicações. Para ilustrar o método, os exemplos mostrados são modelos bidimensionais (2D), uma vez que a maior parte das aplicações envolve dispositivos planos. Estes atuadores são fabricados usando corrosão química em chapas de cobre abaixo de 200 μm de espessura através do método de litografia. Técnica de corrosão química tem um baixo custo e permite-nos fabricar diversos protótipos para testes. Esta técnica pode ser facilmente utilizada no LNLS (Laboratório Nacional de Luz Síncrotron – Campinas). Análise experimental destes protótipos são procedidas para medição de deslocamentos usando uma Probe Station. Como trabalho futuro, estes protótipos serão construídos em escala de MEMS.Flextensional Piezoelectric Actuators consist of a flexible structure actuated by piezoelectric ceramics (or a stack of piezoceramics). The flexible structure connected to the piezoceramic must generate displacements and forces in different specified points of the domain, according to a specific direction. These actuators are applied to precision mechanic applications such as microelectromechanical systems (MEMS), cell manipulators, laser interferometers, nanotechnology equipment, microsurgery equipment, nanopositioners, scanning probe microscopy, etc. However, due to the fact these actuators essentially consist of a compliant mechanism their design is complex. The compliant structure behaves as a mechanical transform by amplifying and changing the direction of small output displacements generated by piezoceramics (order of nanometer). The flexible structure is designed by distributing flexibility and stiffness in the design domain, which can be archieved by using topology optimization. Therefore, the objective of this work is to implement a systematic method based on topology optimization method to design flextensional piezoelectric actuators. Essentially, the topology optimization method consists of finding the optimal material distribution in a perforated design domain with infinite microvoids. The material in each point can change from void to full material, also assuming intermediate (or composite) material. The implemented topology optimization method is based on the SIMP (Simple Isotropic Material with Penalization) material model. The optimization problem is posed as maximization of output displacements (or grabbing forces) in different specified directions and points of the domain. A linear behavior of piezoceramic is considered. By changing the flexibility and stiffness of flexible structure connected to the piezoceramics different types of flextensional piezoelectric actuators can be designed for a desired application. To illustrate the method, examples presented herein are limited to two-dimensional (2D) models once in most part of applications of these actuators they are planar devices. These actuators are manufactured by using chemical corrosion on a 200 um thickness copper plate through lithography method. Chemical corrosion technique has a low cost and it allow us to manufacture several prototypes for testing. For this technique, facilities of the micromachining laboratory of National Sincroton Light Laboratory (LNLS - Campinas) are used. Experimental analysis of these prototypes are conducted by measuring displacements using a probe station. As a future work, these prototypes will be built in a MEMS scale

Design of graded and homogeneous piezoelectric multi-actuators using the topology optimization method.

Microdispositivos piezelétricos tem uma vasta aplicação em mecânica de precisão, como, por exemplo, manipulação de células, microcirurgias, equipamentos de nanotecnologia e principalmente em microeletromecanismos (MEMS). Os microdispositivos piezelétricos considerados nesta tese essencialmente consistem de uma estrutura multi-flexível atuada por duas ou mais piezocerâmicas, que geram deslocamentos e forças em direções e regiões pré-determinadas do domínio, ou seja, a estrutura multi-flexível atua como um transformador mecânico amplificando e alterando os deslocamentos gerados pelas piezocerâmicas nos movimentos de atuação. O desenvolvimento destes microdispositivos piezelétricos em sua grande maioria não utiliza ferramentas sistemáticas e genéricas. A complexidade dos movimentos de atuação torna o desenvolvimento dos microdispositivos piezelétricos complexo, principalmente devido ao surgimento de movimentos indesejados ou acoplados durante a sua atuação. Portanto, é necessário um método sistemático e eficiente como o método de otimização topológica (MOT), que incorpore na sua formulação as principais exigências de projeto dos microdispositivos, como apresentado nesse trabalho. O MOT implementado é baseado na abordagem CAMD (Distribuição Contínua da Distribuição de Material), onde as pseudo-densidades são interpoladas nos nós de cada elemento finito, resultando numa distribuição contínua de material no domínio. Um método adjunto foi implementado para o cálculo das sensibilidades. São consideradas três formulações. A primeira denominada de MAPs (Multi-Atuadores Piezelétricos) considera as regiões piezocerâmicas fixas, otimizando apenas a estrutura multi-flexível no domínio de projeto. Nesta formulação materiais não-piezelétricos (como, por exemplo, Alumínio) e vazio são distribuídos no domínio de projeto, mantendo as regiões piezocerâmicas fixas e homogêneas. Para validar os resultados obtidos com essa formulação foram fabricados protótipos de nanoposicionadores $XY$, que foram caracterizados experimentalmente utilizando técnicas de interferometria laser, considerando excitação quasi-estática. No entanto, essa primeira formulação impõe restrições no problema, limitando a optimalidade da solução obtida pela otimização topológica. Assim, surgiu a necessidade de desenvolver uma segunda formulação, que permite distribuir simultaneamente material não-piezelétrico, piezelétrico e vazio no domínio de projeto, denominada de LOMPs (Localização Ótima do Material Piezelétrico). A formulação dos LOMPs obtém simultaneamente a localização do material piezelétrico na estrutura flexível otimizada pela OT, e inclui também uma variável de projeto para determinar o ângulo ótimo entre as direções de polarização e do campo elétrico. Nesta formulação como as posições dos eletrodos não são conhecidas, ``a priori\'\', é utilizado como abordagem aplicar um campo elétrico constante para determinar a localização do material piezelétrico e conseqüentemente dos eletrodos. Finalmente, foi explorado o conceito de materiais com gradação funcional (MGFs) no projeto dos MAPs. Os MGFs apresentam uma distribuição contínua de materiais na sua microestrutura, não possuindo interface entre os materiais distribuídos, o que possibilita aumentar a vida útil do dispositivo piezelétrico. Assim, foi implementado uma terceira formulação denominada de MAPs MGFs, que permite obter a gradação ótima de materiais piezelétricos e não-piezelétricos no domínio piezocerâmico dos MAPs, conjuntamente com a topologia da estrutura multi-flexível. Essa formulação foi estendida para projetar atuadores bilaminares MGFs. Todas as formulações desenvolvidas utilizam uma função multi-objetivo, que permite controlar a rigidez e a flexibilidade minimizando o movimento acoplado, de cada movimento de atuação. Os exemplos numéricos são limitados a modelos bi-dimensionais, utilizando o estado plano de tensões e deformações mecânicas e elétricas, uma vez que a grande maioria das aplicações dos microdispositivos piezelétricos são bi-dimensionais.Microtools offer significant promise in a wide range of applications such as cell manipulation, microsurgery, nanotechnology processes, and many other fields. The microtools considered in this doctoral thesis essentially consist of a multi-flexible structure actuated by two or more piezoceramic devices that when each piezoceramic is actuated, it generates an output displacement and force at a specified point of the domain and direction. The multi-flexible structure acts as a mechanical transformer by amplifying and changing the direction of the piezoceramic output displacements. Thus, the development of microtools requires the design of actuated flexible structures that can perform complex movements. The development of these microtools is still in the beginning and it can be strongly enhanced by using design tools. In addition, when multiple piezoceramic devices are involved, coupling effects in their movements become critical, especially the appearance of undesired movements, which makes the design task very complex. One way to avoid such undesirable effects is the use of a systematic design method, such as topology optimization, with appropriate formulation of the optimization problem. The topology optimization method implemented is based on the CAMD (Continuous Approximation of Material Distribution) approach where fictitious densities are interpolated at each finite element, providing a continuum material distribution in the domain. The corresponding sensitivity analysis is presented using the adjoint method. Three formulations are considered. The first formulation, called Piezoelectric Multi-Actuators (PMAs), keeps fixed piezoceramic positions in the design domain and only the flexible structure is designed by distributing some non-piezoelectric material (Aluminum, for example). $XY$ Piezoelectric Nanopositioner are manufactured and experimentally analyzed to validate the results of the topology optimization obtained using this formulation. Experimental analyses are conducted using laser interferometry to measure displacement, while considering a quasi-static excitation. However, this first formulation imposes a constraint to the position of piezoelectric material in the optimization problem limiting the optimality of the solution. Thus, the second formulation presented, called LOMPs, allows the simultaneous distribution of non-piezoelectric and piezoelectric material in the design domain, to achieve certain specified actuation movements. The optimization problem is posed as the simultaneous search for an optimal topology of a flexible structure as well as the optimal position of piezoceramics in the design domain and optimal rotation angle of piezoceramic material axes that maximize output displacements or output forces at a specified point of the domain and direction. When the distribution of a non-piezoelectric conductor material and a piezoceramic material is considered in the design domain, the electrode positions are not known ``a priori\'\'. To circumvent this problem, an electric field is applied as electrical excitation. Finally, the concept of functionally graded materials (FGM) is applied to PMAs design. FGMs are special materials that possess continuously graded properties without interfaces which can increase lifetime of piezoelectric devices. Thus, a third formulation is implemented to find the optimum gradation and polarization sign variation of piezoceramic FGMs, while simultaneously optimizing the multi-flexible structural configuration. This formulation is extended to design bimorph type FGM actuators. For all developed formulations, a multi-objective function is defined that controls the stiffness and flexibility, minimizing the coupling movement of each actuated movement. The present examples are limited to two-dimensional models because most part of the applications for such micro-tools are planar devices

Numeric reconstruction of 2D cellular actomyosin network from substrate displacement

Introduction: One of the fundamental structural elements of the cell is the cytoskeleton. Along with myosin, actin microfilaments are responsible for cellular contractions, and their organization may be related to pathological changes in myocardial tissue. Due to the complexity of factors involved, numerical modeling of the cytoskeleton has the potential to contribute to a better understanding of mechanical cues in cellular activities. In this work, a systematic method was developed for the reconstruction of an actomyosin topology based on the displacement exerted by the cell on a flexible substrate. It is an inverse problem which could be considered a phenomenological approach to traction force microscopy (TFM). Methods An actomyosin distribution was found with a topology optimization method (TOM), varying the material density and angle of contraction of each element of the actomyosin domain. The routine was implemented with a linear material model for the bidimensional actomyosin elements and tridimensional substrate. The topology generated minimizes the nodal displacement squared differences between the generated topology and experimental displacement fields obtained by TFM. The structure resulting from TOM was compared to the actin structures observed experimentally with a GFP-attached actin marker. Results The optimized topology reproduced the main features of the experimental actin and its squared displacement differences were 11.24 µm2, 27.5% of the sum of experimental squared nodal displacements (40.87 µm2). Conclusion This approach extends the literature with a model for the actomyosin structure capable of distributing anisotropic material freely, allowing heterogeneous contraction over the cell extension

Integral Piezoactuator System with Optimum Placement of Functionally Graded Material - A Topology Optimization Paradigm

Piezoactuators consist of compliant mechanisms actuated by two or more piezoceramic devices. During the assembling process, such flexible structures are usually bonded to the piezoceramics. The thin bonding layer(s) between the compliant mechanism and the piezoceramic may induce undesirable behavior, including unusual interfacial nonlinearities. This constitutes a drawback of piezoelectric actuators and, in some applications, such as those associated to vibration control and structural health monitoring (e. g., aircraft industry), their use may become either unfeasible or at least limited. A possible solution to this standing problem can be achieved through the functionally graded material concept and consists of developing `integral piezoactuators`, that is those with no bonding layer(s) and whose performance can be improved by tailoring their structural topology and material gradation. Thus, a topology optimization formulation is developed, which allows simultaneous distribution of void and functionally graded piezoelectric materials (including both piezo and non-piezoelectric materials) in the design domain in order to achieve certain specified actuation movements. Two concurrent design problems are considered, that is the optimum design of the piezoceramic property gradation, and the design of the functionally graded structural topology. Two-dimensional piezoactuator designs are investigated because the applications of interest consist of planar devices. Moreover, material gradation is considered in only one direction in order to account for manufacturability issues. To broaden the range of such devices in the field of smart structures, the design of integral Moonie-type functionally graded piezoactuators is provided according to specified performance requirements.FUSP (University of Sao Paulo Foundation)Polytechnic School of University of Sao Paulo (USP, Brazil)[08/51070 - 0]FAPESP (Fundacao de Amparo a Pesquisa do Estado de Sao Paulo)CNPq[s 303689/2009-9]University of Illinois at Urbana-Champaign (UIUC

Design of Multi-Actuated Piezoelectric Mechanisms Using Topology Optimization

Piezoelectric mechanisms o#er significant promise in a wide range of applications such as cell manipulation, microsurgery, nanotechnology processes, and many other fields. The piezoelectric mechanisms considered in this work essentially consist of a multi-flexible structure actuated by two or more piezoceramic devices that must generate di#erent output displacements and forces at di#erent specified points of the domain and in di#erent directions. The multi-flexible structure acts as a mechanical transform by amplifying and changing the direction of the piezoceramic output displacements. However, when multiple piezoceramic devices are involved, coupling e#ects in their movements become critical, especially the appearance of undesired movements, making the design task very complex. One way to avoid undesirable coupling e#ects is to use a systematic design method, such as topology optimization, with appropriate formulation of the optimization problem. In this work, a topology optimization method is developed for the design of multi-actuated piezoelectric mechanisms actuated by multiple piezoceramics, with one aim being to minimize the e#ects of movement coupling. The method is implemented based on the SIMP material model where fictitious densities are interpolated in each finite element, providing a continuum material distribution in the domain. A formulation that allows the simultaneous distribution of conventional and piezoelectric material in the design domain, to achieve certain specified actuation movements is also presented. The optimization algorithm employed is based on sequential linear programming (SLP) and di#erent types of piezoelectric mechanisms required for various applications were designed to demonstrate the usefulness of the proposed method. Examples are ..

Design of pressure vessels using shape optimization: An integrated approach

Previous papers related to the optimization of pressure vessels have considered the optimization of the nozzle independently from the dished end. This approach generates problems such as thickness variation from nozzle to dished end (coupling cylindrical region) and, as a consequence, it reduces the optimality of the final result which may also be influenced by the boundary conditions. Thus, this work discusses shape optimization of axisymmetric pressure vessels considering an integrated approach in which the entire pressure vessel model is used in conjunction with a multi-objective function that aims to minimize the von-Mises mechanical stress from nozzle to head. Representative examples are examined and solutions obtained for the entire vessel considering temperature and pressure loading. It is noteworthy that different shapes from the usual ones are obtained. Even though such different shapes may not be profitable considering present manufacturing processes, they may be competitive for future manufacturing technologies, and contribute to a better understanding of the actual influence of shape in the behavior of pressure vessels. (C) 2011 Elsevier Ltd. All rights reserved.PETROBRAS[0050.0042368.08.4]CNPq[303689/2009-9