МКЭ. Подскажите с чего начать)

[kolo666 255]

Хочу стать асом, начал читать книжку Бате  "Метод конечных элементов". Подскажите, на верном ли я пути) Какие еще книжки стоит почитать?

[MotorManiac 130]

Почитайте литературу, потом переходите к практическим занятиям....

[Борман 2 397]

Найдите любую толковую книгу с основами МКЭ. Потом попробуйте вручную решить задачу с одним 4-х узловым плоским солидом. Начните с аппрксимации перемещений в элементе, постройте функции формы и матрицу жесткости. Решите задачу.

После этого вы будете асом на фоне 95% здесь присутствующих.

Наиболее низкоуровневый доступ к КЭ формулирвке задачи предоставляет классический Ансис. В целом - название пакета не имеет значения.

Удачи.

[kolo666 255]

Найдите любую толковую книгу с основами МКЭ. Потом попробуйте вручную решить задачу с одним 4-х узловым плоским солидом. Начните с аппрксимации перемещений в элементе, постройте функции формы и матрицу жесткости. Решите задачу.

После этого вы будете асом на фоне 95% здесь присутствующих.

Наиболее низкоуровневый доступ к КЭ формулирвке задачи предоставляет классический Ансис. В целом - название пакета не имеет значения.

Удачи.

Спасибо! Будет чем заняться на новогоднии праздники)

[soklakov 1 485]

Потом попробуйте вручную решить задачу с одним 4-х узловым плоским солидом.

Может лучше с плоских ферм начинать?

[Борман 2 397]

Может лучше с плоских ферм начинать?

Не... теряется кухня интегрирования "по шаблону". Но тоже можно, конечно.

 

Так же очень полезно научится вручную приводить давление и объемные силы к узлам.

 

После всего этого очень многое прояснится.

[kolo666 255]

Посоветуйте пожалуйста ТОЛКОВУЮ книгу) или книги.. С практикой, или даже задачками с ответами..если такие есть конечно.

[soklakov 1 485]

Мне, помнится, лекции Маслова понравились. Кажется, здесь их в приличном объеме цитируют.

[Fedor 1 630]

http://pnu.edu.ru/media/filer_public/2013/04/10/6-13_segerlind_1979.pdf   там и коды есть. Можете набрать на с например и порешать. В генераторе ошибочка есть i и j в одном цикле надо поменять местами :) 

[AlexKaz 1 851]

В генераторе ошибочка есть i и j в одном цикле надо поменять местами

//
    // GENERATION OF THE REGION NODE NUMBERS
    //
    KN1=1
    KS1=1
    KN2=NROWS
    KS2=NCOL
    for i=1:4
        NRT=JT(NRG,i)
        if (NRT==0)|(NRT>NRG) then 
            continue
        else
            for j=1:4
                if (JT(NRT,j)==NRG) then NRTS=j
                end
            end
            k=NCOL
            if (i==2)|(i==4) then k=NROWS
            end
            JL=1
            JK=ICOMP(i,NRTS)
            if (JK==-1) then 
                JL=k
            end
            for j=1:k
                if (i==1) then
                    NN(NROWS,j)=NNRB(NRT,NRTS,JL)
                    KN2=NROWS-1
                end
                if (i==2) then
                    NN(j,NCOL)=NNRB(NRT,NRTS,JL)
                    KS2=NCOL-1
                //mprintf('debug NN %d   NRT %d    NRTS %d    JL %d  i %d  j %d\n',NN(j,NCOL),NRT,NRTS,JL,i,j) 
                end
                if (i==3) then
                    NN(1,j)=NNRB(NRT,NRTS,JL)
                    KN1=2
                end
                if (i==4) then
                    NN(j,1)=NNRB(NRT,NRTS,JL)
                    KS1=2
                end
                JL=JL+JK
            end
        end
    end

Случаем не на этом участке?

[AlexKaz 1 851]

После замены

                if (i==2) then
                    NN(j,NCOL)=NN(j,1)//NNRB(NRT,NRTS,JL)
                    KS2=NCOL-1
                end

тестовый пример обрабатывается корректно. На другом разбиении не проверял.

 

ELEMENT GENERATION FOR TORSION PROBLEM

GLOBAL COORDINATES

NUMBER X COORD Y COORD
1 0.000000 0.000000
2 0.125000 0.000000
3 0.250000 0.000000
4 0.400000 0.000000
5 0.500000 0.000000
6 0.500000 0.100000
7 0.500000 0.250000
8 0.500000 0.375000
9 0.500000 0.500000
10 0.425000 0.425000
11 0.350000 0.350000
12 0.300000 0.300000
13 0.250000 0.250000
14 0.375000 0.250000
15 0.200000 0.200000
16 0.150000 0.150000
17 0.075000 0.075000
18 0.250000 0.125000
CONNECTIVETY DATA
REGION SIDE
SIDE 1 2 3 4
1 2 0 0 0
2 0 0 1 3
3 0 2 0 0
*** REGION ****
ROWS, COLUMNS, BOUNDARY NODE NUMBERS
REGION NODE NUMBERS
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
NEL NODE NUMBERS
1 6 7 1 0.300000 0.300000 0.350000 0.318750 0.350000 0.350000
2 7 2 1 0.350000 0.318750 0.387500 0.387500 0.350000 0.350000
3 7 8 2 0.350000 0.318750 0.400000 0.337500 0.387500 0.387500
4 8 3 2 0.400000 0.337500 0.425000 0.425000 0.387500 0.387500
5 8 9 3 0.400000 0.337500 0.450000 0.356250 0.425000 0.425000
6 9 4 3 0.450000 0.356250 0.462500 0.462500 0.425000 0.425000
7 9 10 4 0.450000 0.356250 0.500000 0.375000 0.462500 0.462500
8 10 5 4 0.500000 0.375000 0.500000 0.500000 0.462500 0.462500
9 11 12 6 0.250000 0.250000 0.312500 0.250000 0.300000 0.300000
10 12 7 6 0.312500 0.250000 0.350000 0.318750 0.300000 0.300000
11 12 13 7 0.312500 0.250000 0.375000 0.250000 0.350000 0.318750
12 13 8 7 0.375000 0.250000 0.400000 0.337500 0.350000 0.318750
13 13 14 8 0.375000 0.250000 0.437500 0.250000 0.400000 0.337500
14 14 9 8 0.437500 0.250000 0.450000 0.356250 0.400000 0.337500
15 14 15 9 0.437500 0.250000 0.500000 0.250000 0.450000 0.356250
16 15 10 9 0.500000 0.250000 0.500000 0.375000 0.450000 0.356250
*** REGION ****
ROWS, COLUMNS, BOUNDARY NODE NUMBERS
REGION NODE NUMBERS
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
NEL NODE NUMBERS
17 16 17 12 0.250000 0.187500 0.317188 0.182812 0.312500 0.250000
18 16 12 11 0.250000 0.187500 0.312500 0.250000 0.250000 0.250000
19 17 18 13 0.317188 0.182812 0.381250 0.178125 0.375000 0.250000
20 17 13 12 0.317188 0.182812 0.375000 0.250000 0.312500 0.250000
21 18 19 14 0.381250 0.178125 0.442188 0.173437 0.437500 0.250000
22 18 14 13 0.381250 0.178125 0.437500 0.250000 0.375000 0.250000
23 19 20 15 0.442188 0.173437 0.500000 0.168750 0.500000 0.250000
24 19 15 14 0.442188 0.173437 0.500000 0.250000 0.437500 0.250000
25 21 22 17 0.250000 0.125000 0.321875 0.118750 0.317188 0.182812
26 21 17 16 0.250000 0.125000 0.317188 0.182812 0.250000 0.187500
27 22 23 18 0.321875 0.118750 0.387500 0.112500 0.381250 0.178125
28 22 18 17 0.321875 0.118750 0.381250 0.178125 0.317188 0.182812
29 23 24 19 0.387500 0.112500 0.446875 0.106250 0.442188 0.173437
30 23 19 18 0.387500 0.112500 0.442188 0.173437 0.381250 0.178125
31 24 25 20 0.446875 0.106250 0.500000 0.100000 0.500000 0.168750
32 24 20 19 0.446875 0.106250 0.500000 0.168750 0.442188 0.173437
33 26 27 22 0.250000 0.062500 0.326562 0.057813 0.321875 0.118750
34 26 22 21 0.250000 0.062500 0.321875 0.118750 0.250000 0.125000
35 27 28 23 0.326562 0.057813 0.393750 0.053125 0.387500 0.112500
36 27 23 22 0.326562 0.057813 0.387500 0.112500 0.321875 0.118750
37 28 29 24 0.393750 0.053125 0.451562 0.048438 0.446875 0.106250
38 28 24 23 0.393750 0.053125 0.446875 0.106250 0.387500 0.112500
39 29 30 25 0.451562 0.048438 0.500000 0.043750 0.500000 0.100000
40 29 25 24 0.451562 0.048438 0.500000 0.100000 0.446875 0.106250
41 31 32 27 0.250000 0.000000 0.331250 0.000000 0.326562 0.057813
42 31 27 26 0.250000 0.000000 0.326562 0.057813 0.250000 0.062500
43 32 33 28 0.331250 0.000000 0.400000 0.000000 0.393750 0.053125
44 32 28 27 0.331250 0.000000 0.393750 0.053125 0.326562 0.057813
45 33 34 29 0.400000 0.000000 0.456250 0.000000 0.451562 0.048438
46 33 29 28 0.400000 0.000000 0.451562 0.048438 0.393750 0.053125
47 34 35 30 0.456250 0.000000 0.500000 0.000000 0.500000 0.043750
48 34 30 29 0.456250 0.000000 0.500000 0.043750 0.451562 0.048438
*** REGION ****
ROWS, COLUMNS, BOUNDARY NODE NUMBERS
REGION NODE NUMBERS
36
37
11
38
39
16
40
41
21
42
43
26
44
45
31
NEL NODE NUMBERS
49 38 39 36 0.112500 0.112500 0.181250 0.150000 0.150000 0.150000
50 39 37 36 0.181250 0.150000 0.200000 0.200000 0.150000 0.150000
51 39 16 37 0.181250 0.150000 0.250000 0.187500 0.200000 0.200000
52 16 11 37 0.250000 0.187500 0.250000 0.250000 0.200000 0.200000
53 40 41 38 0.075000 0.075000 0.162500 0.100000 0.112500 0.112500
54 41 39 38 0.162500 0.100000 0.181250 0.150000 0.112500 0.112500
55 41 21 39 0.162500 0.100000 0.250000 0.125000 0.181250 0.150000
56 21 16 39 0.250000 0.125000 0.250000 0.187500 0.181250 0.150000
57 42 43 40 0.037500 0.037500 0.143750 0.050000 0.075000 0.075000
58 43 41 40 0.143750 0.050000 0.162500 0.100000 0.075000 0.075000
59 43 26 41 0.143750 0.050000 0.250000 0.062500 0.162500 0.100000
60 26 21 41 0.250000 0.062500 0.250000 0.125000 0.162500 0.100000
61 44 45 42 0.000000 0.000000 0.125000 0.000000 0.037500 0.037500
62 45 43 42 0.125000 0.000000 0.143750 0.050000 0.037500 0.037500
63 45 31 43 0.125000 0.000000 0.250000 0.000000 0.143750 0.050000
64 31 26 43 0.250000 0.000000 0.250000 0.062500 0.143750 0.050000
BANDWIDTH QUNTITY IS 27 CALCULATED IN ELEMENT 52

[AlexKaz 1 851]

Плюс с заменой 

                if (i==3) then
                    NN(1,j)=NN(NROWS,j)//NNRB(NRT,NRTS,JL)
                    KN1=2
                end

получилось пристроить 4-ю зону под 3-й.

 

И всё же какой-то костыль... У Сегерлинда есть 2-е издание, многое переписано и дополнено, хотя просматривая GRID в ней не заметил.

Изменено 14 февраля 2015 пользователем AlexKaz

[SashaL 5]

А есть ли что менее объемное, пусть и менее подробное по МКЭ, например страниц на 5-10?

[AlexKaz 1 851]

Ну если это не первоапрельская шутка, я ещё в школе распечатывал шпоры по истории, 16 листов в одном. Думаю уложитесь в десяточку =)

У Сегерлинда введение для плоских элементов довольно короткое, если серьёзно.

[Борман 2 397]

И то правда.. в этих МКЭ-шных книжках писанины бестолковой больше чем полезной информации. Нам вот препод лекции по МКЭ читал - все коротенько и понятно.

[AlexKaz 1 851]

в мкэшных книгах примеры разжеваны. а нужна голая теория с лекций - так есть Theory reference ансисовского хэлпа. из одних голых формул талмуд накопился, толще Сегерлинда

[Борман 2 397]

Разжевана балка - это да. А чтоб обычный плоский элементик с кокретными числами - так это мало где найдешь. А ведь если один раз самому проделать ЭТО с плоским элементиком, то все станет ясно.

[AlexKaz 1 851]

И то правда.. в этих МКЭ-шных книжках писанины бестолковой больше чем полезной информации. Нам вот препод лекции по МКЭ читал - все коротенько и понятно.

Подойдёт?

Отчёт 4 - МКЭ, статика, плоская.pdf

[Цветочек 15]

Как решать линейные задача вообще без проблем объяснят в любом вузе, но с тех пор как зенкевич написал свою первую книгу прошло пятьдесят лет. Наибольшая проблема заключается в том, чтобы на основе мизерного времени въехать в современные технологии МКЭ расчетов, на фоне которых L-координаты и сопромат в задачах выглядят позавчерашним днем

Рекомендую книгу на английском, автор Askhar Bhatti.

Изменено 5 января 2016 пользователем Цветочек

[AlexKaz 1 851]

Ну-у, возвращаясь к ТС, есть Бате издание 2010 года, причём на русском. В нём, если копнуть, много что можно нарыть. Но кодов конечно нет. По модалке на русском для начинающих есть разве что Зенкевич (свежие издания куда полнее пятидесятилетних ), Бате (самое старое и самое новое переведённые) и Образцов. В них же есть и нелинейность, куда прибавляется Хечумов. Так что жить можно =)