Forested environments feature a highly complex radiation regime, and solar radiation is hindered from penetrating into the forest by the 3D canopy structure; hence, canopy shortwave radiation varies spatiotemporally, seasonally, and meteorologically, making the radiant flux challenging to both measure and model. Here, we developed a synergetic method using airborne LiDAR data and computer graphics to model the forest canopy and calculate the radiant fluxes of three forest plots (conifer, broadleaf, and mixed). Directional incident solar beams were emitted according to the solar altitude and azimuth angles, and the forest canopy surface was decomposed into triangular elements. A ray tracing algorithm was utilized to simulate the propagation of reflected and transmitted beams within the forest canopy. Our method accurately modeled the solar radiant fluxes and demonstrated good agreement (R2≥0.82) with the plot-scale results of hemispherical photo-based HPEval software and pyranometer measurements. The maximum incident radiant flux appeared in the conifer plot at noon on June 15 due to the largest solar altitude angle (81.21°) and dense clustering of tree crowns; the conifer plot also received the maximum reflected radiant flux (10.91-324.65 kW) due to the higher reflectance of coniferous trees and the better absorption of reflected solar beams. However, the broadleaf plot received more transmitted radiant flux (37.7-226.71 kW) for the trees in the shaded area due to the larger transmittance of broadleaf species. Our method can directly simulate the detailed plot-scale distribution of canopy radiation and is valuable for researching light-dependent biophysiological processes.