Page 3 - Demo

Planar Truss Examplefor ComrelAdd-onRCPConsult, 2023-2025Page 3The first set of statements fromFLIM(1)toDEFFUNC(1)operate with procedure PlanarTrussimplemented in FortranPlanarTrussGFOR.saoas export function that definition placed in source file PlanarTruss_GF.forwritten in the fixed format. In this sourcefile the values from global vector v_XP_and global nameswill be accessed. Note the large 2ndpart of the sourceto create various plots after the line !!*Visualization.!! Code based on Max Ehre work -https://github.com/maxehre/polynomial_surrogates!! Used with his kind permission as version from 2017 that adapted to Strurel.real*8functionPlanarTruss(ier)*bind(C,name='PlanarTruss')implicitreal*8(ah,oz)!! Include file with internal STRUREL definitionsinclude'strurel.hf'!GCC$ ATTRIBUTES DLLEXPORT::PlanarTrussinteger,intent(out)::ier!# ============================================================================!# Planar Truss -Example for FORTRAN!# ============================================================================!# Inputs:!# X = [A1, A2, E1, E2, P1, P2, P3, P4, P5, P6, u_y]!# E1,E2: Youngs modulus of horizontal and inclined bars resp. lognormally distributed!# A1,A2: cross section of horizontal and inclined bars resp. lognormally distributed!# P1P6: negative vertical forces attacking in nodes 813. Gumbel distributed!# u_y: maximal allowed vertical deflection in node 4. -constant!# Output:!# uout: vertical truss deflection at bottom center (N04)!# ============================================================================!# Truss description taken from Lee, S.H., Kwak, B.M. (2006).!# Response surface augmented moment method for efficient reliability analysis.!# Structural Safety 28(3), 261 272.!# https://www.sciencedirect.com/science/article/abs/pii/S0167473005000421!#!# Based on Diego Andr%u00e9s Alvarez Mar%u00edn work -!# https://github.com/diegoandresalvarez/elementosfinitos/tree/master/codigo/repaso_matricial/cercha_2d!# !# N: Nodes!# R: Rods!# ============================================================================!# N13__R4___N12__R8___N11__R12__N10__R16__N09__R20__N08!# /\\/\\/\\/\\/\\/\\!# / \\/ \\/ \\/ \\/ \\/ \\!# R1 R3 R5 R7 R9 R11 R13 R15 R17 R19 R21 R23!# / \\/ \\/ \\/ \\/ \\/ \\!# /___R2___\\/___R6___\\/__R10___\\/__R14___\\/__R18___\\/__R22___\\!# N01 N02 N03 N04 N05 N06 N07!# ||!# \\/!# uout => u_y @ N04!# ============================================================================real*8P(6)equivalence(P,v_XP_(6))!! FEM constantsparameter(nfe=23,ned=2,nnp=13,ndof=2*13,nbd=2*2,nc=3)parameter(PI=3.14159265d0)!! Element, nodes and dofs association!! IEN: connectivity matrix, nfe x nodes bar 1 has nodes 1 and 13, bar 2 has nodes 1 and 2 ...dimensionIEN(ned,nfe)dataIEN/1,13,1,2,13,2,13,12,2,12,2,3,12,3,12,11,3,11,*3,4,11,4,11,10,4,10,4,5,10,5,10,9,5,9,5,6,*9,6,9,8,6,8,6,7,8,7/!! Fixed dofsintegerbc(nc)databc/1,2,14/!! Internal datareal*8LM(nbd,nfe),theta(nfe+1),leng(nfe+1),l1,l2,*Area(nfe+1),E(nfe+1),ke(nfe),K(ndof,ndof),T(nbd,nbd),*KL(nbd,nbd),TT(nbd,nbd),KT(nbd,nbd),KG(nbd,nbd),*fc(ndof),ac(ndof),NN(nfe),ken(nbd),keT(nbd),de(nbd)!! Visualization datareal*8wo,we,ho,he,XY(ned,nnp),v(2,2),sc,fp(ned,nnp),*X(nbd),Y(nbd),XYrot(ned,nbd+1),minA,maxA,ddA,*XYten(ned+1,nbd+1)dataXY/0.d0,0.d0,4.d0,0.d0,8.d0,0.d0,12.d0,0.d0,16.d0,0.d0,*20.d0,0.d0,24.d0,0.d0,22.d0,2.d0,18.d0,2.d0,*14.d0,2.d0,10.d0,2.d0,6.d0,2.d0,2.d0,2.d0/integernfcharacter*200title,statecharacter*7cologicalInit!!---------------------------------------------------------------------------!! Planar truss model!!---------------------------------------------------------------------------!! Global flags to manage plot facilities of engineStrurelPlot=.TRUE.StrurelPlotName=\_GFOR_\StrurelPlotType=\StrurelPlotMode=3!! Vector inputs order from Stochastic ModelA1=v_XP_(1)! cross section of horizontal barsA2=v_XP_(2)! cross section of inclined barsE1=v_XP_(3)! Youngs modulus of horizontal bars