Lbase: Lbase
In ShinsukeSakai0321/LimitState: Package supporting FORM approach

Description Usage Public function Author(s) See Also Examples

Base class for the description of limit state function

1	Arbitrary limit state functions defined by users are inherited from this base class

Standard reference class is R6 class
Public function

function	contents	variable(s)
initialize(n)	initializes the object	n:number of variables
GetN()	returns the number of variables
GetX()	returns the variable list X
SetX(x)	sets variable list x to the variable list X	x: variable list
GetG()	returns calculated g value
SetG(g)	sets calculated g-value to the variable	g:calculated g
GetdGdX()	returns differential variable value list dGdX
SetdGdX(x)	sets differintial variable list dGdX	x: dGdX list

Shinsuke Sakai, Takuyo Kaida

As the application of the following example of GPRSMED, see example of help command LSFM

# Example
#"Probability, Reliability and Statistical Methods in
#  Engineering Design"
# Achintya Haldar  & Sankaran Mahadevan
#  P.218 Table 7.5, 7.6
GPRSMED <- R6::R6Class("GPRSMED",
                   inherit=Lbase,
                   public = list(
                     gcalc = function(){
                       X <- super$GetX()
                       As <- X[1]
                       fy <- X[2]
                       fcd <- X[3]
                       b <- X[4]
                       d <- X[5]
                       eta <- X[6]
                       M <- X[7]
                       g <- As*fy*d*(1.0-eta*As*fy/b/d/fcd)-M
                       super$SetG(g)
                     } ,
                     dGdXcalc = function(){
                       X <- super$GetX()
                       dGdX <- super$GetdGdX()
                       As <- X[1]
                       fy <- X[2]
                       fcd <- X[3]
                       b <- X[4]
                       d <- X[5]
                       eta <- X[6]
                       M <- X[7]
                       dGdX[1] <- fy * d * (1 - eta * As * fy/b/d/fcd) - As * fy * d * (eta * fy/b/d/fcd)
                       dGdX[2] <- As * d * (1 - eta * As * fy/b/d/fcd) - As * fy * d * (eta * As/b/d/fcd)
                       dGdX[3] <- As * fy * d * (eta * As * fy/b/d/fcd^2)
                       dGdX[4] <- As * fy * d * (eta * As * fy/b^2/d/fcd)
                       dGdX[5] <- As * fy * (1 - eta * As * fy/b/d/fcd) + As * fy * d * (eta *
                                                                                                   As * fy/b/d^2/fcd)
                       dGdX[6] <- -(As * fy * d * (As * fy/b/d/fcd))
                       dGdX[7] <- -1
                       super$SetdGdX(dGdX)
                     }
                   )

)