Форум программистов, компьютерный форум, киберфорум
Наши страницы

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
buttermas
0 / 0 / 0
Регистрация: 20.12.2013
Сообщений: 1
#1

метод конечных элементов на си++ - C++

23.12.2013, 20:14. Просмотров 1418. Ответов 2
Метки нет (Все метки)

Всем привет!нужно написатьт прогу для решения одномерной задачи методом конечных элементов. в чем ошибка?
int _tmain(){
setlocale ( LC_ALL, "rus" );
int i=0,j=0;
double h=0.01,N=100,a[100][100],u[100],fi[100],alf[100],bet[100],xi=-1;
a[0][0]=200.0066667;
a[0][1]=-99.99833333;
fi[0]=(0.009999177558*sin(3.141592654*h*(i+1)));
alf[0]=-a[0][1]/a[0][0];
bet[0]=fi[0]/a[0][0];
for(i=1;i<N-1;++i){
a[i][i]=200.0066667;
a[i][i+1]=-99.99833333;
a[i][i-1]=-99.99833333;
fi[i]=0.;
for(j=1;j<200;++j){
xi+=h;
fi[i]+=(sin(3.14*(h*xi+h*i))*(1-abs(xi)));
}
fi[i]*=h*h;
xi=-1;
if(i==98){
a[i+1][i+1]=100.0033333;
a[i+1][i]=-99.99833333;
fi[99]=0.;
for(j=1;j<100;++j){
xi+=h;
fi[99]+=(sin(3.14*(h*xi+h*i))*(1-abs(xi)));
}
fi[99]*=h*h;
}
if(i==1){
alf[i]=-a[i-1][i]/a[i-1][i-1];
bet[i]=fi[i-1]/a[i-1][i-1];
}
alf[i+1]=-a[i][i+1]/(a[i][i-1]*alf[i]+a[i][i]);
bet[i+1]=(fi[i]-a[i][i+1]*bet[i])/(a[i][i-1]*alf[i]+a[i][i]);
if(i==98){
alf[99]=-a[99][98]/a[99][99];
bet[99]=fi[99]/a[99][99];
u[99]=(fi[99])/(a[99][99]);
i=100;
u[100]=0;
}
}

cout<<" u("<<100*h<<")="<<u[100]<<'\n';
cout<<" u("<<99*h<<")="<<u[99]<<'\n';
for(i=98;i>0;--i){
u[i]=alf[i]*u[i+1]+bet[i];
cout<<" u("<<h*i<<")="<<u[i]<<'\n';
if(i==1){
u[i-1]=0;
}
}
cout<<" u("<<0*h<<")="<<u[0]<<'\n';
getch();
return 0;
}
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
23.12.2013, 20:14
Здравствуйте! Я подобрал для вас темы с ответами на вопрос метод конечных элементов на си++ (C++):

Метод конечных элементов - C++
Помогите с алгоритмами по программированию МКЭ. С математическими основами метода я разобрался, но не могу понять с чего начать...

Метод конечных элементов - ищу баг - C++
Всем привет, конкретно застрял на вычислении матрицы жесткости (пока для одного элемента) из этой статьи...

Метод конечных разностей для краевой задачи - C++
Помогите пожалуйста. Надо написать программу, которая решает краевую задачу для ОДУ 2 порядка методом конечных разностей. Никак не могу...

Метод конечных автоматов для списка слов - C++
Дан список слов X1, X2, …, Xk и строка Y длины n. Определить, входит ли одно из слов Xi в строку Y как подстрока за время C∙n (где C –...

Метод Конечных Разностей для уравнения в частных производных - C++
Здравствуйте, требуется ваша помощь в нахождении ошибки в программе. Вычисляю по явной схеме. Выдает какие-то очень большие числа. ...

Метод медиан из трех элементов VS улучшенный быстрый метод сортировки(метод Бентли-Макилроя) - C++
Здравствуйте! Дали весьма интересное задание. Сравнить два вышеуказанных метода сортировки для массива из 10000 элементов, результаты...

2
lpvoid
23.12.2013, 20:17
  #2

Не по теме:

make me unsee it.

0
Hoottie_McGOOB
107 / 105 / 5
Регистрация: 04.10.2013
Сообщений: 231
24.12.2013, 14:56 #3
C++
1
cout<<" u("<<100*h<<")="<<u[100]<<'\n';
как минимум в этой строке выход за пределы массива... размерность 100... элементов 100... а обращение к 101-му элементу.
0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
24.12.2013, 14:56
Привет! Вот еще темы с ответами:

Модуль формирования глобальной С-матрицы для метода конечных элементов - C++
Доброго времени суток! Возможно у кого нибудь имеются уже готовые программы на С++ связанных с методом конечных элементов, очень поможете...

СЛАУ. Метод обратной матрицы, метод Гаусса, метод Крамера, метод Зейделя - C++
Помогите ребят. Не могу построить алгоритмы для этих методов Язык C++

Метод конечных элементов - Численные методы
Здравсвтуйте! Помогите пожалуйста разобраться с методом конечных элементов в двумерном случае с треугольным конечным элементом... ...

Метод Конечных Элементов - Численные методы
Помогите, с понятием метода. Пишу программу расчета деталей на прочность. Пока сделал деление детали на конечные элементы. Делю на...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
3
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.