Influence of interfacial microcracks on the elastic properties of composites

A model is established to quantify the influence of interfacial microcracks on the elastic properties of a particulate composite using a combination of theoretical and finite element analysis. A unique way to construct physical models which could accommodate both crack size and crack density is proposed. Based on energy principles, the influence of a dilute concentration of interfacial microcracks is first studied. The case of a finite concentration of microcracks is solved subsequently by combining the dilute concentration solutions and the differential scheme. Both cases agreed well with existing composite theories for the limiting condition of complete decohesion. The final model predicts the effective elastic properties as functions of both crack size and microcrack density.     

On upgrading the numerics in combustion chemistry codes

A method of updating and reusing legacy FORTRAN codes for combustion simulations is presented using the DAEPACK software package. The procedure is demonstrated on two codes that come with the CHEMKIN-II package, CONP and SENKIN, for the constant-pressure batch reactor simulation. Using DAEPACK generated code, analytical derivative calculations, sparsity pattern information, and hidden discontinuity information can be obtained for the models of interest. This information can be easily integrated with different solvers giving the modeler great flexibility in selecting the best solution procedure. Using the generated code, the CONP code was connected to three different solvers, and the SENKIN code was connected to two different solvers. The effect of model formulation, analytical derivatives, sparsity, and sensitivity equation solution method were analyzed for three large kinetic mechanisms for methane, acetylene, and n-heptane. For the n-heptane model, with 544 species and 2446 reactions, a factor of 10-speed improvement over the original solution procedure was found using analytical derivatives and sparse linear algebra. For sensitivity calculations, for a small number of parameters, a factor of 55 improvement over the original solution procedure was found for the n-heptane problem. Upon closer examination of results, no one method is found to always be superior to other methods, and selection of the appropriate solution procedure requires an examination of the specific kinetic mechanism, which is easily conducted using DAEPACK generated code.     

Analysis of polar compound classes in SRC-II liquids: Comparison of non-aqueous titrametric,i.r. spectrometric and h.p.l.c. methods

Results for various types of polar compounds in SRC-11 coal-derived liquids and other fuels were obtained by one or more of the three analytical methods: non-aqueous titration, i.r. spectroscopy and h.p.l.c. Practical aspects as well as precision, accuracy and assumptions necessary for effective application of each of the methods are discussed. H.p.l.c. is applicable to the widest variety of compound types, is the most rapid, most sensitive and shows the best promise for increased development and improvement. Application and development of these techniques is a logical step toward improving process monitoring, catalyst development, toxicological screening and general fuel analysis.     

Second-order accurate constraint formulation for subdivision finite element simulation of thin shells

We present a new method for enforcing boundary conditions within subdivision finite element simulations of thin shells. The proposed framework is demonstrated to be second-order accurate with respect to increasing refinement in the displacement and energy norm for simply supported, clamped, free and symmetric boundary conditions. Second-order accuracy on the boundary is consistent with the accuracy of subdivision-based approaches for the interior of a body. Our proposed framework is applicable to both triangular and quadrilateral refinement schemes, and does not impose any topological requirements upon the underlying subdivision control mesh. Several examples from an obstacle course of benchmark problems are used to demonstrate the convergence of the scheme. Copyright © 2004 John Wiley & Sons, Ltd.     

Getting a handle on MATLAB graphics

Among engineers and scientists, MATLAB is one of the most popular computational packages. Part of MATLAB's popularity stems from its simple, yet sophisticated, graphics capabilities. While basic plots are relatively easy to obtain, specialized plots require a little more effort to produce. This article introduces MATLAB's Handle Graphics, which provides a mechanism to fully control and customize graphics objects in MATLAB. With a basic understanding of Handle Graphics, users can produce plots that meet the unique needs and quality standards commonly required by the profession, several examples are presented to illustrate the concepts.