The BEM combined with the time-stepping scheme is employed in the nonlinear numerical simulation of the development process of radiation wave profile and hydrodynamic force caused by large amplitude oscillating cylinder. To ensure stability and accuracy of the simulation, special consideration is paid on a series of difficulties in numerical modeling, such as the compatibility of free surface and body surface conditions at the intersection point and the open boundary condition problem. Numerical calculation is carried out for some study cases of different shaped cylinders undergoing large amplitude oscillating motions with distinct modes, the results show. some nonlinear phenomena and indicate the feasibility of the numerical approach. The BEM combined with the time-stepping scheme is employed in the nonlinear numerical simulation of the development process of radiation wave profile and hydrodynamic force caused by large amplitude oscillating cylinder. To ensure stability and accuracy of the simulation, special consideration is paid on a series of difficulties in numerical modeling, such as the compatibility of free surface and body surface conditions at the intersection point and the open boundary condition problem. Numerical calculation is carried out for some study cases of different shaped cylinders undergoing large amplitude oscillating motions with distinct modes, the results show. some nonlinear phenomena and indicate the feasibility of the numerical approach.