耦合分析在ansys中的实现方法与技巧一

耦合作为ansys模拟现实世界非常重要的一种手段，现在越来越受到重视。下面以一个实例来说明。
问题描述：如图1所示，三根杆件几何尺寸以及材料属性均相同，在A点受一集中力作用，试分析三根杆件的轴力。
分析与假设：假设杆件处于线弹性阶段；p=1，l=a=1；由于本问题结果与弹性模量、泊松比无关，所以可设置为任意值；耦合方程：（uy）2+（uy）6 - 2*（uy）4=0。
命令流如下：
/clear,start ！清空内存
/filname,couple anslysis ！定义文件名
/title,bar ！定义标题
fini
/prep7 ！进入前处理
et,1,1 ！定义单元link1
r,1,1 ！定义实常数
mp,ex,1,2e11 ！定义弹性模量
mp,prxy,1,.3 ！定义泊松比
n,1 ！创建节点
n,2,,-1
n,3,1
n,4,1,-1
n,5,2
n,6,2,-1
e,1,2 ！创建单元
egen,3,2,1
d,1,all,,,5,2 ！施加约束
f,2,fy,-1 ！施加集中力
ce,1,,2,uy,1,6,uy,1,4,uy,-2 ！定义耦合方程
allsel,all ！选择所有实体
fini
/solu ！进入求解模块
solve ！执行求解
fini
/post1 ！进入后处理模块
etable,str,ls,1 ！定义单元表
pretable,str ！列表显示单元表
列表显示结果如图2所示。
从图2可以看出一杆轴力为0.83333=5/6，二杆轴力为0.33333=1/3，三杆轴力为-0.16667=-1/6；而理论解为N1=5/6，N2=1/3，N3=-1/6；数值解与理论解吻合的相当好。