pipe288 单元泥面处 轴力发生突变

zhuizhu18

图1：泥面处轴力发生突变（泥面以上有波流力，泥面一下无波流力）
        
图2：不加波流力轴力光滑过渡，轴力为重力引起
      
命令流：

     用pipe288单元建立了一个海上的桩基基础，桩顶高程为2.0m，泥面处为-19.0m，桩底为-21.0m，桩底固接。泥面以上收到波流力的作用，泥面以下不受波流力作用。
    但是为啥泥面处的轴力会发生突变呢？
    按理说波流力产生的是水平方向的力，不应该对轴力有太大影响啊。就算是有影响也不应该在你泥面处发生突变啊。到底是哪个地方出现错误了啊？还望大虾指点啊

     FINISH
/CLEAR
/PREP7
/UNITS,SI
!**********定义单元以及材料属性*****************
ET,1,PIPE288       !pipe288 单元要求将坐标原点建立在水平面位置
KEYOPT,1,3,3       !沿长度方向的形函数  选三次型函数
KEYOPT,1,4,1       !环形应变处理方法  设为薄壁管 为沿厚度为平面应变
MP,EX,1,2.1E11   !定义材料1 钢管的属性
MP,PRXY,1,0.3
MP,DENS,1,7850
MP,DENS,10,1030  !定义材料2 钢管内部液体的密度
MP,DENS,11,1030  !定义外部液体 也就是海水的密度
!**********定义ID*****************
Sec_pile_up_mud=1
Id_wat=2
Id_cur=3
Id_wav=4
Sec_pile_down_mud=10

!**********定义泥面以上桩的管单元的的截面属性 为水动力单元*****************
sectype,Sec_pile_up_mud,pipe,,pile_up
secdata,0.80,0.01,8,,3,10
socean,2     !将ocean id与通过sectype最新定义的的单元截面联系起来
octype,Id_wat,basic,,Id_cur,Id_wav
ocdata,19.0,11,0,-1,1     ! Mat_out_water为外部海水的材料参数编号  考虑浮力作用
                       
octable,,1.2,1.2,,2.0,2.0
octype,Id_cur,curr         !潮流力参数   用octable定义潮流位置 以及 相应位置处的流速作用角度
octable,0,2.0,0,,
octable,-10,2.0,0,,
octable,-19,2.0,0,,
octype,Id_wav,wave
ocdata,1,0,0,0,,,         !波浪较小 所以选择有修正的微波幅理论  并定义波浪作用方向以及 波浪位置类型
octable,1.0,10.0,0
!**********定义泥面以下桩的截面属性  此处不受潮流以及波浪的做用 因此不用设置水动力属性****************
sectype,Sec_pile_down_mud,pipe,,pile_down
secdata,0.80,0.01,8,,3,10
!建立模型  只建立上部结构的模型 桩土相互作用的模型不在此处建立
K,1,0,0,2.0        
K,2,0,0,-19.0
K,3,0,0,-21.0
!创立桩的桩线
l,1,2
LATT,1,,1,,,,Sec_pile_up_mud
ALLSEL
LESIZE,ALL,0.2
LMESH,ALL
!建立下部的桩土   然后耦合连接处的自由度
LSEL,ALL
LSEL,INVE
L,2,3
LATT,1,,1,,,,Sec_pile_down_mud
LESIZE,ALL,0.2
LMESH,ALL
ALLSEL
!施加约束
KSEL,S,LOC,Z,-21.0
NSLK,S
D,ALL,all
ALLSEL
/SOLU
ANTYPE,0
TIME,1
NROPT,FULL
ACEL,,,9.8
SOLVE
/post1
ETABLE,Fx_i,SMISC,1     !单元I点轴力
ETABLE,Fx_j,SMISC,14     !单元J点轴力
plls,Fx_i,Fx_j
                    

sunnygu

我问过技术支持，的确是BUG。尽量减少用这个单元吧
我再想想我当时怎么解决的。

zhuizhu18

不知道这个是这个单元的bug还是自己的命令有问题。

zhuizhu18

sunnygu 发表于 2012-11-17 05:34
我问过技术支持，的确是BUG。尽量减少用这个单元吧
我再想想我当时怎么解决的。 ...

单元bug问题的话，我想可以泥面以下也可以设置ocean参数，只要把ocean的各参数，如波速，流速等设置成小值，这样的话就可以实现轴力连续了。只是浮力好像没起到什么作用，跟手算不考虑浮力的结果差不多。



!**********定义泥面以下桩的截面属性 将水动力参数设置成小值 不考虑波流力的作用
Sec_pile_down_mud=10
Id_wat_1=5
Id_cur_1=6
Id_wav_1=7
sectype,Sec_pile_down_mud,pipe,,pile_down
secdata,0.80,0.01,8,,3,10
socean,Id_wat_1     !将ocean id与通过sectype最新定义的的单元截面联系起来
octype,Id_wat_1,basic,,Id_cur_1,Id_wav_1
ocdata,21.0,11,0,-1,1      ! Mat_out_water为外部海水的材料参数编号 考虑浮力的作用
octable,,1.2,1.2,,2.0,2.0
octype,Id_cur_1,curr       !潮流力参数   用octable定义潮流位置 以及 相应位置处的流速作用角度
octable,0,0.01,0,,
octable,-10,0.01,0,,
octable,-21.0,0.01,0,,     !不考虑波流作用 将此处的流速设置为0就可以了
octype,Id_wav_1,wave
ocdata,1,0,0,0,,,        !波浪较小 所以选择有修正的微波幅理论  并定义波浪作用方向以及波浪位置类型
octable,0.01,10.0,0       !波高必须为一个正值 所以此处必须设置波高

lostangelee

学长太厉害了，完全看不懂。。

坐等树上掉苹果

为什么  使用的ari波 没有多设置几个浪向