Engineering math software and FEA make a potent team for solving tough problems and a way to clearly document techniques for the benefit of others.
What is in this article?:
 A Few Tips for Solving the Unsolvable
 A few tips for Solving the Unsolvable
A few tips for Solving the Unsolvable
Ted Diehl
Engineering Consultant
Du Pont Co.
Wilmington, Del.
(www.dupont.com)
Half the engineering battle in solving problems is coming up with the right tools for the job at hand. For instance, what would you need to model a paper jam in a copier or the shock in a dropped cell phone? One approach uses two complementary analysis tools. One tool is engineering calculation software for applying mathematics, such as Mathcad from Mathsoft Inc., Cambridge, Mass. (www.mathsoft.com). The program lets engineers work naturally using applied math and engineering equations, mixed units, symbolic equations, numerical methods, graphing, text, and pictures, in an easily readable format. The second program is finiteelementanalysis software, such as Abaqus from Abaqus Inc., Pawtucket, R.I. (www.abaqus.com), for generalpurpose and nonlinear simulations. To show how these work together, let's examine a general approach to two different analyses. The first examines paper transport in a copier while the other involves the problem of drop shock in a cellular phone. The idea is not to present detailed solutions but rather a method for attacking problems that can be applied to other engineering problems. A paper modelAnyone using copiers has experienced paper jams and distorted images. Understanding these problems comes down to analyzing the soft rubber rollers that push and pull paper. The problem includes three categories of nonlinear behavior for structural mechanics: materials, kinematics (large deformations), and contact (changing boundary conditions). By observation, foamrubber rollers are repeatedly exposed to relatively large deformations and strains as they turn. Deformations are greater than can be described by Hooke's Law. These large, nonlinear, elastic deformations in the rollers call for a material law based on the theory of hyperelasticity. So the first hurdle to clear is in expanding your understanding of materials from simple Hooke's law, σ_{1} = E ε_{1}, to an OgdenHill hyperelastic strainenergy density function defined by: where W = strain energy density (psi); µ_{n}, α_{n}, and ß_{n} are material constants for an nterm representation; and J are related to strain (ε) through additional equations, and W is related to stress (σ) through various derivatives. This is the recommended material law and equation for nonlinear FEA codes when modeling foamrubber materials. Despite documentation in FEA manuals, it is difficult to use this complicated equation, especially when ones previous experience is limited to Hooke's law. But don't be put off by the expression's complexity. It can be tamed. Several experiments with the material provide enough data for the FEA program to provide a thirdorder curve fit. The FEA program then generates the constants for the equation. Because the engineering math program easily handles symbology, numerical algorithms, graphs, a range of units, and text, users can study such complex equations in more detail and see how different equation parameters affect material response (behavior). After deriving a workable set of material laws (for paper and rubber) we develop a viable approach to model the entire problem using FEA. While nonlinear closedform analysis is not likely to solve the entire systemlevel problem of paper moving in a copier, it does help gain ability and confidence. Studying closedform nonlinear analytical solutions for a nonlinear beam theory called the Elastica helps understand how to model the nonlinear kinematics of paper transported by the foam rubber wheels. (Elastica is the largedeformation, smallstrain beam theory described by Timoshenko and is covered in several texts. Bending a sheet of paper, for instance, produces a large deformation with small strains.) In this case, engineers could type the complex beam equations into Mathcad (just as they appear in the textbook) and ask for numerical solutions. These results can be compared to output from nonlinear Abaqus beam and shell models, and to physical experiments. Shock and fallThe other problem, understanding what happens to a cell phone when it falls to a hard floor, is somewhat more difficult. Typical cellphone housings are often crowded with chips, circuit boards, connectors, glass, foamrubber snubbers, and metal shields. The questions to answer include: What happens during impact? How much does the circuit board deflect? And how badly are critical components shocked? Answers come from a combination of experiments, closedform analysis, and nonlinear FEA simulation. Geometry, materials, contact, and loading are complex, so it is important to ensure models and information processed from them are accurate. A simplified cellphone impact test helps develop basicmodeling approaches along with datacollection and processing methods. To get a handle on the problem, a steel ball was dropped onto an instrumented phone cover, and the acceleration signal was compared to that from a simulated test. FEA models for this type of simulation use explicit dynamics technology, the same analysis method used to simulate crash tests. The first problem here: data not properly collected and manipulated. It is easily corrupted and distorted by aliasing. A simple example of aliasing is the apparent backward motion of wagon wheels in old Western movies. Frames in the film were taken at rates too slow to capture the wheel spoke's accurate locations. So at certain playback speeds, wheels appear to turn backwards, which obviously is not happening. A similar defect often occurs when collecting transient data in simulated drop tests if the data is collected from the FEA model at too slow a sampling rate. Portions of highfrequency noise in the FEA solutions may be recorded as lowfrequency data, and then misinterpreted. Data that is improperly collected cannot, in general, be fixed. If however, the data is collected properly and aliasing is avoided, then filtering can remove noise. In the physical experiment, raw data looks completely different from that collected during an FEA simulation of the same event. FEA results are extremely noisy. One way to remove the noise, assuming all data was collected without aliasing, is to use smoothing  a lowpass filter  that passes low frequencies and rejects high frequencies that make up noise. The problem with explicit dynamic FEA codes, such as Abaqus/Explicit, is that they use an adaptive timeincrement algorithm which gives data a varying time increment throughout the solution. This violates the requirement of digital signalprocessing (DSP) programs that include filtering and other Fourier analysis. Data must have a constant time increment. The math program provided a solution. Its spline function allowed reinterpolating the data to a constanttime increment. Then it was ready for lowpass DSP filtering. When handled properly, the two signals (from experiment and FEA) correlate quite well all across the scale. If the same type of processing is tried in a spreadsheet, several problems arise. First, spreadsheets can handle relatively small data sets and would not be able to open typical FEA data files that contain 100,000 lines of data and more. Second, the spreadsheet does not have sufficient capability to perform the manipulations required. Lastly, most of the calculations done in a spreadsheet would be indecipherable to those not intimately involved. Solving this shock problem called for learning about DSP, data sampling, aliasing, dealing with noisy data, details of Fourier analysis, and more. The details are beyond the scope of this text. But in a nutshell, the effort required several textbooks, a couple short DSP classes, and discussions with other engineers who commonly use DSP technology (acoustic engineers in this case). In particular, Mathcad's WYSIWYG format and general capability helped provide a solid understanding of DSP and a way around FEA's idiosyncrasies. When a project finishes, the details are stored in Mathcad (or engineering calculation and documentation software) for use by others. It's readily understandable, even to those unfamiliar with the software, and can be modified for new tasks by those with a license to it. Most importantly, the engineering math program lets engineers easily create a toolbox of functions that develops over the length of a career to allow processing of closedform equations, FEA, and experimental data in ways that neither approach could do alone. Where math software and FEA work wellAlthough the article describes two solutions that capitalize on the strengths of mathematics software and FEA, there are other complex nonlinear problems that also benefit from similar teamwork. For instance, when dealing with:
