| 13 | 1/1 | 返回列表 |
| 查看: 983 | 回復(fù): 12 | |||
[交流]
求助求助,模擬鉆桿與井壁摩擦
|
|||
|
打算用直接建模法,來建立鉆桿軌跡,然后用pipe59模擬鉆桿,combin40來模擬鉆桿和井壁的摩擦,在鉆桿某個位置產(chǎn)生振動,來得出振動所影響的范圍。但這個振動不是正弦波,所以用瞬態(tài)進行了分析。在下是ansys新手,所以編的apdl一直不收斂,還請各路大神看一下,給出寶貴建議,感激不盡。。。 !建模 finish /clear /PREP7 /ANG,1 /Units,SI !設(shè)置單位為國際單位 *set,zj_pipe,300 !直井段鉆桿的長度 *set,numb_zj_pipe,34 !直井段鉆桿的根數(shù) *set,wq_pipe,100 !彎曲井段鉆桿的長度 *set,numb_wq_pipe,17 !彎曲井段鉆桿的根數(shù) *set,sp_pipe,135 !水平段鉆桿的長度 *set,numb_sp_pipe,18 !水平井鉆桿的根數(shù) *set,drill_collar ,165 !水平段鉆鋌的長度 *set,numb_drill_collar,15 !水平段鉆鋌的根數(shù) *set,R,100 !圓弧半徑為955 *set,pipe_OD,0.127 !鉆桿外徑 *set,pipe_thick,0,0.0092 !鉆桿壁厚 ?這里單位是米么 *set,drill_collar_OD,0.178 !鉆鋌外徑 *set,drill_collar_thick,0,0.0535 !鉆鋌壁厚?這里單位是米么 *set,drilling_fluid_density,1300 !鉆井液密度 *set,drilling_fluid_length,2300 !鉆井液深度 *set,EXs,2.1e11 !彈性模量 *set,NUXYs,0.3 !泊松比 *set,Density,7800 !鉆桿密度 *set,g,9.8 *set,zoomrc,500 !鉆柱系統(tǒng)截面擴大系統(tǒng) *set,CTR10a,0.25 !alpha系數(shù)(質(zhì)量阻尼系數(shù)) *set,CTR10b,0.3 !beta系數(shù)(剛度阻尼系數(shù)) *set,Kspring1,9.85E6 !減震器剛度 *set,Mspring,485 !減震器質(zhì)量 *set,Kspring2,5E6 !減震器剛度2 *set,Mbit,42 !鉆頭質(zhì)量 *set,M,7200 !井口設(shè)備質(zhì)量 *set,Bit_pressure,160000 !鉆壓 /PREP7 !前處理 ET,1,PIPE59 !定義管單元 ET,2,combin40 keyopt,2,1,0 !間隙能力 keyopt,2,3,1 !單元自由度 keyopt,2,4,0 !單元輸出 keyopt,2,6,1 !質(zhì)量位置 R,2,100000,0.12,100,0.1,1000,1000 !R,2,K1,C,M,GAP,FSLIDE,K2 R,1,pipe_OD,pipe_thick,,,drilling_fluid_density,0 !定義1號鉆桿管單元的真實值,鉆桿外徑,鉆桿壁厚,鉆井液密度 Rmore,,,,CTR10a !定義鉆鋌切向阻力系數(shù) TB,WATER,1 TBDATA,3,drilling_fluid_density TBDATA,4,drilling_fluid_length MP,EX,1,EXs !設(shè)置鉆桿的彈性模量 MP,NUXY,1,NUXYs !設(shè)置鉆桿泊松比 MP,DENS,1,Density !設(shè)置鉆桿密度 !建模 /pnum,node,1 N,1,0,0,0 !建立節(jié)點編號為1 N,2,0,-1,0 !建立節(jié)點編號為2 N,35,0,-zj_pipe,0 !設(shè)置鉆桿長度 zicd=300 有33根,從節(jié)點充滿, FILL,2,35,,3 K,1,0,-zj_pipe,0 !按坐標(biāo)值生成一個關(guān)鍵點 K,2,wq_pipe,-zj_pipe,0 K,3,wq_pipe,-(wq_pipe+zj_pipe),0 LARC,1,3,2,R, !由關(guān)鍵點1和3,中心點2,半徑為995生成弧線。 LESIZE,1,,,numb_wq_pipe !直線1劃分165等分 功能:對所選擇的線設(shè)置網(wǎng)格單元大小 LMESH,1 !給弧線1畫網(wǎng)格 TYPE,1 REAL,1 E,1,2 !使節(jié)點2號到3號為單元屬性為1,實常數(shù)為1的鉆桿 EGEN,33,1,-1 !復(fù)制單元屬性,復(fù)制34次,增量為1,最大編號為-1 TYPE,1 REAL,1 E,34,36 !使節(jié)點2號到3號為單元屬性為1,實常數(shù)為1的鉆桿 E,36,38 N,54,wq_pipe,-(wq_pipe+zj_pipe),0 !定義水平段鉆桿的節(jié)點坐標(biāo),根數(shù)為18,單元為管單元,實常數(shù)為R2 N,92,446,-(wq_pipe+zj_pipe),0 FILL,54,92,,55 !將定義好的坐標(biāo)劃分18段,即水平鉆桿的根數(shù) nummrg,node,,,,high !兩個重合的點選擇數(shù)值大的點 TYPE,1 REAL,1 E,38,39 !使節(jié)點2號到3號為單元屬性為1,實常數(shù)為1的鉆桿 EGEN,54,1,-1 !復(fù)制單元屬性,復(fù)制34次,增量為1,最大編號為-1 N,93,wq_pipe,-(wq_pipe+zj_pipe),0 !定義水平段鉆桿的節(jié)點坐標(biāo),根數(shù)為18,單元為管單元,實常數(shù)為R2 N,131,446,-(wq_pipe+zj_pipe),0 FILL,93,131,,94 !將定義好的坐標(biāo)劃分39段,即水平鉆桿的根數(shù) TYPE,2 REAL,2 *DO,i,0,38 E,93+i,54+i *enddo KWPAVE,2 !移動工作平面到關(guān)鍵點2 FINISH /SOL ALLSEL,ALL !選擇所有實體 CSYS,0 !把坐標(biāo)系換成笛卡爾坐標(biāo)系 D,1,ALL D,2,UX,,,34,,UZ,ROTX,ROTY,ROTZ !約束垂直段所有節(jié)點的除Y方向的位移 !約束彎曲段所有節(jié)點的除Y方向的位移 D,54,UY,,,92,,UZ,ROTX,ROTY,ROTZ !保持X方向的位移 D,36,UX,,,,,UZ,ROTX,ROTY,ROTZ D,38,UX,,,53,,UZ,ROTX,ROTY,ROTZ *DO,i,93,131 d,i,ALL *enddo FINISH /solu ACEL,0,-9.8,0 f,1,fy,-100000!鉆壓 pred,off AUTOTS, 1 NSUBST,100,10000000,,on antype,static pstres,on allsel,all,all solve finish !瞬態(tài) /solu ACEL,0,9.8,0 antype,4 TRNOPT,FULL BETAD,0.02 !剛度矩陣乘子 ALPHAD,0.02 !質(zhì)量矩陣乘子 pstres,on ALLSEL,ALL NT=200 dt=0.00025 *do,i,1,NT f,75,fx,300000*sin(2*3.14*dt*i) time,dt*i nsubst,1,,,1 kbc,0 outres,all,all solve *enddo time,dt*NT+2 NSUBST,100,10000000,,on kbc,0 solve finish /post1 etable,F1,smisc,1 !彈簧1上的力 etable,F2,smisc,2 !彈簧2上的力 etable,U1,nmisc,5 !節(jié)點i的位移 etable,U2,nmisc,6 !節(jié)點j的位移 etable,stat,nmisc,1 !單元狀態(tài) nsle nplot pretab,F1,F2,U1,U2,stat |
» 搶金幣啦!回帖就可以得到:
+1/231
+2/224
+1/156
+2/136
+5/125
+1/87
+1/86
+1/86
+1/82
+1/34
+1/23
+1/12
+1/10
+1/9
+1/6
+1/6
+1/6
+1/5
+1/3
+1/3

| 13 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 歡迎211本科同學(xué),過A區(qū)國家線,A區(qū)非偏遠(yuǎn)一本,交叉學(xué)科課題組 +14 | lisimayy 2026-03-04 | 24/1200 |
|
|---|---|---|---|---|
|
[考研] 新疆大學(xué)地質(zhì)與礦業(yè)工程學(xué)院招生 +5 | another12 2026-03-04 | 10/500 |
|
|
[考研] 301求調(diào)劑 +3 | 朝天椒的雙馬尾 2026-03-05 | 3/150 |
|
|
[考研] 328求調(diào)劑 +6 | vuzhdkfjkx 2026-03-04 | 7/350 |
|
|
[考博] 26申博-目前4篇SCI一作 +4 | chen_2024 2026-03-02 | 4/200 |
|
|
[考研] 求調(diào)劑,學(xué)校研究所都可以,材料與化工267分 +6 | wmx1 2026-03-05 | 6/300 |
|
|
[考研] 271求調(diào)劑 +7 | 月色c 2026-03-05 | 8/400 |
|
|
[考研] 0703 東華大學(xué) 理學(xué)化學(xué)均過 a 線,總分 281 求調(diào)劑 +4 | 香香 Lu 2026-03-03 | 8/400 |
|
|
[基金申請]
|
xhuama 2026-03-02 | 17/850 |
|
|
[考研] 一志愿武漢理工大學(xué)-085602-總分296分-求調(diào)劑 +7 | 紫川葡柚 2026-03-04 | 7/350 |
|
|
[考研] 0703化學(xué) 一志愿東北大學(xué) 279分求調(diào)劑 +4 | 聶樂平 2026-03-03 | 4/200 |
|
|
[考研] 281求調(diào)劑 +3 | 我是小小蔥蔥 2026-03-03 | 5/250 |
|
|
[考研]
085600 英一數(shù)二272求調(diào)劑
5+6
|
vida_a 2026-03-01 | 47/2350 |
|
|
[考研] 291求調(diào)劑 +4 | Afy123456 2026-03-03 | 7/350 |
|
|
[考研] 292求調(diào)劑 +3 | sgbl 2026-03-03 | 3/150 |
|
|
[考研] 0805總分292,求調(diào)劑 +12 | 幻想之殤 2026-03-01 | 12/600 |
|
|
[考研] 291求調(diào)劑 +3 | MuoLuo1312 2026-03-02 | 6/300 |
|
|
[考研] 298求調(diào)劑 +10 | 人間唯你是清歡 2026-02-28 | 14/700 |
|
|
[考研] 295求調(diào)劑 +8 | 19171856320 2026-02-28 | 8/400 |
|
|
[考研]
|
LYidhsjabdj 2026-02-28 | 4/200 |
|