[pt] O agronegócio é um importante setor da economia brasileira pelo incremento que proporciona ao PIB nacional. O País exporta vários produtos agropecuários, sendo líder mundial na produção e venda de soja. O mérito deste sucesso é compartilhado pelos diversos elos que formam sua cadeia de produção e distribuição, elevando a agricultura a patamares de competitividade internacional. Apesar disso, o Brasil apresenta escassa infraestrutura logística, principalmente de armazenagem de grãos em fazendas, que representam elevados investimentos. Somado a um modal rodoviário precário, perdas e desperdícios ocorrem em todo o processo, elevando os custos e dificultando a competição no cenário global. Em vista disso, muitos produtores rurais utilizam serviços de cooperativas. Observa-se que tais serviços são muitas vezes excessivamente custosos ou de reduzida transparência. Assim, a despeito dos elevados recursos envolvidos, e da necessidade de integração de novas atividades no processo da pós-colheita, produtores vêm considerando investir em unidades próprias de armazenagem. Foi realizado um estudo de caso em duas fazendas no Mato Grosso do Sul onde se comparou, pelo método de Custeio Baseado em Atividades (ABC), o cenário atual com uma nova estrutura de integração das etapas pós-colheita. Ademais, o referencial teórico proporcionou um levantamento de ganhos potenciais na comercialização. O método inclui ainda uma avaliação estratégico-mercadológica de outros ganhos complementares. Identificou-se que o investimento em unidade armazenadora de grãos própria gera ganhos potenciais sob as formas de redução de custos, maior liberdade nas decisões da pós-colheita, principalmente na comercialização das safras e um potencial incremento de renda com novos serviços. / [en] Agribusiness is an important sector in the Brazilian economy due to its contribution for the gross national domestic product growth. The country exports several agricultural products and is the world leader in soybean sales. The merit of this success is shared by the various members in the supply chain, which lead the local agriculture to high international levels of competitiveness. However, despite these advances, Brazil has poor infrastructure, especially in farms grain storage facilities. This bottleneck, coupled with a poor road system, result in losses and wastes that occur throughout the production process, increasing costs and reducing competitiveness in a global level. As a result, many farmers hire cooperative services. It is observed that such services are often too costly or of limited transparency. Thus, despite the huge resources involved, and the need for integration of new activities in the post-harvest process, producers consider investing in own storage units. A case study was conducted in the state of Mato Grosso do Sul, using the Activity Based Costing (ABC) as method to compare the current situation of two farms using cooperative post-harvest services with a new structure integrating the post-harvest stages. A literature review on conceptual marketing strategies revealed potential gains with this integration process. The method also includes a strategic-marketing evaluation of other additional gains with the integration. It is concluded that investment in own grain storage units creates potential gains for farmers in the forms of cost savings, freedom for post-harvest sales decisions, and an increase in income with new services.

[pt] Muitos problemas reais em modelagem computacional requerem o uso de aproximação de funções. Em alguns casos a função a ser avaliada no computador é muito complexa, portanto seria desejável que ela fosse substituída por uma função mais simples e mais eficiente de ser calculada. Para fazer isso, calcula-se o valor da função escalar f em um conjunto de N pontos {x1, x2, . . . , XN}, onde x(i) (pertence a) R(n), e faz-se uma estimativa dos valores dessa função f em qualquer outro ponto através de um método de interpolação. Um método de interpolação é qualquer procedimento que toma um conjunto de restrições e determina uma boa função que satisfaça essas condições. O método de interpolação de Shepard originalmente calcula o valor estimado dessa função num ponto qualquer x (pertence a) R(N) como uma média ponderada dos valores da função original nas N amostras dadas. Sendo que o peso para cada amostra x(i) é função das potências negativas das distâncias euclidianas entre os pontos x e x(i). Os núcleos K: R(N) × R(N) (EM) R são funções que correspondem ao produto interno no espaço de Hilbert F da imagem dos pontos x e z por uma função phi (conjunto vazio) : R(N) (EM) F, ou seja K(x, z) = < phi (conjunto vazio) (x), phi (conjunto vazio) (z) >. Na prática, as funções núcleos representam implicitamente o mapeamento feito pela função phi (conjunto vazio) , ou seja, se define qual núcleo usar e não qual phi (conjunto vazio) usar. Esse trabalho propõe uma modificação do método de interpolação de Shepard que é uma simples substituição no método original: ao invés de usar a distância euclidiana entre os pontos x e xi sugere-se usar a distância entre as imagens dos pontos x e x(I) por phi (conjunto vazio) no espaço de Hilbert F, que pode ser calculada diretamente com o uso da função núcleo K. Os resultados mostram que essa pequena modificação gera resultados melhores quando comparados com o método de Shepard original. / [en] Several real problem in computational modeling require function approximations. In some cases, the function to be evaluated in the computer is very complex, so it would be nice if this function could be substituted by a simpler and efficient one. To do so, the function f is sampled in a set of N pontos {x1, x2, . . . , xN}, where x(i) (is an element of) R(n), and then an estimate for the value of f in any other point is done by an interpolation method. An interpolation method is any procedure that takes a set of constraints and determines a nice function that satisfies such conditions. The Shepard interpolation method originally calculates the estimate of F(x) for some x (is an element of) R(n) as a weighted mean of the N sampled values of f. The weight for each sample xi is a function of the negative powers of the euclidian distances between the point x and xi. Kernels K : R(n) ×R(n) (IN) R are functions that correspond to an inner product on some Hilbert space F that contains the image of the points x and z by a function phi (the empty set) : R(n) (IN) F, i.e. k(x, z) =< phi (the empty set) (x), phi (the empty set) (z) >. In practice, the kernels represent implicitly the mapping phi (the empty set), i.e. it is more suitable to defines which kernel to use instead of which function phi (the empty set). This work proposes a simple modification on the Shepard interpolation method that is: to substitute the euclidian distance between the points x and xi by a distance between the image of these two point by phi (the empty set) in the Hilbert space F, which can be computed directly with the kernel k. Several tests show that such simple modification has better results when compared to the original method.

[pt] Cilindros fotossensíveis são usados nos processos de impressão e particularmente na impressão eletrofotográfica. O revestimento é aplicado ao cilindro em forma líquida antes de solidificar-se. O líquido é aplicado ao cilindro através de um aplicador de agulha que se translada ao longo da direção axial do cilindro. Durante este processo o cilindro gira em torno de seu próprio eixo levando a uma deposição da tira de líquido sobre a superfície do cilindro, em padrão espiral. Para ajudar a distribuir o líquido lateralmente e assim melhorar a uniformidade da espessura, cada tira de líquido aplicada pela agulha passa por meio de uma lâmina flexível. Este processo leva a um revestimento que apresenta um padrão espiral na espessura da camada revestida que pode causar defeitos no processo eletrofotográfico.O conhecimento, de forma fundamental, do escoamento é vital para a otimização do processo. Um modelo teórico de escoamento de filmes finos sobre superfícies cilíndricas em rotação com uma porta de injeção de líquido em movimento é apresentado. Este modelo é baseado na teoria de lubrificação considerando um filme precursor na frente da linha de contato aparente. A espessura de filme revestido foi obtida através da solução de uma equação diferencial parcial não linear de quarta ordem usando o método de diferenças finitas de segunda ordem. A discretização do tempo foi feita pelo método implícito de Crank- Nicholson. A discretização do sistema a cada passo do tempo leva a um sistema de equações algébricas não lineares que foi resolvido pelo método de Newton. Os resultados mostram como a uniformidade da camada depositada varia com os parâmetros do processo e com as propriedades do líquido. / [en] Photosensitive cylinders are used in printing arts and more particularly in electrophotographic printing (e.g. xerographic copy). The photosensitive coating is applied to the cylinder in liquid form, before it is solidified. The liquid is applied to the cylinder through a needle applicator that translates along the direction of the cylinder axis. The cylinder rotates during this process in order to cover the entire surface. Therefore, the liquid is applied in a spiral pattern. In order to help spreading of the liquid over the cylinder surface and to improve the thickness uniformity, each liquid stream applied from the needle passes under a flexible blade. This process leads to a coating that presents a spiral pattern on the deposited layer thickness, which can cause defects on the electrophotographic process. The complete understanding of the flow is vital to the optimization of the process. A theoretical model of the thin film flow over the surface of a rotating cylinder is presented here. It is based on the lubrication approximation considering a thin precursor film in front of the apparent contact line. The resulting non-linear fourth-order PDE for the film thickness was solved by a second-order finite difference method. The time discretization was done by the implicit Crank-Nicholson scheme. The non-linear algebraic equation at each time step was solved by Newton’s method. The results show how the uniformity of the deposited layer varies with process parameters and liquid properties.

[pt] Neste estudo faz-se uma análise do desempenho de antenas cônicas compactas para operar em largas bandas de freqüência. Esta análise será feita observando o comportamento da impedância da antena. Serão, primeiramente, consideradas as estruturas convencionais formadas por cones e cones sobrepostos próximos a um plano condutor. Posteriormente, será considerada uma estrutura cônica onde o arredondamento das bordas resulta em um comportamento mais uniforme para a impedância de entrada com variações da freqüência, conferindo um aumento na banda para esta configuração. / [en] The discone antenna is well known as an easy to build multioctave broadband antenna with a omnidirectional radiation pattern in the H-plane, It comprises a plane conductive disc element spaced close to and axially aligned with a conductive cone element. For applications in the UHF and microwave ranges, high performance antennas have been designed to operate from 0.5 to 5 GHz with a VSWR of about 3.5:1 or less. If the length is finite, the impedance still strongly dependent on apex angle as long as the cone is longer than about quarter of wavelength and the apex angle is relatively larger. A problem that is experienced with these antenna designs is the relatively large size required to operate at the low frequencies. In applications where utilization of this antenna. In this work, the analysis and design of these antenna is obtained by employing a rigorous formulation of the electromagnetic scattering problem. As a design tool, we employ Method of Moments for the analysis of rotationally symmetric structures excited by TEM mode. To properly account the variations in driven-point impedance with frequency, the coaxial waveguide used to feed the antenna is also included in the analysis. The excitation is simulated by a distribution of equivalent electric and magnetic currents placed inside the coaxial cable, over a cross-section plane. These currents are defined such that only excite the mode TEM towards the cable-antenna junction. These numerical tool is employed in the shaping the metallic surfaces involved in the discone-type antennas in order to obtain more compact structures. Simple solutions can be easily obtained making the center fed cone and the disc element with radially outer edge portions rolled backwardly and away from each other to form donut-shaped configurations. The ruled edges not only reduce the diameter of the antenna but also permit the current to flow around them smoothly and without reflections that degrade the VSWR at low end of the frequency band.

[pt] Estuda-se o espalhamento de ondas eletromagnéticas por condutores elétricos perfeitos com simetria de revolução. A formulação envolve a solução da Equação Integral para o Campo Elétrico (EFIE), aplicável a superfícies fechadas ou abertas, utilizado-se o Método dos Momentos com funções teste iguais às funções base (alternativa de Galerkin). São empregadas três diferentes classes de funções base: locais (triângulos e retângulos), globais (senos e cossenos) e limitadas em banda (funções sampling-like). Os resultados da análise (determinação da densidade de corrente elétrica induzida por uma onda plana axialmente incidente) de diversos tipos de objetos canônicos (esferas, discos e cone-esferas) ao utilizar-se as duas primeiras classes são comparados com aqueles encontráveis na literatura. Em seguida, analisa-se (determinação da corrente induzida, do campo espalhado e do campo irradiado com uma fonte esférica pontual situada no foco) refletores axialmente simétricos (parabolóides e hiperbolóides) com diâmetro entre 10 e 20 (lambda) empregando-se as três classes. É feita uma comparação entre os resultados variando-se o número de funções base por comprimento de onda para cada classe e os resultados obtidos ao empregar-se a Ótica Física (PO). / [en] The problem of electromagnetic scattering from perfectly conducting circularly symmetric bodies is considered. The mathematical formulation herein employed uses the Electrical Field Integral Equation (EFIE), which is applicable to either open or closed surfaces. The solution is effected by the Method of Moments. Test functions are chosen to be equal to basis functions (Galerkin´s alternative). Three different classes of basis functions are studied: local (pulses and triagles), global (sines and cossines) and band-limited ones (sampling-like functions). Results from the analysis (i. E., determination of induced electric currents)of different Kinds of cononical scatters (spheres, discs and cone-spheres) illuminated by an axially incident plane wave, using local and global functions, are compared with results found in the literature. In addition, circularly symmetric reflector antennas (paraboloids and hyperboloids with diameters between 10-20 lambda) are analyzed using the three different classes. The analysis involves the determination of induced currents, scattered and radiated fielsd. A comparision is made among the results obtained when employing different numbers of basis funcions and with the ones derived from Physical Optics (PO) approxomation.

[pt] Os métodos de Brown e Winters são, sem dúvida alguma, os métodos de amortecimento exponencial mais usados para a previsão de séries temporais. Entretanto, ambos podem ser considerados de aplicação limitada, pois ou não admitem componente sazonal (Brown) ou utilizam um modelo linear para a modelagem da tendência (Winters). Apresenta-se aqui uma generalização dos métodos de amortecimento na qual as limitações acima são eliminadas. Em particular, considera-se uma única formulação matemática para o modelo, composto dos termos de tendência (constante, linear ou quadrática) e componente sazonal sob a forma de um conjunto discreto de fatores (aditivos ou multiaplicativos). Fornece-se também uma estimativa para a variância dos erros de previsão, e é proposta uma forma de controle adaptativo para a constante de amortecimento da parte não sazonal. Foi feito um programa de computador que implementa automaticamente o método, inclusive estimando valores iniciais para o processo. Foram geradas e processadas algumas séries para exemplo e análise do desempenho do método. São fornecidas sugestões de pesquisa futura no sentido de possíveis aprimoramentos para o método, mas que demandam maior análise. / [en] The methods of Brown and Winters are, undoubtedly, the most popular exponential smoothing techniques used nowadays. However, both methods have limitations, such as: Brown s method is applicable only to non-seasonal series and Winters use the linear structure as the only possible model for the trend. A generalization of the smoothing methods in which the limitations cited above are eliminated is presented here. In particular, through a unique analytical formulation, the trend model (constant, linear or quadratic) is linked to the seasonal factors (additive or multiplicative). A forecast error variance estimator is provided and the adaptive control of the non-seasonal part smoothing constant is proposed. A computer program was written for automatic implementation of the method. This program also performs initial values estimation for the process initialization. Some series were generated and processed for testing the method performance. Several suggestions are given for future research which may yield, upon further analysis, to method improvement.

[pt] Esta tese apresenta as formulações dos métodos de elementos de contorno convencional e híbrido simplificado para problemas axissimétricos de elasticidade, empregando-se as soluções fundamentais do espaço completo e do semi-espaço. Para problemas de elasticidade axissimétricos no semi-espaço pelos métodos de elementos de contorno, o uso das soluções fundamentais para espaço completo exige a discretização e o truncamento da superfície livre. No entanto, essa discretização é dispensada se as soluções fundamentais empregadas satisfizerem a condição de forças de superfície nulas. Este trabalho apresenta a implementação dessas soluções fundamentais axissimétricas para o espaço completo e o semi-espaço elástico, em termos de integrais do tipo Lipschitz-Hankel. São apresentadas todas as expressões necessárias para o cálculo de resultados em pontos internos e a correta integração numérica das integrais de contorno. Partindo da formulação do espaço completo, mostra-se que é necessária pouca modificação para a implementação da formulação proposta. Esse trabalho também desenvolve a formulação axissimétrica para o método híbrido simplificado dos elementos de contorno, tanto para o espaço completo como para o semi-espaço. Na sua versão original, o uso de propriedades espectrais para a total formulação do problema não era possível para certas configurações topológicas. No entanto, a aplicação de um princípio de contragradiência híbrida às equações do método levou à obtenção de uma nova relação matricial que tornou possível sua total formulação para qualquer topologia, independentemente de propriedades espectrais. A necessidade de se integrar apenas uma matriz e a facilidade de obtenção de resultados em pontos internos tornam o método híbrido simplificado dos elementos de contorno ainda mais vantajoso para problemas axissimétricos. Alguns exemplos numéricos validam as formulações apresentadas. Essa tese é composta por oito capítulos e dois apêndices, como descritos a seguir. O Capítulo 2 trata das soluções fundamentais axissimétricas para o espaço completo e o semi-espaço elástico. As equações governantes para um meio elástico axissimétrico são apresentadas em coordenadas cilíndricas. As soluções fundamentais correspondentes são deduzidas, em termos de integrais do tipo Lipschitz-Hankel, a partir da solução de Muki das equações de equilíbrio de Navier. O Capítulo 3 apresenta o método dos elementos de contorno para problemas axissimétricos no espaço completo e no semi-espaço. A partir das soluções fundamentais apresentadas no Capítulo 2, as equações integrais no contorno são deduzidas, bem como as equações matriciais governantes. Além disso, discute-se a obtenção de uma matriz de rigidez e o cálculo das inversas generalizadas presentes na formulação. As expressões para o cálculo de deslocamentos e tensões no domínio e ao longo do contorno são fornecidas de maneira explícita. O Capítulo 4 apresenta o método híbrido simplificado dos elementos de contorno para problemas axissimétricos no espaço completo e no semi-espaço. Uma nova versão do método é proposta, em que as equações governantes do problema são obtidas a partir de trabalhos virtuais de deslocamentos, uma equação de compatibilidade de deslocamentos e um teorema híbrido de contragradiência. O esquema para o cálculo dos coeficientes indeterminados de U está descrito detalhadamente para o espaço completo, incluindo as soluções analíticas necessárias. A obtenção de uma matriz de rigidez, bem como de deslocamentos e tensões em pontos internos, também é discutida. Bases ortonormais, projetores e inversas generalizadas presentes na formulação são apresentados ao longo do capítulo. O Capítulo 5 apresenta os esquemas numéricos para o cálculo das integrais presentes nos métodos de elementos de contorno convencional e híbrido simplificado aplicados a problemas axissimétricos no espaço completo e no semi-espaço. / [en] This thesis presents the formulation of the conventional and simplified-hybrid boundary element methods for axisymmetric problems, employng fullspace as well as halfspace fundamental solutions. As it is mostly found in the literature on axisymmetric problems in the elastic halfspace, the boundary element formulations make use of fullspace fundamental solutions and insert a mesh discretization of the free surface, with truncation at a reasonable distance from the axis of axisymmetry. This discretization can be circunvented if one employs the fundamental solutions that satisfy in advance the traction free boundary condition on the free surface. This work presents the implementation of these axisymmetric fundamental solutions for both the fullspace and the halfspace, given in terms of integrals of Lipschitz-Hankel type. Explicit equations for post-processing results at internal points are provided, as well as the adequate numerical schemes to evaluate the boundary integrals that arise in the formulation. It is shown that the boundary element method for the halfspace can be easily implemented from existing computation codes for fullspace problems, requiring only a few modifications. This work also addresses the simplified-hybrid boundary element method for the axisymmetric fullspace and halfspace problems. In its original version, the use of spectral properties to completely formulate the method was possible for only strictly non-convex topological configurations. The key contribution of the present developments consisted in the correct application of a hybrid contragradient theorem to derive a simple means of using analytical solutions of the elastic problem in order to substitute for the spectral properties that have been originally proposed. In the simplified-hybrid boundary element method, only one matrix requires integration and the results at internal points can be evaluated directly, which makes the method extremely advantageous for axisymmetric problems. Some numerical examples validate the proposed formulations.

[pt] A presente tese apresenta uma investigação acerca do método da hipótese em Platão. Foram analisadas as passagens do Mênon (86e-87b), do Fédon (101d-e) e da República (VI, 509d-511e), lugares em que o método hipotético sustentaria um discurso científico, inspirado em métodos matemáticos. No estudo, foram utilizados fragmentos de autores antigos que tratam das matemáticas, além de textos de estudiosos e comentadores contemporâneos sobre a natureza da hipótese, a fim de subtraí-los de compreensões filosóficas anacrônicas em relação ao texto platônico. Entre as três passagens metodológicas foram estabelecidas semelhanças e dessemelhanças, para identificarmos até que ponto esse método se constitui ou não num método único. Caso apresente procedimentos diversos, como eles então, verdadeiramente se apresentam? O vocabulário geométrico utilizado por Platão; as disciplinas matemáticas e sua importância no programa de educação platônico; o método do elenchos e o método analítico dos geômetras gregos, também são objetos de estudo da tese. / [en] This thesis is an investigation of Plato s method of hypothesis. With this purpose, the following passages of Plato s dialogues were analyzed, in which the hypothetical method is supposed to sustain a scientific discourse, inspired by mathematical methods: Meno (86e-87b), Phaedo (101d-e) and Republic (VI, 509d- 511e). In this study, we have used fragments of ancient authors on mathematics, besides texts by contemporary scholars and commentators on the nature of the notion of hypothesis as used by Plato, in order to exclude any anachronic philosophical understanding of it. Among the three methodological passages, similarities and dissimilarities were established, in order to verify to which degree this method is or is not one and the same. In case it shows different procedures, what are they, really? The geometrical vocabulary used by Plato, the mathematical disciplines and their importance in the platonic educational program, the socratic elenchus and the analytical method of the Greek geometers are also objects of investigation.

[pt] Este trabalho faz parte de um projeto para a implementação de um programa que possa simular problemas com milhões de graus de liberdade em um computador pessoal. Para isto, combina-se o Método Fast Multipole (FMM) com o Método Expedito dos Elementos de Contorno (EBEM), além de serem utilizados resolvedores iterativos de sistemas de equações. O EBEM é especialmente vantajoso em problemas de complicada topologia, ou que usem funções fundamentais muito complexas. Neste trabalho apresenta-se uma formulação para o Método Fast Multipole (FMM) que pode ser usada para, virtualmente, qualquer função e também para contornos curvos, o que parece ser uma contribuição original. Esta formulação apresenta um formato mais compacto do que as já existentes na literatura, e também pode ser diretamente aplicada a diversos tipos de problemas praticamente sem modificação de sua estrutura básica. É apresentada a validação numérica da formulação proposta. Sua utilização em um contexto do EBEM permite que um programa prescinda de integrações sobre segmentos – mesmo curvos – do contorno quando estes estão distantes do ponto fonte. / [en] This is part of a larger project that aims to develop a program able to simulate problems with millions of degrees of freedom on a personal computer. The Fast Multipole Method (FMM) is combined with the Expedite Boundary Element Method (EBEM) for integration, in the project s final version, with iterative equations solvers. The EBEM is especially advantageous when applied to problems with complicated topology as well as in the case of highly complex fundamental solutions. In this work, a FMM formulation is proposed for the use with virtually any type of fundamental solution and considering curved boundaries, which seems to be an original contribution. This formulation presents a more compact format than the ones shown in the technical literature, and can be directly applied to different kinds of problems without the need of manipulation of its basic structure, being numerically validated for a few applications. Its application in the context of the EBEM leads to the straightforward implementation of higher-order elements for generally curved boundaries that dispenses integration when the boundary segment is relatively far from the source point.

[pt] Este trabalho tem como finalidade a apresentação de um método novo e eficiente para a análise e separação de tensões em materiais compostos (conjugados). O método híbrido consiste no acoplamento da lei ótica para materiais ortotrópicos com as equações de compatibilidade e equilíbrio para corpos anisotrópicos. As soluções destas equações são representadas por funções analíticas de variáveis complexas. O método propõe que se aproximem as funções analíticas por funções polinomiais. Estas funções aproximadas são determinadas através do acoplamento das equações de equilíbrio e compatibilidade com dados experimentais, obtidos da resposta fotoelástica, e com métodos numéricos (Newton-Raphson e mínimos quadrados), possibilitando a determinação dos coeficientes integrantes das funções polinomiais. O método apresenta como principal vantagem a capacidade de separar as tensões, com precisão, em pequenas ou mesmo em grandes regiões de modelos fotoelásticos planos, sem requerer a leitura de franjas isoclínicas. São apresentados excelentes resultados obtidos com a aplicação do método desenvolvido em alguns casos específicos, envolvendo um material composto por uma matriz de resina reforçada por fibras de vidro com arranjo unidirecional. Três casos são objetos de estudo: A) Barra retangular submetida à tração, com a direção do reforço paralela a direção da carga aplicada. B) Barra retangular submetida à tração. com direção do reforço fazendo um ângulo qualquer com a direção da carga aplicada. C) Plano semi-infinito submetido à carga concentrada atuando perpendicularmente à direção do reforço. / [en] The objective of this paper is to present a new and efficient method to analyse and separate stresses in composite materials. The hybrid method associates the stress-optic law for orthotropic materials with the compatibility and equlibrium equations for anisotropic bodies. The solutions of these equatios are represented by analytical functions of complex variables. The method proposes that these analytical functions be approximate by polynomial functions. These aapproximate functions are determined by the association of the compatibility equilibrium equatios with experimental data (obtained from the photoelastic response) and numerical methods (Newton-Raphson and least squares method). The method presents as a major advantage the capacity to separate stresses precisely in small or even in large regions of photoelastic plane models, without needign to collect isoclinic fringe data. The thesis presents excellent results obtained by the application of the developed method in some specific cases, where a fiber glass reinforced epoxy resin is used as a model material, will be shown. Three cases are studied: A) Rectangular bar in tension with the fiber reinforcement axis parallel to the applied load direction B) Rectangular bar in tension with the fiber reinforcement axis making a known angle with the applied load direction C) Semi-infinite plane subjected to a concentrated edge load which is perpendicular to the fiber reinforcement axis.