<< Chapter < Page Chapter >> Page >

Many of the concepts from the single-string case carry over to networks.

Assembling theAndMatrices

We begin by describing the notation of the information represented by our data structures. We denote the i th node to have x y z position vector p i . For each node at p i , there is an associated set of the indices of connected neighbor nodes N i and a set C i containing the physical constants k i j , s i j pertaining to the connection between nodes i and j N i . Since we can divide through by ρ i on both sides of the network wave equation, we can assume without loss of generality that the constant ρ i = 1 , and that any data carried by the density ρ i is now contained in k i j .

Assuming we are given a set of N nodes, along with the x y positions of each node (the z positions are assumed to be 0, such that the web is planar in the x y plane at rest), our first goal is to compute our step sizes h i j and orientation vectors v i j for connections between two nodes i and j N i . To account for Dirichlet boundary conditions, we also create an anchored node for each endpoint. In this implementation, if a node has only one neighbor, we assume it is connected to a pinned endpoint whose position is in the opposite direction but the same distance away as the only neighbor (this is required to calculate an inner product). For a node connected to an endpoint, we append the index 0 to N i .

Given p i , N i , C i , we proceed as follows

  • for each i from 1 , 2 , ... , N
    • for all j N i , calculate h i j = p i - p j 2
    • if the number of elements in any set N i = 1
      • create an“endpoint" node at position p 0 = 2 p i - p j , with step size h i 0 = p i - p 0 2
      • set N i = { N i , 0 } to indicate the i th node is connected to an anchor

In practice, we normalize the positions of our nodes such that the web lies within a box of a desired arbitrary size s . We do this by calculating the maximum distance d max between the anchored endpoints of a web, then scaling the positions p i of every node by the factor s / d max . Since the absolute positions of the nodes don't affect our discretization, we don't need to worry about subtracting off the centroid of all the node positions.

With all our variables now in place, we can now proceed to the actual construction of our discretization matrices. This requires knowing φ i , φ j , a φ i , φ j , and P i j . We deal first with constructing the M matrix, which requires only knowledge of φ i , φ j . Just as in the case of finite element on the single string, most of the basis functions φ i and φ j don't share support and their inner products are zero. However, in addition to calculating inner products of regular hat functions, we need to compute the inner products with joint, generalized and nonsymmetric hat functions as well.

Starting with our M matrix, we only need to calculate the inner product φ i , φ j for j = i and j N i . For the diagonal case j = i , we note that our inner product u , v = i = 1 L 0 i u ( x i ) v ( x i ) d x i needs only be calculated on the support of φ i . For given i ,

φ i , φ i = j N i p i p i + h j φ i 2 ( x ) d x i = j N i h i j 3

The last part is a generalization of our inner product for a uniform grid on a single string. For the off-diagonal case j N i , j i , the inner product is analogous to the single-string case.

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, The art of the pfug. OpenStax CNX. Jun 05, 2013 Download for free at http://cnx.org/content/col10523/1.34
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'The art of the pfug' conversation and receive update notifications?

Ask