施加非法向的面力的方法(from bbs.tongji.edu.com)

qzhang22

Q: 如何施加非法向的面力，比如作用在斜杆上的竖向均布荷载。
A：在任意面施加任意方向任意变化的压力
在某些特殊的应用场合，可能需要在结构件的某个面上施加某个坐标方向的随坐标位置变
化的压力载荷，当然，这在一定程度上可以通过ANSYS表面效应单元实现。如果利用ANSYS
的参数化设计语言，也可以非常完美地实现此功能，下面通过一个小例子描述此方法。

!!!在执行如下加载命令之前,请务必用选择命令asel将需要加载的几何面选择出来
!!!
finish
/prep7
et,500,shell63
press=100e6
amesh,all
esla,s
nsla,s,1
! 如果载荷的反向是一个特殊坐标系的方向,可在此建立局部坐标系,并将
! 所有节点坐标系旋转到局部坐标系下.
*get,enmax,elem,,num,max
dofsel,s,fx,fy,fz
fcum,add !!!将力的施加方式设置为"累加",而不是缺省的"替代"
*do,i,1,enmax
*if,esel,eq,1,then
*get,ae,elem,i,area !此命令用单元真实面积，如用投影面积，请用下几条命令
! *get,ae,elem,i,aproj,x !此命令用单元X投影面积，如用真实面积，请用上一条命令

! *get,ae,elem,i,aproj,y !此命令用单元Y投影面积
! *get,ae,elem,i,aproj,z !此命令用单元Z投影面积
xe=centrx !单元中心X坐标(用于求解压力值)
ye=centry !单元中心Y坐标(用于求解压力值)
ze=centrz !单元中心Z坐标(用于求解压力值)
! 下面输入压力随坐标变化的公式,本例的压力随X和Y坐标线性变化.
p_e=(xe-10)*press+(ye-5)*press
f_tot=p_e*ae
esel,s,elem,,i
nsle,s,corner
*get,nn,node,,count
f_n=f_tot/nn
*do,j,1,nn
f,nelem(i,j),fx,f_n !压力的作用方向为X方向
! f,nelem(i,j),fy,f_n !压力的作用方向为Y方向
! f,nelem(i,j),fz,f_n !压力的作用方向为Z方向
*enddo
*endif
esla,s
*enddo
aclear,all
fcum,repl !!!将力的施加方式还原为缺省的"替代"
dofsel,all
allsel


--


;32m※ 来源:．同舟共济站 http://bbs.tongji.edu.cn [FROM: 192.168.16.146]m


  





[本篇全文] [回复本文] [本篇作者: zhenghzh ] [本篇人气: 7]
发信人: zhenghzh (小郑), 信区: FEA
标  题: Re: 施加非法向的面力的方法
发信站: 同舟共济站 (Wed May 15 17:14:44 2002)

其中有一段应为
*do,i,1,enmax
*if,esel(i),eq,1,then！！！
*get,ae,elem,i,area !此命令用单元真实面积，如用投影面积，请用下几条命令
! *get,ae,elem,i,aproj,x !此命令用单元X投影面积，如用真实面积，请用上一条命令


! *get,ae,elem,i,aproj,y !此命令用单元Y投影面积
! *get,ae,elem,i,aproj,z !此命令用单元Z投影面积
xe=centrx(i) !单元中心X坐标(用于求解压力值)！！！
ye=centry(i) !单元中心Y坐标(用于求解压力值)！！！
ze=centrz(i) !单元中心Z坐标(用于求解压力值)！！！
--------
原文中少了(i)，我去有限元联盟看了一下，发现它用一个电灯泡来放在(i)的位置上，

一共有四处。

这个做法我还没有试过，APDL的功能确实很强大。
--
海阔凭鱼跃，
天高任鸟飞。
;32m※ 来源:．同舟共济站 http://bbs.tongji.edu.cn [FROM: 192.168.46.151]m


  





[本篇全文] [回复本文] [本篇作者: zhenghzh ] [本篇人气: 10]
发信人: zhenghzh (小郑), 信区: FEA
标  题: Re: 施加非法向的面力的方法
发信站: 同舟共济站 (Wed May 15 21:12:21 2002)

刚刚试着用这个方法施加面上的分布荷载，效果还可以。
这个模块基本上不用怎么改，就可以放进去用。
只是还有点疑惑，一开始用shell63对所选的面划分网格，
是不是就是一种辅助手段？我还是第一次碰到。

还有就是他是把力作用在节点上面，由于我用的是对称结构，
荷载也是轴对称，且只取了其中的一片，这样一来，
在对称面上的节点上的力就只有其他非对称面上的节点上的一半。
不知道ansys在处理对称的时候，会不会把对称面上的力也对称一下？
而且这个也比较难检验，现在还不敢确定对不对。

--
海阔凭鱼跃，
天高任鸟飞。
;31m※ 来源:．同舟共济站 http://bbs.tongji.edu.cn [FROM: 192.168.46.151]m


  





[本篇全文] [回复本文] [本篇作者: gxn ] [本篇人气: 2]
发信人: gxn (老农), 信区: FEA
标  题: Re: 施加非法向的面力的方法
发信站: 同舟共济站 (2002年05月22日02:53:47 星期三), 站内信件

xe=centrx(i) !单元中心X坐标(用于求解压力值)！！！
这几句开始ansys不认识，我后来用了*get,xe,centrx,i来代替了


--
我是农民我怕谁
m;36m※ 来源:·同舟共济站 bbs.tongji.edu.cn·[FROM: 211.156.24.13]m

dqtznn

谢谢