静态分析后进行瞬态分析，如何设置

MILAN3

问题：重力引起一块薄板初始变形，一脉冲载荷作用于薄板，引起薄板振动，脉冲作用很短时间后撤销，重力一直存在

分析需求：先静态分析，求出重力引起的初始变形，后进行瞬态分析，求脉冲载荷和重力共同作用下的振动，瞬态分析需在初始变形的基础上开始

请高手讲一下这种分析如何设置？

－－－－－－－－－－－－－－－－－－－－

下面是个人的两种分析方案，好像都不正确

方案一：先进行静态计算，后进行瞬态计算，结果表明瞬态分析开始时的变形不是静分析的初始变形

/SOLU

ANTYPE,STATIC        !静态计算

ACEL,0,GY,GZ        !施加重力,一直延续至后续的瞬态计算中

SOLVE

FINI

/SOLU

ANTYPE,TRANS        !瞬态计算

TIME,1E-9                !脉冲载荷开始时间

AUTOTS,ON               

KBC,1                       

F,LOAD_CRY_FRO,FZ,-1*(FORCE/LOAD_NODE_NUM)       

LSWRITE,1

TIME,4E-9                !脉冲载荷撤销时刻，共持续了4纳秒

FDELE,LOAD_CRY_FRO,FZ

LSWRITE,2

TIME,2                        !重力作用下振动时间，持续2秒针

LSWRITE,3

LSSOLVE,1,2,3

－－－－－－－－－－－－－－－－－

方案二：进行瞬态计算，通过关闭和打开瞬态效应开关考虑重力效应，结果表明没算出振动效应

/SOLU

ANTYPE,STRANS

TIMINT,OFF                !关闭瞬态效应

TIME,1E-3                !小作用时间

ACEL,0,GY,GZ                !施加重力

NSUBST,2                !两个子步

KBC,1                        !阶跃加载

SOLVE       

TIMINT,ON                !打开瞬态效应       

TIME,(1E-3)+(1E-9)        !脉冲载荷开始时间

AUTOTS,ON               

KBC,1                       

F,LOAD_CRY_FRO,FZ,-1*(FORCE/LOAD_NODE_NUM)       

LSWRITE,1

TIME,(1E-3)+(4E-9)        !脉冲载荷撤销时刻，共持续了4纳秒

FDELE,LOAD_CRY_FRO,FZ

LSWRITE,2

TIME,(1E-3)+2                !重力作用下振动时间，持续2秒

LSWRITE,3

LSSOLVE,1,2,3

htbbzzg

应该采用第二种方法。不过命令流中的脉冲载荷不知道是哪一种脉冲？看起来似乎是矩形脉冲。如果是其它脉冲方式，需要进行说明和相应的修改。
另外，你这段命令流好像有一些问题，我修改了一下，见下面，其中主要修改的地方用 !*** 做了标注：

/SOLU

ANTYPE,STRANS

TIMINT,OFF                !关闭瞬态效应

TIME,1E-3                !小作用时间

ACEL,0,GY,GZ                !施加重力

NSUBST,2                !两个子步

KBC,1                        !阶跃加载

SOLVE
      

TIMINT,ON                !打开瞬态效应      

NSUBST,5， 100， 5                !*** 至少 5 个子步

TIME,(1E-3)+(1E-9)        !脉冲载荷开始时间

AUTOTS,ON               

KBC,0                     !***  斜坡方式加载                       

F,LOAD_CRY_FRO,FZ,-1*(FORCE/LOAD_NODE_NUM)      

LSWRITE,1


TIME,(1E-3)+(4E-9)        !*** 脉冲载荷开始撤销时刻，撤销前，脉冲载荷保持不变 4 纳秒，为矩形脉冲

LSWRITE,2


FDELE,LOAD_CRY_FRO,FZ     !*** 撤销脉冲载荷，仍为斜坡方式加载。如果要使用阶梯方式，此处需要改用用  kbc, 1, 下一步再改回 kbc,0

TIME,(1E-3)+(5E-9)        !*** 脉冲载荷撤销到0，用时 1 纳秒

LSWRITE,3                 !*** 增加一个撤销脉冲载荷的载荷步


TIME,(1E-3)+2                !重力作用下振动时间，持续2秒

LSWRITE,4

LSSOLVE,1,4,1             !***  原来的命令用法错误

没有验证过，你试试看。有问题你可以自己改改，或进一步讨论。

MILAN3

htbbzzg 发表于 2023-12-17 09:13
应该采用第二种方法。不过命令流中的脉冲载荷不知道是哪一种脉冲？看起来似乎是矩形脉冲。如果是其它脉冲方 ...

非常感谢。
载荷是矩形脉冲。
按您的修改我重新编一下程序，有了结果后我再来回复一下。

MILAN3

htbbzzg 发表于 2023-12-17 09:13
应该采用第二种方法。不过命令流中的脉冲载荷不知道是哪一种脉冲？看起来似乎是矩形脉冲。如果是其它脉冲方 ...

按您的意见修改了计算命令流，计算出了重力引起的初始变形，也算出了后续的振动效果。

还有一个问题请教，关于重力的作用时长，TIMINT,OFF状态下重力作用多长时间比较合适？会不会有这种情形，重力作用时间过长或过短，会使重力引起的效果增加或减弱，即命令流中的（1E-3）秒这一时长如何设置一个合理的值？

htbbzzg

重力计算部分，由于关闭了瞬态效应，所以是一个静力分析，理论上与使用的时间步长是没有关系的。至于设置为多长合适，要看你的具体问题。一般来说，取你的瞬态分析的第一个步长的 1% 或更小一点，至少两个子步即可。不会影响后续的瞬态计算。