exact energy and momentum conserving algorithms for general models in nonlinear elasticity implicit time integration schemes that inherit the conservation laws of total energy linear and angular momentum are considered for initial boundary value problems in finite deformation elastodynamics conserving schemes are constructed for general hyperelastic material models both compressible and incompressible and are formulated in a way that is independent of spatial discretization three numerical examples for ogden type material models implemented using a finite element discretization in space are given to illustrate the performance of the proposed schemes these examples show that relative to the standard implicit mid point rule the conserving schemes exhibit superior numerical stability properties without a compromise in accuracy

a new dissipative time stepping algorithm for frictional contact problems formulation and analysis this paper presents a new time stepping algorithm for frictional contact problems that exhibits unconditional positive energy dissipation more specifically the proposed scheme preserves a priori stability estimates of the continuum problem for both frictionless and frictional contact leading to improved numerical stability properties in particular for the normal contact component the algorithm exhibits full energy conservation between released states while the energy does not increase over its initial value due to the enforcement of the normal contact constraint during persistent contact a penalty regularization is considered to this purpose a new regularization of the stick conditions is considered for the frictional part the new scheme is shown rigorously to exhibit positive energy dissipation like the continuum physical problem in this frictional case coulomb friction is assumed complete analyses of these considerations as well as a detailed description of their finite element implementation are included in the general finite deformation range representative numerical simulations are presented to assess the performance of the newly proposed methods

variational and projection methods for the volume constraint in finite deformation elasto plasticity this paper focuses on the treatment of volume constraints which in the context of elasto plasticity typically arise as a result of assuming volume preserving plastic flow projection methods based on the modification of the discrete gradient operator b often proposed on an ad hoc basis are systematically obtained in the variational context furnished by a three field hu washizu principle the fully nonlinear formulation proposed here is based on a local multiplicative split of the deformation gradient into volume preserving and dilatational parts without relying on rate forms of the weak form of momentum balance this approach fits naturally in a formulation of plasticity based on the multiplicative decomposition of the deformation gradient and enables one to exactly enforce the condition of volume preserving plastic flow within the framework proposed in this paper rate forms and incrementally objective algorithms are entirely bypassed

universal approximation capability of cascade correlation for structures cascade correlation cc constitutes a training method for neural networks that determines the weights as well as the neural architecture during training various extensions of cc to structured data have been proposed recurrent cascade correlation rcc for sequences recursive cascade correlation reccc for tree structures with limited fan out and contextual recursive cascade correlation creccc for rooted directed positional acyclic graphs dpags with limited fan in and fan out we show that these models possess the universal approximation property in the following sense given a probability measure p on the input set every measurable function from sequences into a real vector space can be approximated by a sigmoidal rcc up to any desired degree of accuracy up to inputs of arbitrary small probability every measurable function from tree structures with limited fan out into a real vector space can be approximated by a sigmoidal reccc with multiplicative neurons up to any desired degree of accuracy up to inputs of arbitrary small probability for sigmoidal creccc networks with multiplicative neurons we show the universal approximation capability for functions on an important subset of all dpags with limited fan in and fan out for which a specific linear representation yields unique codes we give one sufficient structural condition for the latter property which can easily be tested the enumeration of ingoing and outgoing edges should be compatible this property can be fulfilled for every dpag with fan in and fan out two via reenumeration of children and parents and for larger fan in and fan out via an expansion of the fan in and fan out and reenumeration of children and parents in addition the result can be generalized to the case of input output isomorphic transductions of structures thus creccc networks constitute the first neural models for which the universal approximation capability of functions involving fairly general acyclic graph structures is proved

similarity invariant signatures for partially occluded planar shapes a methodology is described for associating local invariant signature functions to smooth planar curves in order to enable their translation rotation and scale invariant recognition from arbitrarily clipped portions the suggested framework incorporates previous approaches based on locating inflections curvature extrema breakpoints and other singular points on planar object boundaries and provides a systematic way of deriving novel invariant signature functions based on curvature or cumulative turn angle of curves these new signatures allow the specification of arbitrarily dense feature points on smooth curves whose locations are invariant under similarity transformations the results are useful for detecting and recognizing partially occluded planar objects a key task in low level robot vision

energy dissipative momentum conserving time stepping algorithms for finite strain multiplicative plasticity abstract this paper presents a new energy dissipative momentum conserving algorithm for multiplicative models of finite strain plasticity the proposed scheme preserves exactly the conservation laws of linear and angular momenta and leads to the exact energy dissipation between the two computed solutions in a given time step in particular the energy is exactly conserved if the step is elastic thus recovering existing energy momentum conserving schemes in the purely elastic range extensions accommodating a strictly non negative numerical energy dissipation in the high frequency to handle the numerical stiffness of the problem are also discussed briefly these conservation dissipation properties are proven rigorously in the very nonlinear setting of multiplicative finite strain plasticity in fact the analysis drives the design of the new return mapping algorithm proposed to integrate the plastic evolution equations so the above dissipation conservation properties hold these analyses account for the spatial discretization of the problem in the context of the finite element method several representative numerical simulations are presented to illustrate the performance of the new algorithm

on automatic threshold selection for polygonal approximations of digital curves polygonal approximation is a very common representation of digital curves a polygonal approximation depends on a parameter e which is the error value in this paper we present a method for an automatic selection of the error value e let e be a polygonal approximation of the original curve with an error value e we define a set of function ns e s s s such that for a given value of s ns e is the number of edges that contain at least s vertices in e the time complexity for computing the set of functions ns e s ss s is almost linear in n the number of vertices in in this paper we analyse the ns graph and show that for adequate values of s a wide plateau is expected to appear at the top of the graph this plateau corresponds to a stable state in the multi scale representation of e we show that the functions ns s s are a statistical representation of some kind of scale space image

formulation and analysis of conserving algorithms for frictionless dynamic contact impact problems this paper presents the formulation of conserving time stepping algorithms for frictionless dynamic contact of solids a new class of finite element methods is proposed for the solution of these problems that exhibit the same conservation laws as the underlying continuum dynamical system the proposed methods are based on a penalty regularization of the constrained contact problem and lead to full conservation of the total energy of the system including the regularization penalty potential during persistent contact and restoration of the original energy upon release both linear and angular momenta are conserved by the scheme furthermore the newly developed methods have the ability to enforce the associated constraints in the velocity besides the impenetrability constraint in the displacements while preserving the conservation restoration properties of the final numerical scheme a modification of these schemes is described that assures positive energy dissipation if desired even in the highly nonlinear setting of contact impact problems leading to contact schemes with high frequency energy dissipation representative numerical simulations are presented illustrating the performance of the proposed numerical schemes

a first order hyperbolic framework for large strain computational solid dynamics an upwind cell centred total lagrangian scheme this paper builds on recent work developed by the authors for the numerical analysis of large strain solid dynamics by introducing an upwind cell centred hexahedral finite volume framework implemented within the open source code openfoam http www openfoam com http www openfoam com in lee gil and bonet 1 a first order hyperbolic system of conservation laws was introduced in terms of the linear momentum and the deformation gradient tensor of the system leading to excellent behaviour in two dimensional bending dominated nearly incompressible scenarios the main aim of this paper is the extension of this algorithm into three dimensions its tailor made implementation into openfoam and the enhancement of the formulation with three key novelties first the introduction of two different strategies in order to ensure the satisfaction of the underlying involutions of the system that is that the deformation gradient tensor must be curl free throughout the deformation process second the use of a discrete angular momentum projection algorithm and a monolithic total variation diminishing runge kutta time integrator combined in order to guarantee the conservation of angular momentum third and for comparison purposes an adapted total lagrangian version of the hyperelastic glace nodal scheme of kluth and despr es 2 is presented a series of challenging numerical examples are examined in order to assess the robustness and accuracy of the proposed algorithm benchmarking it against an ample spectrum of alternative numerical strategies developed by the authors in recent publications

shape representation by multiscale contour approximation an approach is presented for deriving qualitative descriptions of contours containing structures at different unknown scales the descriptions are in terms of straight arcs curved arcs with sign of curvature corners and points delimiting the arcs inflexion points and transitions from straight to curved furthermore the tangents at these points are derived the approach is based on the construction of a hierarchic family of polygons having the scale space property of causality structure can only disappear as scale goes from fine to coarse using the principle that structures that are stable over scale represent significant properties the features of the descriptive representations are then derived

fast decomposition of digital curves into polygons using the haar transform the perimeters of features in digital images are decomposed into sets of straight lines using the haar transform the haar basis functions preserve sharp changes in direction along the perimeter while still allowing smoothing the coefficients of the haar transform are simply related to the error between the polygonal decomposition and the true perimeter the haar transform can be calculated efficiently with a computational overhead linear in the number of data points

design of energy conserving algorithms for frictionless dynamic contact problems this paper proposes a formulation of dynamic contact problems which enables exact algorithmic conservation of linear momentum angular momentum and energy in finite element simulations it is seen that a lagrange multiplier enforcement of an appropriate contact rate constraint produces these conservation properties a related method is presented in which a penalty regularization of the aforementioned rate constraint is utilized this penalty method sacrifices the energy conservation property but is dissipative under all conditions of changing contact so that the global algorithm remains stable notably it is also shown that augmented lagrangian iteration utilizing this penalty kernel reproduces the energy conserving i e lagrange multiplier solution to any desired degree of accuracy the result is a robust stable method even in the context of large deformations as is shown by some representative numerical examples in particular the ability of the formulation to produce accurate results where more traditional integration schemes fail is emphasized by the numerical simulations 1997 by john wiley sons ltd

an algebraic framework to represent finite state machines in single layer recurrent neural networks in this paper we present an algebraic framework to represent finite state machines fsms in single layer recurrent neural networks slrnns which unifies and generalizes some of the previous proposals this framework is based on the formulation of both the state transition function and the output function of an fsm as a linear system of equations and it permits an analytical explanation of the representational capabilities of first order and higher order slrnns the framework can be used to insert symbolic knowledge in rnns prior to learning from examples and to keep this knowledge while training the network this approach is valid for a wide range of activation functions whenever some stability conditions are met the framework has already been used in practice in a hybrid method for grammatical inference reported elsewhere sanfeliu and alquezar 1994

energy consistent algorithms for dynamic finite deformation plasticity abstract this paper addresses the theoretical development and numerical implementation of energy consistent algorithms for dynamic elastoplasticity emphasizing finite strain constitutive formulations so that unconditional stability of the algorithms is assured even in the fully nonlinear regime the key concept behind energy consistency is the requirement that the discretized system obey an a priori stability estimate which may be derived in general using the first and second laws of thermodynamics this approach to computational dynamic plasticity differs from typical application of traditional algorithms such as newmark or hilber hughes taylor methods where local time integration schemes for plasticity laws are developed somewhat independently from the global time integration scheme for the equations of motion without explicit consideration of thermodynamical restrictions two algorithms based on both additive and multiplicative finite deformation plasticity model are formulated within the energy consistent framework both algorithms possess the desirable feature of nonlinear stability of previous energy momentum algorithms for elastodynamics

a hyperelastic based large strain elasto plastic constitutive formulation with combined isotropic kinematic hardening using the logarithmic stress and strain measures summary this paper addresses the formulation of a set of constitutive equations for finite deformation metal plasticity the combined isotropic kinematic hardening model of the infinitesimal theory of plasticity is extended to the large strain range on the basis of three main assumptions i the formulation is hyperelastic based ii the stress strain law preserves the elastic constants of the infinitesimal theory but is written in terms of the hencky strain tensor and its elastic work conjugate stress tensor and iii the multiplicative decomposition of the deformation gradient is adopted since no stress rates are present the formulation is of course numerically objective in the time integration it is shown that the model gives adequate physical behaviour and comparison is made with an equivalent constitutive model based on the additive decomposition of the strain tensor

a new energy and momentum conserving algorithm for the non linear dynamics of shells a numerical time integration scheme for the dynamics of non linear elastic shells is presented that simultaneously and independent of the time step size inherits exactly the conservation laws of total linear total angular momentum as well as total energy the proposed technique generalizes to non linear shells recent work of the authors on non linear elastodynamics and is ideally suited for long term large scale simulations the algorithm is second order accurate and can be immediately extended with no modification to a fourth order accurate scheme the property of exact energy conservation induces a strong notion of non linear numerical stability which manifests itself in actual simulations

scale based detection of corners of planar curves a technique for detecting and localizing corners of planar curves is proposed the technique is based on gaussian scale space which consists of the maxima of absolute curvature of the boundary function presented at all scales the scale space of isolated simple and double corners is first analyzed to investigate the behavior of scale space due to smoothing and interactions between two adjacent corners the analysis shows that the resulting scale space contains line patterns that either persist terminate or merge with a neighboring line next the scale space is transformed into a tree that provides simple but concise representation of corners at multiple scales finally a multiple scale corner detection scheme is developed using a coarse to fine tree parsing technique the parsing scheme is based on a stability criterion that states that the presence of a corner must concur with a curvature maximum observable at a majority of scales experiments were performed to show that the scale space corner detector is reliable for objects with multiple size features and noisy boundaries and compares favorably with other corner detectors tested

non linear dynamics of three dimensional rods exact energy and momentum conserving algorithms the long term dynamic response of non linear geometrically exact rods under going finite extension shear and bending accompanied by large overall motions is addressed in detail the central objective is the design of unconditionally stable time stepping algorithms which exactly preserve fundamental constants of the motion such as the total linear momentum the total angular momentum and for the hamiltonian case the total energy this objective is accomplished in two steps first a class of algorithms is introduced which conserves linear and angular momentum this result holds independently of the definition of the algorithmic stress resultants second an algorithmic counterpart of the elastic constitutive equations is developed such that the law of conservation of total energy is exactly preserved conventional schemes exhibiting no numerical dissipation symplectic algorithms in particular are shown to lead to unstable solutions when the high frequencies are not resolved compared to conventional schemes there is little if any additional computational cost involved in the proposed class of energy momentum methods the excellent performance of the new algorithm in comparison to other standard schemes is demonstrated in several numerical simulations

a new solution procedure for application of energy conserving algorithms to general constitutive models in nonlinear elastodynamics abstract this paper develops a novel solution strategy for temporal discretization of initial value problems in nonlinear elastodynamics enabling exact energy and momentum conservation for a spatially discretized problem this new energy momentum algorithm fits into the general framework established by simo and tarnow z angew math phys 43 1992 757 however in contrast to the presentation given in that reference the current approach extends the critical stress update formula to encompass generic stored energy functions for the hyperelastic continuum in a manner which also allows for quadratically convergent newton raphson solution of the nonlinear algorithmic equations numerical examples using a neo hookean material model are given to illustrate the excellent performance of the proposed scheme not only is the new algorithm guaranteed to be second order accurate and unconditionally stable in the fully nonlinear regime as are other energy momentum strategies but also the approach advocated is demonstrated to be more robust and versatile than its predecessors

improved implicit integrators for transient impact problems geometric admissibility within the conserving framework the value of energy and momentum conserving algorithms has been well established for the analysis of highly non linear systems including those characterized by the nonsmooth non linearities of an impact event this work proposes an improved integration scheme for frictionless dynamic contact seeking to preserve the stability properties of exact energy and momentum conservation without the heretofore unavoidable compromise of violating geometric admissibility as established by the contact constraints the physically motivated introduction of a discrete contact velocity provides an algorithmic framework that ensures exact conservation locally while remaining independent of the choice of constraint treatment thus making full conservation equally possible in conjunction with a penalty regularization as with an exact lagrange multiplier enforcement the discrete velocity effects are incorporated as a post convergence update to the system velocities and thus have no direct effect on the non linear solution of the displacement equilibrium equation the result is a robust implicit algorithmic treatment of dynamic frictionless impact appropriate for large deformations and fully conservative for a range of geometric constraints copyright 2001 john wiley sons ltd

algorithms for coupled problems that preserve symmetries and the laws of thermodynamics part i monolithic integrators and their application to finite strain thermoelasticity abstract we propose discretizations for a general class of nonlinear coupled thermomechanical problems of evolution the most salient feature of the new methods is that they rigorously preserve the two laws of thermodynamics as well as the symmetries of the systems they model to formulate such methods we exploit the geometric structure afforded by the generic formalism of non equilibrium thermodynamics and we follow a systematic methodology that results in discrete evolution equations which mimic the generic structure as an illustration a complete discretization of finite strain thermoelasticity is presented using finite elements in space and a monolithic integrator in time simulations are provided which demonstrate the conservation features of the algorithm as well as its remarkable robustness

formulation of implicit finite element methods for multiplicative finite deformation plasticity some constitutive and computational aspects of finite deformation plasticity are discussed attention is restricted to multiplicative theories of plasticity in which the deformation gradients are assumed to be decomposable into elastic and plastic terms it is shown by way of consistent linearization of momentum balance that geometric terms arise which are associated with the motion of the intermediate configuration and which in general render the tangent operator non symmetric even for associated plastic flow both explicit i e no equilibrium iteration and implicit finite element formulations are considered an assumed strain formulation is used to accommodate the near incompressibility associated with fully developed isochoric plastic flow as an example of explicit integration the rate tangent modulus method is reviewed in some detail an implicit scheme is derived for which the consistent tangents resulting in quadratic convergence of the equilibrium iterations can be written out in closed form for arbitrary material models all the geometrical terms associated with the motion of the intermediate configuration and the treatment of incompressibility are given explicitly examples of application to void growth and coalescence and to crack tip blunting are developed which illustrate the performance of the implicit method

provably unconditionally stable second order time accurate mixed variational methods for phase field models we introduce provably unconditionally stable mixed variational methods for phase field models our formulation is based on a mixed finite element method for space discretization and a new second order accurate time integration algorithm the fully discrete formulation inherits the main characteristics of conserved phase dynamics namely mass conservation and nonlinear stability with respect to the free energy we illustrate the theory with the cahn hilliard equation but our method may be applied to other phase field models we also propose an adaptive time stepping version of the new time integration method we present some numerical examples that show the accuracy stability and robustness of the new method

matching of partially occluded planar curves abstract an algorithm for curves partial matching up to a similarity transformation is presented the matching is based on the total curvature graph which is the direction of the first derivative of the curve no higher derivatives are required the algorithm uses dynamic programming it is not dependent on reference or feature points the base step consists of a comparison of the total curvatures of neighborhoods the neighborhoods are dynamically chosen according to local characteristics of the curve in the recursive step ratios of lengths of matched curve segments are compared in order to extend the matches the matching is not sensitive to local deviations of the curves the algorithm is fast and simple its worst case time complexity is o mn where m and n are the number of points in the two digital curves experimental results are presented

an energy momentum conserving algorithm using the variational formulation of visco plastic updates in this paper we use the variational formulation of elasto plastic updates proposed by ortiz and stainier comput methods appl mech eng 1999 171 419 444 in the context of consistent time integration schemes we show that such a formulation is well suited to obtain a general expression of energy momentum conserving algorithms moreover we present numerical examples that illustrate the efficiency of our developments copyright 2005 john wiley sons ltd

on the detection of dominant points on digital curves a parallel algorithm is presented for detecting dominant points on a digital closed curve the procedure requires no input parameter and remains reliable even when features of multiple sizes are present on the digital curve the procedure first determines the region of support for each point based on its local properties then computes measures of relative significance e g curvature of each point and finally detects dominant points by a process of nonmaximum suppression this procedure leads to the observation that the performance of dominant points detection depends not only on the accuracy of the measure of significance but also on the precise determination of the region of support this solves the fundamental problem of scale factor selection encountered in various dominant point detection algorithms the inherent nature of scale space filtering in the procedure is addressed and the performance of the procedure is compared to those of several other dominant point detection algorithms using a number of examples

assumed strain finite element methods for conserving temporal integrations in non linear solid dynamics this paper presents a new class of assumed strain finite elements to use in combination with general energy momentum conserving time stepping algorithms so that these conservation properties in time are preserved by the fully discretized system in space and time the case of interest corresponds to nearly incompressible material responses in the fully non linear finite strain elastic and elastoplastic ranges the new elements consider the classical scaling of the deformation gradient with an assumed jacobian its determinant defined locally through a weighted averaging procedure at the element level the key aspect of the newly proposed formulation is the definition of the associated linearized strain operator or b bar operator the developments presented here start by identifying the conditions that this discrete operator must satisfy for the fully discrete system in time and space to inherit exactly the conservation laws of linear and angular momenta and the conservation dissipation law of energy for elastic and inelastic problems respectively care is also taken of the preservation of the relative equilibria and the corresponding group motions associated with the momentum conservation laws and characterized by purely rotational and translational motions superimposed to the equilibrium deformed configuration with these developments at hand a new general b bar operator is introduced that satisfies these conditions the new operator not only accounts for the spatial interpolations e g bilinear displacements with piece wise constant volume but also depends on the discrete structure of the equations in time the aforementioned conservation dissipation properties of energy and momenta are then proven to hold rigorously for the final numerical schemes unconditionally of the time step size and the material model elastic or elastoplastic different finite elements are considered in this framework including quadrilateral and triangular elements for plane problems and brick elements for three dimensional problems several representative numerical simulations are presented involving in particular the use of energy dissipating momentum conserving time stepping schemes recently developed by the author and co workers for general finite strain elastoplasticity in order to illustrate the properties of the new finite elements including these conservation dissipation properties in time and their locking free response in the quasi incompressible case copyright 2007 john wiley sons ltd

a general framework for the numerical solution of problems in finite elasto plasticity this article discusses a general framework for the analysis of initial boundary value problems of rate independent finite elasto plasticity based on the theory of green and naghdi a constitutive model is developed within the context of the above theory employing generalized measures of lagrangian strain and work conjugate measures of stress computational implications of the proposed formulation are discussed in conjunction with an implicit time integrator for the differential algebraic equations of plastic flow representative numerical simulations demonstrate the applicability and predictive capacity of the model in the presence of large plastic deformations

optimum uniform piecewise linear approximation of planar curves two dimensional digital curves are often uniformly approximated by polygons or piecewise linear curves several algorithms have been proposed in the literature to find such curves we present an algorithm that finds a piecewise linear curve with the minimal number of segments required to approximate a curve within a uniform error with fixed initial and final points we compare our optimal algorithm to several suboptimal algorithms with respect to the number of linear segments required in the approximation and the execution time of the algorithm

an explicit solution for implicit time stepping in multiplicative finite strain viscoelasticity abstract we consider the numerical treatment of one of the most popular finite strain models of the viscoelastic maxwell body this model is based on the multiplicative decomposition of the deformation gradient combined with neo hookean hyperelastic relations between stresses and elastic strains the evolution equation is six dimensional and describes an incompressible flow such that the volume changes are purely elastic for the corresponding local initial value problem a fully implicit integration procedure is considered and a simple explicit update formula is derived thus no local iterative procedure is required which makes the numerical scheme more robust and efficient the resulting integration algorithm is unconditionally stable and first order accurate the incompressibility constraint of the inelastic flow is exactly preserved a rigorous proof of the symmetry of the consistent tangent operator is provided moreover some properties of the numerical solution like invariance under the change of the reference configuration and positive energy dissipation within a time step are discussed numerical tests show that in terms of accuracy the proposed integration algorithm is equivalent to the classical implicit scheme based on the exponential mapping finally in order to check the stability of the algorithm numerically a representative initial boundary value problem involving finite viscoelastic deformations is considered a fem solution of the representative problem using msc marc is presented

volume preserving energy momentum schemes for isochoric multiplicative plasticity abstract this paper presents a new energy dissipative momentum conserving algorithm for multiplicative finite strain plasticity f f e f p that also preserves exactly the plastic volume for isochoric plastic models the new algorithm exhibits exactly the conservation laws of linear and angular momentum of the underlying physical problem as well as its energy evolution a strictly positive energy dissipation in fact the exact energy dissipation is obtained by design during plastic steps while enforcing the plastic consistency i e the yield condition on the final stress appearing in the equations of motion exact energy conservation is attained in particular during elastic steps the aforementioned preservation of the plastic volume is obtained by a new treatment of the geometric structure behind the considered multiplicative models of finite strain plasticity namely we present a new approximation of the reference and elastic metrics whose contractions with the incremental total and elastic strains lead exactly to the increment of the total and elastic natural volumetric strains respectively the elastic metric is defined in the intermediate configuration defined itself by the proper discrete approximation in time of the plastic deformation gradient f p the new algorithm extends to the plastic range existing energy momentum conserving schemes for nonlinear elastic problems but incorporating a new modified elastic stress formula consistent with this new geometric setting the inherited conservation laws of momenta and especially the non negative character of the energy dissipation leads to an improved performance over existing more classical schemes showing numerical instabilities in the considered highly nonlinear geometric setting of large deformations and strains several numerical simulations are presented illustrating these properties

using simple decomposition for smoothing and feature point detection of noisy digital curves this correspondence presents an algorithm for smoothed polygonal approximation of noisy digital planar curves and feature point detection the resulting smoothed polygonal representation preserves the signs of the curvature function of the curve the algorithm is based on a simple decomposition of noisy digital curves into a minimal number of convex and concave sections the location of each separation point is optimized yielding the minimal possible distance between the smoothed approximation and the original curve curve points within a convex concave section are discarded if their angle signs do not agree with the section sign and if the resulted deviations from the curve are less than a threshold spl epsi which is derived automatically inflection points are curve points between pairs of convex concave sections and cusps are curve points between pairs of convex convex or concave concave sections corners and points of local minimal curvature are detected by applying the algorithm to respective total curvature graphs the detection of the feature points is based on properties of pairs of sections that are determined in an adaptive manner rather than on properties of single points that are based on a fixed size neighborhood the detection is therefore reliable and robust complexity analysis and experimental results are presented

energy conserving and decaying algorithms in non linear structural dynamics a generalized formulation of the energy momentum methodwill be developed within the framework of the generalizedmethodwhich allows at the same time guaranteed conservation or decay of total energy and controllable numerical dissipation of unwanted high frequency response furthermore the latter algorithm will be extended by the consistently integrated constraints of energy and momentum conservation originally derived for the constraint energy momentum algorithm the goal of this general approach of implicit energyconserving and decaying time integration schemes is to compare these algorithms on the basis of an equivalent notation by the means of an overall algorithmic design and hence to investigate their numerical properties numerical stability and controllable numerical dissipation of high frequencies will be studied in application to non linear structural dynamics among the methods considered will be the newmark method the classical methods the energy momentum methodwith and without numerical dissipation the constraint energymomentum algorithm and the constraint energy method copyright 1999 john wiley sons ltd

a first order energy dissipative momentum conserving scheme for elasto plasticity using the variational updates formulation abstract in a previous paper l noels l stainier j p ponthot an energy momentum conserving algorithm using the variational formulation of visco plastic updates int j numer methods engrg 65 2006 904 942 the authors demonstrated the efficiency of the variational formulation of elasto plastic updates to develop energy momentum conserving time integration algorithms indeed within such a framework the stress tensor always derives from an incremental potential even when plastic behavior is considered therefore the verification of the conservation of energy in the non linear range can easily be demonstrated the sum of the reversible stored energy and irreversible dissipated energy exactly corresponds to the work of the external forces applied to the structure although this formulation was shown to be accurate and robust the introduction of numerical dissipation for high frequency numerical modes can be necessary to simulate complex phenomena in this work we propose a modification of the variational updates framework to introduce this numerical property leading to a new energy dissipative momentum conserving time integration algorithm for elasto plasticity

constraint energy momentum algorithm and its application to non linear dynamics of shells abstract a modified time integration method is developed for non linear structural dynamics combining controllable algorithmic dissipation of higher modes and at the same time conservation of energy as well as linear and angular momentum although these features seem to be actually inconsistent numerical stability with large time steps for dynamic buckling and snap through problems as well as for dynamical systems with smooth solutions is guaranteed arbitrary implicit one step time integration schemes with newmark approximations are suitable for use as basic algorithms for the proposed method the desired conservation attributes are obtained if the basic algorithm is augmented by energy and momentum constraints the potential of the suggested algorithm is demonstrated by selected applications to non linear dynamics of shell structures

a vertex centred finite volume jameson schmidt turkel jst algorithm for a mixed conservation formulation in solid dynamics a vertex centred finite volume algorithm is presented for the numerical simulation of fast transient dynamics problems involving large deformations a mixed formulation based upon the use of the linear momentum the deformation gradient tensor and the total energy as conservation variables is discretised in space using linear triangles and tetrahedra in two dimensional and three dimensional computations respectively the scheme is implemented using central differences for the evaluation of the interface fluxes in conjunction with the jameson schmidt turkel jst artificial dissipation term the discretisation in time is performed by using a total variational diminishing tvd two stage runge kutta time integrator the jst algorithm is adapted in order to ensure the preservation of linear and angular momenta the framework results in a low order computationally efficient solver for solid dynamics which proves to be very competitive in nearly incompressible scenarios and bending dominated applications

the curvature primal sketch in this paper we introduce a novel representation of the significant changes in curvature along the bounding contour of planar shape we call the representation the curvature primal sketch because of the close analogy to the primal sketch representation advocated by marr for describing significant intensity changes we define a set of primitive parameterized curvature discontinuities and derive expressions for their convolutions with the first and second derivatives of a gaussian we describe an implemented algorithm that computes the curvature primal sketch by matching the multiscale convolutions of a shape and illustrate its performance on a set of tool shapes several applications of the representation are sketched

generalized energy momentum method for non linear adaptive shell dynamics in the present study a generalization of the energy momentum method denoted by generalized energy momentum method applied to the non linear dynamics of shells will be developed within the framework of the generalized a method this algorithmic environment contains the unconditionally stable energy momentum method and its numerically damped version as well as the classical newmark and a methods as special cases in order to control the size of the time steps of the integration scheme with respect to accuracy and eaciency an adaptive time stepping procedure based on local a posteriori error estimation will be improved for non linear dynamical systems and applied to the proposed class of algorithms the spatial discretization is realized by an eight noded finite shell element of reissner mindlin type including an extensible shell director field permitting the application of three dimensional material laws the original formulation of this finite element will be developed for non linear dynamic analysis and adapted for the employment within the introduced energy conserving decaying time integration scheme 1999 elsevier science s a all rights reserved

energy momentum consistent finite element discretization of dynamic finite viscoelasticity this paper is concerned with energy momentum consistent time discretizations of dynamic finite viscoelasticity energy consistency means that the total energy is conserved or dissipated by the fully discretized system in agreement with the laws of thermodynamics the discretization is energy momentum consistent if also momentum maps are conserved when group motions are superimposed to deformations the performed approximation is based on a three field formulation in which the deformation field the velocity field and a strain like viscous internal variable field are treated as independent quantities the new non linear viscous evolution equation satisfies a non negative viscous dissipation not only in the continuous case but also in the fully discretized system the initial boundary value problem is discretized by using finite elements in space and time thereby the temporal approximation is performed prior to the spatial approximation in order to preserve the stress objectivity for finite rotation increments incremental objectivity although the present approach makes possible to design schemes of arbitrary order the focus is on finite elements relying on linear lagrange polynomials for the sake of clearness the discrete energy momentum consistency is based on the collocation property and an enhanced second piola kirchhoff stress tensor the obtained coupled non linear algebraic equations are consistently linearized the corresponding iterative solution procedure is associated with newly proposed convergence criteria which take the discrete energy consistency into account the iterative solution procedure is therefore not complicated by different scalings in the independent variables since the motion of the element is taken into account for solving the viscous evolution equation representative numerical simulations with various boundary conditions show the superior stability of the new time integration algorithm in comparison with the ordinary midpoint rule both the quasi rigid deformations during a free flight and large deformations arising in a dynamic tensile test are considered copyright 2009 john wiley sons ltd

constructive learning of recurrent neural networks limitations of recurrent cascade correlation and a simple solution it is often difficult to predict the optimal neural network size for a particular application constructive or destructive methods that add or subtract neurons layers connections etc might offer a solution to this problem we prove that one method recurrent cascade correlation due to its topology has fundamental limitations in representation and thus in its learning capabilities it cannot represent with monotone i e sigmoid and hard threshold activation functions certain finite state automata we give a preliminary approach on how to get around these limitations by devising a simple constructive training method that adds neurons during training while still preserving the powerful fully recurrent structure we illustrate this approach by simulations which learn many examples of regular grammars that the recurrent cascade correlation method is unable to learn

extraction of trend lines and extrema from multiscale curves abstract multiscale planar curves are curves conveying information at several levels of detail to extract polygonal approximations of such curves at different resolutions we propose a parallel method based on connecting locally computed centroids the implementation makes use of a new pyramidal data structure the chain pyramid coarser and coarser approximations are obtained through a hierarchy of increasingly coarse tessellations significant extrema of the curves are localized at the highest resolution by segmenting at a natural scale and employing the structure of the chain pyramid for down projection

a first order hyperbolic framework for large strain computational solid dynamics part i total lagrangian isothermal elasticity abstract this paper introduces a new computational framework for the analysis of large strain fast solid dynamics the paper builds upon previous work published by the authors gil et al 2014 48 where a first order system of hyperbolic equations is introduced for the simulation of isothermal elastic materials in terms of the linear momentum the deformation gradient and its jacobian as unknown variables in this work the formulation is further enhanced with four key novelties first the use of a new geometric conservation law for the co factor of the deformation leads to an enhanced mixed formulation advantageous in those scenarios where the co factor plays a dominant role second the use of polyconvex strain energy functionals enables the definition of generalised convex entropy functions and associated entropy fluxes for solid dynamics problems moreover the introduction of suitable conjugate entropy variables enables the derivation of a symmetric system of hyperbolic equations dual of that expressed in terms of conservation variables third the new use of a tensor cross product 61 greatly facilitates the algebraic manipulations of expressions involving the co factor of the deformation fourth the development of a stabilised petrov galerkin framework is presented for both systems of hyperbolic equations that is when expressed in terms of either conservation or entropy variables as an example a polyconvex mooney rivlin material is used and for completeness the eigen structure of the resulting system of equations is studied to guarantee the existence of real wave speeds finally a series of numerical examples is presented in order to assess the robustness and accuracy of the new mixed methodology benchmarking it against an ample spectrum of alternative numerical strategies including implicit multi field fraeijs de veubeke hu washizu variational type approaches and explicit cell and vertex centred finite volume schemes

a stabilised petrov galerkin formulation for linear tetrahedral elements in compressible nearly incompressible and truly incompressible fast dynamics a mixed second order stabilised petrov galerkin finite element framework was recently introduced by the authors lee et al 2014 46 the new mixed formulation written as a system of conservation laws for the linear momentum and the deformation gradient performs extremely well in bending dominated scenarios even when linear tetrahedral elements are used yielding equal order of convergence for displacements and stresses in this paper this formulation is further enhanced for nearly and truly incompressible deformations with three key novelties first a new conservation law for the jacobian of the deformation is added into the system providing extra flexibility to the scheme second a variationally consistent petrov galerkin stabilisation methodology is derived third an adapted fractional step method is presented for both incompressible and nearly incompressible materials in the context of nonlinear elastodynamics for completeness and ease of understanding these three improvements are presented both in small and large strain regimes studying the eigenstructure of the resulting systems a series of numerical examples are presented in order to demonstrate the robustness of the enhanced methodology with respect to the work previously published by the authors

a consistent time fe method for large strain elasto plasto dynamics in the present paper energy consistent momentum conserving time stepping schemes for geometrically nonlinear multiplicative elasto plasto dynamics referred to as ecmc cg methods are developed in this context the entire discretisation procedure of the global equations of motion is based on galerkin methods in space and time whereby special emphasis is placed on the desired conservation properties related to the approximation of time integrals to guarantee energy consistency also for higher order finite elements in time a thermodynamically motivated enhanced stress tensor which defines a non standard quadrature rule is derived for elasto plastic material behaviour concerning the integration of the local evolution equations in time exemplarily a well established first order accurate exponential update is applied and adapted according to the global time stepping scheme based on linear finite elements in time the performance of the resulting scheme is demonstrated by means of several representative numerical examples whereby it is shown that global energy consistency is guaranteed exactly within the calculation accuracy by applying the proposed concepts

evaluation of interest point detectors many different low level feature detectors exist and it is widely agreed that the evaluation of detectors is important in this paper we introduce two evaluation criteria for interest points repeatability rate and information content repeatability rate evaluates the geometric stability under different transformations information content measures the distinctiveness of features different interest point detectors are compared using these two criteria we determine which detector gives the best results and show that it satisfies the criteria well

an energy consistent material point method for dynamic finite deformation plasticity energy consistency for the material point method mpm is examined for thermodynamically consistent hyperelastic plastic materials it is shown that mpm can be formulated with implicit three field variational finite element algorithms which dissipate energy and conserve momentum for that class of material models with a consistent mass matrix the resulting overall numerical method inherits the energy dissipative and momentum conserving properties of the mesh solution thus the proposed mpm algorithm satisfies by construction a time discrete form of the second law of thermodynamics properties of the method are illustrated in numerical examples

motion estimation of two dimensional objects based on the straight line hough transform a new approach abstract in this paper we propose a novel technique for estimating motion of two dimensional objects the technique involves representing objects uniquely by straight line approximations of the boundary p and estimating motion from shifts in the p space the straight line approximation of the boundary is obtained by detecting dominant points during contour tracing and then approximating the pixels between adjacent dominant points along the contour by a straight line using the hough transform we have developed a computationally efficient way of detecting dominant points in a single pass during contour tracing the proposed object representation enables us to determine the translation and rotation with significantly fewer computations compared with existing methods our method can also detect the motion parameters correctly for partially overlapping objects

a framework for finite strain elastoplasticity based on maximum plastic dissipation and the multiplicative decomposition part ii computational aspects computational aspects of the formulation of elastoplasticity proposed in part i of this work employing a multiplicative decomposition and based on the principle of maximum plastic dissipation are examined in detail in particular the implementation of the proposed hyperelastic j2 flow theory is shown to reduce to a trivial modification of the classical radial return algorithm which is amenable to exact linearization this results in a closed form expression for the consistent elastoplastic tangent moduli the algorithmic treatment of the general case is considered subsequently it is emphasized that objective rates and incrementally objective algorithms play no role in the present approach a number of numerical experiments are presented that illustrate the excellent performance of the proposed formulation

thermodynamically consistent time stepping algorithms for non linear thermomechanical systems we present the basic theory for developing novel monolithic and staggered time stepping algorithms for general non linear coupled thermomechanical problems the proposed methods are thermodynamically consistent in the sense that their solutions rigorously comply with the two laws of thermodynamics for isolated systems they preserve the total energy and the entropy never decreases furthermore if the governing equations of the problem have symmetries the proposed integrators preserve them too the formulation of such methods is based on two ideas expressing the evolution equation in the so called general equations for non equilibrium reversible irreversible coupling format and enforcing from their inception certain directionality and degeneracy conditions on the discrete vector fields the new methods can be considered as an extension of the energy momentum integration algorithms to coupled thermomechanical problems to which they reduce in the purely hamiltonian case in the article the new ideas are applied to a simple coupled problem a double thermoelastic pendulum with symmetry numerical simulations verify the qualitative features of the proposed methods and illustrate their excellent numerical stability which stems precisely from their ability to preserve the structure of the evolution equations they discretize copyright 2009 john wiley sons ltd