locking free isogeometric collocation methods for spatial timoshenko rods in this work we present the application of isogeometric collocation techniques to the solution of spatial timoshenko rods the strong form equations of the problem are presented in both displacement based and mixed formulations and are discretized via nurbs based isogeometric collocation several numerical experiments are reported to test the accuracy and efficiency of the considered methods as well as their applicability to problems of practical interest in particular it is shown that mixed collocation schemes are locking free independently of the choice of the polynomial degrees for the unknown fields such an important property is also analytically proven

an isogeometric collocation approach for bernoulli euler beams and kirchhoff plates abstract in this paper iga collocation methods are for the first time introduced for the solution of thin structural problems described by the bernoulli euler beam and kirchhoff plate models in particular a precise description of the proposed methods of the relevant implementation details and of the strategy to efficiently deal with different combinations of boundary conditions is given finally several numerical experiments confirm that the proposed formulations represent an efficient and geometrically flexible tool for the simulation of thin structures

quadratically convergent direct calculation of critical points for 3d structures undergoing finite rotations in this work we present the implementation details of a quadratically converging newton method based algorithm for direct computation of instability points for 3d structures undergoing finite rotations the structural model chosen for illustration is the 3d geometrically exact beam the proposed algorithm makes use of an extended system where equilibrium equations are supplemented with the loss of stability condition which roughly doubles the total number of equations nonetheless the latter requires only an insignificant increase in computational cost due to judicious use of the bordering algorithm for computing the solution the main thrust of our work is directed towards a careful development of linearized forms of the governing equations employed by newton s method the corresponding results are presented both in material and spatial versions a set of numerical examples is used to illustrate a very satisfying performance of the proposed algorithm

on the fundamental domain of affine springer fibers let g be a connected reductive algebraic group over an algebraically closed field k gamma in mathfrak g k epsilon a semisimple regular element we introduce a fundamental domain f gamma for the affine springer fibers mathscr x gamma we show that the purity conjecture of mathscr x gamma is equivalent to that of f gamma via the arthur kottwitz reduction we then concentrate on the unramified affine springer fibers for the group mathrm gl d it turns out that their fundamental domains behave nicely with respect to the root valuation of gamma we formulate a rationality conjecture about a generating series of their poincare polynomials and study them in detail for the group mathrm gl 3 in particular we pave them in affine spaces and we prove the rationality conjecture

on a non linear formulation for curved timoshenko beam elements considering large displacement rotation increments summary an incremental total lagrangian formulation for curved beam elements that includes the effect of large rotation increments is developed a complete and symmetric tangent stiffness matrix is obtained and the numerical results show in general an improvement over the standard formulation where the assumption of infinitesimal rotation increments is made in the derivation of the tangent stiffness matrix

on the parametrization of finite rotations in computational mechanics a classification of concepts with application to smooth shells this paper concerns the computational treatment of large rotations with application to the finite element discretisation of smooth shells formulations based on various rotational parametrizations are reviewed and classified with respect to their update structure category ii is based on total rotational degrees of freedom leading to an additive update structure whereas category i relies on linearized rotational degrees of freedom leading to a multiplicative update structure based on this classification new formulations are developed among them are the rodrigues formula applied within category ii with additive update structure for the two components of the rotation vector and two successive elementary rotations in the sense of euler angles applied within category i yielding a singularity free formulation the numerical examples confirm that every rotational parametrization when applied within category ii leads to singularities whereas the application within category i enables the calculation of overall rotations unrestricted in size without any singularity

an isogeometric implicit g1 mixed finite element for kirchhoff space rods abstract in the paper we present a multi patch isogeometric mixed element with implicit geometric continuity at the ends for the analysis of non polar kirchhoff space rods rotations are introduced at the ends of the element analogously to the hermitian interpolation performing a change of basis for the configuration space a mixed method in which the internal forces are interpolated by means of an l 2 projection leading to a b formulation is used to cope with the membrane flexural locking that affects the solution the efficiency of the method in avoiding locking and instabilities is investigated the formulation adopted has the advantage with respect to other formulations that a non full global stiffness matrix is obtained this approach appears very appealing for efficient analyses of general assemblies of 3d non polar elements

isogeometric collocation neumann boundary conditions and contact abstract isogeometric collocation methods have been proposed recently and their accuracy and efficiency demonstrated for elastostatics and explicit dynamics this paper addresses two important aspects in the development of the isogeometric collocation technology namely the imposition of neumann boundary conditions and the enforcement of contact constraints which are both treated within the same framework it is shown that the strong imposition of neumann boundary conditions may lead to a significant loss of accuracy in some situations in particular when non uniform meshes are used two possible remedies are proposed to restore the desired level of accuracy while keeping the computational cost virtually unchanged i e a hybrid collocation galerkin approach and an enhanced collocation ec approach a frictionless contact formulation suitable for the collocation framework is further proposed and shown to pass the contact patch test to machine precision when combined with the ec approach the formulation is shown to deliver accurate results and to perform robustly also for highly non uniform meshes for all the collocation formulations contact pressures are greater than or equal to zero pointwise in contrast with standard lagrange finite elements

isogeometric collocation for elastostatics and explicit dynamics we extend the development of collocation methods within the framework of isogeometric analysis iga to multi patch nurbs configurations various boundary and patch interface conditions and explicit dynamic analysis the methods developed are higher order accurate stable with no hourglass modes and efficient in that they require a minimum number of quadrature evaluations the combination of these attributes has not been obtained previously within standard finite element analysis

critical study of newmark scheme on manifold of finite rotations in this paper we show that the usual material form of the newmark time stepping scheme of finite rotations is only a simplified version of the correct formula this is because spin material rotation vectors angular velocity vectors and angular accelerations vectors belong to different tangential vector spaces of a manifold at separate time moments we give corrected newmark time stepping schemes for material description and for spatial description

nonlinear isogeometric spatial bernoulli beam abstract a new element formulation of a geometrically nonlinear spatial curved beam assuming bernoulli theory including torsion without warping is proposed the element formulation is derived directly from the 3d continuum by means of nonlinear kinematics thus accounting for large displacements the geometric description of the proposed element is adapted from the spatial rod of greco and cuomo 2013 and extended to a nonlinear element formulation the proposed formulation can handle arbitrary orientations of the cross section along the beam in this publication nurbs are used as basis functions for discretization since they can easily provide the required c 1 continuity between elements the presented element formulation has four degrees of freedom three for displacements and one for the rotation around the center line in order to prove the accuracy of the developed spatial bernoulli beam several numerical examples are presented and compared to analytic solutions and other element formulations

non semisimple and complex gaugings of n 16 supergravity maximal and non maximal supergravities in three dimensions allow for a large variety of semisimple chern simons gauge groups in this paper we analyze non semisimple and complex gauge groups that satisfy the pertinent consistency relations for a maximal n 16 gauged supergravity to exist we give a general procedure how to generate non semisimple gauge groups from known admissible semisimple gauge groups by a singular boost within e8 8 examples include the theories with gauge group so 8 t28 that describe the reduction of iia iib supergravity on the seven sphere in addition we exhibit two strange embeddings of the complex gauge group open image in new window into real e8 8 and prove that both can be consistently gauged we discuss the structure of the associated scalar potentials as well as their relation to those of d 4 gauged supergravities

a consistent co rotational formulation for non linear three dimensional beam elements abstract this paper describes a co rotational formulation for three dimensional beams in which both the internal force vector and tangent stiffness matrix are consistently derived from the adopted strain measures the latter relate to standard beam theory but are embedded in a continuously rotating frame a set of numerical examples show that the element provides an excellent numerical performance

selective and reduced numerical integrations for nurbs based isogeometric analysis we propose a new approach to construct selective and reduced integration rules for isogeometric analysis based on nurbs elements the notion of an approximation space that approximates the target space is introduced we explore the use of various approximation spaces associated with optimal patch wise numerical quadratures that exactly integrate the polynomials in approximation spaces with the minimum number of quadrature points patch rules exploit the higher continuity of spline basis functions the tendency of smooth spline functions to exhibit numerical locking in nearly incompressible problems when using a full gauss legendre quadrature is alleviated with selective or reduced integration stability and accuracy of the schemes are examined analyzing the discrete spectrum in a generalized eigenvalue problem we propose a local algorithm which is robust and computationally efficient to compute element by element the quadrature points and weights in patch rules the performance of the methods is assessed on several numerical examples in two dimensional elasticity and reissner mindlin shell structures

isogeometric collocation methods for the reissner mindlin plate problem abstract within the general framework of isogeometric methods collocation schemes have been recently proposed as a viable and promising low cost alternative to standard isogeometric galerkin approaches in this paper isogeometric collocation methods for the numerical approximation of reissner mindlin plate problems are proposed for the first time locking free primal and mixed formulations are herein considered and the potential of isogeometric collocation as a geometrically flexible and computationally efficient simulation tool for shear deformable plates is shown through the solution of several numerical tests

single variable formulations and isogeometric discretizations for shear deformable beams abstract we present numerical formulations of timoshenko beams with only one unknown the bending displacement and it is shown that all variables of the beam problem can be expressed in terms of it and its derivatives we develop strong and weak forms of the problem the strong form of the problem involves the fourth derivative of the bending displacement whereas the symmetric weak form involves somewhat surprisingly third and second derivatives based on these we develop isogeometric collocation and galerkin formulations that are completely locking free and involve only half the degrees of freedom compared to standard timoshenko beam formulations several numerical tests are presented to demonstrate the performance of the proposed formulations

on the role of frame invariance in structural mechanics models at finite rotations in this work we re examine so called geometrically exact models for structures such as beams shells or solids with independent rotation field with respect to invariance under superposed rigid body motion a special attention is given to clarifying the issues pertaining to the finite element implementation which guarantees that the invariance of the continuum problem is preserved by its discrete approximation several numerical simulations dealing with finite rotations of structural models are presented in order to further illustrate and confirm the given theoretical considerations

consistent quaternion interpolation for objective finite element approximation of geometrically exact beam we explore an isoparametric interpolation of total quaternion for geometrically consistent strain objective and path independent finite element solutions of the geometrically exact beam this interpolation is a variant of the broader class known as slerp the equivalence between the proposed interpolation and that of relative rotation is shown without any recourse to local bijection between quaternions and rotations we show that for a two noded beam element the use of relative rotation is not mandatory for attaining consistency cum objectivity and an appropriate interpolation of total rotation variables is sufficient the interpolation of total quaternion which is computationally more efficient than the one based on local rotations converts nodal rotation vectors to quaternions and interpolates them in a manner consistent with the character of the rotation manifold this interpolation unlike the additive interpolation of total rotation corresponds to a geodesic on the rotation manifold for beam elements with more than two nodes however a consistent extension of the proposed quaternion interpolation is difficult alternatively a quaternion based procedure involving interpolation of relative rotations is proposed for such higher order elements we also briefly discuss a strategy for the removal of possible singularity in the interpolation of quaternions proposed in i romero the interpolation of rotations and its application to finite element models of geometrically exact rods comput mech 34 2004 121 133 the strain objectivity and path independence of solutions are justified theoretically and then demonstrated through numerical experiments this study being focused only on the interpolation of rotations uses a standard finite element discretization as adopted by simo and vu quoc j c simo l vu quoc a three dimensional finite rod model part ii computational aspects comput methods appl mech engrg 58 1986 79 116 the rotation update is achieved via quaternion multiplication followed by the extraction of the rotation vector nodal rotations are stored in terms of rotation vectors and no secondary storages are required

b spline interpolation of kirchhoff love space rods the paper deals with the isogeometric analysis via b splines of space rods under kirchhoff love hypotheses the approach was used by gontier and vollmer 12 for developing a plane curve element within the framework of the timoshenko rod model but they adopted only one patch to represent entirely the geometry of the rod furthermore the authors developed their theory only for plane elements in this work we develop an isogeometric approach for the numerical analysis of the 3d kirchhoff love rod theory we use b splines and bezier interpolations and we show that they are able to attain very good accuracy for rod structures particularly for developing a 3d exact curve element with geometric torsion the paper presents an original parametrization of the geometric torsion that proves to be very effective the use of b splines allows to avoid discontinuities on the geometrical quantities and particularly on the normal fields so that even relatively low order interpolation functions are able to yield accurate results

on the differentiation of the rodrigues formula and its significance for the vector like parameterization of reissner simo beam theory in this paper we present a systematic way of differentiating up to the second directional derivative i the rodrigues formula and ii the spin rotation vector variation relationship to achieve this goal several trigonometric functions are grouped into a family of scalar quantities which can be expressed in terms of a single power series these results are then applied to the vector like parameterization of reissner simo beam theory enabling a straightforward derivation and leading to a clearer formulation in particular and in contrast with previous formulations a relatively compact and obviously symmetric form of the tangent operator is obtained the paper also discusses several relevant issues concerning a beam finite element implementation and concludes with the presentation of a few selected illustrative examples copyright 2002 john wiley sons ltd

computational aspects of vector like parametrization of three dimensional finite rotations theoretical and computational aspects of vector like parametrization of three dimensional finite rotations which uses only three rotation parameters are examined in detail in this work the relationship of the proposed parametrization with the intrinsic representation of finite rotations via an orthogonal matrix is clearly identified careful considerations of the consistent linearization procedure pertinent to the proposed parametrization of finite rotations are presented for the chosen model problem of reissner s non linear beam theory pertaining details of numerical implementation are discussed for the simplest choice of the finite element interpolations for a 2 node three dimensional beam element a number of numerical simulations in three dimensional finite rotation analysis are presented in order to illustrate the proposed approach

geometrically exact 3d beam theory implementation of a strain invariant finite element for statics and dynamics geometrically exact 3d beam theory has been used as a basis for development of a variety of finite element formulations it has recently become apparent that the important requirement of objectivity of adopted strain measures although provided by the theory itself does not automatically extend to a finite element formulation in this paper we present a new finite element formulation of the geometrically exact 3d beam theory specifically designed to preserve the objectivity of the adopted strain measures in order to do so the current local rotations are interpolated in a manner similar to that adopted in co rotational approaches however no approximations typical for co rotational approaches are introduced into the procedure so in contrast to co rotational formulations the present formulation fully preserves the geometric exactness of the theory a range of numerical examples serves to illustrate the problem and to assess the formulation

on the geometrically exact beam model a consistent effective and simple derivation from three dimensional finite elasticity the goal of the present work is to emphasize some intrinsic features of the three dimensional beam model originally proposed by simo simo j c 1985 a finite strain beam formulation the three dimensional dynamic problem part i comput methods appl mech eng 49 55 70 and to derive the model equations for the general finite deformation case as well as for the finite deformation small strain case in a consistent but simpler way with respect to what it is generally done in the literature

an isogeometric collocation method using superconvergent points abstract we develop an iga collocation method modified by collocating at points other than the standard greville abscissae the method is related to orthogonal collocation used for solving differential equations and to the superconvergence theory therefore we refer to this method as super collocation iga sc by carefully choosing the collocation points it can be seen that the iga sc converges in the first derivative energy norms at rates similar to that of the galerkin solution this is different from the collocation at greville abscissae iga c where the convergence in energy norm for odd polynomial degrees is typically suboptimal the method is tested on 1d 2d and 3d numerical examples in which it is compared to iga c and galerkin s method iga g the comparison includes a detailed cost vs accuracy analysis which shows an improved efficiency of the proposed method in particular for odd polynomial degrees

interpolation of rotational variables in nonlinear dynamics of 3d beams the formulation of dynamic procedures for three dimensional 3 d beams requires extensive use of the algebra pertaining to the non linear character of the rotation group in space the corresponding extraction procedure to obtain the rotations that span a time step has certain limitations which can have a detrimental eect on the overall stability of a time integration scheme the paper describes two algorithms for the dynamics of 3 d beams which dier in their manifestation of the above limitation the rst has already been described in the literature and involves the interpolation of iterative rotations while an alternative formulation which eliminates the above eect by design requires interpolation of incremental rotations theoretical arguments are backed by numerical results similarities between the conventional and new formulation are pointed out and are shown to be big enough to enable easy transformation of the conventional formulation into the new one 1998 john wiley sons ltd

improved numerical integration for locking treatment in isogeometric structural elements part i beams a general mathematical framework is proposed in this paper to define new quadrature rules in the context of bb spline nurbs based isogeometric analysis high order continuity across the elements within a patch turned out to have higher accuracy than c0c0 finite elements as well as a better time efficiency unfortunately a maximum regularity accentuates the shear and membrane locking in thick structural elements the improved selective reduced integration schemes are given for uni dimensional beam problems with basis functions of order two and three and can be easily extended to higher orders the resulting bb spline nurbs finite elements are free from membrane and transverse shear locking moreover no zero energy modes are generated the performance of the approach is evaluated on the classical test of a cantilever beam subjected to a distributed moment and compared to lagrange under integrated finite elements

objectivity of strain measures in the geometrically exact three dimensional beam theory and its finite element implementation the paper discusses the issue of discretization of the strainconfiguration relationships in the geometrically exact theory of threedimensional 3d beams which has been at the heart of most recent

on finite element implementation of geometrically nonlinear reissner s beam theory three dimensional curved beam elements abstract finite element implementation of the three dimensional finite strain beam theory of reissner is considered in this work in contrast with some earlier works on the subject discussed here are the beam elements whose reference axes are arbitrary space curved lines we have shown that an improved representation of curved reference geometry significantly increases the accuracy of the results however it also makes the choice of non locking finite element interpolations somewhat more delicate a hierarchical displacement interpolation proposed here is proved to be capable of eliminating both shear and membrane locking phenomena a very satisfying non locking performance is demonstrated for a set of problems in nonlinear elasto statics

studies of refinement and continuity in isogeometric structural analysis we investigate the effects of smoothness of basis functions on solution accuracy within the isogeometric analysis framework we consider two simple one dimensional structural eigenvalue problems and two static shell boundary value problems modeled with trivariate nurbs solids we also develop a local refinement strategy that we utilize in one of the shell analyses we find that increased smoothness that is the k method leads to a significant increase in accuracy for the problems of structural vibrations over the classical c0c0 continuous p method whereas a judicious insertion of c0c0 continuous surfaces about singularities in a mesh otherwise generated by the k method usually outperforms a mesh in which all basis functions attain their maximum level of smoothness we conclude that the potential for the k method is high but smoothness is an issue that is not well understood due to the historical dominance of c0c0 continuous finite elements and therefore further studies are warranted

numerical efficiency locking and unlocking of nurbs finite elements abstract the strategy of using finite elements with nurbs shape functions for approximation of both geometry and displacements isogeometric approach is investigated from the point of view of finite element technology convergence rates are compared to those of classical finite element approaches utilizing standard lagrange shape functions moreover typical locking phenomena are examined it is found that higher order inter element continuity within the nurbs approach results in identical convergence rates but smaller absolute errors compared to c 0 continuous approaches however nurbs finite elements suffer from the same locking problems as finite elements using lagrange shape functions the discrete shear gap dsg method a general framework for formulation of locking free elements is applied to develop a new class of nurbs finite elements the resulting nurbs dsg elements are absolutely free from locking and preserve the property of improved accuracy compared with standard locking free finite elements the method is exemplified for the timoshenko beam model but may be applied to more general cases

accurate efficient and iso geometrically flexible collocation methods for phase field models we propose new collocation methods for phase field models our algorithms are based on isogeometric analysis a new technology that makes use of functions from computational geometry such as for example non uniform rational b splines nurbs nurbs exhibit excellent approximability and controllable global smoothness and can represent exactly most geometries encapsulated in computer aided design cad models these attributes permitted us to derive accurate efficient and geometrically flexible collocation methods for phase field models the performance of our method is demonstrated by several numerical examples of phase separation modeled by the cahn hilliard equation we feel that our method successfully combines the geometrical flexibility of finite elements with the accuracy and simplicity of pseudo spectral collocation methods and is a viable alternative to classical collocation methods

locking free isogeometric collocation formulation for three dimensional geometrically exact shear deformable beams with arbitrary initial curvature abstract this paper presents a displacement based and a mixed isogeometric collocation formulations for three dimensional geometrically exact shear deformable beams with arbitrarily curved initial geometry strong forms of both formulations are consistently linearized with respect to the underlying configuration manifold high efficiency is achieved by using a geometrically consistent kinematic model in which incremental rotations are parametrized through the spatial rotation vector several numerical tests demonstrate that the proposed mixed formulation is totally locking free for any choice of approximation degree mixed formulation also exhibits higher robustness in the incremental iterative algorithm elimination of shear locking and low computational cost assured by the combination of isogeometric collocation with a consistent kinematic model result in a very efficient accurate and geometrically powerful method able to handle geometrically complex beam problems

isogeometric collocation for large deformation elasticity and frictional contact problems abstract isogeometric collocation methods have been recently proposed as an alternative to standard galerkin approaches as they provide a significant reduction in computational cost for higher order discretizations in this work we explore the application of isogeometric collocation to large deformation elasticity and frictional contact problems we first derive the non linear governing equations for the elasticity problem with finite deformation kinematics and provide details on their consistent linearization some numerical examples demonstrate the performance of collocation in its basic and enhanced versions differing by the enforcement of neumann boundary conditions for problems with strong singularities enhanced collocation is shown to outperform basic collocation and to lead to a spatial convergence behavior very similar to galerkin whereas for weaker or no singularities enhanced and basic collocation may give very similar results a large deformation contact formulation is subsequently developed and tested in the frictional setting where collocation confirms the excellent performance already obtained for the frictionless case finally it is shown that the contact formulation in the collocation framework passes the contact patch test to machine precision in a three dimensional setting with arbitrarily inclined non matching discretizations thus outperforming most of the available contact formulations and all those with pointwise enforcement of the contact constraints

an excursion into large rotations the present discourse develops an enlarged exploration of the matrix formulation of finite rotations in space initiated in 1 it is shown how a consistent but subtle matrix calculus inevitably leads to a number of elegant expressions for the transformation or rotation matrix t appertaining to a rotation about an arbitrary axis also analysed is the case of multiple rotations about fixed or follower axes particular attention is paid to an explicit derivation of a single compound rotation vector equivalent to two consecutive arbitrary rotations this theme is discussed in some detail for a number of cases semitangential rotations for which commutativity holds first proposed in 2 3 are also considered furthermore an elementary geometrical analysis of large rotations is also given finally we deduce in an appendix using a judicious reformulation of quarternions the compound pseudovector representing the combined effect of n rotations in the author s opinion the present approach appears preferable to a pure vectorial scheme and even more so to an indicial formulation and is computationally more convenient

avoiding shear locking for the timoshenko beam problem via isogeometric collocation methods in this work we study isogeometric collocation methods for the timoshenko beam problem considering both mixed and displacement based formulations in particular we show that locking free solutions are obtained for mixed methods independently on the approximation degrees selected for the unknown fields moreover several numerical tests are provided in order to support our theoretical results and to show the good behavior and the flexibility of isogeometric collocation methods in this context

efficient quadrature for nurbs based isogeometric analysis abstract we initiate the study of efficient quadrature rules for nurbs based isogeometric analysis a rule of thumb emerges the half point rule indicating that optimal rules involve a number of points roughly equal to half the number of degrees of freedom or equivalently half the number of basis functions of the space under consideration the half point rule is independent of the polynomial order of the basis efficient rules require taking into account the precise smoothness of basis functions across element boundaries several rules of practical interest are obtained and a numerical procedure for determining efficient rules is presented we compare the cost of quadrature for typical situations arising in structural mechanics and fluid dynamics the new rules represent improvements over those used previously in isogeometric analysis

total lagrangian reissner s geometrically exact beam element without singularities in this paper we introduce a new reissner s geometrically exact beam element which is based on a total lagrangian updating procedure the element has the rotation vector as the dependent variable and the singularity problems at the rotation angle 2 and its multiples are passed by the change of parametrization on the rotation manifold the beam formulation has several benefits such as all the unknown vectors belong to the same tangential vector space no need for secondary storage variables the path independence in the static case any standard time integration algorithm may be used and the symmetric stiffness copyright 2006 john wiley sons ltd

an extended isogeometric thin shell analysis based on kirchhoff love theory abstract an extended isogeometric element formulation xiga for analysis of through the thickness cracks in thin shell structures is developed the discretization is based on non uniform rational b splines nurbs the proposed xiga formulation can reproduce the singular field near the crack tip and the discontinuities across the crack it is based on the kirchhoff love theory where c 1 continuity of the displacement field is required this condition is satisfied by the nurbs basis functions hence the formulation eliminates the need of rotational degrees of freedom or the discretization of the director field facilitating the enrichment strategy the performance and validity of the formulation is tested by several benchmark examples

isogeometric collocation cost comparison with galerkin methods and extension to adaptive hierarchical nurbs discretizations abstract we compare isogeometric collocation with isogeometric galerkin and standard c 0 finite element methods with respect to the cost of forming the matrix and residual vector the cost of direct and iterative solvers the accuracy versus degrees of freedom and the accuracy versus computing time on this basis we show that isogeometric collocation has the potential to increase the computational efficiency of isogeometric analysis and to outperform both isogeometric galerkin and standard c 0 finite element methods when a specified level of accuracy is to be achieved with minimum computational cost we then explore an adaptive isogeometric collocation method that is based on local hierarchical refinement of nurbs basis functions and collocation points derived from the corresponding multi level greville abscissae we introduce the concept of weighted collocation that can be consistently developed from the weighted residual form and the two scale relation of b splines using weighted collocation in the transition regions between hierarchical levels we are able to reliably handle coincident collocation points that naturally occur for multi level greville abscissae the resulting method combines the favorable properties of isogeometric collocation and hierarchical refinement in terms of computational efficiency local adaptivity robustness and straightforward implementation which we illustrate by numerical examples in one two and three dimensions

an objective finite element approximation of the kinematics of geometrically exact rods and its use in the formulation of an energy momentum conserving scheme in dynamics we present in this paper a new finite element formulation of geometrically exact rod models in the three dimensional dynamic elastic range the proposed formulation leads to an objective or frame indifferent under superposed rigid body motions approximation of the strain measures of the rod involving finite rotations of the director frame in contrast with some existing formulations this goal is accomplished through a direct finite element interpolation of the director fields defining the motion of the rod s cross section furthermore the proposed framework allows the development of time stepping algorithms that preserve the conservation laws of the underlying continuum hamiltonian system the conservation laws of linear and angular momenta are inherited by construction leading to an improved approximation of the rod s dynamics several numerical simulations are presented illustrating these properties copyright 2002 john wiley sons ltd

isogeometric large deformation frictionless contact using t splines t spline based isogeometric analysis is applied to frictionless contact problems between deformable bodies in the context of large deformations the continuum is discretized with cubic t splines and cubic nurbs a gauss point to surface formulation is combined with the penalty method to treat the contact constraints in the discretized setting it is demonstrated that analysis suitable t splines coupled with local refinement accurately approximate contact pressures with far fewer degrees of freedom than nurbs both two and three dimensional examples are presented additionally all t spline analysis models are generated using commercially available t spline modeling software without intermediate mesh generation or geometry clean up steps

on the choice of finite rotation parameters in this work we discuss some aspects of the three dimensional finite rotations pertinent to the formulation and computational treatment of the geometrically exact structural theories among various possibilities to parameterize the finite rotations special attention is dedicated to a choice featuring an incremental rotation vector some computational aspects pertinent to the implementation of the newton iterative scheme and the newmark time stepping algorithm applied to solving these problems are examined representative numerical simulations are presented in order to illustrate the performance of the proposed formulation

a simple algorithm for obtaining nearly optimal quadrature rules for nurbs based isogeometric analysis we develop new quadrature rules for isogeometric analysis based on the solution of a local nonlinear problem a simple and robust algorithm is developed to determine the rules which are exact for important b spline spaces of uniform and geometrically stretched knot spacings we consider both periodic and open knot vector configurations and illustrate the efficiency of the rules on selected boundary value problems we find that the rules are almost optimally efficient but much easier to obtain than optimal rules which require the solution of global nonlinear problems that are often ill posed