Logo Море(!) аналитической информации!
IT-консалтинг Software Engineering Программирование СУБД Безопасность Internet Сети Операционные системы Hardware
Архив форумов ЦИТФорума
Море(!) вопросов - Море(!) ответов
 
 FAQFAQ   ПоискПоиск   ПользователиПользователи   ГруппыГруппы   РегистрацияРегистрация 
 ПрофильПрофиль   Войти и проверить личные сообщенияВойти и проверить личные сообщения   ВходВход 
Как правильно задавать вопросы

В чем проблема??

 
Перейти:  
Этот форум закрыт, вы не можете писать новые сообщения и редактировать старые.   Эта тема закрыта, вы не можете писать ответы и редактировать сообщения.    Список форумов Архив форумов ЦИТФорума -> Программирование
Предыдущая тема :: Следующая тема  
Автор Сообщение
Zinjvi



Зарегистрирован: 06.04.2012
Сообщения: 1

СообщениеДобавлено: Пт Апр 06 2012 23:14    Заголовок сообщения: В чем проблема?? Ответить с цитатой

Выдает ошибку. мол, память не может быть записана...
помогите пожалуйста



Код:
#include<iostream.h>
#include<math.h>
#include<stdlib.h>

//Глобальные  переменные

float Epsil_1,Epsil_2,PlotnEnerg,PlotnEnergU,Epsil_X,Epsil_Y,Gamma_XY,G;

//Объявление функций

float VnUsM(float ,float , float ,float );
float VnUsQ(float ,float , float ,float );

float KasNapr(float , float ,float ,float ,float ,float );
float NormNapr(float , float ,float ,float ,float ,float );
float GlNapr(float , float ,float ,float );
float Energ(float ,float ,float ,float );



//Главный код

float main()
{   
float L,H, B,dSt,dPolki;
float q1,q2,E, my,Rbt,Rb,x; // все в см. и кг.
float N,J;
L=900;
H=150;
B=200;
dSt=20;
dPolki=20;

N=2;
J=10;

q1=10;
q2=40;
E=325000;
my=0.2;
Rbt=12;
Rb=170;

float* VektVnUsM;
VektVnUsM = (float*)malloc(N*sizeof(float));

float* VektVnUsQ;
VektVnUsQ = (float*)malloc(N*sizeof(float));

float* MasPlotnEnerg;
MasPlotnEnerg = (float*)malloc(N*sizeof(float));

float* MasPlotnEnergU;
MasPlotnEnergU = (float*)malloc(N*sizeof(float));


float* V;
V = (float*)malloc(N*J*sizeof(float));


//Создадим массивы
int i,k;
float **MassivEnerg = new float*[N];
for (i=0;i<N;i++) MassivEnerg = new float*[J];

float **MassivEnergU = new float*[N];
for (i=0;i<N;i++) MassivEnergU = new float*[J];







//Пробный
int t,r;
int dim1 = 3, dim2 = 5;
int **pArr = new int*[dim1];
for (t=0;t<dim1;t++) pArr[t]=new int[dim2];
for (t=0;i<dim1;t++) for (r=0;r<dim2;r++) pArr[t][r]=i+r;
for (t=0;i<dim1;t++) for (r=0;r<dim2;r++) cout << pArr[t][r] <<";";










for (i=0;i<N;i++)
{
   x=L*(i+1)/N;
   VektVnUsM[i]=VnUsM(q1,q2,L,x);
   VektVnUsQ[i]=VnUsQ(q1,q2,L,x);
   //cout << x<< "="<<VektVnUsM[i];

   for (k=0;k<J;k++)
   {
      int y=0.5*H*(k+1)/J;

      //Посчитаем в сечении i и в волокне k величины: Epsil_1,Epsil_2,Epsil_X,Epsil_Y,Gamma_XY,G
      //и сохраним их в глобальных переменных
      GlNapr(my,NormNapr(VektVnUsM[i], dSt,dPolki,B,H,y),KasNapr( VektVnUsQ[i],  dSt, dPolki, B, H, y),E);
      
      //Посчитаем в сечении i и в волокне k величины: PlotnEnerg,PlotnEnergU (воспользовавшись из
      //глобадьных переменных величинами Epsil_1,Epsil_2,Epsil_X,Epsil_Y,Gamma_XY,G)
      //и сохраним их в глобальных переменных
      Energ( E, Rb, Rbt, my);
      
      //Присвоим Плотности єнергии деформации в сечении i и в волокне k.
      //Значения возьмем из глобальной пременной
   

      //   MasPlotnEnerg[i]=PlotnEnerg;


   //   MassivEnerg[i][k]=i+k;
   //   MassivEnergU[i][k]=PlotnEnergU;

   }

}

//for (i=0;i<N;i++) for (k=0;k<J;k++) cout << MassivEnergU [i][k] <<" ; " << '\n';


//for (i=0;i<N;i++)
//{   cout << '\n';
//   for (k=0;k<J;k++)
//   {
//      cout << MasPlotnEnerg[k] <<" ; " ;
//   }
//}



return 0;   
}

// Момент

float VnUsM(float q1,float q2, float L,float x)
{
   float gamma, ksi, M;
   gamma = q1/q2;
   ksi = x/L;
   M = q2*L*L*ksi*ksi*(3*gamma+ksi*(1-gamma))/6;
   return M;
}

// Поперечная сила

float VnUsQ(float q1,float q2, float L,float x)
{
   float gamma, ksi, Q;
   gamma = q1/q2;
   ksi = x/L;
   Q = 0.5*q2*L*ksi*(2*gamma+ksi*(1-gamma));
   return Q;
}




// Касательные напряжения

float KasNapr(float VnUsQ, float dSt,float dPolki,float B,float H,float y)
{   
   float I,S,alfa,betta,fi;


   if(H/2>y && (H-2*dPolki)/2<y)   
   {
      S=(B*H*H)/8-(B*y*y)/2;
   }
   else
   {
      S=(B*dPolki*H)/2-(B*dPolki*dPolki)/2+(dSt*H*H)/8-(dSt*H*dPolki)/2+(dSt*dPolki*dPolki)/2-(dSt*y*y)/2;
   }

   alfa=dSt/B;
   betta=dPolki/H;
   fi=alfa+betta*(1-alfa)*(8*betta*betta-12*betta+6);
   I=(fi*B*H*H*H)/12;
   
   return (VnUsQ*S)/(dSt*I);
}


//Нормальніе напряжения

float NormNapr(float VnUsM, float dSt,float dPolki,float B,float H,float y)
{
   float I,alfa,betta,fi,W;

   alfa=dSt/B;
   betta=dPolki/H;
   fi=alfa+betta*(1-alfa)*(8*betta*betta-12*betta+6);
   I=(fi*B*H*H*H)/12;
   
   W=I/y;

return VnUsM/W;
}

//Главные напряжения

float GlNapr(float my, float NormNapr,float KasNapr,float E)
{
   float p, q;
   // !!! Используются глобальные переменные Epsil_1,Epsil_2,Epsil_X,Epsil_Y,Gamma_XY,G


   Epsil_X=(1-my*my)*NormNapr/E;
   Epsil_Y=-(1+my)*my*NormNapr/E;

   G=E/(2*(1+my));

   Gamma_XY=KasNapr/G;

   p=-(Epsil_X+Epsil_Y);
   q = Epsil_X*Epsil_X-0.25*Gamma_XY*Gamma_XY;
   Epsil_1 = -(p/2)+sqrtf((p*p/4)-q);
   Epsil_2 = -(p/2)-sqrtf((p*p/4)-q);

   if (Epsil_1 > Epsil_2)
   {
   }
   else
   {
   Epsil_1=Epsil_2;
   Epsil_2=Epsil_1;
   }
   return 0;
}

//Текущая и предельная плотность потенциальной энергии деформации

float Energ(float E,float Rb,float Rbt,float my)
{
   float NadLode, Ebcu, Ebtu, Ebshu,Teta ;
   // !!! Используются глобальные переменные Epsil_1,Epsil_2,PlotnEnerg,PlotnEnergU,Epsil_X,Epsil_Y,Gamma_XY,G
   NadLode=(2*Epsil_2-Epsil_1)/Epsil_1;
   
   Ebcu=0.5*Rb*Rb/E;
   Ebtu=0.5*Rbt*Rbt/E;
   Ebshu=0.625*Rbt*Rb/E;

   PlotnEnergU = 0.5*NadLode*NadLode*((NadLode+1)*Ebcu-(NadLode-1)*Ebtu)-(1-NadLode*NadLode)*Ebshu;
   
   Teta=Epsil_X+Epsil_Y;
   
   PlotnEnerg=G*(Epsil_X*Epsil_X+Epsil_Y*Epsil_Y+(my*Teta*Teta)/(1-2*my)+0.5*Gamma_XY*Gamma_XY);   
   
   
   return 0; 
}
[/code]
Вернуться к началу
Посмотреть профиль Отправить личное сообщение
Показать сообщения:   
Этот форум закрыт, вы не можете писать новые сообщения и редактировать старые.   Эта тема закрыта, вы не можете писать ответы и редактировать сообщения.    Список форумов Архив форумов ЦИТФорума -> Программирование Часовой пояс: GMT + 4
Страница 1 из 1

 
Перейти:  
Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете голосовать в опросах


Powered by phpBB © 2001, 2002 phpBB Group
Русская поддержка phpBB

 

IT-консалтинг Software Engineering Программирование СУБД Безопасность Internet Сети Операционные системы Hardware

Информация для рекламодателей PR-акции, размещение рекламы — adv@citforum.ru,
тел. +7 495 6608306, ICQ 232284597
Пресс-релизы — pr@citforum.ru
Послать комментарий
Информация для авторов
This Web server launched on February 24, 1997
Copyright © 1997-2000 CIT, © 2001-2006 CIT Forum
Внимание! Любой из материалов, опубликованных на этом сервере, не может быть воспроизведен в какой бы то ни было форме и какими бы то ни было средствами без письменного разрешения владельцев авторских прав. Подробнее...