24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者izaye将赠送您 10 个金币

izaye

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 6.8
帖子: 10
在线: 15.3小时
虫号: 7974210
注册: 2018-02-08
专业: 岩土与基础工程

[求助] 关于abaqus子程序中精度的问题？

大家能不能帮忙看一我这个abaqus里的vuamp子程序的精度是否正确，如果不正确的话该如何声明变量？纯新手，这两天查资料查的很多概念都蒙圈了，求大佬帮忙
   SUBROUTINE VUAMP(
   *    ampName, time, ampValueOld, dt, nprops, props, nSvars,
   *    svars, lFlagsInfo, nSensor, sensorValues, sensorNames,
   *    jSensorLookUpTable,
   *    AmpValueNew,
   *    lFlagsDefine,
   *    AmpDerivative, AmpSecDerivative, AmpIncIntegral)

   INCLUDE 'VABA_PARAM.INC'

C    time indices
   parameter (iStepTime       = 1,
   *          iTotalTime    = 2,
   *          nTime          = 2)
C    flags passed in for information
   parameter (iInitialization = 1,
   *          iRegularInc    = 2,
   *          ikStep          = 3,
   *          nFlagsInfo       = 3)
C    optional flags to be defined
   parameter (iComputeDeriv    = 1,
   *          iComputeSecDeriv  = 2,
   *          iComputeInteg    = 3,
   *          iStopAnalysis    = 4,
   *          iConcludeStep    = 5,
   *          nFlagsDefine    = 5)
   dimension time(nTime), lFlagsInfo(nFlagsInfo),
   *       lFlagsDefine(nFlagsDefine),
   *       sensorValues(nSensor),
   *       props(nprops),
   *       sVars(nSvars)

   character*80 sensorNames(nSensor)
   character*80 ampName
   dimension jSensorLookUpTable(*)

      tim=time(iStepTime)
      if (ampName(1:6).EQ.'sinsin') then
      r=mod(tim,5.76)
      if(r.LE.0.25) then
         ampValueNew=1000000+7580*sind(8*pi*tim)
      else
         ampValueNew=0
      end if
      end if


      return
      end

回复此楼

» 猜你喜欢

1楼 2020-08-27 10:30:27

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 izaye 的主题更新

返回列表

不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)

普通表情龙兔虎猫