Q为问A为回复
--------------------------------------------------------------------------------
Q：SFA的用法：哪位知道其用法，最好有一个例子。
A：SFA的所有用法全在帮助文件里。
SFA, AREA, LKEY, Lab, VALUE, VALUE2
1、area是你要加载荷的面，可以直接指定面的编号，也可以先用asel命令，选取然后用all参数即可。
2、lkey则是同element有关的载荷关键号，不同的element会有不同的意义。
3、lab是用来指定载荷的意义，比如说要加载压力，就用PRES，要加载热流量，就用HFLUX，要加载对流，就用对流热交换，就用CONV。
4、value就很简单了，就是知载荷的大小。针对某些载荷和某些单元而言，一个value不足以表达得给出，另外一个值来。这时就用到了value2。例如在加载热交换载荷CONV时，往往还得给出环境温度（或者是流体温度），这个温度就用value2给出。
SFA,3,CONV,14235,25
这个命令的意思就是在第3个面上加上一个热交换载荷，交换系数是14235，环境温度是25度。
--------------------------------------------------------------------------------
Q：关于正交异性板有谁知道可以用什么单元模拟两个方向刚度不一样的板如何在ansys中实现。
Amatetial props->orthotropic
A：据我了解ansys的材料属性基本上都可以在xyz三个方向分别定义你仔细找一找应该没有问题。
--------------------------------------------------------------------------------
Q：我这里有两个问题向请教一下：
1、对于平面钢架结构来说，比如说有一根杆是斜的，我想对它施加x或是y方向的均部载荷，而不是垂直于杆的均部载荷。
2、对于一个面如何施加剪应力？
A：
1、梁元施加的荷载是针对单元的各个方向的，例如垂直、切向等，而施加整体坐标下的X或Y方向的荷载，只好分解后施加。
2、查看单元的面力规定方向（数字）即可施加可查HELP。
--------------------------------------------------------------------------------
Q：有一个二维模型大部分用原型法(solid modeling)建的模(存在keypoint,line,area等元素)，模型上某个小部分用直接法(direct generation)建的模(即其中没有keypoint,line,area,只有node和element)，请问我怎样才能overlap这两个部分？
A：曾经有人讲过这方面的事情(E文)，但我忘记了。至少有一点是肯定的，即用原形法建立模型划分好单元后，可以接着用直接法继续定义结点和单元是一体的FEA模型不需要任何操作。
A：如果实体模型已经剖分，然后再在这个基础上定义节点，element，当然是没有问题的，不过直接把没有剖分的实体模型和节点，element做布尔运算，恐怕很不妥吧。我个人认为布尔运算不能使用与节点和单元。
--------------------------------------------------------------------------------
Q：怎么对三维实体加绕轴转动的约束，就是只能绕轴转动？
A：如果在柱坐标下让所有节点的UX,UZ自由度为0，是不是就可以了。不过这里的自由度的方向是节点坐标系的方向，所以你得先把节点得坐标系定义为柱坐标系。
A：你也可以试一试ROTX，ROTY，ROTZ这三个自由度，说不定更方便。
A：ansys中还有一种惯性载荷。其中可以定义角速度和角加速度。命令是CGOMGA和DCGOMG你可以仔细琢磨琢磨。
--------------------------------------------------------------------------------
Q：如果一个线既不是水平，也不是垂直，而是与水平线相交成一定的角度，如何加与之平行的均布力？
A：查看单元图形，不同受力的方向定义为不同数字，选择切向就可以了。
--------------------------------------------------------------------------------
Q：怎么画一个圆柱面？
A：
1、建立圆柱体，仅删除体留下的是圆柱面。
2、建立圆弧线，沿一定的路径拖拉也可以建立圆柱面。
--------------------------------------------------------------------------------
Q：请问ansys中的.full文件中如何看刚度矩阵和质量矩阵？
A：you can use the following steps to extract structural matrices:
1. run a modal analysis to get a file.full
2. generate a temp folder
3. copy binlib.dll, rdwrt.bat, rdfull.f, binlib.lib to temp folder (these files can be found in your ansys system)
4. rename rdfull.f to userprog.f
5. at command prompt, run: rdrwrt userprog
6. at command prompt, execute: userprog.exe
7. you can get MASS.MATRIX and STIFFNESS.MATRIX.
A：jobname.full是二进制文件，它是以标准文件头作为第一条记录，以后的若干条记录都有固定的模式，在Interface Guide中对jobname.full的格式做了详细的介绍。而且ANSYS提供了一系列的调用方法能访问各二进制文件，不是很难的事。
建议你仔细研读一下bintst.f，好象还有几个重要的子程序是
binrd.f
binwrt.f
--------------------------------------------------------------------------------
Q：我现在正在学习用ANSYS进行三维流场分析，但在计算完毕后，我不知道怎么看截面流场矢量图，例如对称面上，或者垂直于空气气流流向的截面上。
A：
1、将模型采用正视图或侧视。
2、用SELECT ELEMENT 菜单选box or circle,先用选取再用不选取就可得想看的截面。
--------------------------------------------------------------------------------
Q：如何在屏幕上不显示ANSYS标记
A：/plopts,vers,0
--------------------------------------------------------------------------------
Q：请问如何在后处理图形中加坐标轴，以及注释？
A：加注释可采用utility menu>plotctrl>annotate..然后在屏上点，填写即可。假如是纯粹的注释性，建议你形成bmp文件，然后到画笔中爱干啥干啥！
Q：如何形成bmp文件,我自己也试过用utility menu>plotctrl>annotate>create 2-D annotation,但是写出来的字是倾斜的,该如何处理？
AUtility menu>plotctrls>hardcopy>to file
在图上直接注解时，先点一点，然后拉出直线则字为正。另外，可以设置字体的大小、字体颜色等。
--------------------------------------------------------------------------------
Q：已将预应力筋按LINK8单元划分30等分，各单元的温降也已计算完毕，现在想为这30个单元以读入数据文件方法加温，可是总是只加第一个温度到所有单元。我已经分别用了ARRAY和TABLE 都不行，请指教。
A：你看看这个自己修改搞定如何？
/prep7
et,1,link8
k
k,,100
k,,200,50
l,1,2
l,2,3
lesize,all, , ,10
lmesh,all
以上是建模试验
-------------
!向线上施加温度的初步例题
!从数据文件读入数据并施加
!从线的一端到另一端
!选择线可采用其它方式数组大小可另定
*dim,tem,,20 !定义数组
*vread,tem(1),wd,txt !从文件wd.txt读入
(20f5.1)
csys,0
lineno=2
allsel,all
lsel,s,,,lineno !选择线
nsll,r,1 !选择其中的结点
*get,k1,line,lineno,kp,1 !得到KP1
*get,ncount,node,0,count !得到结点总数
*get,nmax,node,0,num,max !得到最大结点号
*get,nmin,node,0,num,min !得到最小结点号
nnum=node(kx(k1),ky(k1),kz(k1)) !得到一端的结点号
*do,loopno,1,ncount,1 !循环结点总数
bf,nnum,temp,tem(loopno) !施加温度值
*if,nnum,eq,nmin,then !得到紧邻下一结点号
nnum=ndnext(nmin)+1
*elseif,nnum,eq,nmax,then
nnum=ndnext(nmin)
*else
nnum=ndnext(nnum)
24
*endif
*enddo !循环结束标志
allsel,all
bflis,all
!----------------------
以下是wd.txt文件内容另存为文件
10.0 20.0 30.0 40.0 50.0 60.0 70.0 71.0 72.0 73.0 74.0 75.0
10.0 20.0 30.0 40.0 50.0 60.0 70.0 71.0 72.0 73.0 74.0 75.0
--------------------------------------------------------------------------------Q如何加密求得的等值线,我现在看到的太稀疏了
A设为等值线方式/device,vector,1
定义标记疏密/clabel,1,no1
设根数:/cont,1,no2
其中no1,no2是自己定义的试试疏密满意为止.
--------------------------------------------------------------------------------QA!/nopr
!向线上施加荷载的宏
!需要四个参数
!lineno:线号
!p1:一端的荷载值
!p2:另端的荷载值
!prekey:荷载方向号
!------------------
!从对话框输入上述参数
multipro,'start',4 !对话框输入数值
*cset,1,3,lineno,'the line number'
*cset,4,6,p1,'the pressure at one end'
*cset,7,9,p2,'the pressure at other end'
*cset,10,12,prekey,'the pressure key',1
multipro,'end'
csys
!确定选择集(该线及其单元)并得到线长
lsel,s,,,lineno !选择线
esll,r !选择此线中的单元
*get,linecd,line,lineno,leng !得到线长
*get,qdh,line,lineno,kp,1 !得到KP1
!处理荷载(对选择集而言)
*get,ecount,elem,0,count !得到单元总数
*get,enum,elem,0,num,min !得到单元的最小序号
*do,loopno,1,ECOUNT,1 !循环单元总数
*get,i,elem,ENUM,node,1 !得到当前单元的i结点号
*get,j,elem,ENUM,node,2 !得到当前单元的j结点号
dx=nx-kx(qdh) !求I结点到KP1的距离
dy=ny-ky(qdh)
25
dz=nz-kz(qdh)
inloc=sqrt(dx*dx+dy*dy+dz*dz)
dx=nx(j)-kx(qdh) !求J结点到KP1的距离
dy=ny(j)-ky(qdh)
dz=nz(j)-kz(qdh)
jnloc=sqrt(dx*dx+dy*dy+dz*dz)
ipress=(p2-p1)/linecd*inloc+p1 !计算I结点荷载
jpress=(p2-p1)/linecd*jnloc+p1 !计算j结点荷载
sfbeam,enum,prekey,pres,ipress,jpress !施加该单元荷载
*get,ENUM,elem,ENUM,nxth !从选择集中得到下一个单元号
*enddo !循环结束标志
ALLSEL,ALL
!宏结束
--------------------------------------------------------------------------------
Q:建立实体模型后，计算完毕，要配筋时有无方法直接将应力转换为内力，计算配筋和裂缝？我知道梁单元有输出，那实体呢？
A：梁板才有M、Q。本来就是简化的概念，在ANsys中用beam、shell单元都可直接输出内力，用所谓"sequence number"，solid没有。
--------------------------------------------------------------------------------
Q：ANSYS里面好像只有STRESS 结果，能不能察看内力值？
A：是可以查看诸如轴力、弯矩等结果的，你要详细看一下单元的输出数据（在帮助中）。
A：当然可以了，每种单元都有关于输出数据的定义。例如，beam4想知道y方向的弯矩，查到与之对应的是SMISC2、SMISC8；在general postproc－element table－define table－add。在弹出的对话框中选择by sequence num，定义后在plot或者list element table中查看。
--------------------------------------------------------------------------------
Q：不知道杆单元有没有内力值？还有，我要分析一个筏板基础，用弹簧单元模拟土体，不知道弹簧单元如何加到模型当中，是先离散筏板，还是先加弹簧单元？如果先离散筏板，如果单元不规则，我如何计入弹簧单元的弹黄刚度（弹黄刚度 KN/m=弹黄刚度KN/m3 * 单元面积m2）？
A：杆单元有内力值，查help看你所用的单元的element table分析筏板基础，可选link10单元模拟土弹簧，设定单元只能受压。
A：记住有一点，查到的力和弯矩对应的是单元的局部坐标系。
Q：请问
1、ansys中的Drcker-Prager屈服准则用在土的弹塑性模拟中精度怎么样？
2、有没有人采用广西大学秦荣的QR及样条函数法进行过桩土共同作用的分析？
3、ansys有没有无限元单元和岩石节理元？
A：D-P模型虽然没有帽盖，因为与Mohr-C模型的密切关系。还是可以模拟土的很多特点。但是国内用的人很少Duncan-Chueng，Modified Cambridged 等是在国内流行的。对于实际工程，除了合适的模型外，参数的获得是至关重要的。参数的差异可能大于由于模型带来的差异。你的意思，好象是模型分析，自己定参数，对吗？我认为，限于定性研究D-P模型应该可以的。但是精度，我觉得这里没有精度概念，因为你是理想化的研究。
--------------------------------------------------------------------------------
