找回密码
 注册
Simdroid-非首页
查看: 574|回复: 9

在瞬态分析中加载正弦位移载荷

[复制链接]
发表于 2006-5-9 22:56:27 | 显示全部楼层 |阅读模式 来自 湖南株洲
如题,如何在瞬态分析中加载有预应力的正弦位移载荷(S=S0+S'sin(wt))呢?使用apply/function吗?

谁有相关的命令流可以借鉴一下吗?
谢谢!
发表于 2006-5-10 07:46:59 | 显示全部楼层 来自 北京
Simdroid开发平台
命令流太乱,直接说吧:在parameter/function/edit里定义一个函数
如a+b*sin(w*TIME),然后save成一个文件,
在function/read里读进来,填写变量名S,坐标系,系数a、b、w的数据,就行了,生成一个类似于表的东西,加载和的时候加这个表S就行了!

*SET,_FNCNAME,'S'   
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*SET,_FNC_C1(1),10  
*SET,_FNC_C2(1),10  
*SET,_FNC_C3(1),10  
*SET,_FNCCSYS,0
! /INPUT,..\..\a.func,,,1   
*DIM,%_FNCNAME%,TABLE,6,5,1,,,,%_FNCCSYS%   
!   
! Begin of equation: a+b*sin(w*{TIME})  
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 1, 19, 3, 1
*SET,%_FNCNAME%(0,2,1), 0.0, -1, 9, 1, -1, 0, 0
*SET,%_FNCNAME%(0,3,1),   0, -2, 0, 1, 18, 3, -1
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, 17, 1, -2
*SET,%_FNCNAME%(0,5,1), 0.0, 99, 0, 1, -1, 0, 0
! End of equation: a+b*sin(w*{TIME})
!-->
 楼主| 发表于 2006-5-10 14:13:09 | 显示全部楼层 来自 湖南株洲
可是我这样做出来的并非正弦载荷呢? 附件fun2 result是节点1的位移随时间变化曲线(应该是正弦线的),fun2 array是数组s vs.time,按理说这两种曲线应该表示的是一样的结果吧,为什么现实会不一致呢?

说明:fun2 result是在post26中得到的,fun2 array是在plot/array parameters中得到的(*VPLOT,time,s, , , , , , , ,)
*DEL,s  
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*SET,_FNCNAME,'s'   
! /INPUT,.\APDL\fun2.func   
*DIM,%_FNCNAME%,TABLE,6,12,1
!   
! Begin of equation: 0.5/1000*sin(2*3.14*5*{TIME})  
*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, 0.5, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1000, 0, 0, -1  
*SET,%_FNCNAME%(0,3,1),   0, -3, 0, 1, -1, 4, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 2, 0, 0, 0  
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 3.14, 0, 0, -1  
*SET,%_FNCNAME%(0,6,1), 0.0, -4, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,7,1), 0.0, -1, 0, 5, 0, 0, -4
*SET,%_FNCNAME%(0,8,1), 0.0, -2, 0, 1, -4, 3, -1
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 1, -2, 3, 1
*SET,%_FNCNAME%(0,10,1), 0.0, -1, 9, 1, -1, 0, 0
*SET,%_FNCNAME%(0,11,1), 0.0, -2, 0, 1, -3, 3, -1   
*SET,%_FNCNAME%(0,12,1), 0.0, 99, 0, 1, -2, 0, 0
! End of equation: 0.5/1000*sin(2*3.14*5*{TIME})
!-->
!*  
*GET,s,NODE, ,U,Z   
*TREAD,S, , ,, ,
/solu
ANTYPE,4
TRNOPT,FULL
LUMPM,0
NSUBST,10,0,0   
PSTRES,1
TIME,1  
/GO
D,1, , %S% , , , ,UZ,
solve
/post26
NSOL,2,1,U,Z, UZ_2  
XVAR,1  
PLVAR,2,
 楼主| 发表于 2006-5-10 14:17:41 | 显示全部楼层 来自 湖南株洲
附件

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
发表于 2006-5-11 10:38:43 | 显示全部楼层 来自 黑龙江哈尔滨
C*** START DEFINE THE VOLTA AND VOLTB
*AFUN,RAD
*SET,PI,3.1415926

*SET,RFREQ,20000   !!!!!!!!!
*SET,T_P,1/RFREQ   !ZHOU QI (S)

*SET,VOLTARMS,100   !!!!!ROOT MEAN SQUARE VALUE OF VOLTA
*SET,VOLTBRMS,100   !!!!!!ROOT MEAN SQURAE VALUE OF VOLTB
*SET,VLTA,VOLTARMS*SQRT(2)      !PEAK VALUE OF VOLTA
*SET,VLTB,VOLTBRMS*SQRT(2)      !PEAK VALUE OF VOLTB

*SET,PH_D_D,90          !!!!!!!!!!PHASE DIFFERENCE OF A AND B IN DEGREE
*SET,PH_D,PH_D_D/180*PI !PHASE DIFFERENCE IN RADIAS
*SET,T_S,PH_D_D/360*T_P !PHASE DIFFERENCE TIME SPAN

*SET,PRD_N,20         !PERIOD NUMBER OF SIN INPUT SIGNAL
           
*SET,PRT_N,72         !PARTATION NUMBER IN ONE CRICLE
*SET,T_ICR,T_P/PRT_N    ! TIME INCREMENT

*SET,LST_TIME,T_P*PRD_N  !SETS THE TIME FOR A LOAD STEP:TIME,LST_TIME

*SET,LSTN_1C,16       !LOAD SUBSTEP NUMBER IN ONE CIRCLE
*SET,LSTN_W,LSTN_1C*PRD_N          !THE WHOLE LOAD SUBSTEP NUMBER




*DIM,V_A,TABLE,1+PRT_N*PRD_N,1,1
*SET,V_A(0,1,1),0
*SET,V_A(1,0,1),0
*SET,V_A(1,1,1),0
*DO,I,1,PRD_N
  *DO,J,1,PRT_N
    *SET,V_A((I-1)*PRT_N+J+1,0,1),(I-1)*T_P+J*T_ICR
  *ENDDO
   
  *DO,K,1,PRT_N
    *SET,V_A((I-1)*PRT_N+K+1,1,1),VLTA*SIN((2*PI/T_P)*K*T_ICR)
  *ENDDO
*ENDDO



*DIM,V_B,TABLE,1+PRT_N*PRD_N,1,1
*SET,V_B(0,1,1),0
*SET,V_B(1,0,1),0
*SET,V_B(1,1,1),VLTB*SIN(PH_D)
*DO,I,1,PRD_N
  *DO,J,1,PRT_N
    *SET,V_B((I-1)*PRT_N+J+1,0,1),(I-1)*T_P+J*T_ICR
  *ENDDO
   
  *DO,K,1,PRT_N
    *SET,V_B((I-1)*PRT_N+K+1,1,1),VLTB*SIN((2*PI/T_P)*(K*T_ICR+T_S))
  *ENDDO
*ENDDO


*VPLOT,V_A(1,0,1),V_A(1,1,1)  

*VPLOT,V_B(1,0,1),V_B(1,1,1)




这个是两相相差90度的简谐激励信号,你可以参考一下。
发表于 2006-5-13 08:54:19 | 显示全部楼层 来自 北京

回复 #3 joannall 的帖子

我生成的是正弦曲线,你的怎么不是呢?
重新加一遍表载荷,时间步长取短点,输出每一个时间步的结果!!
发表于 2006-5-13 09:03:02 | 显示全部楼层 来自 北京

正弦位移

位移-时间曲线

[ 本帖最后由 gogouter 于 2006-5-13 09:05 编辑 ]

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
发表于 2010-4-13 17:01:47 | 显示全部楼层 来自 吉林长春
好像对我很有帮助!
回复 不支持

使用道具 举报

发表于 2011-1-20 19:34:54 | 显示全部楼层 来自 辽宁沈阳
确实对我很大帮助。
回复 不支持

使用道具 举报

发表于 2011-2-28 12:32:18 | 显示全部楼层 来自 山西太原
请教一下各位大侠,如果要在网架上加载地震波该使用说明命令呢
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-28 18:06 , Processed in 0.046047 second(s), 13 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表