Step Function Representation of Solid Models and Application to Mesh Free Engineering Analysis

2005 ◽  
Vol 128 (1) ◽  
pp. 46-56 ◽  
Author(s):  
Ashok V. Kumar ◽  
Jongho Lee

Numerical methods for solving boundary value problems that do not require generation of mesh to approximate the analysis domain have been referred to as mesh-free methods. While many of these are “mesh less” methods that do not have connectivity between nodes, a subset of these methods uses a structured mesh or grid for the analysis that does not conform to the geometry of the domain of analysis. Instead the geometry is represented using implicit equations. In this paper we present a method for constructing step functions of solids whose boundaries are represented using implicit equations. Step functions can be used to compute volume integrals over the solid that are needed for mesh free analysis. The step function of the solid has a unit value within the solid and zero outside. A level set of this step function can then be defined as the boundary of the solid. Boolean operators are defined in this paper that enable step functions of half-spaces and primitives to be combined to construct a single step function for more complex solids. Application of step functions to analysis using nonconforming mesh is illustrated.

Author(s):  
Ashok V. Kumar ◽  
Jongho Lee

Implicit equations of curves and surfaces have been shown to be useful for constructing solutions for boundary value problems such that the boundary conditions are satisfied exactly. This application has generated interest in constructing solid models where the geometry is represented using implicit equations rather than parametric equations. In this paper we present a method for constructing step functions of solids that have a unit value within the solid and zero outside. A level set of this step function can then be defined as the boundary of the solid. This step function can be used not only to apply boundary conditions but also to compute volume integrals over the solid. Methods for combining step functions of solid primitives using ordinary and regularized Boolean operations to construct step functions of the Boolean result are also presented.


Author(s):  
Ashok V. Kumar ◽  
Jongho Lee ◽  
Ravi Burla

In traditional solid modeling the boundaries of the solid are represented using parametric equations. Even though the application of implicit equations has also been explored, they have not been widely used. Interest has been rekindled recently due to application of implicit equations to mesh free engineering analysis. In this paper, an implicit representation scheme for solids is presented where the boundaries of primitive solids are defined using implicit equation of surfaces. To ensure that the equations are axis independent, the characteristic functions for the implicit equations are defined by interpolating within hexahedral elements. Primitive solids are defined by sweeping closed 2D profiles. The boundaries of these profiles are defined using implicit equations of curves. Implicit equations can be used for constructing “step function” of the primitives and their Boolean combinations. The step functions of a solid has a unit value inside the solid and zero outside and can be used for computing volume integrals needed for mesh free analysis.


Author(s):  
Linxia Gu ◽  
Ashok V. Kumar

One of the main advantages of meshless methods is that it eliminates the mesh generation, but it is still necessary to place nodes with controlled spacing variation on the boundary and within the domain. However, due to lack of connectivity between nodes it is more difficult to interpolate the field variables and impose boundary conditions. In this paper, a mesh free method is presented for analysis using a structured grid that does not conform to the geometry of the domain. The geometry of the domain is independent of the structured grid and is represented using implicit equations. The implicit equations of the boundaries can be used to construct solution structures that satisfy boundary conditions exactly even though the nodes of the grid are not on the boundaries of the domain. The solution structures are constructed using the implicit equations of the boundary together with a piece-wise interpolation over the structured grid. The implicit equations are also used to construct step function of solid such that its value is equal to unity inside the solid and zero outside. The step function of the solid is used for volume integrations needed for the analysis. The traditional weak form for Poisson’s equation is modified by using this solution structure to eliminate the surface integration terms. The accuracy and implementation of the present mesh free method is illustrated for two-dimensional heat conduction problems governed by Poisson’s equation. Satisfactory results are obtained when compared with analytical results and results from commercial finite element software.


Author(s):  
Teije de Jong

AbstractIn this series of papers I attempt to provide an answer to the question how the Babylonian scholars arrived at their mathematical theory of planetary motion. Papers I and II were devoted to system A theory of the outer planets and of the planet Venus. In this third and last paper I will study system A theory of the planet Mercury. Our knowledge of the Babylonian theory of Mercury is at present based on twelve Ephemerides and seven Procedure Texts. Three computational systems of Mercury are known, all of system A. System A1 is represented by nine Ephemerides covering the years 190 BC to 100 BC and system A2 by two Ephemerides covering the years 310 to 290 BC. System A3 is known from a Procedure Text and from Text M, an Ephemeris of the last evening visibility of Mercury for the years 424 to 403 BC. From an analysis of the Babylonian observations of Mercury preserved in the Astronomical Diaries and Planetary Texts we find: (1) that dates on which Mercury reaches its stationary points are not recorded, (2) that Normal Star observations on or near dates of first and last appearance of Mercury are rare (about once every twenty observations), and (3) that about one out of every seven pairs of first and last appearances is recorded as “omitted” when Mercury remains invisible due to a combination of the low inclination of its orbit to the horizon and the attenuation by atmospheric extinction. To be able to study the way in which the Babylonian scholars constructed their system A models of Mercury from the available observational material I have created a database of synthetic observations by computing the dates and zodiacal longitudes of all first and last appearances and of all stationary points of Mercury in Babylon between 450 and 50 BC. Of the data required for the construction of an ephemeris synodic time intervals Δt can be directly derived from observed dates but zodiacal longitudes and synodic arcs Δλ must be determined in some other way. Because for Mercury positions with respect to Normal Stars can only rarely be determined at its first or last appearance I propose that the Babylonian scholars used the relation Δλ = Δt −3;39,40, which follows from the period relations, to compute synodic arcs of Mercury from the observed synodic time intervals. An additional difficulty in the construction of System A step functions is that most amplitudes are larger than the associated zone lengths so that in the computation of the longitudes of the synodic phases of Mercury quite often two zone boundaries are crossed. This complication makes it difficult to understand how the Babylonian scholars managed to construct System A models for Mercury that fitted the observations so well because it requires an excessive amount of computational effort to find the best possible step function in a complicated trial and error fitting process with four or five free parameters. To circumvent this difficulty I propose that the Babylonian scholars used an alternative more direct method to fit System A-type models to the observational data of Mercury. This alternative method is based on the fact that after three synodic intervals Mercury returns to a position in the sky which is on average only 17.4° less in longitude. Using reduced amplitudes of about 14°–25° but keeping the same zone boundaries, the computation of what I will call 3-synarc system A models of Mercury is significantly simplified. A full ephemeris of a synodic phase of Mercury can then be composed by combining three columns of longitudes computed with 3-synarc step functions, each column starting with a longitude of Mercury one synodic event apart. Confirmation that this method was indeed used by the Babylonian astronomers comes from Text M (BM 36551+), a very early ephemeris of the last appearances in the evening of Mercury from 424 to 403 BC, computed in three columns according to System A3. Based on an analysis of Text M I suggest that around 400 BC the initial approach in system A modelling of Mercury may have been directed towards choosing “nice” sexagesimal numbers for the amplitudes of the system A step functions while in the later final models, dating from around 300 BC onwards, more emphasis was put on selecting numerical values for the amplitudes such that they were related by simple ratios. The fact that different ephemeris periods were used for each of the four synodic phases of Mercury in the later models may be related to the selection of a best fitting set of System A step function amplitudes for each synodic phase.


Author(s):  
Shinya IMAYASU ◽  
Matthias ROTHLIN ◽  
Mansur AKBARI ◽  
Nikolas SCHAAL ◽  
Konrad WEGENER

Author(s):  
S. F. Lukomskii ◽  
G. S. Berdnikov

We consider a class of (N, M)-elementary step functions on the p-adic Vilenkin group. We prove that (N, M)-elementary step function generates a MRA on p-adic Vilenkin group if and only if it is generated by a special N-valid rooted tree on the set of vertices {0,1,…p - 1} with the vector (0,…,0) ∈ ℤN as a root.


Sign in / Sign up

Export Citation Format

Share Document