自由度耦合一例：梁单元和壳单元的连接

Rainyboy

1. !**********************
   
2. !耦合自由度的尝试代码
   
3. !范雨  2009-4-8
   
4. !国际单位制
   
5. !**********************
   
6. FINISH
   
7. /CLEAR
   
8. /TITLE,TEST APDL FOR COUPLING
   
9. 
   
10. /PREP7
    
11. !**********************
    
12. !参数表
    
13. !**********************
    
14. *SET,R_SHELL,0.5!圆柱半径
    
15. *SET,L_SHELL,2!圆柱长度
    
16. *SET,LN_MESH,32!轴向分网
    
17. *SET,RN_MESH,8!周向分网
    
18. *SET,TH,0.01!壳单元厚度
    
19. *SET,MEX,2.06E11!弹性模量
    
20. *SET,MPRXY,0.3!弹性模量
    
21. *SET,MDENS,7850!密度
    
22. *SET,H_BEAM,0.1!梁界面参数：高
    
23. *SET,B_BEAM,0.1!梁界面参数：宽
    
24. *set,IZ_BEAM,(B_BEAM)*(H_BEAM**3)/12!绕Z的惯性矩
    
25. !**********************
    
26. !相关设置
    
27. !**********************
    
28. !设置材料
    
29. MP,EX,1,MEX
    
30. MP,PRXY,1,MPRXY
    
31. MP,DENS,1,MENDS
    
32. !设置壳单元
    
33. ET,1,SHELL63
    
34. R,1,TH
    
35. !设置BEAM单元
    
36. ET,2,BEAM3
    
37. R,2,H_BEAM*B_BEAM,IZ_BEAM,H_BEAM
    
38. !显示SHELL厚度
    
39. /ESHAPE,1
    
40. !以箭头显示压力
    
41. /PSF,PRES,NORM,2,1
    
42. !**********************
    
43. !柱面->几何->分网
    
44. !**********************
    
45. M,1
    
46. TYPE,1
    
47. REAL,1
    
48. !建立柱面
    
49. K,1,0,0,0
    
50. K,2,0,L_SHELL,0
    
51. L,1,2
    
52. CIRCLE,1,R_SHELL,2
    
53. ADRAG,ALL,,,,,,1
    
54. LDELE,1
    
55. !设置分网数目
    
56. ALLSEL,ALL
    
57. LSEL,U,LOC,Y,0
    
58. LSEL,U,LOC,Y,L_SHELL
    
59. LESIZE,ALL,,,LN_MESH
    
60. ALLSEL,ALL
    
61. LSEL,R,LOC,Y,0
    
62. LSEL,A,LOC,Y,L_SHELL
    
63. LESIZE,ALL,,,RN_MESH
    
64. !对圆柱侧面映射分网
    
65. ALLSEL,ALL
    
66. ASEL,U,LOC,Y,0
    
67. ASEL,U,LOC,Y,L_SHELL
    
68. AMESH,ALL
    
69. LESIZE,ALL,,,LN_MESH
    
70. !**********************
    
71. !梁->几何->分网
    
72. !**********************
    
73. !建立辅助直线
    
74. K,100,R_SHELL,0,0
    
75. K,101,R_SHELL,L_SHELL,0
    
76. L,100,101
    
77. K,102,-1*R_SHELL,0,0
    
78. K,103,-1*R_SHELL,L_SHELL,0
    
79. L,102,103
    
80. !设置分网参数
    
81. TYPE,2
    
82. REAL,2
    
83. ALLSEL,ALL
    
84. LSLA,U
    
85. LESIZE,ALL,,,LN_MESH
    
86. !分网
    
87. LMESH,ALL
    
88. !**********************
    
89. !自动耦合自由度
    
90. !**********************
    
91. ALLSEL,ALL
    
92. CPINTF,ALL,0.0001,
    
93. !**********************
    
94. !两端固支：全约束
    
95. !**********************
    
96. NSEL,S,LOC,Y,0
    
97. NSEL,R,LOC,X,-R_SHELL
    
98. D,ALL,ALL,0
    
99. NSEL,S,LOC,Y,L_SHELL
    
100. NSEL,R,LOC,X,-R_SHELL
     
101. D,ALL,ALL,0
     
102. !**********************
     
103. !加载BEAM上的压力载荷
     
104. !**********************
     
105. ESEL,S,TYPE,,2
     
106. ESEL,R,ELEM,,1025,1056
     
107. SFBEAM,ALL,,PRES,-5E5
     
108. ALLSEL,ALL
     
109. !求解
     
110. /SOLU
     
111. SOLVE
     
112. !后处理
     
113. /POST1
     
114. PLNSOL,S,EQV,0,,
     

复制代码

壳体用SHELL53单元建模，上下有加强的梁，用BEAM3建模，它们之间的在共用节点处一一耦合，通过CPINT命令实现。


网格：


位移解（可以看出耦合自由的效果体现了出来）：



应力解：

mkhatebm

恩，命令流很清晰，关键在于耦合自由度那个命令吧，建议用红色标注

Rainyboy

回复 2 # mkhatebm 的帖子

恩，这只是当时做的一个尝试。

zp4655958

耦合自由度  MPC184  约束方程

奋斗中中

Rainyboy 发表于 2011-8-18 11:46

                               
登录/注册后可看大图

回复 2 # mkhatebm 的帖子

恩，这只是当时做的一个尝试。

能不能给个详细的GUI操作啊，新手很多看不懂~~跪求，谢谢了

Rainyboy

奋斗中中 发表于 2012-12-20 16:11

                               
登录/注册后可看大图

能不能给个详细的GUI操作啊，新手很多看不懂~~跪求，谢谢了

迟早要用APDL的啊，想看GUI的话，可以把在帮助里面找相应的命令，会讲到在哪个菜单里面

mxlzhenzhu

Rainyboy 发表于 2012-12-20 22:16

                               
登录/注册后可看大图

迟早要用APDL的啊，想看GUI的话，可以把在帮助里面找相应的命令，会讲到在哪个菜单里面

院长，可不可以给点MPC方面的资料？

比如在一些非协调的时候，如何编辑MPC方程呢？

Rainyboy

mxlzhenzhu 发表于 2013-11-5 13:54

                               
登录/注册后可看大图

院长，可不可以给点MPC方面的资料？

比如在一些非协调的时候，如何编辑MPC方程呢？

这方面的了解有限，ANSYS不是有MPC184单元来处理一些常见的耦合吗？
你说的非协调是指？

非常老妖怪

同求！院长！