| 3 | 1/1 | 返回列表 |
| 查看: 742 | 回復(fù): 2 | ||
proejhy木蟲 (小有名氣)
|
[求助]
為什么三輥卷板的這個例子在我的機(jī)器里總也不收斂呢? 已有1人參與
|
|
非線性分析綜合應(yīng)用實例—鋼板卷制成圓筒 /CLEAR !清除數(shù)據(jù)庫,新建文件 /FILNAM, EXAMPLE20 !指定工作名 /CONFIG, NRES,2000 !設(shè)置最大子步數(shù) !前處理 /PREP7 !進(jìn)入前處理器 /PNUM, VOLU,ON !打開體號 ET,1,SHELL43 !選擇單元類型,殼單元用于劃分上下輥 ET,2,SOLID186 !實體單元用于劃分鋼板 MP,EX,1,2E11 !定義材料模型1的彈性模量 MP,DENS,1,7800 !定義材料模型1的密度 MP,NUXY,1,0.3 !定義材料模型1的泊松比 MP,EX,2,2E11 !定義材料模型2的彈性模量 MP,DENS,2,7800 !定義材料模型2的密度 MP,NUXY,2,0.3 !定義材料模型2的泊松比 TB,BKIN,2,1 TBTEMP,0 TBDATA, ,240E6,2E8 !定義材料模型2的屈服極限、切向模量 R,1,0.02 !定義實常數(shù)1,殼單元的厚度 CYLIND,0.38/2,0,0.3,2.3,0,360 !創(chuàng)建上輥,圓柱體 CYLIND,0.38/4,0,0.3,0,0,360 CYLIND,0.38/4,0,2.6,2.3,0,360 VGLUE,ALL !粘接 VSEL,NONE WPOFF,0.54/2,-0.4 !偏移工作平面原點(diǎn)到右下輥中心 CYLIND,0.3/2,0,0.3,2.3,0,360 !創(chuàng)建右下輥,圓柱體 CYLIND,0.3/4,0,0.3,0,0,360 CYLIND,0.3/4,0,2.6,2.3,0,360 VGLUE,ALL VSEL,NONE WPOFF,- 0. 54 !偏移工作平面原點(diǎn)到左下輥中心 CYLIND,0.3/2,0,0.3,2.3,0,360 !創(chuàng)建左下輥,圓柱體 CYLIND,0.3/4,0,0.3,0,0,360 CYLIND,0.3/4,0,2.6,2.3,0,360 VGLUE,ALL ALLS !選擇所有實體 VDELE,ALL !刪除體,但保留面 HPTCREATE,AREA,5,,COORD !在上輥端面中心處創(chuàng)建硬點(diǎn) HPTCREATE,AREA,10,,COORD,0,0,2.6 HPTCREATE,AREA,34,,COORD,-0.54/2,-0.4,2.6 !在左下輥端面中心處創(chuàng)建硬點(diǎn) HPTCREATE,AREA,29,,COORD,-0.54/2,-0.4,0 HPTCREATE,AREA,22,,COORD,0.54/2,-0.4,2.6 !在右下輥端面中心處創(chuàng)建硬點(diǎn) HPTCREATE,AREA,17,,COORD,0.54/2,-0.4,0 MAT,1 !為上下輥劃分單元指定屬性,材料模型 TYPE,1 !單元類型 REAL,1 !實常數(shù) ESIZE,0.05 !指定單元邊長度 SMRTSIZE,8 !指定智能尺寸級別 MSHAPE,1 !指定單元形狀為三角形 AMESH,ALL !對所有面劃分單元 BLOCK,0.7,-1.0,0.3/2+.02,0.3/2+.02+0.03,0.5,2.1 !創(chuàng)建鋼板,塊 LESIZE,96,,,1 !指定線96(鋼板厚度)被劃分為1段 MAT,2 !為鋼板劃分單元指定屬性 TYPE,2 MSHKEY, 1 !映射網(wǎng)格 MSHAPE,0 !指定單元形狀為六面體 ESIZE,0.05 !指定單元邊長度 VMESH,ALL !對塊劃分網(wǎng)格 ET,4,TARGE170 !指定單元類型,用于創(chuàng)建接觸對 ET,5,CONTA174 KEYOPT,5,12,1 !設(shè)置單元5接觸表面無滑動 R,10,,,0.1 !定義實常數(shù),KFN=0.1 R,11,,,0.1 R,12,,,0.1 ASEL,S,,,40 !在鋼板上表面和上輥表面建立接觸對 NSLA,S,1 NSEL,U,LOC,Z,0.8,1.8 MAT,3 REAL,10 TYPE,4 ESURF ALLS ASEL,S,,,3,4,1 NSLA,S,1 NSEL,U,LOC,Z,0.9,1.7 NSEL,U,LOC,Z,0.3,0.4 NSEL,U,LOC,Z,2.2,2.3 MAT,3 REAL,10 TYPE,5 ESURF ALLS ASEL,S,,,39 !在鋼板下表面和左下輥表面建立接觸對 NSLA,S,1 NSEL,U,LOC,Z,0.8,1.8 MAT,4 REAL,11 TYPE,4 ESURF ALLS ASEL,S,,,27,28,1 NSLA,S,1 NSEL,U,LOC,Z,0.9,1.7 NSEL,U,LOC,Z,0.3,0.4 NSEL,U,LOC,Z,2.2,2.3 MAT,4 REAL,11 TYPE,5 ESURF ALLS ASEL,S,,,39 ! 在鋼板下表面和右下輥表面建立接觸對 NSLA,S,1 NSEL,U,LOC,Z,0.8,1.8 MAT,4 REAL,12 TYPE,4 ESURF ALLS ASEL,S,,,15,16,1 NSLA,S,1 NSEL,U,LOC,Z,0.9,1.7 NSEL,U,LOC,Z,0.3,0.4 NSEL,U,LOC,Z,2.2,2.3 MAT,4 REAL,12 TYPE,5 ESURF ALLS FINI !首先模擬下輥不動,上輥下降 /SOLU !進(jìn)入求解器 D,ALL,UZ !在所有節(jié)點(diǎn)上施加約束 D,ALL,ROTX D,ALL,ROTY NSEL,S,LOC,X !選擇上輥軸線上節(jié)點(diǎn) NSEL,R,LOC,Y D,ALL,UY,-0.08 !施加約束,上輥下降0.08 m D,ALL,UX D,ALL,UZ ALLS CSYS,4 !激活工作平面坐標(biāo)系 ASEL,S,LOC,Y,-2,0.3/2+0.01 !選擇下輥上所有面 NSLA,S,1 !選擇下輥上所有節(jié)點(diǎn) D,ALL,ALL !約束掉所有自由度 ALLS ANTYPE,TRANS !瞬態(tài)分析 NEQIT,100 !指定最大迭代次數(shù) LNSRCH,ON !打開線性搜索 NLGEOM,ON !打開大變形選項 TIME,1 !指定載荷步時間 AUTOT,ON !打開自動載荷步長 NSUBST,30,80,20 !指定子步數(shù)目 KBC,0 !斜坡載荷 OUTRES,ALL,ALL !輸出所有子步所有項目的結(jié)果 SOLVE !解答 !以下模擬下輥轉(zhuǎn)動,驅(qū)動鋼板前進(jìn) CSYS,4 ASEL,S,LOC,Y,-2,0.3/2+0.01 !選擇下輥上所有面 NSLA,S,1 !選擇下輥上所有節(jié)點(diǎn) DDELE, ALL,ALL !刪除掉上一載荷步施加在下輥上的所有約束 ALLS D,ALL,UZ !在所有節(jié)點(diǎn)上施加約束 D,ALL,ROTX D,ALL,ROTY NSEL,S,LOC,X,0 !選擇左下輥軸線上所有節(jié)點(diǎn) NSEL,R,LOC,Y D,ALL,ROTZ,-3.1415926 !轉(zhuǎn)動半周 D,ALL,UX !約束掉移動 D,ALL,UY ALLS WPOFF,0.54 !偏移工作平面 NSEL,S,LOC,X,0 !選擇右下輥軸線上所有節(jié)點(diǎn) NSEL,R,LOC,Y D,ALL,UX !約束掉移動 D,ALL,UY D,ALL,ROTZ,-3.1415926 !轉(zhuǎn)動半周 TIME,2 ALLS KBC,0 NSUBST,200,1000,100 SOLVE !解答 !模擬卸載 ACEL, 0,9.8 !施加重力加速度 CSYS,0 !激活全球直角坐標(biāo)系 NSEL,S,LOC,X !選擇上輥軸線上所有節(jié)點(diǎn) NSEL,R,LOC,Y D,ALL,UY,-0.06 !上輥向上移動,卸載 D,ALL,UX D,ALL,UZ D,ALL,ROTX D,ALL,ROTY ALLS KBC,0 NSUBST,20,300,10 TIME,2.1 SOLVE FINI !查看結(jié)果 /POST1 SET,LAST !讀最后載荷子步計算結(jié)果 PLNS,U,SUM ANTIME,50,0.5, ,1,2,0,2.1 !用動畫查看卷制圓筒過程中的變形情況 FINI |
金蟲 (小有名氣)
|
這個例子我也做過,本來就不收斂,一直在收斂線上下震蕩,不過還是能看到結(jié)果,準(zhǔn)確度不敢恭維。書上例子很多都是這樣的。 發(fā)自小木蟲Android客戶端 |

木蟲 (小有名氣)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +6 | 李李不服輸 2026-03-26 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +6 | 景隨風(fēng) 2026-03-29 | 9/450 |
|
|
[考研] 309求調(diào)劑 +7 | 誰不是少年 2026-03-29 | 7/350 |
|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
[考研] 各位老師好,我的一志愿為北京科技大學(xué)085601材料專碩 +8 | Koxui 2026-03-28 | 8/400 |
|
|
[考研] 329求調(diào)劑 +7 | 星野? 2026-03-26 | 7/350 |
|
|
[考研] 394求調(diào)劑 +3 | 好事多磨靜候佳?/a> 2026-03-26 | 5/250 |
|
|
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +3 | nhj_ 2026-03-25 | 3/150 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
[考研] 085404求調(diào)劑,總分309,本科經(jīng)歷較為豐富 +4 | 來財aa 2026-03-25 | 4/200 |
|
|
[考研] 307求調(diào)劑 +8 | 超級伊昂大王 2026-03-24 | 9/450 |
|
|
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分,英語已過四六級,六級521,求調(diào)劑 +4 | 美味可樂雞翅 2026-03-26 | 4/200 |
|
|
[考研] 312求調(diào)劑 +9 | 上岸吧ZJY 2026-03-22 | 13/650 |
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
[考研] 生物技術(shù)與工程 +3 | 1294608413 2026-03-25 | 4/200 |
|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|
[考研] 一志愿南航材料專317分求調(diào)劑 +5 | 炸呀炸呀炸薯條 2026-03-23 | 5/250 |
|
|
[考研]
|
13659058978 2026-03-24 | 4/200 |
|