The objective of this paper is to examine the utility of direct, numerical solution procedures, such as finite difference or finite element methods, for partial differential equations in chaotic dynamics. In the past, procedures for solving such equations to detect chaos have utilized Galerkin approximations which reduce the partial differential equations to a set of truncated, nonlinear ordinary differential equations. This paper will demonstrate that a finite difference solution is equivalent to a Galerkin solution, and that the finite difference method is more powerful in that it may be applied to problems for which the Galerkin approximations would be difficult, if not impossible to use. In particular, a nonlinear partial differential equation which models a slender, Euler-Bernoulli beam in compression issolvedto investigate chaotic motions under periodic transverse forcing. The equation, cast as a system of firstorder partial differential equations is directly solved by an explicit finite difference scheme. The numerical solutions are shown to be the same as the solutions of an ordinary differential equation approximating the first mode vibration of the buckled beam. Then rigid stops of finite length are incorporated into the model to demonstrate a problem in which the Galerkin procedure is not applicable. The finite difference method, however, can be used to study the stop problem with prescribed restrictions over a selected subdomain of the beam. Results obtained are briefly discussed. The end result is a more general solution technique applicable to problems in chaotic dynamics.