[Akim_1989]

Здравствуйте, коллеги. Пишу программу для интерполяции геологии между тремя скважинами. Есть 3 скважины. Для каждой указана последовательность слоёв.(у слоёв есть данные о толщине, отметки верха и низа).
Первое, что пришло в голову-это генерировать для каждого одинакового слоя по трём точкам уравнение плоскости, которая совпадает с подошвой слоя. В итоге получаем набор плоскостей, из которого можно найти отметку подошвы каждого слоя и построить скважину внутри этого треугольника.
Этот метод очень хорошо работает, когда количество слоёв в скважинах и их последовательность совпадает. А что делать в сложных случаях:

- количество слоёв разное
- последовательность разная
- есть повторяющиеся слои.

Помогите с алгоритмом построения в данном случае. Если у кого-нибудь есть опыт в данной теме, то буду очень благодарен за помощь.

[trir]

1. мощность слоя = 0
2., 3.
a : a : a
b : b : 0
0 : a1 : a1
c : c : c
d : d :d

[DEM]

Я разбирался с данной проблемой...
Надо сразу ставить задачу корректно, иначе заплутаете...
Надо делать сперва триангуляцию по сетке скважин и по сетке точек, где вы собираетесь определять характеристики грунтов....

[Akim_1989]

DEM,С триангуляцией уже разобрался. Если есть триангуляция по скважинам, то зачем ещё по точкам делать?
trir, то есть нужно выбрать скважину с бОльшим кол-вом слоёв, а потом остальные дополнить нулевыми слоями?

[trir]

точнее надо строить литологическую колонку на всю плошадку (все скважины)

[DEM]

Ну понимаешь, нам же нужен результат под какой то точкой а не скважиной.
Я предполагал определять именно несущую способность грунта, а потом интерполировать ее между точками.
К сожалению сейчас занят и не могу заняться этой задачей....

[Владимир.]

Попробуйте нарисовать это самостоятельно на бумаге ручками. А потом пробуйте автоматизировать

[Brandashmыg]

Строю в архикаде объёмными элементами каждый ИГЭ и после делаю разрез где надо.

[Akim_1989]

Нашёл статью http://blog.simmakers.ru/methodika-p...heskoi-modeli/. Но пока не допёр, как реализовать.

[DEM]

Алгоритм примерно думаю можно продумать...
Но придется не хило по кумекать....
Может быть стоит немного другим путем все таки пойти, и определять именно несущую способность...
Тут наверное все таки лучше строить уравнения поверхностей, а не тупо делать разрезы.
Я думаю что наверное проще сделать некую триангуляцию сперва а потом строить уравнения плоскостей.
И уже из плоскостей собирать объемные элементы...

[Владимир.]

Сообщение от DEM Я думаю что наверное проще сделать некую триангуляцию сперва а потом строить уравнения плоскостей.

В этом и проблема у ТС. Как соединить поверхностью подошву ИГЭ, если он то есть, то его нет, а то вообще, два на скважину?

[swell{d}]

всё просто. достаточно набор слоёв во всех скважинах сделать одинаковым, а менять только толщины этих слоёв. причём эта задача должна лечь на пользователя, задающего исходные данные. программа сама не должна принимать таких решений

[Akim_1989]

Сообщение от swell{d} всё просто. достаточно набор слоёв во всех скважинах сделать одинаковым, а менять только толщины этих слоёв. причём эта задача должна лечь на пользователя, задающего исходные данные. программа сама не должна принимать таких решений

Это самый простой путь в плане алгоритма. Но в плане задания исходных данных сложный) Это не наш метод

[swell{d}]

Offtop: двигайтесь от простого к сложному

[trir]

Нажмите для просмотра

http://www.youtube.com/watch?v=fIEe2PxTSJA

[nikе]

В моей программе VarC1 я к примеру использовал ABOS интерполяцию, и ввод исходных данных осуществлял в файле xls.

[DEM]

Сообщение от nikе В моей программе VarC1 я к примеру использовал ABOS интерполяцию, и ввод исходных данных осуществлял в файле xls.

Хммм, вид симпатичный построенной модели, но видны ошибки триангуляции.
Все таки как мне кажется триангуляция, должна происходить на более равные элементы.

[nikе]

Те огрехи, которые видны вызваны ошибкой программы VarC1 которую я устранил. Объемники я строил только для контроля самой интерполяции геологии, в Лире объемная модель грунта не наглядна, поэтому еще пришлось простенький 3D вьювер написать, он постоянно крутит модель вокруг центральной вертикальной оси, более интерактивный вьювер мне писать не интересно. Все же основная функция этой программы именно генерация переменного и интерполированного коэф. постели, а вспомогательная к примеру объемники сгенерировать, или вычислить объем слоя грунта, или просадочные воронки промоделировать. Даже зимой на ней пытался смоделировать ленточный свайный фундамент как сплошной плитный ростверк...

----- добавлено через ~2 мин. -----
А и самое главное в этой программе ввод исходных данных в табличном процессоре, а то в самопальных таблицах типа которые в Лира ГРунт и других программах можно ногу сломать.

[Akim_1989]

Крутая программа) Думаю пока буду фиктивные слои вводить вручную. Потом может додумаюсь как улучшить)))

[nikе]

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

[Akim_1989]

Написал программу. Сравнил с ФОКом и Лирой ГРУНТ. С ФОКом очень близко результаты (видимо они тоже плоскости используют). С ГРУНТом разница в толщинах слоёв до 1м доходит. Видимо у них граница слоя - это поверхность, но странно, что отметки естественного рельефа точно совпадают (Может косякнул при вводе данных)

[DEM]

nikе
В каком виде у тебя получается уравнение поверхности????
Какую библиотеку использовал???

[nikе]

->DEM
http://m.dressler.sweb.cz/AOSIM.pdf
http://m.dressler.sweb.cz/ABOS.htm
там подробно расписана ABOS интерполяция, дается сравнительный анализ эффективности, применимости и точности других методов ещё частичный перевод на русиш где-то есть. Думал сначала библиотеку GsTL применить, но ума не хватило разобраться что к чему. Уравнения поверхности никакого нет, ABOS - метод последовательного натяжения и сглаживания. Исходные данные - массив точек и прямоугольная область, которая дробится сеткой по оси Х и У и в узлах сетки вычисляются высоты поверхности.

[DEM]

nikе
Сразу анекдот вспомнился....

Сидят двое чукчей(Ч) на берегу Финского залива. Всплывает американская подводная лодка. Оттуда высовывается американский солдат и спрашивает у чукчей: - не видели, где русская подводная лодка? (Ч): - ушла на Северо-Запад! Американцы уплывают. Через некоторое время всплывает русская подводка(Р) с тем же вопросом: (Р)- Американскую подводную лодку не видали? (Ч)- Ушла на Северо-Восток! (Р)- Ты не мудри, ты пальцем покажи!

Я нашел и на русском инфу...
http://xn--c1asr.xn--p1ai/wp-content...D1%8F-ABOS.pdf
Меня просто интересует, в каком она виде получается все таки, или в виде численного уравнения, или в виде координат вершин конечных элементов.

[nikе]

По массиву случайно расположенных точек получаем прямоугольный массив точек(координаты вершин прямоугольной сетки). А потом уже можем получить координату произвольной точки, интерполируя по прямоугольнику с известными координатами вершин.

[DEM]

nikе
ЖАль конечно, что не численном виде получается уравнение поверхности.....

[nikе]

Самое главное чтобы метод работал, а когда нам попадется более совершенный метод, мы и из него извлечем пользу.