| 查看: 720 | 回复: 0 | ||
| 【悬赏金币】回答本帖问题，作者izaye将赠送您 10 个金币 | ||
[求助]
关于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 |
回复此楼» 猜你喜欢
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有27人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有4人回复
-
UV压敏胶开发
已经有6人回复
-
有没有H口的？有收到消息的吗？
已经有3人回复
-
十年后又回来了，论文投稿求助
已经有3人回复
-
纯娱乐，不喜欢勿喷
已经有10人回复
-
面上提前没消息，有中的吗
已经有18人回复
-
娱乐
已经有4人回复
-
系统今天又提示维护了，估计离放榜不远了
已经有16人回复
-
一个有机合成实验室都需要哪些设备？
已经有10人回复