Iron-Gallium alloys (Galfenol) are promising transducer materials that combine high magnetostriction, desirable mechanical properties, high permeability, and a wide operational temperature range. Most of all, the material is capable of operating under tensile stress, and is relatively resistant to shock. These materials are generally characterized using a solid, cylindrically-shaped specimen under controlled compressive stress and magnetization conditions. Because the magnetostriction strongly depends on both the applied stress and magnetization, the characterization of the material is usually conducted under controlled conditions so each parameter is varied independently of the other. However, in a real application the applied stress and magnetization will not be maintained constant during operation. Even though the controlled characterization measurement gives insight into standard material properties, usage of this data in an application, while possible, is not straight forward. This study presents an engineering modeling methodology for magnetostrictive materials based on a piezo-electric governing equation. This model suggests phenomenological, nonlinear, three-dimensional functions for strain and magnetic flux density responses as functions of applied stress and magnetic field. Load line performances as a function of maximum magnetic field input were simulated based on the model. To verify the modeling performance, a polycrystalline magnetostrictive rod (Fe-Ga alloy, Galfenol) was characterized under compressive loads using a dead-weight test setup, with strain gages on the rod and a magnetic field driving coil around the sample. The magnetic flux density through the Galfenol rod was measured with a sensing coil; the compressive loads were measured using a load cell on the bottom of the Galfenol rod. The experimental results are compared with the simulation results using the suggested model, showing good agreement.