各位前辈,本人接触ansys不久,最近在学习用ansys分析:轴对称模型施加非轴对称载荷的问题。 我做了一个例题,加载范围在-6°到6°之间。 遇到了许多问题,希望前辈给我一些指导。
1、将载荷用傅里叶级数展开后,每一项作为一个载荷步,如何使用命令流来加载函数,如何定义加载范围?
(因傅里叶展开后有n项,其中每一项为:(2e7)/(n*{PI})*sin(n*6)*cos(n*{Y}),其中{Y}用来表示θ。) (附上的命令流只加了六项,这不是好办法,所以希望用到循环语句,希望前辈给予指导。) 2、傅里叶级数展开第一项为常数项,加载范围是不是0到360°? 3、傅里叶级数其他关于θ的项加载范围是不是-6到6°?
望前辈不吝赐教。
下面附上命令流
/prep7 !创建点
et,1,plane25
mp,ex,1,2.1e11
mp,nuxy,1,0.3
k,1,100 !关键点
k,2,420,
k,3,420,100
k,4,400,100
k,5,400,15
k,6,120,15
k,7,120,100
k,8,100,100
a,1,2,3,4,5,6,7,8 !创建面
wpoffs,,15 !偏移工作平面 Y轴方向
wprota,,90 !旋转工作平面 绕Y轴90
asbw,1 !用工作平面 分割1面
wpoffs,120
wprota,,,90
asbw,4
wpoffs,,,280
asbw,5
LESIZE,7,,,4 !布点网格尺寸大小分布 7号线分成四段
LESIZE,11,,,4
LESIZE,2,,,4
LESIZE,3,,,4
LESIZE,10,,,4
LESIZE,17,,,4
LESIZE,4,,,17
LESIZE,9,,,17
LESIZE,6,,,17
LESIZE,12,,,17
LESIZE,13,,,3
LESIZE,8,,,3
LESIZE,1,,,3
LESIZE,14,,,3
LESIZE,5,,,56
LESIZE,16,,,56
ASEL,ALL !全选
AATT,1,,1 !给所选的面设置划分网格单元属性。材料号,实常数值,单元号
MSHAPE,0,2D !分割单元 四边形 2D
AMESH,ALL !根据面上指定的网格生成一个二维映射网格
/SOL
DL,12,,ALL !模型内径全约束
DL,13,,ALL
wpoffs,,,-280 !调整工作平面与坐标面重合
wprota,,,-90
wpoffs,-120
wprota,,-90
wpoffs,,-15
EPLOTE !显示所有选择的单元网格
CSYS,5 !激活一个预先已定义的坐标系统 Y轴为轴的柱坐标系
allsel,all
MODE,0,1
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNCCSYS
*SET,_FNCNAME,'f0'
*SET,_FNCCSYS,5
! /INPUT,f0.func,,,1
*DIM,%_FNCNAME%,TABLE,6,4,1,,,,%_FNCCSYS%
!
! Begin of equation: (1e7)/30
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 1e7, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 30, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 1, -3, 0, 1, -1, 4, -2
*SET,%_FNCNAME%(0,4,1), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: (1e7)/30
!-->
FLST,2,2,4,ORDE,2
FITEM,2,9
FITEM,2,14
/GO
!*
!*
SFL,9,PRES, %F0%
SFL,14,PRES, %F0%
LSWRITE,1,
allsel,all
MODE,1,1
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNCCSYS
*SET,_FNCNAME,'f1'
*SET,_FNCCSYS,5
! /INPUT,f1.func,,,1
*DIM,%_FNCNAME%,TABLE,6,8,2,,,,%_FNCCSYS%
!
! Begin of equation: {Y}
*SET,%_FNCNAME%(0,0,1), -6, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, 99, 0, 1, 3, 0, 0
*SET,%_FNCNAME%(0,2,1), 0
*SET,%_FNCNAME%(0,3,1), 1
*SET,%_FNCNAME%(0,4,1), 0
*SET,%_FNCNAME%(0,5,1), 0
*SET,%_FNCNAME%(0,6,1), 0
*SET,%_FNCNAME%(0,7,1), 0
*SET,%_FNCNAME%(0,8,1), 0
! End of equation: {Y}
!
! Begin of equation: (2e7)/{PI}*sin(6)*cos({Y})
*SET,%_FNCNAME%(0,0,2), 6, -999
*SET,%_FNCNAME%(2,0,2), 0.0
*SET,%_FNCNAME%(3,0,2), 0.0
*SET,%_FNCNAME%(4,0,2), 0.0
*SET,%_FNCNAME%(5,0,2), 0.0
*SET,%_FNCNAME%(6,0,2), 0.0
*SET,%_FNCNAME%(0,1,2), 1.0, -1, 0, 2e7, 0, 0, 0
*SET,%_FNCNAME%(0,2,2), 0.0, -2, 0, 3.14159265358979310, 0, 0, -1
*SET,%_FNCNAME%(0,3,2), 1, -3, 0, 1, -1, 4, -2
*SET,%_FNCNAME%(0,4,2), 0.0, -1, 9, 6, 0, 0, 0
*SET,%_FNCNAME%(0,5,2), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,6,2), 0.0, -1, 10, 1, 3, 0, 0
*SET,%_FNCNAME%(0,7,2), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,8,2), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: (2e7)/{PI}*sin(6)*cos({Y})
!-->
FLST,2,2,4,ORDE,2
FITEM,2,9
FITEM,2,14
/GO
!*
!*
SFL,9,PRES, %F1%
SFL,14,PRES, %F1%
LSWRITE,2,
allsel,all
MODE,2,1
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNCCSYS
*SET,_FNCNAME,'f2'
*SET,_FNCCSYS,5
! /INPUT,f2.func,,,1
*DIM,%_FNCNAME%,TABLE,6,15,2,,,,%_FNCCSYS%
!
! Begin of equation: {Y}
*SET,%_FNCNAME%(0,0,1), -6, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, 99, 0, 1, 3, 0, 0
*SET,%_FNCNAME%(0,2,1), 0
*SET,%_FNCNAME%(0,3,1), 1
*SET,%_FNCNAME%(0,4,1), 0
*SET,%_FNCNAME%(0,5,1), 0
*SET,%_FNCNAME%(0,6,1), 0
*SET,%_FNCNAME%(0,7,1), 0
*SET,%_FNCNAME%(0,8,1), 0
*SET,%_FNCNAME%(0,9,1), 0
*SET,%_FNCNAME%(0,10,1), 0
*SET,%_FNCNAME%(0,11,1), 0
*SET,%_FNCNAME%(0,12,1), 0
*SET,%_FNCNAME%(0,13,1), 0
*SET,%_FNCNAME%(0,14,1), 0
*SET,%_FNCNAME%(0,15,1), 0
! End of equation: {Y}
!
! Begin of equation: (2e7)/(2*{PI})*sin(2*6)*cos(2*{Y})
*SET,%_FNCNAME%(0,0,2), 6, -999
*SET,%_FNCNAME%(2,0,2), 0.0
*SET,%_FNCNAME%(3,0,2), 0.0
*SET,%_FNCNAME%(4,0,2), 0.0
*SET,%_FNCNAME%(5,0,2), 0.0
*SET,%_FNCNAME%(6,0,2), 0.0
*SET,%_FNCNAME%(0,1,2), 1.0, -1, 0, 2, 0, 0, 0
*SET,%_FNCNAME%(0,2,2), 0.0, -2, 0, 3.14159265358979310, 0, 0, -1
*SET,%_FNCNAME%(0,3,2), 1, -3, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,4,2), 0.0, -1, 0, 2e7, 0, 0, -3
*SET,%_FNCNAME%(0,5,2), 0.0, -2, 0, 1, -1, 4, -3
*SET,%_FNCNAME%(0,6,2), 0.0, -1, 0, 2, 0, 0, 0
*SET,%_FNCNAME%(0,7,2), 0.0, -3, 0, 6, 0, 0, -1
*SET,%_FNCNAME%(0,8,2), 0.0, -4, 0, 1, -1, 3, -3
*SET,%_FNCNAME%(0,9,2), 0.0, -1, 9, 1, -4, 0, 0
*SET,%_FNCNAME%(0,10,2), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,11,2), 0.0, -1, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,12,2), 0.0, -2, 0, 1, -1, 3, 3
*SET,%_FNCNAME%(0,13,2), 0.0, -1, 10, 1, -2, 0, 0
*SET,%_FNCNAME%(0,14,2), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,15,2), 0.0, 99, 0, 1, -2, 0, 0
! End of equation: (2e7)/(2*{PI})*sin(2*6)*cos(2*{Y})
!-->
FLST,2,2,4,ORDE,2
FITEM,2,9
FITEM,2,14
/GO
!*
!*
SFL,9,PRES, %F2%
SFL,14,PRES, %F2%
LSWRITE,3,
allsel,all
outpr,all,all !控制结果的输出
lssolve,1,3,1 !读入并求解多个载荷步
|