Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / САПР изделий из композиционных материалов. Моделирование процессов деформирования и разрушения в среде ANSYS

.pdf
Скачиваний:
24
Добавлен:
19.11.2023
Размер:
68.45 Mб
Скачать

Рис. П2.1. Расчетная схема

Последовательность действий через ГПИ:

– описать скалярные переменные, описывающие геометрию конст-

рукции; меню:

 

 

 

Utility Menu

Parameters

Scalar Parameters...

 

Ввести в строке ввода: pi = 3.1416 Accept

Dr = 10 Accept

Ar =

=60 Accept

B2r = 40 Accept Jpol = pi*Dr**4/32 Accept

Close

– задать тип КЭ; меню:

 

 

 

Main Menu

Preprocessor

Element Type

Add/Edit/Delete

Add

Beam 3D elastic 4;

– задать геометрические характеристики сечения (двутавр 20а): площадь поперечного сечения (AREA = pi*dr**2/4), момент инерции (IZZ = jpol/2) и полярный момент инерции (IXX = jpol) балок в форме вещест-

венных констант (рис. П2.2); меню:

 

 

 

 

Main Menu

Preprocessor

Real Constants

Add/Edit/Delete

Add ...

– описать свойства материала; модуль упругости и модуль сдвига

(ÌÏà); ìåíþ:

 

 

 

 

 

 

 

Main Menu

Preprocessor

Material Props

Constant Isotropic

OK

EX2e

6

 

 

 

 

 

 

GXY8e

5

 

 

 

 

 

 

Для более поздних версий ANSYS соответственно:

 

Main Menu

Preprocessor

Material Props

Material Models

Structural

Linear

Elastic

Orthotropic;

 

 

 

– создать узлы по координатам на плоскости; меню:

 

Main Menu

Preprocessor

 

Create

Nodes

In Active CS

NODE

Node Number

 

1

 

 

 

X,Y,Z

Location in active CS

0

0

0

Apply

NODE

Node Number

 

2

 

 

 

X,Y,Z

Location in active CS

Ar

0

0

Apply

181

Рис. П2.2. Диалоговое окно описания вещественных констант

NODE

Node Number

3

 

 

 

 

X,Y,Z

Location in active CS Ar+B2r

0

0

Apply

NODE

Node Number

4

 

 

 

 

X,Y,Z

Location in active CS Ar+B2r*2

0

0

Apply

NODE

Node Number

5

 

 

 

 

X,Y,Z Location in active CS

(Ar+B2r)*2

0

0

OK

– создать линейные элементы по парам узлов в последовательности:

1 2; 2 3; 3 4; 4 5; ìåíþ:

 

 

 

 

 

Main Menu

Preprocessor

Create

Elements

 

Auto Numbered

Thru Nodes

выбрать узлы 1 и 2

Apply; выбрать узлы 2 и 3

Apply è ò.ä.

 

 

 

 

 

 

182

– задать условия закрепления: заделка в узлах 1 и 5 (ограничение по перемещениям и углу поворота); меню:

Main Menu Preprocessor (Solution) (Define) Loads Apply Displacements On Nodes выбрать узлы 1 и 5 All DOF OK;

задать условия нагружения крутящие моменты в узлах 2 и 4; меню:

Main Menu Preprocessor (Solution) (Define) Loads Apply Force/Moment On Nodes выбрать узлы 2 и 4 OK в диалоговом окне в поле Lab выбрать MX, в поле VALUE внести значение: 8e4 OK;

запустить решатель; меню:

Main Menu Solution Solve Current LS;

– просмотреть опорные реакции; варианты меню:

Main Menu General Postproc List Results Reaction Solu MX; Utility Menu List Results Reaction Solution MX;

– просмотреть угол закручивания; меню:

Main Menu General Postproc List Results Nodal Solution DOF solution Rotation ROTX;

Utility Menu List Results Nodal Solution DOF solution

Rotation ROTX. Соответствующий данному диалоговому способу командный метод приведен в листинге П2.1.

/com,

ПРОГРАММА РЕШЕНИЯ ЗАДАЧИ СОПРОТИВЛЕНИЯ

/com,

МАТЕРИАЛОВ ПО ТЕМЕ: “КРУЧЕНИЕ”

/com,

 

 

/com,

Сборник задач по сопротивлению материалов /

/com,

Под ред. В.К. Качурина. М.: Наука, 1970.

/com,

(по задаче 4.33.

Ñ.94.)

/com,

 

 

/com,

Итоговые результаты записываются

/com,

в стандартный файл вывода *.out

 

/com,

 

 

/prep7

! Входим в препроцессор

 

 

! Все данные в системе СГС

 

pi=3.1416

 

dr=10

! dr - диаметр стержня

 

ar=60

! ar, b2r - геометрические параметры

b2r=40

 

 

183

jpol=pi*dr**4/32

! константа для полярного момента

 

! инерции

et,1,beam4

! балочные КЭ ВЕАМ4

r,1,pi*dr**2/4,jpol/2 ! REAL CONST :

!площадь поперечного сечения AREA=pi*dr**2/4 ,

!момент инерции IZZ=jpol/2

rmore„jpol

! полярный момент инерции

 

! (8-ая REAL-константа для ВЕАМ4)

Mp,ex,1,2e6 ! модуль Юнга ЕХ=2е6

mp,gxy,1,8å5

! модуль сдвига GXY=8е5

n,1

! Определяем узлы по координатам

n,2,ar

 

n,3,ar+b2r

 

n,4,ar+b2r*2

 

n,5,(ar+b2r)*2

 

*do,i,1,4

! Определяем элементы ВЕАМ4

e,i,i+1

 

*enddo

 

finish

 

/solu

! Входим в решатель

antype,static

! Статический тип анализа

d,1,all,0„5,4

! Задаем условия закрепления в узлах

f,2,mx,8å4

! Задаем крутящие моменты в узлах

f,4,mx,8å4

 

outpr,basic,1

! Вывод основных результатов в стан-

 

! дартный файл

solve

! Решаем СЛАУ МКЭ

finish

 

/prep7

! Вариант с КЭ РIРЕ16

et,1,pipe16

! КЭ - трубки РIРЕ16

r,1,dr,dr/2

! REAL CONST :

 

! внешний диаметр сечения OD=dr,

 

! толщина стенок трубы TKWALL=dr/2

 

! (сплошная труба)

finish

 

/solu

 

solve

 

184

finish /exit

Листинг П2.1. Текст программы расчета стержня на кручение

В результате выполнения программы в файле jobname.out (jobname имя создаваемых рабочих файлов) среди различной выводной информации найдем значения реактивных моментов МХ = 8е4 (кг·см) и угол закручивания ROTX = 0.611155е 2 (рад), что совпадает с ответом к задаче.

Повторить решение, заменив КЭ BEAM4 на КЭ PIPE16, используя добавление КЭ в список типов КЭ и переназначение атрибутов всех КЭ (см. главу 3).

Приложение 3

Расчет ферм [16]

В качестве примера расчета фермы рассмотрим задачу ¹ 4.67 из задачника [17]. В данной задаче требуется определить опорные реакции и усилия в стержнях фермы, показанной на рис. П3.1 вместе с действующими в узлах силами.

Отнесем ферму к декартовой системe координат 0XY, и по рисунку, приведенному в задачнике, рассчитаем координаты ее узлов. На рис. П3.1 узлы фермы обозначены метками n1,…,n5. Эти узлы имеют следующие координаты (в метрах): n1 (0; 0), n2 (2; 0), n3 (4; 0), n4 (2.5;3/2), n5 (1;3).

Как известно, ферма есть конструкция, состоящая из стержней, соединенных между собой шарнирами (узлами фермы). Поскольку силы, действующие на ферму, прилагаются в шарнирах, то все стержни фермы испытывают только осевые усилия растяжения-сжатия. Следовательно, подходящими КЭ для стержней ферм являются стержневые элементы LINK1 в случае плоских ферм и LINK8 в случае пространственных ферм. Узлы конечно-элементной модели будут тогда совпадать с узлами фермы, а каждый стержень будет отдельным КЭ.

Рис. П3.1. Расчетная схема конструкции

Приведенная на рис. П3.1 ферма с количеством стержней å = 7 и количеством узлов n = 5 является статически определимой, так как для нее выполняется условие статической определимости плоских ферм: å = 2n3. Данную задачу можно решать методами теоретической механи-

186

ки, и тогда не существенны никакие определяющие параметры стержней, кроме их длин. Однако для КЭ LINK1 (или LINK8) требуется задать, по крайней мере, одно материальное свойство (модуль Юнга ЕХ) и одну вещественную константу КЭ (площадь поперечного сечения АRЕА). Примем произвольно, что для всех стержней фермы ЕХ = 2·1011 (Í/ì2) è AREA = 1·104 2). Отметим, что в рассматриваемой статически определимой задаче значения этих параметров не будут влиять на итоговые искомые величины, подлежащие определению.

Последовательность действий через ГПИ:

– описать скалярные переменные F_R (имя файла для вывода резуль-

татов расчета) и S3 (значение

3); ìåíþ:

 

 

Utility Menu Parameters

Scalar Parameters...

 

Ввести в строке ввода: F_R=’SMs1’

Accept

S3=SQRT(3)

Accept Close;

 

 

 

задать тип КЭ; меню:

Main Menu Preprocessor Element Type Add/Edit/Delete Add Link 2D spar 1;

– задать площадь поперечного сечения стержней в форме вещест-

венной константы; меню:

 

 

 

 

 

Main Menu

Preprocessor

Real Constants Add/Edit/Delete

Add

AREA

1.e4;

 

 

 

 

 

– описать свойства материала; меню:

 

 

 

Main Menu

Preprocessor

Material Props ... ÅÕ

2å11.

 

Для более поздних версий ANSYS соответственно:

 

Main Menu

Preprocessor

Material Props

Material Models

Structural

 

Linear Elastic

Isotropic;

 

 

 

– создать узлы по координатам на плоскости; меню:

 

Main Menu

Preprocessor

Create

Nodes

In Active CS;

 

– создать линейные элементы по парам узлов в последовательности:

1 2; 2 3; 3 4; 4 5; 5 1; 5 2; 2 4; ìåíþ:

 

 

 

Main Menu

Preprocessor

Create

Elements

Auto Numbered

Thru

Nodes выбрать узлы 1 и 2

Apply; выбрать узлы 2 и 3 Apply è ò.ä.

Созданные элементы будут иметь следующие

атрибуты

(другие

не описаны): номер типа элемента ТУРЕ=1 (LINK1); номер вещественной константы REAL=1 (AREA); номер набора материальных свойств ÌÀÒ=1 (ÅÕ);

187

– задать условия закрепления: подвижная шарнирная опора в узле 1 (ограничение по перемещениям по оси Y) и неподвижный шарнир в уз-

ле 3 (ограничение по перемещениям по осям Õ è Y); ìåíþ:

 

Main Menu

Preprocessor

(Solution)

(Define) Loads

Apply

Displacements

On Nodes

выбрать узлы 1 и 3

UY Apply âû-

брать узел 3

UX

OK;

 

 

 

 

– задать условия нагружения описать усилия, приложенные к уз-

ëàì; ìåíþ:

 

 

 

 

 

 

Main Menu

Preprocessor

(Solution)

(Define) Loads

Apply

Force/Moment

On Nodes

выбрать узлы 1 и 3

Apply

в диалого-

вом окне выбрать FY и в поле Value внести значение: 1å3

Apply

выбрать узлы 4 и 5

Apply

в диалоговом окне выбрать FY и в поле

Value внести значение: 2å3

OK;

 

 

 

– запустить решатель; меню:

 

 

 

Main Menu

Solution

Solve Current LS;

 

 

для записи результатов в файл или окно вывода используйте меню:

Utility Menu File Switch Output to File;

Utility Menu File Switch Output to Output Window;

просмотреть опорные реакции; варианты меню:

Main MenuGeneral PostprocList ResultsReaction SoluAll items; Utility Menu List Results Reaction Solution All items;

– просмотреть усилия в стержнях; варианты меню:

Main Menu General Postproc List Results Element Solution LineElem results;

Utility Menu List Results Element Solution LineElem results. Соответствующий данному диалоговому способу командный метод

(листинг П3.1):

F_R=’SMs1’

! Определяем скалярные переменные

S3=SQRT(3)

 

/PREP7

! Входим в препроцессор

ET,1,LINK1

! Задаем тип КЭ

R,1,1.å-4

! Задаем площадь сечения

ÌÐ,ÅÕ,1,2Å11

! Задаем модуль упругости

N,1,0,0

! Создаем узлы

N,2,2,0

 

188

N,3,4,0

 

N,4,2.5,S3/2

 

N,5,1,S3

 

Å,1,2

! Создаем элементы по узлам

Å,2,3

 

Å,3,4

 

Å,4,5

 

Å,5,1

 

Å,5,2

 

Å,2,4

 

FINISH

! Выходим из препроцессора

/SOLU

! Входим в решатель

ÀNÒÓÐÅ,SÒÀÒIÑ

! Тип анализа - статический

 

! Определяем шарнирные опоры

 

! в узлах

D,1,UY,0„3,2

! UY=0 в узлах 1 и 3

D,3,UX

! UX=0 â óçëå Ç

 

! Задаем силы в узлах

F,1,FY,-1ÅÇ„3,2

! FY=-1еЗ в узлах 1 и 3

F,4,FY,-2ÅÇ„5

! FY=-2еЗ в узлах 4 и 5

SOLVE

! Запуск решателя

FINISH

! Выходим из решателя

/ÐÎSÒ1

! Входим в постпроцессор

/OUTPUT,F_R,txt

! Направляем вывод в файл SMs1.txt

PRRSOL

! Печатаем опорные реакции

PRESOL,SMISC,1

! Печатаем усилия в стержнях

/OUTPUT

! Закрываем файл для записи

FINISH

! Выходим из постпроцессора

Листинг П3.1. Текст программы расчета фермы

Перечень команд удобно записать в файл текстового вида построчно (рекомендуется файл с расширением *.txt) и считать его с помощью меню:

Utility Menu File Read Input from...

Сравнение полученных результатов с ответами к задаче ¹ 4.67 [10] (см. табл. П3.1 и П3.2) позволяет оценить правильность работы программы и точность расчетов ферм с помощью ANSYS.

189

 

 

 

Ò à á ë è ö à Ï 3 . 1

Реакции в узлах

 

 

 

 

 

 

 

Номер узла

1

 

3

Реакции в узлах (кН)

Ответ из задачника

3.25

 

2.75

 

Решение ANSYS

3.25

 

2.75

Ò à á ë è ö à Ï 3 . 2

Усилия в стержнях

Номер стержня

1

2

3

4

5

6

7

Усилия (кН)

Ответ из задачника

1.3

3.03

3.5

2.5

2.6

1.73

1.73

 

Решение ANSYS

1.299

3.0311

3.5

2.5

2.581

1.732

1.7321