[DEM]

Родилась идея, о создании бесплатного аналога программ Арабат, Кристал, Вест и отдельных модулей НормКада....
Данный проект можно сделать с открытыми исходниками а разработку его вести на основе ГитХаба или Google Code....
В качестве языка программирования предлагаю принять VB.net т.к. он все таки победил в теме....
Можно было конечно принять матерный, но к сожалению, тогда придется проект делать закрытым, и исходники хранить на флешках в защищенных от доступа детей местах...
Прошу высказывать свои пожелания и предложения к данному проекту...

PS. Прошу всех кто согласен участвовать в проекте отписаться в личку....

Список участников:
DEM
bargool
RrRR
ябс
Cosinus
swell{d}
Akim_1989
Djur
Ratmir
Yasnogor
P1@t0n
GradRusStroy
gabidullin

PS. Всем заинтересованным, во включении каких то программ, прошу заполнить форму, чтобы можно было добавлять программу....

[Meknotek]

Сообщение от Бахил Meknotek, ты зачем дровишек подкидываешь? Сейчас ещё на 4-х страницах будут пинать твой питон.

А чем еще заниматься? Пока люди делом занимаются (zamtmn, AY), хоть змею попинаем. Один пашет, семеро питоном машут...

[Akim_1989]

В питоне можно найти библиотеки на любой вкус и цвет. Хочешь графику, хочешь алгебру, хочешь мкэ. Свой велосипед редко приходится создавать. Простой язык, мощный функционал - что ещё нужно инженеру. Мои потребности в автоматизации повседневной рутины удовлетворяет полностью.

[Meknotek]

Сообщение от Akim_1989 В питоне можно найти библиотеки на любой вкус и цвет. Хочешь графику, хочешь алгебру, хочешь мкэ. Свой велосипед редко приходится создавать. Простой язык, мощный функционал - что ещё нужно инженеру. Мои потребности в автоматизации повседневной рутины удовлетворяет полностью.

А свои библиотеки есть? У меня в процессе (завершены на 1%) пара-тройка своих. Для еще большего удобства пользования

[Akim_1989]

Библиотек нет. Программы есть. Есть для анализа РСУ в стержнях, в фермах, есть наработки по фундаментам(классы грунтов, скважин, фундаментов, определение осадки фундамента, с просадкой ещё не разобрался, есть вопросы) И ещё куча всяких задумок, которые ещё не воплотились в жизнь...

[Chardash]

Добрый вечер.
Сделал программу по прототипу из методички выше выложенной ранее, все работает, но не совсем корректно выводит таблицу, может есть еще ошибки. Еще раз спасибо тем, кто помог с вопросами. Буду рад, если кого то заинтересуют языки С++ и/или Fortran.
зы Если кто-то заметит еще ошибки, просьба указать.

Код:

[Выделить все]

 
/*
Программа COLUMN, при выполнении которой определяются параметры армирования ряда равнопрочных базовому варианту серии 1.020-1/83
новых вариантов армирования
Исходные данные см. список литературы, 1.
Использованная литература:
1. "Программирование инженерных задач на Fortran" Рак Н.А., Смех В.И., Щербак С.Б.
2. "Сравнительный курс языков программирования Fortran и С" Григорьев С.А., Ишанов С.А.
*/

#include <iostream>
#include <fstream>
#include <cmath>
#include <iomanip>
#include <conio.h>

using namespace std;

// Переменные
int  nbd;				//Количество вариантов прочности бетона
int  i;					//Переменная цикла по размеру поперечного сечения
int  j;					//Переменная цикла по прочности бетона
int  m;					//Координата точки интерполяции по гибкости
int  n;					//Координата точки интерполяции по соотношению
//длительной части нагрузки к полной
float  dmin;			//Минимальный диаметр арматуры, мм
float  Asmin;			 //Площадь арматуры при минимальном диаметре, мм2

//Параметры базового варианта колонны

float  nln;				//Отношение длительной нагрузки к полной
float  lr;				//Расчетная длина колонны, мм
float  bb;				//Размер поперечного сечения, мм
float  Rbb;				//Призменная прочность бетона, МПа
float  Asb;				//Площадь арматуры, мм2
float  Rsb;				//Расчетное сопротивление арматуры, МПа
float  lamb;			//Гибкость колонны
float  msb;				//Коэффициент продольного армирования
float  ncb;				//Несущая способность, кН

//Параметры проектных вариантов колонны

float  bd[2] = { 300, 400 };			//Массив размеров поперечного сечения, мм
float  Rbd[8];			//Массив призменных прочностей бетона, МПа
float  Rsd;				//Расчетное сопротивление арматуры, МПа
float  lamd;			//Гибкость колонны
float  Asd;				//Площадь арматуры, мм2
float  ncd;				//Несущая способность, кН
float  ncd1;			//Окончательная несущая способность, кН
float  msd;				//Коэффициент продольного армирования
float  pn;				//Параметр интерполяции по вертикали
float  pl;				//Параметр интерполяции по горизонтали
float  Fi;				//Функция для нахождения коэффициента Fi
float  Fib;				//Величина коэффициента Fib
float  Fisb;			//Величина коэффициента Fisb
float  p, t, r;			//Промежуточные коэффициенты
float  k1, k2;			//Промежуточные коэффициенты

float lamt[8] = { 6, 8, 10, 12, 14, 16, 18, 20 };
float nlnt[3] = { 0, 0.5, 1.0 };

float Fibt[3][8] = { .93, .92, .91, .90, .89, .88, .86, .84,
.92, .91, .90, .89, .86, .82, .78, .72,
.92, .91, .89, .86, .82, .76, .69, .61 };

float Fisbt[3][8] = { .93, .92, .91, .90, .89, .88, .86, .84,
.92, .92, .91, .89, .88, .86, .83, .79,
.92, .91, .90, .89, .87, .84, .79, .74 };


/* Подпрограмма для определения параметров интерполяции
по табл. 26 и 27 */

float point(int &m, int &n, float &pn, float &pl)
{
	for (m = 0; m < 3; m++)
	{
		if (nln <= nlnt[m]) break;
	}
	if (m == 0) { m = 1; pn = 0; }
	else
		pn = (nln - nlnt[m - 1]) / (nlnt[m] - nlnt[m - 1]);

	for (n = 0; n < 8; n++)
	{
		if (lamb <= lamt[n]) break;
	}
	if (n == 0) { n = 1, pl = 0; }
	else
		pl = (lamb - lamt[n - 1]) / (lamt[n] - lamt[n - 1]);

	return 0;
}

/* Функция для определения коэффициента продольного изгиба
Fib по табл. 26 */

float t26(int &m, int &n, float &pn, float &pl)
{
	Fib = (1 - pn) * (1 - pl) * Fibt[m - 1][n - 1] + (1 - pn) * pl * Fibt[m - 1][n] +
		pn * (1 - pl) * Fibt[m][n - 1] + pn * pl * Fibt[m][n];

	return 0;
}

/* Функция для определения коэффициента продольного изгиба
Fisb по табл. 27 */

float t27(int &m, int &n, float &pn, float &pl)
{
	Fisb = (1 - pn) * (1 - pl) * Fisbt[m - 1][n - 1] + (1 - pn) * pl * Fisbt[m - 1][n] +
		pn * (1 - pl) * Fisbt[m][n - 1] + pn * pl * Fisbt[m][n];

	return 0;
}

/* Функция для определения коэффициента продольного изгиба
Fi по табл. 26 и 27а */

float Fit(float Rsb, float Rbb, float msb, float Fib, float Fisb)
{
	float alfs;

	alfs = Rsb * msb / Rbb;
	Fi = Fib + 2 * (Fisb - Fib) * alfs;
	if (Fi > Fisb) { Fi = Fisb ;}

	return 0;
}

int main()
{
	setlocale(LC_ALL, "Russian");

	//Ввод исходных данных из внешнего файла

										
	// Создание потоков двух классов
	ifstream fin;
	ofstream fout;
	// Открытие файлов в соответствующих режимах и связывание их с потоками
	fin.open("column.dat");
	if (!fin.is_open())
	{
		cout << "Создайте файл с данными column.dat (см. методичку)";
		_getch();
		return -1;
	}

	fin >> nln >> lr >> bb >> Rbb >> Asb >> Rsb; 
	fin >> dmin;
	fin >> nbd;
	for (int i(0); i < nbd; i++)
	{
		fin >> Rbd[i];
	}
	fin >> Rsd;

	// Открытие файла результатов вычислений

	fout.open("column.res");
	fout << "Параметры базового расчета" << endl;
	lamb = lr / bb;

	// Проверка возможности использования метода расчета

	if (lamb > 20)
	{
		cout << "Гибкость базовой колонны больше 20" << endl;
		_getch();
		return -1;
	}
	Asmin = 4. * atan(1.0) * pow(dmin, 2);
	msb = Asb / (pow(bb, 2) - Asb);

	//Определение коэффициента Fi для базового варианта

	point(m, n, pn, pl);
	t26(m, n, pn, pl);
	t27(m, n, pn, pl);
	Fit(Rsb, Rbb, msb, Fib, Fisb);

	// Определение несущей способности базового варианта

	ncb = Fi*((pow(bb, 2) - Asb)*Rbb + Rsb*Asb);

	// Вывод результатов в файл результатов для базового варианта

	fout.fill('-');
	fout << "+" << setw(25) << "+" << setw(25) << "+" << endl;
	fout.fill(' ');
	fout << "|" << setw(11) << "Lr = " << lr << " мм " << setw(6) << "|" << setw(12) << "Fib = " << Fib << setw(9) << "|" << endl;
	fout << "+------------------------+------------------------+\n";
	fout << "|" << setw(11) << "bb = " << bb << " мм; " << " lamb = " << lamb << "; " << " Rsb = " << Rsb << " МПа " << setw(5)
		<< "|" << endl;
	fout << "+------------------------+------------------------+\n";
	fout << "|    Rbb, МПа  |    Fi   |   Asb, mm2 |  Ncb, kN  |\n";
	fout << "+------------------------+------------------------+\n";
	fout << "|" << setw(9) << Rbb << setw(6) << "|" << setw(7) << Fi << setw(3) << "|" << setw(8) << Asb << setw(5) << "|"
		<< setw(10) << ncb * 1.e-3 << setw(2) << "|" << endl;
	fout << "+------------------------+------------------------+\n";


	// Расчет проектных вариантов(ПВ) колонн начало цикла по поперечному размеру

	fout << "Параметры проектных вариантов" << endl;

	for (i = 0; i < 2; i++)  // Цикл по поперечному размеру ПВ
	{
		lamd = lr / bd[i];

		// Проверка возможности использования метода расчета

		if (lamb > 20)
		{
			cout << " Bd = " << bd[i] << "Гибкость базовой колонны больше 20" << endl;
			_getch();
			return -1;
		}

		// Определение коэффициентов Fib Fisb для базового варианта
		// Начало цикла по прочности бетона

		for (int j(0); j < nbd; j++)
		{
			msd = Asmin / (pow(bd[i], 2) - Asmin);

			// Определение коэффициента Fi для минимального армирования

			Fi = Fit(Rsd, Rbd[j], msd, Fib, Fisb);

			// Несущая способность при минимальном армировании

			ncd = Fi*((pow(bd[i], 2) - Asmin) * Rbd[j] + Rsd * Asmin);
			Asd = Asmin;

			if (ncd < ncb) {                                   // расчетное армирование
				p = Rsd / Rbd[j];
				Asd = (ncb / Fisb - pow(bd[i], 2) * Rbd[j]) / (Rsd - Rbd[j]);
				msd = Asd / (pow(bd[i], 2) * Rbd[j]) / (Rsd - Rbd[j]);

				if (msd > 0.5 / p || lamd == 6) goto L10;            // случай Fi = Fisb
				t = ncb / (Fib * Rbd[j] * pow(bd[i], 2));
				r = 2 * (Fisb / Fib - 1);
				k1 = 0.5 * ((p * r - 1) + (p - 1) + t) / ((p * r - 1) * (p - 1));
				k2 = (1 - t) / ((p * r - 1) * (p - 1));
				if (k1 > 0.0)
					Asd = pow(bd[i], 2) * (-k1 + sqrt(pow(k1, 2) - k2));
				else
					Asd = pow(bd[i], 2) * (-k1 - sqrt(pow(k1, 2) - k2));
				msd = Asd / (pow(bd[i], 2) - Asd);
			}
		L10:
			if (j == 1) {
				fout.fill('-');
				fout << "+" << setw(25) << "+" << setw(25) << "+" << endl;
				fout.fill(' ');
					fout << "|" << setw(11) << "BD = " << bd[i] << " мм " << setw(7) << "|" << setw(12) << "Fib = " << Fib << setw(9) << "|" << endl;
					fout << "|" << setw(14) << "Lamd = " << lamd << setw(10) << "|" << setw(12) << "Fisb = " << Fisb << setw(9) << "|" << endl;
					fout << "+--------------+---------+------------+-----------+\n";
					fout << "|    Rbd, МПа  |    Fi   |   Asd, mm2 |  Ncd1, kN |\n";
					fout << "+--------------+---------+------------+-----------+\n";
				}

				Fit(Rsd, Rbd[j], msd, Fib, Fisb);
				ncd1 = Fi * ((pow(bd[i], 2) - Asd) * Rbd[j] + Rsd * Asd);
				fout << "|" << setw(9) << Rbd[j] << setw(6) << "|" << setw(7) << Fi << setw(3) << "|" << setw(10) << Asd << setw(3) << "|"
					<< setw(9) << ncd1 * 1.e-3 << setw(3) << "|" << endl;
				}
			fout << "+------------------------+------------------------+\n";

		} // Конец цикла по прочности бетона НВ
	fin.close();
	fout.close();
	return 0;
}

[gomer]

привет из 80х =)

[Chardash]

gomer, отчет легко поправить, вот я сделал вывод в html. По моему, нормально, современно

Изменить нужно несколько строк в коде.

Код:

[Выделить все]

 
fout << "<TABLE CELLSPACING=\"0\" COLS=\"4\" BORDER=\"1\"><COLGROUP SPAN=\"4\" WIDTH=\"113\"></COLGROUP>" << endl;
fout << "<TR><TD STYLE=\" \" COLSPAN=2 HEIGHT=\"26\" ALIGN=\"CENTER\" VALIGN=MIDDLE BGCOLOR=\"#83CAFF\"><FONT SIZE=3>Lr =  " << lr << " мм</FONT></TD>" << endl;
....

Единственное, согласен, без кнопки неудобно. Но и это поправимо.

[Chardash]

В коде где-то ошибка, несколько минусовых площадей в конце отчета, изменил в 172 строке на PI все равно, площадь отрицательная, убрал переходы goto. Все равно пока минус. В исходном коде фортран как и в онлайн калькуляторе все нормально, в программе на С++ пока что не нашел причины. Добавил онлайн калькулятор COLUMN там же, по ссылке выше сообщением. Красивый интерфейс уже делать не буду, для практических задач данный калькулятор возможно не всем подойдет. Но, если нужно, сделаю.
Картинка Copyrights http://xkcd.ru & http://xkcd.com.

[zamtmn]

Chardash
>>По моему, нормально, современно
Не, несовременно. Если для изменения "отчета" нужно ковырять программу - это никак не нормально.

[Chardash]

Offtop: zamtmn, согласен, мне уже тоже не нравится. Оставлю как есть для получения результата в читаемом виде.

[Chardash]

Еще немного доработал калькулятор этой прекрасной программы.
Почти закончил дескоптное приложение (на убунте), осталось добавить генерацию отчета в пдф и сам отчет изменить. Хорошо, если успею.
Еще раз всем спасибо за внимание.

ps справку добавлю еще.

[DEM]

Собственно я сейчас занимаюсь связкой GMSH и Скад.
Мне сейчас API SCADa на NET платформу переводят...

[Chardash]

ясно. Мне больше интересен Code Aster - Salome. В OpenFoam тоже много кода, но больше по механике сплошных сред и методу конечных объемов, очень хотелось бы поработать и с ними. Из плюсов то, что много кода, самому писать матрицы долго бы пришлось, а тут уже есть, останется разобраться . Разбираться нужно... Будет время и возможность как говорится

[DEM]

Собственно говоря программа предполагает следующее...
1. Помощь в подготовке файла geo.
1.1 Удаление дублей вершин
1.2 Удаление дублей отрезков
1.3 Включение отрезков (Line{3, 7, 11} In Surface{76}
1.4 Разбивка отрезков Суфрайсов если на них лежит отрезок другого Суфрайса т.к. (Line{3, 7, 11} In Surface{76} в данном случае не работает
2. Полуавтоматическое создание групп до и после разбития сетки.
3. Экспорт в СКАД расчетной модели, на этом этапе экспорт можно сделать и в Code Aster - Salome и в Робот и в Лиру.

Плюс к этой же программе, можно сделать дополнение по распечатке результатов СКАДа, есть уже наработки по полуавтоматической распечатке результатов, что раз 20 быстрее чем в СКАДе 11, на и раз в 5 быстрее чем в 21-й версии. Правда вот самая большая сложность будет со стержнями, там придется по думать....

[DEM]

В принципе сейчас возникла ситуация, есть набор инструментов по получению данных из расчетных программ....
Из СКАДа и Робота можно, получать данные.
А дальше вот вопрос, то ли обрабатывать так же в Vb.net то ли считать в какой то другой программе и формировать там же отчёты.
пока если честно склоняюсь к расчету в том же Vb.net....
Хотя вот подходящих инструментов для создания отчетов в Vb пока не искал.
Возможно приму самый простой вариант, создание html...

[gesper]

Сообщение от DEM Возможно приму самый простой вариант, создание html...

(И кто то крикнул из ветвей.....)
Сделай два варианта. HTML, с видом на свое усмотрение, и xml файл таблицей с данными, которые потом можно любой другой софтиной прочитать. Разницы между ними не много будет.

[DEM]

Ну да вариант нормальный.
Я уже делал интересные отчеты в Scilab.
Можно в нем сделать отдельный модуль для обработки таких отчётов.

[veb86]

Всем привет, Только узнал про данную тему и решил тоже поучаствовать. Читать 70 страниц некогда, расскажите на чем по решали то? разработкой занимаются? какие успехи? на чем пишем?

[DEM]

Сообщение от veb86 Всем привет, Только узнал про данную тему и решил тоже поучаствовать. Читать 70 страниц некогда, расскажите на чем по решали то? разработкой занимаются? какие успехи? на чем пишем?

Не а не занимаемся, ща пока другими делами занимаюсь.
Но задел у меня например есть по Scilab.

[veb86]

Жаль что работы не стартанули, я вот пытаюсь копать в область MultiCAD написание инструментов под бесплатный НаноКАД, плюсы на лицо - это где бы я в будущем не работал я его смогу скачать и работать, а не клянчить денег на СОФТ. Плюс .Net обладает множеством возможностей! Это мои мысли для написания общего бесплатного софта