TOPREFERAT.COM.KZ - Қазақша рефераттар

войти на сайт

вход на сайт

Логин: :
Пароль :

Забыл пароль Регистрация

Крип құбылысының термофлуктуациялық моделі диплом жұмысы




Крип құбылысының термофлуктуациялық моделі диплом жұмысы
0
Раздел: Математика | Автор: Админ | Дата: 2-05-2015, 10:00
Загрузок: 389







МАЗМҰНЫ - www.topreferat.com.kz

Кіріспе..............................................................................................................4

Кернеу релаксациясын моделдеу...............................................................5

1.1 Феноменологиялық модель......................................................................5

1.2 Кинетикалық модель.................................................................................10

1.3 Молекулалық құрылым бойынша есеп...................................................11

1.4 Есептеулер нәтижесі.................................................................................12

2. Крип құбылысының термофлуктуациялық мо делін жасау...................13

2.1 Амофты қабаттағы крип құбылысының моделін жасау........................13

2.2 Сандық есептеу алгоритмі........................................................................16

2.3 Сандық моделдеу нәтижелері...................................................................17

3. Күйреу деңгейін анықтау............................................................................17

Қорытынды.......................................................................................................18

Қолданылған әдебиеттер.................................................................................19

Қосымша 1.........................................................................................................20

Қосымша 2.........................................................................................................26




Жұмыс түрі: Дипломдық жұмыс
Жұмыс көлемі: - бет
Пәні: Математика

-----------------------------------------------------------------------------------

ДИПЛОМДЫҚ ЖҰМЫСТЫҢ ҚЫСҚАРТЫЛҒАН МӘТІНІ



Мазмұны

Кіріспе..............................................................................................................4

Кернеу релаксациясын моделдеу...............................................................5

1.1 Феноменологиялық модель......................................................................5

1.2 Кинетикалық модель.................................................................................10

1.3 Молекулалық құрылым бойынша есеп...................................................11

1.4 Есептеулер нәтижесі.................................................................................12

2. Крип құбылысының термофлуктуациялық мо делін жасау...................13

2.1 Амофты қабаттағы крип құбылысының моделін жасау........................13

2.2 Сандық есептеу алгоритмі........................................................................16

2.3 Сандық моделдеу нәтижелері...................................................................17

3. Күйреу деңгейін анықтау............................................................................17

Қорытынды.......................................................................................................18

Қолданылған әдебиеттер.................................................................................19

Қосымша 1.........................................................................................................20

Қосымша 2.........................................................................................................26

Реферат

Бұл «Крип құбылысының термофлуктуациялық моделі» деп аталатын магистрлік дисертиция
Қазіргі уақытта конструкциялық полимерлердің кең қолданыста болуына байланысты оларды
Бұл жұмыста кейбір аспектлерге қатысты жаңа болып табылатын криптің
1 қосымшада фортран тілінде жазылған модель теңдеулерін шешіп, күйреу
2 қосымшада магистрлік жұмыстың негізгі нәтижелерінің бірі болып табылатын
Жалпы диссертациялық жұмыстағы негізгі терминдер мыналар: крип, кернеу, температура,
Кіріспе

Жүктелген қатты денелердің бұзылуы туралы қазіргі ілім (оны болжамдау
Ұзақ жүктеу жағдайларында полимерлердің бұзылуына арналған зерттеулер қазіргі уақытта
Термофлуктуациялық бұзылу теориясы негізінде [4-7] – де жылжу жағдайларда
Микроскопикалық мөлшерлерден және желілерде жүктемені кездейсоқ бөлінуіне байланысты материалдың
Термофлуктациялық теорияның шеңберінде тізбек микродеформациясынан және көршілес түйірлердің өзара
Негізгі жеңілдететін жағдай [10] – 4 жүйелі түйірлерден құрылған
Осы жұмыста қорытындыда [10] және «Регель элементі» түрінде аморфты
Полимерлердің механикалық қасиеттері бойынша елеулі ерекшеленетін кристаллды және аморфты
Кернеу релаксациясы кезінде бұзылу аморфты фазада болады және ұзақтығы
Полимерлердің механикалық қасиеттері көптеген қатынастарда металлдардың қасиеттерінен ерекшеленеді. Осы
Поликонденсациялық полимерлерді су шығару арқылы бифункционалды молекулалардың өзара әрекет
Осының бәрі сызықты полимерлердің мысалы, мұнда мономерлы звенолары үзілмейтін
Басты тізбектің бірнеше жерлерінен екінші тізбектер басталуы сияқты тізбектің
Бірінші тарауда релаксацияның негізгі модельдері, екінші тарауда -
1. Кернеу релаксацияны моделдеу

1 Феноменологиялық модельдер

«Кернеу релаксациясы» термині денеде оның тұрақты деформациясы кезінде механикалық
Модельдік ұсыныстар негізінде [1-6, 16] дифференциалды теңдеу алу.

Аргументтің ағымдағы мағынасын ғана емес, сондай-ақ деформалаудың бүкіл тарихын
Жылжудың техникалық теориясы негізінде деформация мен кернеу, олардың өзгеру
Феноменологиялық теориялар, оның ішінде жылжу теориясы суреттелген процесстерді толық
Ең алдымен, теория болжамдау және алдына алу көрудің құралы
Мұралықтың сызықтық теориясы. Мұралықтың сызықты теориясы механикалық полимерлік
(1.1)

мұнда с – тұрақты өлшем, және
(1.1) шартын біртектілік шарты, ал (1.2) - сызықты суперпозиция
(1.2)

(1.3)

мұнда Е - серпімділіктің кенеттік модулі,
R ал функциясы
Белгілі бір функцияның негізінде басқа функцияны белгілеуге мүмкіндік беретін
Бұдан басқа, [3]-те (1.3) теңдеуді жазудың басқа нысаны бар,
(1.4),

(1.3) және (1.4) теңдеулердің айырмашылығы релаксация өзегінің өлшемділігіне байланысты.

Осы теңдеулер бойынша, гер де интегралды теңдеудің (1.4) өзгеі
Теңдеуді дифференциялай (1.4) отырып, R(t) релаксация жылдамдығына пропорционалды екенін
(1.5)

Сонымен, R(t) функциясын белгілеу үшін сараптау түрде белгіленген функциялардың
Вольтер [10] типтегі сызықты уақытша операторларды енгізетін болсақ, Больцманның
(1.6)

Мұнда, функциясы
(1.7)

Осы теңдеулер сыртынан Гука заңына ұқсас, бірақ та Работнов
Жылжу және релаксация функциясының негізінде [4] мынадай теңдеу шығаруға
(1.8)

Бұл үйірткі типтегі интегралдармен ара қатынасты белгілейтін Стилтьес үйірткі
(1.6) және (1.7) теңдеуінен (1.3) және (1.4) интегралды
(1.9)

белгілеп, (1.9) теңдеу шығарамыз.

(1.8) теңдеуді бөлшектеп ықпалдастыра отырып, мынаны шығарамыз:

(1.10)

белгілеп (1.11) теңдеуін шығарамыз.

Кинетикалық модель

1. Созылудың тұрақты деформациясы жағдайларында полимердің аморфты қабатының бұзылуын
-ұзындығы бойынша тізбектерді бөлу тығыздығы, -
Егер де ұзындық бірлігі ретінде аморфтық қабаттың ұзындығын алатын
(1.12)

мұнда аморфтық қабаттың деформациясы болып табылады.

P( l, t) функциясы [7-9] келтірілген теңдеуге бағынады:

(1.13)

коэффициенті W( )=1 мүмкіндігіне жауап беретін тізбектер үзілуінің
[10]-ға сәйкес, 4 түйнектерден тұратын қысқа учаскедегі тізбек үзілуінің
(1.14)

[10]-мен салыстырғанда біз W( ) формуланы деформация терминдарында
Мұнда, - молекула осінің бағытында
( =0 деформация тізбектің тепе-теңдік күйіне жауап береді,
(1.15)

4-түйіндерден құрылған тізбек учаскесінің деформаланбаған ұзындығын арқылы
(1.16)

lmin мөлшерін белгілейік.

мөлшерінің мағынасы бойынша тізбектің минималды ұзындығы теңдікті қанағаттандырады:

Осыдан және (1.13)-ден екені шығады. Деформацияланбаған
(1.17)

(1.16) теңдеу (1.13), (1,14), (1.15), (1.17) бірге
байланысты жеке жағдайды қарастырайық.

Мұндай жағдайда (1.16) жеңіл түрде ықпалдасады:

(1.18)

Жоғарыда біз кристаллиттерден тізбекті шығару нәтижесін елемедік. Сонымен қатар,
Бұзылу функциясын бұзылған тізбектер санының қатынасы арқылы олардың
(1.19)

Тізбекті тарту күші деформациясына үйлесімді деп санайық. Онда
(1.20)

мұнда f0–материалдың константы ретінде қарастыратын күштің өлшемділік көлемі болып
F(t) әрекеті әр l сайын t
(1.19)-(1.20) формулалар аморфтық қабаттағы бәсеңдейтін күшті және бұзылу динамикасын
Соңғы уақытта бірыңғай температура кезінде полимерлерді эксплуатациялау практикасында статикалық
Статикалық қажу құбылысы беріктіктің сынды сипаты туралы ұғынымды қайта
Кинетикалық тұжырымға сәйкес бұзылу мен ұзақ беріктікті жиналудың кейбір
Тұрақты созылу кернеу кезінде жылжымалылық процессінің барысында ұзақ беріктікті
(1.21)

мұнда, - тұрақты созылу кернеу салынған
Мақталы бауды беріктіккке сынау кезінде бұзылғанға дейін уақыттың санды
(1.22)

мұнда, -
Ұзақ мерзміділіктің температурадан тәуелділігін зерттеу қортындылары Бусқа статикалық бұзылулардың
(1.23)

Мұнда, R - Больцманның тұрақтылығы, U -
Ұзақ беріктіктің және қатты денелердің бұзылу процесстерінің, оның ішінде
(1.24)

мұнда, t0 - тұрақты, (10-11–10-13 c) атомдардың жылу
Т – абсолютті температура, U0 – активацияның энергиясы, өлшемі
(1.20) формулаға сәйкес тұрақты кернеу кезінде жарты логарифмикалық координаттарда
Малекулалық қүрылым бойынша есеп

Молекулалық құрылымның параметрлерінің бірі - мөлшері,
- аморфтық және кристалдық фазалардың орташа деформациясы, [16]-да
Сонымен, (а өлшемін макроскопикалық деформациясымен байланыстырғаннан кейін, (1.19)-(1.20) ара
Егер де аморфтық қабат ұзындығының тегіс еместігін елемесек, онда
[13] қорытындыларында негізделе отырып, субмикротрещиналардың сын концентрациясына жауап беретін
Сонда, ((t) = (кр теңдеу ашық емес түрде
үшін формуланы талдай отырып, ұзақ мерзімдік макроскопикалық деформацияға
Релаксация қисық негізінде есеп.

Полимерлер механикасында практикалық есептерде сараптау түрде шығатын қисық жылжулар
Осылар арқылы изохронды қисықтардың жиынтығы құрылады. Изохронды қисықтар мен
(1.12/-(1.14) тапсырма шешуінің итерационды әдісінің схемасы.

q(t)=1+(а(t ) арқылы белгілейік, сонда q(t)
(1.25)

1. (t, ti=i (t уақыт бойынша
(i(l)= ((l, ti.), qi=q(ti ), (0(l)= ((l, 0.), q0=1+(0
2. (i(l) берілгенмен (1.13) қанағаттанатын
3. (1.14) -ден q =qi
4. с i=i+1 2 қадамына
Әр итерация кезінде бұзылудың мағынасын бақылаймыз.

(1.26)

Ұқсастықты тексеру үшін уақыт бойынша есеп жүргізіп, қорытындыларды салыстыру
Есеп айырылысу бағдарламасының қисынды модульдары:

1. Барлық физикалық параметрлерді белгілеу.

2. Ықпалдасу нүктелер санының есеп айырылысу параметрлерін, уақыттың ақырғы
3. (0(l)= ((l, 0.), q0=1+(0 функцияларының бастапқы мәнін
4. (5) теңдеуді шешуге арналған бағдарламаша.

5. q i-ді есептеу

( i(l) есептеу

Келесі итерацияға 5,6 қайталау.

Бұзулықты есептеу. Сын мағына асқаннан кейін итерацияны тоқтату.

Есеп айырылысуды аяқтауға берілген уақытқа жеткеннен кейін, итерацияны тоқтату.

30% дейін уақытты есептеу.

11. бұзылулар

12. Мәліметтердің қортындысы.

Есептеулер нәтижесі

Релаксация құбылысы үшін ұзақ өмірліктің есептері бойынша қорытындылар келтірейік
[10] –да ( t ) есеп
Есептерде ( = 0,5 , ал бұзылуға дейінгі уақыт
(0 мағынасы полиэтилен үшін бағамен келісілген және
Бірқатар параметрлер мағынасының белгісіздігіне қарамастан, есептеу қорытындылары теориялық ұзақ
Мүмкін, бұл жоғары температура кезінде полимердің аморфизация деңгейі өсіп,
2. Крип құбылысының термофлуктуациялық мо делін жасау

2.1 Амофты қабаттағы крип құбылысының моделін жасау

[5] жұмысында [4]-ші жұмыстың нәтижесін қолдана отырып, аморфты
Бұл жұмыста [5] жұмысы негізінде бірөсті созу жағдайындағы
Алынған материалдың құрылысы [5,7-10] жұмыстардағыдай аморфты-кристалды деп қарастырайық. Мұндағы
(2.1)

(2.2)

Осындай теңдеулер [7-9]-жұмыстарында қарастырылған.

коэффициенті ([1]) ықтималдыққа жауап беретін
–бойлық деформациясы шынжырлардың ұзындығына бірқалыпты таралған деп қарастырайық. Онда
(2.3)

Мұндағы – бір элементар аймақтың үзілу ықтималдығы,
(2.4)

Мұнда, - молекула осінің бағытында
(2.4) формуласы кезкелген энергиамен және өзара әсерлесетін, U0
Аморфты қабатша түрақты
(2.5)

мұндағы күштің өлшемдік мөлшері, оны материалдың
Кернеудің тұрақтылық шарты, тасушы шынжырлардың бұзылу ықтималдығы туралы
(2.6)

(2.7)

(2.8)

[5] жұмысындағыдай мұнда да біз кристалиттерден шынжырлардың созылуын ескермейміз,
(2.6)-(2.8) теңдеулерін тасушы шынжырлардың ұзындығы бойынша бстапқы таралуымен толықтыруымыз
(2.9)

[5]-те алынған релаксация қисықтары бастапқы күйге дейінгі лездік деформацияға
Жүктеудің бастапқы тарихын қарастырмай-ақ және сандық есептеу шартының біртүрлілігі
(2.10)

С коэффициенті нормаланған коэффициент:

, мұндағы берілген модельдің параметрлері. [5,10]
шарты (2.5) формуласымен біріге отырып, нөлдік деформация кезінде материалда
Бастапқы таралуды осылай таңдап алу ішкі күшке
Осылайша (2.10)-да есептелген шынжырлардың ұзындығы бойынша бастапқы
2.2 Сандық есептеу алгоритмі

Теңдеуде екі ((l,t), (a(t) белгісіз функциалар және бұзылуға
Мұндай тендеулерді шешудің жалпы теориалық әдістері жоқ. Сол үшін
Белгісіз функциаларды уақыттың ti=ti-1 +(ti, i=0,1,2,3,…тізбектелген кезеңдерінде анықтап аламыз.
1 қадам. (i=(a(ti) белгісіз болсын (1) тендеуінен (a(t)(
.

2 қадам. (2.6)-ға t=ti+1 –ді және алынған ((l,ti+1)-ді
3 қадам. (a(ti+1)=(* деп аламыз. i=i+1, t=ti+1 бастап 1-ші
есептеуде (ti мәні (6)-да кернеудің түрақтылық шарты δ қатынасты
.

ω(t) мәні критикалық бүзылу мәніне жеткенде есептеу тоқталады.

Берілген алгоритм (ti уақыты бойынша аз қадаммен алынған
2.3 Сандық модельдеу нәтижелері

(2.6-2.9) моделі бойынша алғашқы есептеулер жүргізілді және батапқы деформациа
көрсетеді. 3 суретте бастапқы уақыттағы және 30% бұзылу
Модельдің анықталмаған параметрлерінің бірі ω0 өлшемі. [5]-тегі есептеу бойынша
3. Күйреу деңгейін анықтау

Төзімділікті анықтау үшін біз критикалық бұзылуға жеткенге дейінгі
Практикада ([1,11,13]) крип жағдайында төзімділіктің келесі түрдегі формулалары τ=F(σ)exp(U0/kT),
Ln(τ)≈Const(1/T)+f(σ).

Көп жағдайларда ([11,13]), ln(τ) графиктері не бір нүктеден
Формуласымен апроксимацияланса, онда ln(τ)-ның (σ/Т)-ге тәуелділігі мына түрде
ln(τ) )=-γ(σ/T)+f(σ, Т).

Мұндай тәуелділік графиктері теріс иілумен параллель түзулер тобын құруы
Қорытынды

Крипті сипаттаудың ұсынылған модельі теориялық жағынан статистикалық физика заңдарына
ҚОЛДАНЫЛҒАН ӘДЕБИЕТТЕР

В.А.Петров, А.Я.Башкарев, В.И. Веттегрень. Физические основы прогнозирования долговечности конструкционных
Журков С.Н., Нарзуллаев Б.Н. Временная зависимость прочности твердых тел.
Журков С.Н., Томашевский Э.Е.. Исследование прочности твердых тел. //
Муканова Б.Г. Оценка долговечности ориентированной полимерной цепи при релаксации
Муканова Б.Г., Джунисбеков О.Т. Расчет поврежденности и релаксирующего усилия
Бокшицкий М.Н. Длительная прочность полимеров // М.:Химия, 1978-308б.

Регель В.Р., Слуцкер А.И. Томашевский Э.Е. Кинетическая природа прочности
Зайцев М.Г. Связь микро- и макропараметров разрушения в модели
Губанов А.И. Кинематика разрушения полимера // Механика композитных материалов,
Губанов А.И., Кособукин В.А. Оценка распределения проходных полимерных цепей
Ратнер С.Б. Механическое разрушение пластмасс как процесс деструкции полимеров.
Муканова Б.Г., Андрющенко О.В. Моделирование скоростных зависимостей механических свойств
Бартенев Г.М. Прочность и механизм разрушения полимеров. //
14.Аскадский А.А. Деформация полимеров.// М.: Химия, 1973, -448б

15.Бокшицкий М.Н. Длительная прочность полимеров.// М.:Химия, 1978-308б.

16.Крегер А.Ф., Янсон Ю.О. О построении единого спектра времен
17.Бартенев Г.М., Сандитов Д.С. Релаксационные процессы в стеклообразных системах.//
18.Муканова Б.Г. Жумагул Г.О. Моделирование ползучести на основе термофлуктуационной
1 ҚОСЫМША

Программа

MODULE modCr

!******** global parameters*******************

real*8 u0,tem,beta,emi,ema,eps_end,eps_step

real*8 l0,delta,Lmin,Lmax,tmax,psiCr,hl,t0,C0

real*8 nChain

!omega = 4450000000000#, calculated for 100000 sec.

! gaz's universal const }

real*8,parameter :: pi=3.14159265358979,omega = 4.45E+12, R = 8.314

integer*4 NT

integer*4, parameter :: n=40000

real*8 ro_l(0:n),f_l(0:n),ro_l1(0:n),f_l1(0:n)

contains

subroutine initPar

u0=105. !(Kdj/mol)

tem=343.!(Kelvin)

beta=0.3 !layer amorfous degree

nChain=40. !elements on chaun

emi=-0.235 !minimal deformation

ema=0.47 !maximal deformation

l0=1.06 !center of chain's length

delta=0.06

eps_end=0.2

eps_step=0.00002

Lmin=1.

Lmax=1.24

psiCr=0.7!critical damaging

hl=(Lmax-Lmin)/n

end subroutine ! initPar

real*8 function w(w0,n)

real*8 w0,n

If (Abs(w0) > 0.000001) Then

w=1-(1 - w0)**n

Else

w=n*w0*(1-(n-1)/2*w0*(1-(n-2)/3*w0*(1-(n-3)/4*w0)))

End if

end function

real*8 function w0(e)

real*8 a,e, e1,e2,w1,w00

a=Exp(u0*1000./tem/R);

e2=-emi; e1=0.0;

If (E w0=0.0

return

end if

If (E w00 = (E - emi) / (E - e1)

w00 = (w00 * w00 - 3) * a
w00 = w00 + 1

w0 = 1 / w00

return

end if

If (E w00 = (ema - E); w00 = w00 *
w1 = (ema - emi); w1 = w1 *
w1 = (w1 - 4.5 * w00) * a

w00 = w00 * 1.5 * a * a

w1 = w1 + 4.5 * (E - e2)
w1 = w1 + w00

w00 = w00 / w1; w0 = 1 -
return

end if

w0 = 1.

End function

END MODULE !modCr

PROGRAM PrMy

use DFlib

use modCr

integer*4 i,k,nPrm,ntime,ntimemax

real*8 eps,epsmacro,sigma0,x,t,l,S1,S2,S3,S4,epsl,psitmp

real*8 delta_T,delta_tlog,F0,F1,F2,delta_cntrl

real*8 deform(0:100000),psi(0:100000)

logical L_

L_=.true.

nPrm=0

ntimeTmp=100000

!delta_tlog=1./ntimeTmp

delta_t=1./ntimeTmp

tmax=5. !decimal log for time

ntimemax=50000 ntimeTmp=2000

if (nPrm .eq. 0) then

call initPar

print *, 'enter epsilon_start'

read *, epsmacro

else

OPEN(10,FILE='.param.txt')

read (10,101) u0,tem,beta,nChain,emi,ema,ntimemax,sigma0

read (10,102) l0,delta,Lmin,Lmax,psiCr,eps_end,eps_step

read (10,103) epsmacro,ntime,t_Tmp

ntimeTmp=ntime+ntimeTmp

end if

hl=(Lmax-Lmin)/n

l=Lmin

t0= 0.000001!

eps=epsmacro/beta

eps_step=eps_step/beta

if (nPrm .eq. 0) then

!start distribution of chain's lenth

l=Lmin

do i=0,n

if (l > 1. + eps) then

epsl=0.

else

epsl=(1.+eps-l)/l

end if

x = (l - l0)/delta

ro_l(i)= dExp(-x * x)/Sqrt(pi)/delta

ro_l(i)=ro_l(i)*dexp(-w(w0(epsl),nChain*l)*omega*t0)

f_l(i)= ro_l(i)* epsl

l=l+hl;

end do

!calculate normalize factor

S1=0.

do i=1,n-1,2

S1=ro_l(i)+S1

end do

S2=0.

do i=2,n-2,2

S2=ro_l(i)+ S2

end do

C0=(ro_l(0)+ro_l(n)+ S1*4.+S2*2.)*hl/3.

! end calculate normalize factor

!normalize start distribution

do i=0,n

ro_l(i)= ro_l(i)/C0

f_l(i)= f_l(i)/C0

enddo

!calculate start TENSION

S3=0.

do i=1,n-1,2

S3=f_l(i)+S3

end do

S4=0.

do i=2,n-2,2

S4=f_l(i)+ S4

end do

sigma0=(f_l(0)+f_l(n)+ S3*4.+S4*2.)*hl/3.

deform(0)=eps

t=0.;ntime=0;

t=t0

!t=t0*dexp(delta_tlog)

else

read (10,106) (ro_l(i),i=0,n)

read (10,106) (f_l(i),i=0,n)

close(10)

t=t_Tmp

deform(ntime)=eps

endif

k=0

!delta_tlog=dlog(1.D+01)*tmax/ntimemax

!delta_t=dlog(1.D+01)*tmax/ntimemax

delta_eps=eps_step

if (nPrm .eq. 0) then

OPEN(10,FILE='eps.txt')

else

OPEN(10,FILE='eps.txt',ACCESS='APPEND')

endif

do1: DO WHILE(deform(ntime)ntime=ntime+1

eps=deform(ntime-1)

call RO_F(eps,delta_t,F0,psitmp)

delta_cntrl=dabs((F0-sigma0)/sigma0)

do while (delta_cntrl>0.0004)

delta_t=delta_t/2.

call RO_F(eps,delta_t,F0,psitmp)

delta_cntrl=dabs((F0-sigma0)/sigma0)

enddo

ro_l(:)=ro_l1;f_l(:)=f_l1(:);

eps=deform(ntime-1)+delta_eps

F1=F0

do while(F1call RO_F(eps,delta_t,F1,psitmp)

eps=eps+delta_eps

end do

deform(ntime)=deform(ntime-1)+(eps-delta_eps-deform(ntime-1))*(sigma0-F0)/(F1-F0)

call RO_F(deform(ntime),delta_t,F1,psi(ntime))

delta_cntrl=dabs((F1-sigma0)/sigma0)

do while (delta_cntrl>0.0004)

delta_t=delta_t/2.

call RO_F(deform(ntime),delta_t,F1,psi(ntime))

delta_cntrl=dabs((F1-sigma0)/sigma0)

enddo

do while (delta_cntrldelta_t=delta_t*2.

call RO_F(deform(ntime),delta_t,F1,psi(ntime))

delta_cntrl=dabs((F1-sigma0)/sigma0)

enddo

write( 10, '(f14.6,2x,f14.6)') t,deform(ntime)*beta

t=t+delta_t

!t=t*dexp(delta_tlog)

!if (ntime>=ntimemax) exit

if (mod(ntime,200)==0) then

print *, 'f1=',f1,(F1-sigma0)/sigma0,'psi=',psi(ntime)

write( *, 105) deform(ntime)*beta,ntime,t,delta_t

endif

if (psi(ntime)print *, ' Damage, PsiCr= ',psi(ntime)

print *, 'f1=',f1,(F1-sigma0)/sigma0,'psi=',psi(ntime)

write( *, 105) deform(ntime)*beta,ntime,t,delta_t

exit do1

endif

END DO do1

close(10)

101 FORMAT(' U0=', f8.2,' Temp=', f6.1, ' beta=', f4.2,
102 FORMAT(' l0=',f7.3,' delta=' f7.3,' Lmin=',f7.3,' Lmax=' ,f7.3,'
103 FORMAT(' epsmacro=',f14.8,' ntime=',i12,' t=',f14.8)

106 format(' ',20(e14.8))

104 FORMAT(100(E14.6))

105 FORMAT(' def=',F12.6,' ntime=',I8,' t=',F12.4,'
OPEN(11,FILE='param.txt')

write (11,101) u0,tem,beta,nChain,emi,ema,ntimemax,sigma0

write (11,102) l0,delta,Lmin,Lmax,psiCr,eps_end,eps_step

write(11,103) deform(ntime)*beta,ntime,t

write(11,106) (ro_l(i),i=0,n)

write(11,106) (f_l(i),i=0,n)

close(11)

print *, 'Calculated'

END

SUBROUTINE RO_F(eps,t_,F,PSI)

use Modcr, only: f_l,ro_l,f_l1,ro_l1,hl,n,w,w0,nChain,omega,Lmin

real*8 l,q,S1,S2,S3,S4,t_,F,PSI,epsl,eps,tmp1,tmp2,tmp3

integer i

l=Lmin

do i=0,n

if (l > 1. + eps) then

epsl=0.

else

epsl=(1.+eps-l)/l

end if

tmp1=omega*t_;tmp2=w0(epsl);tmp3=w(tmp2,nChain*l)*tmp1

tmp1=dexp(tmp3)

ro_l1(i)=ro_l(i)*dexp(-w(w0(epsl),nChain*l)*omega*t_)

f_l1(i)=ro_l1(i)* epsl

l=l+hl

end do

S1=0.; S3=0.

do i=1,n-1,2

S1=ro_l1(i)+S1; S3=f_l1(i)+S3

end do

S2=0.; S4=0.

do i=2,n-2,2

S2=ro_l1(i)+ S2;S4=f_l1(i)+ S4

end do

F=(f_l1(0)+f_l1(n)+ S3*4.+S4*2.)*hl/3.

PSI=(ro_l1(0)+ro_l1(n)+ S1*4.+S2*2.)*hl/3.

END SUBROUTINE

2 ҚОСЫМША

Суреттерге түсініктеме:

Сурет 1. Әр түрлі температура крип деформациасы Деформация
Сурет 2. Әр түрлі бастапқы Деформациялар үшін крип деформациасы:
Сурет 3. (1) бастапқы уқыттағы тасушы шынжырлардын таралу тығыздығы
Сурет 4. (1) бастапқы уақыттағы шынжырлардың бойымен күштердің таралуы
Сурет 5. Күйреу деңгейінің логарифмінің әр түрлі бастапқы деформациядағы
Сурет 6. Күйреу деңгейі логарифмінің σ0/T-ге ε0-нің мына 0.045(1);
4







Написать комментарий
Имя:*
E-Mail:
Полужирный Наклонный текст Подчеркнутый текст Зачеркнутый текст | Выравнивание по левому краю По центру Выравнивание по правому краю | Вставка смайликов Выбор цвета | Скрытый текст Вставка цитаты Преобразовать выбранный текст из транслитерации в кириллицу Вставка спойлера
Введите код: *


Бұл сайтта Қазақстанның түкпір-түкпірінен жиналған қазақ тіліндегі рефераттар мен курстық және дипломдық жұмыстар ұсынылған. Қазіргі таңда www.topreferat.com.kz сайтының қазақ тіліндегі жұмыстар базасы бүкіл интернеттегі ең үлкен база болып табылады! Біздің базадағы жұмыстар саны 15000-нан асады. Біз бұл жетістікпен тоқтап қалмаймыз! Біз базамызды одан әрі толықтырамыз.
» » Крип құбылысының термофлуктуациялық моделі диплом жұмысы

© 2011-2016 Скачать бесплатно на topreferat.com.kz курсовые, дипломные и рефераты на телефон, на планшет и на компьютер.
При копировании материала активная ссылка на источник обязательна.


Мнение посетителей:
 

После 9 класса Вы:

Пойду в 10, 11, закончу школу полностью
Пойду в Колледж
Пойду в ПТУ
Пойду работать
Снова пойду в 9 класс

 
 
Похожие:
  • Полимерлердің қаттылығын анықтаудағы механикалық қасиеттері және оның ерекш ...
  • Алгебралық теңдеулер жүйесін шешудің кейбір әдістері диплом жұмысы
  • Алгебралық теңдеулер жүйесін құрудың үшбұрышты шекті элементтер әдісі дипло ...
  • Стерженьде серпімді толқындардың таралуы диплом жұмысы
  • Физикалық процесстерді математикалық моделдеу диплом жұмысы
  • Жылу өткізгіштік теориясының негіздері диплом жұмысы
  • Сызықтық теңдеулер жүйесін шешу диплом жұмысы
  • Жоғары дәрежелі теңдеулер диплом жұмысы
  • Электр тізбектер теориясы курстық жұмыс
  • Дифференциалдық теңдеулерді оқытудың әдістемесі курстық жұмыс
  • Қуысты ортада газдың белгіленбеген қозғалысы курстық жұмыс
  • Өзара индуктивтілігі бар тізбектер курстық жұмыс
  • Модуль және оның қасиеттері курстық жұмыс
  • Математиканы тереңдетіп оқытудағы туындының алгебралық қолданылуы курстық ж ...
  • Контурлық токтар әдісі курстық жұмыс
  • Жүктелген параболалық теңдеуді коэффициент арқылы басқару курстық жұмыс
  • Дифференциалдық теңдеулерді мектепте оқыту курстық жұмыс
  • Анықталған интегралды жуықтап шешу курстық жұмыс
  • Полимерлердің жылуфизикалық қасиеттері реферат
  • Теңдеулер жүйесі реферат
  • </