方板的强度分析.doc

YUNQI (云起)我现在做一个方板的强度分析的题，但碰到一个问题像请教各位： 由于板上有均布孔（400个）采用1/4模型还是有100个，这给网可划分带来很多麻烦-产生网格畸变（是这样说的么？） 事实上，我现在还只是在算16个孔的情况，（我老师希望我先试一试计算机能算到多少个孔）。但是，在求解16个空的时候，采用smart6.0网格划分就会在求解时出现警告：警告多少个单元的jacobi值超出了软件默认的值-30，之后，求解就会失败， 问题一：请问出现这种警告应该如何优化网格单元？如何划分网格？ （若采用人工划分网格，如何划分孔的边缘即一个圆（直径：15mm）最好分几份？esize=?才好呢？我试着分了几次都不成功，还是求解失败！） 我自己想了一个办法：将模型分割将有圆孔的一块分割开来，成为两个实体，分别划分网格，之后用add加起来，但是，出现：“单元已划分，不能删除” 问题二：请问这样的话，是否会对求解产生影响即会使得求解的结果不正确？ 问题三：采用1/4模型（含100个孔）应该如何划分网格，才能使计算机（奔4）求解不至于失败？ 模型尺寸：1000mm*1000mm*100mm 圆孔直径：15 （400个）受均布力=3e6N/受力面积 谢谢！bigbird圓孔是否規則排列.如果是.應該不難 作一個圓孔.再COPY成100個應該就可以. 以圓來說.圓弧分六份就很準了.再細當然是可以啦 下例是以2D為例子.3D長出即可. 先COPY再長出. M_c(circle).M_s.M_b是mesh參數. 其他改一改參數就知道了 _ finish /clear,start /prep7 et,1,42 mp,ex,10e6 w=1 l=1 r=0.5 M_c=6 M_s=6 M_b=M_c/2 k, k,w k,-w kgen,2,1,3,l kgen,2,1,3,-l l,1,2 l,2,5 l,5,4 l,4,1 l,4,6 l,6,3 l,3,1 l,3,9 l,9,7 l,7,1 l,7,8 l,8,2 a,1,2,5,4 a,1,4,6,3 a,7,1,3,9 a,7,8,2,1 k,100 circle,1,0.5,10,360 lsel,13,16 lesize,all,M_c al,all allsel asba,all,5 lsel,loc,x,0 lsel,a,loc,y,0 lesize,all,M_s allsel esize,M_b lccat,2,3 lccat,5,6 lccat,8,9 lccat,11,12 amesh,allYUNQI (云起)谢谢你，bigbird!可是我的问题还是没有解决，真的又要麻烦你了，谢谢你！现在我把图贴上来，这是我用sweep网格划分的结果，（我老师认为这种网格挺好）但是现在的问题是，我无法加载和加约束载荷是应该加在400*400的这个小区域上的，可是用了SWEEP之后，这个表面上只有一个面积A1（整个表面）呀，请问如何才能让均布载荷只加在指定的400*400的区域上呢？ 另外，我的约束应该加在底面500*500-450*450的一个外围面积上，同样底面也只有一个面积呀！ 我想过在A1上再造一个面积400*400的，可是，没有网格划分载荷是加不上去的，而且还有那么多圆孔。而且，我总觉得问题不会是这样麻烦的。一定有一种简单的处理方法。可我一直没有找到。 bigbird兄，我想问一下，是否sweep网格划分真的只能有一个源面？对于形状不规则的实体，在一个平面上存在几个面积的情况就不能用sweep划分么？不知你能否讲一讲这种网格划分方法的操作步骤？我找过一些书，但都没有找着具体的操作方法。 另外我想，对于我这个实体是否可以考虑用混合网格划分就是在圆孔处用一种方式，在周边用另一种方式划分网格？ 但我还不知如何用这种混合网格划分曾经有一本书上提到一句，但那书没有给出应该怎样操作，你能讲一讲么？ 最近为这个问题，高的焦头烂额，总是求解失败，而且，别人都说是网格划分的问题，所以，哎！还有100个孔更不会做了！惨！bigbird.在YUNQI還沒把圖傳上來時. 先把我做的傳上來.100個不是問題. 可是求解時內存要夠大就是FOR.zbh.163是SHELLL.我用45做的100個FOR.YUNQI花了一個不算短的時間把你的問題做了一下. 先看看你是否要採用.光四個孔就這麼密了.100! 還有你問到關於加載問題.個人建議可畫分為如圖對於你所問的圓柱開孔問題.個人覺得已是進階範例. 故將以積分1做限制. finish /clear,start /prep7 et,1,42 ET,2,45 mp,ex,10e6 w_r=20 w_2r=40 w_sr=5 bh=60 sh=30 hh=10 M_s=6 M_P=6 M_b=M_P/2 M_k=M_b/2 M_h=2 K, CIRCLE,1,w_sr CIRCLE,1,w_r CIRCLE,1,w_2r L,1,2,M_P L,1,3,M_P L,1,4,M_P L,1,5,M_P A,1,2,3 A,4,1,3 A,5,1,4 A,5,2,1 L,2,6,M_b L,3,7,M_b L,4,8,M_b L,5,9,M_b A,2,6,7,3 A,2,5,9,6 A,5,4,8,9 A,4,3,7,8 L,7,11,M_k L,8,12,M_k L,9,13,M_k L,6,10,M_k A,6,7,11,10 A,9,6,10,13 A,8,9,13,12 A,7,8,12,11 LSEL,U,9,24 LESIZE,ALL,M_s ALLSELTYPE,2 VEXT,ALL,-bh LSEL,LOC,Z,0 LSEL,R,LOC,X,-w_r,w_r ASLL,1 VEXT,ALL,sh-hh ALLSEL CSYS,1 LSEL,LOC,Z,sh-hh LSEL,R,LOC,X,w_sr,w_r ASLL,1 CSYS,0 VEXT,ALL,hh ALLSEL ESIZE,M_h VMESH,ALL圖不見了!再來一次還有這是你的圖.個人覺得很不錯了. 只要按照我所建議的劃分法就應該可以了.YUNQI (云起)谢谢你！太麻烦你了！我要好好想想，把你的建议想明白，该如何做。你的省略号，*_* 要好几天才能看懂噢！&bigbird,你好！我参照你的提示（不知理解的对不对噢？*_*)我现在是这样做的：将这个模型的受载面所在的实体用divide命令分割出来，然后又将外围的实体分割成两个矩形实体，之后将它们用vglue,all命令粘贴起来，这样就可以对他们分别进行分网，均用vsweep,volume,suorce,target命令手动选取源面和目标面，现在已经进行到做100个孔，但是我自己有一个疑问：这样的分割实体之后又粘贴它们（divide-vglue)是否可行？因为自己没有理论解，所以也就不知道所计算出来的结果是否正确？而且，我在分析36个孔时，发现它与解的趋势不相吻合，其他几个都还好。 自这个题目传上来，你就没歇着，一直在帮我下办法，真的太谢谢你了，不知如何来表达这分谢意呢？真的希望自己有一天也能和你一样厉害，向你一样热情的去帮助所有需要帮助的人，让这里成为我们所有学习ansys及其它软件人的家园！我把我的命令流贴上来你看看好么？又让你受累啦！谢谢！他有一点想懒*的裹脚带又臭又长！还是不传好了！宗旨就是先分成三个实体，然后就又将它们全粘贴起来，要不太浪费你的眼睛了！嘻嘻！bigbird.嗯.首先.建模直覺上沒問題.加油喔. 第二. 我不太了解你所謂的趨勢不吻合是什麼意思? 可能原因: 1.網格太粗. 2.邊界條件問題.這可能不是你的原因 先試看看在說吧. 如果還有問題先把有問題的圖PO上來YUNQI (云起)这是我的应力图，但是，你看在粘贴缝上生成了过大的，于是记不相吻合的应力最大的盈利不应该出想在接缝处呀，对么？应该在模型的正中心呀，是不是？你看看这个图。 哎！还是有问题噢！bigbird师兄（可以这样称呼你么？）这是我的网格划分结果，你看看！好么？所谓趋势不稳合，意思是：在我看来，随着受力面上孔的增多，受力面积不是在变小么，那均布力不是在增大么？ 那样的话，在模型的正中心所受的力（等效应力不过现在看来，原来最大的等效应力原来不在模型正中心，哎！那就难说了！可是当时没看出来）但是，我记录下来的从16个孔到100个孔的SEQVSMX的值并不是想我所想的那样应力值在逐渐增大。看来，问题还是出在接缝处呀！怎办？！bigbird.叫我bigbird就好.师兄感覺很老氣.! 針對你這問題 直覺上.不太正確.理由如下 1.這是對偁model. 2.施加負載也是對偁的. 所以.所做出來的也應該是對偁才對. 你是兩邊都施負載嗎?YUNQI (云起)你所说得对称负载是什么意思？我只是在对称面上加了对称的约束，你认为不对么？那对称负载如何加？你看看我这个图.bigbird.對偁負載是.施加均佈力的部分也是對偁的 check 1.你是否是用dsym.作約束的. 2.先以2D模擬一次.看到底是否趨勢吻合.YUNQI (云起)我没有用dysm,但我从来没有见过这个命令，你能说的详细一点么？给我说一个例子好么？说真的，我老师催我。我做的好累噢，今天都快把我说出眼泪来了！bigbird. DSYM commend.用來做對偁約束用.詳細用法看HELP.很簡單的.跟Dcommend一樣. 嗯.!一時之間也想不出有何例子可供說明.你先試試看好了. 哈哈.作research本來就是這樣.別擔心.不要急.一急就更做不出來了YUNQI (云起)我怎么会还有问题了！和我老师讨论了一下，原来约束是应该该已一下，可是还是有一问题：你看这两幅应力图第一副使用最简单的smart3.0网个划分的求解结果，第二副是用divide-glue后用sweep划分的应力图，但是比较这两幅图，为何第二副在约束点有一个最大的应力致点？我总也想不明白呀，你说问题可能出现在事么地方？另一副是：bigbird.這兩張圖趨勢大致相同.可是感覺有singularity. Please check model. 1.兩體積是否有不共面情形. 2.IF 1.沒問題.將mesh再調細一些.jizhenhua00请问bigbird为什么要云lccat，我不用lccat则不能分出你的网格。我在建体网格时，用vext，all 却提示有line concatenate，so cant continue你能否写出后面完整的体网格划分命令bigbird.這是關鍵.好讓你們練習一下.想想why?.為什麼我可以做呢 我只把一個commend給刪除.再試試吧. _ lccat.只要超過四個邊.MAPPED MESH便無法使用.這COMMEND只是將兩個邊.讓ANSYS在MESH時看成一個邊.這樣懂嗎.stree 请问bigbird对称copy怎么做?你先生成一个孔,然后可以生成一百个孔,是怎么实现的YUNQI (云起)先生成一个孔，然后用copy命令就可以生成的，不过需将copy下的kinc定义为0，这样在拷贝时孔上的网格一并copy过去了，你 试一试。bigbird!.終於有些時間看些文章了. 我是用*DO*ENDDO做的. 好處在於.參數設定要做幾個孔.就產生幾個. 只要改一下下就好stree请问yunqi：copy的具体命令是什么，它的菜单方式怎么实现请问bigbird 如果用循环的话，是不是花费的时间比复制多，因为我看资料说，能用复制最好先用复制bigbird.嗯.其實我覺得還好. 我那100個孔也用了幾秒而已(solid model_no mesh). 如果包含mesh. 因為格點數量的關係可能會久一些. 我不太懂能用复制最好先用复制的意思是指哪方面.YUNQI (云起)好久没有来了，谢谢大家一直关注这个问题，对于stree所问的问题用copy或*do*enddo,其实差不多，我都用过，但我自己认为copy更加移懂易些