[求助]单元耦合问题

yinyejun · 发表于 2006-6-6 20:16

请教楼主，exam5中是从下面的部分开始耦合么？我不大懂命令流，有些问题想请教 /view,1,1,1,1 /ang,1 eplot asel,s,loc,x,0 以下三条命令什么意思呀？ da,all,all allsel,all cpintf,ux 以下三条命令什么意思呀？ cpintf,uy cpintf,uz ce,1,0,1226,rotz,-abs(ny(197)-ny(742)),197,ux,1,742,ux,-1 1226，197，742这个节点号是怎么知道的？ ce,2,0,1276,rotz,-abs(ny(201)-ny(971)),201,ux,1,971,ux,-1 ce,3,0,1277,rotz,-abs(ny(205)-ny(969)),205,ux,1,969,ux,-1 ce,4,0,1278,rotz,-abs(ny(209)-ny(967)),209,ux,1,967,ux,-1 ce,5,0,1279,rotz,-abs(ny(213)-ny(965)),213,ux,1,965,ux,-1 ce,6,0,1170,rotz,-abs(ny(194)-ny(701)),194,ux,1,701,ux,-1 ce,7,0,1226,rotx,-abs(ny(197)-ny(742)),197,uz,1,742,uz,-1 ce,8,0,1276,rotx,-abs(ny(201)-ny(971)),201,uz,1,971,uz,-1 ce,9,0,1277,rotx,-abs(ny(205)-ny(969)),205,uz,1,969,uz,-1 ce,10,0,1278,rotx,-abs(ny(209)-ny(967)),209,uz,1,967,uz,-1 ce,11,0,1279,rotx,-abs(ny(213)-ny(965)),213,uz,1,965,uz,-1 ce,12,0,1170,rotx,-abs(ny(194)-ny(701)),194,uz,1,701,uz,-1 nsel,s,loc,x,120 以下几条什么意思呀？ *get,nhzs,node,,count f,all,fy,-3/nhzs allsel,all [em06]

djaity · 发表于 2006-6-6 20:55

临近区耦合(adjacent rejgion) 先选择网格划分较好的模型上的节点,在选则另一模型上的单元,然后选择菜单弹出此对话框既可实现耦合

sysh320 · 发表于 2006-6-7 21:38

楼主的例子很有用的！

yakexi · 发表于 2006-6-7 21:50

<DIV class=quote>以下是引用usefully在2006-6-6 13:22:50的发言： 
梁壳体连接情况的讨论 1 按“杆梁壳体”的顺序，只要后一种单元的自由度完全包容了前一种单元的自由度，则有公用结点即可，不需要约束方程。例如： 杆与梁、壳、体有公用结点即可，不需要写约束方程； 梁与壳有公用结点即可，不需要写约束方程； 梁与体则要同位置的不同结点，需要耦合自由度和约束方程； 壳与体则要同位置的不同结点，需要耦合自由度和约束方程； 2 壳梁自由度数目相同，自由度也相同，尽管壳的rotz是虚的自由度，也不妨碍二者之间的关系，这有点类同于梁与杆的关系。 3 尽管可以采用耦合自由度和约束方程，但建议尽量不同时采用多种单元于一个结构中，除非你对结果的正确性有十足的把握。 4 当然，采用约束方程可能存在应力集中点，不必在意此点的应力。 5 我自认为是正确的，希望各位大侠批评指正。 6 为说明上述说法的正确性，这里提供有5个小例子。例1是全“壳单元”，例2是“梁壳单元”；例3是全“体单元”，例4是“体梁单元”，例5是“体壳单元”。 7 运行于ansys6.1下，三月雨提供。 !********************************************** !梁壳的耦合问题小算例 !采用壳单元时exam1 /prep7 et,1,shell63 mp,ex,1,3e5 mp,prxy,1,0.0 r,1,1.0 wprota,0,90 blc4,,,120,10 aesize,all,2 mshape,0,2d mshkey,1 amesh,all /view,1,1,1,1 /ang,1 eplot nsel,s,loc,x,0 d,all,all nsel,s,loc,x,120 *get,nhzs,node,,count f,all,fy,-3/nhzs allsel,all /solu solve !********************************************** !梁壳的耦合问题小算例 !采用梁壳单元时exam2 /prep7 et,1,shell63 et,2,beam4 mp,ex,1,3e5 mp,prxy,1,0.0 r,1,1.0 r,2,10.0,10/12.0,1000/12.0,10.0,1.0 wprota,0,90 blc4,,,60,10 wpoff,,5 wprota,,90 asbw,all k,100,120,0,5 ksel,s,loc,x,60,120 ksel,r,loc,z,5 *get,kp1,kp,,num,min kp2=kpnext(kp1) l,kp1,kp2 ksel,all wpcsys,-1 lsel,s,loc,x,60,120 latt,1,2,2 lesize,all,,,10 lmesh,all asel,all aatt,1,1,1 aesize,all,2 mshape,0,2d mshkey,1 amesh,all allsel,all /view,1,1,1,1 /ang,1 eplot finish !------------------------ /solu nsel,s,loc,x,0 d,all,all nsel,s,loc,x,120 f,all,fy,-3.0 allsel,all solve /post1 etable,zl1,smisc,1 etable,zl2,smisc,7 etable,jly1,smisc,2 etable,jly2,smisc,8 etable,jlz1,smisc,3 etable,jlz2,smisc,9 etable,mx1,smisc,4 etable,mx2,smisc,10 etable,my1,smisc,5 etable,my2,smisc,11 etable,mz1,smisc,6 etable,mz2,smisc,12 !********************************************** !梁壳体的耦合问题小算例 !全部采用体单元时exam3 /prep7 et,1,solid95 mp,ex,1,3e5 mp,prxy,1,0.0 r,1 blc4,,,20,7,10 blc4,20,3,100,1,10 vglue,all wpoff,0,3 wprota,0,90 vsbw,all wpoff,0,0,-1 vsbw,all wpstyle /view,1,1,1,1 /ang,1 vplot esize,1 mshape,0,2d mshkey,1 vmesh,all finish /solu asel,s,loc,x,0 da,all,all allsel,all nsel,s,loc,x,119.6,120 nsel,r,loc,y,4 *get,nhzs,node,,count f,all,fy,-3/nhzs allsel,all solve !********************************************** !梁壳体的耦合问题小算例 !采用体单元和梁单元时exam4 /prep7 et,1,solid95 et,2,beam4 mp,ex,1,3e5 mp,prxy,1,0.0 r,1 r,2,10.0,10/12.0,1000/12.0,10.0,1.0 blc4,,,20,7,10 wpoff,0,3.5 wprota,0,90 vsbw,all wpoff,0,5 wprota,0,90 vsbw,all wpcsys,-1 k,100,20,3.5,5 k,101,120,3.5,5 l,100,101 lsel,s,loc,x,21,130 latt,1,2,2 lesize,all,,,10 lmesh,all vsel,all vatt,1,1,1 esize,1 mshape,0,2d mshkey,1 vmesh,all allsel,all /view,1,1,1,1 /ang,1 eplot finish !------------------------ /solu asel,s,loc,x,0 da,all,all allsel,all fk,101,fy,-3.0 cp,1,ux,1,21 cp,2,uy,1,21 cp,3,uz,1,21 ce,1,0,626,ux,1,2328,ux,-1,1,roty,-abs(nz(626)-nz(2328)) ce,2,0,67,ux,1,4283,ux,-1,1,rotz,-abs(ny(67)-ny(4283)) ce,3,0,67,uz,1,4283,uz,-1,1,rotx,-abs(ny(67)-ny(4283)) allsel,all solve finish !------------------------ /post1 etable,zl1,smisc,1 etable,zl2,smisc,7 etable,jly1,smisc,2 etable,jly2,smisc,8 etable,jlz1,smisc,3 etable,jlz2,smisc,9 etable,mx1,smisc,4 etable,mx2,smisc,10 etable,my1,smisc,5 etable,my2,smisc,11 etable,mz1,smisc,6 etable,mz2,smisc,12 !********************************************** !梁壳体的耦合问题小算例 !采用体壳单元时exam5 /prep7 et,1,solid95 et,2,shell63 mp,ex,1,3e5 mp,prxy,1,0.0 r,1 r,2,1.0 blc4,,,20,7,10 wpoff,0,3.5 wprota,0,90 vsbw,all wpoff,20 blc4,,,100,10 wpcsys,-1 vsel,all vatt,1,1,1 esize,2 mshape,0,2d mshkey,1 vmesh,all asel,s,loc,x,21,120 aatt,1,2,2 aesize,all,2 mshape,0,2d mshkey,1 amesh,all allsel,all /view,1,1,1,1 /ang,1 eplot asel,s,loc,x,0 da,all,all allsel,all cpintf,ux cpintf,uy cpintf,uz ce,1,0,1226,rotz,-abs(ny(197)-ny(742)),197,ux,1,742,ux,-1 ce,2,0,1276,rotz,-abs(ny(201)-ny(971)),201,ux,1,971,ux,-1 ce,3,0,1277,rotz,-abs(ny(205)-ny(969)),205,ux,1,969,ux,-1 ce,4,0,1278,rotz,-abs(ny(209)-ny(967)),209,ux,1,967,ux,-1 ce,5,0,1279,rotz,-abs(ny(213)-ny(965)),213,ux,1,965,ux,-1 ce,6,0,1170,rotz,-abs(ny(194)-ny(701)),194,ux,1,701,ux,-1 ce,7,0,1226,rotx,-abs(ny(197)-ny(742)),197,uz,1,742,uz,-1 ce,8,0,1276,rotx,-abs(ny(201)-ny(971)),201,uz,1,971,uz,-1 ce,9,0,1277,rotx,-abs(ny(205)-ny(969)),205,uz,1,969,uz,-1 ce,10,0,1278,rotx,-abs(ny(209)-ny(967)),209,uz,1,967,uz,-1 ce,11,0,1279,rotx,-abs(ny(213)-ny(965)),213,uz,1,965,uz,-1 ce,12,0,1170,rotx,-abs(ny(194)-ny(701)),194,uz,1,701,uz,-1 nsel,s,loc,x,120 *get,nhzs,node,,count f,all,fy,-3/nhzs allsel,all /solu solve</DIV>
呵呵，首先楼主总结的不错，不过我还是觉得具体问题具体分析比较好。 比如一个一个板下面有4个柱子可以设为杆或梁，这样的分析老兄作过没，单约束结点是可以但是计算结果是错误的，应力集中太大，如果复杂结构都这么考虑，相关效应就来了。

hhh3836 · 发表于 2006-6-8 18:36

不耦合直接算，得出的结果可信度有多少呢？ 我遇到过类似的问题，当时就没耦合直接算的，结果我除了感觉变形 有些大外，其它的好像没太大的疑点。 大家对这个问题都是怎么处理的啊？[em32][em32]

多情清秋 · 发表于 2006-6-11 10:54

usefully加威望2点，zglecsi hhh3836 yinyejun AaronSpark linqus 渊源1983 yakexi djaity sysh320加威望1点
多情清秋 06.6.11

usefully · 发表于 2006-6-12 20:18

ANSYS的每个位移约束方程可以包含任意个自由度项。
在用CE命令定义位移约束方程时，每一行输入三个自由度项，其余的自由度项按续行处理加入该方程。续行时，同样要用CE命令，只是位移约束方程号和常数项空白即可。
例如：
ce,1,0,1,ux,1,1,uy,1,2,ux,-1 ce,,,2,uy,-1,3,roty,0.5,4,ux,0.1 ce,,,,5,rotx,0.1
它表示位移约束方程1为:
1*Ux(N1)+1*Uy(N1)-1*Ux(N2)-1*Uy(N2)+0.5*Roty(N3)+0.1*Ux(N4)+0.1*Rotx(N5)=0

yinyejun · 发表于 2006-6-12 20:23

多谢楼主！

yinyejun · 发表于 2006-6-12 20:29

请教楼主，请问你的关于壳单元和体单元的耦合问题具体是怎么解决的？我的模型还没有耦合上呢！愁死人了！ 谢谢了！

havvn · 发表于 2006-6-13 19:46

《有限元法实用教程》刘桂容（新加坡）这本书上有一章是论述单元耦合的但是实现方法是用abaqus 我也是ansys入门级的人，水平太低！呵呵期待高手！

账号		自动登录	找回密码
密码			我要加入

[求助]单元耦合问题

回复：（usefully）回复：（usefully）[求助]单元耦合...

回复：（usefully）[求助]单元耦合问题

回复：（yinyejun）请教楼主，exam5中是从下面的部分开...

多谢楼主

请教楼主