Computational modeling for irrigated agriculture planning. Part I: general description and linear programming

Linear programming models are effective tools to support initial or periodic planning of agricultural enterprises, requiring, however, technical coefficients that can be determined using computer simulation models. This paper, presented in two parts, deals with the development, application and tests of a methodology and of a computational modeling tool to support planning of irrigated agriculture activities. Part I aimed at the development and application, including sensitivity analysis, of a multiyear linear programming model to optimize the financial return and water use, at farm level for Jaíba irrigation scheme, Minas Gerais State, Brazil, using data on crop irrigation requirement and yield, obtained from previous simulation with MCID model. The linear programming model outputted a crop pattern to which a maximum total net present value of R$ 372,723.00 for the four years period, was obtained. Constraints on monthly water availability, labor, land and production were critical in the optimal solution. In relation to the water use optimization, it was verified that an expressive reductions on the irrigation requirements may be achieved by small reductions on the maximum total net present value.

Comparison of maps of spatial variability of soil resistance to penetration constructed with and without covariables using a spatial linear model

A study about the spatial variability of data of soil resistance to penetration (RSP) was conducted at layers 0.0-0.1 m, 0.1-0.2 m and 0.2-0.3 m depth, using the statistical methods in univariate forms, i.e., using traditional geostatistics, forming thematic maps by ordinary kriging for each layer of the study. It was analyzed the RSP in layer 0.2-0.3 m depth through a spatial linear model (SLM), which considered the layers 0.0-0.1 m and 0.1-0.2 m in depth as covariable, obtaining an estimation model and a thematic map by universal kriging. The thematic maps of the RSP at layer 0.2-0.3 m depth, constructed by both methods, were compared using measures of accuracy obtained from the construction of the matrix of errors and confusion matrix. There are similarities between the thematic maps. All maps showed that the RSP is higher in the north region.

Impact of elevated temperature scenarios on potato leaf development

The objective of this study was to simulate the impact of elevated temperature scenarios on leaf development of potato in Santa Maria, RS, Brazil. Leaf appearance was estimated using a multiplicative model that has a non-linear temperature response function which calculates the daily leaf appearance rate (LAR, leaves day-1) and the accumulated number of leaves (LN) from crop emergence to the appearance of the upper last leaf. Leaf appearance was estimated during 100 years in the following scenarios: current climate, +1 °C, +2 °C, +3 °C, +4 °C e +5 °C. The LAR model was estimated with coefficients of the Asterix cultivar in five emergence dates and in two growing seasons (Fall and Spring). Variable of interest was the duration (days) of the crop emergence to the appearance of the final leaf number (EM-FLN) phase. Statistical analysis was performed assuming a three-factorial experiment, with main effects being climate scenarios, growing seasons, and emergence dates in a completely randomized design using years (one hundred) as replications. The results showed that warmer scenarios lead to an increase, in the fall, and a decrease, in the spring growing season, in the duration of the leaf appearance phase, indicating high vulnerability and complexity of the response of potato crop grown in a Subtropical environment to climate change.

Correlação linear e espacial entre produtividade de Brachiaria brizantha, densidade do solo e porosidade total em função do sistema de manejo do solo

A densidade e a porosidade são propriedades físicas do solo que são alteradas em função do sistema de manejo utilizado, com consequente influência sobre a produtividade das culturas. A produtividade da planta forrageira, Brachiaria brizantha, em função da densidade do solo e da porosidade total foi analisada em dois sistemas de manejo de solo, em experimento conduzido no segundo semestre de 2009, no município de Santa Teresa, no Estado do Espírito Santo. O objetivo foi estudar a variabilidade e as correlações lineares e espaciais entre os atributos da planta e do solo, visando a selecionar um indicador da qualidade física do solo de boa representatividade para produtividade de forragem. Marcaram-se duas parcelas de 40 m por 50 m a cada 5 m, em duas direções, resultando em um reticulado retangular de 99 pontos, em cada um dos sistemas utilizados: preparo convencional e plantio direto. Os atributos estudados, além de não terem variado aleatoriamente, apresentaram variabilidade dos dados entre média e alta, e seguiram padrões espaciais bem definidos, com alcance entre 20,3 e 24,2 m. Por sua vez, a correlação linear entre o atributo da planta e os do solo, em função do elevado número de observações, foi baixa. A melhor correlação para produtividade de matéria seca foi com a densidade do solo na profundidade de 0,0 - 0,15 m, independentemente do sistema de manejo do solo, indicando que a produtividade e a densidade do solo são inversamente proporcionais. Portanto, a densidade do solo avaliada na camada de 0,0 - 0,15 m apresentou-se como satisfatório indicador da qualidade física do solo, quando se considerou a produtividade da forrageira.

Influência local em um modelo espacial linear da produtividade da soja utilizando distribuição t-Student

O uso das ferramentas da geoestatística, aliadas à agricultura de precisão permitem o acompanhamento das áreas agrícolas produtoras de soja, estabelecendo as relações de dependência espacial entre os pontos amostrados. A modelagem da estrutura de variabilidade espacial possibilita a construção de mapas temáticos dos atributos estudados, utilizando como método de interpolação a krigagem. Porém, a presença de valores atípicos entre os elementos amostrais pode influenciar na construção e interpretação desses mapas. A distribuição de probabilidades t-Student tem sido utilizada na tentativa de diminuir a influência dos valores atípicos durante a estimativa dos parâmetros de dependência espacial, por ter caudas mais pesadas que a distribuição normal. A detecção dos valores influentes na área em estudo, por meio da análise de diagnósticos de influência local, confere maior confiabilidade na utilização dos mapas gerados, corroborando a aplicação de insumos. Deste modo, o objetivo deste trabalho foi aplicar as técnicas de influência local em dados espacialmente referenciados, com os modelos de perturbação aditiva e utilizando a matriz escala, considerando a distribuição t-Student n-variada. Foi utilizado um modelo espacial linear para o estudo de dados da produtividade da soja em função da altura média de plantas e do número médio de vagens por planta. As técnicas de influência local foram eficientes para detectar pontos que influenciam na escolha do modelo geoestatístico, nas estimativas dos parâmetros e na construção do mapa temático.

Anastomose esôfago-jejunal látero-lateral com grampeador linear cortante complementada por sutura manual

One of the most difficult procedures in digestive-tract surgery is esophago-jejunal anastomosis following total gastrectomy. Cost/benefit analysis of this procedure justifies the use of mechanical staplers, in spite of their high cost. A technical variant of the side-to-side esophago-jejunal anastomosis is presented, which incorporates the use of a cutting linear stapler. Technical maneuvers are easy to perform, the cost of the cutting linear stapler is smaller than the circular ones, the amplitude of the anastomosis is wider and the likelihood of fistulae is smaller when compared to other techniques. The side-to-side esophago-jejunal anastomosis with the cutting linear stapler is always complemented by a manual suture.

Predictors of fluid intravasation during operative hysteroscopy: a preplanned prospective observational study with 200 cases

PURPOSE: To verify the predictors of intravasation rate during hysteroscopy.METHODS: Prospective observational study (Canadian Task Force classification II-1). All cases (n=200 women; 22 to 86 years old) were treated in an operating room setting. Considering respective bag overfill to calculate water balance, we tested two multiple linear regression models: one for total intravasation (mL) and the other for absorption rate (mL.min-1). The predictors tested (independent variables) were energy (mono/bipolar), tube patency (with/without tubal ligation), hysterometry (cm), age≤50 years, body surface area (m2), surgical complexity (with/without myomectomy) and duration (min).RESULTS: Mean intravasation was significantly higher when myomectomy was performed (442±616 versus 223±332 mL; p<0.01). In the proposed multiple linear regression models for total intravasation (adjusted R2=0.44; p<0.01), the only significant predictors were myomectomy and duration (p<0.01).In the proposed model for intravasation rate (R2=0.39; p<0.01), only myomectomy and hysterometry were significant predictors (p=0.02 and p<0.01, respectively).CONCLUSIONS: Not only myomectomy but also hysterometry were significant predictors of intravasation rate during operative hysteroscopy.

Linear-quadratic optimal control of descriptor systems

In recent years the analysis and synthesis of (mechanical) control systems in descriptor form has been established. This general description of dynamical systems is important for many applications in mechanics and mechatronics, in electrical and electronic engineering, and in chemical engineering as well. This contribution deals with linear mechanical descriptor systems and its control design with respect to a quadratic performance criterion. Here, the notion of properness plays an important role whether the standard Riccati approach can be applied as usual or not. Properness and non-properness distinguish between the cases if the descriptor system is exclusively governed by the control input or by its higher-order time-derivatives additionally. In the unusual case of non-proper systems a quite different problem of optimal control design has to be considered. Both cases will be solved completely.

BEM/FEM non-linear model applied to transient analysis with viscous damping

In this article a two-dimensional transient boundary element formulation based on the mass matrix approach is discussed. The implicit formulation of the method to deal with elastoplastic analysis is considered, as well as the way to deal with viscous damping effects. The time integration processes are based on the Newmark rhoand Houbolt methods, while the domain integrals for mass, elastoplastic and damping effects are carried out by the well known cell approximation technique. The boundary element algebraic relations are also coupled with finite element frame relations to solve stiffened domains. Some examples to illustrate the accuracy and efficiency of the proposed formulation are also presented.

Non-linear analysis of laminated metal matrix composites by an integrated micro/macro-mechanical model

The present paper describes an integrated micro/macro mechanical study of the elastic-viscoplastic behavior of unidirectional metal matrix composites (MMC). The micromechanical analysis of the elastic moduli is based on the Composites Cylinder Assemblage model (CCA) with comparisons also draw with a Representative Unit Cell (RUC) technique. These "homogenization" techniques are later incorporated into the Vanishing Fiber Diameter (VFD) model and a new formulation is proposed. The concept of a smeared element procedure is employed in conjunction with two different versions of the Bodner and Partom elastic-viscoplastic constitutive model for the associated macroscopic analysis. The formulations developed are also compared against experimental and analytical results available in the literature.

Heat Source Estimation with the Conjugate Gradient Method in Inverse Linear Diffusive Problems

In this work, we present the solution of a class of linear inverse heat conduction problems for the estimation of unknown heat source terms, with no prior information of the functional forms of timewise and spatial dependence of the source strength, using the conjugate gradient method with an adjoint problem. After describing the mathematical formulation of a general direct problem and the procedure for the solution of the inverse problem, we show applications to three transient heat transfer problems: a one-dimensional cylindrical problem; a two-dimensional cylindrical problem; and a one-dimensional problem with two plates.

Non-Linear Unsteady Aerodynamic Response Approximation Using Multi-Layer Functionals

Non-linear functional representation of the aerodynamic response provides a convenient mathematical model for motion-induced unsteady transonic aerodynamic loads response, that accounts for both complex non-linearities and time-history effects. A recent development, based on functional approximation theory, has established a novel functional form; namely, the multi-layer functional. For a large class of non-linear dynamic systems, such multi-layer functional representations can be realised via finite impulse response (FIR) neural networks. Identification of an appropriate FIR neural network model is facilitated by means of a supervised training process in which a limited sample of system input-output data sets is presented to the temporal neural network. The present work describes a procedure for the systematic identification of parameterised neural network models of motion-induced unsteady transonic aerodynamic loads response. The training process is based on a conventional genetic algorithm to optimise the network architecture, combined with a simplified random search algorithm to update weight and bias values. Application of the scheme to representative transonic aerodynamic loads response data for a bidimensional airfoil executing finite-amplitude motion in transonic flow is used to demonstrate the feasibility of the approach. The approach is shown to furnish a satisfactory generalisation property to different motion histories over a range of Mach numbers in the transonic regime.

On non-ideal and non-linear portal frame dynamics analysis using bogoliubov averaging method

We apply the Bogoliubov Averaging Method to the study of the vibrations of an elastic foundation, forced by a Non-ideal energy source. The considered model consists of a portal plane frame with quadratic nonlinearities, with internal resonance 1:2, supporting a direct current motor with limited power. The non-ideal excitation is in primary resonance in the order of one-half with the second mode frequency. The results of the averaging method, plotted in time evolution curve and phase diagrams are compared to those obtained by numerically integrating of the original differential equations. The presence of the saturation phenomenon is verified by analytical procedures.

Understanding the complexity in low dimensional systems

Complex System is any system that presents involved behavior, and is hard to be modeled by using the reductionist approach of successive subdivision, searching for ''elementary'' constituents. Nature provides us with plenty of examples of these systems, in fields as diverse as biology, chemistry, geology, physics, and fluid mechanics, and engineering. What happens, in general, is that for these systems we have a situation where a large number of both attracting and unstable chaotic sets coexist. As a result, we can have a rich and varied dynamical behavior, where many competing behaviors coexist. In this work, we present and discuss simple mechanical systems that are nice paradigms of Complex System, when they are subjected to random external noise. We argue that systems with few degrees of freedom can present the same complex behavior under quite general conditions.