选择球坐标系的奇异点位置 cscir是什么意思？

xiaowei331

/filname,schwedler
/title,analysis of spherical reticulated-shell
/prep7
et,1,beam4
r,1,1382,2093500,2093500,114,114,,
rmore,4187000,,,,,
mp,dens,1,7.85e-6
mp,ex,1,2.06e5
mp,nuxy,1,0.3
local,11,2,0,0,0
cscir,11,1
n,1,20,0,30
n,10,20,180,30
n,18,20,340,30
fill,1,10
fill,10,18
ngen,6,18,1,18,1,0,0,10
n,109,20,0,90
*do,i,1,91,18
  e,i,i+1
  egen,17,1,i,i,1
e,i+17,i
*enddo
e,1,19
egen,5,18,109,109,1
egen,18,1,109,113,1
*do,i,91,108,1
e,i,109
*enddo
*do,i,1,73,18
e,i,i+19
*enddo
egen,17,1,217,221,1
*do,i,18,90,18
e,i,i+1
*enddo
csys,0
nsel,s,loc,z,9.9,10.1
d,all,,,,,,ux,uy,uz,,

大侠帮帮忙  我想问一下这里面的CSCIR，11，1是什么意思   
它的解释是“选择球坐标系的奇异点位置”    可是我不明白的是这儿的奇异点是什么意思      以及这个CSCIR究竟有什么用   谢谢了

[ 本帖最后由 无水1324 于 2007-11-5 21:06 编辑 ]

luobinhan001

正巧我也碰到了这个问题，google了一下，找到了详细说明。
出自http://pera.e-works.net.cn/document/200708/article1756_2.htm


在圆柱坐标系(例如SLKCN=1)中给定梯度时，还应记住以下几点。第一，SLZER以度表示，SLOPE以载荷大小/度表示。第二，应遵循下列两个规则：
规则1:设置CSCIR（为了控制坐标系奇异点位置) 使待加载的表面不通过坐标系奇异点。
规则2:选择SLZER，使之与CSCIR设置一致。即：如果奇异点在180°处[CSCIR,KCN,0]，SLZER应在 ±180°之间；如果奇点在0°处[CSCIR,KCN,1]，SLZER应在0° -360° 之间。
下例说明为什么要给出这些规则。对图2-10所示位于局部柱坐标系11的半圆壳，将对半圆壳施加一个作用外部的楔形压力，压力从-90°位置的400逐渐变化到90°位置的580。缺省情况下，奇异点位于柱坐标系中的180°，因此，壳的坐标范围从-90° to +90°。下列命令将用于施加所需的压力载荷。
SFGRAD,PRES,11,Y,-90,1 ! 使压力沿C.S. 11的θ方向线性变化
! 指定压力作用于-90度，斜率为1个单位/度
SF,ALL,PRES,400! 在所有被选择的节点压力:
!在-90度为400., 在0度为490., 在+90度为580.
在-90°，压力值为400 (指定), 以1个单位/度的斜率增加，在0度位置增加到490，在+90度位置增加到580。
对于SLZER，可能会诱导用户使用270°而不是-90°。
SFGRAD,PRES,11,Y,270,1 ! 使压力沿C.S. 11的θ方向线性变化
! 指定压力作用于270度，斜率为1个单位/度
SF,ALL,PRES,400! 在所有被选择的节点压力:
!在-90度为400., 在0度为490., 在+90度为580.

                               
登录/注册后可看大图

图2-10圆柱壳上线性增加的载荷

然而，如图2-11左边所示，这可能导致所施加的逐渐变化载荷与要求的载荷值不同。这是因为奇点仍位于180°（坐标范围从-90° 到 +90°，而SLZER不在-180°到+180°范围内）。结果，程序将使用270°位置处的载荷400，和1个单位/度的斜率计算得到：施加于+90°位置处的载荷为220，施加于90°位置处的载荷为130，施加于-90°位置处的载荷为40。依照第二个规则，则可避免这种情况的发生。即当奇点在180°位置时，选择SLZER在±180°之间；当奇点在0°位置时，选择SLZER在0°到360°之间；

                               
登录/注册后可看大图

图2-11违背规则2（左）和规则1（右）

假设将奇点位置改变到0°，因此满足第二条规则（270°在0°到360°之间）。但壳的上半部分，节点坐标的范围在0°到+90°之间，而壳的下半部分，节点坐标的范围在270°到360°之间，待加载的表面过奇点，违背规则1。
CSCIR,11,1! 将奇点改变到零度
SFGRAD,PRES,11,Y,270,1 ! C.S. 11.的θ方向施加线性变化的压力
压力作用于270度位置，斜率为1单位/度，SF,ALL,PRES,400! 所有被选择的节点压力：
!在270度位置为400，在360度位置为490!在0度位置为130
程序将使用270度位置的载荷400和1单位/度的斜率计算得到：施加于270度位置的载荷值为400，360度位置的载荷为490，90度位置的载荷为220，0度位置的载荷为130。违背规则1在逐渐变化的载荷上将产生1个奇点，如图2-11右所示。由于节点离散化，实际施加的载荷在奇点将不会发生如图所示的剧烈变化。反而，在所示的情况下，0度处节点的载荷为130，下一个沿顺时针方向的节点（比如在358度位置）的载荷为488。
注：SFGRAD指定对所有随后的载荷施加命令都起作用。要去除该指定，仅需发一个不带变元的SFGRAD命令。而且在读取载荷步文件时，如果SFGRAD指定起作用，那么程序在读取文件前将删除该指定。
大变形作用会显著改变节点位置。基于节点位置的SFGRAD斜率和载荷值计算不会根据这些变化作相应更新。如果需要这种功能，使用表面3加载的SURF19或SURF153命令，或使用表面4加载的SURF22或SURF154命令