Stiff equation and

The main subroutine for evaluation of the elemental stiffness equations and load vectors. 

After the substitution for T from Equation (2.68), dx from Equation (2.70) and global derivatives of shape functions from Equation (2.71) into the elemental stiffness equation (2,.55) we obtain, for the equation corresponding to N[... 

Note that in contrast to the example shown in Section 2.2,2 the element stiffness equation obtained for this problem is not symmetric. After the substitution for the shape functions and algebraic manipulations... 

After the assembly and insertion of the boundai conditions the following set of global stiffness equations is derived... 

Therefore the second-order derivative of/ appearing in the original form of / is replaced by a term involving first-order derivatives of w and/plus a boundary term. The boundary terms are, normally, cancelled out through the assembly of the elemental stiffness equations over the common nodes on the shared interior element sides and only appear on the outside boundaries of the solution domain. However, as is shown later in this chapter, the appropriate treatment of these integrals along the outside boundaries of the flow domain depends on the prescribed boundary conditions. 

In this section the governing Stokes flow equations in Cartesian, polar and axisymmetric coordinate systems are presented. The equations given in two-dimensional Cartesian coordinate systems are used to outline the derivation of the elemental stiffness equations (i.e. the working equations) of various finite element schemes. 

With respect to the selected elemental and global orders of node numbering the elemental stiffness equations for elements ei, eu and em in Figure 6.4 are expressed as... 

Frontal solution requires very intricate bookkeeping for tracking coefficients and making sure that all of the stiffness equations have been assembled and fully reduced. The process time requirement in frontal solvers is hence larger than a straightforward band solver for equal size problems. 

GAUSSP. Gives Gauss point coordinates and weights required in the numerical integration of the members of the elemental stiffness equations. 

SOLVER Assembles elemental stiffness equations into a banded global matrix, imposes boundary conditions and solves the set of banded equations using the LU decomposition method (Gerald and Wheatley, 1984). SOLVER calls the following 4 subroutines. ... 

ASEMB - reads and assembles the elemental stiffness equations in a banded form, as is illustrated in Figure 7.1 to minimize computer memory requirements. 

THIS SUBROUTINE ASSEIIBLES AND SOLVES GLOBAL STIFFNESS EQUATIONS... 

Tsai and Pagano [2-7] ingeniously recast the stiffness transformation equations to enable ready understanding of the consequences of rotating a lamina in a laminate. By use of various trigonometric identities between sin and cos to powers and sin and cos of multiples of the angle, the transformed reduced stiffnesses. Equation (2.85), can be written as... 

For laminates that are symmetric in both geometry and material properties about the middle surface, the general stiffness equations. Equation (4.24), simplify considerably. That symmetry has the form such that for each pair of equal-thickness laminae (1) both laminae are of the same material properties and principal material direction orientations, i.e., both laminae have the same (Qjjlk and (2) if one lamina is a certain distance above the middle surface, then the other lamina is the same distance below the middle surface. A single layer that straddles the middle surface can be considered a pair of half-thickness laminae that satisfies the symmetry requirement (note that such a lamina is inherently symmetric about the middle surface). ... 

Antisymmetric cross-ply laminates were described in Section 4.3.3 and found to have extensional stiffnesses A. , A. 2, A22 = A.. , and Agg bending-extension coupling stiffnesses B., and 822 =-Bn and bending stiffnesses D., D.,2, 822 = and Dgg. The new terms here in comparison to a speciaily orthotro Dic iaminate are B.,and 822- Because of this coupiing, the three equiiibrium differentiai equations are coupied ... 

Equations (2.22) and (2.23) become indeterminate if ks = k. Special forms are needed for the analytical solution of a set of consecutive, first-order reactions whenever a rate constant is repeated. The derivation of the solution can be repeated for the special case or L Hospital s rule can be applied to the general solution. As a practical matter, identical rate constants are rare, except for multifunctional molecules where reactions at physically different but chemically similar sites can have the same rate constant. Polymerizations are an important example. Numerical solutions to the governing set of simultaneous ODEs have no difficulty with repeated rate constants, but such solutions can become computationally challenging when the rate constants differ greatly in magnitude. Table 2.1 provides a dramatic example of reactions that lead to stiff equations. A method for finding analytical approximations to stiff equations is described in the next section. 

Equation (8.22) for a(0) is also special because, due to symmetry, there is only one adjacent point, a(l). The overall set may be solved by any desired method. Euler s method is discussed below and is illustrated in Example 8.5. There are a great variety of commercial and freeware packages available for solving simultaneous ODEs. Most of them even work. Packages designed for stiff equations are best. The stiffness arises from the fact that VJJ) becomes very small near the tube waU. There are also software packages that will handle the discretization automatically. 

---