upgeom_upcood_悬索找形分析
如题，⽬前⼤家对初始缺陷命令upgeom,factor,LSTEP,SB STEP,FNAME,EXT中参数factor的理解，有两种看法：（1）factor 是施加载节点坐标上的因⼦，如factor是0.5，则施加位移值的⼀半到节点坐标上，这也是王新敏⽼师在《ANSYS⼯程数值分析》⼀书中对该参数的解释，⼤部分⼈也这么认为；（2）factor是位移值，单位与分析中采⽤的单位⼀致，如采⽤的单位为m，则当factor为0.5时，会对结构中产⽣的最⼤位移处施加0.5m的位移。
我通过实例认为第⼆种看法是对的，第⼀种看法也算对吧，但是需要做些补充。
大学生恋爱心理
为了验证以上两种看法的正确性，我从王新敏⽼师《ANSYS⼯程数值分析》的例⼦ex7.1c命令流如下，我在命令流的最后加了⼀⾏施加缺陷命令：单位为N,m
!EX7.1C 两端铰⽀柱特征值屈曲分析!
!BEAM188/189单元finish
/clear
/FILENAME,EX7-1
/prep7
b=0.03
h=0.05
l=3
e=2.1e11
et,1,beam189
mp,ex,1,e
mp,prxy,1,0.3 ctype,1,beam,rect cdata,b,h
k,1
k,2,,l
k,10,0,l/2,l/2
l,1,2
dk,1,ux,,,,uy,uz,roty dk,2,ux,,,,uz,roty!多施加了UY约束
latt,1,,1,,10,,1 lesize,all,,,20
lmesh,all
finish
/solu!第⼀次进⼊求解器fk,2,fy,-1000!（此荷载原为-1，我改为1000）pstres,on
solve
finish
/solu!第⼆次进⼊求解器进⾏弹性屈曲分析antype,buckle bucopt,lanb,1!模态改为⼀阶
mxpand,1!改为⼀阶outres,all,all
solve
finish
/POST1
SET,LIST
PLDISP
/PREP7
UPGEOM,1.5,,,EX7-1,RS T!施加了初始缺陷，不是1.5倍就是1.5m？
当第⼀次求解完毕后得到的节点位移变形，见图1（通过后处理general post>plot
月子可以喝牛奶吗results>contour
plot>nodal solu>dof solu>displacement
vector sum得到），之后进⾏了屈曲分析，直⾄施加完初始缺陷后，⽤GUI 中的replot重新显⽰了图形，见图2，⼜再次查看了节点位移变形图，见图3。图1的“最⼤值”（红⾊部分为0.952e-05），图3的“最⼤值”为0.955。取消呼叫转移怎么设置
这根杆件长为3m，⽽图2的跨中最⼤位移很像1.5m，因此我这么认为：有⼈曾说ansys对屈曲后的最⼤位移归⼀化，即最⼤位移不管数值时多少，都会显⽰为1，我认为factor实际上就是位移值，只不过做了归⼀化后，最⼤位移变为1，施加的位移是因⼦factor乘以这个最⼤位移1，其施加缺陷时，这个factor就相当于数值，不知道我理解的对否？这个例⼦中，第⼀次加载位移为0.952e-5，加缺陷后位移变为0.955，显然不是1.5倍。
当然这也不是说第⼀种解释不对，因为最⼤位移归⼀化了，所以这两种说法可以都是对的。只不过便于我这样的初学者理解，这两种说法应做些详细的补充。此外我还做了另外⼀个例⼦的施加初始缺陷的验证，这个例⼦是王新敏⽼师书中EX8.8，命令流如下：
!EX8.8 中⼼受压铰接柱的⼏何⾮线性分析
finish
/clear
/filname,colu
/prep7
aa=100.0
ai=10000/12.0
l0=1000
em=2e5
pcr=acos(-1)*acos(-1)* em*ai/l0/l0
et,1,beam3考研英语和六级英语哪个难
mp,ex,1,em
mp,prxy,1,0
r,1,aa,ai,10
k,1,0
k,2,0,l0/2
k,3,0,l0
l,1,2
l,2,3
lesize,all,,,20
lmesh,all
node1=node(0,l0,0) finish
/solu
dk,1,ux,,,,uy
dk,3,ux
fk,3,fy,-pcr
pstres,on
solve
finish
/solu
antype,1
bucopt,lanb,1 mxpand,1,,,1
solve
finish
/prep7
upgeom,1,,,colu,rst!施加初始缺陷，我将其改为1000后发现位移好像就是1000mm
finish
EX8.8在最后施加初始缺陷时王⽼师这么注明：“施加100%缺陷，柱中1mm”，分析中单位是mm。我后来将“1”改为#2 2010-9-17 12:35
#3 2010-9-17 12:37
#4 2010-9-17 12:39
#5 2010-9-17 17:19 #6 2010-9-17 17:28
求助（单索迭代求解问题）
本⼈现在在研究单索的⼏何⾮线性
静⼒分析，预应⼒⼏何⾮线性模态
分析，单索的预应⼒谐响应分析，
单索的预应⼒瞬态分析，遇到了⼀
些问题，特向各位⽼师请教：
单索的⼏何⾮线性静⼒分析
（1）。有预应⼒的拉索，在什么
情况下，要考虑⾃重引起跨中垂度
对结构刚度的影响，以及对索⼒的
影响？考虑⾃重影响时，是⽤命令（acel，，1）还是（acel，，9.8 ），
西加一笔是什么字
这个问题三⽉⾬⽼师的<
程结构数值分析>>p469⾯⽤的是（acel，，1）。⽽518⾯⽤的是（acel，，9.8 ）为何如此？
（2）。如果我已知道（预张⼒）
索⼒，由于⾃重的影响，预张⼒应
该不等于索⼒？且索⼒是变化的，
预张⼒应该等于索端单元的索⼒？
（3）。已知道预张⼒，如何迭代
求解所得初始构型，为何我的迭代
程序，怎么计算都不能使索端单元
的索⼒收敛到给定的预张⼒，（其中
初应变，我采⽤初应⼒/E/A?，我的
计算程序请看附件，急切希望各位
热⼼的⽼师给予指导，万分感激！）（4）
预应⼒⼏何⾮线性模态分析
当我完成了结构的⼏何⾮线性
静⼒分析之后，就要求解在⼤变形
的情况下，单索的⾃振频率等，⽤
于瞬态分析时计算阻尼项。
（1）我如何才能把通过⼏何⾮线性静⼒分析，迭代得到的单索的初始位置和初始应⼒，以及刚度矩
器欲难量
阵等信息，引⼊到我的预应⼒⼏何⾮线性模态分析中，如何才能写⼊这些信息，求助？或者，我直接在完成找形⼯作之后，再次进⾏静⼒分析（打开pstres，on，写⼊emat ⽂件，ematwrite，yes）之后就是进⾏预应⼒模态分析（利⽤psolve 命令求解），这样可⾏否？为什么我运⾏的时候出现错误？（由于模态分析是⾮线性的，会不会有冲突，这个问题我⼀直在思考，希望专家给予指导）
冬季施工最后⼀个问题，如果我的单索是位于⽔中（深海中），在计算静⼒的时候，如何考虑浮⼒？是直接采⽤惯性荷载还是把浮⼒当成外载荷，静⼒求解之后，再施加上去？？⽔的阻尼对单索的模态分析，到底有多⼤影响，能否忽略？里脊肉的做法
我的单索找形程序，我采⽤的是直接迭代法（<>p469的⽅法。）
我的迭代程序如下：
