35復(fù)合材料三維重構(gòu)及力學(xué)性能的有限元數(shù)值分析
本文關(guān)鍵詞:復(fù)合材料三維重構(gòu)及力學(xué)性能的有限元數(shù)值分析,由筆耕文化傳播整理發(fā)布。
山東大學(xué)博十學(xué)位論文;‰,=菇;t螄,…擴p;同樣求出p-1次曲線的節(jié)點矢量甜川={“f-。唬剩铮玻辏毂牐常ǎ矗┥暾垉(nèi)存空間,用于存放待選;forf=p+ltondo;a=M占1+(1一pg,)+(u7一甜一/7-1;b=u7+∥,+(≯一一u7)=(1-per)+;玩《--inputknotclosettou7;以蟊∈(Cl,6)島=玩elsesf=
山東大學(xué)博十學(xué)位論文
‰,=菇;t螄,…擴p
同樣求出p-1次曲線的節(jié)點矢量甜川={“f-!yUlP-Io??甜。p+-。I)。
Jo2jl睜3)(4)申請內(nèi)存空間,用于存放待選定的節(jié)點矢量S=p。,置,…,Sn…。),令2…2S,2to’Sn+l2Sn.22…5Sn+p+l。t^。
forf=p+ltondo
a=M占1+(1一pg,)+(u7一甜一/7-1)=(1-per)+妒+∥+z蕩1
b=u7+∥,+(≯一一u7)=(1-per)+u7+p卯+“■
玩《--inputknotclosettou7
以蟊∈(Cl,6)島=玩elsesf=z壚
end’
(5)根據(jù)所求丁,s來插值數(shù)據(jù)點Q(f=0,1,…,”),得到曲線c(“),合并節(jié)點矢量u,s并賦值
給u。
該算法預(yù)先定義了P次理想節(jié)點矢量“,以及p—1次節(jié)點矢量“,.,然后判斷候選節(jié)點矢量中是否有處于柔性間距范圍內(nèi)的節(jié)點,如果有,就選擇最接近于理想節(jié)點的存入應(yīng)用于本次曲線的節(jié)點矢量;否則,理想節(jié)點被選中存入應(yīng)用于本次曲線的節(jié)點矢量;應(yīng)用獲得的節(jié)點矢量進(jìn)行曲線插值,在插值曲線完成后,進(jìn)行節(jié)點合并處理,合并后的候選節(jié)點增加了一些新的理想節(jié)點數(shù)據(jù)。
那么根據(jù)這樣的曲線插值方法,可以得到曲面的插值方法:
(1)選擇序列曲線中具有最多型值點的幾條曲線,將其參數(shù)化,求出節(jié)點矢量,然后將其
平均得到最初的公共節(jié)點矢量;
(2)采用各層的型值點按照上述的曲線插值方法進(jìn)行插值,并更新公共節(jié)點矢量:(3)用最終獲得的公共節(jié)點矢量使各條曲線相容,方法是采用節(jié)點插入算法;
(4)獲取縱向的控制頂點和控制節(jié)點矢量。
5.1.2節(jié)點插入【101l
各層曲線相容化的方法是插入節(jié)點。在采用柔性間距的方法獲取了最終的公共節(jié)點矢量后,就需要進(jìn)行相容化處理。其方法就是插入節(jié)點,相應(yīng)的頂點增加相同的個數(shù)。
設(shè)已知一條k次B樣條曲線:
c(甜)=∑d,N1,k(“)(54)
第5章復(fù)合村抖的幾何葷構(gòu)
其中,B樣條基由節(jié)點矢量
U=po,托,…,%+“l】
完全決定,F(xiàn)在要在曲線定義域某個節(jié)點區(qū)間內(nèi)插入一個節(jié)點“e[ut,%。]:Iuk,l/n!浚谑堑玫叫碌墓(jié)點矢量:
U=【甜o,“l,?一,/dr,甜,“f+l,?一,U月+I+I】
重新編號后為:
U=【瓦,玩,一,巧,Z“,玩+2,…,/一ln+l+2】
這個新的節(jié)點矢量驢決定了一組新的B樣條基冠。(“)“=o,I,…,一+I)。原來的B樣條-.15線就由這組新的B櫸條基與未知新頂點Z(,=0,1,…,n+1)表示出來:
n+l
c(“)=∑Z瓢(“)
』,0
控制頂點增加一個,曲線形狀和連續(xù)性均保持不變.
伯姆1134j的節(jié)點插入算法給出了這些未知的新頂點的計算公式:
t=d,
刁“=叫=(1一口:)t+口:乃+,=I|O^¨—.,t一¨川一r一
咖東‰(=00=。)‘55
“j.h1一Hj“\j。
ai=dH_,=i一,十l,…,打+I
其中r表示所插節(jié)點“在老節(jié)點矢量U中的重復(fù)度,若甜。<U<Ⅳ。,則,=0。若重復(fù)度r<k,則有“;“。=“,一l=¨一“…l。這個算法實質(zhì)上就是DeBoor算法求||}次B樣條曲線上的一點c0),“∈(Ua,U!浚恪,“!康牡谝患夁f推。當(dāng),=0時,它僅僅涉及節(jié)點序列”J-k+l,U,4+2'…“m和控制頂點序列d“,dHm…,d,,生成新頂點Z-I+l互^+2,…,互取代老頂點西m”dH+2,…,dj.J。
5.1.3在逼近條件下的序列輪廓公共節(jié)點矢量的確定
在第三章中我們獲得了一定精度條件下的序列輪廓。由于各層輪廓長度的不同,所以在相同的精度條件下,得到的輪廓的控制點是不同的。但是為了能夠使序列輪廓能夠相容,應(yīng)該使控制點數(shù)相同。所以我們的輪廓的相容步驟為:
(I)選擇一定的控制精度下,輪廓的控制點數(shù)最多的輪廓的控制點數(shù)m+l做為標(biāo)準(zhǔn),其余的輪廓的控制點數(shù)調(diào)整為m+l,其中首末控制點重合,然后按照3.6.2節(jié)的方法重新進(jìn)行輪廓的逼近,這樣得出的序列輪廓的精度大于或等于選定的控制精度。
-83.
(2)求出P次曲線的節(jié)點矢量“,=tug,“,,…,“量:,),因為是封閉輪廓,所以其計算方法和開曲線有明顯的不同,其計算方法為:
①重新確定控制頂點的上界值,r/一Jw+P一,,,其中,為頂點的重復(fù)度,在此為1;②按照如。,=dj(j=o,19?a-9k—r)確定I—r+1個重頂點,得到定義該封閉曲線的控制頂點d,“=0,1,…,Ⅲ,m+l,…,帕;
③計算數(shù)據(jù)點參數(shù)丁=f,o,‘,…,o},計算方法可以是累計弦長參數(shù)化法;
④得到節(jié)點矢量U={U0,%,…,U。訓(xùn)),其中定義域外的節(jié)點為:
{U0,“I,…,Up_l}={/'gn-p+t—l,/./月.p+2—1,…,Ⅳn—l}
{/'/n+2,Ⅳ^+,,…,‰+p+1)={1+Up+l,l+Up+2,…j+”2p)
定義域內(nèi)的節(jié)點為:
t,/p=0
,l—I
鏟÷∑o,j-卜Pi=p+l,p+2,…,n
%+l=1
(3)平均各層輪廓的節(jié)點矢量得到初始的節(jié)點矢量:玎=慨,甄,…,酩咖};
(4)更新節(jié)點矢量:
對于每層的輪廓的參數(shù)T={,o。t…,‘},如果存在若干個整數(shù)k∈【p,蚪】滿足參數(shù)f,∈(“。一日,‰+‘)大于(p一2),那么選擇其中節(jié)點跨度k一““j最大的那一個,在其中加入一個新的節(jié)點Ⅳ=扛“+u/)/2,直至所有的,,∈(‰一l一句,‰+‘)小于等于p一2)為止f137】。最終得到公共節(jié)點矢量u=佤,玩,…,_扁,,其中Ct,占,定義為研=pk—l一“¨l和占,=爿”M~‰I,其中P的取值范圍在【o.3,1.o】o
5.1.4逼近條件下的依據(jù)公共節(jié)點矢量的輪廓
在上一節(jié)獲得了公共的節(jié)點矢量,一般情況下,公共節(jié)點矢量比各層輪廓各自的節(jié)點矢量要大,但是各層輪廓的控制點沒有增加,這樣控制點數(shù)少于公共節(jié)點矢量數(shù),就無法獲得相應(yīng)的輪廓曲線。
文獻(xiàn)【137l針對開曲線的公共節(jié)點矢量數(shù)多于型值點數(shù)的情況給出了一種解決方案。獲取一條用公共節(jié)點矢量和少于其的型值點插值出的B樣條曲線使的曲線的能量最小。設(shè)曲線為:
E夠)=f(口p吖叩ffc”汗)dr(5-6)其中c協(xié)和c’,(f,分別是曲線c和)的一階和二階導(dǎo)數(shù),非負(fù)值口,∥稱為剛性和彈性系斛
第5節(jié)復(fù)合材料的幾何事構(gòu)
數(shù).作為B樣條曲線,C(t)可以被表示為:
C(t1=Ⅳ’X
其中,Ⅳ=LⅣo’,,Ⅳl∥…,%】是B樣條基函數(shù),工=[bo,bz,…,k】是控制頂點,那么在式(5-6)中的能量函數(shù)為:
E(c)=£(石)=X7(orfN'N’7dt+#J'NW。7出k=x7KX(5-7)式中Ⅳ(‘)是B樣條基函數(shù)的后階導(dǎo)數(shù),N‘幻---[瞄,(,),Ⅳ乙(,),…,以,,(f)】。
陣,是一個對稱的帶寬為2p—I的矩陣。其具體形式為:置稱作剛性矩
Kq=g叫:。pq、rN?i,吣+pNt:p∞?N?。tj辦t(5-8)蜀,采用Oaussian求積分的形式獲得.控制頂點可以通過最小化二次函數(shù)獲得;minE(Ⅳ)=xm
(st)AX=P
其中St=subject(5.9)tO,引入Lagrange乘子v=Ivo,vl,…,v。1,,將該最小化問題轉(zhuǎn)化為解下面的線性方程:
一][W]
點矢量的曲線。
新的控制頂點。p吶求解上面的方程就可以得出是曲線的能量最小并且過型值點且節(jié)點矢量為公共節(jié)對于逼近條件下的輪廓曲線,我們的步驟是:根據(jù)獲得的公共節(jié)點矢量,取逼近曲線c(“)上具有參數(shù)值氍(f=o,l,…,最}的點c(瓦)。根據(jù)3.3.2節(jié)的方法反算出輪廓曲線的
5.2構(gòu)建一個顆粒及顆粒的顯示
圖5.2.圖5.7給出了~個顆粒的12層輪廓切面,本節(jié)的控制誤差為
如。ma。x(。rain∥q—c(Ⅳ枷_。r~羅mi。.jq,一c(“)}
其中,。,在最大控制誤差選擇3.2的情況下,其中顆粒的第4和第5層輪廓的控制點數(shù)最多,為16個;第11層輪廓的控制點數(shù)最少為12個,其余輪廓的控制點個數(shù)為13.15不等。將輪廓的控制點數(shù)都統(tǒng)一為16個,得到的控制多邊形和逼近曲線如圖5,8.5.13所示。其誤差如表5.1所示。
山東大學(xué)博I一學(xué)位論文
表5.1統(tǒng)一控制頂點數(shù)為16時的各輪廓控制精度
Table5.1Thetoleranceofthecontourswhenthenumberofcontrolpoints缸16I輪廓l2345
E2.93.172.53.23.2
oo
圖5.2顆粒的第1,2層輪廓.
Fig.5.2Thelst.2ndcontoursofthegrain
oo
圖5.4顆粒的第5,6層輪廓
Fig.5.4The5t11.6thCOrltOUtSofthegrain
(二>(二>
圈5,6顆粒的第9,lO層輪廓
Fig,Y.6The9th。10thcontoursofthegrain6789101l12l2.62.62.52.872.882.452.4Ioo圖5,3顆粒的第3,4層輪廓Fig.5.3The3rd.4thcofltoursofthegrainoo圈5.5顆粒的第7,S層輪廓圈5.7顆粒的第ll。12層輪廓Fig.5.7TheIlth,12th‘;oi2toursofthegrain
下載地址:35復(fù)合材料三維重構(gòu)及力學(xué)性能的有限元數(shù)值分析_圖文.Doc
【】最新搜索
復(fù)合材料三維重構(gòu)及力學(xué)性能的有限元數(shù)值分析_圖文
修改錯別字練習(xí)16
臨床實驗室管理52
(0498)《幼教科研方法》網(wǎng)上作業(yè)題及答案
護(hù)理不安全因素分析及防范71
38利用多個LED控制器實現(xiàn)自動漸進(jìn)調(diào)光
數(shù)據(jù)結(jié)構(gòu)順序表算法38
三億文庫
第1題 橋梁加固工程的交工驗收由()組織驗收
96有機硅耐溫漆1
本文關(guān)鍵詞:復(fù)合材料三維重構(gòu)及力學(xué)性能的有限元數(shù)值分析,由筆耕文化傳播整理發(fā)布。
本文編號:210523
本文鏈接:http://sikaile.net/kejilunwen/cailiaohuaxuelunwen/210523.html