با استفاده از کُد corsika نسخه35/7 و مدل اندرکنش هادرونی انرژی بالای qgsjet-04 تعداد 6000 بهمن گسترده هوایی برای اولیه های آهن و پروتون ومنیزیم و اکسیژن در محدوده انرژی ev1017-1014 و زوایای تابش صفر درجه وارتفاع 110متر(ارتفاع آزمایشkascade) شبیه سازی شده است، با استفاده از نتایج شبیه سازی ویژگی های مهمی از جمله(طیف انرژی ، نسبت باردار بودن به خنثی بودن ، تعداد هادرون ها، بیشترین انرژی هادرون ها، توزیع عرضی چگالی وتوزیع عرضی چگالی انرژی و..) از هادرون های ثانویه به همراه توزیع عرضی چگالی و توزیع عرضی چگالی انرژی برای الکترون ها ومیون ها در اولیه های پروتون ،اکسیژن،آهن ومنیزیم در انرژی های اولیه 1014 تا1017 الکترون ولت مطالعه شده است. مانشان داده ایم که طیف انرژی هادرون ها با افزایش انرژی ثانویه ها کم می شودوبرای اولیه های با جرم بیشتر تعداد کمتری از ثانویه ها را خواهیم داشت وهمچنین نشان دادیم که این طیف از قانون توانیe-?تبعیت می کندو ?(شاخص طیفی) رابه عنوان تابعی از انرژی بدست آوردیم، سپس نتیجه گرفتیم که توزیع عرضی چگالی وچگالی انرژی ثانویه ها با افزایش فاصله از مرکز بهمن کاهش می یابد،واین کاهش برای ثانویه های مختلف(الکترون،میون،هادرون) باجرم هاوانرژی های اولیه مختلف، متفاوت است .